CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

MacCormack's Solution

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

LinkBack Thread Tools Search this Thread Display Modes
Old   April 13, 2016, 15:53
Default MacCormack's Solution
New Member
Join Date: Jun 2015
Posts: 25
Rep Power: 11
yow is on a distinguished road
For my undergrad thesis, I recently implemented the MacCormack's space marching method in Matlab for compressible, inviscid flow over two plate angles to output two shocks in an attempt to verify my fluent solutions. The plate angles are: A1 = 10 deg, and A2 = 11.4 deg. I used the method found in Anderson's book CFD - The Basics with Applications (1995) around pg. 400. I use Courant number of .001 to avoid instabilities and CFL# = .5.

Shock angle and flow properties between the first and second shock match fluent and 1D solutions. Unfortunately, for the second shock, any A2 > 10.2 deg. yields imaginary flow properties behind the shock. Interestingly, the real part of all the imaginary results match fluent and 1D solutions. The 2nd shock is even slightly curved due to the nonuniform flow field, just like FLUENT!

What's the underlying reason for these imaginary results and can I simply ignore the imaginary parts? Any help is appreciated. Thank you.

yow is offline   Reply With Quote


anderson, maccormack, plate, shock

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
Why 3D solid-pore geometry showing diverged solution? Sargam05 OpenFOAM 0 December 3, 2012 15:45
Analytic solution for 2D steady Euler equations jojo81 Main CFD Forum 0 October 15, 2012 12:05
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
solution singularity litonx OpenFOAM Running, Solving & CFD 1 February 21, 2007 01:32
Discussion about Mesh independant solution Seb Main CFD Forum 13 May 22, 2001 13:37

All times are GMT -4. The time now is 23:44.