jrrygg |
November 11, 2012 11:46 |
interFoam case blows up
Hi,
I have made a 3D-case with a mesh imported from Ansys. I want water to enter through a cylinder ("jet") and enter the main domain. I have specified an alpha 1 in negative x-direction, and tried with velocities ranging from 0.01 to 40 m/s. I have also experimented with different options for fvSchemes, fvSolution, max courant-numbers and timesteps with no luck.
The error message I get is:
Code:
Create time
Create mesh for time = 0
PIMPLE: Operating solver in PISO mode
Reading field p_rgh
Reading field alpha1
Reading field U
Reading/calculating face flux field phi
Reading transportProperties
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Reading g
Calculating field g.h
time step continuity errors : sum local = 4.67867e-10, global = -4.67867e-10, cumulative = -4.67867e-10
DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 6.66652e-11, No Iterations 111
time step continuity errors : sum local = 3.11808e-20, global = 5.16e-22, cumulative = -4.67867e-10
Courant Number mean: 2.33626e-08 max: 3.83881e-05
Starting time loop
Courant Number mean: 2.33626e-08 max: 3.83881e-05
Interface Courant Number mean: 0 max: 0
Time = 0.0001
MULES: Solving for alpha1
Phase-1 volume fraction = 2.33933e-10 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 4.67867e-10 Min(alpha1) = 0 Max(alpha1) = 1
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0429735, No Iterations 47
time step continuity errors : sum local = 6.12373e-06, global = -1.9493e-07, cumulative = -1.95398e-07
DICPCG: Solving for p_rgh, Initial residual = 0.0589919, Final residual = 0.00254961, No Iterations 8
time step continuity errors : sum local = 0.000113259, global = 1.82915e-06, cumulative = 1.63376e-06
DICPCG: Solving for p_rgh, Initial residual = 0.0149289, Final residual = 9.11547e-08, No Iterations 73
time step continuity errors : sum local = 4.31244e-09, global = 3.18223e-10, cumulative = 1.63407e-06
ExecutionTime = 0.21 s ClockTime = 1 s
Courant Number mean: 0.00757338 max: 12.278
Interface Courant Number mean: 0 max: 0
Time = 0.0002
MULES: Solving for alpha1
Phase-1 volume fraction = 7.05874e-05 Min(alpha1) = -2.93612e-24 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.000141174 Min(alpha1) = -3.28782e-25 Max(alpha1) = 1
DICPCG: Solving for p_rgh, Initial residual = 0.0166851, Final residual = 0.000693113, No Iterations 59
time step continuity errors : sum local = 3.47964e-05, global = 3.5545e-06, cumulative = 5.18857e-06
DICPCG: Solving for p_rgh, Initial residual = 0.113737, Final residual = 0.00554455, No Iterations 22
time step continuity errors : sum local = 0.000272306, global = 6.33387e-06, cumulative = 1.15224e-05
DICPCG: Solving for p_rgh, Initial residual = 0.0730375, Final residual = 9.70425e-08, No Iterations 115
time step continuity errors : sum local = 5.18939e-09, global = -2.03538e-10, cumulative = 1.15222e-05
ExecutionTime = 0.37 s ClockTime = 1 s
Courant Number mean: 0.0331767 max: 3081.05
Interface Courant Number mean: 7.87481e-07 max: 0.00145705
Time = 0.0003
MULES: Solving for alpha1
Phase-1 volume fraction = 0.000141468 Min(alpha1) = -1.13757e-23 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.000141761 Min(alpha1) = -2.82471e-24 Max(alpha1) = 1
DICPCG: Solving for p_rgh, Initial residual = 0.341486, Final residual = 0.014237, No Iterations 71
time step continuity errors : sum local = 0.00127742, global = -0.000206445, cumulative = -0.000194923
DICPCG: Solving for p_rgh, Initial residual = 0.0948778, Final residual = 0.00473849, No Iterations 23
time step continuity errors : sum local = 0.00587709, global = -0.000116325, cumulative = -0.000311248
DICPCG: Solving for p_rgh, Initial residual = 0.261779, Final residual = 9.31408e-08, No Iterations 160
time step continuity errors : sum local = 1.36897e-07, global = -1.75955e-09, cumulative = -0.00031125
ExecutionTime = 0.54 s ClockTime = 1 s
Courant Number mean: 0.786892 max: 104531
Interface Courant Number mean: 7.91815e-07 max: 0.00161964
Time = 0.0004
MULES: Solving for alpha1
Phase-1 volume fraction = 0.000142055 Min(alpha1) = -2.23505e-22 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.000142348 Min(alpha1) = -3.30942e-22 Max(alpha1) = 1
DICPCG: Solving for p_rgh, Initial residual = 0.721516, Final residual = 0.15886, No Iterations 1001
time step continuity errors : sum local = 1.20394, global = -0.0632404, cumulative = -0.0635517
DICPCG: Solving for p_rgh, Initial residual = 0.0976159, Final residual = 0.0607169, No Iterations 1001
time step continuity errors : sum local = 330.267, global = -19.6465, cumulative = -19.71
DICPCG: Solving for p_rgh, Initial residual = 0.612786, Final residual = 5.75419, No Iterations 1001
time step continuity errors : sum local = 1.06021e+06, global = 110793, cumulative = 110773
ExecutionTime = 1.35 s ClockTime = 2 s
Courant Number mean: 2.9712e+07 max: 1.21991e+12
Interface Courant Number mean: 89.6278 max: 214939
Time = 0.0005
MULES: Solving for alpha1
Phase-1 volume fraction = 0.00014264 Min(alpha1) = -17299.9 Max(alpha1) = 5670.03
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 void Foam::MULES::limit<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double, int, bool) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 Foam::MULES::explicitSolve(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/interFoam"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/interFoam"
Floating point exception (core dumped)
May it be a mesh-problem?
My case:
http://dl.dropbox.com/u/2820596/domain-simple.png
|