[Chintan-21]

Hello Guys

I hope you all having the good weekends!!


I am simulating multi phase flow in the Helical pipe and for this I am using the damBreak case as reference. But I am getting the error as follow.



I have attached the ZIP file of the simulation files (file is double zipped due to size issue)



Thank you for your help in advance




Create time

Create mesh for time = 0

Selecting solver compressibleVoF
Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
properties liquid;
energy sensibleInternalEnergy;
}

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}

No MRF models present

Selecting surfaceTensionModel liquidProperties
Courant Number mean: 0 max: 0
Interface Courant Number mean: 0 max: 0
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
Creating fvModels from "constant/fvModels"

Selecting finite volume model type compressible::VoFTurbulenceDamping
Name: turbulenceDamping

PIMPLE: No convergence criteria found


PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode



Starting time loop

Courant Number mean: 0 max: 0
Interface Courant Number mean: 0 max: 0
deltaT = 1.19904e-05
Time = 1.19904e-05s

smoothSolver: Solving for alpha.water, Initial residual = 0, Final residual = 0, No Iterations 0
Phase-1 volume fraction = 0.5 Min(alpha.water) = 0.5 Max(alpha.water) = 0.5
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.5 Min(alpha.water) = 0.5 Max(alpha.water) = 0.5
smoothSolver: Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 1


--> FOAM FATAL ERROR:

request for surfaceScalarField phi.gas from objectRegistry region0 failed
available objects of type surfaceScalarField are

8
(
rhoPhi
ghf
alphaRhoPhi.water
phi
nHatf
magSf
alphaRhoPhi.air
alphaPhi.water
)


From function const Type& Foam:bjectRegistry::lookupObject(const Foam::word&) const [with Type = Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>]
in file /home/ubuntu/OpenFOAM/OpenFOAM-11/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 211.

FOAM aborting


Best regards
Chintan Prajapati

[Chintan-21]

Hello Guys


I hope everyone is doing well


I am working on the project of multi phase simulation of helical pipe, for this case I am using the multi phase Euler foam, but I am getting the following error



--> FOAM FATAL IO ERROR:
keyword pSat is undefined in dictionary "gas_liquid_inThe_gas/gas_liquid_inThe_gas"

file: gas_liquid_inThe_gas/gas_liquid_inThe_gas from line 115 to line 121.

From function const Foam::entry& Foam::dictionary::lookupEntry(const Foam::word&, bool, bool) const
in file db/dictionary/dictionary.C at line 827.

even after changing the phaseProperties file I am getting the same error




Here below i have pasted my phaseProperties file



/--------------------------------- C++ -----------------------------------\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 10
\\/ M anipulation |
\---------------------------------------------------------------------------/
FoamFile
{
format ascii;
class dictionary;
location "constant";
object phaseProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

type interfaceCompositionPhaseChangeMultiphaseSystem;

phases (gas liquid);

gas
{
type multiComponentPhaseModel;
diameterModel isothermal;
isothermalCoeffs
{
d0 3e-3;
p0 1e5;
}

residualAlpha 1e-6;
}

liquid
{
type multiComponentPhaseModel;
diameterModel constant;
constantCoeffs
{
d 1e-4;
}

residualAlpha 1e-6;
}

blending
{
default
{
type linear;
minFullyContinuousAlpha.gas 0.7;
minPartlyContinuousAlpha.gas 0.3;
minFullyContinuousAlpha.liquid 0.7;
minPartlyContinuousAlpha.liquid 0.3;
}

heatTransfer
{
type linear;
minFullyContinuousAlpha.gas 1;
minPartlyContinuousAlpha.gas 0;
minFullyContinuousAlpha.liquid 1;
minPartlyContinuousAlpha.liquid 0;
}

diffusiveMassTransfer
{
$heatTransfer;
}
}

surfaceTension
{
gas_liquid
{
type constant;
sigma 0.07;
}
}

drag
{
gas_dispersedIn_liquid
{
type SchillerNaumann;
residualRe 1e-3;
}

liquid_dispersedIn_gas
{
type SchillerNaumann;
residualRe 1e-3;
}
}

virtualMass
{
gas_dispersedIn_liquid
{
type constantCoefficient;
Cvm 0.5;
}

liquid_dispersedIn_gas
{
type constantCoefficient;
Cvm 0.5;
}
}

interfaceComposition
{
gas_liquid_inThe_gas
{
type saturated;
species ( water );
Le 1.0;
saturationPressure
{
type ArdenBuck;
pSat 3000; // I have updated pSat here
}
}

gas_liquid_inThe_liquid
{
type Henry;
species ( air );
k ( 1.492e-2 );
Le 1.0;

}
}

heatTransfer
{
gas_dispersedIn_liquid_inThe_gas
{
type spherical;
residualAlpha 1e-4;
}

liquid_dispersedIn_gas_inThe_gas
{
type RanzMarshall;
residualAlpha 1e-4;
}

gas_dispersedIn_liquid_inThe_liquid
{
type RanzMarshall;
residualAlpha 1e-4;
}

liquid_dispersedIn_gas_inThe_liquid
{
type spherical;
residualAlpha 1e-4;
}
}

diffusiveMassTransfer
{
gas_dispersedIn_liquid_inThe_gas
{
type spherical;
Le 1.0;
}

liquid_dispersedIn_gas_inThe_gas
{
type Frossling;
Le 1.0;
}

gas_dispersedIn_liquid_inThe_liquid
{
type Frossling;
Le 1.0;
}

liquid_dispersedIn_gas_inThe_liquid
{
type spherical;
Le 1.0;
}
}

phaseTransfer
{}

lift
{}

wallLubrication
{}

turbulentDispersion
{}

interfaceCompression
{}

// ************************* //




I hope someone can help me with this.



Thank you

Best regards
Chintan