def FinalizeSolutionStep(self): super(ShallowWaterSolver, self).FinalizeSolutionStep() epsilon = max(self.advection_epsilon, self.main_model_part.ProcessInfo[SW.DRY_HEIGHT]) SW.ShallowWaterUtilities().ComputeHeightFromFreeSurface( self.main_model_part) SW.ComputeVelocityProcess(self.main_model_part, 1e-3).Execute() SW.ShallowWaterUtilities().ComputeAccelerations(self.main_model_part)
def FinalizeSolutionStep(self): super().FinalizeSolutionStep() KM.VariableUtils().SetFlag(KM.ACTIVE, True, self.main_model_part.Nodes) KM.VariableUtils().SetFlag(KM.ACTIVE, True, self.main_model_part.Elements) dry_height = self.main_model_part.ProcessInfo.GetValue(SW.DRY_HEIGHT) SW.ShallowWaterUtilities().ComputeFreeSurfaceElevation(self.main_model_part) SW.ShallowWaterUtilities().ResetDryDomain(self.main_model_part, dry_height) SW.ComputeVelocityProcess(self.main_model_part, dry_height).Execute() self._CheckWaterLoss()
def FinalizeSolutionStep(self): super().FinalizeSolutionStep() epsilon = max(self.advection_epsilon, self.GetComputingModelPart().ProcessInfo[SW.DRY_HEIGHT]) SW.ShallowWaterUtilities().ComputeHeightFromFreeSurface( self.GetComputingModelPart()) SW.ComputeVelocityProcess(self.GetComputingModelPart(), epsilon).Execute() SW.ShallowWaterUtilities().ComputeAccelerations( self.GetComputingModelPart())
def FinalizeSolutionStep(self): super(LagrangianShallowWaterSolver, self).FinalizeSolutionStep() KM.VariableUtils().SetFlag(KM.ACTIVE, True, self.GetComputingModelPart().Elements) SW.ShallowWaterUtilities().ResetDryDomain(self.GetComputingModelPart(), self.thickness) SW.ShallowWaterUtilities().ComputeFreeSurfaceElevation( self.GetComputingModelPart()) SW.ComputeVelocityProcess(self.GetComputingModelPart(), 1e-1).Execute() SW.ShallowWaterUtilities().ComputeAccelerations( self.GetComputingModelPart())