CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Update Fluxes

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 29, 2006, 03:43
Default Hi all, I try to build in O
  #1
Member
 
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 8
luca is on a distinguished road
Hi all,

I try to build in OpenFOAM a pressure based algorithm for prediciting flows at all speed. I formulate pressure equation in term of pressure-correction like in original SIMPLE procedure.
After I've solved pressure-correction equation i need to update fluxes in order to satisfy local continuity equation.
I use this empirical formulation after I've corrected volScalarField "rho", "p" and "U":

--
p.correctBoundaryConditions();

phi = fvc::interpolate(rho*U) & mesh.Sf();
phi.boundaryField() = ( rho.boundaryField()*U.boundaryField() ) & mesh.Sf().boundaryField();
--

Is it right?!

This procedure works but I'm not sure about that..


Thanks in advance at all OpenFOAM users
luca is offline   Reply With Quote

Old   March 29, 2006, 10:01
Default In this case, the fluxes you'v
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,761
Rep Power: 21
hjasak will become famous soon enough
In this case, the fluxes you've got are not conservative, unless you are very very careful with all the details on interpolation and consistency between the face flux formula in the pressure equation and the interpolation of rho*U on the face (and the expression above says you're not) :-)

Try using pEqn.flux() - if the p equation has been formed correctly, you will get the consistent and conservative flux from it. For details of what flux() member function does, have a dig through the archive or the source code.

Enjoy,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   March 29, 2006, 10:47
Default Hi Hrv, I know the operator
  #3
Member
 
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 8
luca is on a distinguished road
Hi Hrv,

I know the operator pEqn.flux(), in my case I use fvScalarMatrix "pEqn" to solve pressure correction scalar only. Following the Peric's formulation for compressible flows I need to correct the fluxes also with explicitly density corrector. Does pEqn.flux() work properly in this case!?

I hope you can help me another time

Luca
luca is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Fluxes selen CD-adapco 1 August 24, 2003 16:23
Roe upwind fluxes!! Magesh Main CFD Forum 2 July 17, 2003 11:35
Fluxes validation Mavinakere FLUENT 0 May 22, 2002 23:24
Help need on fluxes Zhu Main CFD Forum 7 August 12, 2001 10:14
help needed on fluxes zhu Phoenics 0 August 7, 2001 21:35


All times are GMT -4. The time now is 14:03.