Indeed!
The "total pressure" is defined in Code_Saturne as being: Ptot = Pred + rho0*g*(z-z0) + P0, where Pred is the so-called reduced pressure (i.e. the resolved pressure). If you want the dynamic pressure, you'll need to add 0.5*rho*U**2. You can do that by adding the usvpst.f90 subroutine in which you can defined your own variables, on a given mesh (ipart = -1 for the volumic mesh).
David
ps: In the formula for the total pressure, I've intentionnally forgotten the reference of the reduced pressure Pred0.
Setting non standard boundary conditions in simple calculation
Forum rules
Please read the forum usage recommendations before posting.
Please read the forum usage recommendations before posting.
Re: Setting non standard boundary conditions in simple calculation
Hi everybody,
I have the same issue as Thomas's. I started a simulation with a velocity inlet, and now, I want a pressure inlet. So I changed the usclim.f90 file, and now it seems Code_Saturne only takes into account the velocity inlet I set for the first time.
Here is my code:
I really don't know what I am doing wrong, so if you have any idea...
Best regards
Stéphane
I have the same issue as Thomas's. I started a simulation with a velocity inlet, and now, I want a pressure inlet. So I changed the usclim.f90 file, and now it seems Code_Saturne only takes into account the velocity inlet I set for the first time.
Here is my code:
Code: Select all
iphas=1
call getfbr('inlet', nlelt, lstelt)
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
iel = ifabor(ifac)
itypfb(ifac,iphas) = ientre
icodcl(ifac,ipr(iphas))=1.d0
rcdocl(ifac,ipr(iphas),1)=380.d0
icodcl(ifac,iu(iphas)) = 3.d0
icodcl(ifac,iv(iphas)) = 3.d0
icodcl(ifac,iw(iphas)) = 3.d0
enddo
Best regards
Stéphane
Re: Setting non standard boundary conditions in simple calculation
For sanity check, Stephane's question has been answered in a different thread (see https://code-saturne.info/products/code-saturne/forums/general-usage/817624553 )
;-)
;-)