[nishant_hull]

Hi Matvej

I have advanced a bit ahead with my last problem. The previous error is sorted out. it was because of the difference in meshing and boundary condition.

However I am still wondering with the question of different implementation of p-equation. I mean..
when we use phi in pEqn (), we use
phi += pEqn.flux();
And when we use phid etc then we use
phi = pEqn.flux();

I am still wondering why?? Any hint please?

regards,
Nishant

[mkraposhin]

Hi, I'm very sorry for late reply (too much work).
Anyway, if this answer is still needed:

when we update velocity (phi = velocity*density*area) in diffusion pressure equation, we mean that dU ~ const * grad(dp). If we are using convection-diffusion equation for pressure: ddt(psi,p) + div(phi_d, p) - div(rho/A*grad(p)), we are solving for mass flux, which satisfies continuity equation.

Another words:
- diffusion equation - correction to mass fluxes (to satisfy continuity) - phi+=pEqn.flux()
- convection-diffusion - mass fluxes, which continuity - phi=pEqn.flux()

I hope, this helps you. I'm really sorry for late post

[Annier]

i am new to OpenFOAM.
i was going through the "OpenFOAM Programming Tutorial" written by Tommaso Lucchini. i modified icoFoam to icoScalarTransportFoam and cavity to cavityScalarTransport in accordance to the directions of the notes.
now i cannot set up the case cavityScalarTransport and cannot run it. would you please guide me with the commands for running this new application.thanking you in advance.