def __init__(self, model, project_parameters, field_utility, fluid_solver, dem_solver, variables_manager):
super().__init__(model, project_parameters, field_utility, fluid_solver, dem_solver, variables_manager)
self.frame_angular_vel = Vector([0, 0, self.project_parameters['frame_of_reference']["angular_velocity_of_frame_Z"].GetDouble()])
self.omega = self.project_parameters['frame_of_reference']["angular_velocity_of_frame_Z"].GetDouble()
candelier_pp.include_lift = PT.RecursiveFindParametersWithCondition(self.project_parameters["properties"],
'vorticity_induced_lift_parameters',
condition=lambda value: not (value['name'].GetString()=='default'))
def __init__(self, project_parameters, model_part):
L2_projection_recoverer.L2ProjectionGradientRecoverer.__init__(self, project_parameters, model_part)
self.element_type = "ComputeGradientPouliot20123D"
self.condition_type = "ComputeLaplacianSimplexCondition3D"
self.FillUpModelPart(self.element_type, self.condition_type)
self.DOFs = (Kratos.VELOCITY_COMPONENT_GRADIENT_X, Kratos.VELOCITY_COMPONENT_GRADIENT_Y, Kratos.VELOCITY_COMPONENT_GRADIENT_Z)
self.AddDofs(self.DOFs)
self.calculate_vorticity = (project_parameters["vorticity_calculation_type"].GetInt() > 0
or PT.RecursiveFindParametersWithCondition(project_parameters["properties"],
'vorticity_induced_lift_parameters'))
def __init__(self, project_parameters, model_part):
Pouliot2012EdgeDerivativesRecoverer.__init__(self, project_parameters, model_part)
self.element_type = self.GetElementType(project_parameters)
self.FillUpModelPart(self.element_type)
self.DOFs = self.GetDofs(project_parameters)
self.AddDofs(self.DOFs)
self.calculate_vorticity = (project_parameters["vorticity_calculation_type"].GetInt() > 0
or PT.RecursiveFindParametersWithCondition(project_parameters["properties"],
'vorticity_induced_lift_parameters'))
self.domain_size = project_parameters["fluid_parameters"]["solver_settings"]["domain_size"].GetInt()
def SetHistoryForceOptions(self):
self.history_force_on = False
self.MAE_parameters = Parameters("{}")
for prop in self.project_parameters["properties"].values(): #TODO: now it only works for one property!
self.history_force_on = (PT.RecursiveFindParametersWithCondition(
self.project_parameters["properties"], 'history_force_parameters',
condition=lambda value: value['name'].GetString() != 'default'))
if self.history_force_on:
self.MAE_parameters = prop["hydrodynamic_law_parameters"]["history_force_parameters"]["mae_parameters"]
break
self.do_use_mae = PT.RecursiveFindTrueBoolInParameters(self.MAE_parameters, 'do_use_mae')
def __init__(self, project_parameters, model_part):
recoverer.DerivativesRecoverer.__init__(self, project_parameters, model_part)
self.dimension = project_parameters["fluid_parameters"]["solver_settings"]["domain_size"].GetInt()
self.model_part = model_part
self.use_lumped_mass_matrix = project_parameters["material_acceleration_calculation_type"].GetInt() == 3
self.recovery_model_part = ModelPart("PostGradientFluidPart")
self.custom_functions_tool = SDEM.CustomFunctionsCalculator3D()
self.calculate_vorticity = (project_parameters["vorticity_calculation_type"].GetInt() > 0
or PT.RecursiveFindParametersWithCondition(project_parameters["properties"],
'vorticity_induced_lift_parameters'))
self.GetFieldUtility()
self.CreateCPluPlusStrategies()
def SetOptions(self, parameters):
self.do_include_history_force = (PT.RecursiveFindParametersWithCondition(
parameters["properties"], 'history_force_parameters',
condition=lambda value: value['name'].GetString() != 'default'))
if self.do_include_history_force: #TODO: extend to multiple properties
for prop in parameters["properties"].values():
self.history_force_parameters = prop["hydrodynamic_law_parameters"]["history_force_parameters"]
break
self.do_backward_coupling = parameters["coupling"]["coupling_level_type"].GetInt() > 1
self.backward_coupling_tools = parameters["coupling"]["backward_coupling"]
def __init__(self, project_parameters, model_part):
recoverer.DerivativesRecoverer.__init__(self, project_parameters,
model_part)
self.model_part = model_part
self.use_lumped_mass_matrix = project_parameters[
"material_acceleration_calculation_type"].GetInt() == 3
self.recovery_model_part = ModelPart("PostGradientFluidPart")
self.custom_functions_tool = SDEM.CustomFunctionsCalculator3D()
self.calculate_vorticity = (
project_parameters["vorticity_calculation_type"].GetInt() > 0
or PT.RecursiveFindParametersWithCondition(
project_parameters["properties"],
'vorticity_induced_lift_parameters'))
if self.use_lumped_mass_matrix:
self.model_part.ProcessInfo[Kratos.COMPUTE_LUMPED_MASS_MATRIX] = 1
else:
self.model_part.ProcessInfo[Kratos.COMPUTE_LUMPED_MASS_MATRIX] = 0
self.CreateCPluPlusStrategies()
def __init__(self, model, varying_parameters=Parameters("{}")):
super().__init__(model, varying_parameters)
self._GetSolver().is_rotating_frame = self.project_parameters[
"frame_of_reference"]["frame_type"].GetInt()
self._GetDEMAnalysis().mdpas_folder_path = os.path.join(
self._GetDEMAnalysis().main_path, 'candelier_tests')
candelier_pp.include_history_force = self.vars_man.do_include_history_force
candelier_pp.include_lift = PT.RecursiveFindParametersWithCondition(
self.project_parameters["properties"],
'vorticity_induced_lift_parameters',
condition=lambda value: not (value['name'].GetString() == 'default'
))
Kratos.Logger.PrintInfo("SwimmingDEM", 'candelier_pp.include_lift ',
candelier_pp.include_lift)
candelier.sim = candelier.AnalyticSimulator(candelier_pp)
self.project_parameters["custom_fluid"][
"fluid_already_calculated"].SetBool(True)
self.project_parameters.AddEmptyValue("load_derivatives").SetBool(
False)
def ConstructListsOfVariablesForCoupling(self, parameters):
# fluid coupling variables
self.coupling_fluid_vars = []
self.coupling_fluid_vars += [Kratos.MATERIAL_ACCELERATION]
self.coupling_fluid_vars += [Fluid.MASS_SOURCE]
self.coupling_fluid_vars += [SDEM.EXACT_VELOCITY]
self.coupling_fluid_vars += [SDEM.VECTORIAL_ERROR]
self.coupling_fluid_vars += [SDEM.ERROR_X]
self.coupling_fluid_vars += [SDEM.ERROR_Y]
self.coupling_fluid_vars += [SDEM.ERROR_Z]
self.coupling_fluid_vars += [SDEM.ERROR_P]
self.coupling_fluid_vars += [SDEM.SCALAR_ERROR]
self.coupling_fluid_vars += [SDEM.EXACT_PRESSURE]
self.coupling_fluid_vars += [
Kratos.KratosGlobals.GetVariable(
parameters["body_force_per_unit_mass_variable_name"].GetString(
))
]
if parameters["custom_fluid"]["fluid_model_type"].GetInt() == 0:
self.coupling_fluid_vars += [SDEM.AVERAGED_FLUID_VELOCITY]
if (parameters["custom_fluid"]["fluid_model_type"].GetInt() == 0
or parameters["coupling"]["coupling_level_type"].GetInt() > 1):
self.coupling_fluid_vars += [Kratos.FLUID_FRACTION]
self.coupling_fluid_vars += [Kratos.FLUID_FRACTION_OLD]
if 'DISPERSE_FRACTION' in self.nodal_results:
self.coupling_fluid_vars += [Kratos.DISPERSE_FRACTION]
if parameters["coupling"]["backward_coupling"][
"filter_velocity_option"].GetBool():
self.coupling_fluid_vars += [Kratos.PARTICLE_VEL_FILTERED]
self.coupling_fluid_vars += [Kratos.TIME_AVERAGED_ARRAY_3]
self.coupling_fluid_vars += [Kratos.PHASE_FRACTION]
if parameters["custom_fluid"]["fluid_model_type"].GetInt() >= 1:
self.coupling_fluid_vars += [Kratos.FLUID_FRACTION_GRADIENT]
self.coupling_fluid_vars += [Kratos.FLUID_FRACTION_RATE]
if parameters["coupling"]["coupling_level_type"].GetInt() >= 1:
self.coupling_fluid_vars += [Kratos.HYDRODYNAMIC_REACTION]
if parameters["coupling"]["coupling_level_type"].GetInt(
) >= 1 and parameters["coupling"]["forward_coupling"][
"time_averaging_type"].GetInt() > 0:
self.coupling_fluid_vars += [Kratos.MEAN_HYDRODYNAMIC_REACTION]
if parameters["non_newtonian_fluid"]["non_newtonian_option"].GetBool():
self.coupling_fluid_vars += [Kratos.POWER_LAW_N]
self.coupling_fluid_vars += [Kratos.POWER_LAW_K]
self.coupling_fluid_vars += [Kratos.YIELD_STRESS]
# self.coupling_fluid_vars += [Kratos.GEL_STRENGTH] # TODO: make specific option for this
if parameters["coupling"]["backward_coupling"][
"viscosity_modification_type"].GetInt():
self.coupling_fluid_vars += [Kratos.VISCOSITY]
if parameters["custom_fluid"]["embedded_option"].GetBool():
self.coupling_fluid_vars += [Kratos.DISTANCE]
# dem coupling variables
self.coupling_dem_vars = []
if parameters["coupling"]["coupling_level_type"].GetInt() > 0:
self.coupling_dem_vars += [Kratos.FLUID_VEL_PROJECTED]
self.coupling_dem_vars += [Kratos.FLUID_ACCEL_PROJECTED]
self.coupling_dem_vars += [Kratos.FLUID_DENSITY_PROJECTED]
self.coupling_dem_vars += [Kratos.FLUID_VISCOSITY_PROJECTED]
self.coupling_dem_vars += [Kratos.HYDRODYNAMIC_FORCE]
self.coupling_dem_vars += [Kratos.HYDRODYNAMIC_MOMENT]
self.coupling_dem_vars += [Kratos.MATERIAL_FLUID_ACCEL_PROJECTED]
self.coupling_dem_vars += [Kratos.FLUID_ACCEL_PROJECTED]
self.coupling_dem_vars += [
Kratos.FLUID_ACCEL_FOLLOWING_PARTICLE_PROJECTED
]
self.coupling_dem_vars += [Kratos.ADDITIONAL_FORCE
] # Here for safety for the moment
if PT.RecursiveFindTrueBoolInParameters(
parameters["properties"], 'do_apply_faxen_corrections'):
self.coupling_dem_vars += [Kratos.FLUID_VEL_LAPL_PROJECTED]
self.coupling_dem_vars += [
Kratos.FLUID_VEL_LAPL_RATE_PROJECTED
]
if self.do_include_history_force:
self.coupling_dem_vars += [Kratos.FLUID_VEL_PROJECTED_RATE]
if parameters["coupling"]["coupling_level_type"].GetInt(
) >= 1 or parameters["custom_fluid"]["fluid_model_type"].GetInt() == 0:
self.coupling_dem_vars += [Kratos.FLUID_FRACTION_PROJECTED]
if (parameters["coupling"]["coupling_level_type"].GetInt() >= 1 and
'FLUID_FRACTION_GRADIENT_PROJECTED' in self.dem_printing_vars):
self.coupling_dem_vars += [
Kratos.FLUID_FRACTION_GRADIENT_PROJECTED
]
if parameters["non_newtonian_fluid"]["non_newtonian_option"].GetBool():
self.coupling_dem_vars += [Kratos.POWER_LAW_N]
self.coupling_dem_vars += [Kratos.POWER_LAW_K]
self.coupling_dem_vars += [Kratos.YIELD_STRESS]
if PT.RecursiveFindParametersWithCondition(
parameters["properties"],
'vorticity_induced_lift_parameters',
condition=lambda value: not (value['name'].GetString() ==
'default')):
self.coupling_dem_vars += [Kratos.FLUID_VORTICITY_PROJECTED]
self.coupling_dem_vars += [Kratos.SHEAR_RATE_PROJECTED]
if parameters["custom_fluid"]["embedded_option"].GetBool():
self.coupling_dem_vars += [Kratos.DISTANCE]
if 'REYNOLDS_NUMBER' in self.dem_nodal_results:
self.coupling_dem_vars += [Kratos.REYNOLDS_NUMBER]
if self.do_backward_coupling:
if parameters["coupling"]["backward_coupling"][
"apply_time_filter_to_fluid_fraction_option"].GetBool():
self.time_filtered_vars += [Kratos.FLUID_FRACTION]
if parameters["coupling"]["backward_coupling"][
"filter_velocity_option"].GetBool():
self.time_filtered_vars += [Kratos.PARTICLE_VEL_FILTERED]
if self.time_filtered_vars:
self.coupling_fluid_vars += [Kratos.TIME_AVERAGED_DOUBLE]
self.coupling_fluid_vars += [Kratos.TIME_AVERAGED_ARRAY_3]
def ConstructListsOfVariables(self, parameters):
# PRINTING VARIABLES
# constructing lists of variables to be printed
self.nodal_results, self.gauss_points_results = [], []
self.fluid_parameters = parameters['fluid_parameters']
if self.fluid_parameters.Has('sdem_output_processes'):
gid_output_options = self.fluid_parameters[
"sdem_output_processes"]["gid_output"][0]["Parameters"]
result_file_configuration = gid_output_options[
"postprocess_parameters"]["result_file_configuration"]
gauss_point_results = result_file_configuration[
"gauss_point_results"]
nodal_variables = self.fluid_parameters["sdem_output_processes"][
"gid_output"][0]["Parameters"]["postprocess_parameters"][
"result_file_configuration"]["nodal_results"]
self.nodal_results = [
nodal_variables[i].GetString()
for i in range(nodal_variables.size())
]
self.gauss_points_results = [
gauss_point_results[i].GetString()
for i in range(gauss_point_results.size())
]
self.ConstructListsOfResultsToPrint(parameters)
# COUPLING VARIABLES
# listing the variables involved in the fluid-particles coupling
if parameters["coupling"]["coupling_level_type"].GetInt():
self.ConstructListsOfVariablesForCoupling(parameters)
# VARIABLES TO ADD
# listing nodal variables to be added to the model parts (memory will be allocated for them)
# fluid variables
self.fluid_vars = []
self.fluid_vars += [Kratos.TORQUE]
self.fluid_vars += self.fluid_printing_vars
self.fluid_vars += self.coupling_fluid_vars
if parameters["pressure_grad_recovery_type"].GetInt() > 0:
self.fluid_vars += [Fluid.RECOVERED_PRESSURE_GRADIENT]
if (parameters["gradient_calculation_type"].GetInt() > 1
or parameters["pressure_grad_recovery_type"].GetInt() > 1 or
parameters["material_acceleration_calculation_type"].GetInt()
== 7 or parameters["laplacian_calculation_type"].GetInt() > 1):
self.fluid_vars += [Kratos.NODAL_WEIGHTS]
if parameters["material_acceleration_calculation_type"].GetInt():
self.fluid_vars += [Kratos.MATERIAL_ACCELERATION]
self.fluid_vars += [Kratos.VELOCITY_COMPONENT_GRADIENT]
if (parameters["material_acceleration_calculation_type"].GetInt()
== 5
or parameters["material_acceleration_calculation_type"].
GetInt() == 6):
if parameters["store_full_gradient_option"].GetBool():
self.fluid_vars += [Kratos.VELOCITY_X_GRADIENT]
self.fluid_vars += [Kratos.VELOCITY_Y_GRADIENT]
self.fluid_vars += [Kratos.VELOCITY_Z_GRADIENT]
if (parameters["vorticity_calculation_type"].GetInt() > 0
or PT.RecursiveFindParametersWithCondition(
parameters["properties"],
'vorticity_induced_lift_parameters')):
self.fluid_vars += [Kratos.VORTICITY]
if parameters["laplacian_calculation_type"].GetInt():
self.fluid_vars += [Kratos.VELOCITY_LAPLACIAN]
if PT.RecursiveFindTrueBoolInParameters(parameters["properties"],
'do_apply_faxen_corrections'):
self.fluid_vars += [Kratos.VELOCITY_LAPLACIAN_RATE]
if parameters["coupling"]["backward_coupling"][
"calculate_diffusivity_option"].GetBool():
self.fluid_vars += [Kratos.CONDUCTIVITY]
# dem variables
self.dem_vars = []
self.dem_vars += self.dem_printing_vars
self.dem_vars += self.coupling_dem_vars
self.dem_vars += [Kratos.BUOYANCY]
self.dem_vars += [Kratos.VELOCITY_OLD]
if self.do_include_history_force:
self.dem_vars += [Kratos.BASSET_FORCE]
if parameters["frame_of_reference"]["frame_type"].GetInt(
) and self.do_include_history_force > 0:
self.dem_vars += [SDEM.DISPLACEMENT_OLD]
self.dem_vars += [Kratos.VELOCITY_OLD_OLD]
if (parameters["custom_dem"]["translational_integration_scheme"].
GetString() in {'Hybrid_Bashforth', 'TerminalVelocityScheme'}
or self.do_include_history_force):
self.dem_vars += [Kratos.VELOCITY_OLD]
self.dem_vars += [Kratos.ADDITIONAL_FORCE_OLD]
self.dem_vars += [Kratos.AUX_VEL]
if parameters["add_each_hydro_force_option"].GetBool():
self.dem_vars += [Kratos.DRAG_FORCE]
self.dem_vars += [Kratos.PARTICLE_SPHERICITY
] # TODO: add only when needed
if (PT.RecursiveFindParametersWithCondition(
parameters["properties"], 'vorticity_induced_lift_parameters')
and parameters["add_each_hydro_force_option"].GetBool()):
self.dem_vars += [Kratos.LIFT_FORCE]
if parameters["add_each_hydro_force_option"].GetBool():
self.dem_vars += [Kratos.VIRTUAL_MASS_FORCE]
# clusters variables
self.clusters_vars = []
# rigid faces variables
self.rigid_faces_vars = [
Kratos.VELOCITY, Kratos.ANGULAR_VELOCITY, Kratos.DISPLACEMENT,
DEM.DELTA_DISPLACEMENT, Kratos.DELTA_ROTATION, DEM.CONTACT_FORCES,
DEM.DEM_PRESSURE, DEM.ELASTIC_FORCES, Kratos.PRESSURE,
DEM.TANGENTIAL_ELASTIC_FORCES, DEM.SHEAR_STRESS, Kratos.NODAL_AREA,
Kratos.VELOCITY_OLD
]
if parameters["custom_fluid"]["embedded_option"].GetBool():
self.rigid_faces_vars += [Kratos.FORCE]
self.rigid_faces_vars += [Kratos.POSITIVE_FACE_PRESSURE]
self.rigid_faces_vars += [Kratos.NEGATIVE_FACE_PRESSURE]
self.fluid_vars += self.rigid_faces_vars
# inlet variables
self.inlet_vars = self.dem_vars
