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

Unfamiliar Error

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 26, 2006, 22:16
Default Hello, I am trying to get a
  #1
Member
 
Shaun Cooper
Join Date: Mar 2009
Posts: 54
Rep Power: 8
coops is on a distinguished road
Hello,

I am trying to get a solver to work for my particular application and I have had a lot of trouble. Usually when there are errors I have a rough idea of what I should try to change to eradicate the error but in this case I have no idea. Here is the error output. Can someone please help me by suggesting what might be the issue.

Starting time loop

Mean and max Courant Numbers = 0 -0
deltaT = 0.00048
Time = 0.00048

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
BICCG: Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0
BICCG: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
BICCG: Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0
BICCG: Solving for T, Initial residual = 5.5626846e-286, Final residual = 5.5626846e-286, No Iterations 0
Foam::error::printStack(Foam:: Ostream&)
Foam::sigSegv::sigSegvHandler(int)
[0xd73420]
std::basic_string<char,>, std::allocator<char> > std::operator+<char,>, std::allocator<char> >(char const*, std::basic_string<char,>, std::allocator<char> > const&)
Foam::tmp<foam::fvmatrix<double> > Foam::fvm::ddt<double>(Foam::GeometricField<double ,> const&, Foam::GeometricField<double,>&)
CaseSolver [0x805aa43]
__libc_start_main
__gxx_personality_v0
Segmentation fault

Thanks

Shaun
coops is offline   Reply With Quote

Old   September 27, 2006, 03:21
Default Compile your app with -g -O0 -
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
Compile your app with -g -O0 -DFULLDEBUG, have it dump core and do a 'where' inside gdb.

You have a segmentation fault (see 'man 7 signal') so are probably accessing outside allocated memory.
mattijs is offline   Reply With Quote

Old   October 2, 2006, 03:43
Default I have changed the aachenBomb
  #3
New Member
 
Markus Graf
Join Date: Mar 2009
Posts: 6
Rep Power: 8
graf_m is on a distinguished road
I have changed the aachenBomb case for simulating spray penetration and ignition under cold start conditions. Therefore I built a mesh consisting of many small blocks with lots of internal faces according to the nozzle and glowplug geometry. I have added two more patches but still use the same boundary conditions as in the original case. blockMesh works fine but whenever the solver needs the vertices of the mesh it throws nans out.
Other functions like checking if the injector position is inside the field are working. I checked some parts of the code and found some comments in makeMagSf where it says that producing nans was corrected by implementing subroutine to check if surface area equals zero. But then I got stuck because I am not so familiar with the code yet.

Has Somebody an idea what could be wrong?

Markus

/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.3 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : dieselFoam /home/markus/OpenFOAM/markus-1.3/run/dieselColdStart aachenBomb
Date : Oct 01 2006
Time : 12:13:03
Host : linux
PID : 4899
Root : /home/markus/OpenFOAM/markus-1.3/run/dieselColdStart
Case : aachenBomb
Nprocs : 1
Create time

Create mesh for time = 0


Reading thermophysicalProperties
Selecting thermodynamics package hMixtureThermo<reactingmixture>
Selecting chemistryReader chemkinReader
Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting turbulence model kEpsilon
Creating field DpDt

Constructing chemical mechanism
Selecting ODE solver SIBS
chemistryModel::chemistryModel: Number of species = 5 and reactions = 1

Reading environmentalProperties
Reading combustion properties

Constructing Spray
Selecting injectorType commonRailInjector
injectionPressureProfile_.size() = 2, massFlowRateProfile_.size() = 23
end constructor. in commonRail
Selecting atomizationModel off
Selecting dragModel standardDragModel
Selecting evaporationModel standardEvaporationModel
Selecting heatTransferModel RanzMarshall
Selecting wallModel reflect
Selecting breakupModel ReitzKHRT
Selecting collisionModel off
Selecting dispersionModel off
Selecting injectorModel hollowConeInjector
Selecting pdfType RosinRammler
Average Velocity for injector 0: 375.808 m/s, injection pressure = 500 bar
Constructing three dimensional spray injection.
Mean and max Courant Numbers = nan nan

Starting time loop

Mean and max Courant Numbers = nan nan
deltaT = 2.94118e-06
Time = 2.94118e-06

Evolving Spray
Solving chemistry
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
BICCG: Solving for Ux: solution singularity
BICCG: Solving for Uy: solution singularity
BICCG: Solving for Uz: solution singularity
BICCG: Solving for C7H16: solution singularity
BICCG: Solving for O2: solution singularity
BICCG: Solving for CO2: solution singularity
BICCG: Solving for H2O: solution singularity
BICCG: Solving for h: solution singularity
ICCG: Solving for p: solution singularity
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = nan, global = nan, cumulative = nan
ICCG: Solving for p: solution singularity
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = nan, global = nan, cumulative = nan
BICCG: Solving for epsilon: solution singularity
BICCG: Solving for k: solution singularity

Number of parcels in system | 1
Injected liquid mass....... | 0.00107761 mg
Liquid Mass in system...... | 0.00107761 mg
SMD, Dmax.................. | 82.4948 mu, 82.4948 mu
Added gas mass = nan mg
Evaporation Continuity Error| nan mg
ExecutionTime = 19.78 s ClockTime = 20 s

Mean and max Courant Numbers = nan nan
deltaT = 3.36134e-06
Time = 6.30252e-06

Evolving Spray
interpolationCellPointFace<type>::interpolate(cons t vector&, const label nCell) const : search failed; using closest cellFace value
cell number 1040 position (nan nan nan)
interpolationCellPointFace<type>::interpolate(cons t vector&, const label nCell) const : search failed; using closest cellFace value
cell number 903 position (nan nan nan)
interpolationCellPointFace<type>::interpolate(cons t vector&, const label nCell) const : search failed; using closest cellFace value
cell number 1040 position (nan nan nan)
interpolationCellPointFace<type>::interpolate(cons t vector&, const label nCell) const : search failed; using closest cellFace value
cell number 903 position (nan nan nan)
interpolationCellPointFace<type>::interpolate(cons t vector&, const label nCell) const : search failed; using closest cellFace value
cell number 1040 position (nan nan nan)
interpolationCellPointFace<type>::interpolate(cons t vector&, const label nCell) const : search failed; using closest cellFace value
cell number 903 position (nan nan nan)
graf_m is offline   Reply With Quote

Old   October 3, 2006, 03:34
Default I guess you just didn't mentio
  #4
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
I guess you just didn't mention it: you DID check the mesh with checkMesh?
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   October 3, 2006, 09:46
Default I did check the mesh with chec
  #5
New Member
 
Markus Graf
Join Date: Mar 2009
Posts: 6
Rep Power: 8
graf_m is on a distinguished road
I did check the mesh with checkMesh and it was Ok.
graf_m is offline   Reply With Quote

Old   October 3, 2006, 10:10
Default OK. Sorry for asking. Anoth
  #6
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
OK. Sorry for asking.

Another thing I usually do when having problems with a new case for a complex solver and I'm unsure whether the problem is in the case: I try the solver the complex solver was derived from (or a solver with a subset of the functionality) on that case. In your case that would be one of the compressible solvers. If you still have problems then propably the problem is with the definition of the fluid boundary conditions and not with the spray (if that solver works the inverse is not necessarily true, I'm afraid, but the propability is high).

BTW: the interior field for p is NOT 0 I hope (zero never works for me with compressible solvers ;) ). Just making sure. Sorry for asking.

PS: I'm sorry I can't be more specific, but I never worked with dieselFoam (only used it to extract ideas)
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   October 3, 2006, 12:38
Default I will follow up your advice b
  #7
New Member
 
Markus Graf
Join Date: Mar 2009
Posts: 6
Rep Power: 8
graf_m is on a distinguished road
I will follow up your advice but I doubt that I will be successful because the nan problem starts with calculating the courant number which is needed by other compressible solvers as well. My boundary conditions are simple: zero flow field, pressure field 40 bar, T field 523 K, zero gradients for T,p and zero value for U at walls.

We will see what happens
graf_m is offline   Reply With Quote

Old   October 3, 2006, 13:54
Default You were right. The solver rho
  #8
New Member
 
Markus Graf
Join Date: Mar 2009
Posts: 6
Rep Power: 8
graf_m is on a distinguished road
You were right. The solver rhoTurbFoam does the job. The problem obviously occurs within dieselFoam. Hm ... my bc for dieselFoam should also be Ok. Anyway, thanx so far ...
graf_m 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
Error: Floating point error: invalid number fpingqian FLUENT 4 February 8, 2012 02:20
Errors when Compiling UDF: error C2040/error C2099 Julian K. FLUENT 1 December 21, 2008 01:23
"Error: Floating point error: invalid number" MI Kim FLUENT 2 January 4, 2007 11:00
Fatal error error writing to tmp No space left on device maka OpenFOAM Installation 2 April 3, 2006 08:48
Error: Floating point error: invalid number Bob FLUENT 3 June 3, 2005 18:11


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