|
[Sponsors] |
June 9, 2011, 10:35 |
Using nanoscale with rhoPisoTwinParcelFoam
|
#1 |
New Member
Join Date: Apr 2011
Posts: 10
Rep Power: 14 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|