fluctuating pressure levels
Posted: 23 Dec 2009, 18:42
Hi there
I’ll be grateful if somebody gives me some information about the following part of FlowSolve subroutine. A book, paper or any other theoretical reference will be ideal for me. Furthermore, what is LocalNodes?
!------------------------------------------------------------------------------
! This hack is needed cause of the fluctuating pressure levels
!------------------------------------------------------------------------------
n = NSDOFs * LocalNodes
IF ( NonlinearRelax /= 1.0d0 ) THEN
IF ( CompressibilityModel == Incompressible ) THEN
s = FlowSolution(NSDOFs)
FlowSolution(NSDOFs:n:NSDOFs) = FlowSolution(NSDOFs:n:NSDOFs) - s
PSolution(NSDOFs:n:NSDOFs) = PSolution(NSDOFs:n:NSDOFs) - PSolution(NSDOFs)
END IF
FlowSolution(1:n) = (1 - NonlinearRelax) * PSolution(1:n) + &
NonlinearRelax * FlowSolution(1:n)
IF ( CompressibilityModel == Incompressible ) THEN
FlowSolution(NSDOFs:n:NSDOFs) = FlowSolution(NSDOFs:n:NSDOFs) + s
END IF
END IF
Thanks,
Diana Carlson, NY
I’ll be grateful if somebody gives me some information about the following part of FlowSolve subroutine. A book, paper or any other theoretical reference will be ideal for me. Furthermore, what is LocalNodes?
!------------------------------------------------------------------------------
! This hack is needed cause of the fluctuating pressure levels
!------------------------------------------------------------------------------
n = NSDOFs * LocalNodes
IF ( NonlinearRelax /= 1.0d0 ) THEN
IF ( CompressibilityModel == Incompressible ) THEN
s = FlowSolution(NSDOFs)
FlowSolution(NSDOFs:n:NSDOFs) = FlowSolution(NSDOFs:n:NSDOFs) - s
PSolution(NSDOFs:n:NSDOFs) = PSolution(NSDOFs:n:NSDOFs) - PSolution(NSDOFs)
END IF
FlowSolution(1:n) = (1 - NonlinearRelax) * PSolution(1:n) + &
NonlinearRelax * FlowSolution(1:n)
IF ( CompressibilityModel == Incompressible ) THEN
FlowSolution(NSDOFs:n:NSDOFs) = FlowSolution(NSDOFs:n:NSDOFs) + s
END IF
END IF
Thanks,
Diana Carlson, NY