ADBoundaryFlux3EqnGhostBase.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 "ADBoundaryFlux3EqnGhostBase.h"
#include "ADNumericalFlux3EqnBase.h"

#include "libmesh/elem.h"

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

  params.addClassDescription("Computes boundary fluxes for the 1-D, variable-area Euler equations "
                             "using a numerical flux user object and a ghost cell solution");

  params.addRequiredParam<UserObjectName>("numerical_flux", "Name of numerical flux user object");
  params.addRequiredParam<Real>("normal", "Outward normal");

  return params;
}

ADBoundaryFlux3EqnGhostBase::ADBoundaryFlux3EqnGhostBase(const InputParameters & parameters)
  : ADBoundaryFluxBase(parameters),
    _numerical_flux(getUserObject<ADNumericalFlux3EqnBase>("numerical_flux")),
    _normal(getParam<Real>("normal"))
{
}

void
ADBoundaryFlux3EqnGhostBase::calcFlux(unsigned int iside,
                                      dof_id_type ielem,
                                      const std::vector<ADReal> & U1,
                                      const RealVectorValue & /*normal*/,
                                      std::vector<ADReal> & flux) const
{
  const Elem * elem = _subproblem.mesh().elemPtr(ielem);
  const Elem * side_elem = elem->build_side_ptr(iside).release();
  const Point side_center_point = side_elem->vertex_average();
  delete side_elem;

  const std::vector<ADReal> U2 = getGhostCellSolution(U1, side_center_point);
  flux = _numerical_flux.getFlux(iside, ielem, true, U1, U2, _normal);
}