CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   sonicFOAM for an airfoil error (http://www.cfd-online.com/Forums/openfoam/79887-sonicfoam-airfoil-error.html)

inginer September 7, 2010 13:26

sonicFOAM for an airfoil error
 
Hello,

I'm trying to do a sonicFoam case in OF 1.7.0 with a mesh imported from fluent. it is an airfoil 2d with M=0.8.
after building the case i try to run it and i receive the following error


openfoam@linux-bkii:~/OpenFOAM/openfoam-1.7.0/run/projects/cmyairfoil> sonicFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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-21131bcbd876
Exec : sonicFoam
Date : Sep 07 2010
Time : 16:10:12
Host : linux-bkii
PID : 8245
Case : /home/openfoam/OpenFOAM/openfoam-1.7.0/run/projects/cmyairfoil
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<sutherlandTransport<specieT hermo<hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model SpalartAllmaras
SpalartAllmarasCoeffs
{
sigmaNut 0.66666;
kappa 0.41;
Prt 1;
Cb1 0.1355;
Cb2 0.622;
Cw2 0.3;
Cw3 2;
Cv1 7.1;
Cv2 5;
}


Starting time loop

Time = 2.13e-09

Courant Number mean: 5.01425e-06 max: 0.00279144
#0 Foam::error::printStack(Foam::Ostream&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::operator/(Foam::UList<double> const&, Foam::UList<double> const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linuxGccDPOpt/libOpenFOAM.so"
#5 Foam::diagonalSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linuxGccDPOpt/libOpenFOAM.so"
#6 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linuxGccDPOpt/libfiniteVolume.so"
#7 Foam::lduMatrix::solverPerformance Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) in "/home/openfoam/OpenFOAM/openfoam-1.7.0/applications/bin/linuxGccDPOpt/sonicFoam"
#8 main in "/home/openfoam/OpenFOAM/openfoam-1.7.0/applications/bin/linuxGccDPOpt/sonicFoam"
#9 __libc_start_main in "/lib/libc.so.6"
#10 _start at /usr/src/packages/BUILD/glibc-2.10.1/csu/../sysdeps/i386/elf/start.S:122
Floating point exception
openfoam@linux-bkii:~/OpenFOAM/openfoam-1.7.0/run/projects/cmyairfoil>


anybody receive a similar error? what i did wrong???
Thk you


Ovidiu

bruxellois September 7, 2010 18:39

Hi,

I receive a similar error when trying to run sonicFoam to my case.

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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-279cc8e8233b
Exec : sonicFoam
Date : Sep 08 2010
Time : 00:32:20
Host :
PID : 2229
Case :
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>>>>>
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam170/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam170/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::eConstThermo<Foam:: perfectGas> > > > >::calculate() in "/opt/openfoam170/lib/linuxGccDPOpt/libbasicThermophysicalModels.so"
#4 Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::eConstThermo<Foam:: perfectGas> > > > >::ePsiThermo(Foam::fvMesh const&) in "/opt/openfoam170/lib/linuxGccDPOpt/libbasicThermophysicalModels.so"
#5 Foam::basicPsiThermo::addfvMeshConstructorToTable< Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::eConstThermo<Foam:: perfectGas> > > > > >::New(Foam::fvMesh const&) in "/opt/openfoam170/lib/linuxGccDPOpt/libbasicThermophysicalModels.so"
#6 Foam::basicPsiThermo::New(Foam::fvMesh const&) in "/opt/openfoam170/lib/linuxGccDPOpt/libbasicThermophysicalModels.so"
#7
in "/opt/openfoam170/applications/bin/linuxGccDPOpt/sonicFoam".

Any idea please.

inginer September 9, 2010 11:13

Hy Bruxellois,

Did you manage to repare the error. I try to do the same case but in incompressible, simpleFoam. the code is running but the results were very strange. The velocity U is from 0 to +infinity and the pressure is from -infinity to +infinity; this is totally wrong, but i can see the error in the code.

Maybe OF is not capable to solve transonic case ...

Best Regards,
Ovidiu

kiran September 10, 2010 15:56

Hi
Check with your boundary conditions, turbuence model you are using

Regards
Ambilpur Kiran

inginer September 13, 2010 11:34

Hello kiran,

my turbulence model is Spalart-Allmaras for compressible flow, the temp is 300K and velocity 277.67 m/s (M=0.8).
i check the boundary condition. i receive the same error. i don't have any ideas in what could be wrong.

Best Regards,
Ovidiu

kiran September 13, 2010 14:02

attach your case file

regards
Ambilpur Kiran

inginer September 14, 2010 05:31

Hello Kiran,

bellow you will find a link where you can download the case.

http://www.2shared.com/file/fVbh4wcU/comairfoil.html

Best Regards,
Ovidiu


All times are GMT -4. The time now is 04:51.