[su_junwei]

Hi Foamers
The following thread discussed the implementation of compressibleInterFoam.
http://www.cfd-online.com/Forums/ope...interfoam.html

I am confused with the implementation of dgdt in the file pEqn.H
In this thread, we can get.
dgdt = { ( psi2 / rho2) - (psi1 / rho1) } DpDt
However, in the implementation, it gives

dgdt = (pos(alpha2)*psi2/rho2 - pos(alpha1)*(psi1/rho1))*(p_rghEqnComp & p_rgh);

Is this right? shoud it be
dgdt = (psi2/rho2 - psi1/rho1))*(p_rghEqnComp & p_rgh);

where are alpha1 and alpha2 from?

would you please give me some tips?

Best Regards, Junwei

[su_junwei]

I got it. pos is just a sign function.

[krishtej23]

Originally Posted by su_junwei

I got it. pos is just a sign function.

Hi,

I am trying to run a simulation using compressibleinterfoam application. I changed the alpha1 to alphaair, alphawater and alphadrop by creating each file separately. But when I try to run the code I get an error telling:


/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : compressibleInterFoam
Date : May 07 2014
Time : 17:32:26
Host : "mountaineer"
PID : 15228
Case : /auto/scratch/mdinc/ottawa_simulations/new_corrected_B_cases/3D_simulations/earth_3b_of
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g

PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Reading field alpha1



--> FOAM FATAL IO ERROR:
cannot find file

file: /auto/scratch/mdinc/ottawa_simulations/new_corrected_B_cases/3D_simulations/earth_3b_of/0/alpha1 at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 73.

FOAM exiting


It is appreciable if you can help me in solving this.
Thank you.