#include <NavierStokesSolverBlocks.hpp>

Collaboration diagram for NavierStokesSolverBlocks:

Private Member Functions
void	buildGraphs ()
	Build the graphs. More...

void	updateBCHandler (bcPtr_Type &bc)
	Update the bc handler. More...

void	setSolversOptions (const Teuchos::ParameterList &solversOptions)
	Set options preconditioner. More...

void	applyBoundaryConditionsResidual (vector_Type &r_u, const Real &time=0.0)
	Apply Bc to the residual vector. More...

Private Attributes
commPtr_Type	M_comm

dataFile_Type	M_dataFile

std::string	M_uOrder

std::string	M_pOrder

std::shared_ptr< FESpace< mesh_Type, map_Type > >	M_velocityFESpace

std::shared_ptr< FESpace< mesh_Type, map_Type > >	M_pressureFESpace

std::shared_ptr< FESpace< mesh_Type, map_Type > >	M_velocityFESpaceScalar

std::shared_ptr< ETFESpace_velocity >	M_fespaceUETA

std::shared_ptr< ETFESpace_pressure >	M_fespacePETA

std::shared_ptr< ETFESpace_pressure >	M_fespaceUETA_scalar

bool	M_stiffStrain

graphPtr_Type	M_Mu_graph

graphPtr_Type	M_Btranspose_graph

graphPtr_Type	M_B_graph

graphPtr_Type	M_C_graph

graphPtr_Type	M_A_graph

graphPtr_Type	M_F_graph

graphPtr_Type	M_Jacobian_graph

graphPtr_Type	M_Mp_graph

graphPtr_Type	M_Fp_graph

matrixPtr_Type	M_Mu

matrixPtr_Type	M_Btranspose

matrixPtr_Type	M_B

matrixPtr_Type	M_C

matrixPtr_Type	M_A

matrixPtr_Type	M_F

matrixPtr_Type	M_Jacobian

matrixPtr_Type	M_Mp

matrixPtr_Type	M_Fp

matrixPtr_Type	M_block00

matrixPtr_Type	M_block01

matrixPtr_Type	M_block10

matrixPtr_Type	M_block11

matrixPtr_Type	M_block00_noBC

matrixPtr_Type	M_block01_noBC

matrixPtr_Type	M_block10_noBC

matrixPtr_Type	M_block11_noBC

vectorPtr_Type	M_rhs_noBC

vectorPtr_Type	M_forces

std::shared_ptr< LifeV::Operators::NavierStokesOperator >	M_operLoads

vectorPtr_Type	M_uExtrapolated

vectorPtr_Type	M_rhs

vectorPtr_Type	M_rhs_pressure

vectorPtr_Type	M_velocity

vectorPtr_Type	M_pressure

vectorPtr_Type	M_velocityRhs

vectorPtr_Type	M_velocityExtrapolated

std::shared_ptr< map_Type >	M_monolithicMap

vectorPtr_Type	M_solution

vectorPtr_Type	M_residual_u

vectorPtr_Type	M_residual_p

vectorPtr_Type	M_velocity_old_newton

vectorPtr_Type	M_pressure_old_newton

vectorPtr_Type	M_pressure_extrapolated

Displayer	M_displayer
	Displayer to print in parallel (only PID 0 will print) More...

Real	M_alpha
	Reals. More...

Real	M_timeStep

Real	M_time

bool	M_useGraph
	Booleans. More...

bool	M_graphIsBuilt

bool	M_graphPCDisBuilt

std::shared_ptr< LifeV::Operators::NavierStokesOperator >	M_oper

std::shared_ptr< LifeV::Operators::NavierStokesPreconditionerOperator >	M_prec

std::shared_ptr< Operators::InvertibleOperator >	M_invOper

parameterListPtr_Type	M_pListLinSolver

bool	M_fullyImplicit

bcPtr_Type	M_bcPCD

bool	M_steady

Real	M_relativeTolerance

Real	M_absoluteTolerance

Real	M_etaMax

Int	M_nonLinearLineSearch

UInt	M_maxiterNonlinear

std::ofstream	M_out_res

bcPtr_Type	M_bc

bcPtr_Type	M_bc_sol

bcPtr_Type	M_bc_res_essential

bcPtr_Type	M_bc_res_natural

Real	M_density

Real	M_viscosity

std::shared_ptr< Stabilization >	M_stabilization

bool	M_useStabilization

std::string	M_stabilizationType

bool	M_nonconforming

bool	M_imposeWeakBC

UInt	M_flagWeakBC

Real	M_meshSizeWeakBC

bool	M_computeAerodynamicLoads

std::string	M_methodAerodynamicLoads

UInt	M_flagBody

bool	M_penalizeReverseFlow

UInt	M_flagPenalizeReverseFlow

bool	M_solve_blocks

std::shared_ptr< PostProcessingBoundary< mesh_Type > >	M_postProcessing
	Postprocessing class. More...

bool	M_useFastAssembly

std::shared_ptr< FastAssemblerNS >	M_fastAssembler

Real	M_orderBDF

UInt	M_orderVel

matrixPtr_Type	M_block00_weakBC

matrixPtr_Type	M_block01_weakBC

typedef RegionMesh< LinearTetra >	mesh_Type

typedef std::shared_ptr< mesh_Type >	meshPtr_Type

typedef MapEpetra	map_Type

typedef std::shared_ptr< map_Type >	mapPtr_Type

typedef MatrixEpetra< Real >	matrix_Type

typedef std::shared_ptr< matrix_Type >	matrixPtr_Type

typedef VectorEpetra	vector_Type

typedef std::shared_ptr< vector_Type >	vectorPtr_Type

typedef GetPot	dataFile_Type

typedef std::shared_ptr< dataFile_Type >	dataFilePtr_Type

typedef Epetra_Comm	comm_Type

typedef std::shared_ptr< comm_Type >	commPtr_Type

typedef Epetra_FECrsGraph	graph_Type

typedef std::shared_ptr< Epetra_FECrsGraph >	graphPtr_Type

typedef ETFESpace< mesh_Type, map_Type, 3, 3 >	ETFESpace_velocity

typedef ETFESpace< mesh_Type, map_Type, 3, 1 >	ETFESpace_pressure

typedef std::shared_ptr< BCHandler >	bcPtr_Type

typedef Teuchos::ParameterList	parameterList_Type

typedef std::shared_ptr< parameterList_Type >	parameterListPtr_Type

typedef std::shared_ptr< Epetra_Operator >	invOpPtr_Type

typedef LifeV::Preconditioner	basePrec_Type

typedef std::shared_ptr< basePrec_Type >	basePrecPtr_Type

typedef LifeV::PreconditionerIfpack	prec_Type

typedef std::shared_ptr< prec_Type >	precPtr_Type

typedef Teuchos::RCP< Teuchos::ParameterList >	parameterListRCP_Type

Constructor and Destructor
	NavierStokesSolverBlocks (const dataFile_Type dataFile, const commPtr_Type &communicator)

	~NavierStokesSolverBlocks ()

Methods
void	setup (const meshPtr_Type &mesh, const int &id_domain=36)
	Setup for the solver. More...

void	buildSystem ()
	Assemble constant terms. More...

void	updateSystem (const vectorPtr_Type &u_star, const vectorPtr_Type &rhs_velocity)
	Update the convective term and the right hand side. More...

void	iterate (bcPtr_Type &bc, const Real &time)
	Solve the current timestep, provided the BC. More...

void	iterate (bcPtr_Type &bc, const Real &time, const vectorPtr_Type &velocities)
	Solve the current timestep, provided the BC and a vector of velocities (used only for for aorta example) More...

void	iterate_steady ()
	Solve the steady Navier-Stokes equations. More...

void	iterate_nonlinear (const Real &time)
	Solve the Navier-Stokes equations at a certain time. More...

void	evalResidual (vector_Type &residual, const vector_Type &solution, const UInt iter_newton)
	Evaluate the residual of the NS problem. More...

void	solveJac (vector_Type &increment, const vector_Type &residual, const Real linearRelTol)
	Evaluate the residual of the NS problem. More...

void	updateJacobian (const vector_Type &u_k)
	Update the Jacobian matrix, only the term associated with the linearization of the convective term. More...

void	applyBoundaryConditions (bcPtr_Type &bc, const Real &time)
	Apply the BCs semi-implicit case. More...

void	applyBoundaryConditions (bcPtr_Type &bc, const Real &time, const vectorPtr_Type &velocities)
	Apply the BCs semi-implicit case (example aorta) More...

void	applyBoundaryConditionsJacobian (bcPtr_Type &bc)
	Apply the BCs on the Jacobian matrix. More...

void	applyBoundaryConditionsSolution (const Real &time)
	Apply the BCs to the solution. More...

void	updateConvectiveTerm (const vectorPtr_Type &velocity)
	Update the convective term. More...

VectorSmall< 2 >	computeForces (BCHandler &bcHandlerDrag, BCHandler &bcHandlerLift)
	Compute Aerodynamic Forces. More...

void	solveTimeStep ()
	Solve a time step. More...

void	computeForcesNonLinear (vectorPtr_Type &force, const vectorPtr_Type &solution)
	Compute Aerodynamic Forces - nonlinear case. More...

void	evaluateResidual (const vectorPtr_Type &convective_velocity, const vectorPtr_Type &velocity_km1, const vectorPtr_Type &pressure_km1, const vectorPtr_Type &rhs_velocity, vectorPtr_Type &residual)
	Evaluates the fluid residual in FSI simulations. More...

void	evaluateResidual (const vectorPtr_Type &convective_velocity, const vectorPtr_Type &velocity_km1, const vectorPtr_Type &pressure_km1, const vectorPtr_Type &rhs_velocity, vectorPtr_Type &residualVelocity, vectorPtr_Type &residualPressure)
	Evaluates the fluid residual in FSI simulations. More...

void	updateStabilization (const vector_Type &convective_velocity_previous_newton_step, const vector_Type &velocity_previous_newton_step, const vector_Type &pressure_previous_newton_step, const vector_Type &velocity_rhs)
	Update stabilization terms. More...

void	integrateForces (const vectorPtr_Type &velocity, const vectorPtr_Type &pressure)
	Computation of forces. More...

void	preprocessBoundary (const Real &nx, const Real &ny, const Real &nz, BCHandler &bc, Real &Q_hat, const vectorPtr_Type &Phi_h, const UInt flag, vectorPtr_Type &Phi_h_flag, vectorPtr_Type &V_hat_x, vectorPtr_Type &V_hat_y, vectorPtr_Type &V_hat_z)
	Additional method used for pre-processing on non-circular boundaries (used only in example aorta) More...

void	solveLaplacian (const UInt &flag, bcPtr_Type &bc_laplacian, vectorPtr_Type &laplacianSolution)
	Additional method used for pre-processing on non-circular boundaries (used only in example aorta) More...

void	setupPostProc ()
	Additional setup for postprocessing on boundaries. More...

Real	flux (const markerID_Type &flag, const vector_Type &velocity)
	Compute flux through a boundary face of the domain. More...

Real	area (const markerID_Type &flag)
	Compute area of a boundary face of the domain. More...

Real	pres (const markerID_Type &flag, const vector_Type &pressure)
	Compute mean pressure at a boundary face of the domain. More...

Vector	geometricCenter (const markerID_Type &flag)
	Compute center of a boundary face of the domain. More...

Vector	normal (const markerID_Type &flag)
	Compute normal of a boundary face of the domain. More...

Set methods
void	setAlpha (const Real &alpha)
	Set coefficient associated to the time discretization scheme. More...

void	setRhsVelocity (const vectorPtr_Type &vel_rhs)
	Set the rhs vector, velocity component. More...

void	setExtrapolatedPressure (const vectorPtr_Type &pressure_extrapolated)
	Set the extrapolated pressure vector (used by semi-implicit VMS-LES stabilization) More...

void	setExportFineScaleVelocity (ExporterHDF5< mesh_Type > &exporter, const int &numElementsTotal)
	Set if one needs to export the fine scale component (when using LES models) More...

void	setParameters ()
	Set parameters of the solver. More...

void	setTimeStep (const Real &dt)
	Set time step. More...

void	setBoundaryConditions (const bcPtr_Type &bc)
	Set time step. More...

Get methods
const std::shared_ptr< FESpace< mesh_Type, map_Type > > &	uFESpace () const
	Get the velocity FE space. More...

const std::shared_ptr< FESpace< mesh_Type, map_Type > > &	pFESpace () const
	Get the velocity FE space. More...

const std::shared_ptr< FESpace< mesh_Type, map_Type > > &	uFESpace_scalar () const
	Get the velocity FE space. More...

void	updateVelocity (vectorPtr_Type &velocity)
	Get the velocity vector. More...

void	updatePressure (vectorPtr_Type &pressure)
	Get the pressure vector. More...

matrixPtr_Type const &	getF () const
	Get the F block. More...

matrixPtr_Type const &	getJacobian () const
	Get the jacobian block. More...

matrixPtr_Type const &	getBtranspose () const
	Get the Btranspose block. More...

matrixPtr_Type const &	getB () const
	Get the B block. More...

vectorPtr_Type const &	getRhs () const
	Get the rhs. More...

vectorPtr_Type const &	getRhsPressure () const
	Get the rhs vector (pressure part) More...

matrixPtr_Type const &	Mu () const

Real	density () const

Real	viscosity () const

matrixPtr_Type const &	block00 () const
	Get the (0,0) block. More...

matrixPtr_Type const &	block01 () const
	Get the (0,1) block. More...

matrixPtr_Type const &	block10 () const
	Get the (1,0) block. More...

matrixPtr_Type const &	block11 () const
	Get the (1,1) More...

bool	useStabilization () const
	Get if using a stabilization. More...

void	assembleInterfaceMass (matrixPtr_Type &mass_interface, const mapPtr_Type &interface_map, markerID_Type interfaceFlag, const vectorPtr_Type &numerationInterface, const UInt &offset)

Displayer const &	getDisplayer () const
	Get the displayer. More...

matrixPtr_Type	block00_noBC () const
	Get the (0,0) block without BCs. More...

matrixPtr_Type	block01_noBC () const
	Get the (0,1) block without BCs. More...

vectorPtr_Type	rhs_noBC () const
	Get the rhs without bcs applied. More...

vectorPtr_Type	getForces () const
	Get the forces. More...

Detailed Description

Definition at line 108 of file NavierStokesSolverBlocks.hpp.

Member Typedef Documentation

◆ mesh_Type

typedef RegionMesh<LinearTetra> mesh_Type

Definition at line 116 of file NavierStokesSolverBlocks.hpp.

◆ meshPtr_Type

typedef std::shared_ptr<mesh_Type> meshPtr_Type

Definition at line 117 of file NavierStokesSolverBlocks.hpp.

◆ map_Type

typedef MapEpetra map_Type

Definition at line 119 of file NavierStokesSolverBlocks.hpp.

◆ mapPtr_Type

typedef std::shared_ptr<map_Type> mapPtr_Type

Definition at line 120 of file NavierStokesSolverBlocks.hpp.

◆ matrix_Type

typedef MatrixEpetra<Real> matrix_Type

Definition at line 122 of file NavierStokesSolverBlocks.hpp.

◆ matrixPtr_Type

typedef std::shared_ptr<matrix_Type> matrixPtr_Type

Definition at line 123 of file NavierStokesSolverBlocks.hpp.

◆ vector_Type

typedef VectorEpetra vector_Type

Definition at line 125 of file NavierStokesSolverBlocks.hpp.

◆ vectorPtr_Type

typedef std::shared_ptr<vector_Type> vectorPtr_Type

Definition at line 126 of file NavierStokesSolverBlocks.hpp.

◆ dataFile_Type

typedef GetPot dataFile_Type

Definition at line 128 of file NavierStokesSolverBlocks.hpp.

◆ dataFilePtr_Type

typedef std::shared_ptr<dataFile_Type> dataFilePtr_Type

Definition at line 129 of file NavierStokesSolverBlocks.hpp.

◆ comm_Type

typedef Epetra_Comm comm_Type

Definition at line 131 of file NavierStokesSolverBlocks.hpp.

◆ commPtr_Type

typedef std::shared_ptr< comm_Type > commPtr_Type

Definition at line 132 of file NavierStokesSolverBlocks.hpp.

◆ graph_Type

typedef Epetra_FECrsGraph graph_Type

Definition at line 134 of file NavierStokesSolverBlocks.hpp.

◆ graphPtr_Type

typedef std::shared_ptr<Epetra_FECrsGraph> graphPtr_Type

Definition at line 135 of file NavierStokesSolverBlocks.hpp.

◆ ETFESpace_velocity

typedef ETFESpace<mesh_Type, map_Type, 3, 3 > ETFESpace_velocity

Definition at line 137 of file NavierStokesSolverBlocks.hpp.

◆ ETFESpace_pressure

typedef ETFESpace<mesh_Type, map_Type, 3, 1 > ETFESpace_pressure

Definition at line 138 of file NavierStokesSolverBlocks.hpp.

◆ bcPtr_Type

typedef std::shared_ptr<BCHandler> bcPtr_Type

Definition at line 140 of file NavierStokesSolverBlocks.hpp.

◆ parameterList_Type

typedef Teuchos::ParameterList parameterList_Type

Definition at line 142 of file NavierStokesSolverBlocks.hpp.

◆ parameterListPtr_Type

typedef std::shared_ptr<parameterList_Type> parameterListPtr_Type

Definition at line 143 of file NavierStokesSolverBlocks.hpp.

◆ invOpPtr_Type

typedef std::shared_ptr<Epetra_Operator> invOpPtr_Type

Definition at line 145 of file NavierStokesSolverBlocks.hpp.

◆ basePrec_Type

typedef LifeV::Preconditioner basePrec_Type

Definition at line 147 of file NavierStokesSolverBlocks.hpp.

◆ basePrecPtr_Type

typedef std::shared_ptr<basePrec_Type> basePrecPtr_Type

Definition at line 148 of file NavierStokesSolverBlocks.hpp.

◆ prec_Type

typedef LifeV::PreconditionerIfpack prec_Type

Definition at line 149 of file NavierStokesSolverBlocks.hpp.

◆ precPtr_Type

typedef std::shared_ptr<prec_Type> precPtr_Type

Definition at line 150 of file NavierStokesSolverBlocks.hpp.

◆ parameterListRCP_Type

typedef Teuchos::RCP< Teuchos::ParameterList > parameterListRCP_Type

Definition at line 151 of file NavierStokesSolverBlocks.hpp.

Constructor & Destructor Documentation

◆ NavierStokesSolverBlocks()

NavierStokesSolverBlocks	(	const dataFile_Type	dataFile,
		const commPtr_Type &	communicator
	)

Definition at line 35 of file NavierStokesSolverBlocks.cpp.

◆ ~NavierStokesSolverBlocks()

~NavierStokesSolverBlocks ( )

Definition at line 64 of file NavierStokesSolverBlocks.cpp.

Member Function Documentation

◆ setup()

void setup	(	const meshPtr_Type &	mesh,
		const int &	id_domain = `36`
	)

Setup for the solver.

Parameters

mesh	mesh used
id_domain	used just in case of multiple fluids

Definition at line 77 of file NavierStokesSolverBlocks.cpp.

◆ buildSystem()

void buildSystem ( )

Assemble constant terms.

Definition at line 356 of file NavierStokesSolverBlocks.cpp.

Here is the caller graph for this function:

◆ updateSystem()

void updateSystem	(	const vectorPtr_Type &	u_star,
		const vectorPtr_Type &	rhs_velocity
	)

Update the convective term and the right hand side.

Parameters

u_star	extrapolate velocity
rhs_velocity	right hand side, velocity block

Definition at line 495 of file NavierStokesSolverBlocks.cpp.

Here is the caller graph for this function:

◆ iterate() [1/2]

void iterate	(	bcPtr_Type &	bc,
		const Real &	time
	)

Solve the current timestep, provided the BC.

Parameters

bc	boundary conditions
time	current time

Definition at line 726 of file NavierStokesSolverBlocks.cpp.

◆ iterate() [2/2]

void iterate	(	bcPtr_Type &	bc,
		const Real &	time,
		const vectorPtr_Type &	velocities
	)

Solve the current timestep, provided the BC and a vector of velocities (used only for for aorta example)

Parameters

bc	boundary conditions
time	current time
velocities	vector of velocity (used for imposing velocity on noncircular inflows)

Definition at line 713 of file NavierStokesSolverBlocks.cpp.

◆ iterate_steady()

void iterate_steady ( )

Solve the steady Navier-Stokes equations.

Definition at line 899 of file NavierStokesSolverBlocks.cpp.

Here is the caller graph for this function:

◆ iterate_nonlinear()

void iterate_nonlinear ( const Real & time )

Solve the Navier-Stokes equations at a certain time.

Parameters

time	current time

Definition at line 873 of file NavierStokesSolverBlocks.cpp.

Here is the caller graph for this function:

◆ evalResidual()

void evalResidual	(	vector_Type &	residual,
		const vector_Type &	solution,
		const UInt	iter_newton
	)

Evaluate the residual of the NS problem.

Parameters

residual	residual vector
solution	solution vector
iter_newton	current newton iteration

Definition at line 1168 of file NavierStokesSolverBlocks.cpp.

◆ solveJac()

void solveJac	(	vector_Type &	increment,
		const vector_Type &	residual,
		const Real	linearRelTol
	)

Evaluate the residual of the NS problem.

Parameters

increment	increment vector
residual	residual vector
linearRelTol	tolerance of the linar solver

Definition at line 1508 of file NavierStokesSolverBlocks.cpp.

◆ updateJacobian()

void updateJacobian ( const vector_Type & u_k )

Update the Jacobian matrix, only the term associated with the linearization of the convective term.

Parameters

u_k	velocity previous newton step

Definition at line 1396 of file NavierStokesSolverBlocks.cpp.

◆ applyBoundaryConditions() [1/2]

void applyBoundaryConditions	(	bcPtr_Type &	bc,
		const Real &	time
	)

Apply the BCs semi-implicit case.

Parameters

bc	boundary conditions of the problem
time	current time

Definition at line 667 of file NavierStokesSolverBlocks.cpp.

◆ applyBoundaryConditions() [2/2]

void applyBoundaryConditions	(	bcPtr_Type &	bc,
		const Real &	time,
		const vectorPtr_Type &	velocities
	)

Apply the BCs semi-implicit case (example aorta)

Parameters

bc	boundary conditions of the problem
time	current time
velocities	vector of velocities to be imposed

Definition at line 702 of file NavierStokesSolverBlocks.cpp.

◆ applyBoundaryConditionsJacobian()

void applyBoundaryConditionsJacobian ( bcPtr_Type & bc )

Apply the BCs on the Jacobian matrix.

Parameters

bc	boundary conditions of the problem

Definition at line 1584 of file NavierStokesSolverBlocks.cpp.

◆ applyBoundaryConditionsSolution()

void applyBoundaryConditionsSolution ( const Real & time )

Apply the BCs to the solution.

Parameters

time	current time

Definition at line 915 of file NavierStokesSolverBlocks.cpp.

◆ updateConvectiveTerm()

void updateConvectiveTerm ( const vectorPtr_Type & velocity )

Update the convective term.

Parameters

velocity velocity vector

Definition at line 1357 of file NavierStokesSolverBlocks.cpp.

◆ computeForces()

VectorSmall< 2 > computeForces	(	BCHandler &	bcHandlerDrag,
		BCHandler &	bcHandlerLift
	)

Compute Aerodynamic Forces.

Parameters

bcHandlerDrag	bc to be used for drag computation
bcHandlerLift	bc to be used for lift computation

Definition at line 845 of file NavierStokesSolverBlocks.cpp.

◆ solveTimeStep()

void solveTimeStep ( )

Solve a time step.

Definition at line 740 of file NavierStokesSolverBlocks.cpp.

◆ computeForcesNonLinear()

void computeForcesNonLinear	(	vectorPtr_Type &	force,
		const vectorPtr_Type &	solution
	)

Compute Aerodynamic Forces - nonlinear case.

Parameters

force	vector containing the forces
solution	vector with the solution

Definition at line 1100 of file NavierStokesSolverBlocks.cpp.

◆ evaluateResidual() [1/2]

void evaluateResidual	(	const vectorPtr_Type &	convective_velocity,
		const vectorPtr_Type &	velocity_km1,
		const vectorPtr_Type &	pressure_km1,
		const vectorPtr_Type &	rhs_velocity,
		vectorPtr_Type &	residual
	)

Evaluates the fluid residual in FSI simulations.

Parameters

convective_velocity	difference between fluid velocity and mesh velocity at previous Newton iteration
velocity_km1	fluid velocity at previous Newton iteration
pressure_km1	fluid pressure at previous Newton iteration
pressure_km1	terms of bdf associated to the rhs
residual	residual that will be assembled

Definition at line 999 of file NavierStokesSolverBlocks.cpp.

◆ evaluateResidual() [2/2]

void evaluateResidual	(	const vectorPtr_Type &	convective_velocity,
		const vectorPtr_Type &	velocity_km1,
		const vectorPtr_Type &	pressure_km1,
		const vectorPtr_Type &	rhs_velocity,
		vectorPtr_Type &	residualVelocity,
		vectorPtr_Type &	residualPressure
	)

Evaluates the fluid residual in FSI simulations.

Parameters

convective_velocity	difference between fluid velocity and mesh velocity at previous Newton iteration
velocity_km1	fluid velocity at previous Newton iteration
pressure_km1	fluid pressure at previous Newton iteration
pressure_km1	terms of bdf associated to the rhs
residualVelocity	residual associated to the fluid momentum equation
residualPressure	residual associated to the fluid continuity equation

Definition at line 926 of file NavierStokesSolverBlocks.cpp.

◆ updateStabilization()

void updateStabilization	(	const vector_Type &	convective_velocity_previous_newton_step,
		const vector_Type &	velocity_previous_newton_step,
		const vector_Type &	pressure_previous_newton_step,
		const vector_Type &	velocity_rhs
	)

Update stabilization terms.

Parameters

convective_velocity_previous_newton_step	difference between fluid velocity and mesh velocity at previous Newton iteration
velocity_previous_newton_step	fluid velocity at previous Newton iteration
pressure_previous_newton_step	fluid pressure at previous Newton iteration
velocity_rhs	block of the velocity of the rhs

Definition at line 1428 of file NavierStokesSolverBlocks.cpp.

◆ integrateForces()

void integrateForces	(	const vectorPtr_Type &	velocity,
		const vectorPtr_Type &	pressure
	)

Computation of forces.

Parameters

velocity	velocity vector
pressure	pressure vector

◆ preprocessBoundary()

void preprocessBoundary	(	const Real &	nx,
		const Real &	ny,
		const Real &	nz,
		BCHandler &	bc,
		Real &	Q_hat,
		const vectorPtr_Type &	Phi_h,
		const UInt	flag,
		vectorPtr_Type &	Phi_h_flag,
		vectorPtr_Type &	V_hat_x,
		vectorPtr_Type &	V_hat_y,
		vectorPtr_Type &	V_hat_z
	)

Additional method used for pre-processing on non-circular boundaries (used only in example aorta)

Parameters

nx	x component normal vector to an outflow face
ny	y component normal vector to an outflow face
nz	z component normal vector to an outflow face
bc	boundary conditions
Q_hat	reference flowrate to be imposed
Phi_h	solution laplacian
flag	flag outflow
Phi_h_flag	solution laplacian on the outflow flag
V_hat_x	reference velocity, x component
V_hat_y	reference velocity, y component
V_hat_z	reference velocity, z component

Definition at line 1591 of file NavierStokesSolverBlocks.cpp.

◆ solveLaplacian()

void solveLaplacian	(	const UInt &	flag,
		bcPtr_Type &	bc_laplacian,
		vectorPtr_Type &	laplacianSolution
	)

Additional method used for pre-processing on non-circular boundaries (used only in example aorta)

Parameters

flag	flag of the outflow face
bc_laplacian	boundary conditions laplacian
laplacianSolution	solution laplacian

Definition at line 1640 of file NavierStokesSolverBlocks.cpp.

◆ setupPostProc()

void setupPostProc ( )

Additional setup for postprocessing on boundaries.

Definition at line 1712 of file NavierStokesSolverBlocks.cpp.

Here is the caller graph for this function:

◆ flux()

Real flux	(	const markerID_Type &	flag,
		const vector_Type &	velocity
	)

Compute flux through a boundary face of the domain.

Parameters

flag	flag of the face
velocity	velocity vector

Definition at line 1724 of file NavierStokesSolverBlocks.cpp.

◆ area()

Real area ( const markerID_Type & flag )

Compute area of a boundary face of the domain.

Parameters

flag	flag of the outflow face

Definition at line 1731 of file NavierStokesSolverBlocks.cpp.

◆ pres()

Real pres	(	const markerID_Type &	flag,
		const vector_Type &	pressure
	)

Compute mean pressure at a boundary face of the domain.

Parameters

flag	flag of the outflow face
pressure	pressure vector

Definition at line 1749 of file NavierStokesSolverBlocks.cpp.

◆ geometricCenter()

Vector geometricCenter ( const markerID_Type & flag )

Compute center of a boundary face of the domain.

Parameters

flag	flag of the face

Definition at line 1737 of file NavierStokesSolverBlocks.cpp.

◆ normal()

Vector normal ( const markerID_Type & flag )

Compute normal of a boundary face of the domain.

Parameters

flag	flag of the face

Definition at line 1743 of file NavierStokesSolverBlocks.cpp.

◆ setAlpha()

void setAlpha ( const Real & alpha )

inline

Set coefficient associated to the time discretization scheme.

Parameters

alpha coefficient BDF scheme in front of u^{n+1}

Definition at line 400 of file NavierStokesSolverBlocks.hpp.

Here is the caller graph for this function:

◆ setRhsVelocity()

void setRhsVelocity ( const vectorPtr_Type & vel_rhs )

inline

Set the rhs vector, velocity component.

Parameters

vel_rhs rhs vector, velocity component

Definition at line 413 of file NavierStokesSolverBlocks.hpp.

Here is the caller graph for this function:

◆ setExtrapolatedPressure()

void setExtrapolatedPressure ( const vectorPtr_Type & pressure_extrapolated )

inline

Set the extrapolated pressure vector (used by semi-implicit VMS-LES stabilization)

Parameters

pressure_extrapolated pressure extrapolated vector

Definition at line 422 of file NavierStokesSolverBlocks.hpp.

Here is the caller graph for this function:

◆ setExportFineScaleVelocity()

void setExportFineScaleVelocity	(	ExporterHDF5< mesh_Type > &	exporter,
		const int &	numElementsTotal
	)

Set if one needs to export the fine scale component (when using LES models)

Parameters

exporter	exporter
numElementsTotal	number of total elements (tetrahedral elements)

Definition at line 206 of file NavierStokesSolverBlocks.cpp.

◆ setParameters()

void setParameters ( )

Set parameters of the solver.

Definition at line 68 of file NavierStokesSolverBlocks.cpp.

Here is the caller graph for this function:

◆ setTimeStep()

void setTimeStep ( const Real & dt )

inline

Set time step.

Parameters

dt	time step size

Definition at line 438 of file NavierStokesSolverBlocks.hpp.

Here is the caller graph for this function:

◆ setBoundaryConditions()

void setBoundaryConditions ( const bcPtr_Type & bc )

Set time step.

Parameters

bc	boundary conditions of the problem

Definition at line 1756 of file NavierStokesSolverBlocks.cpp.

◆ uFESpace()

const std::shared_ptr<FESpace<mesh_Type, map_Type> >& uFESpace ( ) const

inline

Get the velocity FE space.

Definition at line 455 of file NavierStokesSolverBlocks.hpp.

◆ pFESpace()

const std::shared_ptr<FESpace<mesh_Type, map_Type> >& pFESpace ( ) const

inline

Get the velocity FE space.

Definition at line 461 of file NavierStokesSolverBlocks.hpp.

◆ uFESpace_scalar()

const std::shared_ptr<FESpace<mesh_Type, map_Type> >& uFESpace_scalar ( ) const

inline

Get the velocity FE space.

Definition at line 467 of file NavierStokesSolverBlocks.hpp.

◆ updateVelocity()

void updateVelocity ( vectorPtr_Type & velocity )

inline

Get the velocity vector.

Definition at line 473 of file NavierStokesSolverBlocks.hpp.

Here is the caller graph for this function:

◆ updatePressure()

void updatePressure ( vectorPtr_Type & pressure )

inline

Get the pressure vector.

Definition at line 479 of file NavierStokesSolverBlocks.hpp.

Here is the caller graph for this function:

◆ getF()

matrixPtr_Type const& getF ( ) const

inline

Get the F block.

Definition at line 485 of file NavierStokesSolverBlocks.hpp.

◆ getJacobian()

matrixPtr_Type const& getJacobian ( ) const

inline

Get the jacobian block.

Definition at line 491 of file NavierStokesSolverBlocks.hpp.

◆ getBtranspose()

matrixPtr_Type const& getBtranspose ( ) const

inline

Get the Btranspose block.

Definition at line 497 of file NavierStokesSolverBlocks.hpp.

◆ getB()

matrixPtr_Type const& getB ( ) const

inline

Get the B block.

Definition at line 503 of file NavierStokesSolverBlocks.hpp.

◆ getRhs()

vectorPtr_Type const& getRhs ( ) const

inline

Get the rhs.

Definition at line 509 of file NavierStokesSolverBlocks.hpp.

◆ getRhsPressure()

vectorPtr_Type const& getRhsPressure ( ) const

inline

Get the rhs vector (pressure part)

Definition at line 515 of file NavierStokesSolverBlocks.hpp.

◆ Mu()

matrixPtr_Type const& Mu ( ) const

inline

Definition at line 520 of file NavierStokesSolverBlocks.hpp.

◆ density()

Real density ( ) const

inline

Definition at line 525 of file NavierStokesSolverBlocks.hpp.

◆ viscosity()

Real viscosity ( ) const

inline

Definition at line 530 of file NavierStokesSolverBlocks.hpp.

◆ block00()

matrixPtr_Type const& block00 ( ) const

inline

Get the (0,0) block.

Definition at line 536 of file NavierStokesSolverBlocks.hpp.

◆ block01()

matrixPtr_Type const& block01 ( ) const

inline

Get the (0,1) block.

Definition at line 542 of file NavierStokesSolverBlocks.hpp.

◆ block10()

matrixPtr_Type const& block10 ( ) const

inline

Get the (1,0) block.

Definition at line 548 of file NavierStokesSolverBlocks.hpp.

◆ block11()

matrixPtr_Type const& block11 ( ) const

inline

Get the (1,1)

Definition at line 554 of file NavierStokesSolverBlocks.hpp.

◆ useStabilization()

bool useStabilization ( ) const

inline

Get if using a stabilization.

Definition at line 560 of file NavierStokesSolverBlocks.hpp.

◆ assembleInterfaceMass()

void assembleInterfaceMass	(	matrixPtr_Type &	mass_interface,
		const mapPtr_Type &	interface_map,
		markerID_Type	interfaceFlag,
		const vectorPtr_Type &	numerationInterface,
		const UInt &	offset
	)