Setting non standard boundary conditions in simple calculation

[David Monfort]

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.

[Stéphane Blanchet]

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:

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 

I really don't know what I am doing wrong, so if you have any idea...

Best regards

Stéphane

[David Monfort]

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 )
;-)