# Numerical Solution of the Blunt Body in Supersonic Flow Problem

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

 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 Sponsored Links 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 Navier-Stokes? 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. gomnam15010 likes this.

 February 16, 2012, 11:28 #2 New Member   Join Date: Feb 2012 Posts: 13 Rep Power: 7 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/cgi-bin/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 time-dependent technique. I am also trying to solve this problem. I am using the MacCormarck's predictor-corrector time marching technique. gomnam15010 and aerosjc like this.

 June 5, 2012, 12:54 #3 New Member   Kathrin Join Date: Jun 2012 Posts: 9 Rep Power: 7 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 gomnam15010 likes this.

 June 5, 2012, 15:33 #4 New Member   Join Date: Feb 2012 Posts: 13 Rep Power: 7 Hi Kathrin, I haven't solved the blunt body problem either, however I've solved the Prandtl-Meyer 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...on-wave-solver Hope this helps. gomnam15010 likes this.

 June 5, 2012, 15:57 #5 New Member   Kathrin Join Date: Jun 2012 Posts: 9 Rep Power: 7 Hi Ivan, coincidentally I've seen your Prandtl-code 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: McGraw-Hill, 1989.) But I hoped that someone already accomplished the EULER equation as a MATLAB-Code 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: http://www.cfd-online.com/Forums/mai...swatitsch.html ?? Best regard Kathrin

 June 5, 2012, 18:09 #6 New Member   Join Date: Feb 2012 Posts: 13 Rep Power: 7 Hi Kathrin, Yes, the Prandtl-Meyer 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 Mach-number-independence 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: 7 Hi Ivan, OK, that explains why the Prandtl-Code 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: 7 Thank you very much for the documents Kathrin, I have already downloaded them. As I said, as soon as I finish the course (3-4 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: 7 Hello Ivan, I just started to develop the space-marching technique for EULER in the purely supersonic case. Have you already made any progesses with the time-dependent 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: 7 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: 7 Hello Ivan, I wrote the code for a steady (!) EULER-case and used the space-marching MacCormack-method 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: 7 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: 7
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
Attached Files
 EULER steady hypersonic.zip (20.5 KB, 43 views)

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: 7 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: 8 Hi Kathrin I want to calculate the aerodynamic forces for re-entry 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: 7 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 Euler-code 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: 8 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: 4 Hey, I am trying to implement a Mac Cormack technique to solve the blunt body problem using finite-differences 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: 3
Quote:
 Originally Posted by andres17 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 Navier-Stokes? 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.
Hello, I am interested to solve blunt body problem with the help of fdm. I am trying to obtain conditions at initial time for this problem. Please suggest me that how to find these conditions at initial time.

 April 26, 2016, 15:29 #20 New Member   Join Date: Dec 2014 Posts: 3 Rep Power: 4 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 time-step and insert artificial viscosity <- This is NOT recommended. Good luck!

 Tags blunt body, cfd, matlab, supersonic flow

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post BigFrank FLUENT 6 April 24, 2014 15:13 Ruben Main CFD Forum 43 May 7, 2011 03:32 jan CFX 1 July 31, 2007 19:44 diaw Main CFD Forum 104 February 16, 2006 06:44 Abhijit Tilak Main CFD Forum 6 February 5, 1999 02:16