diff --git a/modules/navier_stokes/src/physics/WCNSFVFlowPhysics.C b/modules/navier_stokes/src/physics/WCNSFVFlowPhysics.C index 307147cb578e..d0dbc25920e3 100644 --- a/modules/navier_stokes/src/physics/WCNSFVFlowPhysics.C +++ b/modules/navier_stokes/src/physics/WCNSFVFlowPhysics.C @@ -34,8 +34,17 @@ WCNSFVFlowPhysics::validParams() true, "Whether the time derivative term should contribute to the Rhie Chow coefficients. This adds " "stabilization, but makes the solution dependent on the time step size"); - params.addParamNamesToGroup("time_derivative_contributes_to_RC_coefficients characteristic_speed", - "Numerical scheme"); + params.addParam( + "time_derivative_rc_coef_fixed_dt", + "Fixed dt used to make Rhie Chow coefficients independent of the time step before restarts"); + params.addParam("time_derivative_rc_coef_fixed_dt_start", + "Time after which to start leveraging the fixed dt to make Rhie Chow " + "coefficients independent of the time step before restarts"); + params.addParamNamesToGroup("characteristic_speed", "Numerical scheme"); + params.addParamNamesToGroup("time_derivative_contributes_to_RC_coefficients " + "time_derivative_rc_coef_fixed_dt " + "time_derivative_rc_coef_fixed_dt_start", + "Restart"); // Used for flow mixtures, where one phase is solid / not moving under the action of gravity params.addParam( @@ -375,6 +384,12 @@ WCNSFVFlowPhysics::addINSMomentumTimeKernels() params.set("rhie_chow_user_object") = rhieChowUOName(); params.set("contribute_to_rc") = getParam("time_derivative_contributes_to_RC_coefficients"); + if (isParamValid("time_derivative_rc_coef_fixed_dt")) + params.set("fixed_dt_contribution_to_RC") = + getParam("time_derivative_rc_coef_fixed_dt"); + if (isParamValid("time_derivative_rc_coef_fixed_dt_start")) + params.set("start_fixed_dt_contribution_to_RC") = + getParam("time_derivative_rc_coef_fixed_dt_start"); for (const auto d : make_range(dimension())) {