TotalEnergyConvectiveFlux.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "TotalEnergyConvectiveFlux.h"

registerMooseObject("NavierStokesApp", TotalEnergyConvectiveFlux);

InputParameters
TotalEnergyConvectiveFlux::validParams()
{
  InputParameters params = Kernel::validParams();

  params.addRequiredCoupledVar("rho", "density");
  params.addRequiredCoupledVar("rho_u", "rho*u");
  params.addCoupledVar("rho_v", "rho*v");
  params.addCoupledVar("rho_w", "rho*w");
  params.addRequiredCoupledVar("enthalpy", "Enthalpy");
  params.addClassDescription(
      "Implements the advection term for the Navier Stokes energy equation.");

  return params;
}

TotalEnergyConvectiveFlux::TotalEnergyConvectiveFlux(const InputParameters & parameters)
  : Kernel(parameters),
    _rho(coupledValue("rho")),
    _rho_u(coupledValue("rho_u")),
    _rho_v(isCoupled("rho_v") ? coupledValue("rho_v") : _zero),
    _rho_w(isCoupled("rho_w") ? coupledValue("rho_w") : _zero),
    _enthalpy(coupledValue("enthalpy"))
{
}

Real
TotalEnergyConvectiveFlux::computeQpResidual()
{
  RealVectorValue rho_u_vec(_rho_u[_qp], _rho_v[_qp], _rho_w[_qp]);
  return -_enthalpy[_qp] * (rho_u_vec * _grad_test[_i][_qp]);
}

Real
TotalEnergyConvectiveFlux::computeQpJacobian()
{
  return 0.;
}