|
[Sponsors] |
July 27, 2007, 04:33 |
Hi everybody,
From icoFoam:
|
#1 |
Senior Member
Thomas Jung
Join Date: Mar 2009
Posts: 102
Rep Power: 17 |
Hi everybody,
From icoFoam: fvVectorMatrix UEqn ( fvm::ddt(U) + fvm::div(phi, U) - fvm::laplacian(nu, U) ); solve(UEqn == -fvc::grad(p)); // --- PISO loop for (int corr=0; corr<nCorr; corr++) { volScalarField rUA = 1.0/UEqn.A(); U = rUA*UEqn.H(); phi = (fvc::interpolate(U) & mesh.Sf()) + fvc::ddtPhiCorr(rUA, U, phi); adjustPhi(phi, U, p); for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++) { fvScalarMatrix pEqn ( fvm::laplacian(rUA, p) == fvc::div(phi) ); pEqn.setReference(pRefCell, pRefValue); pEqn.solve(); if (nonOrth == nNonOrthCorr) { phi -= pEqn.flux(); } } # include "continuityErrs.H" U -= rUA*fvc::grad(p); U.correctBoundaryConditions(); } Question: Inside the loop UEqn is not beeig re-build. Should rUA, or UEqn.A() and UEqn.H() not remain the same? Why then is this beein repeated all the time? volScalarField rUA = 1.0/UEqn.A(); U = rUA*UEqn.H(); Thank you very much for any hint! Thomas |
|
July 27, 2007, 04:52 |
hi,
you should search the f
|
#2 |
Senior Member
Stephan Gerber
Join Date: Mar 2009
Location: Germany
Posts: 118
Rep Power: 17 |
hi,
you should search the forum - this question was already answered. greeting stephan |
|
July 27, 2007, 05:25 |
finally found it, sorry...
To
|
#3 |
Senior Member
Thomas Jung
Join Date: Mar 2009
Posts: 102
Rep Power: 17 |
finally found it, sorry...
To whom it may interest,too: By Hrvoje Jasak on Thursday, August 04, 2005 - 12:40 pm: Edit Post Not so stupid: UEqn.H() uses the current value of U to evaluate itself and that's where the predicted value gets into the game. |
|
July 27, 2007, 07:02 |
hi,
no problem - nice if i
|
#4 |
Senior Member
Stephan Gerber
Join Date: Mar 2009
Location: Germany
Posts: 118
Rep Power: 17 |
hi,
no problem - nice if i could help you ... greetings stephan |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
PISO Details | Dr B.M. Smith (Smith) | OpenFOAM Running, Solving & CFD | 17 | November 17, 2009 12:31 |
Understanding k from kepsilon | markh83 | OpenFOAM Post-Processing | 3 | December 5, 2008 04:42 |
Incompressible boussiesq flow pressure BC piso and continuity problems | tehache | OpenFOAM Running, Solving & CFD | 4 | November 6, 2008 04:08 |
LUSGS understanding | Gustavo | Main CFD Forum | 0 | February 14, 2007 11:12 |
Understanding implemented PISO | stephan | OpenFOAM Running, Solving & CFD | 0 | January 28, 2006 12:11 |