|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Adriana Vargas
Join Date: Jul 2011
Posts: 8
Rep Power: 15 ![]() |
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: ![]() #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 ![]() |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Member
David GISEN
Join Date: Jul 2009
Location: Germany
Posts: 70
Rep Power: 17 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Adriana Vargas
Join Date: Jul 2011
Posts: 8
Rep Power: 15 ![]() |
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, |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Chemkin Forums? PSR.f program help? | johnjohnmulley | Main CFD Forum | 2 | August 5, 2013 08:53 |
error in program: | mani.rajendiran | Fluent UDF and Scheme Programming | 0 | March 22, 2012 12:25 |
Which opensource and free CAD program? | Xwang | Main CFD Forum | 1 | October 24, 2010 18:29 |
Update boundary conditions calculated by an external program | CedricVH | OpenFOAM | 2 | January 15, 2010 12:55 |
Where can I write my C program "UDF" in Fluent | Chimwemwe Banda | FLUENT | 2 | December 8, 2004 10:07 |