Model class for Extended Navier-Stokes-Korteweg mathematical system of equations (single phase). More...

#include <ModExtendedNavierStokesKorteweg1D.h>

Inherits Model.

Collaboration diagram for ModExtendedNavierStokesKorteweg1D:

Public Member Functions
	ModExtendedNavierStokesKorteweg1D (const int &numbTransports)
	ExtendedNavierStokesKorteweg1D model constructor.
	ModExtendedNavierStokesKorteweg1D (const int &numbTransports, const double &lambda, const double &beta, const double &tau, const double &eps)
	ExtendedNavierStokesKorteweg1D model constructor.
	ModExtendedNavierStokesKorteweg1D (const std::string &name, const int &numbTransports, const double &lambda, const double &beta, const double &tau, const double &eps)
	ExtendedNavierStokesKorteweg1D model constructor.
	~ModExtendedNavierStokesKorteweg1D () override
void	allocateCons (Flux **cons) override
	Allocate conservative variable arrays.
void	allocatePhase (Phase **phase) override
	Instanciate fluid phase variable.
void	allocateMixture (Mixture **mixture) override
	Instanciate mixture variable.
void	fulfillState (Phase *, Mixture ) override
void	fulfillStateResume (Phase *, Mixture ) override
void	initializeAugmentedVariables (Cell *cell) override
	Complete the augmented variables (such as the ones of Navier-Stokes-Korteweg model).
void	solveRiemannIntern (Cell &cellLeft, Cell &cellRight, const double &dxLeft, const double &dxRight, double &dtMax, std::vector< double > &=DEFAULT_VEC_INTERFACE_DATA) const override
	Cell to cell Riemann solver.
void	solveRiemannWall (Cell &cellLeft, const double &dxLeft, double &dtMax, std::vector< double > &boundData) const override
	Wall half Riemann solver.
void	solveRiemannNullFlux () const override
	No flux half Riemann solver (return null flux to use with 1D geometry with smooth varying cross section).
void	reverseProjection (const Coord normal, const Coord tangent, const Coord binormal) const override
	Flux reverse projection in the absolute Cartesian coordinate system.
double	selectScalar (Phase *phases, Mixture , Transport *transports, Variable nameVariable, int num=0) const override
	Select a specific scalar variable.
const double &	getSM () override
	Return the local fluid velocity.
const Coord &	getVelocity (const Cell *cell) const override
	Return the fluid velocity of the corresponding cell.
Coord &	getVelocity (Cell *cell) override
const std::string &	whoAmI () const override
virtual void	computeMaxWaveSpeed (const double &rhoL, const double &rhoR, const double &uL, const double &uR, const double &rhoSurEtaL, const double &rhoSurEtaR, const double &vecPxL, const double &vecPxR, const double &kappaL, const double &kappaR, const double &kappaPrimeL, const double &kappaPrimeR, const double &kappaSecondL, const double &kappaSecondR, const double &c2L, const double &c2R, double &sL, double &sR) const
	Compute maximal wave speed using Davis approximation and eigenvalues of hyperbolic equations.
Public Member Functions inherited from Model
	Model (const std::string &name, const int &numbTransports)
	Generic model constructor.
virtual	~Model ()
virtual void	allocatePhaseSolid (Phase **)
	Instanciate solid phase variable.
virtual void	allocatePhaseGradient (GradPhase **)
	Instanciate gradient phase variable.
virtual void	allocatePhaseSolidGradient (GradPhase **)
	Instanciate gradient solid-phase variable.
virtual void	allocateMixtureGradient (GradMixture **)
	Instanciate gradient mixture variable.
void	allocateEos (Cell &cell) const
	Associate equations of state.
void	initializeRelaxation (Cell *cell) const
	Initialize the theoritical critical pressure of the fluid (only required for PTMu relax).
virtual void	solveRiemannInternMRF (Cell &, Cell &, const double &, const double &, double &, const Coord &, const Coord &, const Coord &, const Coord &, const Coord &) const
	Cell to cell Riemann solver + compute fluxBuffMRF for MRF interface.
virtual void	solveRiemannSymmetry (Cell &cellLeft, const double &dxLeft, double &dtMax, std::vector< double > &boundData) const
	Symmetry half Riemann solver.
virtual void	solveRiemannMovingWall (Cell &, const double &, double &, double, std::vector< double > &) const
	Moving wall half Riemann solver.
virtual void	solveRiemannInletTank (Cell &, const double &, double &, const double , const double , const double &, const double &, std::vector< double > &) const
	Inlet tank half Riemann solver.
virtual void	solveRiemannInletInjStagState (Cell &, const double &, double &, const double, const double , const double , const double *, std::vector< double > &) const
	Inlet injection using stagnation state half Riemann solver.
virtual void	solveRiemannInletInjTemp (Cell &, const double &, double &, const double, const double , const double , std::vector< double > &) const
	Inlet injection using temperature half Riemann solver.
virtual void	solveRiemannOutletPressure (Cell &, const double &, double &, const double, std::vector< double > &) const
	Outlet at imposed pressure half Riemann solver.
virtual void	solveRiemannOutletMassflow (Cell &, const double &, double &, const double, std::vector< double > &) const
	Outlet at imposed massflow half Riemann solver.
virtual void	solveRiemannTransportIntern (Cell &, Cell &)
	Cell to cell Riemann solver for transport equations.
virtual void	solveRiemannTransportWall ()
	Wall half Riemann solver for transport equations.
virtual void	solveRiemannTransportMovingWall (Cell &, double)
	Moving wall half Riemann solver for transport equations.
virtual void	solveRiemannTransportInletInjStagState (Cell &, double *)
	Inlet injection using stagnation state half Riemann solver for transport equations.
virtual void	solveRiemannTransportInletTank (Cell &, double *)
	Inlet tank half Riemann solver for transport equations.
virtual void	solveRiemannTransportOutletPressure (Cell &, double *)
	Outlet at imposed pressure half Riemann solver for transport equations.
void	relaxations (Cell *cell, const double &dt, Prim type=vecPhases) const
	Run all the relaxations.
virtual void	correctionEnergy (Cell *, Prim=vecPhases) const
	Energy correction after relaxations (only for PUEq model).
virtual void	lowMachSoundSpeed (double &, const double &, double &, const double &=Errors::defaultDouble, double &=Tools::uselessDouble) const
virtual void	addNonConsMrfFlux (Phase **)
virtual void	reverseProjectionMrfFlux (const Coord, const Coord, const Coord) const
virtual void	readCompactionParameters (tinyxml2::XMLElement *, std::string, const int &)
	Read compaction parameters.
virtual bool	compaction (const int &) const
	Return if compaction is activated or not.
virtual const double &	getDensityReference (const int &) const
	Return the compaction reference density of phase k.
virtual double	computeEnergyCompaction (const int &, const double &, const double &, const double &) const
	Return the compaction energy of the fluid/solid.
virtual double	computeCompactionPressure_dkappadxi (const int &, const double &, const double &, const double &) const
	Compute the compaction pressure term including the derivative d kappa_k / d xi_k.
virtual double	computeCompactionPressure_dkappadalpha (const int &, const double &, const double &, const double &) const
	Compute the compaction pressure term including the derivative d kappa_k / d alpha_k.
virtual double	computeCompactionPlasticity (const int &, const double &, const double &, const double &, const double &, double &) const
	Compute the compaction plasticity term of the fluid/solid for relaxation.
virtual void	readSolidParameters (tinyxml2::XMLElement *, std::string, const int &)
	Read solid parameters.
virtual const double *	getDensityReference () const
	Return the reference density for each phase.
virtual const double *	getElasticParameterA () const
	Return the elastic parameter a of the one-parameter model for each phase.
virtual const double *	getShearModulus () const
	Return the shear modulus for each phase.
virtual const double *	getChi () const
	Return the chi modulus for each phase.
virtual const double *	getElasticityLimit () const
	Return the limit of elasticity for each phase.
virtual double	computeElasticEnergy (const int &, const Tensor &) const
	Compute the elastic energy of the solid.
virtual void	computeElasticEnergyAndStressTensor (const int &, const Tensor &, const double &, const double &, double &, Tensor &) const
	Compute the elastic energy and stress tensor of the solid.
virtual double	computeSquareLongitudinalWaveSpeed (const int &, const Phase &)
	Compute the square longitudinal wave speed of the solid.
virtual const Coord &	getVectorP (const Cell *) const
virtual Coord &	getVectorP (Cell *)
std::vector< Relaxation * > *	getRelaxations ()
void	printInfo () const
virtual void	setSmoothCrossSection1d (const bool &)
virtual const bool &	isSmoothCrossSection1d () const
const bool &	getLowMach () const
void	setLowMach (const bool &lowMach)
virtual void	setMachRefMin (const double &machRefMin)
virtual void	setPhaseSlope (Phase &, int)
	Set the slope values for the extrapolation within the cell.
virtual void	setSlopeDistance (const double &)
	Set the geometric distance over which the in-cell slopes extrapolate (cell width along the face normal).

Static Private Attributes
static const std::string	NAME = "EXTENDEDNAVIERSTOKESKORTEWEG1D"

Additional Inherited Members
Protected Attributes inherited from Model
std::string	m_name
	Name of the required model.
std::vector< Relaxation * >	m_relaxations
	Vector of relaxation procedure.
bool	m_lowMach
	Low-Mach preconditioning (default: false).
double	m_machRefMin
	Minimum Mach number limit for L-M preconditionning used when local Mach number is below this value (default: 0.01).
bool	m_smoothCrossSection1d
	1D geometry with smooth cross section variation (default: false)

Detailed Description

Model class for Extended Navier-Stokes-Korteweg mathematical system of equations (single phase).

Constructor & Destructor Documentation

◆ ModExtendedNavierStokesKorteweg1D() [1/3]

ModExtendedNavierStokesKorteweg1D::ModExtendedNavierStokesKorteweg1D ( const int & numbTransports )

ExtendedNavierStokesKorteweg1D model constructor.

Parameters

numbTransports number of additional transport equations

Here is the call graph for this function:

◆ ModExtendedNavierStokesKorteweg1D() [2/3]

ModExtendedNavierStokesKorteweg1D::ModExtendedNavierStokesKorteweg1D	(	const int &	numbTransports,
		const double &	lambda,
		const double &	beta,
		const double &	tau,
		const double &	eps )

ExtendedNavierStokesKorteweg1D model constructor.

Parameters

numbTransports	number of additional transport equations
lambda	parameter lambda of Navier-Stokes-Korteweg equations
beta	parameter beta of Navier-Stokes-Korteweg equations
tau	small parameter controlling dissipartion of Navier-Stokes-Korteweg equations
eps	small parameter controlling dispersion of Navier-Stokes-Korteweg equations

Here is the call graph for this function:

◆ ModExtendedNavierStokesKorteweg1D() [3/3]

ModExtendedNavierStokesKorteweg1D::ModExtendedNavierStokesKorteweg1D	(	const std::string &	name,
		const int &	numbTransports,
		const double &	lambda,
		const double &	beta,
		const double &	tau,
		const double &	eps )

ExtendedNavierStokesKorteweg1D model constructor.

Parameters

name	model name
numbTransports	number of additional transport equations
lambda	parameter lambda of Navier-Stokes-Korteweg equations
beta	parameter beta of Navier-Stokes-Korteweg equations
tau	small parameter controlling dissipartion of Navier-Stokes-Korteweg equations
eps	small parameter controlling dispersion of Navier-Stokes-Korteweg equations

Here is the call graph for this function:

◆ ~ModExtendedNavierStokesKorteweg1D()

ModExtendedNavierStokesKorteweg1D::~ModExtendedNavierStokesKorteweg1D ( )

override

Member Function Documentation

◆ allocateCons()

void ModExtendedNavierStokesKorteweg1D::allocateCons ( Flux ** )

overridevirtual

Allocate conservative variable arrays.

Parameters

cons	conservative variable array to allocate

Reimplemented from Model.

◆ allocateMixture()

void ModExtendedNavierStokesKorteweg1D::allocateMixture ( Mixture ** )

overridevirtual

Instanciate mixture variable.

Parameters

mixture mixture to instanciate

Reimplemented from Model.

◆ allocatePhase()

void ModExtendedNavierStokesKorteweg1D::allocatePhase ( Phase ** )

overridevirtual

Instanciate fluid phase variable.

Parameters

phase phase to instanciate

Reimplemented from Model.

◆ computeMaxWaveSpeed()

void ModExtendedNavierStokesKorteweg1D::computeMaxWaveSpeed	(	const double &	rhoL,
		const double &	rhoR,
		const double &	uL,
		const double &	uR,
		const double &	rhoSurEtaL,
		const double &	rhoSurEtaR,
		const double &	vecPxL,
		const double &	vecPxR,
		const double &	kappaL,
		const double &	kappaR,
		const double &	kappaPrimeL,
		const double &	kappaPrimeR,
		const double &	kappaSecondL,
		const double &	kappaSecondR,
		const double &	c2L,
		const double &	c2R,
		double &	sL,
		double &	sR ) const

virtual

Compute maximal wave speed using Davis approximation and eigenvalues of hyperbolic equations.

Parameters

cellLeft	cell left
cellRight	cell right
rhoL	density left
rhoR	density right
uL	velocity left
uR	velocity right
etaL	eta left
etaR	eta right
vecPxL	vecPxL left
vecPxR	vecPxR right
kappaL	kappaL left
kappaR	kappaR right
kappaPrimeL	kappaPrimeL left
kappaPrimeR	kappaPrimeR right
kappaSecondL	kappaSecondL left
kappaSecondR	kappaSecondR right
c2L	squareSoundSpeedL left
c2R	squareSoundSpeedR right
sL	RP left characteristic
sR	RP right characteristic

Here is the caller graph for this function:

◆ fulfillState()

void ModExtendedNavierStokesKorteweg1D::fulfillState	(	Phase **	phases,
		Mixture *	)

overridevirtual

Does nothing for this model

Reimplemented from Model.

Here is the call graph for this function:

◆ fulfillStateResume()

void ModExtendedNavierStokesKorteweg1D::fulfillStateResume	(	Phase **	,
		Mixture *	)

inlineoverridevirtual

Does nothing for this model

Reimplemented from Model.

◆ getSM()

const double & ModExtendedNavierStokesKorteweg1D::getSM ( )

inlineoverridevirtual

Return the local fluid velocity.

Returns: the velocity solution of the local Riemann problem

Reimplemented from Model.

Here is the call graph for this function:

◆ getVelocity() [1/2]

Coord & ModExtendedNavierStokesKorteweg1D::getVelocity ( Cell * cell )

inlineoverridevirtual

Reimplemented from Model.

Here is the call graph for this function:

◆ getVelocity() [2/2]

const Coord & ModExtendedNavierStokesKorteweg1D::getVelocity ( const Cell * ) const

inlineoverridevirtual

Return the fluid velocity of the corresponding cell.

Parameters

cell	pointer to corresponding cell

Returns: velocity

Reimplemented from Model.

Here is the call graph for this function:

Here is the caller graph for this function:

◆ initializeAugmentedVariables()

void ModExtendedNavierStokesKorteweg1D::initializeAugmentedVariables ( Cell * cell )

overridevirtual

Complete the augmented variables (such as the ones of Navier-Stokes-Korteweg model).

Parameters

cell cell

Reimplemented from Model.

Here is the call graph for this function:

◆ reverseProjection()

void ModExtendedNavierStokesKorteweg1D::reverseProjection	(	const Coord	,
		const Coord	,
		const Coord	) const

overridevirtual

Flux reverse projection in the absolute Cartesian coordinate system.

Parameters

normal	normal vector associated to the cell interface
tangent	tangent vector associated to the cell interface
binormal	binormal vector associated to the cell interface

Reimplemented from Model.

◆ selectScalar()

double ModExtendedNavierStokesKorteweg1D::selectScalar	(	Phase **	phases,
		Mixture *	,
		Transport *	transports,
		Variable	nameVariable,
		int	num = 0 ) const

overridevirtual

Select a specific scalar variable.

Parameters

phases	phases array variables
mixture	mixture variables
vecTransports	vector of transport variables
nameVariables	Name of the variable to select
numPhases	Phases number's

Reimplemented from Model.

Here is the call graph for this function:

◆ solveRiemannIntern()

void ModExtendedNavierStokesKorteweg1D::solveRiemannIntern	(	Cell &	,
		Cell &	,
		const double &	,
		const double &	,
		double &	,
		std::vector< double > &	= DEFAULT_VEC_INTERFACE_DATA ) const

overridevirtual

Cell to cell Riemann solver.

Parameters

cellLeft	left cell
cellRight	right cell
dxLeft	left characteristic lenght
dxRight	right characteristic lenght
dtMax	maximum explicit time step
boundData	boundary dataset used for output

Reimplemented from Model.

Here is the call graph for this function:

◆ solveRiemannNullFlux()

void ModExtendedNavierStokesKorteweg1D::solveRiemannNullFlux ( ) const

overridevirtual

No flux half Riemann solver (return null flux to use with 1D geometry with smooth varying cross section).

Reimplemented from Model.

◆ solveRiemannWall()

void ModExtendedNavierStokesKorteweg1D::solveRiemannWall	(	Cell &	,
		const double &	,
		double &	,
		std::vector< double > &	) const

overridevirtual

Wall half Riemann solver.

Parameters

cellLeft	left cell
dxLeft	left characteristic lenght
dtMax	maximum explicit time step
boundData	boundary dataset used for output

Reimplemented from Model.

Here is the call graph for this function:

◆ whoAmI()

const std::string & ModExtendedNavierStokesKorteweg1D::whoAmI ( ) const

inlineoverridevirtual

Reimplemented from Model.

Member Data Documentation

◆ NAME

const std::string ModExtendedNavierStokesKorteweg1D::NAME = "EXTENDEDNAVIERSTOKESKORTEWEG1D"

staticprivate

The documentation for this class was generated from the following files:

/home/runner/work/ECOGEN/ECOGEN/doc_sources/src/Models/ExtendedNavierStokesKorteweg1D/ModExtendedNavierStokesKorteweg1D.h
/home/runner/work/ECOGEN/ECOGEN/doc_sources/src/Models/ExtendedNavierStokesKorteweg1D/ModExtendedNavierStokesKorteweg1D.cpp

Public Member Functions

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ModExtendedNavierStokesKorteweg1D() [1/3]

◆ ModExtendedNavierStokesKorteweg1D() [2/3]

◆ ModExtendedNavierStokesKorteweg1D() [3/3]

◆ ~ModExtendedNavierStokesKorteweg1D()

Member Function Documentation

◆ allocateCons()

◆ allocateMixture()

◆ allocatePhase()

◆ computeMaxWaveSpeed()

◆ fulfillState()

◆ fulfillStateResume()

◆ getSM()

◆ getVelocity() [1/2]

◆ getVelocity() [2/2]

◆ initializeAugmentedVariables()

◆ reverseProjection()

◆ selectScalar()

◆ solveRiemannIntern()

◆ solveRiemannNullFlux()

◆ solveRiemannWall()

◆ whoAmI()

Member Data Documentation

◆ NAME