INSFVMushyPorousFrictionFunctorMaterial

Computes the mushy zone porous resistance for solidification/melting problems.

This material computes the Darcy_coefficient and the Forchheimer_coefficient for a solidification problem. The model uses a mushy-zone approach to compute the friction coefficients.

The coefficients are defined as follows:

D=μK,D = \frac{\mu}{K} \,,

where the permeability KK is defined as follows:

K=fl3fs2Fkcs,K = \frac{f_l^3}{f_s^2 F_k c_s} \,,

where flf_l and fsf_s are the liquid and solid fraction, the switching function is defined as Fk=0.5atan(s(fsfs,crit)/π)F_k = 0.5 * \operatorname{atan}(s (f_s - f_{s,crit}) / \pi) with s=100s=100 and fs,crit=0.27f_{s,crit}=0.27 and cs=c/(λ)2c_s = c / (\lambda)^2 with c=180c = 180 and the dendrite spacing is λ\lambda. By default, we set λ=104\lambda = 10^{-4}. However, this coefficient may be tuned by the user if experimental data is available.

The Forchheimer coefficient is defined as follows:

F=CFρlK,F = C_F \frac{\rho_l}{\sqrt{K}} \,,

where the Ergun coefficient is CF=0.55C_F = 0.55 and ρl\rho_l is the density of the liquid phase.

commentnote

This material is compatible with PINSFVMomentumFriction for the non standard (simplified) formulation which multiplies the incoming Forchheimer coefficient by the velocity magnitude times velocity component; it is incompatible with the standard Forchheimer formulation which also multiplies with density over two. Thus the user should set the parameter "standard_friction_formulation = false"

Input Parameters

  • liquid_fractionLiquid Fraction Functor. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:Liquid Fraction Functor. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • muThe liquid dynamic viscosity. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:The liquid dynamic viscosity. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • rho_lThe liquid density (not the mixture one). A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:The liquid density (not the mixture one). A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

Required Parameters

  • Darcy_coef_nameDarcy_coefficientName of the Darcy friction coefficient. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    Default:Darcy_coefficient

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:Name of the Darcy friction coefficient. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • Forchheimer_coef_nameForchheimer_coefficientName of the Forchheimer friction coefficient. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    Default:Forchheimer_coefficient

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:Name of the Forchheimer friction coefficient. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • blockThe list of blocks (ids or names) that this object will be applied

    C++ Type:std::vector<SubdomainName>

    Unit:(no unit assumed)

    Controllable:No

    Description:The list of blocks (ids or names) that this object will be applied

  • boundaryThe list of boundaries (ids or names) from the mesh where this object applies

    C++ Type:std::vector<BoundaryName>

    Unit:(no unit assumed)

    Controllable:No

    Description:The list of boundaries (ids or names) from the mesh where this object applies

  • constant_onNONEWhen ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped

    Default:NONE

    C++ Type:MooseEnum

    Unit:(no unit assumed)

    Options:NONE, ELEMENT, SUBDOMAIN

    Controllable:No

    Description:When ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped

  • declare_suffixAn optional suffix parameter that can be appended to any declared properties. The suffix will be prepended with a '_' character.

    C++ Type:MaterialPropertyName

    Unit:(no unit assumed)

    Controllable:No

    Description:An optional suffix parameter that can be appended to any declared properties. The suffix will be prepended with a '_' character.

  • dendrite_spacing_scaling1e-4The dendrite spacing scaling. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    Default:1e-4

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:The dendrite spacing scaling. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • execute_onALWAYSThe list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.

    Default:ALWAYS

    C++ Type:ExecFlagEnum

    Unit:(no unit assumed)

    Options:NONE, INITIAL, LINEAR, NONLINEAR_CONVERGENCE, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM, ALWAYS

    Controllable:No

    Description:The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.

  • prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.

    C++ Type:MaterialPropertyName

    Unit:(no unit assumed)

    Controllable:No

    Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.

  • use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

    Default:False

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

Optional Parameters

  • control_tagsAdds user-defined labels for accessing object parameters via control logic.

    C++ Type:std::vector<std::string>

    Unit:(no unit assumed)

    Controllable:No

    Description:Adds user-defined labels for accessing object parameters via control logic.

  • enableTrueSet the enabled status of the MooseObject.

    Default:True

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:Yes

    Description:Set the enabled status of the MooseObject.

  • implicitTrueDetermines whether this object is calculated using an implicit or explicit form

    Default:True

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:Determines whether this object is calculated using an implicit or explicit form

  • seed0The seed for the master random number generator

    Default:0

    C++ Type:unsigned int

    Unit:(no unit assumed)

    Controllable:No

    Description:The seed for the master random number generator

Advanced Parameters

  • output_propertiesList of material properties, from this material, to output (outputs must also be defined to an output type)

    C++ Type:std::vector<std::string>

    Unit:(no unit assumed)

    Controllable:No

    Description:List of material properties, from this material, to output (outputs must also be defined to an output type)

  • outputsnone Vector of output names where you would like to restrict the output of variables(s) associated with this object

    Default:none

    C++ Type:std::vector<OutputName>

    Unit:(no unit assumed)

    Controllable:No

    Description:Vector of output names where you would like to restrict the output of variables(s) associated with this object

Outputs Parameters