MassDiffusionSpeciesGasMixDGKernel.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 "MassDiffusionSpeciesGasMixDGKernel.h"

registerMooseObject("ThermalHydraulicsApp", MassDiffusionSpeciesGasMixDGKernel);

InputParameters
MassDiffusionSpeciesGasMixDGKernel::validParams()
{
  InputParameters params = MassDiffusionBaseGasMixDGKernel::validParams();

  params.addClassDescription("Adds mass diffusion to the species equation for FlowChannelGasMix.");

  return params;
}

MassDiffusionSpeciesGasMixDGKernel::MassDiffusionSpeciesGasMixDGKernel(
    const InputParameters & parameters)
  : MassDiffusionBaseGasMixDGKernel(parameters)
{
}

ADReal
MassDiffusionSpeciesGasMixDGKernel::computeQpFlux() const
{
  Real dx, dx_side;
  computePositionChanges(dx, dx_side);

  const ADReal rho = linearlyInterpolate(_rho_elem[_qp], _rho_neig[_qp], dx, dx_side);
  const ADReal D = linearlyInterpolate(_D_elem[_qp], _D_neig[_qp], dx, dx_side);
  const ADReal dxi_dx = computeGradient(_mass_fraction_elem[_qp], _mass_fraction_neig[_qp], dx);

  return -rho * D * dxi_dx;
}