|
[Sponsors] |
March 29, 2006, 04:43 |
Hi all,
I try to build in O
|
#1 |
Member
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 17 |
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 |
|
March 29, 2006, 11:01 |
In this case, the fluxes you'v
|
#2 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,907
Rep Power: 33 |
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 |
|
March 29, 2006, 11:47 |
Hi Hrv,
I know the operator
|
#3 |
Member
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 17 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Fluxes | selen | Siemens | 1 | August 24, 2003 17:23 |
Roe upwind fluxes!! | Magesh | Main CFD Forum | 2 | July 17, 2003 12:35 |
Fluxes validation | Mavinakere | FLUENT | 0 | May 23, 2002 00:24 |
Help need on fluxes | Zhu | Main CFD Forum | 7 | August 12, 2001 11:14 |
help needed on fluxes | zhu | Phoenics | 0 | August 7, 2001 22:35 |