Writing a basic complete solver: myVOFFoam, VOF/homogeneous mixture solver for two phases. p-V coupling for a basic mixture solver with gravitational forces is explained before coding. A case is set and then run.
Hello professor; I reproduced your code as you explained it. Then I tried to adopt another method as shown below. The two methods do not give identical results. Do you have an explanation? Thank you. ..... ..... ..... fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U) - fvm::laplacian(mu, U) - rho*g ); if (piso.momentumPredictor()) { solve(UEqn == -fvc::grad(p)); } // --- PISO loop while (piso.correct()) { volScalarField rAU(1.0/UEqn.A()); volVectorField HbyA(constrainHbyA(rAU*UEqn.H(), U, p)); surfaceScalarField phiHbyA ( "phiHbyA", fvc::flux(HbyA) + fvc::interpolate(rho*rAU)*fvc::ddtCorr(U, phi) ); ..... ..... .....
@@SantiagoMarquezD Yes you're right prof, the gravitational force term is included in the H operator. What I've found is that the results of the two approaches start out perfectly identical, but over time they diverge little by little. My question: in your experience, which of the two approaches is the most reliable? Thank you once again.