FVAdvectedUpwind

Overview

This object provides first-order upwind interpolation for an advected quantity on a finite-volume face. The upwind state is selected using the sign of the face mass flux ( $var element = document.getElementById("moose-equation-b6883932-2deb-4e71-96ea-d539a896285c");katex.render("m_f", element, {displayMode:false,throwOnError:false});$ ), yielding a monotone, Total Variation Diminishing (TVD) discretization that is robust near discontinuities at the cost of increased numerical diffusion (Moukalled et al. (2016), Harten (1997)).

Let $var element = document.getElementById("moose-equation-46901b6c-b7c3-4561-943a-0f02cee373bb");katex.render("\\phi_C", element, {displayMode:false,throwOnError:false});$ and $var element = document.getElementById("moose-equation-e9400b0c-3044-4238-83bf-87a624d7c96e");katex.render("\\phi_N", element, {displayMode:false,throwOnError:false});$ denote the cell-centered values on the element and neighbor sides of a face. Define the upwind value

var element = document.getElementById("moose-equation-9a0c5303-d4a0-4e3d-8b06-5f235936a042");katex.render("\\phi_U = \\begin{cases} \\phi_C & \\text{if } m_f \\ge 0,\\\\ \\phi_N & \\text{if } m_f < 0. \\end{cases}", element, {displayMode:true,throwOnError:false});

Then the interpolated face value is simply

var element = document.getElementById("moose-equation-8d3ef826-8d01-4b74-ad6b-3c4d614d9b2d");katex.render("\\phi_f = \\phi_U.", element, {displayMode:true,throwOnError:false});

In matrix form this corresponds to weights $var element = document.getElementById("moose-equation-4a253d7d-206e-40fd-9e5f-b055c222dd9e");katex.render("(w_C, w_N)\\in\\{(1,0),(0,1)\\}", element, {displayMode:false,throwOnError:false});$ chosen by the mass flux sign, so no gradients are required. This is a common baseline scheme for advection and is often used as the low-order component in higher-order schemes with limiting or deferred correction (see Limiters).

Example Syntax

Declare the interpolation method in [FVInterpolationMethods]:

[FVInterpolationMethods<<<{"href": "../../syntax/FVInterpolationMethods/index.html"}>>>]
  [upwind]
    type = FVAdvectedUpwind<<<{"description": "Upwind interpolation for advected quantities using the face mass flux sign.", "href": "FVAdvectedUpwind.html"}>>>
  []
[]

Use it in a linear FV advection kernel via "advected_interp_method_name":

[LinearFVKernels<<<{"href": "../../syntax/LinearFVKernels/index.html"}>>>]
  [advection]
    type = LinearFVAdvection<<<{"description": "Represents the matrix and right hand side contributions of an advection term in a partial differential equation.", "href": "../linearfvkernels/LinearFVAdvection.html"}>>>
    variable<<<{"description": "The name of the variable whose linear system this object contributes to"}>>> = u
    velocity<<<{"description": "Constant advection velocity"}>>> = "0.5 0 0"
    advected_interp_method_name<<<{"description": "Name of the FVInterpolationMethod to use for the advected quantity."}>>> = upwind
  []
[]

Input 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:No
Description:Set the enabled status of the MooseObject.

References

Ami Harten. High resolution schemes for hyperbolic conservation laws. Journal of computational physics, 135(2):260–278, 1997.

@article{harten1997,
    author = "Harten, Ami",
    title = "High resolution schemes for hyperbolic conservation laws",
    journal = "Journal of computational physics",
    volume = "135",
    number = "2",
    pages = "260--278",
    year = "1997",
    publisher = "Elsevier"
}

Fadl Moukalled, L Mangani, Marwan Darwish, and others. The finite volume method in computational fluid dynamics. Volume 6. Springer, 2016.

@book{moukalled2016finite,
    author = "Moukalled, Fadl and Mangani, L and Darwish, Marwan and others",
    title = "The finite volume method in computational fluid dynamics",
    volume = "6",
    year = "2016",
    publisher = "Springer"
}

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

[Mesh]
  [gmg]
    type = GeneratedMeshGenerator
    dim = 1
    nx = 2
  []
[]

[Problem]
  linear_sys_names = 'u_sys'
[]

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

[FVInterpolationMethods]
  [upwind]
    type = FVAdvectedUpwind
  []
  [average]
    type = FVGeometricAverage
  []
  [muscl_venkat]
    type = FVAdvectedVenkatakrishnanDeferredCorrection
    deferred_correction_factor = 1.0
  []
  [nvd_vanleer]
    type = FVAdvectedVanLeerWeightBased
    blending_factor = 1.0
  []
  [nvd_minmod]
    type = FVAdvectedMinmodWeightBased
    blending_factor = 1.0
  []
[]

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

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

[Functions]
  [source_func]
    type = ParsedFunction
    expression = '0.5*x'
  []
  [analytic_solution]
    type = ParsedFunction
    expression = '0.5+0.5*x*x'
  []
[]

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

[Convergence]
  [linear]
    type = IterationCountConvergence
    max_iterations = 1
    converge_at_max_iterations = true
  []
[]

[Executioner]
  type = Steady
  system_names = u_sys
  l_tol = 1e-10
  multi_system_fixed_point=true
  multi_system_fixed_point_convergence=linear
  multi_system_fixed_point_relaxation_factor = 1.0
  petsc_options_iname = '-pc_type -pc_factor_shift_type -pc_factor_shift_amount -pc_factor_mat_solver_type -mat_mumps_icntl_14'
  petsc_options_value = 'lu       NONZERO               1e-12                     mumps                    50'
[]

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

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

[Mesh]
  [gmg]
    type = GeneratedMeshGenerator
    dim = 1
    nx = 2
  []
[]

[Problem]
  linear_sys_names = 'u_sys'
[]

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

[FVInterpolationMethods]
  [upwind]
    type = FVAdvectedUpwind
  []
  [average]
    type = FVGeometricAverage
  []
  [muscl_venkat]
    type = FVAdvectedVenkatakrishnanDeferredCorrection
    deferred_correction_factor = 1.0
  []
  [nvd_vanleer]
    type = FVAdvectedVanLeerWeightBased
    blending_factor = 1.0
  []
  [nvd_minmod]
    type = FVAdvectedMinmodWeightBased
    blending_factor = 1.0
  []
[]

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

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

[Functions]
  [source_func]
    type = ParsedFunction
    expression = '0.5*x'
  []
  [analytic_solution]
    type = ParsedFunction
    expression = '0.5+0.5*x*x'
  []
[]

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

[Convergence]
  [linear]
    type = IterationCountConvergence
    max_iterations = 1
    converge_at_max_iterations = true
  []
[]

[Executioner]
  type = Steady
  system_names = u_sys
  l_tol = 1e-10
  multi_system_fixed_point=true
  multi_system_fixed_point_convergence=linear
  multi_system_fixed_point_relaxation_factor = 1.0
  petsc_options_iname = '-pc_type -pc_factor_shift_type -pc_factor_shift_amount -pc_factor_mat_solver_type -mat_mumps_icntl_14'
  petsc_options_value = 'lu       NONZERO               1e-12                     mumps                    50'
[]

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

Overview
Example Syntax
Input Parameters
References