Class to solve the linear conservation law equation in 3-D by a MUSCL Finite Volume scheme on tetrahedra. More...

Inheritance diagram for LCL3DT:

Public Types
enum	Method { FIRST_ORDER_METHOD = 0, MULTI_SLOPE_Q_METHOD = 1, MULTI_SLOPE_M_METHOD = 2, M_FORCE_WORD = 0xFFFFFFFF }
	Enumeration for flux choice. More...
enum	Limiter { MINMOD_LIMITER = 0, VANLEER_LIMITER = 1, SUPERBEE_LIMITER = 2, VANALBADA_LIMITER = 3, MAX_LIMITER = 4, L_FORCE_WORD = 0xFFFFFFFF }
	Enumeration of flux limiting methods. More...
enum	SolverType { ROE_SOLVER = 0, VFROE_SOLVER = 1, LF_SOLVER = 2, RUSANOV_SOLVER = 3, HLL_SOLVER = 4, HLLC_SOLVER = 5, MAX_SOLVER = 6, FORCE_WORD = 0xFFFFFFFF }
	Enumeration of various solvers for the Riemann problem. More...
Public Member Functions
	LCL3DT (Mesh &m)
	Constructor using mesh.
	LCL3DT (Mesh &m, Vect< double > &U)
	Constructor using mesh and initial field.
	~LCL3DT ()
	Destructor.
void	setInitialCondition (Vect< double > &u)
	Set elementwise initial condition.
void	setInitialCondition (double u)
	Set a constant initial condition.
void	setReconstruction ()
	Reconstruct flux using Muscl scheme.
double	runOneTimeStep ()
	Run one time step.
void	setBC (double u)
	Set Dirichlet boundary condition. Assign a constant value u to all boundary sides.
void	setBC (const Side &sd, double u)
	Set Dirichlet boundary condition.
void	setBC (int code, double u)
	Set Dirichlet boundary condition.
void	setVelocity (const Vect< double > &v)
	Set convection velocity.
void	setVelocity (const LocalVect< double, 3 > &v)
	Set (constant) convection velocity.
void	setReferenceLength (double dx)
	Assign reference length value.
double	getReferenceLength () const
	Return reference length.
void	Forward (const Vect< double > &Flux, Vect< double > &Field)
	Computation of the primal variable n->n+1.
double	getMinimumFaceArea () const
	Return minimum area of faces in the mesh.
double	getMinimumElementVolume () const
	Return minimum volume of elements in the mesh.
double	getMaximumFaceArea () const
	Return maximum area of faces in the mesh.
double	getMaximumElementVolume () const
	Return maximum volume of elements in the mesh.
double	getMeanFaceArea () const
	Return mean area of faces in the mesh.
double	getMeanElementVolume () const
	Return mean volume of elements in the mesh.
double	getMinimumEdgeLength () const
	Return minimum length of edges in the mesh.
double	getMinimumVolumebyArea () const
	Return minimum volume by area in the mesh.
double	getMaximumEdgeLength () const
	Return maximum length of edges in the mesh.
double	getTauLim () const
	Return value of tau lim.
double	getComega () const
	Return value of Comega.
void	setbetalim (double bl)
	Assign value of beta lim.
void	setTimeStep (double dt)
	Assign time step value.
double	getTimeStep () const
	Return time step value.
void	setCFL (double CFL)
	Assign CFL value.
double	getCFL () const
	Return CFL value.
Mesh &	getMesh () const
	Return reference to Mesh instance.
void	setVerbose (int v)
	Set verbosity parameter.
bool	setReconstruction (const Vect< double > &U, Vect< double > &LU, Vect< double > &RU, size_t dof)
	Function to reconstruct by the Muscl method.
void	setMethod (const Method &s)
	Choose a flux solver.
void	setSolidZoneCode (int c)
	Choose a code for solid zone.
bool	getSolidZone () const
	Return flag for presence of solid zones.
int	getSolidZoneCode () const
	Return code of solid zone, 0 if this one is not present.
void	setLimiter (Limiter l)
	Choose a flux limiter.

Detailed Description

Class to solve the linear conservation law equation in 3-D by a MUSCL Finite Volume scheme on tetrahedra.

Member Enumeration Documentation

enum Method [inherited]

Enumeration for flux choice.

Enumerator:

FIRST_ORDER_METHOD	First Order upwind method
MULTI_SLOPE_Q_METHOD	Multislope Q method
MULTI_SLOPE_M_METHOD	Multislope M method
M_FORCE_WORD	for compatibity reasons, gcc may return a warning anyway

enum Limiter [inherited]

Enumeration of flux limiting methods.

Enumerator:

MINMOD_LIMITER	MinMod limiter
VANLEER_LIMITER	Van Leer limiter
SUPERBEE_LIMITER	Superbee limiter
VANALBADA_LIMITER	Van Albada limiter
MAX_LIMITER	Max limiter
L_FORCE_WORD	for compatibity reason, gcc may return a warning anyway

enum SolverType [inherited]

Enumeration of various solvers for the Riemann problem.

Enumerator:

ROE_SOLVER	Roe solver
VFROE_SOLVER	Finite Volume Roe solver
LF_SOLVER	LF solver
RUSANOV_SOLVER	Rusanov solver
HLL_SOLVER	HLL solver
HLLC_SOLVER	HLLC solver
MAX_SOLVER	Max solver
FORCE_WORD	for compatibity reason, gcc may return a warning anyway

Constructor & Destructor Documentation

LCL3DT	(	Mesh &	m,
		Vect< double > &	U
	)

Constructor using mesh and initial field.

Parameters:

[in]	m	Reference to Mesh instance
[in]	U	Vector containing initial (elementwise) solution

Member Function Documentation

void setInitialCondition ( Vect< double > & u )

Set elementwise initial condition.

Parameters:

[in] u Vect instance containing initial condition values

void setInitialCondition ( double u )

Set a constant initial condition.

Parameters:

[in] u Value of initial condition to assign to all elements

void setBC	(	const Side &	sd,
		double	u
	)

Set Dirichlet boundary condition.

Assign a constant value to a side

Parameters:

[in]	sd	Side to which value is prescibed
[in]	u	Value to prescribe

void setBC	(	int	code,
		double	u
	)

Set Dirichlet boundary condition.

Assign a constant value sides with a given code

Parameters:

[in]	code	Code of sides to which value is prescibed
[in]	u	Value to prescribe

void setVelocity ( const Vect< double > & v )

Set convection velocity.

Parameters:

[in] v Vect instance containing velocity

void setVelocity ( const LocalVect< double, 3 > & v )

Set (constant) convection velocity.

Parameters:

[in] v Vector containing constant velocity to prescribe

void Forward	(	const Vect< double > &	Flux,
		Vect< double > &	Field
	)

Computation of the primal variable n->n+1.

Vector Flux contains elementwise fluxes issued from the Riemann problem, calculated with, as left element, getNeighborElement(1) and right element getNeighborElement(2) if getNeighborElement(2) doesn't exist, we are on a boundary and we prescribe a symmetry condition

void setTimeStep ( double dt ) [inherited]

Assign time step value.

Parameters:

[in] dt Time step value

Reimplemented in ICPG3DT.

void setCFL ( double CFL ) [inherited]

Assign CFL value.

Parameters:

[in] CFL Value of CFL

Reimplemented in ICPG3DT.

void setVerbose ( int v ) [inherited]

Set verbosity parameter.

Parameters:

[in] v Value of verbosity parameter

bool setReconstruction	(	const Vect< double > &	U,
		Vect< double > &	LU,
		Vect< double > &	RU,
		size_t	dof
	)		`[inherited]`

Function to reconstruct by the Muscl method.

Parameters:

[in]	U	Field to reconstruct
[out]	LU	Left gradient vector
[out]	RU	Right gradient vector
[in]	dof	Label of dof to reconstruct

void setMethod ( const Method & s ) [inherited]

Choose a flux solver.

Parameters:

[in] s Solver to choose

void setLimiter ( Limiter l ) [inherited]

Choose a flux limiter.

Parameters:

[in] l Limiter to choose

Public Types

Public Member Functions

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation