multiphaseEulerFoam FOAM FATAL IO ERROR

qutadah.r · December 7, 2021, 07:25

Hi everyone, I am having the following error when running a two-phase flow in OpenFOAM.

The original case works perfectly and doesnt have any issues, I come across this issue when i change some values of alpha or U in the boundary conditions...
* I am using collated filehandler and parallel processing using two nodes, 48 processors each.

Code:

--> FOAM FATAL IO ERROR:
Expected a ')' while reading binaryBlock, found on line 21 the word 'e▒L▒3'

file: /rwthfs/rz/cluster/hpcwork/jg235250/SetupsMarkus/EE2F-MUSIG/06_3D_Pumpe_M1_v6_Qtot696_IGVF_5_SST_OF9_EE2F_MUSIG/0.288173000000247/U.water at line 21.

    From function Foam::Istream &Foam::Istream::readEnd(const char *)
    in file db/IOstreams/IOstreams/Istream.C at line 109.

Any help would be appreciated! Thanks!

qutadah.r · December 7, 2021, 13:25

so i dont know what caused this still, i was using nano editor to edit the binary boundary files... I just copied the original case again and used vi editor and it worked.... Even though i did the same and edited the binary files with nano and got this error .....

If someone knows if nano editor messes up the binary files somehow please let me kindly know, thanks

dakaezekoye · December 8, 2021, 12:56

I am experiencing almost similar issues, in my case I swapped out the cavitating bullet geometry with a 3D hydrofoil and I tried to carefully change the boundary conditions. The mesh was ok but when I ran the command "interPhaseChangeFoam" I got the error below
Please can anyone help!

engr@engr-dave:~/foam/dave-4.0/run/cavitatingBullet$ of8
engr@engr-dave:~/foam/dave-4.0/run/cavitatingBullet$ interPhaseChangeFoam
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 8
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 8-1c9b5879390b
Exec : interPhaseChangeFoam
Date : Dec 09 2021
Time : 00:16:41
Host : "engr-dave"
PID : 13384
I/O : uncollated
Case : /home/engr/foam/dave-4.0/run/cavitatingBullet
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

PIMPLE: No convergence criteria found

PIMPLE: No corrector convergence criteria found
Calculations will do 3 corrections

PIMPLE: Operating solver in transient mode with 3 outer correctors

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Creating phaseChangeTwoPhaseMixture

Selecting phaseChange model SchnerrSauer
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
bounding k, min: 0 max: 0.06 average: 0.06
bounding epsilon, min: 0 max: 0.0495 average: 0.0495
RAS
{
model kEpsilon;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

Reading g

Reading hRef
Calculating field g.h

No finite volume options present
GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.40042e-07, global = 5.21498e-22, cumulative = 5.21498e-22
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5 Foam::RASModels::kEpsilon<Foam::IncompressibleMome ntumTransportModel<Foam::transportModel> >::correctNut() at ??:?
#6 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/interPhaseChangeFoam"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/interPhaseChangeFoam"
Floating point exception (core dumped)
engr@engr-dave:~/foam/dave-4.0/run/cavitatingBullet$

night-hawk · December 8, 2021, 21:48

Quote:

Originally Posted by dakaezekoye

I am experiencing almost similar issues, in my case I swapped out the cavitating bullet geometry with a 3D hydrofoil and I tried to carefully change the boundary conditions. The mesh was ok but when I ran the command "interPhaseChangeFoam" I got the error below
Please can anyone help!

engr@engr-dave:~/foam/dave-4.0/run/cavitatingBullet$ of8
engr@engr-dave:~/foam/dave-4.0/run/cavitatingBullet$ interPhaseChangeFoam
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 8
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 8-1c9b5879390b
Exec : interPhaseChangeFoam
Date : Dec 09 2021
Time : 00:16:41
Host : "engr-dave"
PID : 13384
I/O : uncollated
Case : /home/engr/foam/dave-4.0/run/cavitatingBullet
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

PIMPLE: No convergence criteria found

PIMPLE: No corrector convergence criteria found
Calculations will do 3 corrections

PIMPLE: Operating solver in transient mode with 3 outer correctors

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Creating phaseChangeTwoPhaseMixture

Selecting phaseChange model SchnerrSauer
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
bounding k, min: 0 max: 0.06 average: 0.06
bounding epsilon, min: 0 max: 0.0495 average: 0.0495
RAS
{
model kEpsilon;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

Reading g

Reading hRef
Calculating field g.h

No finite volume options present
GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.40042e-07, global = 5.21498e-22, cumulative = 5.21498e-22
#0 Foam::error:

rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field&, Foam::UList const&, Foam::UList const&) at ??:?
#4 Foam::tmp Foam:

perator/(Foam::tmp const&, Foam::GeometricField const&) at ??:?
#5 Foam::RASModels::kEpsilon::correctNut() at ??:?
#6 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/interPhaseChangeFoam"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/interPhaseChangeFoam"
Floating point exception (core dumped)
engr@engr-dave:~/foam/dave-4.0/run/cavitatingBullet$

Try changing nOuterCorrectors nNonOrthonormalCorrectors.

You can also try changing schemes.

dakaezekoye · December 9, 2021, 02:46

Quote:

Originally Posted by night-hawk

Try changing nOuterCorrectors nNonOrthonormalCorrectors.

You can also try changing schemes.

Thank you night-hawk, I will try your suggestions ASAP

dakaezekoye · December 9, 2021, 07:12

Quote:

Originally Posted by night-hawk

Try changing nOuterCorrectors nNonOrthonormalCorrectors.

You can also try changing schemes.

I have tried to increase and decrease the value of the nOuterCorrectors nNonOrthogonalCorrectors, I have also tried to change several values in the fvscheme file but nothing is working.

I also tried to calculate using the 2D version of the same hydrofoil and the calculation started running but after a while, it terminated with an error message as shown below.

Please I need your help with this problem

Courant Number mean: 1.31656e-06 max: 5.03315
Interface Courant Number mean: 3.78911e-07 max: 5.03315
deltaT = 1.85645e-46
Time = 3.31999e-06

PIMPLE: Iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 7.9757e-08, Final residual = 2.59373e-12, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113975 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113975 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.784898, Final residual = 0.0776456, No Iterations 6
smoothSolver: Solving for Uy, Initial residual = 0.784338, Final residual = 0.0728061, No Iterations 8
GAMG: Solving for p_rgh, Initial residual = 2.50221e-33, Final residual = 2.50221e-33, No Iterations 0
PIMPLE: Iteration 2
smoothSolver: Solving for alpha.water, Initial residual = 1.83132e-08, Final residual = 1.94545e-12, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113974 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113974 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.673668, Final residual = 0.0639442, No Iterations 9
smoothSolver: Solving for Uy, Initial residual = 0.627473, Final residual = 0.0565589, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 2.08326e-32, Final residual = 2.08326e-32, No Iterations 0
PIMPLE: Iteration 3
smoothSolver: Solving for alpha.water, Initial residual = 2.30372e-08, Final residual = 3.42331e-12, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113974 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113974 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.575274, Final residual = 2.32922e-06, No Iterations 1000
smoothSolver: Solving for Uy, Initial residual = 0.649715, Final residual = 2.30958e-06, No Iterations 1000
GAMG: Solving for p_rgh, Initial residual = 6.85021e-34, Final residual = 6.85021e-34, No Iterations 0
smoothSolver: Solving for epsilon, Initial residual = 6.91275e-18, Final residual = 9.62875e-19, No Iterations 1
bounding epsilon, min: -1.24782e+114 max: 2.56052e+128 average: 6.87117e+124
smoothSolver: Solving for k, Initial residual = 0.938039, Final residual = 1.12015e-05, No Iterations 1000
bounding k, min: -3.96646e+78 max: 8.53537e+82 average: 3.25647e+80
ExecutionTime = 473.36 s ClockTime = 473 s

Courant Number mean: 4.39296e-07 max: 18.162
Interface Courant Number mean: 6.7485e-08 max: 1.3656
deltaT = 5.11082e-47
Time = 3.31999e-06

PIMPLE: Iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 1.24025e-08, Final residual = 3.75624e-14, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.873234, Final residual = 0.044149, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.956389, Final residual = 0.0385667, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 2.12784e-32, Final residual = 2.12784e-32, No Iterations 0
PIMPLE: Iteration 2
smoothSolver: Solving for alpha.water, Initial residual = 5.95103e-09, Final residual = 1.06362e-12, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.885566, Final residual = 0.0729296, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.951138, Final residual = 0.0617451, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 1.44958e-31, Final residual = 1.44958e-31, No Iterations 0
PIMPLE: Iteration 3
smoothSolver: Solving for alpha.water, Initial residual = 4.63717e-07, Final residual = 2.58962e-11, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.905716, Final residual = 9.8713e-07, No Iterations 210
smoothSolver: Solving for Uy, Initial residual = 0.942901, Final residual = 9.99183e-07, No Iterations 169
GAMG: Solving for p_rgh, Initial residual = 1.77193e-31, Final residual = 1.77193e-31, No Iterations 0
smoothSolver: Solving for epsilon, Initial residual = 9.18298e-13, Final residual = 1.07957e-18, No Iterations 1
bounding epsilon, min: -2.89875e+109 max: 1.91737e+128 average: 7.25083e+124
smoothSolver: Solving for k, Initial residual = 0.960219, Final residual = 0.000947622, No Iterations 1000
bounding k, min: -4.85718e+78 max: 2.00365e+87 average: 9.3082e+82
ExecutionTime = 473.82 s ClockTime = 473 s

Courant Number mean: 1.11795e-05 max: 760.236
Interface Courant Number mean: 1.34447e-08 max: 0.326501
deltaT = 3.36133e-49
Time = 3.31999e-06

PIMPLE: Iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 3.14423e-10, Final residual = 2.10051e-16, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.898913, Final residual = 0.0463407, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.965174, Final residual = 0.0949937, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 2.03761e-30, Final residual = 2.03761e-30, No Iterations 0
PIMPLE: Iteration 2
smoothSolver: Solving for alpha.water, Initial residual = 3.8785e-10, Final residual = 2.14871e-15, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.945497, Final residual = 0.0425314, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.890285, Final residual = 0.0789811, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 2.32038e-29, Final residual = 2.32038e-29, No Iterations 0
PIMPLE: Iteration 3
smoothSolver: Solving for alpha.water, Initial residual = 1.54941e-08, Final residual = 8.2065e-13, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.980694, Final residual = 9.85966e-07, No Iterations 621
smoothSolver: Solving for Uy, Initial residual = 0.845164, Final residual = 9.83645e-07, No Iterations 598
GAMG: Solving for p_rgh, Initial residual = 6.60654e-29, Final residual = 6.60654e-29, No Iterations 0
smoothSolver: Solving for epsilon, Initial residual = 1.13019e-14, Final residual = 2.9905e-16, No Iterations 1
bounding epsilon, min: -8.95102e+115 max: 1.29507e+128 average: 6.2465e+124
smoothSolver: Solving for k, Initial residual = 0.89142, Final residual = 0.0248855, No Iterations 1000
bounding k, min: -1.02608e+81 max: 1.19168e+87 average: 3.5906e+84
ExecutionTime = 474.5 s ClockTime = 474 s

Courant Number mean: 2.93316e-05 max: 543.774
Interface Courant Number mean: 1.269e-06 max: 145.594
deltaT = 3.09075e-51
Time = 3.31999e-06

PIMPLE: Iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 4.16811e-09, Final residual = 9.74614e-13, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.986337, Final residual = 0.0408722, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.99739, Final residual = 0.055142, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 6.89508e-26, Final residual = 6.89508e-26, No Iterations 0
PIMPLE: Iteration 2
smoothSolver: Solving for alpha.water, Initial residual = 2.20691e-07, Final residual = 3.93015e-12, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.984175, Final residual = 0.0316829, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.998156, Final residual = 0.0649702, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 1.65663e-25, Final residual = 1.65663e-25, No Iterations 0
PIMPLE: Iteration 3
smoothSolver: Solving for alpha.water, Initial residual = 0.000573257, Final residual = 1.26673e-09, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.989864, Final residual = 1.63168e-05, No Iterations 1000
smoothSolver: Solving for Uy, Initial residual = 0.992148, Final residual = 8.81591e-06, No Iterations 1000
GAMG: Solving for p_rgh, Initial residual = 1.41589e-24, Final residual = 1.41589e-24, No Iterations 0
smoothSolver: Solving for epsilon, Initial residual = 1.9448e-08, Final residual = 2.42767e-09, No Iterations 1
bounding epsilon, min: -4.72803e+115 max: 1.29219e+128 average: 6.24286e+124
smoothSolver: Solving for k, Initial residual = 0.90641, Final residual = 0.0565075, No Iterations 1000
bounding k, min: -1.87534e+83 max: 2.70993e+92 average: 6.16761e+88
ExecutionTime = 475.36 s ClockTime = 475 s

Courant Number mean: 0.00342074 max: 57888.8
Interface Courant Number mean: 0.00157151 max: 56837.4
deltaT = 1.08757e-55
Time = 3.31999e-06

PIMPLE: Iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 2.46886e-07, Final residual = 5.77144e-13, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.989243, Final residual = 0.0532337, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.97756, Final residual = 0.0346156, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 3.312e-24, Final residual = 3.312e-24, No Iterations 0
PIMPLE: Iteration 2
smoothSolver: Solving for alpha.water, Initial residual = 5.29372e-08, Final residual = 7.05084e-14, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.972908, Final residual = 0.0578711, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.953992, Final residual = 0.0306631, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 1.10713e-22, Final residual = 1.10713e-22, No Iterations 0
PIMPLE: Iteration 3
smoothSolver: Solving for alpha.water, Initial residual = 3.09189e-08, Final residual = 5.15276e-12, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.993079, Final residual = 1.65713e-05, No Iterations 1000
smoothSolver: Solving for Uy, Initial residual = 0.925153, Final residual = 0.000585797, No Iterations 1000
GAMG: Solving for p_rgh, Initial residual = 1.49259e-22, Final residual = 1.49259e-22, No Iterations 0
smoothSolver: Solving for epsilon, Initial residual = 6.01779e-07, Final residual = 2.68609e-13, No Iterations 1
bounding epsilon, min: -3.35361e+120 max: 9.43312e+131 average: 2.17325e+128
smoothSolver: Solving for k, Initial residual = 0.83036, Final residual = 0.0011249, No Iterations 1000
bounding k, min: -3.71259e+83 max: 2.37214e+102 average: 3.15691e+98
ExecutionTime = 476.21 s ClockTime = 476 s

Courant Number mean: 4.45473e-05 max: 266.226
Interface Courant Number mean: 2.69846e-07 max: 4.90987
deltaT = 2.04258e-57
Time = 3.31999e-06

PIMPLE: Iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 2.05394e-09, Final residual = 2.378e-16, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0624997, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0624997, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 5.95873e-21, Final residual = 5.95873e-21, No Iterations 0
PIMPLE: Iteration 2
smoothSolver: Solving for alpha.water, Initial residual = 3.14317e-09, Final residual = 4.02205e-13, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.98147, Final residual = 0.0315776, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.999539, Final residual = 0.031261, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 1.29708e-20, Final residual = 1.29708e-20, No Iterations 0
PIMPLE: Iteration 3
smoothSolver: Solving for alpha.water, Initial residual = 1.29734e-08, Final residual = 2.68178e-12, No Iterations 1
Phase-1 volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.963473 Min(alpha.water) = -0.113971 Max(alpha.water) = 1.9541
smoothSolver: Solving for Ux, Initial residual = 0.919007, Final residual = 9.24707e-07, No Iterations 12
smoothSolver: Solving for Uy, Initial residual = 0.70392, Final residual = 6.86551e-07, No Iterations 11
GAMG: Solving for p_rgh, Initial residual = 2.05474e-19, Final residual = 2.05474e-19, No Iterations 0
#0 Foam::error:

rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::tmp<Foam:

imensionedField<double, Foam::volMesh> > Foam:

perator*<Foam::volMesh>(Foam::tmp<Foam:

im ensionedField<double, Foam::volMesh> > const&, Foam:

imensionedField<double, Foam::volMesh> const&) in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/interPhaseChangeFoam"
#5 Foam::RASModels::kEpsilon<Foam::IncompressibleMome ntumTransportModel<Foam::transportModel> >::correct() at ??:?
#6 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/interPhaseChangeFoam"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/interPhaseChangeFoam"
Floating point exception (core dumped)

night-hawk · December 9, 2021, 20:03

There is probably issue with boundary condition. Can you upload your case folder.

dakaezekoye · December 9, 2021, 21:34

Quote:

Originally Posted by night-hawk

There is probably issue with boundary condition. Can you upload your case folder.

I have uploaded the zero and system folders in zipped format, kindly find attached

night-hawk · December 9, 2021, 22:13

I will check. I will let you know if I am able to made any progress.

dakaezekoye · December 10, 2021, 01:08

Quote:

Originally Posted by night-hawk

I will check. I will let you know if I am able to made any progress.

Ok, Thank you

dakaezekoye · December 10, 2021, 09:53

Night-hawk I removed the polymesh folder just like you suggested and the file size dropped significantly. I have uploaded the 3D case folder.

Secondly, as you can see from the screenshot attached, the 2D case is currently calculating thanks to your correction but so far I have not seen a single calculation folder been generated by OpenFOAM in my case folder. It seems abnormal. Please what can I do about this?

night-hawk · December 10, 2021, 20:18

Your write interval seems don't match with deltaT.

December 7, 2021, 07:25	multiphaseEulerFoam FOAM FATAL IO ERROR	#1
qutadah.r Senior Member qutadah Join Date: Jun 2021 Location: USA Posts: 101 Rep Power: 5	Hi everyone, I am having the following error when running a two-phase flow in OpenFOAM. The original case works perfectly and doesnt have any issues, I come across this issue when i change some values of alpha or U in the boundary conditions... * I am using collated filehandler and parallel processing using two nodes, 48 processors each. Code: --> FOAM FATAL IO ERROR: Expected a ')' while reading binaryBlock, found on line 21 the word 'e▒L▒3' file: /rwthfs/rz/cluster/hpcwork/jg235250/SetupsMarkus/EE2F-MUSIG/06_3D_Pumpe_M1_v6_Qtot696_IGVF_5_SST_OF9_EE2F_MUSIG/0.288173000000247/U.water at line 21. From function Foam::Istream &Foam::Istream::readEnd(const char *) in file db/IOstreams/IOstreams/Istream.C at line 109. Any help would be appreciated! Thanks!

Thread Tools	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[OpenFOAM] ParaView command in Foam-extend-4.1	mitu_94	ParaView	0	March 4, 2021 13:46
Mesquite - Adaptive mesh refinement / coarsening?	philippose	OpenFOAM Running, Solving & CFD	94	January 27, 2016 09:40
Compiling dynamicTopoFvMesh for OpenFOAM 2.1.x	Saxwax	OpenFOAM Installation	25	November 29, 2013 05:34
checking the system setup and Qt version	vivek070176	OpenFOAM Installation	22	June 1, 2010 12:34
Version 15 on Mac OS X	gschaider	OpenFOAM Installation	113	December 2, 2009 10:23

December 7, 2021, 13:25		#2
qutadah.r Senior Member qutadah Join Date: Jun 2021 Location: USA Posts: 101 Rep Power: 5	so i dont know what caused this still, i was using nano editor to edit the binary boundary files... I just copied the original case again and used vi editor and it worked.... Even though i did the same and edited the binary files with nano and got this error ..... If someone knows if nano editor messes up the binary files somehow please let me kindly know, thanks

December 9, 2021, 20:03		#7
night-hawk Member Ashutosh Join Date: Jul 2021 Location: India Posts: 76 Rep Power: 5	There is probably issue with boundary condition. Can you upload your case folder.

December 9, 2021, 22:13		#9
night-hawk Member Ashutosh Join Date: Jul 2021 Location: India Posts: 76 Rep Power: 5	I will check. I will let you know if I am able to made any progress.

December 10, 2021, 20:18		#12
night-hawk Member Ashutosh Join Date: Jul 2021 Location: India Posts: 76 Rep Power: 5	Your write interval seems don't match with deltaT.