CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

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

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   October 30, 2012, 15:07
Default update program v1.7 to v2.1.1..help!
  #1
New Member
 
Adriana Vargas
Join Date: Jul 2011
Posts: 8
Rep Power: 6
adriana.nino is on a distinguished road
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:rintStack(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
Attached Files
File Type: c fields2.C (6.2 KB, 7 views)
adriana.nino is offline   Reply With Quote

Old   October 31, 2012, 10:06
Default
  #2
Member
 
David GISEN
Join Date: Jul 2009
Location: Germany
Posts: 37
Rep Power: 8
David* is on a distinguished road
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
David* is offline   Reply With Quote

Old   October 31, 2012, 11:38
Default
  #3
New Member
 
Adriana Vargas
Join Date: Jul 2011
Posts: 8
Rep Power: 6
adriana.nino is on a distinguished road
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,
adriana.nino is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Chemkin Forums? PSR.f program help? johnjohnmulley Main CFD Forum 2 August 5, 2013 07: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 17: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


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