Negative temperature with twoPhaseEulerFoam

tankerp · May 13, 2018, 09:20

Hi !!!

I am currently using the twoPhaseEulerFoam solver to study spray formation with air and water.

I initially simulated a single phase flow with rhoPimpleFoam(since the air flow is compressible). The simulation converged, so I moved onto the two phase simulations.

With twoPhaseEulerFoam, my simulations continuously result in negative temperatures for air. Does anyone have any suggestions for this problem? I have attached the error log below:

fieldAverage fieldAverage1 write:
Calculating averages

Courant Number mean: 9.677663e-05 max: 0.07496172
Max Ur Courant Number = 0.5125739
deltaT = 4.898478e-09
Time = 0.0006120702

PIMPLE: iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1
Constructing momentum equations
smoothSolver: Solving for e.air, Initial residual = 1.196049e-05, Final residual = 5.437223e-10, No Iterations 1
smoothSolver: Solving for e.water, Initial residual = 1.687357e-06, Final residual = 1.802551e-12, No Iterations 1
min T.air 4.285799
min T.water 227.2014
GAMG: Solving for p_rgh, Initial residual = 3.477521e-05, Final residual = 4.661494e-16, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 4.661153e-16, Final residual = 4.661153e-16, No Iterations 0
GAMG: Solving for p_rgh, Initial residual = 2.884329e-05, Final residual = 4.12967e-16, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 4.130639e-16, Final residual = 4.130639e-16, No Iterations 0
PIMPLE: iteration 2
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1
Constructing momentum equations
smoothSolver: Solving for e.air, Initial residual = 3.511325e-06, Final residual = 2.756874e-10, No Iterations 1
smoothSolver: Solving for e.water, Initial residual = 2.847783e-05, Final residual = 1.862502e-11, No Iterations 1
min T.air 0.1135057
min T.water 224.8915
GAMG: Solving for p_rgh, Initial residual = 2.938289e-05, Final residual = 4.642091e-16, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 4.643053e-16, Final residual = 4.643053e-16, No Iterations 0
GAMG: Solving for p_rgh, Initial residual = 2.892439e-05, Final residual = 4.064227e-16, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 4.065217e-16, Final residual = 4.065217e-16, No Iterations 0
smoothSolver: Solving for epsilonm, Initial residual = 0.00254583, Final residual = 2.93189e-08, No Iterations 1
smoothSolver: Solving for km, Initial residual = 0.0004375, Final residual = 6.476975e-09, No Iterations 1
ExecutionTime = 1122.68 s

fieldAverage fieldAverage1 write:
Calculating averages

Courant Number mean: 9.440734e-05 max: 0.07565087
Max Ur Courant Number = 0.5126871
deltaT = 4.777259e-09
Time = 0.000612075

PIMPLE: iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1
Constructing momentum equations
smoothSolver: Solving for e.air, Initial residual = 2.587203e-05, Final residual = 5.203098e-08, No Iterations 1
smoothSolver: Solving for e.water, Initial residual = 1.830379e-06, Final residual = 1.922872e-12, No Iterations 1
min T.air -0.5202558
min T.water 224.9178
GAMG: Solving for p_rgh, Initial residual = 5.205457e-05, Final residual = 4.613364e-16, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 4.608789e-16, Final residual = 4.608789e-16, No Iterations 0
GAMG: Solving for p_rgh, Initial residual = 2.082419e-05, Final residual = 4.038371e-16, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 4.039299e-16, Final residual = 4.039299e-16, No Iterations 0
PIMPLE: iteration 2
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1
Constructing momentum equations
[7] #0 Foam::error:

rintStack(Foam::Ostream&) at ??:?
[7] #1 Foam::sigFpe::sigHandler(int) at ??:?
[7] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[7] #3 ? in "/lib/x86_64-linux-gnu/libm.so.6"
[7] #4 pow in "/lib/x86_64-linux-gnu/libm.so.6"
[7] #5 Foam:

ow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:?
[7] #6 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:

ow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) at ??:?
[7] #7 Foam::dragModels::SchillerNaumann::CdRe() const at ??:?
[7] #8 Foam::dragModel::Ki() const at ??:?
[7] #9 Foam::dragModel::K() const at ??:?
[7] #10 Foam::BlendedInterfacialModel<Foam::dragModel>::K( ) const at ??:?
[7] #11 Foam::twoPhaseSystem::Kd() const at ??:?
[7] #12 ? at ??:?
[7] #13 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[7] #14 ? at ??:?
[carnage:23562] *** Process received signal ***
[carnage:23562] Signal: Floating point exception (8)
[carnage:23562] Signal code: (-6)
[carnage:23562] Failing at address: 0x3e900005c0a
[carnage:23562] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f917d3bc4b0]
[carnage:23562] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f917d3bc428]
[carnage:23562] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f917d3bc4b0]
[carnage:23562] [ 3] /lib/x86_64-linux-gnu/libm.so.6(+0x12601)[0x7f917d979601]
[carnage:23562] [ 4] /lib/x86_64-linux-gnu/libm.so.6(pow+0x14)[0x7f917d988e54]
[carnage:23562] [ 5] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libOpenFOAM.so(_ZN4Foam3powERNS_5FieldIdEERKNS_5UL istIdEERKd+0x38)[0x7f917e7f2418]
[carnage:23562] [ 6] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libcompressibleTurbulenceModels.so(_ZN4Foam3powINS _12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14Geometr icFieldIdT_T0_EEEERKS7_RKNS_11dimensionedIdEE+0x23 7)[0x7f9182a582b7]
[carnage:23562] [ 7] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam10dragModels15SchillerNaumann4CdReEv+0x1e0)[0x7f9181b5d1c0]
[carnage:23562] [ 8] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam9dragModel2KiEv+0x145)[0x7f9181b21455]
[carnage:23562] [ 9] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam9dragModel1KEv+0x2e)[0x7f9181b2123e]
[carnage:23562] [10] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam23Blende dInterfacialModelINS_9dragModelEE1KEv+0x4be)[0x7f9181f69ace]
[carnage:23562] [11] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam14twoPha seSystem2KdEv+0x32)[0x7f9181f45532]
[carnage:23562] [12] twoPhaseEulerFoam[0x43a9fb]
[carnage:23562] [13] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f917d3a7830]
[carnage:23562] [14] twoPhaseEulerFoam[0x44e3b9]
[carnage:23562] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 7 with PID 23562 on node carnage exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

tankerp · May 13, 2018, 09:24

I can provide more information if required!

ed_bar · May 15, 2018, 23:43

You've got a divide by zero in there somewhere. Have you copied across a tutorial case and changed mesh and BC/ICs for your conditions? If so, you will need to determine if the discretisation methods and solution controls are appropriate for your case. How many outer correctors do you have? How's your mesh (hint: checkMesh). If the non-orthogonality is above 0.5, then you will need to throw in some nNonOrthogonalCorrectors.

There are many things that could have gone wrong. Please reply with your case attached, or link to the case on github if possible.

May 13, 2018, 09:20	Negative temperature with twoPhaseEulerFoam	#1
tankerp New Member krishnan Join Date: Feb 2018 Posts: 16 Rep Power: 9	Hi !!! I am currently using the twoPhaseEulerFoam solver to study spray formation with air and water. I initially simulated a single phase flow with rhoPimpleFoam(since the air flow is compressible). The simulation converged, so I moved onto the two phase simulations. With twoPhaseEulerFoam, my simulations continuously result in negative temperatures for air. Does anyone have any suggestions for this problem? I have attached the error log below: fieldAverage fieldAverage1 write: Calculating averages Courant Number mean: 9.677663e-05 max: 0.07496172 Max Ur Courant Number = 0.5125739 deltaT = 4.898478e-09 Time = 0.0006120702 PIMPLE: iteration 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 1.196049e-05, Final residual = 5.437223e-10, No Iterations 1 smoothSolver: Solving for e.water, Initial residual = 1.687357e-06, Final residual = 1.802551e-12, No Iterations 1 min T.air 4.285799 min T.water 227.2014 GAMG: Solving for p_rgh, Initial residual = 3.477521e-05, Final residual = 4.661494e-16, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 4.661153e-16, Final residual = 4.661153e-16, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 2.884329e-05, Final residual = 4.12967e-16, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 4.130639e-16, Final residual = 4.130639e-16, No Iterations 0 PIMPLE: iteration 2 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 3.511325e-06, Final residual = 2.756874e-10, No Iterations 1 smoothSolver: Solving for e.water, Initial residual = 2.847783e-05, Final residual = 1.862502e-11, No Iterations 1 min T.air 0.1135057 min T.water 224.8915 GAMG: Solving for p_rgh, Initial residual = 2.938289e-05, Final residual = 4.642091e-16, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 4.643053e-16, Final residual = 4.643053e-16, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 2.892439e-05, Final residual = 4.064227e-16, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 4.065217e-16, Final residual = 4.065217e-16, No Iterations 0 smoothSolver: Solving for epsilonm, Initial residual = 0.00254583, Final residual = 2.93189e-08, No Iterations 1 smoothSolver: Solving for km, Initial residual = 0.0004375, Final residual = 6.476975e-09, No Iterations 1 ExecutionTime = 1122.68 s fieldAverage fieldAverage1 write: Calculating averages Courant Number mean: 9.440734e-05 max: 0.07565087 Max Ur Courant Number = 0.5126871 deltaT = 4.777259e-09 Time = 0.000612075 PIMPLE: iteration 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 2.587203e-05, Final residual = 5.203098e-08, No Iterations 1 smoothSolver: Solving for e.water, Initial residual = 1.830379e-06, Final residual = 1.922872e-12, No Iterations 1 min T.air -0.5202558 min T.water 224.9178 GAMG: Solving for p_rgh, Initial residual = 5.205457e-05, Final residual = 4.613364e-16, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 4.608789e-16, Final residual = 4.608789e-16, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 2.082419e-05, Final residual = 4.038371e-16, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 4.039299e-16, Final residual = 4.039299e-16, No Iterations 0 PIMPLE: iteration 2 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.9950781 Min(alpha.air) = 0 Max(alpha.air) = 1 Constructing momentum equations [7] #0 Foam::error:rintStack(Foam::Ostream&) at ??:? [7] #1 Foam::sigFpe::sigHandler(int) at ??:? [7] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" [7] #3 ? in "/lib/x86_64-linux-gnu/libm.so.6" [7] #4 pow in "/lib/x86_64-linux-gnu/libm.so.6" [7] #5 Foam:ow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:? [7] #6 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:ow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) at ??:? [7] #7 Foam::dragModels::SchillerNaumann::CdRe() const at ??:? [7] #8 Foam::dragModel::Ki() const at ??:? [7] #9 Foam::dragModel::K() const at ??:? [7] #10 Foam::BlendedInterfacialModel<Foam::dragModel>::K( ) const at ??:? [7] #11 Foam::twoPhaseSystem::Kd() const at ??:? [7] #12 ? at ??:? [7] #13 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" [7] #14 ? at ??:? [carnage:23562] * Process received signal * [carnage:23562] Signal: Floating point exception (8) [carnage:23562] Signal code: (-6) [carnage:23562] Failing at address: 0x3e900005c0a [carnage:23562] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f917d3bc4b0] [carnage:23562] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f917d3bc428] [carnage:23562] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f917d3bc4b0] [carnage:23562] [ 3] /lib/x86_64-linux-gnu/libm.so.6(+0x12601)[0x7f917d979601] [carnage:23562] [ 4] /lib/x86_64-linux-gnu/libm.so.6(pow+0x14)[0x7f917d988e54] [carnage:23562] [ 5] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libOpenFOAM.so(_ZN4Foam3powERNS_5FieldIdEERKNS_5UL istIdEERKd+0x38)[0x7f917e7f2418] [carnage:23562] [ 6] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libcompressibleTurbulenceModels.so(_ZN4Foam3powINS _12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14Geometr icFieldIdT_T0_EEEERKS7_RKNS_11dimensionedIdEE+0x23 7)[0x7f9182a582b7] [carnage:23562] [ 7] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam10dragModels15SchillerNaumann4CdReEv+0x1e0)[0x7f9181b5d1c0] [carnage:23562] [ 8] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam9dragModel2KiEv+0x145)[0x7f9181b21455] [carnage:23562] [ 9] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam9dragModel1KEv+0x2e)[0x7f9181b2123e] [carnage:23562] [10] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam23Blende dInterfacialModelINS_9dragModelEE1KEv+0x4be)[0x7f9181f69ace] [carnage:23562] [11] /home/rhoheth/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt64Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam14twoPha seSystem2KdEv+0x32)[0x7f9181f45532] [carnage:23562] [12] twoPhaseEulerFoam[0x43a9fb] [carnage:23562] [13] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f917d3a7830] [carnage:23562] [14] twoPhaseEulerFoam[0x44e3b9] [carnage:23562] * End of error message * -------------------------------------------------------------------------- mpirun noticed that process rank 7 with PID 23562 on node carnage exited on signal 8 (Floating point exception). --------------------------------------------------------------------------

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
whats the cause of error?	immortality	OpenFOAM Running, Solving & CFD	13	March 24, 2021 07:15
Error Prevention: Negative Temperature Throws Fatal Error in Fluid Sim	ben_	OpenFOAM Running, Solving & CFD	1	August 28, 2017 12:56
is internalField(U) equivalent to zeroGradient?	immortality	OpenFOAM Running, Solving & CFD	7	March 29, 2013 01:27
negative fluid temperature	JeanPierre	FLOW-3D	2	January 3, 2012 16:24
Negative temperature	Julie	Siemens	7	September 13, 2004 08:39

May 13, 2018, 09:24		#2
tankerp New Member krishnan Join Date: Feb 2018 Posts: 16 Rep Power: 9	I can provide more information if required!

May 15, 2018, 23:43		#3
ed_bar New Member Ed Barry Join Date: May 2018 Posts: 4 Rep Power: 9	You've got a divide by zero in there somewhere. Have you copied across a tutorial case and changed mesh and BC/ICs for your conditions? If so, you will need to determine if the discretisation methods and solution controls are appropriate for your case. How many outer correctors do you have? How's your mesh (hint: checkMesh). If the non-orthogonality is above 0.5, then you will need to throw in some nNonOrthogonalCorrectors. There are many things that could have gone wrong. Please reply with your case attached, or link to the case on github if possible.