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/)
-   -   Two strange problems! (https://www.cfd-online.com/Forums/openfoam-solving/218191-two-strange-problems.html)

NewKid June 11, 2019 15:04

Two strange problems!
 
Hi, Foamers, I encountered two strange problems while solving two cases.
Problem #1: After solving one case from 0 to 2e-4 sec, I had to continue to solve the case from 2e-4 to 1e-3 sec. Therefore, I changed the controlDict in 'system' folder as follows:
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application sprayFoam;

startFrom startTime;

startTime 0.0002;//0.000276;//0;//3.4e-5;//0.00324;//0.00027;

stopAt endTime;

endTime 1e-3;//0.000279;//0.1;//630=1.5873e-04s

deltaT 1e-9;

writeControl timeStep;//adjustableRunTime;

writeInterval 1000;//1000;//1e-04;

purgeWrite 0;

writeFormat binary;

writePrecision 6;

writeCompression uncompressed;
Anyway, once I started the sprayFoam, it gave me the error messages as follows:
Starting time loop

Courant Number mean: 0.000916657 max: 0.0340124
Time = 0.000200001


Solving 3-D cloud sprayCloud


--> FOAM FATAL ERROR:
Integrating outside table 0..0.0002 not implemented. t1:0.000199967 t2:0.000200001

From function linearInterpolationWeights::integrationWeights(..)
in file interpolations/interpolationWeights/linearInterpolationWeights/linearInterpolationWeights.C at line 180.

FOAM exiting

I have no ideas why it happens, since I did it early to continue my computation from an existing time.

Problem #2: I constructed another case which shared the same mesh, boundary conditions and almost everything as the case in Problem #1, the only difference between the two cases is that the property 'coupled' was set true in 'sprayCloudProperties' file in 'constant' folder while it was set as false in the case in problem #1. Once the deltaT was set the same value as that in the first case, it would not take too long before 'Floating point exception' occurred. As the deltaT decreased, the time before error may be a little longer.

Courant Number mean: 1.83702e-06 max: 0.11739
Time = 1.739e-06


Solving 3-D cloud sprayCloud
Cloud: sprayCloud
Current number of parcels = 17
Current mass in system = 4.97306e-05
Linear momentum = (-3.50334e-05 -0.000157153 0.0042042)
|Linear momentum| = 0.00420728
Linear kinetic energy = 0.179701
model1:
number of parcels added = 17
mass introduced = 4.97309e-05
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Temperature min/max = 287.495, 424.183
Mass transfer phase change = 2.30384e-10
D10, D32, Dmax (mu) = 65.1829, 91.4678, 152.159
Liquid penetration 95% mass (m) = 0.000139409

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
smoothSolver: Solving for Ux, Initial residual = 2.90473e-05, Final residual = 3.64878e-11, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.00010294, Final residual = 2.61305e-09, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 3.49818e-05, Final residual = 6.6992e-11, No Iterations 1
DILUPBiCG: Solving for C8H18, Initial residual = 2.94076e-06, Final residual = 1.79667e-10, No Iterations 1
DILUPBiCG: Solving for O2, Initial residual = 2.94074e-06, Final residual = 7.62067e-11, No Iterations 1
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for h, Initial residual = 2.95332e-05, Final residual = 1.60485e-10, No Iterations 1
T gas min/max = 216.37, 1600
GAMG: Solving for p, Initial residual = 2.82251e-05, Final residual = 2.26802e-18, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.45359e-17, global = 1.87181e-18, cumulative = 4.86682e-09
rho max/min : 9.77452 0
GAMG: Solving for p, Initial residual = 2.2221e-09, Final residual = 2.2221e-09, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.79589e-10, global = 4.7143e-13, cumulative = 4.86729e-09
rho max/min : 9.77452 0
ExecutionTime = 530.29 s ClockTime = 535 s

Courant Number mean: 1.83696e-06 max: 0.0965713
Time = 1.7392e-06


Solving 3-D cloud sprayCloud
Cloud: sprayCloud
Current number of parcels = 17
Current mass in system = 4.97306e-05
Linear momentum = (-3.50334e-05 -0.000157153 0.0042042)
|Linear momentum| = 0.00420728
Linear kinetic energy = 0.179701
model1:
number of parcels added = 17
mass introduced = 4.97309e-05
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Temperature min/max = 287.495, 424.183
Mass transfer phase change = 2.30411e-10
D10, D32, Dmax (mu) = 65.1829, 91.4677, 152.159
Liquid penetration 95% mass (m) = 0.000139426

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
smoothSolver: Solving for Ux, Initial residual = 2.9045e-05, Final residual = 3.65374e-11, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.000105027, Final residual = 2.67707e-09, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 3.49899e-05, Final residual = 6.88286e-11, No Iterations 1
DILUPBiCG: Solving for C8H18, Initial residual = 2.95174e-06, Final residual = 1.55185e-10, No Iterations 1
DILUPBiCG: Solving for O2, Initial residual = 2.95171e-06, Final residual = 6.59306e-11, No Iterations 1
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
yhrun: error: cn10680: task 39: Floating point exception
yhrun: First task exited 60s ago
yhrun: tasks 0-38,40-179: running
yhrun: task 39: exited abnormally
yhrun: Terminating job step 11564550.0
slurmd[cn10677]: *** STEP 11564550.0 KILLED AT 2019-06-12T02:26:41 WITH SIGNAL 9 ***
yhrun: Job step aborted: Waiting up to 2 seconds for job step to finish.
slurmd[cn10677]: *** STEP 11564550.0 KILLED AT 2019-06-12T02:26:41 WITH SIGNAL 9 ***
yhrun: error: cn11673: tasks 144-155: Killed

Courant number is not big before the error.
Any clues? Thanks in advance!

NewKid June 11, 2019 15:52

Problem #1 has been solved!
I forget to change the 'flowRateProfile table' in the 'sprayCloudProperties' file in 'constant' folder.

Anything about Problem #2 is welcomed!


All times are GMT -4. The time now is 22:08.