NSMomentumInviscidBC.h

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
 
#pragma once
 
#include "NSIntegratedBC.h"
#include "NSPressureDerivs.h"
 
// Forward Declarations
 
class NSMomentumInviscidBC : public NSIntegratedBC
{
public:
  static InputParameters validParams();
 
  NSMomentumInviscidBC(const InputParameters & parameters);
 
protected:
  // Which spatial component of the momentum equations (0,1, or 2) is this
  // kernel applied in?
  const unsigned _component;
 
  // An object for computing pressure derivatives.
  // Constructed via a reference to ourself
  NSPressureDerivs<NSMomentumInviscidBC> _pressure_derivs;
 
  // Declare ourselves friend to the helper class.
  template <class U>
  friend class NSPressureDerivs;
 
  // These functions can be mix-n-matched by derived classes to implement
  // any of the following boundary conditions:
  // .) Fully unspecified (both (rho*u)(u.n) and p computed implicitly, is this valid?)
  // .) Specified pressure/unspecified (rho*u)(u.n)
  // .) Unspecified pressure/specified (rho*u)(u.n)
  // .) Fully specified (both pressure and (rho*u)(u.n) given, this may not be physically
  // meaningful?)
 
  // Depending on the passed-in value, will compute the residual for either a specified
  // pressure value or the residual at the current value of the pressure.
  Real pressureQpResidualHelper(Real pressure);
 
  // If the pressure is fixed, the Jacobian of the pressure term is zero, otherwise
  // we return the Jacobian value for the passed-in variable number.
  Real pressureQpJacobianHelper(unsigned var_number);
 
  // Depending on the passed-in vector, will compute the residual for either a specified
  // value of (rho*u)(u.n) or the residual at the current value of (rho*u)(u.n).
  // The passed-in value is the _component'th entry of the (rho*u)(u.n) vector.
  Real convectiveQpResidualHelper(Real rhou_udotn);
 
  // If the value of (rho*u)(u.n) is fixed, the Jacobian of the
  // convective term is zero, otherwise we return the correct value
  // based on the passed-in variable number.
  Real convectiveQpJacobianHelper(unsigned var_number);
};