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 "INSFVMomentumGravity.h"
      11             : #include "NS.h"
      12             : 
      13             : registerMooseObject("NavierStokesApp", INSFVMomentumGravity);
      14             : 
      15             : InputParameters
      16        2601 : INSFVMomentumGravity::validParams()
      17             : {
      18        2601 :   InputParameters params = FVElementalKernel::validParams();
      19        2601 :   params += INSFVMomentumResidualObject::validParams();
      20        2601 :   params.addClassDescription(
      21             :       "Computes a body force due to gravity in Rhie-Chow based simulations.");
      22        5202 :   params.addRequiredParam<RealVectorValue>("gravity", "Direction of the gravity vector");
      23        2601 :   params.addParam<MooseFunctorName>(NS::density, NS::density, "The value for the density");
      24        2601 :   return params;
      25           0 : }
      26             : 
      27         584 : INSFVMomentumGravity::INSFVMomentumGravity(const InputParameters & params)
      28             :   : FVElementalKernel(params),
      29             :     INSFVMomentumResidualObject(*this),
      30         584 :     _gravity(getParam<RealVectorValue>("gravity")),
      31        1168 :     _rho(getFunctor<ADReal>(NS::density))
      32             : {
      33         584 : }
      34             : 
      35             : ADReal
      36     4197116 : INSFVMomentumGravity::computeQpResidual()
      37             : {
      38     8394232 :   return -_rho(makeElemArg(_current_elem), determineState()) * _gravity(_index);
      39             : }