|
[Sponsors] |
May 18, 2016, 20:33 |
barotropicCavitatingFoam
|
#1 |
New Member
|
Hi Foamers,
I am new to OpenFOAM and I have a specific question regarding barotropicCavitatingFoam in Foam extend. I haven't seen many discussion regarding this solver in this forum. Can somebody explain to me about the pEqn.H of the solver? the code looks like this Code:
{ surfaceScalarField rUAf = 1.0/fvc::interpolate(UEqn.A()); U = UEqn.H()/UEqn.A(); phi = fvc::interpolate(U) & mesh.Sf(); // Bug fix: must change name of phi on copy to keep objectRegistry happy // HJ, 7/Nov/2010 surfaceScalarField phiU ( "phiU", phi ); surfaceScalarField phip = fvc::interpolate(psiByRho)*phi; for(int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++) { fvScalarMatrix pEqn ( fvm::laplacian(rUAf, p) ); solve ( psiByRho*fvm::ddt(p) + fvm::div(phip, p) - fvm::Sp(fvc::div(phip), p) + fvc::div(phi) - pEqn ); if (nonOrth == nNonOrthCorr) { phi -= pEqn.flux(); } } U += fvc::reconstruct(phi - phiU); U.correctBoundaryConditions(); } How is it different from the cavitatingFoam Solver? |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|