
[Sponsors] 
September 15, 2010, 15:43 
Numerical Solution of the Blunt Body in Supersonic Flow Problem

#1 
New Member
Noel Andres
Join Date: Sep 2010
Posts: 1
Rep Power: 0 
Hi,
I want to do a CFD code to solve the blunt body problem in supersonic flow. I'll probably pick a sphere in supersonic flow. I want to find the shock wave location and the static pressure at the surface of the sphere. This will be a 2D problem. Are the Euler equations enough for this, or do I need to go NavierStokes? Also, since I don't know beforehand the shape and location of the bow shock wave, the code has to find it itself. Will the Euler equations handle that? I will be using the finite difference approach. I'll be coding this in Matlab since it is the language I know. I will post the code in the Matlab Central page. Any help with the questions I have are appreciated. Thanks. 

February 16, 2012, 11:28 

#2 
New Member
Join Date: Feb 2012
Posts: 13
Rep Power: 6 
Hi, since you seem to be solving an inviscid flow problem, the Euler equations should be enough. Take a look at this document: http://www.dtic.mil/cgibin/GetTRDoc...f&AD=AD0846913. There you will find the proper equations and transformation.
The best way to handle the detached shock wave in the finite difference method is to use the timedependent technique. I am also trying to solve this problem. I am using the MacCormarck's predictorcorrector time marching technique. 

June 5, 2012, 12:54 

#3 
New Member
Kathrin
Join Date: Jun 2012
Posts: 9
Rep Power: 6 
Hi andres17,
I'm also very interested in your program with the EULER equations in Matlab. Is it possible to have a look at the code?? Best regards Kathrin 

June 5, 2012, 15:33 

#4 
New Member
Join Date: Feb 2012
Posts: 13
Rep Power: 6 
Hi Kathrin,
I haven't solved the blunt body problem either, however I've solved the PrandtlMeyer expansion wave problem form Anderson's book (chapter 8). It is based on Euler equations. You can find it here: http://www.mathworks.es/matlabcentra...onwavesolver Hope this helps. 

June 5, 2012, 15:57 

#5 
New Member
Kathrin
Join Date: Jun 2012
Posts: 9
Rep Power: 6 
Hi Ivan,
coincidentally I've seen your Prandtlcode working, yesterday. But for large angles (i. e. 20°) the iteration just stopped at about 250 iterations. Do you know why? Nevertheless, for comparisons I need the mentoined algorithm from andres17. Of course, I can do it my own with the deteiled desciption in chapter 5 (Anderson, John David: Hypersonic and High Temperature Gas Dynamics. New York: McGrawHill, 1989.) But I hoped that someone already accomplished the EULER equation as a MATLABCode and shares it... By the way, I've an other problem, where I don't find a solution  even in the case of a simple steamline. Can you help me here: Machnumberindependence Oswatitsch ?? Best regard Kathrin 

June 5, 2012, 18:09 

#6 
New Member
Join Date: Feb 2012
Posts: 13
Rep Power: 6 
Hi Kathrin,
Yes, the PrandtlMeyer code is very sensitive to angle variations. It is a very primitive implementation and does not work as a general solver. In fact, the computation blows up most of the times when changing the parameters, as you may have seen. It is my first CFD code, but is very good way to understand the physics of Euler equations and the MacCormack's predictor corrector technique. Anyway, at the end of chapter 8 of Anderson's CFD book (pag. 410), he says that when the leading egde of the expansion wave exits the computational domain through the upper boundary, and not through the right boundary, the boundary conditions of the upper boundary must be changed. That's why it doesn't work for large angles, because the boundary condition I set was simply the specification of uniform conditions equal to those in the upstream uniform flow. Maybe in a future I will try to improve the code. Regarding the blunt body problem, I have read the explanation in chapter 5 of Hypersonic and... but I still have some doubts about how to start the problem. I will take a look at it more closely this summer. I have read your post about Machnumberindependence Oswatitsch but I don't know Crocco's theorem and I have never implemented a time dependent MacCormack method, so, by now, I'm not able to help you with it. But as I said before, I will try to solve the blunt body in a few days. If it works I will inform you. Just a last thing, do you know where could I find this document? A time dependent computational method for blunt body flows, AIAA Journal, Vol 4, No. 12. Dec, 1966. Moretti, G. and Abbett, M. I think it describes all the necessary points to solve the problem. Best regards Ivan 

June 6, 2012, 02:05 

#7 
New Member
Kathrin
Join Date: Jun 2012
Posts: 9
Rep Power: 6 
Hi Ivan,
OK, that explains why the PrandtlCode only holds for certain angles. But the code looks really good. To your question: I found this paper from Abbett concerning the algothm https://www.unibw.de/mb/fakultaet/we...download/down2 and that concerning different boundary condions at the wall https://www.unibw.de/mb/fakultaet/we...download/down1 I can leave these files uploaded for about 24 hours. Could you please tell me when you got them?? You're right, I think the desciption in the paper, you were looking for, is quite helpful. As you said, it would be very nice if you would tell me, when you got a solution for this problem. Greetings Kathrin 

June 6, 2012, 05:12 

#8 
New Member
Join Date: Feb 2012
Posts: 13
Rep Power: 6 
Thank you very much for the documents Kathrin, I have already downloaded them.
As I said, as soon as I finish the course (34 weeks) I will start with the blunt body problem. I will tell you about it then. Best regards, Ivan 

August 15, 2012, 04:11 

#9 
New Member
Kathrin
Join Date: Jun 2012
Posts: 9
Rep Power: 6 
Hello Ivan,
I just started to develop the spacemarching technique for EULER in the purely supersonic case. Have you already made any progesses with the timedependent technique for sub and supersonic flows? Best regards Kathrin 

September 16, 2012, 11:44 

#10 
New Member
Join Date: Feb 2012
Posts: 13
Rep Power: 6 
Hi Kathrin,
I'm sorry for not answering until today, I've been on holiday and I couldn't send any message. I have made several progresses with the numerical solution of the blunt body problem. I've wirten a Matlab code followin both the documents we talked about and Ansderson's books. However, I have not been able to obtain a converged solution of the problem. After a huge number of time steps the calculation blows up. Anyway, the code works pretty well and I believe in a few days I will manage to run it until the end. Have obtained any solution as well? Greetings, Ivan 

September 25, 2012, 07:16 

#11 
New Member
Kathrin
Join Date: Jun 2012
Posts: 9
Rep Power: 6 
Hello Ivan,
I wrote the code for a steady (!) EULERcase and used the spacemarching MacCormackmethod from Anderson. However, I got nearly the same problem: the calculations blow up, too, and I was unable to find the mistake, although the code is rather simple. :( 

September 25, 2012, 13:00 

#12 
New Member
Join Date: Feb 2012
Posts: 13
Rep Power: 6 
Hi,
I am still working on the time marching solution, but I haven't managed to obtain a steady solution yet, however the results are improving little by little. Which space marching problem are you solving? The one in Anderson's hypersonic book? I can have a look at it if you want Best regards Iván 

September 25, 2012, 14:05 

#13 
New Member
Kathrin
Join Date: Jun 2012
Posts: 9
Rep Power: 6 
Hi Ivan,
that sounds good. How do I create an attachement so that only you can open it or which is only visible for a certain period of time? It's written on MATLAB R2012 A. Greetings Last edited by Kathrin; September 25, 2012 at 14:37. Reason: ... and here the attachement 

September 25, 2012, 16:44 

#14 
New Member
Join Date: Feb 2012
Posts: 13
Rep Power: 6 
OK, thanks, I will take a look at it and try to find out any error, however it will be difficult.
Good luck Iván 

November 24, 2013, 02:36 

#15 
New Member
Amir
Join Date: Mar 2011
Posts: 9
Rep Power: 7 
Hi Kathrin
I want to calculate the aerodynamic forces for reentry capsule, can i use your code "EULER steady hypersonic" for this purpose? If it is possible put your email address for me due future calls. Last edited by gomnam15010; November 24, 2013 at 02:37. Reason: it is mistake 

November 24, 2013, 08:06 

#16 
New Member
Kathrin
Join Date: Jun 2012
Posts: 9
Rep Power: 6 
Hi gomnam,
you can use this program as a very simple way of solving your problem but this code is just valid for completely supersonic fields. And it was one one of my first trials. So, there are some mistakes in it. By now, for my Phd. I developed a new and very fast Eulercode that works for both sub and supersonic regions behind the shock, including chemical reactions. Best regards 

November 24, 2013, 11:12 

#17 
New Member
Amir
Join Date: Mar 2011
Posts: 9
Rep Power: 7 
Hi Kathrin again
How is Arrangement of implementation subprograms in "EULER steady hypersonic" zipfile? and How can run these subprograms as one program in matlab? its 2D solver or 3D solver? thanks for your attention 

February 28, 2015, 16:20 

#18 
New Member
Join Date: Dec 2014
Posts: 3
Rep Power: 3 
Hey,
I am trying to implement a Mac Cormack technique to solve the blunt body problem using finitedifferences but it doesn't work for for than 30 iterations before temperature becomes negative somewhere. (specifically at the third node from the wall). Anyone who can give me ideas on what can/is causing this error? 

April 26, 2016, 12:59 

#19  
New Member
Iqrar Raza
Join Date: Mar 2016
Posts: 9
Rep Power: 2 
Quote:


April 26, 2016, 15:29 

#20 
New Member
Join Date: Dec 2014
Posts: 3
Rep Power: 3 
Take your initial condition to be the free stream value at the Mach number you are solving for. In essence you are inserting the body at a given instant in the free stream flow immediately.
Your Euler code should 'generate' the shock in front of the body provided you have the appropriate Upwind scheme and you are marching in time. < Recommended If you are using central differences, keep an eye on your timestep and insert artificial viscosity < This is NOT recommended. Good luck! 

Tags 
blunt body, cfd, matlab, supersonic flow 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
back pressure at exit for supersonic flow  BigFrank  FLUENT  6  April 24, 2014 15:13 
Steady state  Ruben  Main CFD Forum  43  May 7, 2011 03:32 
demo free flow blunt body in cfx ansys 11  jan  CFX  1  July 31, 2007 19:44 
Can 'shock waves' occur in viscous fluid flows?  diaw  Main CFD Forum  104  February 16, 2006 06:44 
Wall functions  Abhijit Tilak  Main CFD Forum  6  February 5, 1999 02:16 