RhieChowMassFlux

Computes H/A and 1/A together with face mass fluxes for segregated momentum-pressure equations using linear systems.

Overview

This object is responsible for generating the following fields for a SIMPLE-type segregated solver:

  • (inverse of the matrix diagonal) which is used as a diffusivity for the pressure equation. This field is stored in a face-based functor, so face values are easy to access but cell values need to be reconstructed. This is mainly used in the pressure Poisson equation where only face values are queried.

  • whose divergence is used as a source in the pressure Poisson equation. This field is also stored in a face-based functor, so face values are easy to access, but cell-center values need to be reconstructed.

  • which is the Rhie-Chow corrected face mass flux. This is also stored in a face-based functor, so face values are easy to access, but cell-center values need to be reconstructed.

Besides these capabilities, this user object is also responsible for reconstructing cell velocities at the end of the pressure corrector step. For more information on these fields and processes, we suggest visiting SIMPLE.

The object enables the computation of the standard (SIMPLE) or consistent (SIMPLEC) momentum projection matrix () and neighbour face flux () vector via "pressure_projection_method". In general, SIMPLEC will be stable with higher relaxation factors for pressure than SIMPLE. This is particularly useful in problems with slow-converging pressure fields, such as those with high Reynolds numbers, complex geometries, viscous flows in narrow channels, multiphase flows, problems with rapidly varying thermophysical properties, and, in general, when using high-resolution grids.

Input Parameters

  • p_diffusion_kernelThe diffusion kernel acting on the pressure.

    C++ Type:std::string

    Controllable:No

    Description:The diffusion kernel acting on the pressure.

  • pressureThe pressure variable.

    C++ Type:VariableName

    Unit:(no unit assumed)

    Controllable:No

    Description:The pressure variable.

  • rhoDensity functor. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:Density functor. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.

  • uThe x-component of velocity

    C++ Type:VariableName

    Unit:(no unit assumed)

    Controllable:No

    Description:The x-component of velocity

Required Parameters

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

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

    Controllable:No

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

  • pressure_projection_methodstandardThe method to use in the pressure projection for Ainv - standard (SIMPLE) or consistent (SIMPLEC)

    Default:standard

    C++ Type:MooseEnum

    Options:standard, consistent

    Controllable:No

    Description:The method to use in the pressure projection for Ainv - standard (SIMPLE) or consistent (SIMPLEC)

  • vThe y-component of velocity

    C++ Type:VariableName

    Unit:(no unit assumed)

    Controllable:No

    Description:The y-component of velocity

  • wThe z-component of velocity

    C++ Type:VariableName

    Unit:(no unit assumed)

    Controllable:No

    Description:The z-component of velocity

Optional Parameters

  • allow_duplicate_execution_on_initialFalseIn the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).

    Default:False

    C++ Type:bool

    Controllable:No

    Description:In the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).

  • execution_order_group0Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.

    Default:0

    C++ Type:int

    Controllable:No

    Description:Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.

  • force_postauxFalseForces the UserObject to be executed in POSTAUX

    Default:False

    C++ Type:bool

    Controllable:No

    Description:Forces the UserObject to be executed in POSTAUX

  • force_preauxFalseForces the UserObject to be executed in PREAUX

    Default:False

    C++ Type:bool

    Controllable:No

    Description:Forces the UserObject to be executed in PREAUX

  • force_preicFalseForces the UserObject to be executed in PREIC during initial setup

    Default:False

    C++ Type:bool

    Controllable:No

    Description:Forces the UserObject to be executed in PREIC during initial setup

Execution Scheduling Parameters

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

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

    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

    Controllable:Yes

    Description:Set the enabled status of the MooseObject.

  • use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

    Default:False

    C++ Type:bool

    Controllable:No

    Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

  • 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

    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.

Material Property Retrieval Parameters

Input Files