Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "PINSFVMomentumBoussinesq.h"
      11             : #include "PINSFVSuperficialVelocityVariable.h"
      12             : #include "NS.h"
      13             : 
      14             : registerMooseObject("NavierStokesApp", PINSFVMomentumBoussinesq);
      15             : 
      16             : InputParameters
      17         213 : PINSFVMomentumBoussinesq::validParams()
      18             : {
      19         213 :   InputParameters params = INSFVMomentumBoussinesq::validParams();
      20         213 :   params.addClassDescription(
      21             :       "Computes a body force for natural convection buoyancy in porous media: eps alpha (T-T_0)");
      22         213 :   params.addRequiredParam<MooseFunctorName>(NS::porosity, "Porosity auxiliary variable");
      23             : 
      24         213 :   return params;
      25           0 : }
      26             : 
      27         124 : PINSFVMomentumBoussinesq::PINSFVMomentumBoussinesq(const InputParameters & params)
      28         124 :   : INSFVMomentumBoussinesq(params), _eps(getFunctor<ADReal>(NS::porosity))
      29             : {
      30         124 :   if (!dynamic_cast<PINSFVSuperficialVelocityVariable *>(&_var))
      31           0 :     mooseError("PINSFVMomentumBoussinesq may only be used with a superficial velocity "
      32             :                "variable, of variable type PINSFVSuperficialVelocityVariable.");
      33         124 : }
      34             : 
      35             : ADReal
      36      676672 : PINSFVMomentumBoussinesq::computeQpResidual()
      37             : {
      38      676672 :   return _eps(makeElemArg(_current_elem), determineState()) *
      39     1353344 :          INSFVMomentumBoussinesq::computeQpResidual();
      40             : }