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

Using nanoscale with rhoPisoTwinParcelFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 9, 2011, 09:35
Default Using nanoscale with rhoPisoTwinParcelFoam
  #1
New Member
 
Join Date: Apr 2011
Posts: 10
Rep Power: 6
jabecker is on a distinguished road
Hello Foamers,

I posted earlier in regards to this topic, but now another problem has arisen. I keep getting a "floating point exception" error and I cannot get rid of it. I have done some serious searching and digging on the forum and the Internet and I have found nothing conclusive.

So again, my case is a cube filled with air that has a length, width and height of 1e-06m. I use manual injection to place stationary particles of SnO2 inside of the box . The closed box contains a heated bottom plate at 573K, but the air and particles all start at room temperature of 298K. I want to investigate the heat transfer through the particles and the air. I am using OF 1.7.0

The exact error I get is:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.7.0 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.7.0-113391ee57bd
Exec : rhoPisoTwinParcelFoam
Date : Jun 09 2011
Time : 14:18:40
Host : alfred
PID : 5756
Case : /home/becker/tutorials/lagrangian/rhoPisoTwinParcelFoam/NANOsimplifiedSiwek
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g
Reading thermophysical properties

Selecting thermodynamics package hsPsiThermo<pureMixture<sutherlandTransport<specie Thermo<janafThermo<perfectGas>>>>>

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type laminar
Creating field DpDt

Constructing thermoCloud1
--> FOAM Warning :
From function Cloud<ParticleType>::initCloud(const bool checkClass)
in file /opt/openfoam170/src/lagrangian/basic/lnInclude/CloudIO.C at line 124
Cannot read particle positions file
"/home/becker/tutorials/lagrangian/rhoPisoTwinParcelFoam/NANOsimplifiedSiwek/0/lagrangian/thermoCloud1"
assuming the initial cloud contains 0 particles.
Selecting DispersionModel none
Selecting DragModel none
Selecting InjectionModel ManualInjection
Constructing 3-D injection
Selecting pdfType RosinRammler
Selecting PatchInteractionModel StandardWallInteraction
Selecting PostProcessingModel none
Selecting U IntegrationScheme Euler
Selecting HeatTransferModel RanzMarshall
Selecting T IntegrationScheme Analytical
Constructing kinematicCloud1
--> FOAM Warning :
From function Cloud<ParticleType>::initCloud(const bool checkClass)
in file /opt/openfoam170/src/lagrangian/basic/lnInclude/CloudIO.C at line 124
Cannot read particle positions file
"/home/becker/tutorials/lagrangian/rhoPisoTwinParcelFoam/NANOsimplifiedSiwek/0/lagrangian/kinematicCloud1"
assuming the initial cloud contains 0 particles.
Selecting DispersionModel none
Selecting DragModel SphereDrag
Selecting InjectionModel ManualInjection
Constructing 3-D injection
Selecting pdfType RosinRammler
Selecting PatchInteractionModel StandardWallInteraction
Selecting PostProcessingModel none
Selecting U IntegrationScheme Euler
Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 Foam::Time::adjustDeltaT() in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#4
in "/home/becker/OpenFOAM/becker-1.7.0/applications/bin/linux64GccDPOpt/rhoPisoTwinParcelFoam"
#5 __libc_start_main in "/lib/libc.so.6"
#6
in "/home/becker/OpenFOAM/becker-1.7.0/applications/bin/linux64GccDPOpt/rhoPisoTwinParcelFoam"
Floating point exception

So my C++ skills are minimal, but from the error I assumed there was something wrong with the DeltaT. I varied this greatly and still the error remained. I then tried making the mesh smaller and that did not work as well. I tried changing the maxCo from 0.2 to 0.9... I tried also changing the Cp and mu values to be temperature dependent and still nothing. I saw a suggestion about checking the initial conditions, but the program runs perfectly fine with it is 1m instead of 1e-06m, so it does not make sense to me that this could be the problem.

The only suggestion I found that I have not tried is to use lower order differencing. This did not make sense to me and I am not sure what it means. Maybe someone can help me in this regard.

As you can see, I am getting very frustrated and running out of ideas. So please, if you have any other suggestions I would GREATLY appreciated them. I can send you any files that you need in order to help fix the problem.

Thank you in advance for your help and assistance in the matter. I really look forward to hearing from someone,

Jen
jabecker 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



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