Using nanoscale with rhoPisoTwinParcelFoam
 

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::printStack(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