CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   twoPhaseEulerFoam confined plunging jet simulation failed (https://www.cfd-online.com/Forums/openfoam-solving/228012-twophaseeulerfoam-confined-plunging-jet-simulation-failed.html)

ves June 17, 2020 04:53

twoPhaseEulerFoam confined plunging jet simulation failed
 
Hello, Foamers!
My confined plunging jet simulation in twoPhaseEulerFoam failed after 3-rd PIMPLE iteration on the first time step:
Boundary condition

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object epsilonm;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -3 0 0 0 0];

internalField uniform 1.5e-4;

boundaryField
{

opening
{
type fixedValue;
value $internalField;
}

inlet
{
type fixedValue;
value $internalField;
}

outlet
{
type inletOutlet;
phi phim;
inletValue $internalField;
value $internalField;
}

wall
{
type zeroGradient;
value $internalField;
}

wallinner
{
type zeroGradient;
value $internalField;
}
}

dimensions [0 2 -2 0 0 0 0];

internalField uniform 3.75e-5;

boundaryField
{

opening
{
type fixedValue;
value $internalField;
}

inlet
{
type fixedValue;
value $internalField;
}

outlet
{
type inletOutlet;
phi phim;
inletValue $internalField;
value $internalField;
}

wall
{
type zeroGradient;
value $internalField;
}

wallinner
{
type zeroGradient;
value $internalField;
}

}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 1e5;

boundaryField
{

opening
{
type calculated;
value $internalField;
}

inlet
{
type calculated;
value $internalField;
}

outlet
{
type calculated;
value $internalField;
}

wall
{
type calculated;
value $internalField;
}

wallinner
{
type calculated;
value $internalField;
}

}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{


opening
{
type totalPressure;
p0 uniform 0;
U U.air;
phi phi.air;
value uniform 0;
}

inlet
{
type fixedFluxPressure;
value uniform 0;
}

outlet
{
type fixedFluxPressure;
value prghPressure;
}

wall
{
type fixedFluxPressure;
value uniform 0;
}

wallinner
{
type fixedFluxPressure;
value uniform 0;
}

}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object Theta;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0.0;

boundaryField
{
inlet
{
type fixedValue;
value uniform 1.0e-7;
}

outlet
{
type inletOutlet;
inletValue uniform 1.0e-7;
value uniform 1.0e-7;
}

wall
{
type zeroGradient;
}

wallinner
{
type zeroGradient;
}

opening
{
type fixedValue;
value uniform 1.0e-7;
}
}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object alphat.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{

opening
{
type calculated;
value $internalField;
}

inlet
{
type calculated;
value $internalField;
}

outlet
{
type calculated;
value $internalField;
}

wall
{
type compressible::alphatWallFunction;
Prt 0.85;
value $internalField;
}

wallinner
{
type compressible::alphatWallFunction;
Prt 0.85;
value $internalField;
}


}


FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object epsilon.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -3 0 0 0 0];

internalField uniform 1.5e-4;

boundaryField
{
opening
{
type fixedValue;
value $internalField;
}

inlet
{
type fixedValue;
value uniform 0.000371;
}

outlet
{
type inletOutlet;
phi phi.air;
inletValue $internalField;
value $internalField;
}

wall
{
type epsilonWallFunction;
value $internalField;
}

wallinner
{
type epsilonWallFunction;
value $internalField;
}
}


FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object k.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 3.75e-5;

boundaryField
{

opening
{
type fixedValue;
value $internalField;
}

inlet
{
type fixedValue;
value $internalField;
}

outlet
{
type inletOutlet;
phi phi.air;
inletValue $internalField;
value $internalField;
}

wall
{
type kqRWallFunction;
value $internalField;
}

wallinner
{
type kqRWallFunction;
value $internalField;
}

}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object kappai.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 -1 0 0 0 0 0];

internalField uniform 2000;

boundaryField
{
wall
{
type zeroGradient;
}
wallinner
{
type zeroGradient;
}
opening
{
type inletOutlet;
phi phi.air;
inletValue $internalField;
value $internalField;
}
inlet
{
type zeroGradient;
}

outlet
{
type zeroGradient;
}


}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -1 0 0 0 0];

internalField uniform 1e-8;

boundaryField
{

opening
{
type calculated;
value uniform 0;
}

inlet
{
type fixedValue;
value uniform 1.3416e-05;
}

outlet
{
type calculated;
value $internalField;
}

wall
{
type nutkWallFunction;
value uniform 0;
}

wallinner
{
type nutkWallFunction;
value uniform 0;
}

}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 300;

boundaryField
{

opening
{
type fixedValue;
value $internalField;
}

inlet
{
type fixedValue;
value $internalField;
}

outlet
{
type inletOutlet;
phi phi.air;
inletValue $internalField;
value $internalField;
}


wall
{
type zeroGradient;
}

wallinner
{
type zeroGradient;
}
}

FoamFile
{
version 2.0;
format binary;
class volVectorField;
object U.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{

opening
{
//type fluxCorrectedVelocity;
type pressureInletOutletVelocity;
value uniform (0 0 0);
}

inlet
{
type fixedValue;
value uniform (0 0 0);
}
outlet
{
type pressureInletOutletVelocity;
phi phi.air;
value $internalField;
}

wall
{
type noSlip;

}

wallinner
{
type noSlip;
}
}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object alphat.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{

opening
{
type calculated;
value $internalField;
}

inlet
{
type calculated;
value $internalField;
}

outlet
{
type calculated;
value $internalField;
}

wall
{
type compressible::alphatWallFunction;
Prt 0.85;
value $internalField;
}

wallinner
{
type compressible::alphatWallFunction;
Prt 0.85;
value $internalField;
}


}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object epsilon.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -3 0 0 0 0];

internalField uniform 1.5e-4;

boundaryField
{
opening
{
type fixedValue;
value $internalField;
}

inlet
{
type fixedValue;
value uniform 0.000371;
}

outlet
{
type inletOutlet;
phi phi.water;
inletValue $internalField;
value $internalField;
}

wall
{
type epsilonWallFunction;
value $internalField;
}

wallinner
{
type epsilonWallFunction;
value $internalField;
}
}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object k.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 3.75e-5;

boundaryField
{

opening
{
type fixedValue;
value $internalField;
}

inlet
{
type fixedValue;
value $internalField;
}

outlet
{
type inletOutlet;
phi phi.water;
inletValue $internalField;
value $internalField;
}

wall
{
type kqRWallFunction;
value $internalField;
}

wallinner
{
type kqRWallFunction;
value $internalField;
}

}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object nut.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -1 0 0 0 0];

internalField uniform 1e-8;

boundaryField
{

opening
{
type calculated;
value uniform 0;
}

inlet
{
type fixedValue;
value uniform 1.3416e-05;
}

outlet
{
type calculated;
value $internalField;
}

wall
{
type nutkWallFunction;
value uniform 0;
}

wallinner
{
type nutkWallFunction;
value uniform 0;
}

}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 300;

boundaryField
{

opening
{
type fixedValue;
value $internalField;
}

inlet
{
type fixedValue;
value $internalField;
}

outlet
{
type inletOutlet;
phi phi.air;
inletValue $internalField;
value $internalField;
}


wall
{
type zeroGradient;
}

wallinner
{
type zeroGradient;
}
}

FoamFile
{
version 2.0;
format binary;
class volVectorField;
object U.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
opening
{
//type fluxCorrectedVelocity;
type pressureInletOutletVelocity;
value uniform (0 0 0);
}

inlet
{
type flowRateInletVelocity;
massFlowRate 0.205;
extrapolateProfile yes;
rho rho;
rhoInlet 1.0;
value uniform (0 0 0);
}

outlet
{
type flowRateInletVelocity;
massFlowRate 0.205;
extrapolateProfile yes;
rho rho;
rhoInlet 1.0;
value uniform (0 0 0);
}

wall
{
type noSlip;
}

wallinner
{
type noSlip;
}

}

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object setFieldsDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

defaultFieldValues
(
volScalarFieldValue alpha.air 1
volScalarFieldValue alpha.water 0
);

regions
(
boxToCell
{
box (-0.04 0 -0.04) (0.04 1.0 0.04);
fieldValues
(
volScalarFieldValue alpha.air 0
volScalarFieldValue alpha.water 1
);
}
);

ves June 17, 2020 04:55

solver settings and error
 
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
alpha.air
{
nAlphaCorr 1;
nAlphaSubCycles 2;
MULESCorr yes; // Switches on semi-implicit MULES
nLimiterIter 8; // Number of MULES iterations over the limiter
alphaApplyPrevCorr true;
}

p_rgh
{
solver GAMG;
smoother DIC;
tolerance 1e-8;
relTol 0;
}

p_rghFinal
{
$p_rgh;
relTol 0;
}

"U.*"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-7;
relTol 0;
minIter 1;
}

"e.*"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-7;
relTol 0;
minIter 1;
}

"(k|epsilon|Theta).*"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-7;
relTol 0;
minIter 1;
}
}

PIMPLE
{
nOuterCorrectors 3;
nCorrectors 1;
nNonOrthogonalCorrectors 1;
}

relaxationFactors
{
equations
{
".*" 0.1;
}
}






Create time

Create mesh for time = 0


PIMPLE: no residual control data found. Calculations will employ 3 corrector loops


Reading g

Reading hRef
Creating twoPhaseSystem

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

Calculating face flux field phi.air
Selecting diameterModel for phase air: IATE
Selecting turbulence model type RAS
Selecting RAS turbulence model mixtureKEpsilon
RAS
{
RASModel mixtureKEpsilon;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 1.92;
Cp 0.25;
sigmak 1;
sigmaEps 1.3;
}

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

Calculating face flux field phi.water
Selecting diameterModel for phase water: constant
Selecting turbulence model type RAS
Selecting RAS turbulence model mixtureKEpsilon
RAS
{
RASModel mixtureKEpsilon;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 1.92;
Cp 0.25;
sigmak 1;
sigmaEps 1.3;
}

Selecting default blending method: hyperbolic
Selecting aspectRatioModel for (air in water): constant
Selecting aspectRatioModel for (water in air): constant
Selecting dragModel for (air and water): segregated
Selecting swarmCorrection for (air and water): none
Selecting dragModel for (air in water): SchillerNaumann
Selecting swarmCorrection for (air in water): Tomiyama
Selecting dragModel for (water in air): SchillerNaumann
Selecting swarmCorrection for (water in air): Tomiyama
Selecting virtualMassModel for (air in water): constantCoefficient
Selecting virtualMassModel for (water in air): constantCoefficient
Selecting heatTransferModel for (air in water): RanzMarshall
Selecting heatTransferModel for (water in air): RanzMarshall
Selecting liftModel for (air in water): LegendreMagnaudet
Selecting liftModel for (water in air): LegendreMagnaudet
Selecting wallLubricationModel for (air in water): Frank
Selecting wallLubricationModel for (water in air): Frank
Selecting turbulentDispersionModel for (air in water): constantCoefficient
Selecting turbulentDispersionModel for (water in air): constantCoefficient
Calculating field g.h

Reading field p_rgh

Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 0 max: 0
Max Ur Courant Number = 0
Calculating field DDtU1 and DDtU2


Starting time loop

fieldAverage fieldAverage1:
Restarting averaging for fields:
U.air: starting averaging at time 0
U.water: starting averaging at time 0
alpha.air: starting averaging at time 0
p: starting averaging at time 0

Courant Number mean: 0 max: 0
Max Ur Courant Number = 0
deltaT = 6e-08
Time = 6e-08

PIMPLE: iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.09078137948 Min(alpha.air) = 0 Max(alpha.air) = 1
smoothSolver: Solving for kappai.air, Initial residual = 1, Final residual = 0.001264025234, No Iterations 1000
Constructing momentum equations
smoothSolver: Solving for e.air, Initial residual = 4.792804359e-05, Final residual = 4.211162828e-05, No Iterations 1000
smoothSolver: Solving for e.water, Initial residual = 0.4524787668, Final residual = 0.4138812702, No Iterations 1000
min T.air 300
min T.water 300
Selecting patchDistMethod meshWave
GAMG: Solving for p_rgh, Initial residual = 0.001385065939, Final residual = 2.795003837e-15, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 9.083520552e-16, Final residual = 9.083520552e-16, No Iterations 0
PIMPLE: iteration 2
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.09077876671 Min(alpha.air) = -11.41765647 Max(alpha.air) = 1
smoothSolver: Solving for kappai.air, Initial residual = 0.9999999892, Final residual = 3.872285928e-08, No Iterations 1
Constructing momentum equations
smoothSolver: Solving for e.air, Initial residual = 0.9999978023, Final residual = 1.312056595e-10, No Iterations 2
smoothSolver: Solving for e.water, Initial residual = 0.9999999985, Final residual = 7.404774966e-10, No Iterations 2
min T.air 299.2787486
min T.water 274.0068754
GAMG: Solving for p_rgh, Initial residual = 0.0001938205006, Final residual = 8.062402264e-14, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1.396477599e-16, Final residual = 1.396477599e-16, No Iterations 0
PIMPLE: iteration 3
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.09077953577 Min(alpha.air) = -2.486427675 Max(alpha.air) = 1

ves June 17, 2020 05:08

solvers setting from cfd-online and errors
 
with this settings https://www.cfd-online.com/Forums/op...tml#post607189


FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
"alpha.*"
{
nAlphaCorr 1;
nAlphaSubCycles 2;

smoothLimiter 0.1;

implicitPhasePressure yes;
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-5;
relTol 0;
minIter 1;
}

p_rgh
{
solver GAMG;
smoother DIC;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 500;
agglomerator faceAreaPair;
mergeLevels 1;
tolerance 1e-5;
relTol 0;
}

p_rghFinal
{
$p_rgh;
relTol 0;
}

"U.*"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-5;
relTol 0;
minIter 1;
}

"(h|e).*"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-5;
relTol 0;
minIter 1;
}

"Theta.*"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-5;
relTol 0;
minIter 1;
}

"(k|epsilon).*"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-5;
relTol 0;
minIter 1;
}

"epsilonFinal"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-08;
relTol 0;
}
}

PIMPLE
{
nOuterCorrectors 3;
nCorrectors 1;
nNonOrthogonalCorrectors 1;

pRefCell 0;
pRefValue 0;
residualControl
{
U
{
tolerance 1e-5;
relTol 0;
}
p
{
tolerance 5e-4;
relTol 0;
}
}
}

relaxationFactors
{
equations
{
".*" 1;
}
}



Courant Number mean: 0 max: 0
Max Ur Courant Number = 0
deltaT = 6e-08
Time = 6e-08

PIMPLE: iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.09078137948 Min(alpha.air) = 0 Max(alpha.air) = 1
smoothSolver: Solving for kappai.air, Initial residual = 1, Final residual = 0.0009141209896, No Iterations 1000
Constructing momentum equations
smoothSolver: Solving for e.air, Initial residual = 4.819361816e-05, Final residual = 4.231010327e-05, No Iterations 1000
smoothSolver: Solving for e.water, Initial residual = 0.4591847899, Final residual = 0.4243305836, No Iterations 1000
min T.air 300
min T.water 300
Selecting patchDistMethod meshWave
GAMG: Solving for p_rgh, Initial residual = 0.00148599746, Final residual = 6.291118185e-13, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 2.053537248e-13, Final residual = 2.053537248e-13, No Iterations 0
PIMPLE: iteration 2
MULES: Solving for alpha.air
MULES: Solving for alpha.air


--> FOAM FATAL ERROR:
cannot be called for a calculatedFvPatchField
on patch outlet of field alpha.air in file "/mnt/c/Users/admin/tutorial/confjet/confjettfef/0/alpha.air"
You are probably trying to solve for a field with a default boundary condition.

From function Foam::tmp<Foam::Field<Type> > Foam::calculatedFvPatchField<Type>::gradientIntern alCoeffs() const [with Type = double]
in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 188.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::calculatedFvPatchField<double>::gradientInte rnalCoeffs() const at ??:?
#3 Foam::fv::gaussLaplacianScheme<double, Foam::SymmTensor<double> >::fvmLaplacianUncorrected(Foam::GeometricField<do uble, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#4 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5 Foam::twoPhaseSystem::solve() at ??:?
#6 ? at ??:?
#7 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#8 ? at ??:?
Aborted (core dumped)


All times are GMT -4. The time now is 20:33.