[Horus]

Hello,

I had successfully simulated a simple using icoFoam (the geometry was created with CATIA V5, meshing with netgen, mesh conversion with netgenNeutralToFoam).

Now I'm trying to adapt that case for sonicFoam. I've used the sonicFoam/laminar/forwardStep example as an exemplar, almost all settings are identical. The entire case is available at http://xgm.de/nozzle.zip

When running sonicFoam it crashes:

florian@horus:~/OpenFOAM/run/nozzle> sonicFoam 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.6                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 1.6-f802ff2d6c5a                                                      
Exec   : sonicFoam                                                             
Date   : Feb 19 2010                                                           
Time   : 17:36:39                                                              
Host   : horus                                                                 
PID    : 5568                                                                  
Case   : /home/florian/OpenFOAM/run/nozzle                                     
nProcs : 1                                                                     
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).             

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

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package ePsiThermo<pureMixture<constTransport<specieThermo<eConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type laminar

Starting time loop

Time = 0.002

#0  Foam::error::printStack(Foam::Ostream&) in "/home/florian/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1  Foam::sigFpe::sigFpeHandler(int) in "/home/florian/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2  ?? in "/lib64/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/florian/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#4  void Foam::divide<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/florian/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/sonicFoam"
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator/<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/home/florian/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/sonicFoam"
#6  main in "/home/florian/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/sonicFoam"
#7  __libc_start_main in "/lib64/libc.so.6"
#8  _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception

I hope anyone got an idea what could be wrong.

Thanks!

Florian

[ziooscuro]

are your BC ok??