For the Full Potential Equation of Transonic Flow*
MATHEMATICS OF COMPUTATION VOLUME 44, NUMBER 169 JANUARY 1985, PAGES 1-29 Entropy Condition Satisfying Approximations for the Full Potential Equation of Transonic Flow* By Stanley Osher, Mohamed Hafez and Woodrow Whitlow, Jr. Abstract. We shall present a new class of conservative difference approximations for the steady full potential equation. They are, in general, easier to program than the usual density biasing algorithms, and in fact, differ only slightly from them. We prove rigorously that these new schemes satisfy a new discrete "entropy inequality", which rules out expansion shocks, and that they have sharp, steady, discrete shocks. A key tool in our analysis is the construction of an "entropy inequality" for the full potential equation itself. We conclude by presenting results of some numerical experiments using our new schemes. I. Introduction. The full potential equation is a common model for describing supersonic and subsonic flow close to the speed of sound. The flow is assumed to be that of a perfect gas, and the assumptions of irrotationality and constant entropy are made. The resulting equation is a single nonlinear partial differential equation of second order, which changes type from hyperbolic to elliptic, as the flow goes from supersonic to subsonic. Flows with a supersonic component generally have solutions with shocks, so the conservation form of the equation is important. This formulation, (FP), is one of three conservative formulations used for inviscid transonic flows. The other two are transonic small-disturbance equation, (TSD), and Euler equation, (EU), which is the exact inviscid formulation. The FP formulation is the most efficient of the three in terms of accuracy-to-cost ratio for a wide range of inviscid transonic flow applications for real geometries.
[Show full text]