Add squirrel tests

.github/workflows/tests.yml

@@ -1,7 +1,7 @@
 # Adapted from the RACCOON project's test workflow:
 # https://github.com/hugary1995/raccoon/blob/devel/.github/workflows/tests.yml
 
-name: Tests
+name: Squirrel Continuous Integration (CI) Testing
 
 on:
   # Triggers the workflow on pushes to the master branch

.gitignore

@@ -152,3 +152,6 @@ squirrel.yaml
 
 # VSCode dotfolder
 .vscode
+
+# JIT files
+**/.jitcache/*

include/bcs/PostprocessorVelocityFunctionTemperatureInflowBC.h

@@ -0,0 +1,24 @@
+#pragma once
+
+#include "PostprocessorVelocityFunctionInflowBC.h"
+#include "JvarMapInterface.h"
+#include "DerivativeMaterialInterface.h"
+
+class PostprocessorVelocityFunctionTemperatureInflowBC
+    : public DerivativeMaterialInterface<JvarMapIntegratedBCInterface<PostprocessorVelocityFunctionInflowBC>>
+{
+public:
+  PostprocessorVelocityFunctionTemperatureInflowBC(const InputParameters & parameters);
+
+  static InputParameters validParams();
+
+protected:
+  virtual void initialSetup() override;
+  virtual Real computeQpResidual() override;
+  virtual Real computeQpJacobian() override;
+
+  const MaterialProperty<Real> & _rho;
+  const MaterialProperty<Real> & _d_rho_d_u;
+  const MaterialProperty<Real> & _cp;
+  const MaterialProperty<Real> & _d_cp_d_u;
+};

src/bcs/DiffusiveFluxBC.C

@@ -27,12 +27,12 @@ DiffusiveFluxBC::initialSetup()
 Real
 DiffusiveFluxBC::computeQpResidual()
 {
-  return _normals[_qp] * -_D[_qp] * _grad_u[_qp];
+  return -_test[_i][_qp] * _normals[_qp] * _D[_qp] * _grad_u[_qp];
 }
 
 Real
 DiffusiveFluxBC::computeQpJacobian()
 {
-  return -_normals[_qp] *
+  return -_test[_i][_qp] * _normals[_qp] *
          (_d_D_d_u[_qp] * _phi[_j][_qp] * _grad_u[_qp] + _D[_qp] * _grad_phi[_j][_qp]);
 }

src/bcs/PostprocessorTemperatureInflowBC.C

@@ -6,6 +6,8 @@ InputParameters
 PostprocessorTemperatureInflowBC::validParams()
 {
   InputParameters params = PostprocessorInflowBC::validParams();
+  params.addParam<MaterialPropertyName>("rho", "rho", "Density material property name");
+  params.addParam<MaterialPropertyName>("cp", "cp", "Specific heat material property name");
   return params;
 }
 

src/bcs/PostprocessorVelocityFunctionTemperatureInflowBC.C

@@ -0,0 +1,43 @@
+#include "PostprocessorVelocityFunctionTemperatureInflowBC.h"
+
+registerMooseObject("SquirrelApp", PostprocessorVelocityFunctionTemperatureInflowBC);
+
+InputParameters
+PostprocessorVelocityFunctionTemperatureInflowBC::validParams()
+{
+  InputParameters params = PostprocessorVelocityFunctionInflowBC::validParams();
+  params.addParam<MaterialPropertyName>("rho", "rho", "Density material property name");
+  params.addParam<MaterialPropertyName>("cp", "cp", "Specific heat material property name");
+  return params;
+}
+
+PostprocessorVelocityFunctionTemperatureInflowBC::PostprocessorVelocityFunctionTemperatureInflowBC(
+    const InputParameters & parameters)
+  : DerivativeMaterialInterface<JvarMapIntegratedBCInterface<PostprocessorVelocityFunctionInflowBC>>(parameters),
+    _rho(getMaterialProperty<Real>("rho")),
+    _d_rho_d_u(getMaterialPropertyDerivative<Real>("rho", _var.name())),
+    _cp(getMaterialProperty<Real>("cp")),
+    _d_cp_d_u(getMaterialPropertyDerivative<Real>("cp", _var.name()))
+{
+}
+
+void
+PostprocessorVelocityFunctionTemperatureInflowBC::initialSetup()
+{
+  validateNonlinearCoupling<Real>("rho");
+  validateNonlinearCoupling<Real>("cp");
+}
+
+Real
+PostprocessorVelocityFunctionTemperatureInflowBC::computeQpResidual()
+{
+  return _rho[_qp] * _cp[_qp] * PostprocessorVelocityFunctionInflowBC::computeQpResidual();
+}
+
+Real
+PostprocessorVelocityFunctionTemperatureInflowBC::computeQpJacobian()
+{
+  return _rho[_qp] * _cp[_qp] * PostprocessorVelocityFunctionInflowBC::computeQpJacobian() +
+         _d_rho_d_u[_qp] * _phi[_j][_qp] * _cp[_qp] * PostprocessorVelocityFunctionInflowBC::computeQpResidual() +
+         _rho[_qp] * _d_cp_d_u[_qp] * _phi[_j][_qp] * PostprocessorVelocityFunctionInflowBC::computeQpResidual();
+}

src/bcs/VelocityFunctionTemperatureOutflowBC.C

@@ -10,6 +10,8 @@ VelocityFunctionTemperatureOutflowBC::validParams()
   params.addRequiredParam<FunctionName>("vel_x_func", "The x velocity function");
   params.addRequiredParam<FunctionName>("vel_y_func", "The y velocity function");
   params.addRequiredParam<FunctionName>("vel_z_func", "The z velocity function");
+  params.addParam<MaterialPropertyName>("rho", "rho", "Density material property name");
+  params.addParam<MaterialPropertyName>("cp", "cp", "Specific heat material property name");
   return params;
 }
 

src/dgkernels/DGFunctionTemperatureAdvection.C

@@ -7,6 +7,8 @@ DGFunctionTemperatureAdvection::validParams()
 {
   InputParameters params = DGFunctionConvection::validParams();
   params.addClassDescription("DG upwinding for func temp convection");
+  params.addParam<MaterialPropertyName>("rho", "rho", "Density material property name");
+  params.addParam<MaterialPropertyName>("cp", "cp", "Specific heat material property name");
   return params;
 }
 

src/interface_kernels/InterTemperatureAdvection.C

@@ -42,12 +42,12 @@ InterTemperatureAdvection::computeQpResidual(Moose::DGResidualType type)
       if (vdotn >= 0)
         r += vdotn * _u[_qp] * _test[_i][_qp];
       else
-        r += vdotn * _neighbor_value[_qp] * _test[_i][_qp] - _heat_source;
+        r += (vdotn * _neighbor_value[_qp] - _heat_source) * _test[_i][_qp];
       break;
 
     case Moose::Neighbor:
       if (vdotn >= 0)
-        r -= vdotn * _u[_qp] * _test_neighbor[_i][_qp] - _heat_source;
+        r -= (vdotn * _u[_qp] - _heat_source) * _test_neighbor[_i][_qp];
       else
         r -= vdotn * _neighbor_value[_qp] * _test_neighbor[_i][_qp];
       break;

src/kernels/VelocityFunctionTemperatureAdvection.C

@@ -10,6 +10,8 @@ VelocityFunctionTemperatureAdvection::validParams()
   params.addRequiredParam<FunctionName>("vel_x_func", "The x velocity function");
   params.addRequiredParam<FunctionName>("vel_y_func", "The y velocity function");
   params.addRequiredParam<FunctionName>("vel_z_func", "The z velocity function");
+  params.addParam<MaterialPropertyName>("rho", "rho", "Density material property name");
+  params.addParam<MaterialPropertyName>("cp", "cp", "Specific heat material property name");
   return params;
 }
 

tests/auxkernels/density_from_log/density_from_log.i

@@ -0,0 +1,58 @@
+# DensityFromLog exponentiates a log-density field. Here the log-density falls
+# linearly with x, so the recovered density is rho = exp(-x/H).
+
+H = 0.5
+
+[Mesh]
+  [column]
+    type = GeneratedMeshGenerator
+    dim = 1
+    nx = 4
+    xmax = 1
+  []
+[]
+
+[Problem]
+  solve = false
+  kernel_coverage_check = false
+[]
+
+[Variables]
+  [dummy]
+  []
+[]
+
+[AuxVariables]
+  [log_rho]
+  []
+  [rho]
+  []
+[]
+
+[Functions]
+  [log_rho_func]
+    type = ParsedFunction
+    expression = '-x / ${H}'
+  []
+[]
+
+[AuxKernels]
+  [log_rho]
+    type = FunctionAux
+    variable = log_rho
+    function = log_rho_func
+  []
+  [rho]
+    type = DensityFromLog
+    variable = rho
+    density_log = log_rho
+  []
+[]
+
+[Executioner]
+  type = Steady
+[]
+
+[Outputs]
+  exodus = true
+[]

tests/auxkernels/density_from_log/tests

@@ -0,0 +1,16 @@
+[Tests]
+  [density_from_log]
+    type = Exodiff
+    input = 'density_from_log.i'
+    exodiff = 'density_from_log_out.e'
+    requirement = 'The system shall recover a density field by exponentiating a log-density variable.'
+  []
+  [larger_scale_height]
+    type = Exodiff
+    input = 'density_from_log.i'
+    cli_args = 'H=2
+                Outputs/file_base=larger_scale_height_out'
+    exodiff = 'larger_scale_height_out.e'
+    requirement = 'The system shall flatten the recovered density profile as the log-density scale height increases.'
+  []
+[]

tests/auxkernels/function_derivative_aux/function_derivative_aux.i

@@ -0,0 +1,61 @@
+# FunctionDerivativeAux samples a component of a function gradient. For the
+# field u(x) = x^2 the x derivative is 2x, evaluated both at nodes and at
+# element quadrature points.
+
+[Mesh]
+  [rod]
+    type = GeneratedMeshGenerator
+    dim = 1
+    nx = 4
+    xmax = 1
+  []
+[]
+
+[Problem]
+  solve = false
+  kernel_coverage_check = false
+[]
+
+[Variables]
+  [dummy]
+  []
+[]
+
+[Functions]
+  [quadratic]
+    type = ParsedFunction
+    expression = 'x * x'
+  []
+[]
+
+[AuxVariables]
+  [dudx_nodal]
+  []
+  [dudx_elemental]
+    order = CONSTANT
+    family = MONOMIAL
+  []
+[]
+
+[AuxKernels]
+  [nodal]
+    type = FunctionDerivativeAux
+    variable = dudx_nodal
+    function = quadratic
+    component = 1
+  []
+  [elemental]
+    type = FunctionDerivativeAux
+    variable = dudx_elemental
+    function = quadratic
+    component = 1
+  []
+[]
+
+[Executioner]
+  type = Steady
+[]
+
+[Outputs]
+  exodus = true
+[]

tests/auxkernels/function_derivative_aux/tests

@@ -0,0 +1,17 @@
+[Tests]
+  [x_derivative]
+    type = Exodiff
+    input = 'function_derivative_aux.i'
+    exodiff = 'function_derivative_aux_out.e'
+    requirement = 'The system shall evaluate a component of a function gradient at nodes and at element quadrature points.'
+  []
+  [y_derivative]
+    type = Exodiff
+    input = 'function_derivative_aux.i'
+    cli_args = 'AuxKernels/nodal/component=2
+                AuxKernels/elemental/component=2
+                Outputs/file_base=y_derivative_out'
+    exodiff = 'y_derivative_out.e'
+    requirement = 'The system shall return a zero derivative for the direction in which the sampled function is constant.'
+  []
+[]