Currently, the IMOC (Interactive Method of Characteristics) program developed for the University of Queensland in1991 by Jacobs and Gourlay provides solutions only for isentropic, compressible flows. Development of the (MOC) equations required to obtain a solution for 2D non-isentropic flows with shocks is described in this thesis. Starting with the equations of motion for rotational flow from Liepmann and Roshko (1966), a characteristic equations are determined by relating them to Cartesian co-ordinates such that integration over a finite characteristic mesh is possible. Using these equations, unit processes have been developed for an interior point calculation, a wall point, as well as a point along the shock. Combinations of the three unit processes have been implemented in MATLAB to obtain a solution for a test case of a sharp nosed body describe in Zucrow & Hoffman (1977).
Numerical methods required for computer implementation of the equations were identified to incorporate in the MATLAB code. Testing of the MATLAB code against the test case indicated a reasonable solution and verified the equations and unit processes used. These results indicated that there is an error in the convergence of the solution with mesh refinement. Prior to implementation of the equations and processes outlined in this thesis into IMOC, it is recommended that further testing be done to determine the affect of this error, and determine a solution.
To enhance the capabilities of the program, it is proposed that the equations be manipulated to be valid for axi-symmetric flow. Furthermore, the processes describe in this thesis assume a uniform flow before the shock wave, it would broaden the programs applications if this could be manipulated to be valid for non-uniformed pre-shock flow.