Implementation of the DDR scheme for the Stokes problem in curl-curl formulation. More...

Collaboration diagram for DDR_stokes:

Classes
struct	HArDCore3D::StokesNorms
	Structure to store norm components (for velocity, its curl, pressure, and its gradient), as well as Hcurl and Hgrad norms of velocity and pressure. More...

struct	HArDCore3D::Stokes
	Assemble a Stokes problem. More...

Typedefs
typedef Eigen::SparseMatrix< double >	HArDCore3D::Stokes::SystemMatrixType

typedef std::function< Eigen::Vector3d(const Eigen::Vector3d &)>	HArDCore3D::Stokes::ForcingTermType

typedef std::function< Eigen::Vector3d(const Eigen::Vector3d &)>	HArDCore3D::Stokes::VelocityType

typedef std::function< Eigen::Vector3d(const Eigen::Vector3d &)>	HArDCore3D::Stokes::VorticityType

typedef std::function< double(const Eigen::Vector3d &)>	HArDCore3D::Stokes::PressureType

typedef std::function< Eigen::Vector3d(const Eigen::Vector3d &)>	HArDCore3D::Stokes::PressureGradientType

typedef IntegralWeight	HArDCore3D::Stokes::ViscosityType

Functions
	HArDCore3D::StokesNorms::StokesNorms (double norm_u, double norm_curl_u, double norm_p, double norm_grad_p)
	Constructor.

	HArDCore3D::StokesNorms::StokesNorms ()

	HArDCore3D::Stokes::Stokes (const DDRCore &ddrcore, bool use_threads, std::ostream &output=std::cout)
	Constructor.

void	HArDCore3D::Stokes::assembleLinearSystem (const ForcingTermType &f, const VelocityType &u, const VorticityType &omega, const ViscosityType &nu)
	Assemble the global system

size_t	HArDCore3D::Stokes::dimension () const
	Returns the global problem dimension (without Lagrange multiplier, just velocity & pressure)

size_t	HArDCore3D::Stokes::nbSCDOFs_u () const
	Returns the number of statically condensed DOFs (velocity)

size_t	HArDCore3D::Stokes::nbSCDOFs_p () const
	Returns the number of statically condensed DOFs (pressure)

size_t	HArDCore3D::Stokes::nbSCDOFs () const
	Returns the number of statically condensed DOFs (both velocity and pressure)

size_t	HArDCore3D::Stokes::sizeSystem () const
	Returns the size of the statically condensed system (with Lagrange multiplier)

const XGrad &	HArDCore3D::Stokes::xGrad () const
	Returns the space XGrad.

const XCurl &	HArDCore3D::Stokes::xCurl () const
	Returns the space XCurl.

const XDiv &	HArDCore3D::Stokes::xDiv () const
	Returns the space XDiv.

const SystemMatrixType &	HArDCore3D::Stokes::systemMatrix () const
	Returns the linear system matrix.

SystemMatrixType &	HArDCore3D::Stokes::systemMatrix ()
	Returns the linear system matrix.

const Eigen::VectorXd &	HArDCore3D::Stokes::systemVector () const
	Returns the linear system right-hand side vector.

Eigen::VectorXd &	HArDCore3D::Stokes::systemVector ()
	Returns the linear system right-hand side vector.

const SystemMatrixType &	HArDCore3D::Stokes::scMatrix () const
	Returns the static condensation recovery operator.

Eigen::VectorXd &	HArDCore3D::Stokes::scVector ()
	Returns the static condensation rhs.

const double &	HArDCore3D::Stokes::stabilizationParameter () const
	Returns the stabilization parameter.

double &	HArDCore3D::Stokes::stabilizationParameter ()
	Returns the stabilization parameter.

std::vector< StokesNorms >	HArDCore3D::Stokes::computeStokesNorms (const std::vector< Eigen::VectorXd > &list_dofs) const
	Compute the discrete L2 norms, for a family of Eigen::VectorXd representing velocities & pressures, of u, curl u, p and grad p.

std::pair< StokesNorms, StokesNorms >	HArDCore3D::Stokes::computeContinuousErrorsNorms (const Eigen::VectorXd &v, const VelocityType &u, const VorticityType &curl_u, const PressureType &p, const PressureGradientType &grad_p) const
	Compute the continuous Hcurl errors in u and Hgrad error in p (1st component), and same with continuous norms (second component)

double	HArDCore3D::Stokes::evaluateDefectCommutationGradInterp (const PressureType &p, const PressureGradientType &grad_p, const size_t deg_quad_interpolate) const
	Evaluate the defect of commutation G_h (Igrad p) = Icurl (grad p) by computing the norm of the difference (which should be 0).

Variables
double	HArDCore3D::StokesNorms::u
	Norm of velocity.

double	HArDCore3D::StokesNorms::curl_u
	Norm of curl of velocity (vorticity)

double	HArDCore3D::StokesNorms::p
	Norm of pressure.

double	HArDCore3D::StokesNorms::grad_p
	Norm of grad of pressure.

double	HArDCore3D::StokesNorms::hcurl_u
	Hcurl norm of u.

double	HArDCore3D::StokesNorms::hgrad_p
	Hgrad norm of p.

static const double	HArDCore3D::PI = boost::math::constants::pi<double>()

static Stokes::VelocityType	HArDCore3D::trigonometric_u

static Stokes::VorticityType	HArDCore3D::trigonometric_curl_u

static Stokes::PressureType	HArDCore3D::trigonometric_p

static Stokes::PressureGradientType	HArDCore3D::trigonometric_grad_p

static Stokes::ForcingTermType	HArDCore3D::trigonometric_f

static Stokes::ViscosityType	HArDCore3D::trigonometric_nu = Stokes::ViscosityType(1.)

static Stokes::VelocityType	HArDCore3D::linear_u

static Stokes::VorticityType	HArDCore3D::linear_curl_u

static Stokes::PressureType	HArDCore3D::linear_p

static Stokes::PressureGradientType	HArDCore3D::linear_grad_p

static Stokes::ForcingTermType	HArDCore3D::linear_f

static Stokes::ViscosityType	HArDCore3D::linear_nu = Stokes::ViscosityType(1.)

static Stokes::VelocityType	HArDCore3D::field_u

static Stokes::VorticityType	HArDCore3D::field_curl_u

static Stokes::PressureType	HArDCore3D::field_p

static Stokes::PressureGradientType	HArDCore3D::field_grad_p

static Stokes::ForcingTermType	HArDCore3D::field_f

static Stokes::ViscosityType	HArDCore3D::field_nu = Stokes::ViscosityType(1.)

Detailed Description

Implementation of the DDR scheme for the Stokes problem in curl-curl formulation.

Typedef Documentation

◆ ForcingTermType

typedef std::function<Eigen::Vector3d(const Eigen::Vector3d &)> HArDCore3D::Stokes::ForcingTermType

◆ PressureGradientType

typedef std::function<Eigen::Vector3d(const Eigen::Vector3d &)> HArDCore3D::Stokes::PressureGradientType

◆ PressureType

typedef std::function<double(const Eigen::Vector3d &)> HArDCore3D::Stokes::PressureType

◆ SystemMatrixType

typedef Eigen::SparseMatrix<double> HArDCore3D::Stokes::SystemMatrixType

◆ VelocityType

typedef std::function<Eigen::Vector3d(const Eigen::Vector3d &)> HArDCore3D::Stokes::VelocityType

◆ ViscosityType

typedef IntegralWeight HArDCore3D::Stokes::ViscosityType

◆ VorticityType

typedef std::function<Eigen::Vector3d(const Eigen::Vector3d &)> HArDCore3D::Stokes::VorticityType

Function Documentation

◆ assembleLinearSystem()

void Stokes::assembleLinearSystem	(	const ForcingTermType &	f,
		const VelocityType &	u,
		const VorticityType &	omega,
		const ViscosityType &	nu
	)

Assemble the global system

Parameters

f	Forcing term
u	Boundary condition on velocity
omega	Boundary condition on vorticity
nu	Viscosity

◆ computeContinuousErrorsNorms()

std::pair< StokesNorms, StokesNorms > Stokes::computeContinuousErrorsNorms	(	const Eigen::VectorXd &	v,
		const VelocityType &	u,
		const VorticityType &	curl_u,
		const PressureType &	p,
		const PressureGradientType &	grad_p
	)		const

Compute the continuous Hcurl errors in u and Hgrad error in p (1st component), and same with continuous norms (second component)

Parameters

v	The vector, contains both velocity and pressure
u	Continuous velocity
curl_u	Curl of continuous velocity
p	Continuous pressure
grad_p	Grad of continuous pressure

◆ computeStokesNorms()

std::vector< StokesNorms > Stokes::computeStokesNorms ( const std::vector< Eigen::VectorXd > & list_dofs ) const

Compute the discrete L2 norms, for a family of Eigen::VectorXd representing velocities & pressures, of u, curl u, p and grad p.

Parameters

list_dofs List of vectors representing velocities and pressures

◆ dimension()

size_t HArDCore3D::Stokes::dimension ( ) const

inline

Returns the global problem dimension (without Lagrange multiplier, just velocity & pressure)

◆ evaluateDefectCommutationGradInterp()

double Stokes::evaluateDefectCommutationGradInterp	(	const PressureType &	p,
		const PressureGradientType &	grad_p,
		const size_t	deg_quad_interpolate
	)		const

Evaluate the defect of commutation G_h (Igrad p) = Icurl (grad p) by computing the norm of the difference (which should be 0).

Parameters

p	The pressure
grad_p	Gradient of the pressure
deg_quad_interpolate	Degree of quadrature for computing the interpolate of grad p (the interpolate of p is computed using the highest possible degree)

◆ nbSCDOFs()

size_t HArDCore3D::Stokes::nbSCDOFs ( ) const

inline

Returns the number of statically condensed DOFs (both velocity and pressure)

◆ nbSCDOFs_p()

size_t HArDCore3D::Stokes::nbSCDOFs_p ( ) const

inline

Returns the number of statically condensed DOFs (pressure)

◆ nbSCDOFs_u()

size_t HArDCore3D::Stokes::nbSCDOFs_u ( ) const

inline

Returns the number of statically condensed DOFs (velocity)

◆ scMatrix()

const SystemMatrixType & HArDCore3D::Stokes::scMatrix ( ) const

inline

Returns the static condensation recovery operator.

◆ scVector()

Eigen::VectorXd & HArDCore3D::Stokes::scVector ( )

inline

Returns the static condensation rhs.

◆ sizeSystem()

size_t HArDCore3D::Stokes::sizeSystem ( ) const

inline

Returns the size of the statically condensed system (with Lagrange multiplier)

◆ stabilizationParameter() [1/2]

double & HArDCore3D::Stokes::stabilizationParameter ( )

inline

Returns the stabilization parameter.

◆ stabilizationParameter() [2/2]

const double & HArDCore3D::Stokes::stabilizationParameter ( ) const

inline

Returns the stabilization parameter.

◆ Stokes()

Stokes::Stokes	(	const DDRCore &	ddrcore,
		bool	use_threads,
		std::ostream &	output = `std::cout`
	)

Constructor.

Parameters

ddrcore	Core for the DDR space sequence
use_threads	True for parallel execution, false for sequential execution
output	Output stream to print status messages

◆ StokesNorms() [1/2]

HArDCore3D::StokesNorms::StokesNorms ( )

inline

◆ StokesNorms() [2/2]

HArDCore3D::StokesNorms::StokesNorms	(	double	norm_u,
		double	norm_curl_u,
		double	norm_p,
		double	norm_grad_p
	)

inline

Constructor.

◆ systemMatrix() [1/2]

SystemMatrixType & HArDCore3D::Stokes::systemMatrix ( )

inline

Returns the linear system matrix.

◆ systemMatrix() [2/2]

const SystemMatrixType & HArDCore3D::Stokes::systemMatrix ( ) const

inline

Returns the linear system matrix.

◆ systemVector() [1/2]

Eigen::VectorXd & HArDCore3D::Stokes::systemVector ( )

inline

Returns the linear system right-hand side vector.

◆ systemVector() [2/2]

const Eigen::VectorXd & HArDCore3D::Stokes::systemVector ( ) const

inline

Returns the linear system right-hand side vector.

◆ xCurl()

const XCurl & HArDCore3D::Stokes::xCurl ( ) const

inline

Returns the space XCurl.

◆ xDiv()

const XDiv & HArDCore3D::Stokes::xDiv ( ) const

inline

Returns the space XDiv.

◆ xGrad()

const XGrad & HArDCore3D::Stokes::xGrad ( ) const

inline

Returns the space XGrad.

Variable Documentation

◆ curl_u

double HArDCore3D::StokesNorms::curl_u

Norm of curl of velocity (vorticity)

◆ field_curl_u

Stokes::VorticityType HArDCore3D::field_curl_u

static

Initial value:

= [](const Eigen::Vector3d & x) -> Eigen::Vector3d {
                double cux = (PI*PI)*cos(PI*x(0))*cos(PI*x(2))*pow(sin(PI*x(0)),2.0)*pow(sin(PI*x(1)),3.0)*pow(sin(PI*x(2)),2.0)*9.0;
                double cuy = (PI*PI)*cos(PI*x(1))*cos(PI*x(2))*pow(sin(PI*x(0)),3.0)*pow(sin(PI*x(1)),2.0)*pow(sin(PI*x(2)),2.0)*9.0;
                double cuz = (PI*PI)*pow(sin(PI*x(0)),3.0)*pow(sin(PI*x(1)),3.0)*pow(sin(PI*x(2)),3.0)*6.0-(PI*PI)*pow(cos(PI*x(0)),2.0)*sin(PI*x(0))*pow(sin(PI*x(1)),3.0)*pow(sin(PI*x(2)),3.0)*6.0-(PI*PI)*pow(cos(PI*x(1)),2.0)*pow(sin(PI*x(0)),3.0)*sin(PI*x(1))*pow(sin(PI*x(2)),3.0)*6.0;
 
                return Eigen::Vector3d(cux, cuy, cuz);
                }

◆ field_f

Stokes::ForcingTermType HArDCore3D::field_f

static

Initial value:

= [](const Eigen::Vector3d & x) -> Eigen::Vector3d {
              double fx = (PI*PI*PI)*pow(cos(PI*x(1)),3.0)*pow(sin(PI*x(0)),3.0)*pow(sin(PI*x(2)),3.0)*-6.0+(PI*PI*PI)*cos(PI*x(1))*pow(sin(PI*x(0)),3.0)*pow(sin(PI*x(1)),2.0)*pow(sin(PI*x(2)),3.0)*3.9E+1-(PI*PI*PI)*pow(cos(PI*x(0)),2.0)*cos(PI*x(1))*sin(PI*x(0))*pow(sin(PI*x(1)),2.0)*pow(sin(PI*x(2)),3.0)*1.8E+1-(PI*PI*PI)*cos(PI*x(1))*pow(cos(PI*x(2)),2.0)*pow(sin(PI*x(0)),3.0)*pow(sin(PI*x(1)),2.0)*sin(PI*x(2))*1.8E+1;
 
              double fy = (PI*PI*PI)*pow(cos(PI*x(0)),3.0)*pow(sin(PI*x(1)),3.0)*pow(sin(PI*x(2)),3.0)*6.0-(PI*PI*PI)*cos(PI*x(0))*pow(sin(PI*x(0)),2.0)*pow(sin(PI*x(1)),3.0)*pow(sin(PI*x(2)),3.0)*3.9E+1+(PI*PI*PI)*cos(PI*x(0))*pow(cos(PI*x(1)),2.0)*pow(sin(PI*x(0)),2.0)*sin(PI*x(1))*pow(sin(PI*x(2)),3.0)*1.8E+1+(PI*PI*PI)*cos(PI*x(0))*pow(cos(PI*x(2)),2.0)*pow(sin(PI*x(0)),2.0)*pow(sin(PI*x(1)),3.0)*sin(PI*x(2))*1.8E+1;
 
              double fz = 0.;
               return Eigen::Vector3d(fx, fy, fz);
             }

◆ field_grad_p

Stokes::PressureGradientType HArDCore3D::field_grad_p

static

Initial value:

= [](const Eigen::Vector3d & x) -> Eigen::Vector3d {
                    return Eigen::Vector3d(0., 0., 0.);
                  }

◆ field_nu

Stokes::ViscosityType HArDCore3D::field_nu = Stokes::ViscosityType(1.)

static

◆ field_p

Stokes::PressureType HArDCore3D::field_p

static

Initial value:

= [](const Eigen::Vector3d & x) -> double {
               return 0.;
             }

◆ field_u

Stokes::VelocityType HArDCore3D::field_u

static

Initial value:

= [](const Eigen::Vector3d & x) -> Eigen::Vector3d {
               double ux = PI*cos(PI*x(1))*pow(sin(PI*x(0)),3.0)*pow(sin(PI*x(1)),2.0)*pow(sin(PI*x(2)),3.0)*3.0;
               double uy = PI*cos(PI*x(0))*pow(sin(PI*x(0)),2.0)*pow(sin(PI*x(1)),3.0)*pow(sin(PI*x(2)),3.0)*(-3.0);
               return Eigen::Vector3d(ux, uy, 0.);
             }

◆ grad_p

double HArDCore3D::StokesNorms::grad_p

Norm of grad of pressure.

◆ hcurl_u

double HArDCore3D::StokesNorms::hcurl_u

Hcurl norm of u.

◆ hgrad_p

double HArDCore3D::StokesNorms::hgrad_p

Hgrad norm of p.

◆ linear_curl_u

Stokes::VorticityType HArDCore3D::linear_curl_u

static

Initial value:

= [](const Eigen::Vector3d & x) -> Eigen::Vector3d {
                    return Eigen::Vector3d::Zero();
                  }

◆ linear_f

Stokes::ForcingTermType HArDCore3D::linear_f

static

Initial value:

= [](const Eigen::Vector3d & x) -> Eigen::Vector3d {
               return linear_grad_p(x);
             }

◆ linear_grad_p

Stokes::PressureGradientType HArDCore3D::linear_grad_p

static

Initial value:

= [](const Eigen::Vector3d & x) -> Eigen::Vector3d {
                    return pressure_scaling * Eigen::Vector3d(1., -1., 0.);
                  }

◆ linear_nu

Stokes::ViscosityType HArDCore3D::linear_nu = Stokes::ViscosityType(1.)

static

◆ linear_p

Stokes::PressureType HArDCore3D::linear_p

static

Initial value:

= [](const Eigen::Vector3d & x) -> double {
               return pressure_scaling * (x(0) - x(1));
             }

◆ linear_u

Stokes::VelocityType HArDCore3D::linear_u

static

Initial value:

= [](const Eigen::Vector3d & x) -> Eigen::Vector3d {
               return Eigen::Vector3d(x(0), -x(1), 0.);
             }

◆ p

double HArDCore3D::StokesNorms::p

Norm of pressure.

◆ PI

const double HArDCore3D::PI = boost::math::constants::pi<double>()

static

◆ trigonometric_curl_u

Stokes::VorticityType HArDCore3D::trigonometric_curl_u

static

Initial value:

= [](const Eigen::Vector3d & x) -> Eigen::Vector3d {
                      return 3. * PI * Eigen::Vector3d(
                                                      cos(2. * PI * x(0)) * sin(2. * PI * x(1)) * sin(2. * PI * x(2)),
                                                      -sin(2. * PI * x(0)) * cos(2. * PI * x(1)) * sin(2. * PI * x(2)),
                                                      0.
                                                      );
                         }

◆ trigonometric_f

Stokes::ForcingTermType HArDCore3D::trigonometric_f

static

Initial value:

= [](const Eigen::Vector3d & x) -> Eigen::Vector3d {
                      return 6. * std::pow(PI, 2) * Eigen::Vector3d(
                                                                    sin(2. * PI * x(0)) * cos(2. * PI * x(1)) * cos(2. * PI * x(2)),
                                                                    cos(2. * PI * x(0)) * sin(2. * PI * x(1)) * cos(2. * PI * x(2)),
                                                                    -2. * cos(2. * PI * x(0)) * cos(2. * PI * x(1)) * sin(2. * PI * x(2))
                                                                    )
                        + trigonometric_grad_p(x);
                    }

◆ trigonometric_grad_p

Stokes::PressureGradientType HArDCore3D::trigonometric_grad_p

static

Initial value:

= [](const Eigen::Vector3d & x) -> Eigen::Vector3d {
                      return pressure_scaling * 2. * PI * Eigen::Vector3d(
                                                      cos(2. * PI * x(0)) * sin(2. * PI * x(1)) * sin(2. * PI * x(2)),
                                                      sin(2. * PI * x(0)) * cos(2. * PI * x(1)) * sin(2. * PI * x(2)),
                                                      sin(2. * PI * x(0)) * sin(2. * PI * x(1)) * cos(2. * PI * x(2))
                                                      );
                         }

◆ trigonometric_nu

Stokes::ViscosityType HArDCore3D::trigonometric_nu = Stokes::ViscosityType(1.)

static

◆ trigonometric_p

Stokes::PressureType HArDCore3D::trigonometric_p

static

Initial value:

= [](const Eigen::Vector3d & x) -> double {
                      return pressure_scaling * sin(2. * PI * x(0)) * sin(2. * PI * x(1)) * sin(2. * PI * x(2));
                    }

◆ trigonometric_u

Stokes::VelocityType HArDCore3D::trigonometric_u

static

Initial value:

= [](const Eigen::Vector3d & x) -> Eigen::Vector3d {
                      return Eigen::Vector3d(
                                             0.5 * sin(2. * PI * x(0)) * cos(2. * PI * x(1)) * cos(2. * PI * x(2)),
                                             0.5 * cos(2. * PI * x(0)) * sin(2. * PI * x(1)) * cos(2. * PI * x(2)),
                                             -cos(2. * PI * x(0)) * cos(2. * PI * x(1)) * sin(2. * PI * x(2))
                                             );
                    }

◆ u

double HArDCore3D::StokesNorms::u

Norm of velocity.

Classes

Typedefs

Functions

Variables

Detailed Description

Typedef Documentation

◆ ForcingTermType

◆ PressureGradientType

◆ PressureType

◆ SystemMatrixType

◆ VelocityType

◆ ViscosityType

◆ VorticityType

Function Documentation

◆ assembleLinearSystem()

◆ computeContinuousErrorsNorms()

◆ computeStokesNorms()

◆ dimension()

◆ evaluateDefectCommutationGradInterp()

◆ nbSCDOFs()

◆ nbSCDOFs_p()

◆ nbSCDOFs_u()

◆ scMatrix()

◆ scVector()

◆ sizeSystem()

◆ stabilizationParameter() [1/2]

◆ stabilizationParameter() [2/2]

◆ Stokes()

◆ StokesNorms() [1/2]

◆ StokesNorms() [2/2]

◆ systemMatrix() [1/2]

◆ systemMatrix() [2/2]

◆ systemVector() [1/2]

◆ systemVector() [2/2]

◆ xCurl()

◆ xDiv()

◆ xGrad()

Variable Documentation

◆ curl_u

◆ field_curl_u

◆ field_f

◆ field_grad_p

◆ field_nu

◆ field_p

◆ field_u

◆ grad_p

◆ hcurl_u

◆ hgrad_p

◆ linear_curl_u

◆ linear_f

◆ linear_grad_p

◆ linear_nu

◆ linear_p

◆ linear_u

◆ p

◆ PI

◆ trigonometric_curl_u

◆ trigonometric_f

◆ trigonometric_grad_p

◆ trigonometric_nu

◆ trigonometric_p

◆ trigonometric_u

◆ u