CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Understanding implemented PISO (

stephan January 28, 2006 12:11

hi, does anybody know why t

does anybody know why the flux in the icoFoam-case is calculated as it is?
there are lines like that:
volScalarField A=UEqn.A();
phi=(fvc::interpolate(U+UphiCoeff*fvc::ddt0(U)/A,"interpolate((H(U)|A(U))")&mesh .Sf())-UphiCoeff*fvc::interpolate(1/A)*fvc::ddt0(phi);

where phi itself is "linearinterpolate(U)&mesh.Sf()"
i know that i need a div-term for the pressure-eqn but what are these other terms in the phi-calculation? for me it looks like adding something in the interpolation and than substracting the same after interpolation.
Am i right, that ddt0(U) gives the old-time-value for U in discretizised form (with coefficient)?
I hope somebody can explain what happens here.
thanks in advance

All times are GMT -4. The time now is 15:29.