SteadyStokes111
Description of Kernelβ
- (1) Equal order interpolation for velocity and pressure field
- (1) Coupled linear system (dof=[u, p])
- (1) Constant material properties
Variational multiscale stabilization method is used, and following options are allowed
- SubscalePressure = True or False
- BoundarySubscale = True or False
- WeakDBC = True or False
- Multiphase = True or False
TODOβ
- Body-force
- gravity
- conservative and non-conservative form.
- Neumann boundary condition
- higher order terms
- Positive and negative sign in case of HOT
- boundary subscale
- pressure subscale
- no-slip boundary condition
- free-slip boundary condition
- inlet boundary condition
- outlet boundary-condition
- add an option to calculate stabilization parameter only at the center of the element
- handle quadrature mapping in general form
- Add reference pressure node and pressure value.
Ideas for optimizationβ
- Constant stabilization parameter
- Better precondition for iterative solver
Methodsβ
ποΈ SteadyStokes111
Description of Kernel
ποΈ AddFluidMaterial
This routine adds Fluid material to the AbstractSteadyStokes.
ποΈ AddPressureDirichletBC
This routine sets the Dirichlet boundary condition for pressure field in AbstractSteadyStokes kernel.
ποΈ AddPressureNeumannBC
This routine sets the Neumann boundary condition for pressure field.
ποΈ AddVelocityDirichletBC
This routine sets the Dirichlet boundary condition for Velocity field in AbstractSteadyStokes kernel.
ποΈ AddVelocityNeumannBC
This routine sets the Neumann boundary condition for Velocity field in AbstractSteadyStokes kernel.
ποΈ ApplyDirichletBC
This routine apply the Dirichlet boundary condition to pressure and velocity field.
ποΈ Assemble
This procedure pointer assembles the problem
ποΈ AssembleRHS
This procedure pointer assembles the right-hand-side vector.
ποΈ AssembleTanmat
This procedure pointer assembles the global tangent matrix.
ποΈ AssembleTanmat_InternalFacet
Assemble tangent matrix internal facet.
ποΈ AssembleTanMat_NeumanFacet
Assemble tangent matrix internal facet.
ποΈ CheckEssentialParam
Check essential parameter.
ποΈ ComputeStabParam
Compute stabilization parameter.
ποΈ Display
Display the content of SteadyStokes221.
ποΈ Export
This routine exports the kernel to HDF5file
ποΈ GetPressureDirichletBCPointer
This routine returns the pointer to Dirichlet boundary condition of pressure field in AbstractSteadyStokes kernel, that is obj%DBCForPressure(dbcNo)%ptr.
ποΈ GetPressureNeumannBCPointer
This routine returns the pointer to Neumann boundary condition of pressure field in AbstractSteadyStokes kernel, that is obj%NBCForPressure(nbcNo)%ptr.
ποΈ GetVelocityDirichletBCPointer
This routine returns the pointer to Dirichlet boundary condition of Velocity field in AbstractSteadyStokes kernel, that is obj%DBCForVelocity(dbcNo)%ptr.
ποΈ GetVelocityNeumannBCPointer
This routine returns the pointer to Neumann boundary condition of Velocity field in AbstractSteadyStokes kernel, that is obj%NBCForVelocity(nbcNo)%ptr.
ποΈ Import
This routine imports the properties of kernel and initiates it.
ποΈ Initiate
This routine initiates the kernel.
ποΈ InitiateFields
This routine initiates the field variable, deferred type.
ποΈ IsConverged
Returns true if steady state is achieved.
ποΈ IsSteadyState
Returns true if steady state is achieved.
ποΈ Run
Run the simulation
ποΈ Set
Set the kernel.
ποΈ SetSteadyStokes111Param
Set SteadyStokes111 parameter.
ποΈ Solve
This subroutine solves the system of linear equation.
ποΈ SteadyStokesUtility
SetMeshData
ποΈ Update
This subroutine Updates the system of linear equation
ποΈ UpdateIteration
This subroutine update the state of the kernel
ποΈ WriteData
This routine writes the data in the hdf5 file.