def __init__(self, model_part, settings):
super().__init__(model_part, settings)
default_settings = Kratos.Parameters(r'''
{
"formulation_name": "fractional_step_k_epsilon",
"incompressible_potential_flow_initialization_settings": {},
"fractional_step_flow_solver_settings": {},
"k_epsilon_solver_settings": {},
"max_iterations": 1
}''')
settings.ValidateAndAssignDefaults(default_settings)
if (not settings[
"incompressible_potential_flow_initialization_settings"].
IsEquivalentTo(Kratos.Parameters("{}"))):
self.incompressible_potential_flow_formulation = IncompressiblePotentialFlowRansFormulation(
model_part, settings[
"incompressible_potential_flow_initialization_settings"])
self.AddRansFormulation(
self.incompressible_potential_flow_formulation)
self.fractional_step_formulation = FractionalStepVelocityPressureRansFormulation(
model_part, settings["fractional_step_flow_solver_settings"])
self.AddRansFormulation(self.fractional_step_formulation)
self.k_epsilon_formulation = KEpsilonRansFormulation(
model_part, settings["k_epsilon_solver_settings"])
self.AddRansFormulation(self.k_epsilon_formulation)
self.SetMaxCouplingIterations(settings["max_iterations"].GetInt())
class MonolithicKEpsilonRansFormulation(RansFormulation):
def __init__(self, model_part, settings):
super().__init__(model_part, settings)
default_settings = Kratos.Parameters(r'''
{
"formulation_name": "monolithic_k_epsilon",
"incompressible_potential_flow_initialization_settings": {},
"monolithic_flow_solver_settings": {},
"k_epsilon_solver_settings": {},
"max_iterations": 1
}''')
settings.ValidateAndAssignDefaults(default_settings)
if (not settings[
"incompressible_potential_flow_initialization_settings"].
IsEquivalentTo(Kratos.Parameters("{}"))):
self.incompressible_potential_flow_formulation = IncompressiblePotentialFlowRansFormulation(
model_part, settings[
"incompressible_potential_flow_initialization_settings"])
self.AddFormulation(self.incompressible_potential_flow_formulation)
self.monolithic_formulation = MonolithicVelocityPressureRansFormulation(
model_part, settings["monolithic_flow_solver_settings"])
self.AddRansFormulation(self.monolithic_formulation)
self.k_epsilon_formulation = KEpsilonRansFormulation(
model_part, settings["k_epsilon_solver_settings"])
self.AddRansFormulation(self.k_epsilon_formulation)
self.SetMaxCouplingIterations(settings["max_iterations"].GetInt())
def SetConstants(self, settings):
self.k_epsilon_formulation.SetConstants(settings)
def Initialize(self):
super().Initialize()
if (self.monolithic_formulation.ElementHasNodalProperties()):
nut_nodal_update_process = KratosRANS.RansNutNodalUpdateProcess(
self.GetBaseModelPart().GetModel(),
self.GetBaseModelPart().Name,
self.k_epsilon_formulation.echo_level)
self.k_epsilon_formulation.AddProcess(nut_nodal_update_process)