Skip to content

Updates for TDGLinearElasticity Kernel #534

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
shishiousan wants to merge 7 commits into
base: vikas-dev
Choose a base branch
from
Open

Updates for TDGLinearElasticity Kernel #534

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
17aa70e
Updates in STVectorFieldClass
shishiousan Dec 13, 2025
eba71c5
Updates in STVectorFieldClass
shishiousan Dec 13, 2025
3c4c2d1
Updates in VectorField Class
shishiousan Dec 13, 2025
bd6da1f
Updates in STVectorField Class
shishiousan Dec 14, 2025
05a15dd
Adding TDGAlgorithm2 Class
shishiousan Dec 14, 2025
4e1a327
Updates in STVectorFields Classes
shishiousan Dec 14, 2025
1e5a2f1
Updates in STVectorFields
shishiousan Dec 14, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
267 changes: 258 additions & 9 deletions src/modules/STVectorField/src/STVectorField_Class.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -19,16 +19,17 @@
! summary: STVector field data type is defined

MODULE STVectorField_Class
USE GlobalData, ONLY: DFP, I4B, LGT, DOF_FMT, NODES_FMT, NodesToDOF
USE BaseType, ONLY: FEVariable_
USE AbstractField_Class, ONLY: AbstractField_
USE AbstractNodeField_Class, ONLY: AbstractNodeField_
USE BaseType, ONLY: FEVariable_
USE DirichletBC_Class, ONLY: DirichletBC_, DirichletBCPointer_
USE ExceptionHandler_Class, ONLY: e
USE HDF5File_Class, ONLY: HDF5File_
USE FEDOF_Class, ONLY: FEDOF_, FEDOFPointer_
USE DirichletBC_Class, ONLY: DirichletBC_, DirichletBCPointer_
USE UserFunction_Class, ONLY: UserFunction_
USE GlobalData, ONLY: DFP, I4B, LGT, DOF_FMT, NODES_FMT, NodesToDOF
USE HDF5File_Class, ONLY: HDF5File_
USE TimeFEDOF_Class, ONLY: TimeFEDOF_, TimeFEDOFPointer_
USE UserFunction_Class, ONLY: UserFunction_
USE VectorField_Class, ONLY: VectorField_

IMPLICIT NONE

Expand Down Expand Up @@ -68,7 +69,10 @@ MODULE STVectorField_Class

! CONSTRUCTOR:
! @ConstructorMethods
PROCEDURE, NON_OVERRIDABLE, PUBLIC, PASS(obj) :: Initiate2 => obj_Initiate2
PROCEDURE, PUBLIC, PASS(obj) :: Initiate2 => obj_Initiate2
!! Initiate by copy

PROCEDURE, PUBLIC, PASS(obj) :: Initiate4 => obj_Initiate4
!! Initiate by copy

PROCEDURE, PUBLIC, PASS(obj) :: DEALLOCATE => obj_Deallocate
Expand Down Expand Up @@ -146,10 +150,11 @@ MODULE STVectorField_Class
Set7, Set8, Set9, Set10, Set11, Set12, Set13, Set14, Set15, Set16

PROCEDURE, PUBLIC, NON_OVERRIDABLE, PASS(obj) :: SetByFunction => obj_SetByFunction

!! Set by function
PROCEDURE, PUBLIC, NON_OVERRIDABLE, PASS(obj) :: SetFromVectorField => &
obj_SetFromVectorField
!! Set by function
PROCEDURE, PUBLIC, NON_OVERRIDABLE, PASS(obj) :: SetToVectorField => &
obj_SetToVectorField

! GET:
! @GetMethods
Expand Down Expand Up @@ -204,6 +209,36 @@ MODULE STVectorField_Class
GENERIC, PUBLIC :: ApplyDirichletBC => ApplyDirichletBC1, &
ApplyDirichletBC2, ApplyDirichletBC3

! SET:
! @BodySourceMethods
PROCEDURE, NON_OVERRIDABLE, PASS(obj) :: &
ApplyBodySource1 => obj_ApplyBodySource1
!! Add contribution of body source to the scalar field
!! body source is given as user function
PROCEDURE, NON_OVERRIDABLE, PASS(obj) :: &
ApplyBodySource2 => obj_ApplyBodySource2
!! Add contribution of body source to the scalar field
!! body source is given external scalar field
PROCEDURE, NON_OVERRIDABLE, PASS(obj) :: &
ApplyBodySource3 => obj_ApplyBodySource3
!! Add contribution of body source to the scalar field
!! body source is given as user function
GENERIC, PUBLIC :: ApplyBodySource => ApplyBodySource1, &
ApplyBodySource2, ApplyBodySource3
!! Generic method for setting body source

! SET:
! @SurfaceNBCMethods
PROCEDURE, PUBLIC, NON_OVERRIDABLE, PASS(obj) :: ApplySurfaceNeumannBC => &
obj_ApplySurfaceNeumannBC
!! Apply Surface neumann boundary condition

! SET:
! @PointNBCMethods
PROCEDURE, PUBLIC, NON_OVERRIDABLE, PASS(obj) :: ApplyPointNeumannBC => &
obj_ApplyPointNeumannBC
!! Apply point neumann boundary condition

END TYPE STVectorField_

!---------------------------------------------------------------------------
Expand All @@ -222,7 +257,7 @@ MODULE STVectorField_Class
! date: 2023-03-29
! summary: Initiate2

INTERFACE STVectorFieldInitiate
INTERFACE
MODULE SUBROUTINE obj_Initiate2( &
obj, obj2, copyFull, copyStructure, usePointer)
CLASS(STVectorField_), INTENT(INOUT) :: obj
Expand All @@ -232,6 +267,93 @@ MODULE SUBROUTINE obj_Initiate2( &
LOGICAL(LGT), OPTIONAL, INTENT(IN) :: copyStructure
LOGICAL(LGT), OPTIONAL, INTENT(IN) :: usePointer
END SUBROUTINE obj_Initiate2
END INTERFACE

INTERFACE STVectorFieldInitiate
MODULE PROCEDURE obj_Initiate2
END INTERFACE STVectorFieldInitiate

!----------------------------------------------------------------------------
!
!----------------------------------------------------------------------------

INTERFACE
MODULE SUBROUTINE obj_Initiate4( &
obj, name, engine, fieldType, storageFMT, comm, local_n, global_n, &
spaceCompo, isSpaceCompo, isSpaceCompoScalar, timeCompo, isTimeCompo, &
isTimeCompoScalar, tPhysicalVarNames, physicalVarNames, &
isPhysicalVarNames, isPhysicalVarNamesScalar, tNodes, isTNodes, &
isTNodesScalar, tSize, fedof, geofedof, timefedof)
CLASS(STVectorField_), INTENT(INOUT) :: obj
CHARACTER(*), INTENT(IN) :: name
!! name of the field, needed
CHARACTER(*), INTENT(IN) :: engine
!! name of the engine, needed
INTEGER(I4B), OPTIONAL, INTENT(IN) :: fieldType
!! field type, default is FIELD_TYPE_NORMAL, needed
!! following options are available FIELD_TYPE_NORMAL
!! FIELD_TYPE_CONSTANT
INTEGER(I4B), OPTIONAL, INTENT(IN) :: storageFMT
!! storage format of the Vector field
!! Not required.
INTEGER(I4B), OPTIONAL, INTENT(IN) :: comm
!! communication group
!! Only needed for parallel environment
INTEGER(I4B), OPTIONAL, INTENT(IN) :: local_n
!! local size of field on each processor
!! Only needed for parallel environment
INTEGER(I4B), OPTIONAL, INTENT(IN) :: global_n
!! global size of field on distributed on processors
!! Only needed for parallel environment
INTEGER(I4B), OPTIONAL, INTENT(IN) :: spaceCompo(:)
!! space components
LOGICAL(LGT), OPTIONAL, INTENT(IN) :: isSpaceCompo
!! if true we will try to access spaceCompo, NOT REQUIRED
LOGICAL(LGT), OPTIONAL, INTENT(IN) :: isSpaceCompoScalar
!! is space component scalar,
!! in this case we only access spaceCompo(1), NOT REQUIRED
INTEGER(I4B), OPTIONAL, INTENT(IN) :: timeCompo(:)
!! Time components, needed for space-time field
LOGICAL(LGT), OPTIONAL, INTENT(IN) :: isTimeCompo
!! if true we will try to access TimeCompo
!! Not required, Not needed
LOGICAL(LGT), OPTIONAL, INTENT(IN) :: isTimeCompoScalar
!! is Time component scalar,
!! in this case we only access TimeCompo(1), Not needed
INTEGER(I4B), OPTIONAL, INTENT(IN) :: tPhysicalVarNames
!! total physical variable names
!! if it is zero, then physicalVarNames will not be written
!! evenif physicalVarNames is present, and isPhysicalVarNames
!! is true, Not required
CHARACTER(*), OPTIONAL, INTENT(IN) :: physicalVarNames(:)
!! Names of the physical variables, NOT REQUIRED
LOGICAL(LGT), OPTIONAL, INTENT(IN) :: isPhysicalVarNames
!! logical variable to check if physicalVarNames is present or not
!! if it is false then physicalVarNames will not be written
!! NOT REQUIRED
LOGICAL(LGT), OPTIONAL, INTENT(IN) :: isPhysicalVarNamesScalar
!! if true then physicalVarNames is scalar
INTEGER(I4B), OPTIONAL, INTENT(IN) :: tNodes(:)
!! total number of nodes in each physical variable
!! Not required
LOGICAL(LGT), OPTIONAL, INTENT(IN) :: isTNodes
!! if true we will try to access tNodes
!! Not required
LOGICAL(LGT), OPTIONAL, INTENT(IN) :: isTNodesScalar
!! is tNodes scalar
!! Not required
INTEGER(I4B), OPTIONAL, INTENT(IN) :: tSize
!! total size of node field
!! not required
CLASS(FEDOF_), TARGET, INTENT(IN) :: fedof, geofedof
!! FEDOF object
CLASS(TimeFEDOF_), OPTIONAL, TARGET, INTENT(IN) :: timefedof
!! TimeFEDOF object
END SUBROUTINE obj_Initiate4
END INTERFACE

INTERFACE STVectorFieldInitiate
MODULE PROCEDURE obj_Initiate4
END INTERFACE STVectorFieldInitiate

!----------------------------------------------------------------------------
Expand Down Expand Up @@ -879,6 +1001,25 @@ MODULE SUBROUTINE obj_SetFromVectorField(obj, VALUE, timeCompo, &
END SUBROUTINE obj_SetFromVectorField
END INTERFACE

!----------------------------------------------------------------------------
! Set@SetMethods
!----------------------------------------------------------------------------

!> author: Shion Shimizu
! date: 2025-12-11
! summary: Set to VectorField_

INTERFACE
MODULE SUBROUTINE obj_SetToVectorField(obj, VALUE, timeCompo, &
scale, addContribution)
CLASS(STVectorField_), INTENT(INOUT) :: obj
CLASS(VectorField_), INTENT(INOUT) :: VALUE
INTEGER(I4B), INTENT(IN) :: timeCompo
REAL(DFP), OPTIONAL, INTENT(IN) :: scale
LOGICAL(LGT), OPTIONAL, INTENT(IN) :: addContribution
END SUBROUTINE obj_SetToVectorField
END INTERFACE

!----------------------------------------------------------------------------
! Set@SetMethods
!----------------------------------------------------------------------------
Expand Down Expand Up @@ -1262,6 +1403,114 @@ MODULE SUBROUTINE obj_ApplyDirichletBC3(obj, times)
END SUBROUTINE obj_ApplyDirichletBC3
END INTERFACE

!----------------------------------------------------------------------------
! ApplyBodySource@NBCMethods
!----------------------------------------------------------------------------

!> author: Shion Shimizu
! date: 2025-12-11
! summary: Apply Body source to ST vector field

INTERFACE
MODULE SUBROUTINE obj_ApplyBodySource1(obj, bodySource, scale, times)
CLASS(STVectorField_), INTENT(INOUT) :: obj
CLASS(UserFunction_), INTENT(INOUT) :: bodySource
!! Body source user function
!! It should be a vector function with
!! total arguments 4 (x, y, z, time)
REAL(DFP), INTENT(IN) :: scale
!! scale for body source
!! obj = obj + scale * bodySource integral
REAL(DFP), INTENT(IN) :: times(:)
!! time, which will be passed to the body source function
END SUBROUTINE obj_ApplyBodySource1
END INTERFACE

!----------------------------------------------------------------------------
! ApplyBodySource@NBCMethods
!----------------------------------------------------------------------------

!> author: Shion Shimizu
! date: 2025-12-11
! summary: Apply Body source to ST vector field

INTERFACE
MODULE SUBROUTINE obj_ApplyBodySource2(obj, bodySource, scale)
CLASS(STVectorField_), INTENT(INOUT) :: obj
!! ST vector field
!! The test function will be corresponding to the obj
CLASS(STVectorField_), INTENT(INOUT) :: bodySource
!! Body source in terms of vector field
REAL(DFP), INTENT(IN) :: scale
!! scale for body source
!! obj = obj + scale * bodySource integral
END SUBROUTINE obj_ApplyBodySource2
END INTERFACE

!----------------------------------------------------------------------------
! ApplyBodySource@NBCMethods
!----------------------------------------------------------------------------

!> author: Shion Shimizu
! date: 2025-12-11
! summary: Apply Body source to ST vector field

INTERFACE
MODULE SUBROUTINE obj_ApplyBodySource3(obj, bodySource, scale, times)
CLASS(STVectorField_), INTENT(INOUT) :: obj
CLASS(UserFunction_), INTENT(INOUT) :: bodySource
!! Body source user function
!! It should be a vector function with
!! total arguments 4 (x, y, z, time)
REAL(DFP), INTENT(IN) :: scale
!! scale for body source
!! obj = obj + scale * bodySource integral
REAL(DFP), INTENT(IN) :: times
!! time, which will be passed to the body source function
!! This time can also represent a quadrature point time
END SUBROUTINE obj_ApplyBodySource3
END INTERFACE

!----------------------------------------------------------------------------
! ApplyNeumannBC@NBCMethods
!----------------------------------------------------------------------------

!> author: Shion Shimizu
! date: 2025-12-11
! summary: Add Contribution of neumann boundary condition

INTERFACE
MODULE SUBROUTINE obj_ApplySurfaceNeumannBC(obj, nbcField, scale, times)
CLASS(STVectorField_), INTENT(INOUT) :: obj
!! Vector field
CLASS(STVectorField_), INTENT(INOUT) :: nbcField
!! Vector field where we will keep the neumann boundary condition
!! extension to the entire domain
REAL(DFP), INTENT(IN) :: scale
!! Scale for neumann boundary condition
REAL(DFP), INTENT(IN) :: times(:)
!! times
END SUBROUTINE obj_ApplySurfaceNeumannBC
END INTERFACE

!----------------------------------------------------------------------------
! ApplyPointNeumannBC@NBCMethods
!----------------------------------------------------------------------------

!> author: Shion Shimizu
! date: 2025-12-11
! summary: Add Contribution of point neumann boundary condition

INTERFACE
MODULE SUBROUTINE obj_ApplyPointNeumannBC(obj, scale, times)
CLASS(STVectorField_), INTENT(INOUT) :: obj
!! Vector field
REAL(DFP), INTENT(IN) :: scale
!! scale for neumann boundary condition
REAL(DFP), INTENT(IN) :: times(:)
END SUBROUTINE obj_ApplyPointNeumannBC
END INTERFACE

!----------------------------------------------------------------------------
!
!----------------------------------------------------------------------------
Expand Down
7 changes: 6 additions & 1 deletion src/modules/TDGAlgorithms/CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,4 +16,9 @@
#

set(src_path "${CMAKE_CURRENT_LIST_DIR}/src/")
target_sources(${PROJECT_NAME} PRIVATE ${src_path}/TDGAlgorithm1_Class.F90)
target_sources(
${PROJECT_NAME}
PRIVATE
${src_path}/TDGAlgorithm1_Class.F90
${src_path}/TDGAlgorithm2_Class.F90
)
Loading