Navier Stokes System Requirements Specification

This template follows INL template TEM-135, "IT System Requirements Specification".

note

This document serves as an addendum to Framework System Requirements Specification and captures information for SRS specific to the Navier Stokes module.

Framework System Requirements Specification

Introduction

System Purpose

The purpose of this software is to allow simulation of fluid flow in regular and porous media. It should be able to determine the pressure, velocity, and fluid temperature fields, as well as the solid temperature field in the case of conjugate heat transfer simulations. It aims to also be able to transport scalar species in the flow fields determined.

System Scope

The MOOSE Navier Stokes module provides numerical discretizations of the Navier Stokes equations to model the flow of fluid through regular and porous media. The equations discretized include the conservation of mass, momentum, energy and of transported scalars / species. It covers a wide range of fluids and of fluid flow regimes. It covers both natural and forced convection regime, and should be able to model conjugate heat transfer between the fluid and solid phases. A number of scalar species can be transported using the velocity field determined from the fluid flow equations. It can be used as a standalone application or can be included in downstream applications interested in modeling multiphysics problems involving fluid flow.

System Overview

System Context

The Navier Stokes module is command-line driven. Like MOOSE, this is typical for a high-performance software that is designed to run across several nodes of a cluster system. As such, all usage of the software is through any standard terminal program generally available on all supported operating systems. Similarly, for the purpose of interacting through the software, there is only a single user, "the user", which interacts with the software through the command-line. The Navier Stokes module does not maintain any back-end database or interact with any system daemons. It is an executable, which may be launched from the command line and writes out various result files as it runs.

Figure 1: Usage of the Navier Stokes module and other MOOSE-based applications.

System Functions

Since the Navier Stokes module is a command-line driven application, all functionality provided in the software is operated through the use of standard UNIX command line flags and the extendable MOOSE input file. The Navier Stokes module is completely extendable so individual design pages should be consulted for specific behaviors of each user-defined object.

User Characteristics

Like MOOSE, there are three kinds of users working on the Navier Stokes module:

Navier Stokes module Developers: These are the core developers of the Navier Stokes module. They are responsible for following and enforcing the software development standards of the module, as well as designing, implementing, and maintaining the software.
Developers: A scientist or engineer that uses the Navier Stokes module alongside MOOSE to build their own application. This user will typically have a background in modeling or simulation techniques (and perhaps numerical analysis) but may only have a limited skillset when it comes to code development using the C++ language. This is the primary focus group of the module. In many cases, these developers will be encouraged to contribute module-appropriate code back to the Navier Stokes module, or to MOOSE itself.
Analysts: These are users that will run the code and perform analysis on the simulations they perform. These users may interact with developers of the system requesting new features and reporting bugs found and will typically make heavy use of the input file format.

Assumptions and Dependencies

The Navier Stokes module is developed using MOOSE and can itself be based on various MOOSE modules, as such the SRS for the Navier Stokes module is dependent upon the files listed at the beginning of this document. Any further assumptions or dependencies are outlined in the remainder of this section.

The Navier Stokes module is designed with the fewest possible constraints on hardware and software. For more context on this point, the Navier Stokes module SRS defers to the framework Assumptions and Dependencies. The Navier Stokes module depends on the fluid properties module for fluid properties. All physics-based assumptions in this module's objects are highlighted in each object's respective documentation pages. A summary of the key assumptions may be found on the Navier Stokes module homepage. The key physics-based assumptions when modeling fluid flow consist of the compressibility, from incompressible to fully compressible, the flow regime, whether turbulent or regular, and whether the media is porous or not, which induces volumetric friction. The key numerical decisions to be made when using the Navier Stokes module are whether to use a finite volume or finite element discretization, and within those which numerical scheme to use, notably for discretizing the advective term.

References

ISO/IEC/IEEE 24765:2010(E). Systems and software engineering—Vocabulary. first edition, December 15 2010.

@manual{ISO-systems-software,
    author = "24765:2010(E), ISO/IEC/IEEE",
    title = "Systems and software engineering---Vocabulary",
    edition = "First",
    year = "2010",
    month = "December 15"
}

ASME NQA-1. ASME NQA-1-2008 with the NQA-1a-2009 addenda: Quality Assurance Requirements for Nuclear Facility Applications. first edition, August 31 2009.

@manual{ASME-NQA-1-2008,
    author = "NQA-1, ASME",
    title = "ASME NQA-1-2008 with the NQA-1a-2009 addenda: Quality Assurance Requirements for Nuclear Facility Applications",
    orginization = "American Socieity of Mechanical Engineers",
    year = "2009",
    month = "August 31",
    edition = "First"
}

Definitions and Acronyms

This section defines, or provides the definition of, all terms and acronyms required to properly understand this specification.

Definitions

Verification: (1) The process of: evaluating a system or component to determine whether the products of a given development phase satisfy the conditions imposed at the start of that phase. (2) Formal proof of program correctness (e.g., requirements, design, implementation reviews, system tests) (24765:2010(E), 2010).

Acronyms

Acronym	Description
INL	Idaho National Laboratory
LGPL	GNU Lesser General Public License
MOOSE	Multiphysics Object Oriented Simulation Environment
NQA-1	Nuclear Quality Assurance Level 1
POSIX	Portable Operating System Interface
SRS	Software Requirement Specification

System Requirements

In general, the following is required for MOOSE-based development:

A POSIX compliant Unix-like operating system. This includes any modern Linux-based operating system (e.g., Ubuntu, Fedora, Rocky, etc.), or a Macintosh machine running either of the last two MacOS releases.

Hardware	Information
CPU Architecture	x86_64, ARM (Apple Silicon)
Memory	8 GB (16 GBs for debug compilation)
Disk Space	30GB

Libraries	Version / Information
GCC	8.5.0 - 12.2.1
LLVM/Clang	10.0.1 - 16.0.6
Intel (ICC/ICX)	Not supported at this time
Python	3.7 - 3.11
Python Packages	packaging pyaml jinja2

Functional Requirements

navier_stokes: Auxkernels
5.1.1The system shall be able to compute the liquid fraction based on solidus and liquidus temperatures.
5.1.2
The system shall be able to compute element thermal Peclet numbers for:
1. finite volume discretizations
2. finite element discretizations
5.1.3
The system shall be able to compute element Reynolds numbers for:
1. finite volume discretizations
2. finite element discretizations

navier_stokes: Finite Element
5.2.1The system shall be able to solve the Euler equations for subsonic flow with a bump in the channel.
5.2.2The system shall be able to solve the incompressible Navier-Stokes equations in an RZ coordinate system while not integrating the pressure term by parts.
5.2.3The system shall be able to solve the incompressible Navier-Stokes equations in an RZ coordinate system while integrating the pressure term by parts.
5.2.4The system shall be able to solve the incompressible Navier-Stokes equations for a high Reynolds number in an RZ coordinate system.
5.2.5The system shall be able to compute an accurate Jacobian for the incompressible Navier-Stokes equations in an RZ coordinate system.
5.2.6The system shall be able to solve the transient incompressible Navier-Stokes equations using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts and reproduce the results of a hand-coded Jacobian implementation.
5.2.7The system shall be able to solve the transient incompressible Navier-Stokes equations using an automatic differentiation, vector variable implementation in an RZ coordinate system while not integrating the pressure term by parts, using a traction form for the viscous term, and using a no-bc boundary condition, and reproduce the results of a hand-coded Jacobian implementation.
5.2.8The system shall be able to solve the steady incompressible Navier-Stokes equations using an automatic differentiation, vector variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a natural outflow boundary condition.
5.2.9The system shall be able to solve the steady incompressible Navier-Stokes equations using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a natural outflow boundary condition
5.2.10The system shall be able to solve the steady incompressible Navier-Stokes equations using an automatic differentiation, vector variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a NoBC outflow boundary condition.
5.2.11The system shall be able to solve the steady incompressible Navier-Stokes equations using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a NoBC outflow boundary condition
5.2.12The system shall be able to solve the steady incompressible Navier-Stokes equations using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a natural outflow boundary condition and reproduce the results of the AD, vector variable implementation.
5.2.13The system shall be able to solve the steady incompressible Navier-Stokes equations using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a natural outflow boundary condition and reproduce the results of the AD, vector variable implementation.
5.2.14The system shall be able to solve the steady incompressible Navier-Stokes equations using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a NoBC outflow boundary condition and reproduce the results of the AD, vector variable implementation.
5.2.15The system shall be able to solve the steady incompressible Navier-Stokes equations using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a NoBC outflow boundary condition and reproduce the results of the AD, vector variable implementation.
5.2.16The system shall be able to solve the steady incompressible Navier-Stokes equations with SUPG and PSPG stabilization and a first order velocity basis using an automatic differentiation, vector variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a natural outflow boundary condition.
5.2.17The system shall be able to solve the steady incompressible Navier-Stokes equations with SUPG and PSPG stabilization and a first order velocity basis using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a natural outflow boundary condition and reproduce the results of the AD, vector variable implementation.
5.2.18The system shall be able to solve the steady incompressible Navier-Stokes equations with SUPG and PSPG stabilization and a first order velocity basis using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a natural outflow boundary condition.
5.2.19The system shall be able to solve the steady incompressible Navier-Stokes equations with SUPG and PSPG stabilization and a first order velocity basis using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a natural outflow boundary condition and reproduce the results of the AD, vector variable implementation.
5.2.20The system shall be able to solve the steady incompressible Navier-Stokes equations with SUPG and PSPG stabilization and a second order velocity basis using an automatic differentiation, vector variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a natural outflow boundary condition.
5.2.21The system shall be able to solve the steady incompressible Navier-Stokes equations with SUPG and PSPG stabilization and a second order velocity basis using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a natural outflow boundary condition.
5.2.22The system shall be able to solve the steady incompressible Navier-Stokes equations with SUPG and PSPG stabilization and a second order velocity basis using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a natural outflow boundary condition.
5.2.23The system shall be able to solve the steady incompressible Navier-Stokes equations with SUPG and PSPG stabilization and a second order velocity basis using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a natural outflow boundary condition.
5.2.24The system shall compute an accurate Jacobian using automatic differentiation when solving the incompressible Navier Stokes equations in an axisymmetric coordinate system with SUPG and PSPG stabilization
5.2.25The system shall be able to solve the steady incompressible Navier-Stokes equations with SUPG and PSPG stabilization and a first order velocity basis using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts, using a traction form for the viscous term, and applying a natural outflow boundary condition.
5.2.26The system shall be able to solve the steady incompressible Navier-Stokes equations with SUPG and PSPG stabilization and a first order velocity basis using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while integrating the pressure term by parts, using a traction form for the viscous term, and applying a natural outflow boundary condition and reproduce the results of the AD, vector variable implementation.
5.2.27The system shall be able to solve the steady incompressible Navier-Stokes equations with SUPG and PSPG stabilization and a first order velocity basis using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts, using a traction form for the viscous term, and applying a natural outflow boundary condition and obtain a perfect Jacobian.
5.2.28
The system shall be able to solve the steady incompressible Navier-Stokes equations in an axisymmetric coordinate system, using a Jacobian computed via automatic differentiation, on a displaced mesh, with the viscous term in
1. traction form
2. laplace form
5.2.29The system shall be able to reproduce the unstabilized conical flow results from a stabilized continuous finite element implementation using generic unstabilized objects.
5.2.30The system shall compute inflow and outflow boundary conditions for advected variables
5.2.31We shall error if the user provides less velocity components than the mesh dimension
5.2.32We shall error if the user provides more than 3 velocity components
5.2.33We shall allow the user to supply more velocity components than the mesh dimension (up to 3 components)
5.2.34The system shall be able to solve two different kernel sets with two different material domains.
5.2.35The system shall be able to solve two different kernel sets within one material domain.
5.2.36The system shall be able to solve one kernel set with two different material domains.
5.2.37The system shall be able to reproduce benchmark results for a Rayleigh number of 1e3.
5.2.38The system shall be able to reproduce benchmark results for a Rayleigh number of 1e4.
5.2.39The system shall be able to reproduce benchmark results for a Rayleigh number of 1e5.
5.2.40The system shall be able to reproduce benchmark results for a Rayleigh number of 1e6.
5.2.41The system shall be able to simulate natural convection by adding the Boussinesq approximation to the incompressible Navier-Stokes equations.
5.2.42The system shall support passing constants for material property names to a kernel which passes its parameters to a material.
5.2.43The system shall be able to solve mass, momentum, and energy incompressible Navier-Stokes equations with multiple threads.
5.2.44The system shall have an accurate Jacobian provided by automatic differentiation when computing the Boussinesq approximation.
5.2.45The system shall be able to support SUPG and PSPG stabilization of the incompressible Navier Stokes equations including the Boussinesq approximation.
5.2.46The system shall be able to solve stablized mass, momentum, and energy incompressible Navier-Stokes equations with multiple threads.
5.2.47The system shall have an accurate Jacobian provided by automatic differentiation when computing the Boussinesq approximation with SUPG and PSPG stabilization.
5.2.48The system shall be able to reproduce results of incompressible Navier-Stokes with Boussinesq approximation using a customized and condensed action syntax.
5.2.49The system shall be able to solve mass, momentum, and energy incompressible Navier-Stokes equations with a custom action syntax using multiple threads.
5.2.50The system shall be able to solve a channel flow problem using a hybrid CG-DG discretization with first Lagrange pressure and first monomial velocity.
5.2.51
The system shall be able to solve a lid driven cavity problem for a Reynolds number of 200 using a hybrid CG-DG scheme in which the pressure is first order Lagrange and the velocity is first order monomial and show
1. accurate results, and
2. a perfect Jacobian.
5.2.52The system shall be able to solve the incompressible Navier-Stokes equations on triangular meshes, using a hybrid CG-DG scheme with first order Lagrange pressure and second order monomial velocity, with Dirichlet boundary conditions for the velocity, and demonstrate third order convergence for the velocity and second order convergence for pressure in the L2 error measure.
5.2.53The system shall be able to solve the incompressible Navier-Stokes equations on triangular meshes, using a hybrid CG-DG scheme with first order Lagrange pressure and first order monomial velocity, with Dirichlet boundary conditions for the velocity, and demonstrate second order convergence for the velocity and first order convergence for pressure in the L2 error measure.
5.2.54The system shall be able to solve the incompressible Navier-Stokes equations using a hybrid CG-DG scheme with first order Lagrange pressure and first order monomial velocity and demonstrate second order convergence for the velocity and first order convergence for pressure in the L2 error measure.
5.2.55The system shall be able to solve the incompressible Navier-Stokes equations using a hybrid CG-DG scheme with first order Lagrange pressure and second order monomial velocity and demonstrate third order convergence for the velocity and second order convergence for pressure in the L2 error measure.
5.2.56The system shall be able to solve the incompressible Navier-Stokes equations using a hybrid CG-DG scheme with first order Lagrange pressure and second order L2 Lagrange velocity and demonstrate third order convergence for the velocity and second order convergence for pressure in the L2 error measure.
5.2.57The system shall be able to solve the incompressible Navier-Stokes equations using a hybrid CG-DG scheme with first order Lagrange pressure and second order L2 hierarchic velocity and demonstrate third order convergence for the velocity and second order convergence for pressure in the L2 error measure.
5.2.58The system shall be able to apply an external force to the incompressible Navier-Stokes momentum equation through a coupled variable.
5.2.59The system shall be able to compute an accurate Jacobian when applying an external force to the incompressible Navier-Stokes momentum equation through a coupled variable.
5.2.60The system shall be able to apply an external force to the incompressible Navier-Stokes momentum equation through a vector function.
5.2.61The system shall be able to compute an accurate Jacobian when applying an external force to the incompressible Navier-Stokes momentum equation through a vector function.
5.2.62The system shall be able to apply an external force to the incompressible Navier-Stokes momentum equation through a coupled variable, with the problem setup through automatic action syntax.
5.2.63The system shall be able to compute an accurate Jacobian when applying an external force to the incompressible Navier-Stokes momentum equation through a coupled variable, with the problem setup through automatic action syntax.
5.2.64The system shall be able to apply an external force to the incompressible Navier-Stokes momentum equation through a vector function, with the problem setup through automatic action syntax.
5.2.65The system shall be able to compute an accurate Jacobian when applying an external force to the incompressible Navier-Stokes momentum equation through a vector function, with the problem setup through automatic action syntax.
5.2.66
The system shall be able to solve the Navier-Stokes equations with a coupled variable force and a gravity force
1. provided through a dedicated object,
2. or through a generic object that can simultaneously add multiple forces through both a coupled variable and a function.
3. The generic object shall also be able to compute the forces solely through multiple coupled variables,
4. or solely through multiple vector functions.
5. The system shall be able to add the generic object through an automatic action syntax and provide two forces either through a coupled variable and a function,
6. two coupled variables,
7. or two functions.
5.2.67
The system shall exhibit global conservation of energy when using the continuous-Galerkin finite element spatial discretization with streamline-upwind Petrov-Galerkin stabilization and with
1. q2q1 elements
2. q1q1 elements and pressure-stabilized Petrov-Galerkin stabilization
5.2.68The system shall be able to model a volumetric heat source and included it in stabilization terms.
5.2.69The system shall be able to build a volumetric heat source model using an automatic action syntax.
5.2.70The system shall be able to model a volumetric heat source with a coupled variable and included it in stabilization terms.
5.2.71The system shall be able to build a volumetric heat source model, provided through a coupled variable, using an automatic action syntax.
5.2.72The system shall be able to solve the incompressible Navier-Stokes equations in a cavity using a hybridized discontinuous Galerkin scheme with broken Lagrange basis and produce second order convergence for all variables.
5.2.73The system shall be able to solve the incompressible Navier-Stokes equations in a cavity using a hybridized discontinuous Galerkin scheme with broken Hierarchic basis and produce second order convergence for all variables.
5.2.74The system shall be able to solve the incompressible Navier-Stokes equations in a channel using a hybridized discontinuous Galerkin scheme with broken Lagrange basis and produce second order convergence for all variables.
5.2.75The system shall be able to solve the incompressible Navier-Stokes equations in a channel using a hybridized discontinuous Galerkin scheme with broken Hierarchic basis and produce second order convergence for all variables.
5.2.76The system shall be able to solve a lid-driven cavity problem using a hybridized discontinuous Galerkin discretization.
5.2.77The system shall be able to solve a channel flow problem using a hybridized discontinuous Galerkin discretization.
5.2.78
The system shall produce a symmetric matrix for a hybridizable discontinuous Galerkin discretization of the Stokes equation for a
1. lid driven cavity problem, and
2. channel flow problem.
5.2.79The system shall error if hybridized discontinuous Galerkin kernels and boundary conditions implement different physics.
5.2.80The system shall be able to model the effect of gravity on incompressible flow using a finite element discretization.
5.2.81The system shall compute accurate Jacobians for the incompressible Navier-Stokes equation.
5.2.82The system shall compute accurate Jacobians for the incompressible Navier-Stokes equation with stabilization.
5.2.83The system shall compute accurate Jacobians for the incompressible Navier-Stokes equation with stabilization with a traction boundary condition.
5.2.84
The system shall be able to solve Jeffery-Hamel flow in a 2D wedge and compare to the analytical solution
1. with pressure Dirichlet boundary conditions
2. and with natural advection boundary conditions.
5.2.85The system shall support solving a steady energy equation and transient momentum equations and apply the correct stabilization.
5.2.86The system shall support solving a steady energy equation and transient momentum equations with correct stabilization and compute a perfect Jacobian.
5.2.87We shall be able to solve a canonical lid-driven problem without stabilization, using mixed order \ finite elements for velocity and pressure.
5.2.88Least squares commutator (LSC) preconditioning shall require only a small number of linear iterations to converge when using Newton with a Reynolds number of unity, performing a transient march to steady-state.
5.2.89Least squares commutator (LSC) preconditioning shall require only a small number of linear iterations to converge when using Picard with a Reynolds number of 500, performing a few timesteps in a transient.
5.2.90Least Squares Commutator (LSC) preconditioning shall require only a small number of linear iterations to converge when using Newton with a Reynolds number of 500, performing a transient march to steady-state.
5.2.91The system shall be able to solve a canonical lid-driven problem using same order variables and the PSPG/SUPG stabilization scheme.
5.2.92
The system shall be able to run scalable Taylor-Hood finite element simulations using a field split preconditioner with least-squares commutator preconditioning for the Schur complement and multigrid for the sub-solves, using a steady executioner with a Reynolds number of unity, and
1. using the velocity mass matrix for scaling within the Least Squares Commutator preconditioner with a component variable implementation, or
2. using the diagonal of the velocity-velocity block for scaling within the Least Squares Commutator preconditioner with a component variable implementation, or
3. using the velocity mass matrix for scaling within the Least Squares Commutator preconditioner with a vector variable implementation, or
4. commuting operators in the style of Olshanskii with a vector variable implementation.
5.2.93The system shall be able to run scalable Taylor-Hood finite element simulations using a field split preconditioner with least-squares commutator preconditioning for the Schur complement and multigrid for the sub-solves, using a steady executioner with a Reynolds number of 500.
5.2.94
The system shall be able to efficiently precondition Taylor-Hood finite elements using a Schur complement field split with the preconditioner for the Schur complement formed from the pressure mass matrix for a
1. a Stokes problem and a
2. a high Reynolds number (1000) Navier-Stokes problem with a grad-div stabilization.
5.2.95We shall be able to reproduce the results from the hand-coded lid-driven simulation using automatic differentiation objects.
5.2.96We shall be able to run lid-dirven simulation using a global mean-zero pressure constraint approach.
5.2.97The Jacobian for the mixed-order INS problem shall be perfect when provided through automatic differentiation.
5.2.98We shall be able to solve the lid-driven problem using equal order shape functions with pressure-stabilized petrov-galerkin stabilization. We shall also demonstrate SUPG stabilization.
5.2.99We shall be able to reproduce the hand-coded stabilized results with automatic differentiation objects.
5.2.100The Jacobian for the automatic differentiation stabilized lid-driven problem shall be perfect.
5.2.101Simulation with equal-order shape functions without pressure stabilization shall be unstable.
5.2.102We shall be able to solve the INS equations using the classical Chorin splitting algorithm.
5.2.103The system shall be able to reproduce unstabilized incompressible Navier-Stokes results with hand-coded Jacobian using a customized and condensed action syntax.
5.2.104The system shall be able to reproduce stabilized incompressible Navier-Stokes results with hand-coded Jacobian using a customized and condensed action syntax.
5.2.105The system shall be able to reproduce unstabilized incompressible Navier-Stokes results with auto-differentiation using a customized and condensed action syntax.
5.2.106The system shall be able to reproduce stabilized incompressible Navier-Stokes results with auto-differentiation using a customized and condensed action syntax.
5.2.107The system shall be able to solve a steady stabilized mass/momentum/energy incompressible Navier-Stokes formulation.
5.2.108The system shall be able to solve a transient stabilized mass/momentum/energy incompressible Navier-Stokes formulation.
5.2.109The system shall be able to solve a steady stabilized mass/momentum/energy incompressible Navier-Stokes formulation with action syntax.
5.2.110The system shall be able to solve a transient stabilized mass/momentum/energy incompressible Navier-Stokes formulation with action syntax.
5.2.111The system shall be able to solve a transient incompressible Navier-Stokes with nonlinear Smagorinsky eddy viscosity.
5.2.112
The system shall be able to apply pressure stabilization using an alpha parameter of 1e-6 on a
1. 4x4,
2. 8x8,
3. 16x16,
4. and 32x32 mesh.
5.2.113
The system shall be able to apply pressure stabilization using an alpha parameter of 1e-3 on a
1. 4x4,
2. 8x8,
3. 16x16,
4. and 32x32 mesh.
5.2.114
The system shall be able to apply pressure stabilization using an alpha parameter of 1e0 on a
1. 4x4,
2. 8x8,
3. 16x16,
4. and 32x32 mesh.
5.2.115
The system shall be able to apply streamline-upwind stabilization using an alpha parameter of 1e-6 on a
1. 4x4,
2. 8x8,
3. 16x16,
4. and 32x32 mesh.
5.2.116
The system shall be able to apply streamline-upwind stabilization using an alpha parameter of 1e-3 on a
1. 4x4,
2. 8x8,
3. 16x16,
4. and 32x32 mesh.
5.2.117
The system shall be able to apply streamline-upwind stabilization using an alpha parameter of 1e0 on a
1. 4x4,
2. 8x8,
3. 16x16,
4. and 32x32 mesh.
5.2.118The system shall be able to solve high Reynolds number incompressible flow problems through use of streamline upwind Petrov-Galerkin stabilization and with a Q2Q1 discretization
5.2.119The system shall be able to solve high Reynolds number incompressible flow problems through use of streamline upwind and pressure stabilized Petrov-Galerkin and with a Q1Q1 discretization
5.2.120The system shall allow MOOSE applications to specify nonzero malloc behavior; for the Navier-Stokes application, new nonzero allocations shall be errors.
5.2.121
The system shall be able to solve for incompressible fluid flowing through a 2D channel driven by pressure inlet and outlet boundary conditions
1. using the kernel formulation,
2. using the action formulation
3. and using a field split preconditioning.
5.2.122
The system shall be able to solve an axisymmetric pipe flow problem using a finite element discretization in which the axis of symmetry is the x-axis, using a Laplace form for the viscous term
1. in which the pressure is constrained with a Dirichlet boundary condition on the outlet, using a Jacobian computed via automatic differentiation, without PSPG and SUPG stabilization,
2. in which the pressure is constrained with a Dirichlet boundary condition on the outlet, using a Jacobian computed via user-provided functions, without PSPG and SUPG stabilization,
3. in which the pressure is constrained with a Dirichlet boundary condition on the outlet, using a Jacobian computed via automatic differentiation, with PSPG and SUPG stabilization,
4. in which the pressure is constrained with a Dirichlet boundary condition on the outlet, using a Jacobian computed via user-provided functions, with PSPG and SUPG stabilization,
5. in which the pressure is constrained with natural boundary conditions for the velocity equations on the outlet, using a Jacobian computed via automatic differentiation, with PSPG and SUPG stabilization,
6. and in which the pressure is constrained with natural boundary conditions for the velocity equations on the outlet, using a Jacobian computed via user-provided functions, with PSPG and SUPG stabilization.
5.2.123
The system shall be able to solve an axisymmetric pipe flow problem using a finite element discretization in which the axis of symmetry is the x-axis, using a traction form for the viscous term
1. in which the pressure is constrained with a Dirichlet boundary condition on the outlet, using a Jacobian computed via automatic differentiation, without PSPG and SUPG stabilization,
2. in which the pressure is constrained with a Dirichlet boundary condition on the outlet, using a Jacobian computed via user-provided functions, without PSPG and SUPG stabilization,
3. in which the pressure is constrained with a Dirichlet boundary condition on the outlet, using a Jacobian computed via automatic differentiation, with PSPG and SUPG stabilization,
4. in which the pressure is constrained with a Dirichlet boundary condition on the outlet, using a Jacobian computed via user-provided functions, with PSPG and SUPG stabilization,
5. in which the pressure is constrained with natural boundary conditions for the velocity equations on the outlet, using a Jacobian computed via automatic differentiation, with PSPG and SUPG stabilization,
6. and in which the pressure is constrained with natural boundary conditions for the velocity equations on the outlet, using a Jacobian computed via user-provided functions, with PSPG and SUPG stabilization.
5.2.124
The system shall be able to solve for an incompressible fluid flowing through a 1D channel with Streamline Upwind Petrov Galerkin stabilization.
1. using the optimal tau stabilization,
2. using the modified tau stabilization,
3. and still satisfy MMS testing in 1D
4. and in 2D.
5.2.125
The system shall be able to solve for incompressible fluid evolving in a corner cavity with Dirichlet boundary conditions on velocity.
1. in 2D
2. and in 2D RZ axisymmetric simulations.
5.2.126
The system shall be able to solve for incompressible fluid flowing through a 2D channel with only inlet velocity boundary conditions
1. with the regular volumetric integration of the momentum pressure term
2. and with the momentum pressure term integrated by parts.
5.2.127
The system shall be able to solve for incompressible fluid flowing through a 2D channel with only inlet velocity boundary conditions with streamline upwind Petrov Galerkin stabilization and a traction form for the viscous term
1. using a hand-coded Jacobian, and
2. using an automatic differentiation computed Jacobian and compute identical results, indicating the traction implementations with second order derivatives for the viscous term in the stabilization term are identical.
5.2.128The system shall be able to model heat transfer from ambient surroundings using a volumetric approximation.
5.2.129The system shall be able to build a simulation, modeling heat transfer from ambient surroundings, using an automated action syntax.
5.2.130The system shall be able to add a incompressible Navier-Stokes energy/temperature equation using an action, but use a temperature variable already added in the input file.
5.2.131The system shall be able to solve a porous flow with pressure drop due to viscous and inertia frictions.
5.2.132The system shall be able to solve a pressure gradient driven porous flow.
5.2.133The system shall be able to solve a porous flow with internal heating source.
5.2.134The system shall be able to solve a porous flow when reverse flow happens.
5.2.135The system shall be able to solve a porous flow using slip-wall boundary condition.

navier_stokes: Finite Volume
5.3.1The system shall be able to use Weller's reconstruction method to build cell vector values from face fluxes with a second order convergence.
5.3.2The system shall be able to solve the 1D Sod shock-tube benchmark problem using an HLLC scheme to compute convective fluxes.
5.3.3The system shall be able to solve the steady Euler equations in a heated channel using Kurganov-Tadmor with linearly reconstructed data with Van-Leer limiting for the convection term and a primitive variable set and show a flat momentum profile
5.3.4The system shall be able to impose boundary advective fluxes for HLLC discretizations that use implicit/interior cell information.
5.3.5The system shall exhibit first order convergence for all variables for the free-flow Euler equations with added artificial diffusion using a HLLC discretization scheme for the advection flux and with specified temperature and momentum at one boundary and specified pressure at another boundary.
5.3.6
The system shall exhibit first order convergence for all variables for the porous Euler equations using a HLLC discretization scheme for the advection flux and with specified temperature and momentum/velocity at one boundary and specified pressure at another boundary with
1. constant porosity
2. varying porosity
5.3.7
The system shall be able to use a primitive variable set and compute intercell fluxes using a Kurganov-Tadmor scheme
1. when using central differencing to interpolate cell center values to faces and display second order convergence
2. when using directional upwinding to interpolate cell center values to faces and display first order convergence
3. when using linear interpolation of cell center values to faces with Van-Leer limiting and display at least second order convergence
5.3.8
The system shall be able to use a conserved variable set and compute intercell fluxes using a Kurganov-Tadmor scheme
1. when using central differencing to interpolate cell center values to faces and display second order convergence
2. when using directional upwinding to interpolate cell center values to faces and display first order convergence
3. when using linear interpolation of cell center values to faces with Van-Leer limiting and display at least second order convergence
5.3.9
The system shall be able to solve a problem with continuously varying porosity provided through a function object, using a primitive variable set, and compute intercell fluxes using a Kurganov-Tadmor (KT) scheme with the KT Method for computing omega
1. when using central differencing to interpolate cell center values to faces and display second order convergence
2. when using directional upwinding to interpolate cell center values to faces and display first order convergence
3. when using linear interpolation of cell center values to faces with Van-Leer limiting and display at least second order convergence
5.3.10
The system shall be able to solve a problem with continuously varying porosity provided through a function object, using a primitive variable set, and compute intercell fluxes using a Kurganov-Tadmor scheme with the Kurganov-Noelle-Petrova method for computing omega
1. when using central differencing to interpolate cell center values to faces and display second order convergence
2. when using directional upwinding to interpolate cell center values to faces and display first order convergence
3. when using linear interpolation of cell center values to faces with Van-Leer limiting and display at least second order convergence
5.3.11
The system shall be able to solve a problem with continuously varying porosity provided through a function object, using a mixed variable set, and compute intercell fluxes using a Kurganov-Tadmor scheme
1. when using central differencing to interpolate cell center values to faces and display second order convergence
2. when using directional upwinding to interpolate cell center values to faces and display first order convergence
5.3.12
The system displays issues when trying to solve hyperbolic equations with sources when using a Godunov method with HLLC approximate Riemann solver on an irregular grid
1. when the source has a cell-centered volumetric discretization
5.3.13
On a regular grid, using a HLLC scheme to calculate inter-cell fluxes, the system shall show, via the momentum variable
1. conservation of mass when no sources are present
2. violation of conservation of mass when sources are present
3. lesser violation of conservation of mass when sources are present and the mesh is refined
5.3.14The system shall demonstrate first order convergence rates for pressure and superficial velocity when using an upwind interpolation for advected quantities in a weakly compressible formulation of the mass and momentum Euler equations.
5.3.15The system shall be able to model subsonic nozzle flow using an HLLC discretization with a specified outlet pressure.
5.3.16The system shall be able to advect a scalar using density and velocity computed through solution of the Euler equations.
5.3.17The system shall be able to run a two-dimensional version of Sod's shocktube problem.
5.3.18The system shall be able to model supersonic nozzle flow using an HLLC advective flux discretization and with inlet boundary conditions based on stagnation temperature and stagnation pressure.
5.3.19
The system shall be able to solve a series of stages of continuous porosity changes with different schemes for computing the convective fluxes assuming piecewise constant data including
1. the Kurganov-Tadmor scheme
2. the HLLC scheme
5.3.20The system shall be able to solve a two-dimensional y-channel problem with frictional drag and a series of porosity jumps smoothed into a continuous porosity function, using the Kurganov-Tadmor scheme for computing intercell convective fluxes with upwind limiting interpolation (e.g. the + cell centroid value is used as the + side value at the face).
5.3.21The system shall be able to solve a two-dimensional y-channel problem using a mixed variable set with frictional drag and a series of porosity jumps smoothed into a continuous porosity function, using the Kurganov-Tadmor scheme for computing intercell convective fluxes with upwind limiting interpolation (e.g. the + cell centroid value is used as the + side value at the face).
5.3.22The system shall support the deferred correction algorithm for transitioning from low-order to high-order representations of the convective flux during a transient simulation.
5.3.23The system shall be able to run a two-dimensional symmetric flow problem with an HLLC discretization for advection.
5.3.24The system shall be able to compute wave speeds for HLLC Riemann solvers.
5.3.25The system shall be able to compute wave speeds for HLLC Riemann solvers in multiple dimensions.
5.3.26The system shall be able to switch Dirichlet boundary conditions on and off over the course of a transient in finite volume contexts.
5.3.27The system shall be able to switch the behavior of velocity and pressure boundary conditions during runtime.
5.3.28The system shall provide a boundary condition to split a constant heat flux according to local values of porosity, using functor material properties.
5.3.29The system shall provide a boundary condition to split a constant heat flux according to domain-averaged values of porosity, using functor material properties.
5.3.30The system shall provide a boundary condition to split a constant heat flux according to local values of thermal conductivity, using functor material properties.
5.3.31The system shall provide a boundary condition to split a constant heat flux according to domain-averaged values of thermal conductivity, using functor material properties.
5.3.32The system shall provide a boundary condition to split a constant heat flux according to local values of effective thermal conductivity, using functor material properties.
5.3.33The system shall provide a boundary condition to split a constant heat flux according to domain-averaged values of effective thermal conductivity, using functor material properties.
5.3.34The system shall provide a boundary condition to split a constant heat flux according to local values of porosity.
5.3.35The system shall provide a boundary condition to split a constant heat flux according to domain-averaged values of porosity.
5.3.36The system shall provide a boundary condition to split a constant heat flux according to local values of thermal conductivity.
5.3.37The system shall provide a boundary condition to split a constant heat flux according to domain-averaged values of thermal conductivity.
5.3.38The system shall provide a boundary condition to split a constant heat flux according to local values of effective thermal conductivity.
5.3.39The system shall provide a boundary condition to split a constant heat flux according to domain-averaged values of effective thermal conductivity.
5.3.40The system shall be able to impose a wall shear stress at the wall according to the algebraic wall function.
5.3.41
The system shall be able to solve for wall-convection with a user-specified heat transfer coefficient
1. for a cavity problem
2. and for a channel problem.
5.3.42The system shall throw an error if the number of momentum inlet types does not match the number of inlet boundaries.
5.3.43The system shall throw an error if the number of momentum inlet functions is not equal to the problem dimension for a fixed-velocity inlet.
5.3.44The system shall throw an error if the number of entries for momentum inlet types does not match the total number of fixed-velocity and fixed-pressure inlet boundaries.
5.3.45The system shall throw an error if the number of momentum outlet types does not match the number of outlet boundaries.
5.3.46The system shall throw an error if the number of pressure outlet functions is not the same the pressure outlet boundaries.
5.3.47The system shall throw an error if the number of momentum wall types does not match the number of wall boundaries.
5.3.48The system shall throw an error if the number of energy inlet types does not match the number of inlet boundaries.
5.3.49The system shall throw an error if the number of passive scalar inlet types does not match the number of inlet boundaries.
5.3.50The system shall throw an error if the number of passive scalar inlet function blocks does not match the number of scalar variables.
5.3.51The system shall throw an error if the number of passive scalar inlet functions does not match the number of inlet boundaries for a specific scalar variable.
5.3.52The system shall throw an error if the number of energy wall types does not match the number of wall boundaries.
5.3.53The system shall throw an error if the number of energy wall functions does not match the number of energy wall types.
5.3.54The system shall throw an error if the number of defined initial conditions is different than the number of created scalar variables.
5.3.55The system shall throw an error if the number of components for the initial velocity is not equal to the dimension or 3.
5.3.56The system shall throw an error if friction correction is requested with no porous medium treatment.
5.3.57The system shall throw an error if consistent scaling is defined without using friction correction
5.3.58The system shall throw an error if porosity smoothing is requested without porous medium treatment.
5.3.59The system shall throw an error if two-term extrapolation is elected for prosity jump faces without enabling the Bernoulli treatment.
5.3.60The system shall throw an error if the user defines an inappropriate number of passive scalar diffusivities
5.3.61The system shall throw an error if the user defines an inappropriate number of passive scalar source functions
5.3.62The system shall throw an error if the user defines an inappropriate number of passive scalar coupled source functions
5.3.63The system shall throw an error if the user defines an inappropriate number of passive scalar coupled source coefficients with regards to the number of scalar equations
5.3.64The system shall throw an error if the user defines an inappropriate number of passive scalar coupled source coefficients with regards to the number of sources
5.3.65The system shall throw an error if the user defines an inappropriate number of passive scalar Schmidt numbers
5.3.66The system shall throw an error if the user supplies a velocity variable which does not exist
5.3.67The system shall throw an error if the user supplies an inappropriate number of externally created velocity components
5.3.68The system shall throw an error if the user supplies unallowed names for the external velocity components in a porous medium setting
5.3.69The system shall throw an error if the user supplies a pressure variable which does not exist
5.3.70The system shall throw an error if the user supplies a fluid temperature variable which does not exist
5.3.71The system shall throw an error if the user requests the currently unimplemented porous flow scalar quantity advection
5.3.72The system shall throw an error if the user supplies vector and scalar thermal conductivities together
5.3.73The system shall throw an error if the user supplies vector thermal conductivity with non-porous treatment
5.3.74The system shall throw an error if the user defines a non-existing block for thermal conductivity
5.3.75The system shall throw an error if there is a mismatch in the number of friction coefficients and the number of friction types
5.3.76
The system shall throw an error if
1. an initial condition is provided for velocities when they are defined outside of the action
2. an initial condition is provided for pressure when it is defined outside of the action
3. an initial condition is provided for temperature when it is defined outside of the action
5.3.77The system shall be able to run the NSFVAction if no block-restriction is defined on a mesh which doesn't have a default block.
5.3.78The system shall throw an error if the block-restrictions of the external variable and the action are different.
5.3.79
The system shall be able to perform a variety of limiting schemes when solving fluid flow equations. These schemes include
1. second-order upwind
2. van Leer
3. min-mod
4. QUICK
5.3.80The system should be able to solve the Navier-Stokes equations with block-restricted variables using the SIMPLE algorithm.
5.3.81The system shall be able to block-restrict all variables in a heated channel simulation with passive scalar advection.
5.3.82The system shall be able to reproduce benchmark results for a Rayleigh number of 1e3 using a finite volume discretization.
5.3.83The system shall be able to reproduce benchmark results for a Rayleigh number of 1e4 using a finite volume discretization.
5.3.84The system shall be able to reproduce benchmark results for a Rayleigh number of 1e5 using a finite volume discretization.
5.3.85The system shall be able to reproduce benchmark results for a Rayleigh number of 1e6 using a finite volume discretization.
5.3.86The system shall be able to reproduce benchmark results for a Rayleigh number of 1e6 using the INSFV actions.
5.3.87The system shall be able to model natural convection using a weakly compressible implementation.
5.3.88The system shall be able to model transient natural convection with a low Rayleigh number using a weakly compressible implementation.
5.3.89The system shall be able to solve incompressible Navier-Stokes channel flow with no-slip boundary conditions on the wall in an axisymmetric coordinate system using an average interpolation scheme for the velocity.
5.3.90The system shall be able to solve incompressible Navier-Stokes channel flow with no-slip boundary conditions on the wall in an axisymmetric coordinate system using a Rhie-Chow interpolation scheme for the velocity.
5.3.91The system shall be able to solve incompressible Navier-Stokes channel flow with free-slip boundary conditions on the wall in an axisymmetric coordinate system using a Rhie-Chow interpolation scheme for the velocity.
5.3.92The system shall be able to solve a diverging channel problem in cylindrical coordinates with no slip boundary conditions.
5.3.93The system shall be able to solve a straight channel problem in cylindrical coordinates using triangular elements with no slip boundary conditions.
5.3.94The system shall be able to solve a straight channel problem in cylindrical coordinates using triangular elements with no slip boundary conditions and the NSFV action syntax.
5.3.95The system shall be able to solve a diverging channel problem in cylindrical coordinates with free slip boundary conditions.
5.3.96The system shall conserve mass when solving a Cartesian channel flow problem with one symmetry boundary condition and one no-slip wall boundary condition.
5.3.97The system shall be able to model linear volumetric friction in a channel.
5.3.98The system shall be able to model quadratic volumetric friction in a channel.
5.3.99The system shall be able to model quadratic volumetric friction with an exponential friction factor correlation in a channel.
5.3.100The system shall be able to build the correct Jacobian using exponential correlation for the friction factor.
5.3.101The system shall be able to model linear volumetric friction in a channel using a custom action syntax.
5.3.102The system shall be able to model quadratic volumetric friction in a channel using a custom action syntax.
5.3.103The system shall be able to model quadratic volumetric friction with an exponential friction factor correlation in a channel a custom action syntax.
5.3.104The system shall be able to solve the steady-state Navier-Stokes problem in a 2D channel using the SIMPLE algorithm using the linear finite volume system.
5.3.105The system shall be able to solve the steady-state Navier-Stokes problem in a 3D channel using the SIMPLE algorithm using the linear finite volume system.
5.3.106The system shall give back the correct results for a channel flow with the SIMPLE algorithm using nonlinear system assembly.
5.3.107The system shall give back the correct results for a channel flow with the SIMPLE algorithm using linear system assembly.
5.3.108The system shall give back the correct results for a channel flow with the SIMPLE algorithm using linear system assembly and the Physics shorthand syntax.
5.3.109The system shall be able to solve the steady-state Navier-Stokes problem in a 2D channel using the SIMPLE algorithm with separating the momentum components into different systems.
5.3.110The system shall be able to solve the steady-state Navier-Stokes problem in a 2D cylindrical channel using the SIMPLE algorithm.
5.3.111The system shall be able to solve the steady-state Navier-Stokes problem in a 2D cylindrical channel with slip boundary conditions using the SIMPLE algorithm.
5.3.112The system shall be able to solve the steady-state Navier-Stokes problem in a 2D cylindrical channel with homogenized friction using the SIMPLE algorithm.
5.3.113The system shall be able to solve the steady-state Navier-Stokes problem together with the energy equation in a 2D channel using the SIMPLE algorithm.
5.3.114The system shall be able to produce a constant solution for temperature in a symmetric 2D channel flow without heat sources using the SIMPLE algorithm.
5.3.115The system shall be able to solve the steady-state Navier-Stokes problem together with the energy equation in a symmetric 2D channel using the SIMPLE algorithm.
5.3.116The system shall be able to solve the steady-state Navier-Stokes problem together with scalar transport equations in a 2D channel using the SIMPLE algorithm.
5.3.117The system shall conserve of a passive scalar while solving using the SIMPLE algorithm.
5.3.118The system shall be able to produce a constant solution for scalar values in a symmetric 2D channel flow without sources using the SIMPLE algorithm.
5.3.119The system shall be able to solve the steady-state Navier-Stokes problem in a 3D channel using the SIMPLE algorithm with separating the momentum components into different systems.
5.3.120The system shall be able to solve the steady-state Navier-Stokes problem together with the energy equation in a 3D channel using the SIMPLE algorithm.
5.3.121The system shall be able to solve the steady-state Navier-Stokes problem together with scalar transport equations in a 3D channel using the SIMPLE algorithm.
5.3.122The system shall be able to solve the steady-state Navier-Stokes problem on nonorthogonal meshes using the SIMPLE algorithm.
5.3.123The system shall be able to model free-slip conditions in a 1D channel; specifically the tangential velocity shall have a uniform value of unity and the pressure shall not change.
5.3.124The system shall be able to model free-slip conditions in a channel; specifically the tangential velocity shall have a uniform value of unity, the normal velocity shall have a uniform value of zero, and the pressure shall not change.
5.3.125The system shall be able to model free-slip conditions in a channel using the NSFV action.
5.3.126The system shall be able to initialize variables added by NSFV action from the mesh file.
5.3.127The system shall be able to restart a transient calculation from a steady calculation and immediately register the steady solution as a solution to the transient problem.
5.3.128
The system shall be able to model no-slip conditions in a channel; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile. This problem shall be solvable using a
1. monolithic preconditioning technique
2. field-split preconditioning technique
5.3.129The system shall be able to model no-slip conditions in a channel using incompressible Navier Stokes action.
5.3.130The system shall be able to transport arbitrary scalar field variables in a fluid flow field.
5.3.131The system shall be able to create and transport scalar field variables in a fluid flow field using the NSFVAction syntax.
5.3.132The system shall be able to create and transport scalar field variables in a fluid flow field using the Navier Stokes weakly-compressible Physics syntax.
5.3.133The system shall be able to transport an externally-generated scalar field variables in a fluid flow field using the NSFVAction syntax.
5.3.134The system shall be able to use flux boundary conditions for the momentum and match results produced by using flux kernels.
5.3.135The system shall be able to extrapolate a pressure value at a fully developed outflow boundary and use a mean pressure approach to eliminate the nullspace for the pressure.
5.3.136The system shall be able to model the effect of Reynolds-averaged parameters on the momentum and passive scalar advection equations using a mixing length model
5.3.137The system shall be able to model the effect of Reynolds-averaged parameters on the momentum and passive scalar advection equations using a mixing length model and the NSFVAction syntax
5.3.138The system shall be able to model the effect of Reynolds-averaged parameters on the momentum and passive scalar advection equations using a mixing length model and the Physics syntax
5.3.139The system shall be able to model the effect of Reynolds-averaged parameters on the momentum and passive scalar advection equations using a mixing length model and show a perfect Jacobian
5.3.140The system shall be able to model ambient volumetric convection in a channel.
5.3.141The system shall be able to model ambient volumetric convection in a channel using the Navier-Stokes action.
5.3.142
The system shall be able to run incompressible Navier-Stokes channel-flow simulations with
1. two-dimensional triangular elements
2. three-dimensional tetrahedral elements
5.3.143The system shall be able to model free-slip conditions in a 3D square channel; specifically the tangential velocity shall have a uniform value of unity and the pressure shall not change.
5.3.144The system shall be able to compute gradients, when extrapolated boundary face values make the gradient computation singular, by catching the singularity error and re-running without doing boundary face value extrapolation.
5.3.145The system shall show a monotone pressure profile in the presence of discontinuous body forces, in this case a transient from free flow to linear drag.
5.3.146The system shall be able to simulate a transient flow in a channel with enthalpy as an advected quantity.
5.3.147The system shall be able to simulate a transient flow in a channel with enthalpy as an advected quantity using the NSFV action syntax.
5.3.148The system shall be able to simulate a transient flow in a channel with enthalpy as an advected quantity using the Navier Stokes weakly compressible Physics syntax.
5.3.149The system shall be able to constrain the pressure nullspace by imposing an average pressure value on a boundary.
5.3.150The system shall report an error if a user specifies a residual object for the momentum equation that is not a incompressible Navier-Stokes finite volume momentum residual object.
5.3.151The system shall report an error if the interpolation object has block restriction different from the nonlinear flow variables.
5.3.152The system shall be able to solve a problem with surface heat transfer between a fluid and a solid block using FVConvectionCorrelationInterface.
5.3.153The system shall be able to solve Jeffery-Hamel flow using the finite volume method in a 2D wedge and compare to the analytical solution with Dirichlet boundary conditions
5.3.154The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity using the SIMPLE algorithm and the linear finite volume system.
5.3.155The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity using the SIMPLE algorithm, the linear finite volume system and the shorthand Physics syntax.
5.3.156The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity using the SIMPLE algorithm.
5.3.157The system shall be able to solve the incompressible Navier-Stokes equations together with the energy equation in a lid-driven cavity using the SIMPLE algorithm.
5.3.158The system shall be able to solve the incompressible buoyant Navier-Stokes equations with the Boussinesq approximation in a lid-driven cavity using the SIMPLE algorithm.
5.3.159The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity using the finite volume method.
5.3.160The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity using the finite volume Navier-Stokes action.
5.3.161The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity using the finite volume Navier-Stokes action and an average pressure pin applied at the end of every time step.
5.3.162The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity using the finite volume Navier-Stokes action with an approximate computation of the Rhie-Chow coefficients.
5.3.163The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity using the finite volume method on a (zero-)displaced mesh.
5.3.164The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity by fixing the point value of the pressure at a certain coordinate.
5.3.165The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity by fixing the point value of the pressure at a certain coordinate using the NSFV action syntax.
5.3.166The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity by fixing the point value of the pressure at a certain coordinate using the Navier-Stokes finite volume action syntax, with a post-treatment of the pressure rather than a Lagrange multiplier-based constraint.
5.3.167The system shall throw an error if the user requests integral value pressure pinning while specifying a point for the pin.
5.3.168The system shall be able to solve an incompressible Navier-Stokes problem with dirichlet boundary conditions for all the normal components of velocity, using the finite volume method, and have a nonsingular system matrix.
5.3.169The system shall be able to compute a perfect Jacobian when solving a lid-driven incompressible Navier-Stokes problem with the finite volume method.
5.3.170The system shall be able to transport scalar quantities using the simultaneously calculated velocity field from the incompressible Navier Stokes equations.
5.3.171The system shall be able to get the same result as the enthalpy transport example using the NSFVAction to set up the run.
5.3.172The system shall be able to get the same result as the enthalpy transport example using the Physics syntax to set up the run.
5.3.173The system shall be able to transport scalar quantities using the simultaneously calculated velocity field from the transient incompressible Navier Stokes equations.
5.3.174The system shall yield a quiescent fluid in an axisymmetric coordinate system with a gravitational force applied and Rhie-Chow interpolation used for the velocity field.
5.3.175The system shall compute an accurate Jacobian when a scaling factor is applied to a scalar variable.
5.3.176
The system shall be able to compute the turbulent viscosity based on the capped mixing length model and store it in a variable when performing
1. transient simulations
2. steady simulations
5.3.177
The system shall be able to calculate the material property comprising the total turbulent viscosity, based on the capped mixing length model when performing
1. transient simulations
2. steady simulations
3. steady simulations with action syntax
4. steady simulations with functions describing the boundary layer parameter, and
5. steady simulations with functions describing the boundary layer parameter with action syntax.
5.3.178The system shall exhibit second order convergence for all variables in a Cartesian, no-slip, channel-flow problem using a Rhie-Chow interpolation, including body forces, and two term boundary expansion for gradient and extrapolated boundary face value computation.
5.3.179The system shall exhibit second order convergence for all variables in a Cartesian, no-slip, channel-flow problem with a symmetry axis using a Rhie-Chow interpolation, including body forces, and two term boundary expansion for gradient and extrapolated boundary face value computation.
5.3.180The system shall exhibit second order convergence for all variables in an axisymmetric, no-slip, channel-flow problem using a Rhie-Chow interpolation, including body forces, and two term boundary expansion for gradient and extrapolated boundary face value computation.
5.3.181The system shall be able to solve a problem with channel-flow like boundary conditions in the coordinate system with an average interpolation for the velocity and demonstrate second order convergence in the velocity variables and first order convergence in the pressure variable.
5.3.182The system shall be able to solve a problem with channel-flow like boundary conditions in the coordinate system with a Rhie-Chow interpolation for the velocity and demonstrate second order convergence in the velocity and pressure variables.
5.3.183The system shall be able to solve the incompressible Navier-Stokes equations in an RZ coordinate system, including energy, using an average interpolation for the velocity, with a mix of Dirichlet and zero-gradient boundary conditions for each variable, and demonstrate second order convergence for each variable other than the pressure which shall demonstrate first order convergence.
5.3.184The system shall be able to solve the incompressible Navier-Stokes equations in an RZ coordinate system, including energy, using a RC interpolation for the velocity, with a mix of Dirichlet and zero-gradient boundary conditions for each variable, and demonstrate second order convergence for each variable.
5.3.185The system shall demonstrate global second order convergence for all variables on a rotated mesh when using an average interpolation for the velocity and a two term Taylor series expansion for face values on non-Dirichlet boundaries.
5.3.186The system shall demonstrate global second order convergence for all variables on a rotated mesh when using an RC interpolation for the velocity and a two term Taylor series expansion for face values on non-Dirichlet boundaries.
5.3.187The system shall demonstrate global second order convergence for velocity variables and first order convergence for the pressure variable on a rotated mesh when using an average interpolation for the velocity and a one term Taylor series expansion for face values on non-Dirichlet boundaries.
5.3.188The system shall demonstrate global second order convergence for all variables on a rotated mesh when using an RC interpolation for the velocity and a one term Taylor series expansion for face values on non-Dirichlet boundaries.
5.3.189The system shall be able to solve the incompressible Navier-Stokes equations in one dimension with prescribed inlet velocity and outlet pressure and implicit zero gradient boundary conditions elsewhere, and demonstrate second order convergence in both velocity and pressure when using an average interpolation scheme for the velocity.
5.3.190The system shall be able to solve the incompressible Navier-Stokes equations in two dimensions with prescribed inlet velocity and outlet pressure, free slip along the walls, and implicit zero gradient boundary conditions elsewhere, and demonstrate second order convergence in both velocity and pressure when using an average interpolation scheme for the velocity.
5.3.191The system shall be able to solve the incompressible Navier-Stokes equations in two dimensions with prescribed inlet velocity and outlet pressure, free slip along the walls, and implicit zero gradient boundary conditions elsewhere, and demonstrate second order convergence in both velocity and pressure when using a Rhie-Chow interpolation scheme for the velocity.
5.3.192The system shall be able to solve the incompressible Navier-Stokes equations in two dimensions with prescribed inlet velocity and outlet pressure, free slip along the walls, and implicit zero gradient boundary conditions elsewhere, and demonstrate second order convergence in both velocity and pressure when using an approximate Rhie-Chow interpolation scheme for the velocity.
5.3.193The system shall demonstrate global second order convergence for all variables when using an RC interpolation for the velocity and a two term Taylor series expansion for face values on non-Dirichlet boundaries.
5.3.194The system shall demonstrate global second order convergence for all variables when using an RC interpolation for the velocity and a one term Taylor series expansion for face values on non-Dirichlet boundaries.
5.3.195The system shall be able to solve the incompressible Navier-Stokes equations, including energy, using an average interpolation for the velocity, with a mix of Dirichlet and zero-gradient boundary conditions for each variable, and demonstrate second order convergence for each variable.
5.3.196The system shall be able to solve the incompressible Navier-Stokes equations, including energy, using a Rhie-Chow interpolation for the velocity, with a mix of Dirichlet and zero-gradient boundary conditions for each variable, and demonstrate second order convergence for each variable.
5.3.197The system shall be able to solve the incompressible Navier-Stokes equations in 2D cylindrical coordinates, using a Rhie-Chow scheme, dirichlet boundary conditions for both variables, and demonstrate second order convergence for the velocity and pressure.
5.3.198The system shall be able to subtract the time rate of change of mesh displacements from the velocity field in order to compute the correct rate of momentum advection in an Abritrary Lagrangian Eulerian simulation.
5.3.199The system shall be able to solve the Navier Stokes equations using the SIMPLE algorithm and obtain second order spatial convergence for velocity and at least first order spatial convergence for pressure on an orthogonal grid.
5.3.200The system shall be able to solve the Navier Stokes equations using the SIMPLE algorithm and obtain second order spatial convergence for velocity and at least first order spatial convergence for pressure on a nonorthogonal grid.
5.3.201The system shall be able to solve the incompressible Navier-Stokes equations on triangular meshes, using a Rhie-Chow scheme and skewness-correction, Dirichlet boundary conditions for the velocity, and demonstrate second order convergence for the velocity and first order convergence for pressure.
5.3.202The system shall be able to solve the incompressible Navier-Stokes equations on triangular meshes, using an approximated Rhie-Chow scheme and skewness-correction, Dirichlet boundary conditions for the velocity, and demonstrate second order convergence for the velocity and first order convergence for pressure.
5.3.203The system shall be able to solve the incompressible Navier-Stokes equations on triangular meshes using skewness-correction and the NSFV action syntax.
5.3.204The system shall be able to automatically expand its ghosting pattern when skew-corrected face gradients are involved in the simulation.
5.3.205The system shall be able to solve the incompressible Navier-Stokes equations using a Rhie-Chow interpolation scheme and produce second order convergence for all variables.
5.3.206The system shall be able to solve the incompressible Navier-Stokes equations using an approximate Rhie-Chow interpolation scheme and produce second order convergence for all variables.
5.3.207The system shall be able to transport passive scalar quantities in an application different from the fluid flow.
5.3.208The system shall be able to transport passive scalar quantities in an application different from the fluid flow, using the NSFV simplified action syntax.
5.3.209The system shall report an error if both average and Rhie Chow velocity interpolation parameters are provided.
5.3.210The system shall report an error if Rhie Chow velocity interpolation is requested but the a coefficients are not provided and cannot be computed from the momemtum equation.
5.3.211
The system shall report an error if the Rhie Chow a coefficients are provided in a manner inconsistent with the dimension of the problem
1. for example with the Y coefficient but not the X one
2. for example with the Z coefficient but not the X one
3. for example if the x coefficient is missing in a >1D problem
4. for example if the y coefficient is missing in a >2D problem
5.3.212The system shall be able to use pressure inlet and outlet boundary conditions to compute open, chimney-type natural circulation problems.
5.3.213The system shall be able to correctly predict the maximum fluid and solid temperatures.
5.3.214The system shall be able to group sidesets logically connected to neighboring subdomains but contacting the fluid subdomains into the fluid flow boundary set, allowing for accurate solution of the fluid flow equations.
5.3.215The system shall maintain the correct flow and maximum speed (no oscillations) with one forward rotating pump.
5.3.216The system shall maintain the correct flow and maximum speed (no oscillations) with one backwards rotating pump.
5.3.217The system shall maintain the correct flow and maximum speed (no oscillations) with two competing pumps.
5.3.218The system shall maintain the correct flow and maximum speed (no oscillations) with a body force correction to Rhie-Chow computed by examining the difference between two different kinds of body force interpolation to faces.
5.3.219The system shall maintain the correct flow and maximum speed (no oscillations) with a body force correction to Rhie-Chow computed by checking for a non-zero body force face interpolation.
5.3.220The system shall be able to solve solidification without advection using the finite volume method.
5.3.221The system shall be able to solve solidification with advection in a pipe using the finite volume method.
5.3.222The system shall be able to conserve energy during solidification.
5.3.223The system must be able to obtain the correct temperature profile over the center line in a solidification problem.
5.3.224
The system shall be able to solve fluid flow problems with k-epsilon turbulence model for a backward facing step,
1. and reach converged results with segregated solvers,
2. and pass debugging checks.
5.3.225
The system shall be able to solve fluid flow problems with k-epsilon turbulence model for a standard channel,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.226The system shall be able to show a perfect Jacobian when solving the flow equations, the energy conservation equation and the k-epsilon turbulence equations.
5.3.227
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with standard wall functions,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.228
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with standard wall functions using bulk wall treatment for the turbulent viscosity,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.229
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with the time scale limiter deactivated,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.230
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with non-equilibrium wall functions,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.231
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with linear wall functions,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.232
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with incremental wall functions,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.233
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with no wall treatment and the first cell in the logarithmic regime,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.234
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with no wall treatment,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.235
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with non-equilibrium bulk treatment,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.236
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with transported energy,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.237
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with temperature wall functions,
1. and reach converged results with segregated solvers.
2. and pass debugging checks,
3. using the Navier Stokes Physics shorthand syntax.
5.3.238
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with linearized temperature wall functions,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.239
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with temperature wall functions in the log-layer regime,
1. and reach converged results with segregated solvers.
2. and pass debugging checks.
5.3.240
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with capped production,
1. and reach converged results with segregated solvers.
2. and pass debugging checks,
3. using the Navier Stokes Physics shorthand syntax.
5.3.241
The system shall be able to run fluid flow with k-epsilon turbulence model in enclosures with standard wall functions,
1. and reach converged results with monolithic nonlinear solvers.
2. and pass debugging checks,
3. using the Navier Stokes Physics shorthand syntax.
5.3.242The system shall be able to solve for flow in a 3D channel while not caching cell gradients.
5.3.243The system shall be able to solve for flow in a 3D channel while caching cell gradients.
5.3.244The system shall be able to compute the turbulent viscosity based on the capped mixing length model.
5.3.245The system shall be able to model flow around a bend with the porous incompressible Navier Stokes equations using a finite volume discretization and an Ergun drag correlation.
5.3.246The system shall be able to model a simple friction-based volumetric flow diode
5.3.247
The system shall be able to control a simple friction-based volumetric flow diode
1. with a simple time-based turn-on criterion,
2. with a simple pressure-based turn-on criterion, and
3. with a simple flow-based turn-on criterion.
5.3.248The system shall be able to compute phase-averaged quantities using a material.
5.3.249The system shall be able to run a flow simulation using both scalar advection and K-epsilon turbulence.
5.3.250
The system shall be able to restart all the flow variable in the shorthand Navier Stokes Physics-syntax
1. with the default initial conditions,
2. with user-defined initial conditions,
3. when performing a regular checkpoint restart, but still obeying the user-defined initial conditions,
4. when performing manual restart from a mesh file, ignoring the default initial conditions.
5.3.251The system should be able to solve the porous Navier-Stokes equations with block-restricted variables using the SIMPLE algorithm.
5.3.252The system shall be able to simulate a porous medium flow with block restriction.
5.3.253
The system shall be able to solve for fluid energy diffusion, advection and convection with the solid phase in a 2D channel
1. with a Cartesian geometry, only modeling the fluid phase,
2. with a Cartesian geometry, only modeling the fluid phase, with objects running on the displaced mesh,
3. in rz geometry,
4. with an effective diffusion coefficient,
5. with an effective diffusion coefficient and NSFVAction syntax,
6. with an effective diffusion coefficient and dirichlet inlet boundary condition for fluid temperature,
5.3.254The system shall be able to solve for fluid energy diffusion, advection and convection with the solid phase in a 2D channel, modeling both fluid and solid temperature.
5.3.255The system shall be able to solve for fluid energy diffusion, advection and convection with the solid phase in a 2D channel using the NSFV action syntax.
5.3.256The system shall be able to solve for fluid energy diffusion, advection and convection with the solid phase in a 2D channel using the WCNSFV Physics syntax.
5.3.257The system shall be able to solve for fluid energy diffusion, advection and block-restricted convection.
5.3.258The system shall be able to solve for fluid energy diffusion, advection and block-restricted convection with the solid phase in a 2D channel using the NSFV action syntax.
5.3.259The system shall be able to solve transient relaxations with fluid energy diffusion, advection and convection with the solid phase in a 2D channel, modeling both fluid and solid temperature.
5.3.260The system shall be able to run transient simulations with fluid energy diffusion, advection and convection with the solid phase in a 2D channel, modeling both fluid and solid temperature using the NSFV action syntax.
5.3.261The system shall be able to solve for fluid energy diffusion, advection and convection with the solid phase in a 2D channel with a Boussinesq approximation for the influence of temperature on density.
5.3.262The system shall be able to solve for fluid energy diffusion, advection and convection with the solid phase in a 2D channel with a Boussinesq approximation using the NSFV action syntax.
5.3.263The system shall be able to model a smooth porosity gradient in a 2D channel.
5.3.264The system shall be able to model a discontinuous porosity jump in a 1D channel with average interpolation of velocity and advected quantity.
5.3.265The system shall be able to model a discontinuous porosity jump in a 1D channel with average interpolation of velocity and upwinding of the advected quantity.
5.3.266The system shall be able to model a discontinuous porosity jump in a 1D channel with Rhie Chow interpolation of velocity and averaging of the advected quantity.
5.3.267The system shall be able to model a discontinuous porosity jump in a 1D channel with Rhie Chow interpolation of velocity and upwinding of the advected quantity.
5.3.268The system shall be able to model a discontinuous porosity jump in a 2D channel.
5.3.269The system shall exhibit a monotonic pressure profile when a number of reconstructions (essentially smoothing) are applied to the porosity.
5.3.270The system shall report an error if the user attempts to create a reconstructed field from a variable.
5.3.271
The system shall be able to automatically impose a change in pressure in line with the Bernoulli equation at porosity jumps in
1. one dimension, using upwinding for advection and a Rhie-Chow interpolation of velocity
2. one dimension, using central differencing for advection and velocity
3. two dimensions, using upwinding for advection and a Rhie-Chow interpolation of velocity
4. two dimensions, using central differencing for advection and velocity
5. two dimensions using compact action syntax, upwinding for advection, and a Rhie-Chow interpolation of velocity
5.3.272The system shall be able to mark what elements have porosity jump faces using an auxiliary kernel.
5.3.273The system shall be able to recognize discontinuities in parsed functions at subdomain changes and produce the correct Bernoulli pressure drop as a consequence.
5.3.274The system shall be able to use two-term expansion to compute the face pressure on the downstream side on a porosity jump face.
5.3.275The system shall be able to recognize discontinuities in functor material properties and produce the correct Bernoulli pressure drop as a consequence.
5.3.276The system shall be able to solve the steady-state porous Navier-Stokes problem in a 2D channel using the SIMPLE algorithm.
5.3.277The system shall be able to solve the steady-state porous Navier-Stokes problem in a 2D channel with slip and symmetry boundary conditions using the SIMPLE algorithm.
5.3.278The system shall be able to solve the steady-state porous Navier-Stokes problem in a 2D channel with friction caused by porous media using the SIMPLE algorithm.
5.3.279The system shall be able to solve the steady-state porous Navier-Stokes problem coupled with both solid and fluid energy equations in a 2D channel using the SIMPLE algorithm.
5.3.280The system shall be able to model free-slip conditions in a porous media channel; specifically the tangential velocity shall have a uniform value of unity, the normal velocity shall have a uniform value of zero, and the pressure shall not change.
5.3.281The system should be able to run a porous medium channel computation using the NSFV action.
5.3.282The system shall be able to model free-slip conditions in a porous media cylindrical channel; specifically the tangential velocity shall have a uniform value of unity, the normal velocity shall have a uniform value of zero, and the pressure shall not change.
5.3.283The system shall be able to model no-slip conditions in a porous media channel; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.
5.3.284The system shall be able to model no-slip conditions in a porous media channel with flow driven by a pressure differential; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.
5.3.285The system shall be able to model no-slip conditions in a porous media channel with flow driven by a pressure differential using the NSFVAction syntax.
5.3.286The system shall be able to model no-slip conditions in a porous media channel with a set mean pressure; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.
5.3.287The system shall be able to model no-slip conditions in a porous media channel using the NSFV action and a set mean pressure.
5.3.288The system shall be able to model no-slip conditions in a porous media channel using an average interpolation for velocity; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.
5.3.289The system shall be able to model no-slip conditions in a porous media channel with a porosity of 1; specifically, it should match a regular INSFV simulation results.
5.3.290The system shall be able to model no-slip conditions in a porous media channel with reflective boundary conditions on one side; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.
5.3.291The system shall be able to model no-slip conditions in a symmetric porous media channel using the NSFV action.
5.3.292The system shall be able to model no-slip conditions in a cylindrical porous media channel with reflective boundary conditions on one side; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.
5.3.293The system shall be able to model porous flow with volumetric friction, using the Darcy and Forchheimer friction models with no slip boundary conditions on the channel walls.
5.3.294The system shall be able to model porous flow with volumetric friction, using the Darcy and Forchheimer friction models and the NSFV action input syntax.
5.3.295The system shall be able to model porous flow with volumetric friction, using the Darcy and Forchheimer friction models with free slip boundary conditions on the channel walls.
5.3.296The system shall be able to model porous flow with volumetric friction, using the Darcy and Forchheimer friction models with free slip boundary conditions using the NSFV action.
5.3.297The system shall be able to model porous flow with block-restricted volumetric friction.
5.3.298The system shall be able to model porous flow with block-restricted volumetric friction using the NSFV action.
5.3.299The system shall be able to solve the porous flow equations in RZ geometry using an integration by parts of the pressure term.
5.3.300The system shall be able to model porous flow with volumetric friction with a friction factor that depends linearly on velocity.
5.3.301The system shall be able to compute the speed, eg the norm of the interstitial velocity, in a porous media incompressible flow problem.
5.3.302The system shall be able to compute the speed in a porous media incompressible flow problem from within the NSFVAction.
5.3.303The system shall be able to solve the incompressible porous flow Navier-Stokes equations using a Rhie-Chow interpolation scheme in a 1D channel with a continuously varying porosity and produce second order convergence for all variables.
5.3.304The system shall be able to solve the incompressible porous flow Navier-Stokes equations using a Rhie-Chow interpolation scheme in a 2D channel with a continuously varying porosity and produce second order convergence for all variables.
5.3.305The system shall show second order convergence for all variables when the porosity is interpolated and reconstructed multiple times and when a correction is applied to the pressure interpolation.
5.3.306The system shall show second order convergence for all variables when using porosity reconstructions and friction corrections together with NSFV action syntax.
5.3.307The system shall show a perfect Jacobian when the porosity is interpolated and reconstructed multiple times and when a correction is applied to the pressure interpolation.
5.3.308The system shall be able to solve the incompressible porous flow Navier-Stokes equations in a 1D channel using a Rhie-Chow interpolation scheme and produce second order convergence for all variables.
5.3.309The system shall be able to solve the incompressible porous flow Euler equations in a 1D channel using a Rhie-Chow interpolation scheme for velocity and upwind interpolation for advected quantities and produce first order convergence for pressure and 1.5 order for velocity.
5.3.310The system shall be able to solve the incompressible porous flow Euler equations in a 1D channel using a Rhie-Chow interpolation scheme for velocity and upwind interpolation for advected quantities and produce first order convergence for pressure and 1.5 order for velocity with using flux boundary conditions for both variables.
5.3.311The system shall be able to solve the incompressible porous flow Navier-Stokes equations in a 2D channel using a Rhie-Chow interpolation scheme and produce second order convergence for all variables.
5.3.312
The system shall be able to connect subdomains within a diffusion problem
1. using interface kernels,
2. using a functor material automatically defined by the physics syntax to combine the thermal diffusivities.
5.3.313The system shall be able to model a heated porous solid phase.
5.3.314The system shall be able to model heated flow in a porous medium with convection between phases.
5.3.315The system shall be able to pull all nonlocal Rhie-Chow coefficient data from the owning process to a given process such that user code may access Rhie-Chow data at arbitrary locations.
5.3.316The system shall be able to model a momentum inlet condition based on mass flow rate for porous weakly-compressible flow.
5.3.317The system shall be able to solve transient relaxations within the weakly compressible approximation, with fluid energy diffusion, advection and convection with the solid phase in a 2D channel, modeling both fluid and solid temperature.
5.3.318The system shall be able to solve transient relaxations within the weakly compressible approximation, using the time derivative of the specific enthalpy for the time derivative.
5.3.319The system shall be able to solve weakly compressible transient problems with the NSFV action syntax.
5.3.320The system shall be able to solve weakly compressible transient problems with flow physics syntax.
5.3.321The system shall be able to solve transient relaxations within the weakly compressible approximation, with fluid energy diffusion, advection and convection with the solid phase in a 2D channel, modeling both fluid and solid temperature and show a perfect Jacobian.
5.3.322The system shall be able to track interfacial area in steady-state, two-phase problems dominated by pressure driven growth.
5.3.323The system shall be able to track interfacial area in steady-state, two-phase problems dominated by turbulence driven growth.
5.3.324The system shall be able to track interfacial area in transient, two-phase problems dominated by pressure driven growth.
5.3.325The system shall give the correct solution for two phase flow in a lid-driven cavity with the mixture model.
5.3.326The system shall give the correct solution for two phase flow in a lid-driven cavity with a mixture model using Navier Stokes physics syntax.
5.3.327The system shall be able solve two phase Rayleigh Bernard convection with the mixture model. The problem has multiple solutions so RunApp is used for testing.
5.3.328The system shall be able to solve two-phase flow in a channel using the mixture drift flux model.
5.3.329The system shall be able to solve two-phase flow in a channel using a mixture drift flux model using Navier Stokes physics syntax.
5.3.330The system shall be able to solve transient two-phase flow in a channel using the mixture drift flux model.
5.3.331The system shall be able to solve two-phase flow in a channel using the mixture advection-slip model.
5.3.332The system shall be able to solve two-phase flow in a channel using amixture advection-slip model using Navier Stokes physics syntax.
5.3.333The system shall report an error if it tries to require quadrature point computations on variable that were deliberately set not to perform quadrature point calculations.
5.3.334The system shall be able to simulate flow in a channel using finite volume methods with an optimization that avoids computations / re-initialization on quadrature points.
5.3.335
The system shall be able to use Dirichlet boundary conditions for specifying inlet conditions in a weakly compressible fluid flow simulation
1. using a velocity postprocessor
2. using a mass flow rate postprocessor
5.3.336
The system shall be able to use flux boundary conditions for specifying inlet conditions in a weakly compressible fluid flow simulation
1. specifying the flux values directly though the flow rates
2. using the inlet velocity to compute the fluxes
3. using the inlet velocity to compute the fluxes with the NSFV action syntax
4. using the inlet mass flow rate to compute the fluxes
5. using the inlet mass flow rate to compute the fluxes with NSFVAction syntax
6. and correctly handle inflow and outflow boundaries
5.3.337
The system shall report an error if
1. the inlet velocity cannot be determined from the mass flow rate due to insufficient parameters
2. the inlet temperature cannot be determined because neither the temperature or the energy flow rate have been provided
3. the inlet temperature cannot be determined from the energy flow rate due to missing fluid velocity information
4. the inlet temperature cannot be determined from the velocity due to insufficient parameters
5. the inlet temperature cannot be determined from the mass flow rate due to insufficient parameters
5.3.338
The system shall return a warning if
1. redundant information is provided for setting the inlet velocity
2. redundant information is provided for setting the inlet temperature
5.3.339
The system shall report an error if
1. only the mass flow rate is provided to compute the inlet mass flux
2. only the mass flow rate is provided to compute the inlet momentum flux
3. the inlet energy cannot be determined because neither the temperature or the energy flow rate have been provided
4. the inlet energy flux cannot be determined because the temperature was provided but neither the inlet mass flow rate or velocity were provided
5. neither the scalar quantity flux or boundary values are specified when attempting to compute the scalar quantity flux
6. only the mass flow rate is provided to compute the inlet scalar quantity flux
7. only the scalar quantity boundary value is provided to compute the inlet scalar quantity flux
5.3.340
The system shall return a warning if
1. redundant information is provided for setting the inlet mass flux
2. redundant information is provided for setting the inlet momentum flux
3. redundant information is provided for setting the inlet energy flux
4. redundant information is provided for setting the inlet scalar flux
5.3.341The system shall throw an error when the mass flux boundary (with velocity) is defined on an internal face without explicitly specifying the direction of the flow.
5.3.342The system shall throw an error when the mass flux boundary is defined with an incorrect direction vector.
5.3.343The system shall throw an error when the momentum flux boundary is defined on an internal face without explicitly specifying the direction of the flow.
5.3.344The system shall throw an error when the momentum flux boundary is defined with an incorrect direction vector.
5.3.345The system shall throw an error when the energy flux boundary (with velocity) is defined on an internal face without explicitly specifying the direction of the flow.
5.3.346The system shall throw an error when the energy flux boundary is defined with an incorrect direction vector.
5.3.347The system shall throw an error when the passive scalar flux boundary is defined on an internal face without explicitly specifying the direction of the flow.
5.3.348The system shall throw an error when the passive scalar flux boundary (with is defined with an incorrect direction vector.
5.3.349
The system shall be able to model a momentum inlet condition based on mass flow rate for porous weakly-compressible flow
1. with the incoming flow being parallel to the surface normal of the inlet;
2. with the incoming flow direction spanning an angle with the normal surface of the inlet;
5.3.350
The system shall be able to model a momentum inlet condition based on velocity magnitude for porous weakly-compressible flow
1. with the incoming flow being parallel to the surface normal of the inlet;
2. with the incoming flow direction spanning an angle with the normal surface of the inlet;
5.3.351The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation.
5.3.352The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation using the NSFV action syntax.
5.3.353The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation using the WCNSFV Physics syntax.
5.3.354The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation and a mixing length turbulence model.
5.3.355The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation and a mixing length turbulence model using the NSFV action syntax.
5.3.356The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation and show a perfect Jacobian.
5.3.357The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation and a mixing length turbulence model and show a perfect Jacobian.
5.3.358The system shall be able to use realistic fluid properties in a weakly compressible flow simulation
5.3.359The system shall be able to output grandeurs, derivatives and non-dimensional quantities from realistic functor fluid properties
5.3.360The system shall be able to neglect the derivatives with regards to nonlinear variables of the density first order time derivative.
5.3.361The system shall be able to use pressure inlet and outlet boundary conditions to compute open, chimney-type natural circulation problems using weakly compressible navier-stokes equations.

navier_stokes: Ics
5.4.1The system shall be able to set initial conditions for fluid flow variables.
5.4.2The system shall be able to set initial conditions for fluid flow variables using functions.
5.4.3The system shall be able to set intial conditions for porous flow variables.

navier_stokes: Postprocessors
5.5.1The system shall be able to compute mass and momentum flow rates at internal and external boundaries of a straight channel with a finite element incompressible Navier Stokes model.
5.5.2The system shall be able to compute mass and momentum flow rates at internal and external boundaries of a diverging channel with a finite element incompressible Navier Stokes model.
5.5.3The system shall be able to compute flow rates and prove mass, momentum and energy conservation at internal and external boundaries of a frictionless heated straight channel with a finite volume incompressible Navier Stokes model.
5.5.4
The system shall be able to compute flow rates and prove mass, momentum and energy conservation at internal and external boundaries of a frictionless heated diverging channel with a finite volume incompressible Navier Stokes model,
1. with a quadrilateral mesh in XY geometry, with mass flow measured using either a variable or material property,
2. with a quadrilateral mesh in RZ geometry,
3. with a triangular mesh in XY geometry,
4. with upwind interpolation of advected quantities,
5. with no-slip boundary conditions, for which momentum and energy will be dissipated at the wall.
6. with uniform refinement near an internal interface.
7. at the very beginning of the simulation, with the initialized velocities
5.5.5The system shall be able to compute flow rates and prove mass, momentum and energy conservation at internal and external boundaries of a frictionless heated straight channel with a finite volume porous media incompressible Navier Stokes model.
5.5.6
The system shall be able to compute flow rates and prove mass, momentum and energy conservation at internal and external boundaries of a frictionless heated diverging channel with a finite volume porous media incompressible Navier Stokes model,
1. with a quadrilateral mesh in XY geometry, with mass flow measured using either a variable or material property,
2. with a quadrilateral mesh in RZ geometry,
3. with upwind interpolation of advected quantities,
4. and with no-slip boundary conditions, for which momentum and energy will be dissipated at the wall.
5.5.7The system shall compute mass-flux weighted averages of quantities over boundaries.
5.5.8
The system shall report an error if
1. a volumetric flow rate is requested at the initialization of the simulation on a boundary internal to the flow domain when using finite volume and Rhie Chow interpolation, as this is not supported
5.5.9The system shall be able to compute the pressure drop in a straight channel with a finite element incompressible Navier Stokes model.
5.5.10
The system shall be able to compute the pressure drop in a diverging channel with a finite element incompressible Navier Stokes model
1. with a regular face pressure evaluation, and
2. with a pressure drop face evaluation weighted by the local velocity.
5.5.11
The system shall be able to compute the pressure drop in a frictionless heated straight channel with a finite volume incompressible Navier Stokes model
1. with a regular face pressure evaluation, and
2. with a pressure drop face evaluation weighted by the local velocity
5.5.12
The system shall be able to compute the pressure drop in a frictionless heated diverging channel with a finite volume incompressible Navier Stokes model,
1. with a quadrilateral mesh in XY geometry, with mass flow measured using either a variable or material property, and
2. with a quadrilateral mesh in RZ geometry.
5.5.13
The system shall report an error in a pressure drop calculation if
1. an upstream boundary for the pressure is not a boundary for the postprocessor,
2. a downstream boundary for the pressure is not a boundary for the postprocessor,
3. a boundary for the postprocessor is not part of either the upstream or downstream pressure evaluation,
4. a downstream boundary is also an upstream boundary for the pressure drop,
5. the weighting functor integrates to 0, and
6. a face interpolation rule is specified for a finite element pressure variable.
5.5.14The system shall be able to compute the Rayleigh number in a natural convection flow simulation

Usability Requirements

No requirements of this type exist for this application, beyond those of its dependencies.

Performance Requirements

No requirements of this type exist for this application, beyond those of its dependencies.

System Interfaces

No requirements of this type exist for this application, beyond those of its dependencies.

System Operations

Human System Integration Requirements

The Navier Stokes module is command line driven and conforms to all standard terminal behaviors. Specific human system interaction accommodations shall be a function of the end-user's terminal. MOOSE (and therefore the Navier Stokes module) does support optional coloring within the terminal's ability to display color, which may be disabled.

Maintainability

The latest working version (defined as the version that passes all tests in the current regression test suite) shall be publicly available at all times through the repository host provider.
Flaws identified in the system shall be reported and tracked in a ticket or issue based system. The technical lead will determine the severity and priority of all reported issues and assign resources at their discretion to resolve identified issues.
The software maintainers will entertain all proposed changes to the system in a timely manner (within two business days).
The core software in its entirety will be made available under the terms of a designated software license. These license terms are outlined in the LICENSE file alongside the Navier Stokes module source code. As a MOOSE physics module, the license for the Navier Stokes module is identical to that of the framework - that is, the LGPL version 2.1 license.

Reliability

The regression test suite will cover at least 77% of all lines of code within the Navier Stokes module at all times. Known regressions will be recorded and tracked (see Maintainability) to an independent and satisfactory resolution.

System Modes and States

MOOSE applications normally run in normal execution mode when an input file is supplied. However, there are a few other modes that can be triggered with various command line flags as indicated here:

Command Line Flag	Description of mode
`-i <input_file>`	Normal execution mode
`--split-mesh <splits>`	Read the mesh block splitting the mesh into two or more pieces for use in a subsequent run
`--use-split`	(implies -i flag) Execute the simulation but use pre-split mesh files instead of the mesh from the input file
`--yaml`	Output all object descriptions and available parameters in YAML format
`--json`	Output all object descriptions and available parameters in JSON format
`--syntax`	Output all registered syntax
`--registry`	Output all known objects and actions
`--registry-hit`	Output all known objects and actions in HIT format
`--mesh-only` (implies -i flag)	Run only the mesh related tasks and output the final mesh that would be used for the simulation
`--start-in-debugger <debugger>`	Start the simulation attached to the supplied debugger

note

The list of system-modes may not be extensive as the system is designed to be extendable to end-user applications. The complete list of command line options for applications can be obtained by running the executable with zero arguments. See the command line usage.

Physical Characteristics

The Navier Stokes module is software only with no associated physical media. See System Requirements for a description of the minimum required hardware necessary for running the Navier Stokes module.

Environmental Conditions

Not Applicable

System Security

MOOSE-based applications such as the Navier Stokes module have no requirements or special needs related to system security. The software is designed to run completely in user-space with no elevated privileges required nor recommended.

Information Management

The core framework and all modules in their entirety will be made publicly available on an appropriate repository hosting site. Day-to-day backups and security services will be provided by the hosting service. More information about MOOSE backups of the public repository on INL-hosted services can be found on the following page: GitHub Backups

Polices and Regulations

MOOSE-based applications must comply with all export control restrictions.

System Life Cycle Sustainment

MOOSE-based development follows various agile methods. The system is continuously built and deployed in a piecemeal fashion since objects within the system are more or less independent. Every new object requires a test, which in turn requires an associated requirement and design description. The Navier Stokes module development team follows the NQA-1 standards.

Packaging, Handling, Shipping and Transportation

No special requirements are needed for packaging or shipping any media containing MOOSE and Navier Stokes module source code. However, some MOOSE-based applications that use the Navier Stokes module may be export-controlled, in which case all export control restrictions must be adhered to when packaging and shipping media.

Verification

The regression test suite will employ several verification tests using comparison against known analytical solutions, the method of manufactured solutions, and convergence rate analysis.

Introduction
System Overview
References
Definitions and Acronyms
System Requirements
Verification