void NavierStokesExplicit::on_initial_conditions_set(InitialConditions& initial_conditions)
{
// Initial condition for the viscosity, defaulting to the molecular viscosity
Handle<ProtoAction> visc_ic(initial_conditions.create_initial_condition("navier_stokes_viscosity", "cf3.solver.ProtoAction"));
visc_ic->set_expression(nodes_expression(nu_eff = nu));
m_viscosity_initial_condition = visc_ic;
m_velocity_initial_condition = initial_conditions.create_initial_condition("navier_stokes_u_solution");
m_pressure_initial_condition = initial_conditions.create_initial_condition("navier_stokes_p_solution");
m_iteration_initial_condition = initial_conditions.create_initial_condition("navier_stokes_explicit_iteration");
}
void NavierStokes::on_initial_conditions_set(InitialConditions& initial_conditions)
{
// Initial condition for the viscosity, defaulting to the molecular viscosity
Handle<ProtoAction> visc_ic(initial_conditions.create_initial_condition("navier_stokes_viscosity", "cf3.solver.ProtoAction"));
visc_ic->set_expression(nodes_expression(nu_eff = nu));
m_initial_conditions = initial_conditions.create_initial_condition(solution_tag());
// Use a proto action to set the linearized_velocity easily
Handle<ProtoAction> lin_vel_ic (initial_conditions.create_initial_condition("linearized_velocity", "cf3.solver.ProtoAction"));
lin_vel_ic->set_expression(nodes_expression(group(u_adv = u, u1 = u, u2 = u, u3 = u)));
}
void NavierStokesManual::on_initial_conditions_set(InitialConditions& initial_conditions)
{
// Use proto to add an initial condition on the viscosity, setting it equal to the molecular viscosity
FieldVariable<0, ScalarField> nu_eff("EffectiveViscosity", "navier_stokes_viscosity");
PhysicsConstant nu("kinematic_viscosity");
Handle<ProtoAction> visc_ic(initial_conditions.create_initial_condition("navier_stokes_viscosity", "cf3.solver.ProtoAction"));
visc_ic->set_expression(nodes_expression(nu_eff = nu));
}
void ScalarAdvection::on_initial_conditions_set ( InitialConditions& initial_conditions )
{
initial_conditions.create_initial_condition(solution_tag());
}
void HeatConductionSteady::on_initial_conditions_set(InitialConditions& initial_conditions)
{
initial_conditions.create_initial_condition(solution_tag());
}
