D:/demona/2.0/manual/doc/2-31set/global.f90

Go to the documentation of this file.

module global

use demonaParameters

implicit none

type setI
        
        integer, allocatable :: tVec(:) !< vector that stores the indices of local dof at that "i" index
        
        integer :: maxT !< size of the unkVec

        integer, allocatable :: endIndex(:)

end type setI


!< type 1
type setJ

        integer, allocatable :: iVec(:) !< vector that stores the indices of "i"s of that "j" index

        integer :: maxI !< size of the iVec
        
        type(setI), allocatable :: i(:) !< array to define reach different dVec's

end type setJ

!< type 2
type setUnk

        integer, allocatable :: jVec(:)  !< vector that stores the indices of "i"s of that "j" index
        
        integer :: maxJ !< size of the jVec 

        type(setJ), allocatable :: j(:) !< array to define reach different dVec's

end type setUnk


! type-1: bc%j%i%unkVec
! type-2: bc%unk%j%iVec
! either type1 or type2 branches are active
! default is type1 - type2 is to be used for u1,u2,u3..., v1,v2,v3 type problems.

type setBC

        ! type-1
        
        integer, allocatable :: jVec(:)  !< vector that stores the indices of "j"s of that "bc" index

        integer :: maxJ !< size of the jVec 

        type(setJ), allocatable :: j(:) !< array to define reach different jVec's

        ! type-2

        integer, allocatable :: unkVec(:)  !< vector that stores the indices of "j"s of that "bc" index

        integer :: maxUnk !< size of the jVec 

        type(setUnk),allocatable :: unk(:) !< array to define reach different uVec's

end type setBC



type setSet

        integer, allocatable :: bcVec(:) !< vector that stores the indices of "bc"s of that set

        type(setBC), allocatable :: bc(:) !< array to define reach different bc's
        
        integer maxBc  !< size of the bcVec 

end type setSet



type dom

        integer, allocatable :: nodeMatrix(:,:,:)

        integer, allocatable :: dofMatrix(:,:,:,:)

        integer, allocatable :: nodeVector(:,:)

        integer, allocatable :: dofVector(:,:)

        type(setSet) set

        integer, allocatable :: setVec(:)

        integer :: allBC

        integer :: minI

        integer :: minJ

        integer :: maxI

        integer :: maxJ

        integer :: nx

        integer :: ny

        integer :: numberOfGhostCells

        integer totalDof

        integer totalNode

        integer setDof !< if picard of set is used, then setDof might be less than expected

        real(8),pointer :: x(:) !< used as myModel%physics(p)%domain(d)%x=>myModel%physics(p)%x(lowerbound:upperbound)

        real(8),pointer :: xPTO(:) !< used as myModel%physics(p)%domain(d)%x=>myModel%physics(p)%x(lowerbound:upperbound)

        real(8),pointer :: xRTO(:) !< used as myModel%physics(p)%domain(d)%x=>myModel%physics(p)%x(lowerbound:upperbound)

        real(8),allocatable :: xCoor(:)
        
        real(8),allocatable :: yCoor(:)

        real(8) relXcoor
        
        real(8) relYcoor

        real(8) Lx

        real(8) Ly

        integer colorStartIndex

        integer colorEndIndex

        integer firstIndex

        integer lastIndex

        integer colorMin

        integer colorMax

        integer,allocatable :: dofVec(:)

        integer,allocatable :: gsSetVec(:)

        integer relativeIndex

        integer relativeIndexNode


end type dom


type phy

        real(8),pointer :: x(:) !< used as myModel%physics(p)%x=>discreteY(lowerbound:upperbound)

        real(8),pointer :: xPTO(:) !< used as myModel%physics(p)%x=>discreteY(lowerbound:upperbound)

        real(8),pointer :: xRTO(:) !< used as myModel%physics(p)%x=>discreteY(lowerbound:upperbound)

        type(dom),allocatable :: domain(:)

        integer :: numberOfDomains=1 !< number of domains. default is 1 (at least one domain is needed of course!)

        integer :: dof

        real(8) :: refXcoor !< x coordinate of the reference point
        
        real(8) :: refYcoor !< y coordinate of the reference point
        
        logical :: realTime=NO !< default setting is steady state. set to YES for transient problems like vof.

        logical :: pseudoTime=NO !< default setting is no pseudo time. set to YES if artificial compressilibity is used or pseudo time continuation

        logical :: setActive

        integer :: dFX

        character(79) :: tecVariables

        integer :: totalDof

        integer :: totalNode

        integer :: setDof

        integer, allocatable :: dVec(:)

        integer numberOfDomainGroups

        real(8),allocatable :: damping(:)

        integer,allocatable :: endIndex(:)

        real(8) :: Lx

        real(8) :: Ly

        integer :: numberOfPhysicalParameters
        
        real(8),allocatable :: physicalParameters(:)
        
        real(8),allocatable :: physicalParametersInitial(:)
        
        real(8),allocatable :: physicalParametersFinal(:)

        real(8) :: discOrder

        real(8),pointer :: evalOrder


end type phy



type dGroup

        integer,allocatable :: dVec(:)

end type dGroup



type pGroup

        type(dGroup),allocatable :: dg(:)

        integer, allocatable :: pVec(:)

end type pGroup


type DemonaModel

        type(phy),allocatable :: physics(:)

        integer :: numberOfPhysics !=1 !< number of physics. default is 1 (at least one physics is needed of course!)

        integer,allocatable :: pVec(:)

        real(8),allocatable :: x(:)

        real(8),allocatable :: tempx(:)

        real(8),allocatable :: xPTO(:)

        real(8),allocatable :: xRTO(:)

        integer totaldof

        integer totalNode

        integer,allocatable :: setVec(:)

        type(pGroup),allocatable :: pg(:)

        integer :: numberOfPhysicsGroups

        integer,allocatable :: endIndex(:)

end type DemonaModel


type(DemonaModel),target :: myModel








integer,target :: fullNx,fullNy,fullNz



! will be physics dependent - but for now it is added standalone
integer :: fV=0
integer :: ghostCells=2




integer,target :: currentSetLength
integer,target :: globalCurrentSetLength


!integer,parameter :: relativeModel=relativeBC
integer :: relativeModel=relativeLex




integer :: globalCorrection=0

logical :: setRuleActive=OFF


integer :: currentPhysics

integer :: currentDomain


! pointers

integer,pointer :: numberOfPhysics

integer,pointer :: numberOfDomains

integer,pointer :: dof

integer,pointer :: dFX

integer,pointer :: numberOfGhostCells

integer,pointer :: pVec(:)

integer,pointer :: dVec(:)

character(len=79) :: TecplotVariables

integer,pointer :: totalDof

integer,pointer :: modelDof

integer,pointer :: physicsDof

integer,pointer :: domainDof

integer,pointer :: totalNode

integer,pointer :: modelNode

integer,pointer :: physicsNode

integer,pointer :: domainNode

integer,pointer :: setVec(:)

integer,pointer :: localSetVec(:)

integer,pointer :: maxPhysicsGroup

integer,pointer :: maxDomainGroup

real(8),pointer :: refXcoor

real(8),pointer :: refYcoor

real(8),pointer :: Lx

real(8),pointer :: Ly

integer,pointer :: firstIndex

integer,pointer :: lastIndex

integer,pointer :: colorMin

integer,pointer :: colorMax

integer,pointer :: setDof

integer,pointer :: relativeIndex

integer,pointer :: numberOfPhysicalParameters


type(setSet), pointer :: discSet


type phyY

        real(8),pointer :: y

end type phyY

type(phyY), allocatable :: setY(:)




real(8),pointer :: y(:) !< used as y=>myModel%x

real(8),pointer :: yPTO(:) !< used as y=>myModel%x

real(8),pointer :: yRTO(:) !< used as y=>myModel%x

real(8),pointer :: tempy(:) !< used as tempy=>myModel%tempx

real(8),pointer :: discreteY(:)  !< used as discreteY=>y or discreteY=>tempy

real(8),pointer :: x(:) !< used as x=>myModel%physics(p)%domain(d)%x

real(8),pointer :: xPTO(:) !< used as x=>myModel%physics(p)%domain(d)%x

real(8),pointer :: xRTO(:) !< used as x=>myModel%physics(p)%domain(d)%x




integer colorFunctionIndex
integer colorValMax



real(8),allocatable :: yBackup(:)


!-----------------------
!< time variables


integer timeIter

real(8) dt


!-----------------------
!< performance variables

integer totalLinearIterations
integer totalNonLinearIterations


real(8), allocatable :: setTol(:)



!try to write a code that prints the default values.

!##############################
!SOLVER PARAMETERS
!##############################
!NEWTON
!------------------------------
logical :: freezedJacobian=off

integer :: maxNewtonIter=20 !< default is 20

integer :: newtonSuccess

logical :: inexactNewton=ON

real(8) :: jacSF=1.1d0   !< safety factor for the size of the jacVal

integer :: jacInterval=1

real(8) :: jacSkipTol=1e6       ! if set to 1.d0 then it is like no skip even the skip is on 
                                                        ! since normF=tol means problem already converged on the previous step
                                                        ! unless stepsize calculation is also perfomed,  then yes it changes !
                                                        ! meaning if f is converged then jac is constant, solve for better deltax only
logical :: jacSkip=off

integer :: maxBacktrackCount=10

logical :: doBacktrack=on

integer :: preconditioner=rp

logical :: bootStrapping=on

real(8) :: bootdamping=1.d0

real(8),target :: jacDiscOrder=lowOrder
!##############################
!PICARD
!------------------------------
real(8),pointer :: damping(:)

integer currentDof

integer :: maxPicardIter=50
!##############################
!GAUSS SEIDEL
!------------------------------
integer :: gsMaxNewton=1

integer :: maxGsIter=1000

logical :: blockGS=no

logical :: gsSolver=off

logical :: accumulating=no

integer :: maxGsPhysics=1
!##############################
!MULTIPLICATIVE SCHWARZ
!------------------------------
integer :: overlap=0

integer :: overlapX !< might be used

integer :: overlapY !< might be used

integer :: maxMSIter=100 !<default value

integer :: xPartition=1 

integer :: yPartition=1

integer :: numberOfMSDomains

logical :: msSolver=on

logical :: msRandom=off

logical :: msArtificialBC=on

real(8) :: msArtificialAlpha(2)  !<--------------------!!!!!!!!!!!!!!
!#######################################################################################
! model related parameters



integer :: maxRealTimeIter=1

integer :: maxPhysicsIter=1

integer :: maxPseudoTimeIter=1

integer :: maxDomainIter=1





integer :: mySolver=1

integer :: myHybrid=1


logical :: readInitialGuess=off

logical :: poissonCompatibility=on



!<dimensionless groups

!real(8), pointer :: Re
!real(8), pointer :: Ra
!real(8), pointer :: Pr

real(8) dtau




!debugging variables
integer test





end module global

---