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 "INSADGravityForce.h"
      11             : #include "INSADObjectTracker.h"
      12             : #include "FEProblemBase.h"
      13             : 
      14             : registerMooseObject("NavierStokesApp", INSADGravityForce);
      15             : 
      16             : InputParameters
      17         304 : INSADGravityForce::validParams()
      18             : {
      19         304 :   InputParameters params = ADVectorKernelValue::validParams();
      20         304 :   params.addClassDescription("Computes a body force due to gravity.");
      21         608 :   params.addRequiredParam<RealVectorValue>("gravity", "Direction of the gravity vector");
      22         304 :   return params;
      23           0 : }
      24             : 
      25         176 : INSADGravityForce::INSADGravityForce(const InputParameters & parameters)
      26             :   : ADVectorKernelValue(parameters),
      27         352 :     _gravity_strong_residual(getADMaterialProperty<RealVectorValue>("gravity_strong_residual"))
      28             : {
      29             :   // Bypass the UserObjectInterface method because it requires a UserObjectName param which we
      30             :   // don't need
      31             :   auto & obj_tracker = const_cast<INSADObjectTracker &>(
      32         176 :       _fe_problem.getUserObject<INSADObjectTracker>("ins_ad_object_tracker"));
      33         352 :   for (const auto block_id : blockIDs())
      34             :   {
      35         176 :     obj_tracker.set("has_gravity", true, block_id);
      36         528 :     obj_tracker.set("gravity", getParam<RealVectorValue>("gravity"), block_id);
      37             :   }
      38         176 : }
      39             : 
      40             : ADRealVectorValue
      41    15527136 : INSADGravityForce::precomputeQpResidual()
      42             : {
      43    15527136 :   return _gravity_strong_residual[_qp];
      44             : }