437619c6e7
git-svn-id: http://svn.sintef.no/trondheim/IFEM/trunk@835 e10b68d5-8a6e-419e-a041-bce267b0401d
124 lines
5.7 KiB
C++
124 lines
5.7 KiB
C++
//==============================================================================
|
|
//!
|
|
//! \file NavierStokesG2MP.h
|
|
//!
|
|
//! \date August 11 2010
|
|
//!
|
|
//! \author Runar Holdahl / SINTEF
|
|
//!
|
|
//! \brief Integrand implementations for Navier-Stokes equations with G2 stabilization
|
|
//! and. Time integration by mid-point method
|
|
//!
|
|
//==============================================================================
|
|
|
|
#ifndef _NAVIER_STOKES_G2_MP_H
|
|
#define _NAVIER_STOKES_G2_MP_H
|
|
|
|
#include "NavierStokesG2.h"
|
|
|
|
/*!
|
|
\brief Class representing the integrand of the pressure stabilized Navier-Stokes problem.
|
|
\details This class supports G2-stabilized Navier-Stokes formulations using equal
|
|
order elements for velocity and pressure. Time integration by the mid-point method.
|
|
*/
|
|
|
|
|
|
class NavierStokesG2MP : public NavierStokesG2
|
|
{
|
|
public:
|
|
|
|
//! \brief The default constructor initializes all pointers to zero
|
|
//! \param[in] n Number of spatial dimensions
|
|
NavierStokesG2MP(short int n, SIM::Formulation form, int itg = 3);
|
|
//! \brief The destructor clears the static work arrays used internally.
|
|
virtual ~NavierStokesG2MP();
|
|
|
|
//! \brief Evaluates the integrand at an interior point.
|
|
//! \param elmInt The local integral object to receive the contributions
|
|
//! \param[in] time Parameters for nonlinear and time-dependent simulations
|
|
//! \param[in] detJW Jacobian determinant times integration point weight
|
|
//! \param[in] N Basis function values
|
|
//! \param[in] dNdX Basis function gradients
|
|
//! \param[in] X Cartesian coordinates of current integration point
|
|
//!
|
|
//! \details This interface is used when \a getIntegrandType returns 1.
|
|
virtual bool evalInt(LocalIntegral*& elmInt,
|
|
const TimeDomain& time, double detJW,
|
|
const Vector& N, const Matrix& dNdX,
|
|
const Vec3& X) const
|
|
{ return false; }
|
|
|
|
//! \brief Evaluates the integrand at an interior point.
|
|
//! \param elmInt The local integral object to receive the contributions
|
|
//! \param[in] detJW Jacobian determinant times integration point weight
|
|
//! \param[in] time Parameters for nonlinear and time-dependent simulations
|
|
//! \param[in] N Basis function values
|
|
//! \param[in] dNdX Basis function gradients
|
|
//! \param[in] d2NdX2 Basis function second derivatives
|
|
//! \param[in] X Cartesian coordinates of current integration point
|
|
//! \param[in] h Characteristic element length
|
|
//!
|
|
//! \details This interface is used when \a getIntegrandType returns 2.
|
|
//! The default implementation forwards to the stationary version.
|
|
//! Reimplement this method for time-dependent or non-linear problems.
|
|
virtual bool evalInt(LocalIntegral*& elmInt,
|
|
const TimeDomain& time, double detJW,
|
|
const Vector& N, const Matrix& dNdX,
|
|
const Matrix3D& d2NdX2, const Vec3& X,
|
|
double h = 0.0) const;
|
|
|
|
//! \brief Evaluates the integrand at an interior point.
|
|
//! \param elmInt The local integral object to receive the contributions
|
|
//! \param[in] detJW Jacobian determinant times integration point weight
|
|
//! \param[in] time Parameters for nonlinear and time-dependent simulations
|
|
//! \param[in] N Basis function values
|
|
//! \param[in] dNdX Basis function gradients
|
|
//! \param[in] Navg Volume-averaged basis function values over the element
|
|
//! \param[in] X Cartesian coordinates of current integration point
|
|
//!
|
|
//! \details This interface is used when \a getIntegrandType returns 3.
|
|
//! The default implementation forwards to the stationary version.
|
|
//! Reimplement this method for time-dependent or non-linear problems.
|
|
virtual bool evalInt(LocalIntegral*& elmInt,
|
|
const TimeDomain& time, double detJW,
|
|
const Vector& N, const Matrix& dNdX,
|
|
const Vector& Navg, const Vec3& X) const
|
|
{ return true; }
|
|
|
|
//! \brief Evaluates the integrand at an interior point.
|
|
//! \param elmInt The local integral object to receive the contributions
|
|
//! \param[in] detJW Jacobian determinant times integration point weight
|
|
//! \param[in] N Basis function values
|
|
//! \param[in] dNdX Basis function gradients
|
|
//! \param[in] X Cartesian coordinates of current integration point
|
|
virtual bool evalInt(LocalIntegral*& elmInt, double detJW,
|
|
const Vector& N, const Matrix& dNdX,
|
|
const Vec3& X) const
|
|
{ return false; }
|
|
//! \brief Evaluates the integrand at an interior point.
|
|
//! \param elmInt The local integral object to receive the contributions
|
|
//! \param[in] detJW Jacobian determinant times integration point weight
|
|
//! \param[in] N Basis function values
|
|
//! \param[in] dNdX Basis function gradients
|
|
//! \param[in] d2NdX2 Basis function second derivatives
|
|
//! \param[in] X Cartesian coordinates of current integration point
|
|
//! \param[in] h characteristic element length
|
|
virtual bool evalInt(LocalIntegral*& elmInt, double detJW,
|
|
const Vector& N, const Matrix& dNdX,
|
|
const Matrix3D& d2NdX2,
|
|
const Vec3& X, double h = 0.0) const;
|
|
//! \brief Evaluates the integrand at an interior point.
|
|
//! \param elmInt The local integral object to receive the contributions
|
|
//! \param[in] detJW Jacobian determinant times integration point weight
|
|
//! \param[in] N Basis function values
|
|
//! \param[in] dNdX Basis function gradients
|
|
//! \param[in] Navg Average value of basis function over element
|
|
//! \param[in] X Cartesian coordinates of current integration point
|
|
virtual bool evalInt(LocalIntegral*& elmInt, double detJW,
|
|
const Vector& N, const Matrix& dNdX,
|
|
const Vector& Navg, const Vec3& X) const
|
|
{ return false; }
|
|
};
|
|
|
|
#endif
|