ADBoundaryFlux3EqnFreeOutflow.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 "ADBoundaryFlux3EqnFreeOutflow.h"
#include "THMIndicesVACE.h"
#include "Numerics.h"

registerMooseObject("ThermalHydraulicsApp", ADBoundaryFlux3EqnFreeOutflow);

InputParameters
ADBoundaryFlux3EqnFreeOutflow::validParams()
{
  InputParameters params = ADBoundaryFluxBase::validParams();

  params.addClassDescription("Computes the outflow boundary flux directly for the 1-D, 1-phase, "
                             "variable-area Euler equations");

  params.addRequiredParam<UserObjectName>("fluid_properties",
                                          "Name of fluid properties user object");

  return params;
}

ADBoundaryFlux3EqnFreeOutflow::ADBoundaryFlux3EqnFreeOutflow(const InputParameters & parameters)
  : ADBoundaryFluxBase(parameters),

    _fp(getUserObject<SinglePhaseFluidProperties>("fluid_properties"))
{
}

void
ADBoundaryFlux3EqnFreeOutflow::calcFlux(unsigned int /*iside*/,
                                        dof_id_type /*ielem*/,
                                        const std::vector<ADReal> & U1,
                                        const RealVectorValue & /*normal*/,
                                        std::vector<ADReal> & flux) const
{
  // extract the solution and area
  const ADReal rhoA1 = U1[THMVACE1D::RHOA];
  const ADReal rhouA1 = U1[THMVACE1D::RHOUA];
  const ADReal rhoEA1 = U1[THMVACE1D::RHOEA];
  const ADReal A1 = U1[THMVACE1D::AREA];

  const ADReal rho1 = rhoA1 / A1;
  const ADReal vel1 = rhouA1 / rhoA1;
  const ADReal v1 = 1.0 / rho1;
  const ADReal e1 = rhoEA1 / rhoA1 - 0.5 * vel1 * vel1;
  const ADReal p1 = _fp.p_from_v_e(v1, e1);

  flux.resize(THMVACE1D::N_FLUX_OUTPUTS);
  flux[THMVACE1D::MASS] = vel1 * rhoA1;
  flux[THMVACE1D::MOMENTUM] = (vel1 * rhouA1 + p1 * A1);
  flux[THMVACE1D::ENERGY] = vel1 * (rhoEA1 + p1 * A1);
}