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:

$var element = document.getElementById("moose-equation-7fc9117b-38ca-4c9c-9bbe-ede9f895008e");katex.render("A^{-1}", element, {displayMode:false,throwOnError:false});$ (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.
$var element = document.getElementById("moose-equation-6c23d36e-0a69-43a8-a3a5-8c2b35dcd5fd");katex.render("A^{-1}H(u)", element, {displayMode:false,throwOnError:false});$ 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.
$var element = document.getElementById("moose-equation-68133a6d-69ca-4a41-90e4-624edbb8ff60");katex.render("(\\rho \\vec{u} \\vec{n})_{RC}", element, {displayMode:false,throwOnError:false});$ 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 ( $var element = document.getElementById("moose-equation-ab2966f1-e4bd-4c05-9691-579081638633");katex.render("A^{-1}", element, {displayMode:false,throwOnError:false});$ ) and neighbour face flux ( $var element = document.getElementById("moose-equation-89344cc9-dfaf-4c73-b48f-6423d8528aea");katex.render("H(u)", element, {displayMode:false,throwOnError:false});$ ) 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.

Reconstructed Pressure Gradient

The reconstructed pressure-gradient option uses the Aguerre reconstruction (Aguerre et al. (2018)) implemented by FVReconstructedPressureGradient. It is useful when the conservative Rhie-Chow face flux is smooth, but the cell-centered velocity still exhibits oscillations because its pressure gradient is not fully consistent with that face flux.

After a pressure correction, the method reconstructs a cell velocity that is compatible with the corrected conservative face flux and then determines the pressure gradient required by the cell momentum balance. The newly reconstructed gradient corrects the velocity immediately so that the cell velocity and continuity-preserving face flux remain consistent. A relaxed version of that gradient is used by the next momentum predictor to avoid introducing an abrupt change into the pressure-velocity coupling. The accepted gradient is carried between time steps and is preserved through time-step retries and restarts.

The following options control the reconstruction:

"gradient_relaxation" controls how strongly the newly reconstructed gradient changes the gradient used by the next momentum solve. The default value of 0.1 is deliberately conservative. Smaller values provide more damping and can improve robustness, but usually require more outer iterations. Values closer to 1 respond more quickly to the latest pressure correction, but can strengthen pressure-velocity oscillations.
"base_gradient_method" selects the ordinary pressure-gradient method used before the first reconstructed gradient is available. The default is green-gauss.

For example, the following tested input selects the reconstructed method through an input-file variable and uses that selection for the pressure variable:

    gradient_method = ${pressure_gradient_method}

[Variables<<<{"href": "../../syntax/Variables/index.html"}>>>]
  [pressure]
    type = MooseLinearVariableFVReal<<<{"description": "Base class for Moose variables. This should never be the terminal object type", "href": "../variables/MooseLinearVariableFV.html"}>>>
    initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = 0
    solver_sys<<<{"description": "If this variable is a solver variable, this is the solver system to which it should be added."}>>> = pressure_system
    gradient_method<<<{"description": "Default gradient computation method to register when a consumer requests gradients from this variable. This may be a built-in method name like 'green-gauss' or 'green-gauss-venkatakrishnan', or the name of an object in [FVGradientMethods]."}>>> = ${pressure_gradient_method}
  []
[]
[FVGradientMethods<<<{"href": "../../syntax/FVGradientMethods/index.html"}>>>]
  [reconstructed]
    type = FVReconstructedPressureGradient<<<{"description": "Reconstructs and relaxes the pressure gradient used for Rhie-Chow momentum coupling.", "href": "../fvgradientmethods/FVReconstructedPressureGradient.html"}>>>
    base_gradient_method<<<{"description": "Gradient method used before Rhie-Chow has computed reconstructed gradients."}>>> = green-gauss
    # This is the relaxation used by the reconstructed-gradient channel regression tests.
    gradient_relaxation<<<{"description": "Relaxation factor applied when updating the reconstructed pressure-coupling gradient."}>>> = 0.1
  []
[]

Use the same reconstructed pressure-gradient definition for every LinearFVMomentumPressure component coupled to this pressure equation. Mixing gradient definitions among velocity components would make the reconstructed cell velocity inconsistent with the coupled momentum balance. Other equations and diagnostic quantities should continue to use an ordinary gradient method.

Pressure Diffusion Interpolation

The "pressure_diffusion_interpolation" parameter selects whether average or harmonic interpolation is used when computing the face values of Ainv, the $var element = document.getElementById("moose-equation-e397a02c-a2ca-4310-8136-e9df219a9afb");katex.render("A^{-1}", element, {displayMode:false,throwOnError:false});$ diffusion tensor in the pressure correction diffusion term. When using Navier Stokes Flow Segregated / WCNSLinearFVFlowPhysics, select this through the pressure_diffusion_interpolation parameter in the Physics block.

Input Parameters

p_diffusion_kernelThe LinearFVPressureCorrectionDiffusion kernel acting on the pressure.
C++ Type:std::string
Controllable:No
Description:The LinearFVPressureCorrectionDiffusion 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_diffusion_interpolationaverageThe face interpolation method for Ainv in the pressure correction diffusion term.
Default:average
C++ Type:MooseEnum
Options:average, harmonic
Controllable:No
Description:The face interpolation method for Ainv in the pressure correction diffusion term.
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

References

Horacio J Aguerre, Cesar I Pairetti, Cesar M Venier, Santiago Márquez Damián, and Norberto M Nigro. An oscillation-free flow solver based on flux reconstruction. Journal of Computational Physics, 365:135–148, 2018.

BibTeX

@article{aguerre2018oscillation,
    author = "Aguerre, Horacio J and Pairetti, Cesar I and Venier, Cesar M and Dami{\'a}n, Santiago M{\'a}rquez and Nigro, Norberto M",
    title = "An oscillation-free flow solver based on flux reconstruction",
    journal = "Journal of Computational Physics",
    volume = "365",
    pages = "135--148",
    year = "2018",
    publisher = "Elsevier"
}

RIS

TY  - JOUR
AU  - Aguerre, Horacio J
AU  - Pairetti, Cesar I
AU  - Venier, Cesar M
AU  - Damián, Santiago Márquez
AU  - Nigro, Norberto M
TI  - An oscillation-free flow solver based on flux reconstruction
JO  - Journal of Computational Physics
PY  - 2018
VL  - 365
PB  - Elsevier
SP  - 135
EP  - 148
ER  -

Plain Text

Horacio J Aguerre, Cesar I Pairetti, Cesar M Venier, Santiago Márquez Damián, and Norberto M Nigro. An oscillation-free flow solver based on flux reconstruction. Journal of Computational Physics, 365:135–148, 2018.

(moose/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/linear-segregated/2d/reconstructed-force-channel.i)

mu = 1
rho = 1
pressure_gradient_method = reconstructed

[Mesh]
  [channel]
    type = CartesianMeshGenerator
    dim = 2
    # Twenty cells per section resolve the interface-localized oscillation without making this
    # two-solve comparison unnecessarily expensive.
    dx = '1 1 1'
    ix = '20 20 20'
    dy = 1
    iy = 3
    subdomain_id = '0 1 2'
  []
  [blocks]
    type = RenameBlockGenerator
    input = channel
    old_block = '0 1 2'
    new_block = 'upstream forced downstream'
  []
[]

[Problem]
  linear_sys_names = 'u_system v_system pressure_system'
  previous_nl_solution_required = true
[]

[UserObjects]
  [rc]
    type = RhieChowMassFlux
    u = vel_x
    v = vel_y
    pressure = pressure
    rho = ${rho}
    p_diffusion_kernel = p_diffusion
  []
[]

[Variables]
  [vel_x]
    type = MooseLinearVariableFVReal
    initial_condition = 1
    solver_sys = u_system
  []
  [vel_y]
    type = MooseLinearVariableFVReal
    initial_condition = 0
    solver_sys = v_system
  []
  [pressure]
    type = MooseLinearVariableFVReal
    initial_condition = 0
    solver_sys = pressure_system
    gradient_method = ${pressure_gradient_method}
  []
[]

[FVGradientMethods]
  [reconstructed]
    type = FVReconstructedPressureGradient
    base_gradient_method = green-gauss
    # This is the relaxation used by the reconstructed-gradient channel regression tests.
    gradient_relaxation = 0.1
  []
[]

[FVInterpolationMethods]
  [upwind]
    type = FVAdvectedUpwind
  []
[]

[LinearFVKernels]
  [u_advection_stress]
    type = LinearWCNSFVMomentumFlux
    variable = vel_x
    advected_interp_method_name = upwind
    mu = ${mu}
    u = vel_x
    v = vel_y
    momentum_component = x
    rhie_chow_user_object = rc
    use_nonorthogonal_correction = false
  []
  [v_advection_stress]
    type = LinearWCNSFVMomentumFlux
    variable = vel_y
    advected_interp_method_name = upwind
    mu = ${mu}
    u = vel_x
    v = vel_y
    momentum_component = y
    rhie_chow_user_object = rc
    use_nonorthogonal_correction = false
  []
  [u_pressure]
    type = LinearFVMomentumPressure
    variable = vel_x
    pressure = pressure
    momentum_component = x
  []
  [v_pressure]
    type = LinearFVMomentumPressure
    variable = vel_y
    pressure = pressure
    momentum_component = y
  []
  [middle_force]
    type = LinearFVSource
    variable = vel_x
    # The large source makes the pressure-gradient jump dominate truncation error and visibly
    # excites the unsupported cell-velocity mode when reconstruction is disabled.
    source_density = 1e4
    block = forced
  []
  [p_diffusion]
    type = LinearFVPressureCorrectionDiffusion
    variable = pressure
    diffusion_tensor = Ainv
    use_nonorthogonal_correction = false
  []
  [HbyA_divergence]
    type = LinearFVDivergence
    variable = pressure
    face_flux = HbyA
    force_boundary_execution = true
  []
[]

[LinearFVBCs]
  [inlet_u]
    type = LinearFVAdvectionDiffusionFunctorDirichletBC
    boundary = left
    variable = vel_x
    functor = 1
  []
  [inlet_v]
    type = LinearFVAdvectionDiffusionFunctorDirichletBC
    boundary = left
    variable = vel_y
    functor = 0
  []
  [symmetry_u]
    type = LinearFVVelocitySymmetryBC
    boundary = 'top bottom'
    variable = vel_x
    momentum_component = x
    u = vel_x
    v = vel_y
  []
  [symmetry_v]
    type = LinearFVVelocitySymmetryBC
    boundary = 'top bottom'
    variable = vel_y
    momentum_component = y
    u = vel_x
    v = vel_y
  []
  [pressure_flux]
    type = LinearFVPressureFluxBC
    boundary = left
    variable = pressure
    HbyA_flux = HbyA
    Ainv = Ainv
    u = vel_x
    v = vel_y
    rho = ${rho}
  []
  [pressure_symmetry]
    type = LinearFVPressureSymmetryBC
    boundary = 'top bottom'
    variable = pressure
    HbyA_flux = HbyA
  []
  [outlet_u]
    type = LinearFVAdvectionDiffusionOutflowBC
    boundary = right
    variable = vel_x
    use_two_term_expansion = false
  []
  [outlet_v]
    type = LinearFVAdvectionDiffusionOutflowBC
    boundary = right
    variable = vel_y
    use_two_term_expansion = false
  []
  [outlet_p]
    type = LinearFVAdvectionDiffusionFunctorDirichletBC
    boundary = right
    variable = pressure
    functor = 0
  []
[]

[Postprocessors]
  [upstream_vel_x]
    type = ElementAverageValue
    variable = vel_x
    block = upstream
    execute_on = FINAL
  []
  [forced_vel_x]
    type = ElementAverageValue
    variable = vel_x
    block = forced
    execute_on = FINAL
  []
  [downstream_vel_x]
    type = ElementAverageValue
    variable = vel_x
    block = downstream
    execute_on = FINAL
  []
[]

[Executioner]
  type = SIMPLE
  momentum_l_abs_tol = 1e-10
  pressure_l_abs_tol = 1e-10
  momentum_l_tol = 0
  pressure_l_tol = 0
  rhie_chow_user_object = rc
  momentum_systems = 'u_system v_system'
  pressure_system = pressure_system
  momentum_equation_relaxation = 0.8
  pressure_variable_relaxation = 0.3
  # Five hundred iterations allow the reconstructed coupling pressure gradient to relax to the same
  # uniform steady field as the standard gradient. The reconstructed momentum residual is
  # cancellation-dominated at the source jumps and does not satisfy the generic normalized monitor.
  num_iterations = 500
  continue_on_max_its = true
  # Relative to the 1e4 source, this requires the ordinary outer residuals to fall by eleven orders
  # of magnitude while remaining above platform-level roundoff stagnation.
  pressure_absolute_tolerance = 1e-7
  momentum_absolute_tolerance = 1e-7
  momentum_petsc_options_iname = '-pc_type'
  momentum_petsc_options_value = lu
  pressure_petsc_options_iname = '-pc_type'
  pressure_petsc_options_value = lu
  print_fields = false
[]

[Outputs]
  [csv]
    type = CSV
    execute_on = FINAL
  []
  exodus = true
[]

(moose/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/linear-segregated/2d/reconstructed-force-channel.i)

mu = 1
rho = 1
pressure_gradient_method = reconstructed

[Mesh]
  [channel]
    type = CartesianMeshGenerator
    dim = 2
    # Twenty cells per section resolve the interface-localized oscillation without making this
    # two-solve comparison unnecessarily expensive.
    dx = '1 1 1'
    ix = '20 20 20'
    dy = 1
    iy = 3
    subdomain_id = '0 1 2'
  []
  [blocks]
    type = RenameBlockGenerator
    input = channel
    old_block = '0 1 2'
    new_block = 'upstream forced downstream'
  []
[]

[Problem]
  linear_sys_names = 'u_system v_system pressure_system'
  previous_nl_solution_required = true
[]

[UserObjects]
  [rc]
    type = RhieChowMassFlux
    u = vel_x
    v = vel_y
    pressure = pressure
    rho = ${rho}
    p_diffusion_kernel = p_diffusion
  []
[]

[Variables]
  [vel_x]
    type = MooseLinearVariableFVReal
    initial_condition = 1
    solver_sys = u_system
  []
  [vel_y]
    type = MooseLinearVariableFVReal
    initial_condition = 0
    solver_sys = v_system
  []
  [pressure]
    type = MooseLinearVariableFVReal
    initial_condition = 0
    solver_sys = pressure_system
    gradient_method = ${pressure_gradient_method}
  []
[]

[FVGradientMethods]
  [reconstructed]
    type = FVReconstructedPressureGradient
    base_gradient_method = green-gauss
    # This is the relaxation used by the reconstructed-gradient channel regression tests.
    gradient_relaxation = 0.1
  []
[]

[FVInterpolationMethods]
  [upwind]
    type = FVAdvectedUpwind
  []
[]

[LinearFVKernels]
  [u_advection_stress]
    type = LinearWCNSFVMomentumFlux
    variable = vel_x
    advected_interp_method_name = upwind
    mu = ${mu}
    u = vel_x
    v = vel_y
    momentum_component = x
    rhie_chow_user_object = rc
    use_nonorthogonal_correction = false
  []
  [v_advection_stress]
    type = LinearWCNSFVMomentumFlux
    variable = vel_y
    advected_interp_method_name = upwind
    mu = ${mu}
    u = vel_x
    v = vel_y
    momentum_component = y
    rhie_chow_user_object = rc
    use_nonorthogonal_correction = false
  []
  [u_pressure]
    type = LinearFVMomentumPressure
    variable = vel_x
    pressure = pressure
    momentum_component = x
  []
  [v_pressure]
    type = LinearFVMomentumPressure
    variable = vel_y
    pressure = pressure
    momentum_component = y
  []
  [middle_force]
    type = LinearFVSource
    variable = vel_x
    # The large source makes the pressure-gradient jump dominate truncation error and visibly
    # excites the unsupported cell-velocity mode when reconstruction is disabled.
    source_density = 1e4
    block = forced
  []
  [p_diffusion]
    type = LinearFVPressureCorrectionDiffusion
    variable = pressure
    diffusion_tensor = Ainv
    use_nonorthogonal_correction = false
  []
  [HbyA_divergence]
    type = LinearFVDivergence
    variable = pressure
    face_flux = HbyA
    force_boundary_execution = true
  []
[]

[LinearFVBCs]
  [inlet_u]
    type = LinearFVAdvectionDiffusionFunctorDirichletBC
    boundary = left
    variable = vel_x
    functor = 1
  []
  [inlet_v]
    type = LinearFVAdvectionDiffusionFunctorDirichletBC
    boundary = left
    variable = vel_y
    functor = 0
  []
  [symmetry_u]
    type = LinearFVVelocitySymmetryBC
    boundary = 'top bottom'
    variable = vel_x
    momentum_component = x
    u = vel_x
    v = vel_y
  []
  [symmetry_v]
    type = LinearFVVelocitySymmetryBC
    boundary = 'top bottom'
    variable = vel_y
    momentum_component = y
    u = vel_x
    v = vel_y
  []
  [pressure_flux]
    type = LinearFVPressureFluxBC
    boundary = left
    variable = pressure
    HbyA_flux = HbyA
    Ainv = Ainv
    u = vel_x
    v = vel_y
    rho = ${rho}
  []
  [pressure_symmetry]
    type = LinearFVPressureSymmetryBC
    boundary = 'top bottom'
    variable = pressure
    HbyA_flux = HbyA
  []
  [outlet_u]
    type = LinearFVAdvectionDiffusionOutflowBC
    boundary = right
    variable = vel_x
    use_two_term_expansion = false
  []
  [outlet_v]
    type = LinearFVAdvectionDiffusionOutflowBC
    boundary = right
    variable = vel_y
    use_two_term_expansion = false
  []
  [outlet_p]
    type = LinearFVAdvectionDiffusionFunctorDirichletBC
    boundary = right
    variable = pressure
    functor = 0
  []
[]

[Postprocessors]
  [upstream_vel_x]
    type = ElementAverageValue
    variable = vel_x
    block = upstream
    execute_on = FINAL
  []
  [forced_vel_x]
    type = ElementAverageValue
    variable = vel_x
    block = forced
    execute_on = FINAL
  []
  [downstream_vel_x]
    type = ElementAverageValue
    variable = vel_x
    block = downstream
    execute_on = FINAL
  []
[]

[Executioner]
  type = SIMPLE
  momentum_l_abs_tol = 1e-10
  pressure_l_abs_tol = 1e-10
  momentum_l_tol = 0
  pressure_l_tol = 0
  rhie_chow_user_object = rc
  momentum_systems = 'u_system v_system'
  pressure_system = pressure_system
  momentum_equation_relaxation = 0.8
  pressure_variable_relaxation = 0.3
  # Five hundred iterations allow the reconstructed coupling pressure gradient to relax to the same
  # uniform steady field as the standard gradient. The reconstructed momentum residual is
  # cancellation-dominated at the source jumps and does not satisfy the generic normalized monitor.
  num_iterations = 500
  continue_on_max_its = true
  # Relative to the 1e4 source, this requires the ordinary outer residuals to fall by eleven orders
  # of magnitude while remaining above platform-level roundoff stagnation.
  pressure_absolute_tolerance = 1e-7
  momentum_absolute_tolerance = 1e-7
  momentum_petsc_options_iname = '-pc_type'
  momentum_petsc_options_value = lu
  pressure_petsc_options_iname = '-pc_type'
  pressure_petsc_options_value = lu
  print_fields = false
[]

[Outputs]
  [csv]
    type = CSV
    execute_on = FINAL
  []
  exodus = true
[]

Overview
Input Parameters
References