[Aditya Patil]

Dear All,

Kindly find attached the case and let me know if someone can help me out with this. Its a Shell and tube heat exchanger, i am using RNGkepsilon model.
I have generated mesh in icem and imported it in OpenFOAM (which took two months for me to figure it out how to import internal walls).

Now when i run the solver, i get the following error. My Master's Thesis depends on this.

It would be great help if someone can help me out with this.

apatil996@gmail.com

[Bloerb]

There is no case attached?

[Aditya Patil]

I am not able to upload the file if you can share your mail id i can share the files with you

[Aditya Patil]

Solving for fluid region water
smoothSolver: Solving for Ux, Initial residual = 0.9790392, Final residual = 0.005736963, No Iterations 5
smoothSolver: Solving for Uy, Initial residual = 0.9308278, Final residual = 0.005755914, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.9411201, Final residual = 0.005565112, No Iterations 4
smoothSolver: Solving for h, Initial residual = 0.5149063, Final residual = 5.30946e-006, No Iterations 1
Min/max T:300 330
GAMG: Solving for p_rgh, Initial residual = 0.05148115, Final residual = 0.0004485215, No Iterations 7
GAMG: Solving for p_rgh, Initial residual = 0.1919115, Final residual = 0.001488724, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 0.019925, Final residual = 0.0001984095, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.00315037, Final residual = 2.630818e-005, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.0005758737, Final residual = 4.178398e-006, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.0001134084, Final residual = 8.907421e-007, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 2.646995e-005, Final residual = 1.930104e-007, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 6.741278e-006, Final residual = 6.422404e-008, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 1.953217e-006, Final residual = 3.998076e-008, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 5.720107e-007, Final residual = 4.862596e-008, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1.797649e-007, Final residual = 1.529353e-008, No Iterations 1
time step continuity errors : sum local = 0.7081896, global = -0.004061983, cumulative = -0.004190197
Min/max rho:998.2 998.2
DILUPBiCG: Solving for epsilon, Initial residual = 0.966629, Final residual = 0.005164755, No Iterations 6
DILUPBiCG: Solving for k, Initial residual = 0.9137956, Final residual = 0.003374522, No Iterations 3
bounding k, min: -0.003801141 max: 3.597947 average: 0.1961555

Solving for solid region solid
DICPCG: Solving for h, Initial residual = 0.9999976, Final residual = 1.355595e-005, No Iterations 4
DICPCG: Solving for h, Initial residual = 0.1632307, Final residual = 3.701048e-005, No Iterations 3
DICPCG: Solving for h, Initial residual = 0.07990981, Final residual = 1.714419e-005, No Iterations 3
DICPCG: Solving for h, Initial residual = 0.05253247, Final residual = 1.017037e-005, No Iterations 3
Min/max T:300 300.0007
ExecutionTime = 77.766 s ClockTime = 78 s

Time = 3


Solving for fluid region gas
smoothSolver: Solving for Ux, Initial residual = 0.4012478, Final residual = 0.002204488, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.4779169, Final residual = 0.00178186, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.4704772, Final residual = 0.001735573, No Iterations 3
smoothSolver: Solving for h, Initial residual = 0.07085427, Final residual = 0.0001470701, No Iterations 1
Min/max T:300 450
GAMG: Solving for p_rgh, Initial residual = 0.4096719, Final residual = 0.002440951, No Iterations 7
GAMG: Solving for p_rgh, Initial residual = 0.1968471, Final residual = 0.000888509, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.01984266, Final residual = 0.0001255817, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.00363527, Final residual = 1.506567e-005, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.0006360669, Final residual = 5.624142e-006, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.0001621235, Final residual = 1.24138e-006, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 6.003564e-005, Final residual = 3.667581e-007, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 2.887354e-005, Final residual = 2.00293e-007, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 1.847279e-005, Final residual = 1.114627e-007, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 1.462804e-005, Final residual = 9.228124e-008, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 1.281861e-005, Final residual = 8.007105e-008, No Iterations 4
time step continuity errors : sum local = 0.0005627474, global = -5.62112e-006, cumulative = -0.004195818
Min/max rho:0.55 0.55
DILUPBiCG: Solving for epsilon, Initial residual = 0.2496675, Final residual = 0.002029951, No Iterations 14
DILUPBiCG: Solving for k, Initial residual = 0.694531, Final residual = 0.006792958, No Iterations 58

Solving for fluid region water
smoothSolver: Solving for Ux, Initial residual = 0.328951, Final residual = 0.001382712, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.6077026, Final residual = 0.002674739, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.5900856, Final residual = 0.00245067, No Iterations 3
smoothSolver: Solving for h, Initial residual = 0.2817345, Final residual = 9.76883e-006, No Iterations 1
Min/max T:300 330
GAMG: Solving for p_rgh, Initial residual = 0.8069214, Final residual = 0.003880355, No Iterations 10
GAMG: Solving for p_rgh, Initial residual = 0.2790827, Final residual = 0.00278334, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 0.02385825, Final residual = 0.000143512, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.003995125, Final residual = 2.444764e-005, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.0008686896, Final residual = 6.135138e-006, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.0001456225, Final residual = 1.005563e-006, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 4.349722e-005, Final residual = 2.500999e-007, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 9.529139e-006, Final residual = 5.844189e-008, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 2.743393e-006, Final residual = 8.732328e-008, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 7.584416e-007, Final residual = 6.50084e-008, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 2.198307e-007, Final residual = 2.503013e-008, No Iterations 1
time step continuity errors : sum local = 0.6008994, global = 0.001049441, cumulative = -0.003146376
Min/max rho:998.2 998.2
DILUPBiCG: Solving for epsilon, Initial residual = 0.6036671, Final residual = 0.001586308, No Iterations 2
bounding epsilon, min: -1639.298 max: 123706.8 average: 220.2834
DILUPBiCG: Solving for k, Initial residual = 0.5545638, Final residual = 0.003996428, No Iterations 25
bounding k, min: -0.2369039 max: 96.80032 average: 0.5651651

Solving for solid region solid
DICPCG: Solving for h, Initial residual = 0.3641707, Final residual = 7.723615e-005, No Iterations 3
DICPCG: Solving for h, Initial residual = 0.132716, Final residual = 2.676192e-005, No Iterations 3
DICPCG: Solving for h, Initial residual = 0.07120789, Final residual = 1.422047e-005, No Iterations 3
DICPCG: Solving for h, Initial residual = 0.04786442, Final residual = 8.93938e-006, No Iterations 3
Min/max T:300 300.0108
ExecutionTime = 100.723 s ClockTime = 101 s

Time = 4


Solving for fluid region gas
smoothSolver: Solving for Ux, Initial residual = 0.997895, Final residual = 0.005140062, No Iterations 5
smoothSolver: Solving for Uy, Initial residual = 0.9951449, Final residual = 0.003863077, No Iterations 5
smoothSolver: Solving for Uz, Initial residual = 0.9957609, Final residual = 0.004041113, No Iterations 5
smoothSolver: Solving for h, Initial residual = 0.9998887, Final residual = 0.00215796, No Iterations 1
Min/max T:300 450
GAMG: Solving for p_rgh, Initial residual = 0.1707758, Final residual = 0.000691169, No Iterations 8
GAMG: Solving for p_rgh, Initial residual = 0.2205068, Final residual = 0.0007038493, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.02429237, Final residual = 0.0001814, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.00381364, Final residual = 2.121552e-005, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.0008595016, Final residual = 7.324344e-006, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.0001714744, Final residual = 9.041996e-007, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 5.990239e-005, Final residual = 3.708448e-007, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 1.419537e-005, Final residual = 1.132744e-007, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 5.840415e-006, Final residual = 5.827181e-008, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 2.05174e-006, Final residual = 6.060779e-008, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 9.098214e-007, Final residual = 2.987628e-008, No Iterations 2
time step continuity errors : sum local = 0.142911, global = 0.01224158, cumulative = 0.009095205
Min/max rho:0.55 0.55
DILUPBiCG: Solving for epsilon, Initial residual = 0.963659, Final residual = 0.009422052, No Iterations 11
bounding epsilon, min: -4.29482e+007 max: 6.201085e+014 average: 3.407446e+012
DILUPBiCG: Solving for k, Initial residual = 0.8746861, Final residual = 0.001735249, No Iterations 1
bounding k, min: -130.283 max: 1357584 average: 13843.08

Solving for fluid region water
smoothSolver: Solving for Ux, Initial residual = 0.4258334, Final residual = 0.002836647, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.5487691, Final residual = 0.003929843, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.5808286, Final residual = 0.003483625, No Iterations 3
smoothSolver: Solving for h, Initial residual = 0.003828602, Final residual = 1.012619e-005, No Iterations 1
Min/max T:300 330
GAMG: Solving for p_rgh, Initial residual = 0.841575, Final residual = 0.006002557, No Iterations 10
GAMG: Solving for p_rgh, Initial residual = 0.9909519, Final residual = 0.008380402, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.1576031, Final residual = 0.001345169, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.02590365, Final residual = 0.0001482172, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.00611356, Final residual = 3.423935e-005, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.001096689, Final residual = 5.206981e-006, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.0002935457, Final residual = 1.731125e-006, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 6.052588e-005, Final residual = 2.527115e-007, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 1.618924e-005, Final residual = 1.075399e-007, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 3.839643e-006, Final residual = 7.358797e-008, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 1.003823e-006, Final residual = 7.899073e-008, No Iterations 1
time step continuity errors : sum local = 4.670205e+010, global = 3.188148e+008, cumulative = 3.188148e+008
Min/max rho:998.2 998.2
DILUPBiCG: Solving for epsilon, Initial residual = 1, Final residual = 4.047027, No Iterations 101
bounding epsilon, min: -7.46515e+017 max: 4.762213e+017 average: 3.265049e+015
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.007218006, No Iterations 47
bounding k, min: -5.993124e+013 max: 2.393794e+020 average: 3.747147e+015

Solving for solid region solid
DICPCG: Solving for h, Initial residual = 0.9855608, Final residual = 2.461337e-005, No Iterations 4
DICPCG: Solving for h, Initial residual = 0.153145, Final residual = 5.186669e-005, No Iterations 3
DICPCG: Solving for h, Initial residual = 0.08033832, Final residual = 2.406222e-005, No Iterations 3
DICPCG: Solving for h, Initial residual = 0.05559123, Final residual = 1.43891e-005, No Iterations 3
Min/max T:300 306.5061
ExecutionTime = 124.02 s ClockTime = 124 s

Time = 5


Solving for fluid region gas
smoothSolver: Solving for Ux, Initial residual = 0.2762095, Final residual = 0.001942218, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.5057565, Final residual = 0.002582315, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.5033928, Final residual = 0.002791919, No Iterations 3
smoothSolver: Solving for h, Initial residual = 0.9881829, Final residual = 0.002039114, No Iterations 1
Min/max T:300 450
GAMG: Solving for p_rgh, Initial residual = 0.9601526, Final residual = 0.008333369, No Iterations 8
GAMG: Solving for p_rgh, Initial residual = 0.2487248, Final residual = 0.0009569042, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.0236621, Final residual = 0.0002168057, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.0045142, Final residual = 2.598943e-005, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 0.001340493, Final residual = 1.241336e-005, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.0002866616, Final residual = 2.806589e-006, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 9.378252e-005, Final residual = 8.509813e-007, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 3.749074e-005, Final residual = 2.650459e-007, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 1.567895e-005, Final residual = 1.513883e-007, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 9.221193e-006, Final residual = 9.212312e-008, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 5.732154e-006, Final residual = 5.995503e-008, No Iterations 4
time step continuity errors : sum local = 0.2473274, global = 0.004257796, cumulative = 3.188148e+008
Min/max rho:0.55 0.55
DILUPBiCG: Solving for epsilon, Initial residual = 0.9998053, Final residual = 1.254621e-007, No Iterations 1
bounding epsilon, min: -1.120723e+013 max: 1.542508e+013 average: 1.128546e+010
DILUPBiCG: Solving for k, Initial residual = 0.6029677, Final residual = 0.004540687, No Iterations 26
bounding k, min: -30040.46 max: 3361104 average: 47591.83

Solving for fluid region water
smoothSolver: Solving for Ux, Initial residual = 0.6442359, Final residual = 0.006100699, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.6184839, Final residual = 0.006119722, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.696461, Final residual = 0.006920729, No Iterations 3
smoothSolver: Solving for h, Initial residual = 0.002719333, Final residual = 9.114461e-006, No Iterations 1
Min/max T:300 450
GAMG: Solving for p_rgh, Initial residual = 0.9199075, Final residual = 0.2172644, No Iterations 1000
GAMG: Solving for p_rgh, Initial residual = 0.9384384, Final residual = 4.203924, No Iterations 1000
GAMG: Solving for p_rgh, Initial residual = 0.9999959, Final residual = 0.00531034, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.4049094, Final residual = 0.003023114, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.138636, Final residual = 0.001346446, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 0.04641499, Final residual = 0.0003863779, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.01776378, Final residual = 0.0001020957, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.004052399, Final residual = 2.448172e-005, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 0.002342693, Final residual = 2.173853e-005, No Iterations 3
GAMG: Solving for p_rgh, Initial residual = 0.0008084469, Final residual = 5.883629e-006, No Iterations 14
GAMG: Solving for p_rgh, Initial residual = 0.0006418857, Final residual = 4.189607e-006, No Iterations 10
time step continuity errors : sum local = 1.712006e+042, global = 6.981974e+036, cumulative = 6.981974e+036
Min/max rho:998.2 998.2
DILUPBiCG: Solving for epsilon, Initial residual = 1, Final residual = 51.83246, No Iterations 101
bounding epsilon, min: -6.045852e+061 max: 1.02814e+062 average: 4.12244e+058
DILUPBiCG: Solving for k, Initial residual = 0.9817559, Final residual = 0.008315693, No Iterations 8
bounding k, min: -1.657296e+059 max: 3.275837e+067 average: 5.043551e+062

Solving for solid region solid
DICPCG: Solving for h, Initial residual = 0.9603932, Final residual = 5.594731e-006, No Iterations 4
DICPCG: Solving for h, Initial residual = 0.1651377, Final residual = 2.231357e-005, No Iterations 3
DICPCG: Solving for h, Initial residual = 0.07915659, Final residual = 1.065266e-005, No Iterations 3
DICPCG: Solving for h, Initial residual = 0.04969039, Final residual = 6.360309e-006, No Iterations 3
Min/max T:300 450
ExecutionTime = 180.901 s ClockTime = 181 s

Time = 6


Solving for fluid region gas
smoothSolver: Solving for Ux, Initial residual = 0.6578743, Final residual = 0.002162361, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.5158366, Final residual = 0.00193017, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.5809276, Final residual = 0.003274813, No Iterations 3
smoothSolver: Solving for h, Initial residual = 0.05025332, Final residual = 0.0001729062, No Iterations 1
Min/max T:300 450
GAMG: Solving for p_rgh, Initial residual = 0.815518, Final residual = 0.004156763, No Iterations 12
GAMG: Solving for p_rgh, Initial residual = 0.5630934, Final residual = 0.002945261, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 0.1379688, Final residual = 0.000791544, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.03394696, Final residual = 0.0002271798, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.01405888, Final residual = 8.895992e-005, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.009398774, Final residual = 6.151751e-005, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 0.01363227, Final residual = 9.886645e-005, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 0.0266589, Final residual = 0.0002343467, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 0.05440896, Final residual = 0.0005054366, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 0.1066077, Final residual = 0.001021793, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 0.1834895, Final residual = 0.001765707, No Iterations 5
time step continuity errors : sum local = 1.271341e+022, global = 2.479593e+021, cumulative = 6.981974e+036
Min/max rho:0.55 0.55
DILUPBiCG: Solving for epsilon, Initial residual = 1, Final residual = nan, No Iterations 101
DILUPBiCG: Solving for k, Initial residual = nan, Final residual = nan, No Iterations 101

Solving for fluid region water
smoothSolver: Solving for Ux, Initial residual = 7.045395e-008, Final residual = 7.045395e-008, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 3.340629e-008, Final residual = 3.340629e-008, No Iterations 0
smoothSolver: Solving for Uz, Initial residual = 1.967e-008, Final residual = 1.967e-008, No Iterations 0
smoothSolver: Solving for h, Initial residual = 0.002853995, Final residual = 1.045642e-005, No Iterations 1
Min/max T:300 450
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = nan, No Iterations 1000
GAMG: Solving for p_rgh, Initial residual = nan, Final residual = nan, No Iterations 1000


--> FOAM FATAL IO ERROR:
wrong token type - expected Scalar, found on line 0 the word 'nan'

file: C:/OpenFOAM/Aadhi-3.0.x/run/EGRC2_OF_Case/system/water/data.solverPerformance.p_rgh at line 0.

From function operator>>(Istream&, Scalar&)
in file lnInclude/Scalar.C at line 93.

FOAM exiting

[piu58]

Your simulation gains turbulence energy in every step. This may be due to an unfortunately formulation of your problem or due to incorrect boundary conditions.
I recommend to write the result every step and look at it with paraView. The regions which give the problems arise then. After that you have to check your model.

[Bloerb]

Hi,

I have looked at the case. A few quick things. Why the total temperature boundary condition instead of fixed value? Why is the flowRate defined on both inlet and outlet. Keep the inletOutlet Condition but delete the flowRate lines at the outlet. You define k as a fixedValue on both inlet and outlet. The boundary conditions are all mixed up.

And one last thing to consider. This is not a good mesh. remeshing with blockMesh would greatly improve your results!

[Aditya Patil]

Dear All,

I am still stuck with the same problem even though i have made the corrections mentioned in earlier comments, i am not able to get the results.
If someone has solved a similar case please contact me so that i can know where i am going wrong.

Regards,
Aditya Patil