OFELI's Logo

An Object Oriented Finite Element Library

Builds finite element arrays for transient incompressible fluid flow using Navier-Stokes equations in 2-D domains. Numerical approximation uses stabilized 3-node triangle finite elements for velocity and pressure. 2nd-order projection scheme is used for time integration. More...

#include <TINS2DT3S.h>

+ Inheritance diagram for TINS2DT3S:

Public Member Functions

 TINS2DT3S ()
 Default Constructor.
 
 TINS2DT3S (Mesh &mesh)
 Constructor using mesh. More...
 
 TINS2DT3S (Mesh &mesh, Vect< real_t > &u)
 Constructor using mesh and velocity. More...
 
 ~TINS2DT3S ()
 Destructor.
 
void setInput (EqDataType opt, Vect< real_t > &u)
 Set equation input data. More...
 
int runOneTimeStep ()
 Run one time step.
 
int run ()
 Run (in the case of one step run)
 
- Public Member Functions inherited from Equa_Fluid< 3, 6, 2, 4 >
 Equa_Fluid ()
 Default constructor. More...
 
virtual ~Equa_Fluid ()
 Destructor.
 
void Reynolds (const real_t &Re)
 Set Reynolds number.
 
void Viscosity (const real_t &visc)
 Set (constant) Viscosity.
 
void Viscosity (const string &exp)
 Set viscosity given by an algebraic expression.
 
void Density (const real_t &dens)
 Set (constant) Viscosity.
 
void Density (const string &exp)
 Set Density given by an algebraic expression.
 
void ThermalExpansion (const real_t *e)
 Set (constant) thermal expansion coefficient.
 
void ThermalExpansion (const string &exp)
 Set thermal expansion coefficient given by an algebraic expression.
 
void setMaterial ()
 Set material properties.
 
- Public Member Functions inherited from Equation< NEN_, NEE_, NSN_, NSE_ >
 Equation ()
 
 Equation (Mesh &mesh)
 Constructor with mesh instance. More...
 
 Equation (Mesh &mesh, Vect< real_t > &u)
 Constructor with mesh instance and solution vector. More...
 
 Equation (Mesh &mesh, Vect< real_t > &u, real_t &init_time, real_t &final_time, real_t &time_step)
 Constructor with mesh instance, matrix and right-hand side. More...
 
 ~Equation ()
 Destructor.
 
void updateBC (const Element &el, const Vect< real_t > &bc)
 Update Right-Hand side by taking into account essential boundary conditions. More...
 
void DiagBC (DOFSupport dof_type=NODE_DOF, int dof=0)
 Update element matrix to impose bc by diagonalization technique. More...
 
void LocalNodeVector (Vect< real_t > &b)
 Localize Element Vector from a Vect instance. More...
 
void ElementNodeVector (const Vect< real_t > &b, LocalVect< real_t, NEE_ > &be)
 Localize Element Vector from a Vect instance. More...
 
void SideNodeVector (const Vect< real_t > &b, LocalVect< real_t, NSE_ > &bs)
 Localize Side Vector from a Vect instance. More...
 
void SideSideVector (const Vect< real_t > &b, std::valarray< real_t > &bs)
 Localize Side Vector from a Vect instance. More...
 
void ElementNodeVectorSingleDOF (const Vect< real_t > &b, LocalVect< real_t, NEN_ > &be)
 Localize Element Vector from a Vect instance. More...
 
void ElementNodeVector (const Vect< real_t > &b, LocalVect< real_t, NEN_ > &be, int dof)
 Localize Element Vector from a Vect instance. More...
 
void ElementSideVector (const Vect< real_t > &b, LocalVect< real_t, NSE_ > &be)
 Localize Element Vector from a Vect instance. More...
 
void ElementVector (const Vect< real_t > &b, DOFSupport dof_type=NODE_DOF, int flag=0)
 Localize Element Vector. More...
 
void SideVector (const Vect< real_t > &b, std::valarray< real_t > &sb)
 Localize Side Vector. More...
 
void ElementNodeCoordinates ()
 Localize coordinates of element nodes. More...
 
void SideNodeCoordinates ()
 Localize coordinates of side nodes. More...
 
void ElementAssembly (Matrix< real_t > *A)
 Assemble element matrix into global one. More...
 
void ElementAssembly (BMatrix< real_t > &A)
 Assemble element matrix into global one. More...
 
void ElementAssembly (SkSMatrix< real_t > &A)
 Assemble element matrix into global one. More...
 
void ElementAssembly (SkMatrix< real_t > &A)
 Assemble element matrix into global one. More...
 
void ElementAssembly (SpMatrix< real_t > &A)
 Assemble element matrix into global one. More...
 
void ElementAssembly (TrMatrix< real_t > &A)
 Assemble element matrix into global one. More...
 
void DGElementAssembly (Matrix< real_t > *A)
 Assemble element matrix into global one for the Discontinuous Galerkin approximation. More...
 
void DGElementAssembly (SkSMatrix< real_t > &A)
 Assemble element matrix into global one for the Discontinuous Galerkin approximation. More...
 
void DGElementAssembly (SkMatrix< real_t > &A)
 Assemble element matrix into global one for the Discontinuous Galerkin approximation. More...
 
void DGElementAssembly (SpMatrix< real_t > &A)
 Assemble element matrix into global one for the Discontinuous Galerkin approximation. More...
 
void DGElementAssembly (TrMatrix< real_t > &A)
 Assemble element matrix into global one for the Discontinuous Galerkin approximation. More...
 
void SideAssembly (Matrix< real_t > *A)
 Assemble side (edge or face) matrix into global one. More...
 
void SideAssembly (SkSMatrix< real_t > &A)
 Assemble side (edge or face) matrix into global one. More...
 
void SideAssembly (SkMatrix< real_t > &A)
 Assemble side (edge or face) matrix into global one. More...
 
void SideAssembly (SpMatrix< real_t > &A)
 Assemble side (edge or face) matrix into global one. More...
 
void ElementAssembly (Vect< real_t > &v)
 Assemble element vector into global one. More...
 
void SideAssembly (Vect< real_t > &v)
 Assemble side (edge or face) vector into global one. More...
 
void AxbAssembly (const Element &el, const Vect< real_t > &x, Vect< real_t > &b)
 Assemble product of element matrix by element vector into global vector. More...
 
void AxbAssembly (const Side &sd, const Vect< real_t > &x, Vect< real_t > &b)
 Assemble product of side matrix by side vector into global vector. More...
 
size_t getNbNodes () const
 Return number of element nodes.
 
size_t getNbEq () const
 Return number of element equations.
 
real_t setMaterialProperty (const string &exp, const string &prop)
 Define a material property by an algebraic expression. More...
 
- Public Member Functions inherited from Equa
 Equa ()
 Default constructor.
 
virtual ~Equa ()
 Destructor.
 
void setMesh (Mesh &m)
 Define mesh and renumber DOFs after removing imposed ones.
 
MeshgetMesh () const
 Return reference to Mesh instance. More...
 
LinearSolver< real_t > & getLinearSolver ()
 Return reference to linear solver instance.
 
Matrix< real_t > * getMatrix () const
 Return pointer to matrix.
 
void setSolver (Iteration ls, Preconditioner pc=IDENT_PREC)
 Choose solver for the linear system. More...
 
void setLinearSolver (Iteration ls, Preconditioner pc=IDENT_PREC)
 Choose solver for the linear system. More...
 
void setMatrixType (int t)
 Choose type of matrix. More...
 
int solveLinearSystem (Matrix< real_t > *A, Vect< real_t > &b, Vect< real_t > &x)
 Solve the linear system with given matrix and right-hand side. More...
 
int solveLinearSystem (Vect< real_t > &b, Vect< real_t > &x)
 Solve the linear system with given right-hand side. More...
 

Detailed Description

Builds finite element arrays for transient incompressible fluid flow using Navier-Stokes equations in 2-D domains. Numerical approximation uses stabilized 3-node triangle finite elements for velocity and pressure. 2nd-order projection scheme is used for time integration.

Author
Rachid Touzani

Constructor & Destructor Documentation

◆ TINS2DT3S() [1/2]

TINS2DT3S ( Mesh mesh)

Constructor using mesh.

Parameters
[in]meshMesh instance

◆ TINS2DT3S() [2/2]

TINS2DT3S ( Mesh mesh,
Vect< real_t > &  u 
)

Constructor using mesh and velocity.

Parameters
[in]meshMesh instance
[in,out]uVect instance containing initial velocity. This vector is updated during computations and will therefore contain velocity at each time step

Member Function Documentation

◆ setInput()

void setInput ( EqDataType  opt,
Vect< real_t > &  u 
)

Set equation input data.

Parameters
[in]optParameter to select type of input (enumerated values)
  • INITIAL_FIELD: Initial temperature
  • BOUNDARY_CONDITION: Boundary condition (Dirichlet)
  • SOURCE: Body force applied to fluid
  • TRACTION: Heat flux (Neumann boundary condition)
  • VELOCITY_FIELD: Velocity vector (for the convection term)
[in]uVector containing input data (Vect instance)





Copyright © 1998-2022 Rachid Touzani