update program v1.7 to v2.1.1..help!
 

Good afternoon, I'm trying to convert a program to resolve DQMOM, It was made in v1.7 and I want to change it to v2.1.1 .. when I run the simulation, it is observed this error:

Starting time loop

Courant Number mean: 1.32874256268e-06 max: 1.33361663878e-06
deltaT = 5.99880023995e-06
Time = 5.9988e-06

DILUPBiCG: Solving for WEIGHT1, Initial residual = 0.999999999998, Final residual = 5.08151305829e-17, No Iterations 1
DILUPBiCG: Solving for ETA1, Initial residual = 1, Final residual = 9.98376403815e-17, No Iterations 1
DILUPBiCG: Solving for WEIGHT2, Initial residual = 0.0397686186206, Final residual = 1.44498674972e-18, No Iterations 1
DILUPBiCG: Solving for ETA2, Initial residual = 1, Final residual = 3.18479856381e-16, No Iterations 1
DILUPBiCG: Solving for WEIGHT3, Initial residual = 1, Final residual = 1.15648228185e-16, No Iterations 1
DILUPBiCG: Solving for ETA3, Initial residual = 1, Final residual = 6.19568860656e-17, No Iterations 1
DILUPBiCG: Solving for WEIGHT4, Initial residual = 1, Final residual = 1.24113473485e-16, No Iterations 1
DILUPBiCG: Solving for ETA4, Initial residual = 1, Final residual = 2.28495440613e-17, No Iterations 1
DILUPBiCG: Solving for WEIGHT1, Initial residual = 0.999999971134, Final residual = 8.99399726894e-17, No Iterations 1
DILUPBiCG: Solving for ETA1, Initial residual = 1, Final residual = 1.195727455e-16, No Iterations 1
DILUPBiCG: Solving for WEIGHT2, Initial residual = 0.0522088348692, Final residual = 1.69431964894e-18, No Iterations 1
DILUPBiCG: Solving for ETA2, Initial residual = 1, Final residual = 2.69543587336e-16, No Iterations 1
DILUPBiCG: Solving for WEIGHT3, Initial residual = 1, Final residual = 2.90282428617e-16, No Iterations 1
DILUPBiCG: Solving for ETA3, Initial residual = 1, Final residual = 1.91281127538e-16, No Iterations 1
DILUPBiCG: Solving for WEIGHT4, Initial residual = 1, Final residual = 1.79580712974e-16, No Iterations 1
DILUPBiCG: Solving for ETA4, Initial residual = 1, Final residual = 3.23398178388e-16, No Iterations 1
Dispersed phase volume fraction = 2.14252085634e+20
Dispersed phase Sauter diameter = 576685.058315
ExecutionTime = 0.02 s


Courant Number mean: 1.59433165799e-06 max: 1.60017996454e-06
deltaT = 7.19752392741e-06
Time = 1.31963e-05

DILUPBiCG: Solving for WEIGHT1, Initial residual = 1, Final residual = 1.95826080093e-16, No Iterations 1
DILUPBiCG: Solving for ETA1, Initial residual = 1, Final residual = 9.30882348876e-17, No Iterations 1
DILUPBiCG: Solving for WEIGHT2, Initial residual = 1, Final residual = 6.97758366553e-16, No Iterations 1
DILUPBiCG: Solving for ETA2, Initial residual = 1, Final residual = 9.2284786077e-17, No Iterations 1
DILUPBiCG: Solving for WEIGHT3, Initial residual = 1, Final residual = 6.49750592699e-16, No Iterations 1
DILUPBiCG: Solving for ETA3, Initial residual = 1, Final residual = 2.74127967709e-16, No Iterations 1
DILUPBiCG: Solving for WEIGHT4, Initial residual = 1, Final residual = 1.8273534406e-16, No Iterations 1
DILUPBiCG: Solving for ETA4, Initial residual = 1, Final residual = 9.03106382163e-17, No Iterations 1
#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 in "/lib/x86_64-linux-gnu/libm.so.6"
#4 in "/lib/x86_64-linux-gnu/libm.so.6"
#5 pow in "/lib/x86_64-linux-gnu/libm.so.6"
#6 Foam::aggregationModels::aggMcCoyMadras::aggregKer nel(int) const in "/home/adriana/OpenFOAM/adriana-2.1.1/platforms/linux64GccDPOpt/lib/libAggregationModels.so"
#7
at ~/Documents/OpenFoam/24_10_2012/pbeDev/fields/fields2/fields2.C:107
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
in "/home/adriana/OpenFOAM/adriana-2.1.1/platforms/linux64GccDPOpt/bin/fields2"
Floating point exception (core dumped)

The error appears when are calculated "Dispersed phase volume fraction" and
"Dispersed phase Sauter diameter". I think that probably is because I isn't updating the program correctly ...

I am attaching the fields2.C program..

Some suggestions??


Thanks ;)

Hello Adriana,
as the first timestep is calculated completely, I would guess that it is a problem of numerical nature and not of your update. "Floating point exception" means that some number gets way too big. Sorry, but my only advice could be to check your BCs and numerical settings!
Cheers,
David

Thanks David, I tried to use different numerical settings but didn't work, I believe that the problem is in the program. If you see, dispersed phase fraction diverges completely in first iteration! :

Dispersed phase volume fraction = 2.14252085634e+20
Dispersed phase Sauter diameter = 0.576685058315
ExecutionTime = 0.02 s

what do you think?

Regards,