CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

hypersonic cfd free code???

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

Like Tree1Likes
  • 1 Post By Solarberiden

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 8, 2010, 08:17
Default hypersonic cfd free code???
  #1
Member
 
inginer's Avatar
 
Ovidiu Michiu
Join Date: Apr 2010
Location: Munich, Germany
Posts: 53
Rep Power: 16
inginer is on a distinguished road
hello forum,

as the title says my question is the following:

does anybody know a hypersonic cfd FREE code?

best,
Ovidiu
inginer is offline   Reply With Quote

Old   December 9, 2010, 13:25
Default
  #2
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 16
nakul is on a distinguished road
Hi,

At present no solver exists for hypersonic flow in OpenFOAM. But you can always make your own. I was thinking of doing it but just couldn't find time.

If you want you may begin working on it and I could help.
nakul is offline   Reply With Quote

Old   December 10, 2010, 08:45
Default
  #3
New Member
 
Nishit Joseph
Join Date: Nov 2010
Posts: 29
Rep Power: 16
JLight is on a distinguished road
Well...none part of the official package ... there is one called aeroFoam that may be of interest to you

http://www.aero.polimi.it/freecase/?...namic_problems

This post may also be of interest (the above link is taken from here)

http://www.cfd-online.com/Forums/ope...-aerofoam.html

But I have not tested it (yet!). I am in the process of developing a solver myself. Hypersonic CFD code is a bit too general ... what type of flow conditions are you interested in.
JLight is offline   Reply With Quote

Old   December 10, 2010, 13:46
Default density based solver
  #4
Member
 
MSR CHANDRA MURTHY
Join Date: Mar 2009
Posts: 33
Rep Power: 17
chandramurthy is on a distinguished road
Hypersonic flows are characterized by strong energy-flow copuling, requiring a coupled solver (density based). Currenly this solver is not avilable in OF. A brief intriduction about this technology is highlighed by Prof. H jasak in the recent opensourcecfd conferece 2010.http://www.opensourcecfd.com/conference2010/proceedings/documents/Session_III/OSCIC2010_Rusche_Jasak.pdf.

alternatively, you can try open source density based code http://www.freecfd.com/
chandramurthy is offline   Reply With Quote

Old   December 16, 2010, 08:44
Default
  #5
Member
 
inginer's Avatar
 
Ovidiu Michiu
Join Date: Apr 2010
Location: Munich, Germany
Posts: 53
Rep Power: 16
inginer is on a distinguished road
hello JLight,

I have only the Mach number, which is 10 and the vehicle is flying on the layer of the stratosphere of the Earth.

thk for aeroFoam.
I have a small problem. Is not working on OF 1.7.0. My supposition is this version doesn´t have this C++ library: Cloud.H and passiveParticle.H

best,
Ovidiu
inginer is offline   Reply With Quote

Old   December 16, 2010, 09:51
Default AeroFoam for OpenFOAM-1.7.x
  #6
New Member
 
Nishit Joseph
Join Date: Nov 2010
Posts: 29
Rep Power: 16
JLight is on a distinguished road
hi Ovidiu
OpenFoam does have the library. You need to edit the file in the Make directory. I have uploaded the necessary file in Make.zip

That was the easy part , aeroFoam still has many errors . I do not have the time to go through it to figure out the logic and correct the code as of now (assuming i am able too ), but i am open to the possibility in the coming week.

Hopefully you will have better luck than me!

Cheers!

JLight

PS: I also uploaded the log file for building AeroFoam to give you an idea of the problems encountered.
Attached Files
File Type: zip Make.zip (468 Bytes, 56 views)
File Type: txt log.txt (8.8 KB, 64 views)
JLight is offline   Reply With Quote

Old   December 16, 2010, 10:11
Default
  #7
Member
 
MSR CHANDRA MURTHY
Join Date: Mar 2009
Posts: 33
Rep Power: 17
chandramurthy is on a distinguished road
Quote:
Originally Posted by inginer View Post

I have only the Mach number, which is 10 and the vehicle is flying on the layer of the stratosphere of the Earth.

Ovidiu
U check the order of Knudsen number, is continuum mechanics valid at such a rarefied atmosphere?
chandramurthy is offline   Reply With Quote

Old   December 25, 2010, 10:56
Default
  #8
New Member
 
rock.senthilkumar's Avatar
 
senthilkumar.R
Join Date: Sep 2010
Location: Ranchi, India
Posts: 20
Rep Power: 16
rock.senthilkumar is on a distinguished road
Hi JLight,
If u completed updation of AeroFoam for OF 1.6 + , please reply with details.
__________________
Senthil Kumar. R
Department of space Engineering and Rocketry,
Birla Institute of Technology,Mesra.
rock.senthilkumar@gmail.com
rock.senthilkumar is offline   Reply With Quote

Old   December 28, 2010, 06:18
Default
  #9
Senior Member
 
Ivan Flaminio Cozza
Join Date: Mar 2009
Location: Torino, Piemonte, Italia
Posts: 210
Rep Power: 18
ivan_cozza is on a distinguished road
Send a message via MSN to ivan_cozza
I suggest you to contact Giulio Romanelli and Elisa Serioli at Politecnico di Milano, the authors of aeroFoam, maybe they could help you. The are very kind guys!

Ivan


Quote:
Originally Posted by inginer View Post
hello JLight,

I have only the Mach number, which is 10 and the vehicle is flying on the layer of the stratosphere of the Earth.

thk for aeroFoam.
I have a small problem. Is not working on OF 1.7.0. My supposition is this version doesn´t have this C++ library: Cloud.H and passiveParticle.H

best,
Ovidiu
ivan_cozza is offline   Reply With Quote

Old   December 28, 2010, 06:43
Default
  #10
New Member
 
rock.senthilkumar's Avatar
 
senthilkumar.R
Join Date: Sep 2010
Location: Ranchi, India
Posts: 20
Rep Power: 16
rock.senthilkumar is on a distinguished road
Hi Ivan,
Thank you for your suggestion.Surely I will contact them.
__________________
Senthil Kumar. R
Department of space Engineering and Rocketry,
Birla Institute of Technology,Mesra.
rock.senthilkumar@gmail.com
rock.senthilkumar is offline   Reply With Quote

Old   March 19, 2011, 17:42
Default Hypersonic CFD code
  #11
New Member
 
Nishit Joseph
Join Date: Nov 2010
Posts: 29
Rep Power: 16
JLight is on a distinguished road
Hi All
Sorry I had gone MIA. I had some personal trouble that are now sorted.

Regarding AeroFoam, I did not get a lot of opportunity to fix it. However, in my attempts I was able to remove all the compiling errors ( to the best I knew was right) by updating the codes.

But I am not sure of a certain section which is the connectivity . In the file "CreateConnectivity.C" @ about line 110

Code:
forAll( mesh.boundaryMesh(), iPatch )
{
    forAll( mesh.boundaryMesh()[iPatch].faceAreas(), ii )
    {
        
        // Input
        id_L = mesh.boundaryMesh()[iPatch].faceCells()[ii];           
        i    = mag(mesh.boundaryMesh()[iPatch].whichFace(id_L));
                 
        // *** ADAPTIVE VERSION *** ( works better with quad/hexa cells )
        if ( adaptiveConnectivitySearch == 1 )
the whichFace() will return -ve indexes from the way it is currently defined. Any ideas why whichFace does that?

More importantly I am not sure what the Authors are trying to do here so I was not sure if it needs 'fixing'.

This however is not reason I get runtime errors! The problem arises in the file "BuldFluxes.C" @ the lines about 57
Code:
    forAll( mesh->Sf(), i )
    {
        // Input 
        
        id_L        = (*mesh).faceOwner()[i];
        id_R        = (*mesh).faceNeighbour()[i];
Not just the above, but all loops that use the mesh->.Sf() . What happens is the loops seems to step over the arrays by 1 all the time (causing an error). I tried working this out but my debugger takes like a hour to step every single time into the code .

Also I would appreciated if someone tried this on a give known working example. I have not done this, may be the results may be different!

Any help will be appreciated.

The files are located here
http://sourceforge.net/projects/aerofoam/

Nishit
JLight is offline   Reply With Quote

Old   June 17, 2011, 04:42
Default
  #12
New Member
 
TCH
Join Date: Jul 2010
Location: Beijing City
Posts: 15
Rep Power: 16
Solarberiden is on a distinguished road
JLight, thanks for you major contribution about the correction and update the AeroFoam 1.7.1 version, I tried some cases in the aeroFoam 1.7.1 and I feel that, the code in CreateConeectivity.C line 110 should be changed form

i = mag(mesh.boundaryMesh()[iPatch].whichFace(id_L));

into

i = mag(mesh.boundaryMesh()[iPatch].start() + ii);

during running the case of OneraM6 and YF-17, the storage of connectivity file routing should be changed from

./Data
./Log

to

../Data
../Log

further feedback about this post is welcome. :-]
aerothermal likes this.
Solarberiden is offline   Reply With Quote

Old   January 18, 2012, 09:16
Default AeroFoam 1.7.1
  #13
Member
 
Join Date: Nov 2011
Posts: 32
Rep Power: 15
Andy_bm is on a distinguished road
Hi.
I have downloaded AeroFoam 1.7.1 from http://sourceforge.net/projects/aerofoam/ and has compiled it.Errors at compilation was not found, but when I have begun test AeroFoam on an example OneraM6,ObliqueShock it doesn't work and gives out the following error:
for ObliqueShock log .o file:
Code:
Create time
Create mesh for time = 0
// using new solver syntax:
cellDisplacement
{
solver PCG;
preconditioner DIC;
tolerance 1e-08;
relTol 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
This is a start 10 % 
19 % 
This is a start 10 % 
19 % 
29 % 
This is a start 10 % 
19 % 
29 % 
This is a start 10 % 
19 % 
29 % 
This is a start 10 % 
19 % 
29 % 
This is a start 9 % 
19 % 
28 % 
This is a start 9 % 
19 % 
28 % 
----------------------------------------
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_i216SIZE84Loop on Internal Faces...
i=0i=1i=2i=3i=4i=5i=6i=7 10 % 
i=8i=9i=10i=11i=12i=13i=14i=15 19 % 
29 % 
39 % 
48 % 
58 % 
67 % 
77 % 
87 % 
96 % 
9 % 
19 % 
29 % 
39 % 
49 % 
59 % 
69 % 
79 % 
89 % 
99 % 
38 % 
47 % 
56 % 
66 % 
75 % 
85 % 
94 % 
9 % 
19 % 
29 % 
39 % 
49 % 
59 % 
69 % 
79 % 
89 % 
99 % 
i=16i=17i=18i=19i=20i=21i=22i=23 29 % 
i=24i=25i=26i=27i=28i=29i=30i=31 38 % 
i=32i=33i=34i=35i=36i=37i=38i=39 48 % 
i=40i=41i=42i=43i=44i=45i=46i=47 57 % 
i=48i=49i=50i=51i=52i=53i=54i=55 67 % 
i=56i=57i=58i=59i=60i=61i=62i=63 76 % 
i=64i=65i=66i=67i=68i=69i=70i=71 86 % 
i=72i=73i=74i=75i=76i=77i=78i=79 95 % 
i=80i=81i=82i=83----------------------------------------
Loop on Boundary Faces...
9 % 
19 % 
29 % 
39 % 
49 % 
59 % 
69 % 
79 % 
89 % 
98 % 
----------------------------------------
========================================
Mesh 
39 % 
48 % 
58 % 
67 % 
77 % 
87 % 
96 % 
9 % 
18 % 
28 % 
37 % 
46 % 
56 % 
65 % 
75 % 
84 % 
94 % 
38 % 
48 % 
57 % 
67 % 
76 % 
86 % 
95 % 
9 % 
19 % 
29 % 
39 % 
49 % 
59 % 
69 % 
79 % 
89 % 
98 % 
38 % 
48 % 
57 % 
67 % 
76 % 
86 % 
95 % 
9 % 
19 % 
29 % 
39 % 
49 % 
59 % 
69 % 
79 % 
89 % 
98 % 
38 % 
48 % 
57 % 
67 % 
76 % 
86 % 
95 % 
9 % 
19 % 
29 % 
39 % 
49 % 
59 % 
69 % 
79 % 
89 % 
98 % 
38 % 
47 % 
56 % 
66 % 
75 % 
85 % 
94 % 
9 % 
19 % 
29 % 
39 % 
49 % 
59 % 
69 % 
79 % 
89 % 
99 % 
========================================
# of Faces = 84
# of Cells = 50
# of BndPatch = 8
0) supersonicInlet @ inlet
1) extrapolatedOutlet @ outlet
2) symmetryPlane @ bottom
3) supersonicInlet @ top
4) empty @ defaultFaces
5) @ procBoundary0to1
6) @ procBoundary0to2
7) @ 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
----------------------------------------
Does everything look fine? [y] to continue [n] otherwise 
Starting time loop...
SIZE84blf fx i=0
blf fx i=1
...........
blf fx i=83
And further the case isn't considered and gives out an error
(Foam::error:rintStack)
For other cases there are same problems.

Thanks.
Andy_bm is offline   Reply With Quote

Old   January 18, 2012, 09:44
Default
  #14
New Member
 
TCH
Join Date: Jul 2010
Location: Beijing City
Posts: 15
Rep Power: 16
Solarberiden is on a distinguished road
Actually I have poured a lot of time on aeroFoam for promote more time integration feature for the code but I could not manage to build more feasible and robust flux schemes on high mach number backward facing step (strong expansion wave problems), I am willing to offer my code but I don't know how. And I would like to provoke some others on proceeding the more powerful aerofoam features. Any body interested of sharing my code?
Solarberiden is offline   Reply With Quote

Old   January 18, 2012, 09:53
Default
  #15
Member
 
Join Date: Nov 2011
Posts: 32
Rep Power: 15
Andy_bm is on a distinguished road
I am interested. what its difference from code Aerofoam, i don't understand? What version openfoam it is constructed? Whether you know how to rectify my error?
Andy_bm is offline   Reply With Quote

Old   January 18, 2012, 10:06
Default
  #16
New Member
 
TCH
Join Date: Jul 2010
Location: Beijing City
Posts: 15
Rep Power: 16
Solarberiden is on a distinguished road
On what version the code was constructed didn't matter, I am just using the code on openfoam 2.0 and it is working well. Questions:
1 - What test case are you just running?
2 - Have you ever tried other test cases aerofoam?
3 - you should first locate where is the "blf fx i=1
...........
blf fx i=83" instructions, it looks strange for me and it does not appear to be right for me.
could you please let me know where you locate the "blf fx i = xxx" instructions.
:-]
Solarberiden is offline   Reply With Quote

Old   January 18, 2012, 10:26
Default
  #17
Member
 
Join Date: Nov 2011
Posts: 32
Rep Power: 15
Andy_bm is on a distinguished road
I tested Aerofoam on obliqueshock case.I don't know how this solver worked,and don't know that means "blf fx i=1" in the log.

I work over a problem stream injection in vacuum. Start Mach numder > 3..4, at present I use rhoCentralFoam on OpenFOAM v2.0.1,heard that aerofoam gives out very good results, and i wanted to compare it.
Andy_bm is offline   Reply With Quote

Old   January 18, 2012, 11:07
Default
  #18
New Member
 
TCH
Join Date: Jul 2010
Location: Beijing City
Posts: 15
Rep Power: 16
Solarberiden is on a distinguished road
please leave your email address and I will send you a copy of the code if I've got the time. Currently, I am not on my working post. I guess you need to wait for quite a while. I will get you informed. :-]
Solarberiden is offline   Reply With Quote

Old   May 31, 2013, 09:49
Default Aerofoam
  #19
Senior Member
 
Nishant
Join Date: Mar 2009
Location: Glasgow, UK
Posts: 166
Rep Power: 17
nishant_hull is on a distinguished road
Hi Solarbiden,

I would be interested in helping to debug the case. I have previously done some development on openfoam. However I am new to Aerofoam and I was attracted towards it for its applicability in supersonic cases. I was also wondering if it can be used in aeroelastic applications? My email is nishant.nitj@gmail.com for any communication.

Thanks and regards,

Nishant
__________________
Thanks and regards,

Nishant
nishant_hull is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
where can i find such free CFD code? euler Main CFD Forum 9 October 14, 2018 07:17
CFD Design...The CFD Future John C. Chien Main CFD Forum 20 November 20, 2015 00:40
Register CFD events and get a free CFD book Jonas Larsson Main CFD Forum 8 April 1, 2005 15:27
CFD code structure (F90) ma Main CFD Forum 4 January 10, 2005 21:47
Looking for a free CFD code a. Main CFD Forum 2 September 22, 2004 04:13


All times are GMT -4. The time now is 06:30.