New densitybased solver AeroFoam
 

Hi all,

As a part of our M. Sc. Thesis at Politecnico di Milano, I and Elisa Serioli developed a new density-based solver for compressible flows called AeroFoam and we tested it on a number of 2D and 3D aerodynamic and aeroelastic benchmark test problems.

Should you be curious, we finally managed to upload a web page presenting the project at www.aero.polimi.it/freecase/. There you can also download the source code, a lot of test cases ready to be run and a quickstart guide.

Let us know if you have any comments, questions, suggestions.

Best regards.

Hi Guilo, thanks for publis
 

Hi Guilo,

thanks for publishing all the stuff, looks great!

In your presentation you mention, that you will perform computations of a Apollo reentry capsule, right?! I am involved in a related project, about boundary layer transition on a Apollo reentry capsule at the Aerodynamics Lab at TU Delft and I will (try to) compute an Apollo reentry capsule in Mach 7.5 with Fluent. Until now I have validated Fluent on a more simple blunt nosed geometry in inviscid supersonic and hypersonic flow.

I am highly interested in collaboration and comparison of results for the Apollo reentry capsule! What do you think about it?!

with best regards,
Florian

Hi Florian, The project you
 

Hi Florian,

The project you are involved in is terrific.
Unfortunately at the moment the capability of solving hypersonic flows is still a ToDo. In fact the flux function Jacobian matrix eigenvectors necessary to build the 1st order (Roe) and 2nd order (LW) numerical fluxes at the moment are written assuming a thermodynamic model of polytropic ideal gas.
To go above M = 5 it should be necessary to modify those accordingly to a thermodynamic model of mixture of reacting real gases at least in thermochemical equilibrium. It is not a too difficult task (we have coded it in the quasi-1D case in Fortran 90) but it requires some time. What do you think? Anyway, do you have the NASA TN D-4688 report?

Best regards.

Hi Giulio, actually I am (f
 

Hi Giulio,

actually I am (fortunatley) not going to include any effects due to chemical reactions or vibrational excitation in my calculations. My plan is to obtain the transition Reynolds number along the capsule front surface with the aid of windtunnel test in a Ludwig tube under cold flow conditions, hence no high temperature effects.

Reynolds numbers(fluent)
+ boundary layer transition zone (windtunnel)
= transition Reynolds number along the capsule front

Indeed, I have the mentioned NASA report just loaded a few weeks before.

best regards

Hi Giulio I am interested i
 

Hi Giulio

I am interested in FSI alanysis so as to couple the results to a reanalysis via an optimization engine acting in the definition of the geometry...

Analyze CFD->Analyze FEA->rank design->perturbate design parameters->redefine Geom->remesh... and it iterates like this

The design perturbation would be the result of a genetic evolution of the population of "feasible designs", and the ranking procedure would be the result of a quadratic or logarithmic penalty function applied to some index of merit in my analysis,,,,

Where can I download the cae2d and the cae3d toolboxes...? is there any link on http://www.aero.polimi.it/freecase/?FreeCASE...

My work has been done around Calculix as a FEA solver... I am not very familiar with Code Aster yet...

Best regards

Alex

Hi Alex, Cool the shape opt
 

Hi Alex,

Cool the shape optimization thing... If the displacements of the boundaries are not too big you could simply deform the mesh, rather than remesh everything, no? In that case you may find useful the showDisplacement utility included in the AeroFoam.tar.gz package.

The CAE2D and CAE3D toolboxes are included in the AeroFoam.tar.gz package. To use them you must remember to change the #define CAE2D 0 to #define CAE2D 1 in the top section of AeroFoam.C.

Best regards,

Giulio

Hi Giulio Thanks... Good j
 

Hi Giulio

Thanks... Good job... :-)

Best regards

Alex

Hi Giulio, Just one questi
 

Hi Giulio,

Just one question, If you have time to answer of course... :-)

And I beg all you Foamers to forgive me for using this post to make this question... but I do not know where to post it.. :-(

I have ( tried ) to make a solver, ant the exe compiles ok.. meaning that it resolves the dependencies, the header files are ok... but when I execute the damn thing, a message of error just pos out... and I do not know wether this is a class grammatics related error, or just a numeric error... I do not know enough c++ to understand what is happening ... The message is...


root@iskandhar:/home/alex/Desktop/wkshp# mhdIcoFoam . mhdDuctCase
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : mhdIcoFoam . mhdDuctCase
Date : Jan 11 2009
Time : 17:52:38
Host : iskandhar
PID : 7327
Root : /home/alex/Desktop/wkshp
Case : mhdDuctCase
Nprocs : 1
Create time

Create mesh for time = 0

Reading transportProperties

Reading field sigma

Reading field ElPot

Reading field A

Calculate electrical current density Je

Calculating magnetic field B

Calculate Lorentz Force F

#0 Foam::error::printStack(Foam:stream&) in "/root/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/root/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xffffe420]
#3 Foam::tmp<foam::geometricfield<foam::crossproduct< foam::vector<double>, Foam::Vector<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::operator^<foam::vector<double>, Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<foam::vector<d ouble>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/root/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/mhdIcoFoam"
#4 main in "/root/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/mhdIcoFoam"
#5 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#6 Foam::regIOobject::readIfModified() in "/root/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/mhdIcoFoam"
Floating point exception (core dumped)
root@iskandhar:/home/alex/Desktop/wkshp#

And I am just defining F=Je^B in the source code... either Je and B are volVectorFields, F is defined in createFields.H as

Info << "Calculate Lorentz Force F \n" << endl;
volVectorField F
(
IOobject
(
"F",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
Je^B
);


Is this syntactically correct...?

I just do not know wether this is a numeric build up error, or a C++ grammatical one...

Could you help me here...?

Thanks in advance

Alex

Hi Giulio... Solved it alre
 

Hi Giulio...

Solved it already... it was just a badly built blockMesh domain... :-)

Thankz anyway.

Regards,
Alex

Hi Giulio Taking ypour Idea
 

Hi Giulio

Taking ypour Idea of a hypersonic reacting flow solver,

which solver do you think would be more adapted to the necessary changes, (arbitrary equation of state, diffusive pressure equation) reactingFoam or rhoSonicFoam... ?

Best regards

Alex

Hi Giulio Taking ypour Idea
 

Hi Giulio

Taking ypour Idea of a hypersonic reacting flow solver,

which solver do you think would be more adapted to the necessary changes, (arbitrary equation of state, diffusive pressure equation) reactingFoam or rhoSonicFoam... ?

Best regards

Alex

Running AeroFoam parallel
 

hi

Thanks in bunch for the wonderful AeroFoam solver.
I have downloaded and compiled with OpenFoam 1.4.1. It is working fine when excuteed serially but I am facing problem when I use parallel option

The command I have used is for Onera Case..

mpirun -np 4 AeroFoam . Onera_CASE -parallel


The error I got is (sorry for such a long output posting)


MPI Pstream initialized with:
floatTransfer : 1
nProcsSimpleSum : 0
scheduledTransfer : 0

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

Exec : AeroFoam . OneraM6_parallel -parallel
[0] Date : Apr 02 2009
[0] Time : 13:10:50
[0] Host : n5
[0] PID : 29450
[1] Date : Apr 02 2009
[1] Time : 13:10:50
[1] Host : n5
[1] PID : 29451
[2] Date : Apr 02 2009
[2] Time : 13:10:50
[2] Host : n5
[2] PID : 29452
[3] Date : Apr 02 2009
[3] Time : 13:10:50
[3] Host : n5
[3] PID : 29453
[3] Root : /HOME/crl_cfd/santosh/comp_sol/OneraM6
[3] Case : OneraM6_parallel
[3] Nprocs : 4
[1] Root : /HOME/crl_cfd/santosh/comp_sol/OneraM6
[1] Case : OneraM6_parallel
[1] Nprocs : 4
[2] Root : /HOME/crl_cfd/santosh/comp_sol/OneraM6
[2] Case : OneraM6_parallel
[2] Nprocs : 4
[0] Root : /HOME/crl_cfd/santosh/comp_sol/OneraM6
[0] Case : OneraM6_parallel
[0] Nprocs : 4
[0] Slaves :
[0] 3
[0] (
[0] n5.29451
[0] n5.29452
[0] n5.29453
[0] )
[0]
Create time

Create mesh for time = 0

========================================
Reading Thermodynamic Properties...
========================================
gamma [-] = 1.4
Pr [-] = 0.72
R [J/KgK] = 287.05
Cp [J/KgK] = 717.625
Cv [J/KgK] = 1004.67
----------------------------------------

========================================
Reading Initial Conditions...
========================================
Reading field p
Reading field T
Reading field U
----------------------------------------
Creating conservative variables at t^(k)
Creating conservative variables at t^(k-1)
Creating conservative variables at t^(k-2)
Creating conservative variables fluxes at t^(k)
Creating primitive variables gradient fields (optional)
----------------------------------------

========================================
Creating Mesh Connectivity...
========================================
Using standard algorithm...
Loop on Internal Faces...
10 %
10 %
10 %
20 %
20 %
20 %
30 %
30 %
30 %
40 %
40 %
40 %
50 %
50 %
50 %
60 %
60 %
60 %
70 %
70 %
70 %
80 %
80 %
80 %
90 %
90 %
90 %
100 %
10 %
20 %
30 %
40 %
50 %
100 %
100 %
----------------------------------------
Loop on Boundary Faces...
10 %
10 %
60 %
20 %
20 %
70 %
30 %
30 %
80 %
40 %
40 %
90 %
50 %
50 %
100 %
60 %
60 %
70 %
70 %
80 %
80 %
90 %
90 %
100 %
100 %
----------------------------------------

========================================
Loading Mesh Connectivity...
========================================
Done!
----------------------------------------

========================================
Mesh
========================================
# of Faces = 162125
# of Cells = 85908
# of BndPatch = 10
0) slip @ wing
1) symmetryPlane @ symmetry
2) extrapolatedOutlet @ lateral
3) extrapolatedOutlet @ top
4) extrapolatedOutlet @ outlet
5) supersonicInlet @ bottom
6) supersonicInlet @ inlet
7) @ procBoundary0to1
8) @ procBoundary0to2
9) @ procBoundary0to3
----------------------------------------

========================================
Numerical Scheme
========================================
flowType = Euler
timeScheme = RK2
CourantMax = 0
CourantSteps = 25
MonotoneFlux = Roe
HighResFlux = LW
fluxLimiter = VL
entropyFix = HH2
extrapolateBC = 1
residualNorm = L1
minResidual = -1
----------------------------------------

Starting time loop...
10 %
20 %
30 %
40 %
50 %
60 %
70 %
80 %
90 %
100 %
10 %
20 %
30 %
40 %
50 %
60 %
70 %
80 %
90 %
100 %
[1] #0 Foam::error::printStack(Foam::Ostream&)[2] #0 Foam::error::printStack(Foam::Ostream&)[3] #0 Foam::error::printStack(Foam::Ostream&)[0] #0 Foam::error::printStack(Foam::Ostream&) in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigFpeHandler(int) in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[2] #1 Foam::sigFpe::sigFpeHandler(int) in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigFpeHandler(int) in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigFpeHandler(int) in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #2 [2] #2 in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #2 in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[0] #2 __restore_rt__restore_rt__restore_rt__restore_rt in "/lib64/tls/libc.so.6"
[1] #3 in "/lib64/tls/libc.so.6"
[2] #3 in "/lib64/tls/libc.so.6"
[3] #3 in "/lib64/tls/libc.so.6"
[0] #3 lwhr_flux(double, double, double, Foam::Vector<double>, Foam::Vector<double>, double, double, double, double, Foam::Vector<double>, Foam::Vector<double>, double, double, Foam::Vector<double>, double, double, double*, Foam::Vector<double>*, double*, Foam::word)lwhr_flux(double, double, double, Foam::Vector<double>, Foam::Vector<double>, double, double, double, double, Foam::Vector<double>, Foam::Vector<double>, double, double, Foam::Vector<double>, double, double, double*, Foam::Vector<double>*, double*, Foam::word)lwhr_flux(double, double, double, Foam::Vector<double>, Foam::Vector<double>, double, double, double, double, Foam::Vector<double>, Foam::Vector<double>, double, double, Foam::Vector<double>, double, double, double*, Foam::Vector<double>*, double*, Foam::word)lwhr_flux(double, double, double, Foam::Vector<double>, Foam::Vector<double>, double, double, double, double, Foam::Vector<double>, Foam::Vector<double>, double, double, Foam::Vector<double>, double, double, double*, Foam::Vector<double>*, double*, Foam::word) in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/AeroFoam"
[1] #4 in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/AeroFoam"
[2] #4 in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/AeroFoam"
[0] #4 in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/AeroFoam"
[3] #4 build_fluxes(Foam::Field<double>*, Foam::word, Foam::dimensioned<double>*, Foam::dimensioned<double>*, Foam::dimensioned<double>*, Foam::fvMesh*, Foam::Field<int>*, Foam::Field<int>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>*build_fluxes(Foam::Field<double >*, Foam::word, Foam::dimensioned<double>*, Foam::dimensioned<double>*, Foam::dimensioned<double>*, Foam::fvMesh*, Foam::Field<int>*, Foam::Field<int>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>*build_fluxes(Foam::Field<double >*, Foam::word, Foam::dimensioned<double>*, Foam::dimensioned<double>*, Foam::dimensioned<double>*, Foam::fvMesh*, Foam::Field<int>*, Foam::Field<int>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>*build_fluxes(Foam::Field<double >*, Foam::word, Foam::dimensioned<double>*, Foam::dimensioned<double>*, Foam::dimensioned<double>*, Foam::fvMesh*, Foam::Field<int>*, Foam::Field<int>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>* at Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>*, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>*, Foam::word, Foam::word, Foam::word, Foam::word, double)
[2] #5 at Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>*, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>*, Foam::word, Foam::word, Foam::word, Foam::word, double)
[0] #5 at Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>*, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>*, Foam::word, Foam::word, Foam::word, Foam::word, double)
[1] #5 at Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>*, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>*, Foam::word, Foam::word, Foam::word, Foam::word, double)
[3] #5 mainmainmainmain in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/AeroFoam"
[2] #6 __libc_start_main in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/AeroFoam"
[0] #6 __libc_start_main in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/AeroFoam"
[1] #6 __libc_start_main in "/lib64/tls/libc.so.6"
[2] #7 in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/AeroFoam"
[3] #6 __libc_start_main in "/lib64/tls/libc.so.6"
[0] #7 in "/lib64/tls/libc.so.6"
[1] #7 __gxx_personality_v0__gxx_personality_v0 in "/lib64/tls/libc.so.6"
[3] #7 __gxx_personality_v0__gxx_personality_v0 in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/AeroFoam"
[n5:29450] *** Process received signal ***
[n5:29450] Signal: Floating point exception (8)
[n5:29450] Signal code: (-6)
[n5:29450] Failing at address: 0x21a0000730a
[n5:29450] [ 0] /lib64/tls/libc.so.6 [0x399e12e2b0]
[n5:29450] [ 1] /lib64/tls/libc.so.6(gsignal+0x3d) [0x399e12e21d]
[n5:29450] [ 2] /lib64/tls/libc.so.6 [0x399e12e2b0]
[n5:29450] [ 3] AeroFoam [0x419e90]
[n5:29450] [ 4] AeroFoam [0x421e80]
[n5:29450] [ 5] AeroFoam [0x43273e]
[n5:29450] [ 6] /lib64/tls/libc.so.6(__libc_start_main+0xdb) [0x399e11c3fb]
[n5:29450] [ 7] AeroFoam(__gxx_personality_v0+0xe2) [0x410aea]
[n5:29450] *** End of error message ***
in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/AeroFoam"
[n5:29452] *** Process received signal ***
[n5:29452] Signal: Floating point exception (8)
[n5:29452] Signal code: (-6)
[n5:29452] Failing at address: 0x21a0000730c
[n5:29452] [ 0] /lib64/tls/libc.so.6 [0x399e12e2b0]
[n5:29452] [ 1] /lib64/tls/libc.so.6(gsignal+0x3d) [0x399e12e21d]
[n5:29452] [ 2] /lib64/tls/libc.so.6 [0x399e12e2b0]
[n5:29452] [ 3] AeroFoam [0x419e90]
[n5:29452] [ 4] AeroFoam [0x421e80]
[n5:29452] [ 5] AeroFoam [0x43273e]
[n5:29452] [ 6] /lib64/tls/libc.so.6(__libc_start_main+0xdb) [0x399e11c3fb]
[n5:29452] [ 7] AeroFoam(__gxx_personality_v0+0xe2) [0x410aea]
[n5:29452] *** End of error message ***
in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/AeroFoam"
[n5:29451] *** Process received signal ***
[n5:29451] Signal: Floating point exception (8)
[n5:29451] Signal code: (-6)
[n5:29451] Failing at address: 0x21a0000730b
[n5:29451] [ 0] /lib64/tls/libc.so.6 [0x399e12e2b0]
[n5:29451] [ 1] /lib64/tls/libc.so.6(gsignal+0x3d) [0x399e12e21d]
[n5:29451] [ 2] /lib64/tls/libc.so.6 [0x399e12e2b0]
[n5:29451] [ 3] AeroFoam [0x419e90]
[n5:29451] [ 4] AeroFoam [0x421e80]
[n5:29451] [ 5] AeroFoam [0x43273e]
[n5:29451] [ 6] /lib64/tls/libc.so.6(__libc_start_main+0xdb) [0x399e11c3fb]
[n5:29451] [ 7] AeroFoam(__gxx_personality_v0+0xe2) [0x410aea]
[n5:29451] *** End of error message ***
in "/scratch/crl_cfd/sagar/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/AeroFoam"
[n5:29453] *** Process received signal ***
[n5:29453] Signal: Floating point exception (8)
[n5:29453] Signal code: (-6)
[n5:29453] Failing at address: 0x21a0000730d
[n5:29453] [ 0] /lib64/tls/libc.so.6 [0x399e12e2b0]
[n5:29453] [ 1] /lib64/tls/libc.so.6(gsignal+0x3d) [0x399e12e21d]
[n5:29453] [ 2] /lib64/tls/libc.so.6 [0x399e12e2b0]
[n5:29453] [ 3] AeroFoam [0x419e90]
[n5:29453] [ 4] AeroFoam [0x421e80]
[n5:29453] [ 5] AeroFoam [0x43273e]
[n5:29453] [ 6] /lib64/tls/libc.so.6(__libc_start_main+0xdb) [0x399e11c3fb]
[n5:29453] [ 7] AeroFoam(__gxx_personality_v0+0xe2) [0x410aea]
[n5:29453] *** End of error message ***
[n5:29447] [0,0,0] ORTE_ERROR_LOG: Timeout in file base/pls_base_orted_cmds.c at line 275
[n5:29447] [0,0,0] ORTE_ERROR_LOG: Timeout in file pls_rsh_module.c at line 1164
[n5:29447] [0,0,0] ORTE_ERROR_LOG: Timeout in file errmgr_hnp.c at line 90
mpirun noticed that job rank 0 with PID 29450 on node n5 exited on signal 8 (Floating point exception).
1 additional process aborted (not shown)
[n5:29447] [0,0,0] ORTE_ERROR_LOG: Timeout in file base/pls_base_orted_cmds.c at line 188
[n5:29447] [0,0,0] ORTE_ERROR_LOG: Timeout in file pls_rsh_module.c at line 1196
--------------------------------------------------------------------------
mpirun was unable to cleanly terminate the daemons for this job. Returned value Timeout instead of ORTE_SUCCESS.

--------------------------------------------------------------------------


Can any body please help me out...

Update on AeroFoam
 

Has there been any updates to AeroFoam? Specifically porting it OF 1.5 and the parallelization issues?

I am also interested to know whether it is possible to use AeroFoam in version 1.5. Has anyone tested it?

I'd also be interested in a version which is compatible with OF1.5

Hi Giulio

I would also be interested in a Version of AeroFoam for OpenFOAM-1.5 or 1.6.x...

Thanks in advance

Alex

Hi Chris, Paul, Julian, Alexandre.

Sorry for the late reply. Yes, there is a AeroFoam version for OpenFOAM 1.5-dev fully parallel with some new options such as Residual Smoothing for steady-state computations. I will try to update the website in the next few days.

Best regards,

Giulio

That sounds very good! Thanks Giulio!

Thanks Giulio

Does it couple with Aster 10.0...? would be awesome... :-)

Thats a gre8 news. I was waiting this to happen.
Thanks a lot Giulio.
Is turbulence also available in this version ?

Is this AeroFoam for v1.5 available yet ? If so where can I find it ?

rhoCentralFoam
 

There is this code rhoCentralFoam, how does it differ from the newer AeroFoam?

AeroFoam for OpenFoam 1.6
 

Hi all,

Is AeroFoam available for OpenFoam 1.6. If so where can i download the source code.

thanks
mecbe

Aerofoam for OpenFoam 1.6
 

hello everybody,

i just wanted to know if AeroFoam is available for OF 1.6 .
thank you

Unfortunately No.

Thanks
mecbe

Hello,

Is it there any version of Aerofoam for OpenFoam 1.7.1?

Regards,

Hi all, i would like to use the 3D toolbox to make a movie of the AGARD wing of the tutorial. I have it activated, but i do not know how to use it.

Can someone tell me please?

Thanks a lot!

Marta

Dear all, i am trying to run the AGARD wing tutorial, but I do not know how to make the movie i saw on the website.
I have already activated the 3D toolbox and data are saved correctly, so everything is fine with the run.

Thanks a lot,

Marta

Problem regarding adding another monoflux scheme in Aerofoam
 

Hi Guilo,

I am trying to add a new monoflux scheme (KFVS) similar to that of AUSM, Roe fluxes, etc.,. I have made correction to include the file in header for Aerofoam.C (#include KFVSFlux.C) file and called successfully from build_fluxes.C, I got no compilation errors when I ran wmake in the AeroFoam folder through Terminal.

But when I use the scheme to run a test case it gives the following error.

Does everything look fine? [y] to continue [n] otherwise y

Starting time loop...

#0 Foam::error:: printStack(Foam::Ostream&) in "/home/xx/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/xx/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0x997400]
#3 kfvs_flux(double, double, double, Foam::Vector<double>, Foam::Vector<double>, double, double, Foam::Vector<double>, double*, Foam::Vector<double>*, double*, double) in "/home/madhu/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/AeroFoam"
#4 build_fluxes(Foam::Field<double>*, Foam::word, Foam::dimensioned<double>*, Foam::dimensioned<double>*, Foam::dimensioned<double>*, Foam::fvMesh*, Foam::Field<int>*, Foam::Field<int>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>*, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>* at Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>*, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>*, Foam::word, Foam::word, Foam::word, Foam::word, double)
#5 main in "/home/xx/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/AeroFoam"
#6 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#7 _start at /build/buildd/eglibc-2.10.1/csu/../sysdeps/i386/elf/start.S:122
Floating point exception


Do I need to make changes in other places to incorporate the file.

Thanks in advance,
Eeshwar.

AeroFoam
 

run open foam
 

hi
I new user open foam and can generate mesh but I can't run it
please help me.

Hi Iman,
Explain how you generated mesh in OpenFOAM. If you generated with Mesh Conversion utilities , please go through the following link http://www.afs.enea.it/funel/CFD/Ope...shconvtest.htm and user guide .

dear Rock.Senthil Kumar
Its good
thanks for help you:)

new version of AeroFOAM for OpenFOAM 1.7.1
 

Attachment 6360

Attachment 6361

enjoy~! this is not my contribution, but I hope it is helpful!:cool:

This is really helpful and I hope work continues in this area of OpenFOAM.

Thanks Solarberiden.

Is it not your contribution? I am sure that you contributed posting here...:)
Do you know the author of the codes? Or a URL where it is located?

Problem with Aerofoam for OF1.7.1
 

Hello everybody,

I'm trying to use the solver AeroFoam with the 1.7.1 version of OF; it's part of an important project for my study. I downloaded the Aerofoam files at :
https://sourceforge.net/projects/aerofoam/files/Source/
And now I want to test this solver with a tutorial case : the oblique shock.
However, (after 'BlockMesh'), when I launch AeroFoam, I have :

-------------------------------------------------------------

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

Create mesh for time = 0

========================================
Reading Thermodynamic Properties...
========================================
gamma [-] = 1.4
Pr [-] = 0.72
R [J/KgK] = 287.05
Cp [J/KgK] = 717.625
Cv [J/KgK] = 1004.67
----------------------------------------

========================================
Reading Initial Conditions...
========================================
Reading field p
Reading field T
Reading field U
----------------------------------------
Creating conservative variables at t^(k)
Creating conservative variables at t^(k-1)
Creating conservative variables at t^(k-2)
Creating conservative variables fluxes at t^(k)
Creating primitive variables gradient fields (optional)
----------------------------------------

========================================
Creating Mesh Connectivity...
========================================
Using adaptive algorithm...
id_LL_i7290SIZE3510Loop on Internal Faces...
i=0i=1i=2i=3i=4i=5i=6i=7i=8i=9i=10i=11i=1
until i=3509
Loop on Boundary Faces...
----------------------------------------

========================================
Mesh
========================================
# of Faces = 3510
# of Cells = 1800
# of BndPatch = 5
0) inlet @ top
1) supersonicInlet @ inlet
2) extrapolatedOutlet @ outlet
3) symmetryPlane @ bottom
4) empty @ frontAndBack
----------------------------------------

========================================
Numerical Scheme
========================================
flowType = Euler
timeScheme = RK2
CourantMax = 0.1
CourantSteps = 25
MonotoneFlux = Roe
HighResFlux = LW
fluxLimiter = VL
entropyFix = HH2
extrapolateBC = 1
residualNorm = L1
minResidual = -1
----------------------------------------
Does everything look fine? [y] to continue [n] otherwise

-------------------------------------------------------------
I type yes and get :

-------------------------------------------------------------
blf fx i=1
...
blf fx i=3509
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3
at Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>*, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>*, Foam::word, Foam::word, Foam::word, Foam::word, double)
#4
in "/home/emilie/OpenFOAM/bruno-1.7.1/applications/bin/linux64GccDPOpt/AeroFoam"
#5 __libc_start_main in "/lib/libc.so.6"
#6
in "/home/emilie/OpenFOAM/bruno-1.7.1/applications/bin/linux64GccDPOpt/AeroFoam"
Segmentation fault

-------------------------------------------------------------

Does anyone has an idea why I have this error ?:confused:

Thank you very much for your help ! :)

Emilie

is there any version can run in parallel？

could you please rebuild your openfoam in debug mode and repost the error message? The final error message you posted didn't bring forwards any location about the error position information, which is a little bit difficult for other "depanner" , LOL.
by the way, the author of the version you are using and downloading from sourceforge JLight2010 had mentioned that
http://www.cfd-online.com/Forums/ope...free-code.html
he is not sure about the code correctness, further discuss would be depending your reply, I have tried some examples including the OneraM6 case and YF-17 case.

Finally, till now it didn't happen, the parallel version