LinearFVAdvectionDiffusionOutflowBC

Description

LinearFVAdvectionDiffusionOutflowBC will contribute to the system matrix and right hand side of a linear finite volume system. The contributions can be derived using the integral of the advective flux over a boundary face ( $var element = document.getElementById("moose-equation-c7b47f4a-29ec-4aba-bc60-f824ea427c78");katex.render("S_b", element, {displayMode:false,throwOnError:false});$ ) of a boundary element:

var element = document.getElementById("moose-equation-73e3b920-2ec2-4524-88a0-86d60063eab7");katex.render("\\int\\limits_{S_b} \\vec{v} \\cdot \\vec{n} u dS \\approx \\vec{v}_f \\cdot \\vec{n} u_f |S_b|~,", element, {displayMode:true,throwOnError:false});

where $var element = document.getElementById("moose-equation-b56dd10e-cbd8-43df-8eea-a93e8fa03529");katex.render("\\vec{v}_f", element, {displayMode:false,throwOnError:false});$ , $var element = document.getElementById("moose-equation-dfeec38c-87e4-4d9f-9872-24801fda2c79");katex.render("\\vec{n}", element, {displayMode:false,throwOnError:false});$ and $var element = document.getElementById("moose-equation-473c07e0-9bcb-488a-9922-50d7a0d99227");katex.render("|S_b|", element, {displayMode:false,throwOnError:false});$ are the outlet face velocity, outward pointing surface vector and the surface area,

The value of $var element = document.getElementById("moose-equation-0017a126-579e-4f9d-ac28-d1494ffab548");katex.render("u_f", element, {displayMode:false,throwOnError:false});$ can be computed two different ways depending on the setting of the "use_two_term_expansion" parameter. When the two-term expansion is enabled the face value is approximated as:

var element = document.getElementById("moose-equation-70de6558-d4d6-420e-8eec-cde1876c18b4");katex.render("u_f = u_C + \\nabla u_C \\cdot \\vec{d}_{Cf}~,", element, {displayMode:true,throwOnError:false});

where $var element = document.getElementById("moose-equation-33bc245d-48a7-41f7-a233-3d7b05f16331");katex.render("u_C", element, {displayMode:false,throwOnError:false});$ and $var element = document.getElementById("moose-equation-3f02305c-2c81-4f42-9cec-dd09f519f95d");katex.render("\\nabla u_C", element, {displayMode:false,throwOnError:false});$ are the solution value and gradient in the boundary cell, while $var element = document.getElementById("moose-equation-aedd7404-8850-4d4f-8d3d-6b0d53a8dc90");katex.render("\\vec{d}_{Cf}", element, {displayMode:false,throwOnError:false});$ is the vector pointing to the face center from the boundary cell centroid. When "use_two_term_expansion" is disabled the following first-order approximation is used:

var element = document.getElementById("moose-equation-140ae1f0-0fb1-466f-bced-c8e7befe7ac6");katex.render("u_f = u_C.", element, {displayMode:true,throwOnError:false});

This boundary condition assumes zero normal gradient contribution to the diffusion terms.

note

This boundary condition should only be used for problems which involve advection and/or diffusion problems.

Example Syntax

[LinearFVBCs]
  [inflow]
    type = LinearFVAdvectionDiffusionFunctorDirichletBC
    variable = u
    boundary = "left top bottom"
    functor = analytic_solution
  []
  [outflow]
    type = LinearFVAdvectionDiffusionOutflowBC
    variable = u
    boundary = "right"
    use_two_term_expansion = false
  []
[]

Input Parameters

boundaryThe list of boundary IDs from the mesh where this object applies
C++ Type:std::vector<BoundaryName>
Unit:(no unit assumed)
Controllable:No
Description:The list of boundary IDs from the mesh where this object applies
variableThe name of the variable that this boundary condition applies to
C++ Type:LinearVariableName
Unit:(no unit assumed)
Controllable:No
Description:The name of the variable that this boundary condition applies to

Required Parameters

use_two_term_expansionFalseIf an approximate linear expansion should be used to compute the face value.
Default:False
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:If an approximate linear expansion should be used to compute the face value.

Optional Parameters

absolute_value_vector_tagsThe tags for the vectors this residual object should fill with the absolute value of the residual contribution
C++ Type:std::vector<TagName>
Unit:(no unit assumed)
Controllable:No
Description:The tags for the vectors this residual object should fill with the absolute value of the residual contribution
extra_matrix_tagsThe extra tags for the matrices this Kernel should fill
C++ Type:std::vector<TagName>
Unit:(no unit assumed)
Controllable:No
Description:The extra tags for the matrices this Kernel should fill
extra_vector_tagsThe extra tags for the vectors this Kernel should fill
C++ Type:std::vector<TagName>
Unit:(no unit assumed)
Controllable:No
Description:The extra tags for the vectors this Kernel should fill
matrix_tagssystemThe tag for the matrices this Kernel should fill
Default:system
C++ Type:MultiMooseEnum
Unit:(no unit assumed)
Options:nontime, system
Controllable:No
Description:The tag for the matrices this Kernel should fill
vector_tagsrhsThe tag for the vectors this Kernel should fill
Default:rhs
C++ Type:MultiMooseEnum
Unit:(no unit assumed)
Options:rhs, time
Controllable:No
Description:The tag for the vectors this Kernel should fill

Tagging 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

Advanced Parameters

(moose/test/tests/linearfvkernels/advection/advection-2d.i)

[Mesh]
  [gmg]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 2
    ny= 1
    ymax = 0.5
  []
[]

[Problem]
  linear_sys_names = 'u_sys'
[]

[Variables]
  [u]
    type = MooseLinearVariableFVReal
    solver_sys = 'u_sys'
    initial_condition = 1.0
  []
[]

[LinearFVKernels]
  [advection]
    type = LinearFVAdvection
    variable = u
    velocity = "0.5 0 0"
    advected_interp_method = upwind
  []
  [source]
    type = LinearFVSource
    variable = u
    source_density = source_func
  []
[]

[LinearFVBCs]
  [inflow]
    type = LinearFVAdvectionDiffusionFunctorDirichletBC
    variable = u
    boundary = "left top bottom"
    functor = analytic_solution
  []
  [outflow]
    type = LinearFVAdvectionDiffusionOutflowBC
    variable = u
    boundary = "right"
    use_two_term_expansion = false
  []
[]

[Functions]
  [source_func]
    type = ParsedFunction
    expression = '0.5*pi*sin(2*y*pi)*cos(x*pi)'
  []
  [analytic_solution]
    type = ParsedFunction
    expression = 'sin(x*pi)*sin(2*y*pi) + 1.5'
  []
[]

[Postprocessors]
  [error]
    type = ElementL2FunctorError
    approximate = u
    exact = analytic_solution
    execute_on = FINAL
  []
  [h]
    type = AverageElementSize
    execute_on = FINAL
  []
[]

[Executioner]
  type = LinearPicardSteady
  linear_systems_to_solve = u_sys
  number_of_iterations = 1
  petsc_options_iname = '-pc_type -pc_factor_shift_type -pc_factor_shift_amount'
  petsc_options_value = 'lu NONZERO               1e-10'
[]

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

Description
Example Syntax
Input Parameters