Optimum plane diffusers in laminar flow

The problem of determining the profile of a plane diffuser (of given upstream width and length) that provides the maximum static pressure rise is solved. Two-dimensional, incompressible, laminar flow governed by the steady-state Navier-Stokes equations is assumed through the diffuser. Recent advances in computational resources and algorithms have made it possible to solve the ‘direct’ problem of determining such a flow through a body of known geometry. In this paper, a set of ‘adjoint’ equations is obtained, the solution to which permits the calculation of the direction and relative magnitude of change in the diffuser profile that leads to a higher pressure rise. The direct as well as the adjoint set of partial differential equations are obtained for Dirichlet-type inflow and outflow conditions. Repeatedly modifying the diffuser geometry with each solution to these two sets of equations with the above boundary conditions would in principle lead to a diffuser with the maximum static pressure rise, also called the optimum diffuser. The optimality condition, that the shear stress all along the wall must vanish for the optimum diffuser, is also recovered from the analysis. It is postulated that the adjoint set of equations continues to hold even if the computationally inconvenient Dirichlet-type outflow boundary condition is replaced by Neumann-type conditions. It is shown that numerical solutions obtained in this fashion do satisfy the optimality condition.

[1]  A. Rizzi,et al.  Computation of flow around wings based on the Euler equations , 1984, Journal of Fluid Mechanics.

[2]  V. N. Vatsa Accurate solutions for transonic viscous flow over finite wings , 1986 .

[3]  G. Volpe,et al.  The design of transonic aerofoils by a well‐posed inverse method , 1986 .

[4]  P. R. Garabedian,et al.  NUMERICAL DESIGN OF TRANSONIC AIRFOILS**The work presented in this paper is supported by the AEC Computing and Applied Mathematics Center, Courant Institute of Mathematical Sciences, New York University, under Contract AT(30-1)-1480 with the U. S. Atomic Energy Commission. , 1971 .

[5]  O. Pironneau On optimum design in fluid mechanics , 1974 .

[6]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[7]  Ephraim M Sparrow,et al.  Experiments on hydrodynamically developing flow in rectangular ducts , 1970 .

[8]  Towards the Minimum Drag on a Body of Given Volume in Slow Viscous Flow , 1971 .

[9]  T. Tranen,et al.  A rapid computer aided transonic airfoil design method , 1974 .

[10]  R. M. Hicks,et al.  Wing Design by Numerical Optimization , 1977 .

[11]  R. J. Goldstein,et al.  Measurement of Laminar Flow Development in a Square Duct Using a Laser-Doppler Flowmeter , 1967 .

[12]  G. McFadden,et al.  Computational Fluid Dynamics of Airfoils and Wings , 1982 .

[13]  P. Henne,et al.  An inverse transonic wing design method , 1980 .

[14]  M. Drela NEWTON SOLUTION OF DIRECT AND INVERSE TRANSONIC , 1985 .

[15]  Antony Jameson,et al.  Aerodynamic design via control theory , 1988, J. Sci. Comput..

[16]  Arthur Rizzi,et al.  Computation of inviscid incompressible flow with rotation , 1985, Journal of Fluid Mechanics.

[17]  M. Giles,et al.  Newton solution of direct and inverse transonic Euler equations , 1985 .

[18]  H. Çabuk,et al.  Shape optimization analysis: first- and second-order necessary conditions , 1990 .

[19]  O. Pironneau On optimum profiles in Stokes flow , 1973, Journal of Fluid Mechanics.

[20]  Roland Glowinski,et al.  On the numerical computation of the minimum-drag profile in laminar flow , 1975, Journal of Fluid Mechanics.

[21]  E. Turkel,et al.  A multistage time-stepping scheme for the Navier-Stokes equations , 1985 .

[22]  Joe F. Thompson,et al.  Automatic numerical generation of body-fitted curvilinear coordinate system for field containing any number of arbitrary two-dimensional bodies , 1974 .