NSMomentumInviscidFluxWithGradP.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 "NSKernel.h"
#include "NSPressureDerivs.h"

// ForwardDeclarations

class NSMomentumInviscidFluxWithGradP : public NSKernel
{
public:
  static InputParameters validParams();

  NSMomentumInviscidFluxWithGradP(const InputParameters & parameters);

protected:
  virtual Real computeQpResidual();
  virtual Real computeQpJacobian();
  virtual Real computeQpOffDiagJacobian(unsigned int jvar);

  // Coupled gradients
  const VariableGradient & _grad_p;

  // Parameters
  const unsigned int _component;

private:
  // Computes the Jacobian contribution due to the pressure term,
  // by summing over the appropriate Hessian row.
  Real pressureQpJacobianHelper(unsigned var_number);

  // Single vector to refer to all gradients.  We have to store
  // pointers since you can't have a vector<Foo&>.  Initialized in
  // the ctor.
  std::vector<const VariableGradient *> _gradU;

  // An object for computing pressure derivatives.
  // Constructed via a reference to ourself
  NSPressureDerivs<NSMomentumInviscidFluxWithGradP> _pressure_derivs;

  // Declare ourselves friend to the helper class.
  template <class U>
  friend class NSPressureDerivs;
};