AddPressureDirichletBC

This routine sets the Dirichlet boundary condition for pressure field in AbstractSteadyStokes_ kernel.

• It also makes the obj%DBCForPressure(dbcNo)
• dbcNo should be lesser than total dirichlet boundary condition for pressure field

Interface

܀ Interface

INTERFACE
  MODULE SUBROUTINE addPressureDirichletBC(obj, dbcNo, param, boundary)
    CLASS(AbstractSteadyStokes_), INTENT(INOUT) :: obj
    INTEGER(I4B), INTENT(IN) :: dbcNo
    !! Dirichlet boundary nunber
    TYPE(ParameterList_), INTENT(IN) :: param
    !! parameter for constructing [DirichletBC_](DirichletBC_.md).
    TYPE(MeshSelection_), INTENT(IN) :: boundary
    !! Boundary region
  END SUBROUTINE addPressureDirichletBC
END INTERFACE

️܀ See example

In this example we will learn how to add pressure and velocity dirichlet boundary condition

PROGRAM main
  USE easifemBase
  USE easifemClasses
  USE easifemMaterials
  USE easifemKernels
  USE SteadyStokes221_Class
  IMPLICIT NONE

Declare variables

TYPE( SteadyStokes221_ ) :: obj
TYPE( ParameterList_ ) :: param
TYPE( DomainPointer_ ) :: domains(2)
CLASS( Domain_ ), POINTER :: dom => NULL()
CHARACTER( LEN = * ), PARAMETER :: domainFileNamePressure="./mesh_tri3.h5"
CHARACTER( LEN = * ), PARAMETER :: domainFileNameVelocity="./mesh_tri6.h5"
TYPE(String) :: filename(2)
TYPE( MeshSelection_ ) :: region

New variables in this example

CLASS( DirichletBC_ ), POINTER :: dbc => NULL()

Set parameters

CALL FPL_INIT(); CALL param%Initiate()

CALL SetSteadyStokes221Param( &
  & param=param, &
  & isConservativeForm=.TRUE., &
  & gravity = [0.0_DFP, -9.8_DFP, 0.0_DFP], &
  & domainFileForPressure = domainFileNamePressure, &
  & domainFileForVelocity = domainFileNameVelocity, &
  & engine="NATIVE_SERIAL", &
  & CoordinateSystem=KERNEL_CARTESIAN, &
  & maxIter = 100, &
  & rtoleranceForPressure = REAL( 1.0E-6, DFP ), &
  & rtoleranceForVelocity = REAL( 1.0E-6, DFP ), &
  & atoleranceForPressure = REAL( 1.0E-6, DFP ), &
  & atoleranceForVelocity = REAL( 1.0E-6, DFP ), &
  & toleranceForSteadyState = REAL( 1.0E-6, DFP ), &
  & tFluidMaterials=1, &
  & tDirichletBCForPressure=1, &
  & tDirichletBCForVelocity=3, &
  & baseInterpolationForSpace="LagrangeInterpolation", &
  & baseContinuityForSpace="H1", &
  & quadratureTypeForSpace="GaussLegendre", &
  & postProcessOpt=1)

Set parameters for linear solver.

CALL SetLinSolverParam( &
  & param=param, &
  & solverName=LIS_GMRES,&
  & preconditionOption=NO_PRECONDITION, &
  & convergenceIn=convergenceInRes, &
  & convergenceType=relativeConvergence, &
  & maxIter=100, &
  & relativeToRHS=.TRUE., &
  & KrylovSubspaceSize=20, &
  & rtol=REAL( 1.0E-10, DFP ), &
  & atol=REAL( 1.0D-10, DFP ) )

Initiate the domain.

filename = [String(domainFileNamePressure), &
  & String(domainFileNameVelocity)]
CALL e%setQuietMode(.TRUE.)
CALL Initiate(domains=domains, filename=filename)

Initiate kernel

CALL e%setQuietMode(.FALSE.)
CALL obj%Initiate(param=param, domains=domains )

Add fluid material to kernel. To do so, we first create an instance of MeshSelection_. Then we add this instance to the kernel.

CALL e%setQuietMode(.TRUE.)
CALL region%Initiate( isSelectionByMeshID=.TRUE. )
CALL region%Add( dim=obj%nsd, meshID=[1] )
CALL SetFluidMaterialParam( &
  & param=param, &
  & name="fluidMaterial", &
  & massDensity=1000.0_DFP, &
  & dynamicViscosity = 0.001_DFP, &
  & stressStrainModel="NewtonianFluidModel" )
CALL SetNewtonianFluidModelParam( &
  & param = param, &
  & dynamicViscosity = 0.001_DFP )

CALL e%setQuietMode(.FALSE.)
CALL obj%AddFluidMaterial( &
  & materialNo=1, &
  & materialName="fluidMaterial", &
  & param=param, &
  & region=region)
CALL region%Deallocate()

Now we show how to add dirichlet boundary condition. To this end first we create an instance of [[MeshSelection_]] to select the region of the mesh. Then we define the dirichlet bonundary condition, and pass these two information to kernel.

Set parameters for dirichlet boundary condition:

CALL SetDirichletBCParam( &
  & param=param, &
  & name="ZeroV1", &
  & idof=1, &
  & nodalValueType=Constant, &
  & useFunction=.FALSE. )

select the mesh region:

#define BOTTOM 1
#define RIGHT 2
#define TOP 3
#define LEFT 4

CALL region%Initiate( isSelectionByMeshID=.TRUE. )
CALL region%Add( dim=obj%nsd-1, meshID=[BOTTOM, RIGHT, LEFT] )
CALL region%Set()

Add dirichlet boundary condition and the region to kernel:

CALL obj%AddVelocityDirichletBC( &
  & dbcNo=1, &
  & param=param, &
  & boundary=region )
dbc => obj%GetVelocityDirichletBCPointer( dbcNo=1 )
CALL dbc%Set( ConstantNodalValue=0.0_DFP )
dbc=>NULL()

AddDirichletBC, V1=U, set parameters for dirichlet boundary condition:

CALL SetDirichletBCParam( &
  & param=param, &
  & name="UpstreamV1", &
  & idof=1, &
  & nodalValueType=Constant, &
  & useFunction=.FALSE. )

select the mesh region:

CALL region%Deallocate()
CALL region%Initiate( isSelectionByMeshID=.TRUE. )
CALL region%Add( dim=obj%nsd-1, meshID=[TOP] )
CALL region%Set()

Add dirichlet boundary condition and the region to kernel:

CALL obj%AddVelocityDirichletBC( &
  & dbcNo=2, &
  & param=param, &
  & boundary=region )
dbc => obj%GetVelocityDirichletBCPointer( dbcNo=2 )
CALL dbc%Set( ConstantNodalValue=0.1_DFP )
dbc=>NULL()
CALL region%Deallocate()

AddDirichletBC, V2=0, Set parameters for dirichlet boundary condition:

CALL SetDirichletBCParam( &
  & param=param, &
  & name="ZeroV2", &
  & idof=2, &
  & nodalValueType=Constant, &
  & useFunction=.FALSE. )

select the mesh region:

CALL region%Deallocate()
CALL region%Initiate( isSelectionByMeshID=.TRUE. )
CALL region%Add( dim=obj%nsd-1, meshID=[BOTTOM, RIGHT, TOP, LEFT] )
CALL region%Set()

Add dirichlet boundary condition and the region to kernel:

CALL obj%AddVelocityDirichletBC( &
  & dbcNo=3, &
  & param=param, &
  & boundary=region )
dbc => obj%GetVelocityDirichletBCPointer( dbcNo=3 )
CALL dbc%Set( ConstantNodalValue=0.0_DFP ); dbc=>NULL()
CALL region%Deallocate()

AddDirichletBC, P=0, Set parameters for dirichlet boundary condition:

CALL SetDirichletBCParam( &
  & param=param, &
  & name="ZeroP", &
  & idof=1, &
  & nodalValueType=Constant, &
  & useFunction=.FALSE. )

select the mesh region:

CALL region%Deallocate()
CALL region%Initiate( isSelectionByMeshID=.TRUE. )
CALL region%Add( dim=obj%nsd-1, meshID=[BOTTOM] )
CALL region%Set()

Add dirichlet boundary condition and the region to kernel:

CALL obj%AddPressureDirichletBC( dbcNo=1, param=param, &
  & boundary=region )
dbc => obj%GetPressureDirichletBCPointer( dbcNo=1 )
CALL dbc%Set( ConstantNodalValue=0.0_DFP ); dbc=>NULL()
CALL region%Deallocate()

CALL obj%Display("")

CALL obj%Deallocate()
CALL param%Deallocate(); CALL FPL_FINALIZE()
END PROGRAM main

↢