RichardsSUPG.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
//  Base class for Richards SUPG
//

#include "GeneralUserObject.h"

class RichardsSUPG : public GeneralUserObject
{
public:
  static InputParameters validParams();

  RichardsSUPG(const InputParameters & parameters);

  void initialize();
  void execute();
  void finalize();

  virtual RealVectorValue velSUPG(RealTensorValue perm,
                                  RealVectorValue gradp,
                                  Real density,
                                  RealVectorValue gravity) const = 0;

  virtual RealTensorValue dvelSUPG_dgradp(RealTensorValue perm) const = 0;

  virtual RealVectorValue
  dvelSUPG_dp(RealTensorValue perm, Real density_prime, RealVectorValue gravity) const = 0;

  virtual RealVectorValue bb(RealVectorValue vel,
                             int dimen,
                             RealVectorValue xi_prime,
                             RealVectorValue eta_prime,
                             RealVectorValue zeta_prime) const = 0;

  virtual RealVectorValue dbb2_dgradp(RealVectorValue vel,
                                      RealTensorValue dvel_dgradp,
                                      RealVectorValue xi_prime,
                                      RealVectorValue eta_prime,
                                      RealVectorValue zeta_prime) const = 0;

  virtual Real dbb2_dp(RealVectorValue vel,
                       RealVectorValue dvel_dp,
                       RealVectorValue xi_prime,
                       RealVectorValue eta_prime,
                       RealVectorValue zeta_prime) const = 0;

  virtual Real tauSUPG(RealVectorValue vel, Real traceperm, RealVectorValue b) const = 0;

  virtual RealVectorValue dtauSUPG_dgradp(RealVectorValue vel,
                                          RealTensorValue dvel_dgradp,
                                          Real traceperm,
                                          RealVectorValue b,
                                          RealVectorValue db2_dgradp) const = 0;

  virtual Real dtauSUPG_dp(RealVectorValue vel,
                           RealVectorValue dvel_dp,
                           Real traceperm,
                           RealVectorValue b,
                           Real db2_dp) const = 0;

  virtual bool SUPG_trivial() const = 0;
};