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 "PINSFEFluidVelocityTimeDerivative.h"
      11             : 
      12             : registerMooseObject("NavierStokesApp", PINSFEFluidVelocityTimeDerivative);
      13             : registerMooseObjectRenamed("NavierStokesApp",
      14             :                            PMFluidVelocityTimeDerivative,
      15             :                            "02/01/2024 00:00",
      16             :                            PINSFEFluidVelocityTimeDerivative);
      17             : 
      18             : InputParameters
      19         228 : PINSFEFluidVelocityTimeDerivative::validParams()
      20             : {
      21         228 :   InputParameters params = TimeDerivative::validParams();
      22         456 :   params.addParam<bool>("conservative_form", false, "if conservative form is used");
      23         456 :   params.addRequiredCoupledVar("pressure", "coupled pressure");
      24         456 :   params.addRequiredCoupledVar("temperature", "coupled temperature");
      25         456 :   params.addRequiredParam<UserObjectName>("eos", "The name of equation of state object to use.");
      26         228 :   params.addClassDescription("Add the transient term for one component of the porous media "
      27             :                              "momentum conservation equation");
      28             : 
      29         228 :   return params;
      30           0 : }
      31             : 
      32         120 : PINSFEFluidVelocityTimeDerivative::PINSFEFluidVelocityTimeDerivative(
      33         120 :     const InputParameters & parameters)
      34             :   : TimeDerivative(parameters),
      35         120 :     _conservative_form(getParam<bool>("conservative_form")),
      36         120 :     _pressure(coupledValue("pressure")),
      37         120 :     _temperature(coupledValue("temperature")),
      38         120 :     _temperature_dot(coupledDot("temperature")),
      39         120 :     _pressure_dot(coupledDot("pressure")),
      40         240 :     _rho(getMaterialProperty<Real>("rho_fluid")),
      41         240 :     _eos(getUserObject<SinglePhaseFluidProperties>("eos"))
      42             : {
      43         120 : }
      44             : 
      45             : Real
      46    78095360 : PINSFEFluidVelocityTimeDerivative::computeQpResidual()
      47             : {
      48    78095360 :   Real res = _rho[_qp] * TimeDerivative::computeQpResidual();
      49    78095360 :   if (_conservative_form)
      50             :   {
      51             :     Real rho, drho_dp, drho_dT;
      52    29862400 :     _eos.rho_from_p_T(_pressure[_qp], _temperature[_qp], rho, drho_dp, drho_dT);
      53    29862400 :     Real drho_dt = drho_dT * _temperature_dot[_qp] + drho_dp * _pressure_dot[_qp];
      54    29862400 :     res += _u[_qp] * drho_dt * _test[_i][_qp];
      55             :   }
      56             : 
      57    78095360 :   return res;
      58             : }
      59             : 
      60             : Real
      61     8657920 : PINSFEFluidVelocityTimeDerivative::computeQpJacobian()
      62             : {
      63     8657920 :   Real jac = _rho[_qp] * TimeDerivative::computeQpJacobian();
      64     8657920 :   if (_conservative_form)
      65             :   {
      66             :     Real rho, drho_dp, drho_dT;
      67      665600 :     _eos.rho_from_p_T(_pressure[_qp], _temperature[_qp], rho, drho_dp, drho_dT);
      68      665600 :     Real drho_dt = drho_dT * _temperature_dot[_qp] + drho_dp * _pressure_dot[_qp];
      69      665600 :     jac += _phi[_j][_qp] * drho_dt * _test[_i][_qp];
      70             :   }
      71             : 
      72     8657920 :   return jac;
      73             : }