[Rachel]

@Hamed : If you are upgrading to OpenSUSE 11.1 for OF-1.6 it is not necessary. I am using SLED 10 till now. paraFoam does not work...

@ Sandy: I could not find interPhaseChangeFoam tutorial in OF-1.6

[Rachel]

Anybody managed to get rhoReactingFoam working ?

I am still having problems with this solver... as mentioned in this thread.
rhoReactingFoam

[chiven]

Hi, guys, how about your progresses on the two phase flow simulation? I am trying to use OpenFOAM to simulate the two phase flow with condensation (vapor to water). I think it is difficult, but very interesting.

Best regards,
Chiven

[sanjibdsharma]

Hi Matvej,

I am new to OpenFoam. Having started with the basics 2 months ago, I am now stuck with modeling of steady-state, liquid-phase, multi-species transport with chemical reaction. Can you provide some input as to how I proceed i.e., which modules should I use.

As far as I understood, I compiled with the following:

simpleFoam solver + basicMixture.H + basicThermo.H + liquidMixture.H + liquid.H + pEqn.H + UEqn.H

I have tried to use YEqn.H, but gave error. How to combine this and proceed ? It will be great if you have any example.

Regards,

Sanjib

[isabel]

I have a doubt about interPhaseChangeFoam.
The momentum equation is:

fvVectorMatrix UEqn
(
fvm::ddt(rho, U)
+ fvm::div(rhoPhi, U)
- fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
- fvm::laplacian(muf, U)
- (fvc::grad(U) & fvc::grad(muf))
);

UEqn.relax();

if (momentumPredictor)
{
solve
(
UEqn
==
fvc::reconstruct
(
(
fvc::interpolate(interface.sigmaK())*fvc::snGrad(g amma)
- ghf*fvc::snGrad(rho)
- fvc::snGrad(pd)
) * mesh.magSf()
)
);
}


Why the red line is ghf*fvc::snGrad(rho) ???
I think it would be ghf*rho

[sandy]

Hi isabel,

Do you understand the blue line?? How to derive it?

fvVectorMatrix UEqn
(
fvm::ddt(rho, U)
+ fvm::div(rhoPhi, U)
- fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
- fvm::laplacian(muf, U)
- (fvc::grad(U) & fvc::grad(muf))
);
UEqn.relax();

Thanks.

[isabel]

I don't know. For me fvc::ddt(rho) + fvc::div(rhoPhi) is the continuity equation, which is equal to zero.

In interFoam (a similar solver), you don't have this blue line:

fvVectorMatrix UEqn
(
fvm::ddt(rho, U)
+ fvm::div(rhoPhi, U)
- fvm::laplacian(muf, U)
- (fvc::grad(U) & fvc::grad(muf))
);

if (momentumPredictor)
{
solve
(
UEqn
==
fvc::reconstruct
(
(
fvc::interpolate(interface.sigmaK())*fvc::snGrad(g amma)
- ghf*fvc::snGrad(rho)
- fvc::snGrad(pd)
) * mesh.magSf()
)
);




But I still don't know why the term ghf*rho doesn't appear and what the term snGgrad(rho) means.

[lucky1729]

@isabel

snGrad(rho) calculates the normal component of gradient to a surface. Momentum equation in differential form requires 'rho*g' in right hand side. So to get that 'ghf'{g*hf} is multiplied by snGrad(rho).

[isabel]

Thanks lucky. But I still think that I don't need grad(rho), so I have modified the UEqn.H as follows. Is this way wrong???

surfaceScalarField muf = twoPhaseProperties.muf();
volVectorField gravityterm = rho*vector(0, -9.81, 0); // this means rho*g

fvVectorMatrix UEqn
(
fvm::ddt(rhoU)
+ fvm::div(rhoPhi, U)
- fvm::laplacian(muf, U)
- (fvc::grad(U) & fvc::grad(nuf))
- gravityterm
);

UEqn.relax();

if (momentumPredictor)
{
solve
(
UEqn
==
fvc::reconstruct
(
fvc::interpolate(interface.sigmaK())*fvc::snGrad(g amma)

- fvc::snGrad(pd)
) * mesh.magSf()
);
}

[isabel]

I already know why the term fvc::snGrad(rho) appears. It is a consequence of the UEqn equation works with dynamic pressure (pd) instead total pressure (p).


But I still have a question:


What is the advantage of working with dynamic pressure (pd) instead total pressure (p)

[haghajani]

Hi Rachel,
Have you got new experiences with rhoReactingFoam, to simulate phase change of a liquid into another medium?
Have you find out sth from counterFlowFlame2D?
Please keep us posted,

Hamed

[haghajani]

Originally Posted by mkraposhin

Here is a link for source code of caviatingPimpleFoam solver for OpenFOAM 1.5

http://www.os-cfd.narod.ru/small_fil...leFoam_1.5.tgz

Dear Matvej,

In cavitatingPimpleFoam/pEqn, Line 17, phid2 = fvc::interpolate(psi) * (phid1 + fvc::ddtPhiCorr(rUA, rho, U, phi) / fvc::interpolate(rho)), the two sides of "+" are inconsistent because of red bold part. Am I right?

Regarding the thermodynamicProperties, beside the items indicated in rascavitatingfoam tutorial, there should be some data for following parameters, according the chosen gammaModel, say Kunz;
CgammaP [0 0 0 0 0]
CgammaA [0 0 0 0 0]
UInf [0 1 -1 0 0] ?
tauInf [0 0 1 0 0] ?
pCoeff1 [0 0 0 0 0] ?
CgammaEff [0 0 0 0 0] ?
Would you please drop a line on their Initiation, values.

I am facing some other incompatible dimension errors, which may share with you here, later.

Best regards,
Hamed

[sandy]

Hi hamed, you think, what is the difference between the interPhaseChangeFoam solver with the cavitatingPimpleFoam solver? Thanks.

[haghajani]

Sandy,

As far as I have understood, cavitatingPimpleFoam is a Homogenous cavitation solver which uses PISO-SIMPLE coupling loop. Homogeneous might means there is no slip velocity between liquid and its cavitated phase. It is compressible, so uses a barotropic EOS, to deal with density. It solves for one phase wich might result in two phase [first phase + vapour/liquid of first Phase]
On the other hand, interPhaseChangeFoam considers the Phase change of, say a liquid as water, into another phase, say Air. so at the end we may have two phases, (liquid and [vapour of that liquid + Air]). and the solver is incompressible here.

Hope it would be a step forward,

Best regards,
Hamed

[neewbie]

hello,

i was pretty happy to find these post. since i´m trying to understand the rasCavtiation-solver in OF1.5. Like usual, i was trying to read through the code and got stuck while reading the rhoEqn.H.

I think my question is: what is this phiv doing? It is some sort of correction? I read
phiv=phi/fvc::interpolate(rho) but i didn´t got the grasp.
I found two parts for correction in pEqn.H: resetphiPatches and resetphivpatches... but i have no clou what they are doing.

Could somebody please be so kind and give me a hint?

Thanks in advance,

neewbie

[mach]

Hi,

Does any body used moment method for aerosol coagulation of size less then 0.1 um in Fluent?.

Does any know about moment method UDF?.

Please let me know

Thanks

[matejfor]

John, sending the same question to this forum as well as to Fluent one is not a good thing. This forum is for OpenFOAM. It will make people annoyed, not attracted to answer your question.

good luck
matej

[nimasam]

hi openFoam former
im going to work on evaporation of incompressible flow, can any body send me a simple model which can work properly with openFoam 1.6 solver (interPhaseChangeFoam) ?
nimasam65@gmail.com

[sandy]

No, I guess interPhaseChangeFoam is not suiltable to your problem, because its mode is about the specified cavitation flows, and they are some empirical formula to defined the phase change, but not really solve equations of the heat and mass transfer. Moreover, this phase change is an isothermal process.

If you need the tutorial about this solver, you can find it in OF 1.6.x

[nimasam]

thanks sandy
is there any solver in openFoam which is appropriate for my simulation?
if no , whats ur suggestion for further progress? i dont know use whether interfoam (direct model) or bubbleFoam (Eulerian _Eulerian model) or interChangePhaseFoam as a basic solver and then add mass and heat transfers equation.
could you give me some hints
thanks