Abstract:
A pressure-recovery device for a supersonic, continuous wave gas laser is capable of recovering the pressure of the supersonic fluid in an efficient manner and in a relatively light and compact assembly. The supersonic diffuser has a substantially rectangular inlet and top and bottom walls that extend from the inlet in the direction of the flow of gas. The top and bottom walls are separated by a height (l 1 ) at the inlet and are connected by two side walls. The side walls are separated by a width (h) at the inlet, wherein the ratio of height to width ranges at the inlet from about 0.001 to about 1. A compression region is formed by the top wall, bottom wall and side walls, wherein said top and bottom walls either are or are not symmetric about a central longitudinal axis of the compression region. The top and bottom walls may each in the shape of a Prandtl-Meyer streamline and the side walls are substantially planar, or the top and bottom walls are each in the shape of a lens analogy streamline and the side walls are substantially planar. A plurality of diffusers may be stacked in series with one another, each of the diffusers being substantially similar in construction.

Description:
CROSS-REFERENCE TO RELATED APPLICATION 
       [0001]    This application is a continuation-in-part of application Ser. No. 10/874,040, filed Jun. 22, 2004. 
     
    
     TECHNICAL FIELD 
       [0002]    This application relates to diffusers and, more particularly, to a design for an efficient supersonic diffuser. 
       BACKGROUND 
       [0003]    Supersonic diffusers have a variety of applications: wind tunnels, ground testing of high altitude rocket engines, engine inlet for supersonic aircraft, and any supersonic device that operates at a static pressure below ambient. One such device is the continuous wave (cw) supersonic chemical laser. While the present invention focuses on this last application, it may be applicable to other supersonic flows where a diffuser is required. 
         [0004]    A supersonic flow, where the pressure is below ambient, occurs in the laser cavity of a cw chemical laser. CW chemical lasers involve a steady, supersonic, low-pressure gaseous flow inside the laser cavity. In fact, all high-performance cw chemical lasers operate supersonically with a laser cavity pressure of a few Torr. For example, a chemical oxygen-iodine laser (COIL) typically operates in the 1 to 10 Torr range, although pressures as high as 20 Torr may be possible. A pressure value above about 4 Torr usually is achieved by adding diluent gas to the singlet oxygen generator (SOG) that drives the COIL device. A diffuser is then used to increase the device&#39;s exit pressure. If this exit pressure is still below ambient, the diffuser is followed by a pumping system that typically consists of mechanical pumps or an ejector system. In specialized cases, chemical pumping may be used. 
         [0005]    The type of high-performance laser under consideration might be mounted on a motorized vehicle, naval vessel, or on an aircraft. It requires a pressure recovery system to increase the pressure of the high-speed, spent laser effluent from its several Torr value to a pressure level slightly above ambient. At sea level, ambient pressure is about 760 Torr, while ambient pressure for an aircraft at a 40,000 foot altitude would be about 141 Torr. In any case, the size and weight of the pressure recovery system is of crucial importance for the viability of the overall laser system. It is thus important that the pressure recovery system be as compact and lightweight as possible for the intended application. 
         [0006]    Diffusers are common devices whose function is to convert as much as possible of the kinetic energy of a fluid, at the inlet of the diffuser, into an increased value for the pressure of the fluid at the device exit. The exit pressure, Pr, is referred to as the recovered pressure. Diffusers come in two categories, depending on whether the inlet flow is subsonic or supersonic. Supersonic diffusers are bulky and generally very inefficient, especially when their inlet Mach number is large. 
         [0007]    Diffusers used for chemical lasers generally consist of a converging supersonic section, followed by a throat region, which is then followed by a slowly diverging subsonic diffuser. The throat region and subsonic section can be lengthy. Most of the pressure increase occurs in the throat region, which is a duct containing a system of oblique shock waves. A subsonic diffuser, in order to avoid boundary-layer separation, increases its cross-sectional area gradually. This subsonic section only provides a modest amount of pressure recovery, typically less than 10% of the overall value of a supersonic/subsonic diffuser. For reasons of compactness, the subsonic portion of the diffuser may be attenuated or even bypassed. 
         [0008]    An oblique shock system generally starts at the diffuser&#39;s inlet and continues into the throat section. Most of the static pressure increase stems from the shock system and not from an isentropic process. By contrast, the stagnation pressure steadily decreases through the shock system. The overall decrease in the stagnation pressure of a diffuser is of crucial importance. At the exit, where the Mach number is small compared to unity, the recovered (static) pressure is essentially the stagnation pressure. A diffuser&#39;s efficiency is, therefore, just the ratio of the exit to inlet stagnation pressures. At a supersonic inlet, the stagnation pressure significantly exceeds the static pressure. For steady, isentropic flow, the stagnation pressure is a constant, and the diffuser&#39;s second-law efficiency is thus 100%. 
         [0009]    The standard pressure recovery system for any supersonic gas flow would be a supersonic/subsonic diffuser as is commonly used with a wind tunnel. Common to such diffusers is a poor efficiency because an oblique shock system is present inside the supersonic section, including the throat section, of the diffuser. The shock system decreases the stagnation pressure of the flow and causes boundary-layer separation. 
         [0010]    Use of a conventional supersonic/subsonic diffuser for a chemical laser has a number of major drawbacks, aside from its poor efficiency. They are bulky, heavy, and for a conventional COIL do not provide nearly enough pressure recovery, even for aircraft operation of a COIL at a 40,000 foot altitude. For a conventional COIL system with diluent in the SOG and a diffuser, an ejector system is still required for an aircraft-based system. Compared to a diffuser, an ejector system may be more compact but is heavier, much heavier, as the laser run time increases. 
         [0011]    Accordingly, there is a need for a supersonic diffuser able to provide greater pressure recovery than existing diffuser designs while remaining relatively compact and light in weight. 
       SUMMARY 
       [0012]    It is a general object of the present invention to provide a pressure-recovery device for a supersonic, continuous wave gas laser that is capable of recovering the pressure of the supersonic fluid in an efficient manner and in a relatively light and compact assembly. This and other objects of the present invention are achieved by providing a supersonic diffuser having a substantially rectangular inlet. Top and bottom walls extend from the inlet in the direction of the flow of gas. The top and bottom walls are separated by a height (l 1 ) at the inlet and are connected by two side walls. The side walls are separated by a width (h) at the inlet, wherein the ratio of height to width ranges at the inlet from about 0.001 to about 1. A compression region is formed by the top wall, bottom wall and side walls, wherein said top and bottom walls either are or are not symmetric about a central longitudinal axis of the compression region. 
         [0013]    According to the preferred embodiment of the present invention, the supersonic flow entering the diffuser is substantially free of intense pressure gradients, such as associated with shock waves. 
         [0014]    According to a preferred embodiment of the present invention, at least a portion of the surfaces of the top, bottom and side walls that contact flow of gas are treated to facilitate the maintenance of a turbulent boundary layer. 
         [0015]    According to a preferred embodiment of the present invention, the top and bottom walls are each in the shape of a Prandtl-Meyer streamline and the side walls are substantially planar. 
         [0016]    According to a preferred embodiment of the present invention, the top and bottom walls are each in the shape of a lens analogy streamline and the side walls are substantially planar. 
         [0017]    According to a preferred embodiment of the present invention, the diffuser is substantially isentropic in operation. 
         [0018]    According to a preferred embodiment of the present invention, the substantially rectangular inlet is in fluid communication with an optical chamber of a chemical laser to receive the nearly uniform supersonic gas flow. 
         [0019]    According to a preferred embodiment of the present invention, the chemical laser is a chemical oxygen-iodine laser (COIL). 
         [0020]    According to a preferred embodiment of the present invention, a duct having a rectangular cross section precedes the rectangular inlet, the duct for inducing turbulent boundary layer flow. 
         [0021]    According to a preferred embodiment of the present invention, a shock holder is disposed in a duct downstream of the compression region, the shock holder for positioning a weak normal shock in the gas flow. 
         [0022]    According to the preferred embodiment of the present invention, a plurality of diffusers are stacked in series with one another, each of the diffusers being substantially similar in construction. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0023]    For a more complete understanding of the present invention, and the advantages thereof, reference is now made to the following descriptions taken in conjunction with the accompanying drawings, in which: 
           [0024]      FIG. 1  is schematic diagram of a chemical laser system including a diffuser and a Brayton cycle pump. 
           [0025]      FIG. 2  is a cross section of a Prandtl-Meyer diffuser. 
           [0026]      FIG. 3  is a diagram showing the shape of the lower wall of a Prandtl-Meyer diffuser. 
           [0027]      FIG. 4  is a diagram showing the transonic region of a diffuser. 
           [0028]      FIG. 5  is a diagram showing the shape of the lower wall of a lens-analogy diffuser. 
           [0029]      FIG. 6  is a diagram schematically depicting the shape of the lower wall of two lens-analogy diffusers according to the present invention in a stacked (serial) relationship or orientation. 
       
    
    
     DETAILED DESCRIPTION 
       [0030]    Shown in  FIG. 1  is a chemical oxygen iodine laser (“COIL”) system. A singlet oxygen generator (“SOG”)  12  generates singlet delta oxygen for the laser. The singlet delta oxygen, along with a diluent (collectively “laser gases”) is forced from the SOG through the laser nozzle  14 . The laser gases flow through an optical cavity  16  at supersonic velocity and at a pressure typically between 1 and 20 Torr. After passing through the optical cavity  16 , the laser gases continue, preferably, to a diffuser  18  and a Brayton cycle pump (“BCP”)  22 . The preferred Brayton Cycle Pump is disclosed in U.S. patent application Ser. No. 10/874,039, filed concurrently herewith and entitled “Laser With Brayton Cycle Pump.” 
         [0031]    The diffuser  18  and BCP  22  each act to raise the pressure of the effluent stream so that, upon exiting the BCP  22 , the static pressure of the effluent will be greater than the ambient pressure and the effluent may be vented. The diffuser  18  is preferably designed according to the lens analogy or the Prandtl-Meyer formulation. 
         [0032]    An ejector, reciprocating pump or some equivalent pressure enhancing device may be used in place of the Brayton cycle pump. In some applications, the diffuser alone may increase the pressure above the ambient level, in which case no further pressure enhancement is necessary. The upstream flow entering the diffuser has a velocity {right arrow over (w)} 1 . The laser gas flow is generated by the laser nozzle  14  ( FIG. 1 ). The gas passes through a nearly constant cross-sectional area laser cavity  16  before entering the diffuser. The supersonic diffuser according to the present invention has two embodiments: one designed according to the Prandtl-Meyer (P-M) model and one according to the lens analogy (L-A) model. In both embodiments, the inlet is rectangular in cross-section, and the side walls of the diffuser are planar and parallel. In the P-M embodiment, the upper and lower or top and bottom walls are not symmetric with respect to each other: they are curved in accordance with a Prandtl-Meyer streamline of the flow through the diffuser. In the L-A embodiment, the top and bottom walls are symmetric with respect to one another and a central axis of the diffuser: each wall is curved in accordance with the lens analogy. 
         [0033]    The determination of the shapes of the sidewalls of both the P-M embodiment and the L-A embodiment are described below. Both embodiments strive to achieve substantially isentropic operation. 
         [0034]      FIG. 2  shows a diffuser  30  using a P-M compression. The diffuser is preferably a rectangular conduit with height l and width h. The width h is constant throughout the length of the diffuser—that is, the side walls are parallel. The height l varies along the length of the diffuser as shown in  FIG. 2  and described below. The height at the inlet  33  is designated l 1  and is generally in the 10 −3  h to h range. 
         [0035]    In order to obtain a nearly isentropic compression, with either a P-M or L-A diffuser, it is imperative that the supersonic flow entering the diffuser be substantially free of any intense pressure gradient, such as associated with shock waves. The preferred “Stable/Unstable Optical Cavity Resonator for Laser,” disclosed in U.S. patent application Ser. No. 10/874,064, filed Jun. 22, 2004, is thus warranted. In this disclosure, the supersonic flow inside the laser cavity is confined on all four sides where the confining walls are designed such that shock waves, or other intense pressure gradients, are not generated. 
         [0036]    There is preferably a short, roughened wall section  32  near the diffuser inlet  33  just before the beginning of the compression region  34  of the diffuser. The roughened wall section  32  can be used, if needed, to trip the boundary layer and helps to ensure a turbulent boundary layer in the diffuser. This approach dispenses with the need for boundary-layer bleed and/or boundary-layer energization that is sometimes used to shorten the length of supersonic diffusers and improve performance. In addition, the walls of the entire diffuser  30  may be roughened. 
         [0037]    In  FIG. 2 , the compression is shown as a centered P-M compression. It is not essential that the compression be centered; it is, however, convenient. The focal point of the compression is above the upper surface  42  and is the origin for the x,y coordinate system shown in  FIG. 2 , where the x axis is parallel to {right arrow over (w)} 1 . The focal point could, alternatively, be located below surface  44 ; in this circumstance, the diffuser would curve downward. 
         [0038]    The compression is bounded by its leading edge (LE) and its trailing edge (TE), which are straight characteristics of the same family. Along straight characteristic, all flow properties, such as the pressure and Mach number, are constant. Thus, on the LE, the Mach number and pressure are M 1  and p 1 , while on the TE they are M 2  and p 2 . At start-up, a normal shock wave must pass through the system. For this to occur, the trailing edge Mach number must exceed unity. For this reason, and to avoid unsteady transonic effects in a compressive flow, M 2  is preferably limited to a minimum value of about 1.4. 
         [0039]    The parameter μ is the Mach angle and is shown for the LE and the TE. Both walls are P-M streamlines; consequently, the flow between the LE and TE is a simple wave region in which the straight characteristics pass through the origin. The flows upstream and downstream of the LE and TE, respectively, are uniform flows. The downstream flow has an angle θ 2  relative to the upstream flow. 
         [0040]      FIG. 3  shows the r,η coordinates that are used for the wall equations. Note that η is measured from the positive y-axis. The angular spread of the compression fan is η 1 -η 2 . All lengths, except h, scale with l 1 . Similarly, pressures can be normalized with p 1  and temperatures with the inlet stagnation temperature, T o . 
         [0041]    A non-dimensional evaluation requires only γ, M 1 , M 2 , and y 1 /l 1 , where γ is the ratio of specific heats, M 1  is the Mach number at the LE and M 2  is the Mach number at the TE. Dimensional results also require values for l 1 , h, p 1 , T o , and R, where p 1  is the static pressure at LE, T o  is the inlet stagnation temperature, and R is the universal gas constant. 
         [0042]    The viscosity of the gas at the inlet, {tilde over (μ)} 1 , only appears in the Reynolds number. The two Mach numbers are the prescribed values across the compression. The y 1  parameter is shown in  FIG. 2 ; it has a negative value that must have an absolute value greater than l 1 . 
         [0043]    The slope (μ 1 ) of LE, measured from the x-axis, is given by 
         [0000]    
       
         
           
             
               μ 
               1 
             
             = 
             
               
                 sin 
                 
                   - 
                   1 
                 
               
                
               
                 ( 
                 
                   1 
                   
                     M 
                     1 
                   
                 
                 ) 
               
             
           
         
       
     
         [0000]    A diffuser inlet Reynolds number, based on the l 1  length, is given by 
         [0000]    
       
         
           
             Re 
             = 
             
               
                 
                   ( 
                   
                     γ 
                     
                       RT 
                       1 
                     
                   
                   ) 
                 
                 
                   1 
                   / 
                   2 
                 
               
                
               
                 
                   ( 
                   
                     
                       pM 
                        
                       
                           
                       
                        
                        
                     
                     
                       μ 
                       ~ 
                     
                   
                   ) 
                 
                 1 
               
             
           
         
       
     
         [0044]    The p o  parameter is the inlet stagnation pressure, while p r  is an estimate for the recovered pressure. The normal shock recovered pressure is p nr , while the l 2 /l 1  quantity is the area ratio across the compression. The Mach number M 3  occurs just downstream of the normal shock; note that 
         [0000]      M 1 &gt;M 2 &gt;1&gt;M 3    
         [0045]    The mass flow rate is {dot over (m)}. 
         [0046]    The derivation is outlined for the a and b walls that are streamlines of a centered P-M compression.  FIG. 3  shows the r,η velocity components, v r  and v n , that satisfy the simple-wave region relation 
         [0000]    
       
         
           
             
               
                 
                   
                     
                        
                       
                         v 
                         r 
                       
                     
                     
                        
                       η 
                     
                   
                   = 
                   
                     
                       v 
                       n 
                     
                     = 
                     
                       - 
                       a 
                     
                   
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where the speed of sound, a, is given by 
         [0000]    
       
         
           
             
               
                 
                   
                     a 
                     2 
                   
                   = 
                   
                     
                       γ 
                        
                       
                         p 
                         ρ 
                       
                     
                     = 
                     
                       
                         
                           
                             γ 
                             - 
                             1 
                           
                           2 
                         
                          
                         
                           q 
                           m 
                           2 
                         
                       
                       - 
                       
                         
                           
                             γ 
                             - 
                             1 
                           
                           2 
                         
                          
                         
                           ( 
                           
                             
                               v 
                               r 
                               2 
                             
                             + 
                             
                               v 
                               n 
                               2 
                             
                           
                           ) 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   2 
                   ) 
                 
               
             
           
         
       
     
         [0000]    The maximum speed is 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       q 
                       m 
                     
                     = 
                     
                       
                         
                           ( 
                           
                             
                               2 
                               
                                 γ 
                                 - 
                                 1 
                               
                             
                              
                             
                               X 
                               1 
                             
                           
                           ) 
                         
                         
                           1 
                           / 
                           2 
                         
                       
                        
                       
                         a 
                         1 
                       
                     
                   
                    
                   
                     
 
                   
                    
                   where 
                    
                   
                     
 
                   
                    
                   
                     
                       X 
                       1 
                     
                     = 
                     
                       1 
                       + 
                       
                         
                           
                             γ 
                             - 
                             1 
                           
                           2 
                         
                          
                         
                           M 
                           1 
                           2 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
         [0000]    Combining the above yields 
         [0000]    
       
         
           
             
               
                 
                   
                     
                        
                       
                         v 
                         r 
                       
                     
                     
                        
                       η 
                     
                   
                   = 
                   
                     
                       - 
                       
                         
                           ( 
                           
                             
                               γ 
                               - 
                               1 
                             
                             
                               γ 
                               + 
                               1 
                             
                           
                           ) 
                         
                         
                           1 
                           / 
                           2 
                         
                       
                     
                      
                     
                       
                         ( 
                         
                           
                             q 
                             m 
                             2 
                           
                           - 
                           
                             v 
                             r 
                             2 
                           
                         
                         ) 
                       
                       
                         1 
                         / 
                         2 
                       
                     
                   
                 
               
               
                 
                   ( 
                   4 
                   ) 
                 
               
             
           
         
       
     
         [0047]    Upon integration, we have 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       v 
                       r 
                     
                     
                       q 
                       m 
                     
                   
                   = 
                   
                     sin 
                      
                     
                         
                     
                      
                     z 
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
           
         
       
     
         [0048]    The z variable is introduced for analytical convenience. It is defined by a linear relation with respect to η 
         [0000]    
       
         
           
             
               
                 
                   z 
                   = 
                   
                     
                       
                         
                           ( 
                           
                             
                               γ 
                               - 
                               1 
                             
                             
                               γ 
                               + 
                               1 
                             
                           
                           ) 
                         
                         
                           1 
                           / 
                           2 
                         
                       
                        
                       
                         ( 
                         
                           
                             η 
                             1 
                           
                           - 
                           η 
                         
                         ) 
                       
                     
                     + 
                     
                       
                         
                           sin 
                           
                             - 
                             1 
                           
                         
                          
                         
                           ( 
                           
                             
                               - 
                               
                                 
                                   γ 
                                   - 
                                   1 
                                 
                                 2 
                               
                             
                              
                             
                               
                                 Z 
                                 1 
                               
                               
                                 X 
                                 1 
                               
                             
                           
                           ) 
                         
                       
                       
                         1 
                         / 
                         2 
                       
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
         [0000]    but is also related to the Mach number by 
         [0000]    
       
         
           
             
               
                 
                   
                     M 
                     2 
                   
                   = 
                   
                     1 
                     + 
                     
                       
                         
                           γ 
                           + 
                           1 
                         
                         
                           γ 
                           - 
                           1 
                         
                       
                        
                       
                         tan 
                         2 
                       
                        
                       z 
                     
                   
                 
               
               
                 
                   ( 
                   7 
                   ) 
                 
               
             
           
         
       
     
         [0049]    The differential equation for the wall shape is 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       1 
                       r 
                     
                      
                     
                       
                          
                         r 
                       
                       
                          
                         η 
                       
                     
                   
                   = 
                   
                     
                       cot 
                        
                       
                           
                       
                        
                       μ 
                     
                     = 
                     
                       
                         
                           ( 
                           
                             
                               M 
                               2 
                             
                             - 
                             1 
                           
                           ) 
                         
                         
                           1 
                           / 
                           2 
                         
                       
                       = 
                       
                         
                           - 
                           
                             
                               ( 
                               
                                 
                                   γ 
                                   + 
                                   1 
                                 
                                 
                                   γ 
                                   - 
                                   1 
                                 
                               
                               ) 
                             
                             
                               1 
                               / 
                               2 
                             
                           
                         
                          
                         tan 
                          
                         
                             
                         
                          
                         z 
                       
                     
                   
                 
               
               
                 
                   ( 
                   8 
                   ) 
                 
               
             
           
         
       
     
         [0000]    Equations (6) and (8) can be integrated, to yield 
         [0000]    
       
         
           
             
               
                 
                   
                     r 
                     a 
                   
                   = 
                   
                     
                       - 
                       
                         M 
                         1 
                       
                     
                      
                     
                       
                         
                           y 
                           1 
                         
                          
                         
                           ( 
                           
                             X 
                             
                               X 
                               1 
                             
                           
                           ) 
                         
                       
                       
                         
                           ( 
                           
                             γ 
                             + 
                             1 
                           
                           ) 
                         
                         / 
                         
                           [ 
                           
                             2 
                              
                             
                               ( 
                               
                                 γ 
                                 - 
                                 1 
                               
                               ) 
                             
                           
                           ] 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     9 
                      
                     a 
                   
                   ) 
                 
               
             
             
               
                 
                   
                     r 
                     b 
                   
                   = 
                   
                     
                       - 
                       
                         
                           M 
                           1 
                         
                          
                         
                           ( 
                           
                             
                               l 
                               1 
                             
                             + 
                             
                               y 
                               1 
                             
                           
                           ) 
                         
                       
                     
                      
                     
                       
                         ( 
                         
                           X 
                           
                             X 
                             1 
                           
                         
                         ) 
                       
                       
                         
                           ( 
                           
                             γ 
                             + 
                             1 
                           
                           ) 
                         
                         / 
                         
                           [ 
                           
                             2 
                              
                             
                               ( 
                               
                                 γ 
                                 - 
                                 1 
                               
                               ) 
                             
                           
                           ] 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     9 
                      
                     b 
                   
                   ) 
                 
               
             
           
         
       
     
         [0050]    (Note that y 1  and l 1 +y 1  are negative.) 
         [0051]    Wall values are evaluated at an arbitrary sequence of Mach numbers, given by 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       M 
                       i 
                     
                     = 
                     
                       
                         M 
                         1 
                       
                       - 
                       
                         
                           ( 
                           
                             
                               
                                 M 
                                 1 
                               
                               - 
                               
                                 M 
                                 2 
                               
                             
                             N 
                           
                           ) 
                         
                          
                         
                           ( 
                           
                             i 
                             - 
                             1 
                           
                           ) 
                         
                       
                     
                   
                   , 
                   
                     i 
                     = 
                     1 
                   
                   , 
                   2 
                   , 
                   … 
                    
                   
                       
                   
                   , 
                   
                     N 
                     + 
                     1 
                   
                 
               
               
                 
                   ( 
                   10 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where i=1 corresponds to points a 1  and b 1 , in  FIG. 2 , and i=N+1 corresponds to points a 2  and b 2 . The η i  values are 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       η 
                       i 
                     
                     = 
                     
                       
                         η 
                         1 
                       
                       - 
                       
                         
                           
                             ( 
                             
                               
                                 γ 
                                 + 
                                 1 
                               
                               
                                 γ 
                                 - 
                                 1 
                               
                             
                             ) 
                           
                           
                             1 
                             / 
                             2 
                           
                         
                          
                         
                           [ 
                           
                             
                               z 
                               i 
                             
                             + 
                             
                               
                                 
                                   sin 
                                   
                                     - 
                                     1 
                                   
                                 
                                  
                                 
                                   ( 
                                   
                                     
                                       
                                         γ 
                                         - 
                                         1 
                                       
                                       2 
                                     
                                      
                                     
                                       
                                         Z 
                                         1 
                                       
                                       
                                         X 
                                         i 
                                       
                                     
                                   
                                   ) 
                                 
                               
                               
                                 1 
                                 / 
                                 2 
                               
                             
                           
                           ] 
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   where 
                    
                   
                     
 
                   
                    
                   
                     X 
                     = 
                     
                       
                         1 
                         + 
                         
                           
                             
                               γ 
                               - 
                               1 
                             
                             2 
                           
                            
                           
                             M 
                             2 
                           
                            
                           
                               
                           
                            
                           and 
                            
                           
                               
                           
                            
                           Z 
                         
                       
                       = 
                       
                         
                           M 
                           2 
                         
                         - 
                         1 
                       
                     
                   
                 
               
               
                 
                   ( 
                   11 
                   ) 
                 
               
             
           
         
       
     
         [0052]    The two wall shapes are then given by the relations 
         [0000]      x ai =r ai  sin η i , y ai =r ai  cos η 1   (12a) 
         [0000]      x bi =r bi  sin η i , y bi =r bi  cos η i   (12b) 
         [0000]    where r ai  and r bi  are provided by Eqs. (9a) and (9b). The walls have similar shapes, since they are Prandtl-Meyer streamlines. The diagonal distance, L, is the length between points a 1  and b 2 . It is used as a measure of the longitudinal size of the diffuser, and is given by 
         [0000]        L =[( x   b2   −x   a1 ) 2 +( y   b2   −y   a1 ) 2 ] 1/2   (13) 
         [0000]    Note that this estimate does not include a length associated with the short throat section. 
         [0053]      FIG. 4  is a closer view of the transonic region of the preferred embodiment of supersonic diffuser according to the present invention. Downstream of the TE characteristic there is a short duct  36  which is the throat, or transonic, region. The cross-sectional area of the duct  36  should gradually increase with a wall slope of about 0.5° to 1°, relative to the direction of gas flow. This compensates for boundary-layer growth and helps stabilize the shock wave. Inside the duct  36  is a porous fibrous material  38  reminiscent of the filter material used in a heating or air-conditioning circuit, but here more open or porous. The fibrous material  38  provides a small, distributed disturbance that acts as a shock holder for a weak normal shock wave. Alternatively, several screens (not shown) can be used as a shock holder. Downstream of the fibrous material  38  is a high Mach number, subsonic flow. Additional pressure recovery, if needed, can be obtained with a subsonic diffuser (not shown), whose length can be shortened, e.g., using slat or vane inserts. 
         [0054]    The diffuser described above is subject to the following constraints for supersonic diffusers. A near vacuum pressure can be used for starting, the boundary layer must be turbulent, and the configuration is designed for a prescribed inlet Mach number. However, in contrast to a conventional supersonic diffuser, the diffuser shown in  FIG. 2  is curved. 
         [0055]      FIG. 5  shows an alternative design for a supersonic diffuser that is based on the lens analogy (L-A). The lens analogy is based on conceptual similarities between nozzles and diffusers in isentropic gas flow and what several lenses can do with an optical beam. 
         [0056]    L-A flow can be complicated by the presence of a limit line. A limit line does not occur in a P-M flow. A limit line was discovered with the hodograph transformation. It is present when the Jacobian of the transformation is zero along a curve that occurs in both the hodograph and physical planes. Because of the Jacobian association, discussion of limit lines has generally been restricted to the two-dimensional hodograph transformation. The fluid acceleration is infinite at a limit line. This interpretation allows the concept to be applied to flows not associated with a hodograph or any other type of transformation. It is the infinite acceleration interpretation that is utilized here. 
         [0057]    In the physical plane, an inviscid flow cannot be continued across a limit line, i.e., a streamline that encounters a limit line will reflect off of it by making a 1800 turn. Thus, flows with a limit line are analytically multi-valued. Limit lines are only present in sonic or supersonic flows. 
         [0058]    A limit line is an analytical phenomenon, somewhat analogous to the crossing of Mach lines of the same family. Neither is physically possible; in both cases a major change in the flow field must occur. Overlapping Mach lines are replaced with a shock wave. Since a limit line only occurs at low supersonic Mach numbers, it can be avoided by terminating the compressive flow, where a limit line can occur, at a slightly higher Mach number, thereby leaving the bulk of the flow field unaltered. In contrast to a P-M diffuser, in a L-A diffuser the M 2  value is not arbitrary. 
         [0059]      FIG. 5  shows the lower half of a symmetric L-A diffuser where M 1 =4, M 2 =2, and the ratio of specific heats, γ, is 1.4. The M 2  value is sufficiently large such that a limit line is not present.  FIG. 5  does not show the upstream laser or the downstream throat section. The throat section is similar to the throat section of the diffuser discussed above, and is shown in  FIG. 3 . The lines A′A and B′B are left-running characteristics, while AB′ and BC′ are right-running characteristics. These characteristics separate regions of different flow types. Upstream of AA′ and downstream of BC′ are uniform flow regions with Mach numbers M 1  and M 2 , where 
         [0000]      M 1 &gt;M 2 &gt;1&gt;M 3    
         [0000]    as before. 
         [0060]    Regions A′AB′A′ and B′BC′B′ are simple wave regions. Region AB′BA is a non-simple wave region, but is easily modeled, since it is a line-source flow region whose (suction) origin is the origin of the x,y coordinate system. 
         [0061]    As is evident in  FIG. 5 , the wall  52  has an inflection at point B′. A straight-wall section that points to the origin of the coordinate system can be inserted at point B′. When this is done, the pressure gradient along the wall is discontinuous at the two points where the straight and curved wall sections intersect. This discontinuity does not occur when the curved wall sections meet at the inflection point. Moreover, the diffuser has a minimum length without straight wall segments. Therefore, a L-A diffuser without straight wall sections is preferred. 
         [0062]    Much of the earlier discussion relating to  FIGS. 2 and 3  also applies here. 
         [0063]    The shape of the wall  52  is parametrically given in terms of the Mach number. For the wall between A′ and B′, 
         [0000]    
       
         
           
             
               
                 
                   x 
                   = 
                   
                     
                       - 
                       
                         
                           l 
                           1 
                         
                         2 
                       
                     
                      
                     
                       1 
                       φ 
                     
                      
                     
                       
                         M 
                         1 
                       
                       M 
                     
                      
                     
                       
                         
                           ( 
                           
                             X 
                             
                               X 
                               1 
                             
                           
                           ) 
                         
                         
                           
                             ( 
                             
                               γ 
                               + 
                               1 
                             
                             ) 
                           
                           / 
                           
                             [ 
                             
                               2 
                                
                               
                                 ( 
                                 
                                   γ 
                                   - 
                                   1 
                                 
                                 ) 
                               
                             
                             ] 
                           
                         
                       
                        
                       
                         [ 
                         
                           
                             cos 
                              
                             
                                 
                             
                              
                             θ 
                           
                           + 
                           
                             
                               ( 
                               
                                 φ 
                                 - 
                                 θ 
                               
                               ) 
                             
                              
                             
                               ( 
                               
                                 
                                   
                                     Z 
                                     
                                       1 
                                       / 
                                       2 
                                     
                                   
                                    
                                   cos 
                                    
                                   
                                       
                                   
                                    
                                   θ 
                                 
                                 - 
                                 
                                   sin 
                                    
                                   
                                       
                                   
                                    
                                   θ 
                                 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     14 
                      
                     a 
                   
                   ) 
                 
               
             
             
               
                 
                   
                     y 
                     = 
                     
                       
                         - 
                         
                           
                             l 
                             1 
                           
                           2 
                         
                       
                        
                       
                         1 
                         φ 
                       
                        
                       
                         
                           M 
                           1 
                         
                         M 
                       
                        
                       
                         
                           
                             ( 
                             
                               X 
                               
                                 X 
                                 1 
                               
                             
                             ) 
                           
                           
                             
                               ( 
                               
                                 γ 
                                 + 
                                 1 
                               
                               ) 
                             
                             / 
                             
                               [ 
                               
                                 2 
                                  
                                 
                                   ( 
                                   
                                     γ 
                                     - 
                                     1 
                                   
                                   ) 
                                 
                               
                               ] 
                             
                           
                         
                          
                         
                           [ 
                           
                             
                               sin 
                                
                               
                                   
                               
                                
                               θ 
                             
                             + 
                             
                               
                                 ( 
                                 
                                   φ 
                                   - 
                                   θ 
                                 
                                 ) 
                               
                                
                               
                                 ( 
                                 
                                   
                                     
                                       Z 
                                       
                                         1 
                                         / 
                                         2 
                                       
                                     
                                      
                                     sin 
                                      
                                     
                                         
                                     
                                      
                                     θ 
                                   
                                   - 
                                   
                                     cos 
                                      
                                     
                                         
                                     
                                      
                                     θ 
                                   
                                 
                                 ) 
                               
                             
                           
                           ] 
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   where 
                 
               
               
                 
                   ( 
                   
                     14 
                      
                     b 
                   
                   ) 
                 
               
             
             
               
                 
                   
                     M 
                     2 
                   
                   &lt; 
                   
                     M 
                     
                       B 
                       ′ 
                     
                   
                   ≤ 
                   M 
                   ≤ 
                   
                     M 
                     1 
                   
                 
               
               
                 
                   ( 
                   
                     15 
                      
                     a 
                   
                   ) 
                 
               
             
             
               
                 
                   θ 
                   = 
                   
                     
                       v 
                       1 
                     
                     - 
                     
                       v 
                        
                       
                         ( 
                         M 
                         ) 
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     15 
                      
                     b 
                   
                   ) 
                 
               
             
             
               
                 
                   
                     φ 
                     ≡ 
                     
                       θ 
                       
                         B 
                         ′ 
                       
                     
                   
                   = 
                   
                     
                       1 
                       2 
                     
                      
                     
                       ( 
                       
                         
                           v 
                           1 
                         
                         - 
                         
                           v 
                           2 
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   
                     15 
                      
                     c 
                   
                   ) 
                 
               
             
           
         
       
     
         [0000]    and v is the Prandtl-Meyer function. M B′  is the implicit solution of 
         [0000]    
       
         
           
             
               
                 
                   
                     v 
                      
                     
                       ( 
                       
                         M 
                         
                           B 
                           ′ 
                         
                       
                       ) 
                     
                   
                   = 
                   
                     
                       1 
                       2 
                     
                      
                     
                       ( 
                       
                         
                           v 
                           1 
                         
                         + 
                         
                           v 
                           2 
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   
                     15 
                      
                     d 
                   
                   ) 
                 
               
             
           
         
       
     
         [0064]    For instance, when M=M 1 , we obtain 
         [0000]      θ 1 =0  (16a) 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       x 
                       1 
                     
                     
                       l 
                       1 
                     
                   
                   = 
                   
                     
                       - 
                       
                         1 
                         2 
                       
                     
                      
                     
                       ( 
                       
                         
                           Z 
                           1 
                           
                             1 
                             / 
                             2 
                           
                         
                         + 
                         
                           1 
                           φ 
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   
                     16 
                      
                     b 
                   
                   ) 
                 
               
             
             
               
                 
                   
                     
                       y 
                       1 
                     
                     
                       l 
                       1 
                     
                   
                   = 
                   
                     - 
                     
                       1 
                       2 
                     
                   
                 
               
               
                 
                   ( 
                   
                     16 
                      
                     c 
                   
                   ) 
                 
               
             
           
         
       
     
         [0000]    where the y 1  value is the vertical coordinate of point A′. 
         [0065]    For the wall between points B′ and C′, we have 
         [0000]    
       
         
           
             
               
                 
                   x 
                   = 
                   
                     
                       
                         l 
                         1 
                       
                       2 
                     
                      
                     
                       1 
                       φ 
                     
                      
                     
                       
                         M 
                         1 
                       
                       M 
                     
                      
                     
                       
                         
                           ( 
                           
                             X 
                             
                               X 
                               1 
                             
                           
                           ) 
                         
                         
                           
                             ( 
                             
                               γ 
                               + 
                               1 
                             
                             ) 
                           
                           / 
                           
                             [ 
                             
                               2 
                                
                               
                                 ( 
                                 
                                   γ 
                                   - 
                                   1 
                                 
                                 ) 
                               
                             
                             ] 
                           
                         
                       
                        
                       
                         [ 
                         
                           
                             
                               - 
                               cos 
                             
                              
                             
                                 
                             
                              
                             θ 
                           
                           + 
                           
                             
                               ( 
                               
                                 φ 
                                 - 
                                 θ 
                               
                               ) 
                             
                              
                             
                               ( 
                               
                                 
                                   sin 
                                    
                                   
                                       
                                   
                                    
                                   θ 
                                 
                                 + 
                                 
                                   
                                     Z 
                                     
                                       1 
                                       / 
                                       2 
                                     
                                   
                                    
                                   cos 
                                    
                                   
                                       
                                   
                                    
                                   θ 
                                 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     17 
                      
                     a 
                   
                   ) 
                 
               
             
             
               
                 
                   
                     y 
                     = 
                     
                       
                         - 
                         
                           
                             l 
                             1 
                           
                           2 
                         
                       
                        
                       
                         1 
                         φ 
                       
                        
                       
                         
                           M 
                           1 
                         
                         M 
                       
                        
                       
                         
                           
                             ( 
                             
                               X 
                               
                                 X 
                                 1 
                               
                             
                             ) 
                           
                           
                             
                               ( 
                               
                                 γ 
                                 + 
                                 1 
                               
                               ) 
                             
                             / 
                             
                               [ 
                               
                                 2 
                                  
                                 
                                   ( 
                                   
                                     γ 
                                     - 
                                     1 
                                   
                                   ) 
                                 
                               
                               ] 
                             
                           
                         
                          
                         
                           [ 
                           
                             
                               sin 
                                
                               
                                   
                               
                                
                               θ 
                             
                             + 
                             
                               
                                 ( 
                                 
                                   φ 
                                   - 
                                   θ 
                                 
                                 ) 
                               
                                
                               
                                 ( 
                                 
                                   
                                     cos 
                                      
                                     
                                         
                                     
                                      
                                     θ 
                                   
                                   + 
                                   
                                     
                                       Z 
                                       
                                         1 
                                         / 
                                         2 
                                       
                                     
                                      
                                     sin 
                                      
                                     
                                         
                                     
                                      
                                     θ 
                                   
                                 
                                 ) 
                               
                             
                           
                           ] 
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   where 
                 
               
               
                 
                   ( 
                   
                     17 
                      
                     b 
                   
                   ) 
                 
               
             
             
               
                 
                   1 
                   &lt; 
                   
                     M 
                     l 
                   
                   &lt; 
                   
                     M 
                     2 
                   
                   ≤ 
                   M 
                   ≤ 
                   
                     M 
                     
                       B 
                       ′ 
                     
                   
                 
               
               
                 
                   ( 
                   
                     18 
                      
                     a 
                   
                   ) 
                 
               
             
             
               
                 
                   θ 
                   = 
                   
                     
                       v 
                        
                       
                         ( 
                         M 
                         ) 
                       
                     
                     - 
                     
                       v 
                       2 
                     
                   
                 
               
               
                 
                   ( 
                   
                     18 
                      
                     b 
                   
                   ) 
                 
               
             
           
         
       
     
         [0000]    and M l  is the limit line Mach number, given by 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         8 
                         
                           γ 
                           + 
                           1 
                         
                       
                        
                       
                         
                           Z 
                           l 
                           
                             3 
                             / 
                             2 
                           
                         
                         
                           M 
                           l 
                           4 
                         
                       
                     
                     + 
                     
                       v 
                        
                       
                         ( 
                         
                           M 
                           l 
                         
                         ) 
                       
                     
                   
                   = 
                   
                     v 
                     1 
                   
                 
               
               
                 
                   ( 
                   
                     18 
                      
                     c 
                   
                   ) 
                 
               
             
           
         
       
     
         [0066]    For instance, when M=M 2 , we obtain 
         [0000]    
       
         
           
             
               
                 
                   
                     θ 
                     2 
                   
                   = 
                   0 
                 
               
               
                 
                   ( 
                   
                     19 
                      
                     a 
                   
                   ) 
                 
               
             
             
               
                 
                   
                     
                       x 
                       2 
                     
                     
                       l 
                       2 
                     
                   
                   = 
                   
                     
                       1 
                       2 
                     
                      
                     
                       ( 
                       
                         
                           Z 
                           2 
                           
                             1 
                             / 
                             2 
                           
                         
                         - 
                         
                           1 
                           φ 
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   
                     19 
                      
                     b 
                   
                   ) 
                 
               
             
             
               
                 
                   
                     
                       y 
                       2 
                     
                     
                       l 
                       2 
                     
                   
                   = 
                   
                     - 
                     
                       1 
                       2 
                     
                   
                 
               
               
                 
                   ( 
                   
                     19 
                      
                     c 
                   
                   ) 
                 
               
             
           
         
       
     
         [0000]    where the y 2  value is the vertical coordinate of point C′. 
         [0067]    Both Eqs. (14) and (17) hold when M=M B′ , as is evident since θ B′  equals φ. The length of the diffuser is given by 
         [0000]    
       
         
           
             
               
                 
                   L 
                   = 
                   
                     
                       1 
                       2 
                     
                      
                     
                       
                         l 
                         1 
                       
                        
                       
                         [ 
                         
                           
                             
                               
                                 M 
                                 1 
                               
                               
                                 M 
                                 2 
                               
                             
                              
                             
                               
                                 ( 
                                 
                                   
                                     X 
                                     2 
                                   
                                   
                                     X 
                                     1 
                                   
                                 
                                 ) 
                               
                               
                                 
                                   ( 
                                   
                                     γ 
                                     + 
                                     1 
                                   
                                   ) 
                                 
                                 / 
                                 
                                   [ 
                                   
                                     2 
                                      
                                     
                                       ( 
                                       
                                         γ 
                                         - 
                                         1 
                                       
                                       ) 
                                     
                                   
                                   ] 
                                 
                               
                             
                              
                             
                               ( 
                               
                                 
                                   Z 
                                   2 
                                   
                                     1 
                                     / 
                                     2 
                                   
                                 
                                 - 
                                 
                                   1 
                                   φ 
                                 
                               
                               ) 
                             
                           
                           + 
                           
                             Z 
                             1 
                             
                               1 
                               / 
                               2 
                             
                           
                           + 
                           
                             1 
                             φ 
                           
                         
                         ] 
                       
                     
                   
                 
               
               
                 
                   ( 
                   20 
                   ) 
                 
               
             
           
         
       
     
         [0000]    Three cases are evaluated for each approach, see Table 1. The specific heat ratio is 1.4, and the M 1  values span a range expected for COIL. 
         [0000]    
       
         
               
             
               
               
               
               
               
               
               
             
           
               
                 TABLE 1 
               
             
             
               
                   
               
               
                 Performance Summary 
               
             
          
           
               
                 Type 
                 M 1   
                 M 2   
                 η r   
                 η nr   
                 L/l 1   
                 
                   
                     
                       
                         
                           d 
                            
                           
                             ( 
                             
                               p 
                               / 
                               
                                 p 
                                 o 
                               
                             
                             ) 
                           
                         
                         
                           
                             d 
                              
                             
                               ( 
                               
                                 s 
                                 / 
                                 
                                   l 
                                   1 
                                 
                               
                               ) 
                             
                           
                           m 
                         
                       
                     
                   
                 
               
               
                   
               
               
                 P-M 
                 2.5 
                 1.4 
                 0.9582 
                 0.4990 
                 3.576 
                 0.4334 
               
               
                   
                 3.5 
                 1.4 
                 0.9582 
                 0.2129 
                 5.228 
                 1.116 
               
               
                   
                 4.5 
                 1.4 
                 0.9582 
                 0.0917 
                 6.755 
                 28.62 
               
               
                 L-A 
                 2.5 
                 1.3 
                 0.9794 
                 0.4990 
                 2.349 
                 6.278 
               
               
                   
                 3.5 
                 1.5 
                 0.9298 
                 0.2129 
                 2.790 
                 13.94 
               
               
                   
                 4.5 
                 1.7 
                 0.8557 
                 0.0917 
                 3.224 
                 27.89 
               
               
                   
               
             
          
         
       
     
         [0068]    The only irreversible process stems from the normal shock whose upstream Mach number is M 2 . Hence η nr  has a constant value for a P-M diffuser. For the L-A diffuser, the limit line Mach number values are 1.26, 1.43, or 1.61, respectively, when M 1  is 2.5, 3.5, or 4.5. The increasing M 2  value causes some fall-off in η r . Nevertheless, both approaches provide excellent pressure recovery performance as compared to a conventional diffuser, whose nominal performance is represented by its normal shock recovered pressure efficiency η nr . When the inlet Mach number is 4.5, the difference in recovered pressure, as provided by η r  and η nr , is about an order-of-magnitude. 
         [0069]    Table 1 shows that L/l 1  favors the L-A approach. At lower inlet Mach numbers, the P-M approach is favored for its smaller value for the maximum value of the adverse pressure gradient. This changes when M 1 =4.5. Note that the gradient is rapidly increasing with M 1 . The rate of increase is less for the L-A diffuser, because M 2  is also increasing. Unfortunately, the large gradient value occurs where the Mach number is smallest, and, therefore, where the supersonic boundary layer is most prone to separate. 
         [0070]    In both approaches, the wall pressure gradient increases from a near-zero value at the inlet to its maximum value at the outlet. Most of the increase occurs near the downstream end of the compression region. This is evident in  FIG. 5 , where the relatively short B′C′ wall section provides much of the Mach number decrease. It is essential that the maximum value of the wall pressure gradient not cause boundary-layer separation. This value increases rapidly with M 1 , but M 1  is fixed by the upstream laser system. The maximum value of the pressure gradient decreases gradually with an increasing M 2  value for the P-M diffuser, and much more rapidly with an L-A diffuser. Recall that increasing M 2  decreases the diffuser&#39;s efficiency. Alternatively, a non-centered P-M compression can be used. 
         [0071]    A particularly effective approach for decreasing the maximum value of the pressure gradient is to use “stacking.” In this approach, two or more diffusers are in series (outlet of first diffuser coupled to the inlet of the second), as shown schematically in  FIG. 6 . Region  2  has a Mach number now intermediate between M 1  and M 3 . Of course, the length of a stacked diffuser exceeds that of a single diffuser operating between M 1  and M 3 . 
         [0072]    The foregoing discussion insures that both approaches can operate without boundary-layer separation. 
         [0073]    A new approach is introduced for the supersonic portion and the downstream throat section of a diffuser. Supersonic compression is considered using a Prandtl-Meyer compression (P-M) or the lens analogy (L-A). In each case, the flow is shock free with an isentropic compression; the engineering logic for this expectation is discussed. A short throat section is used, in both cases, that contains a shock holder for a weak normal shock wave. The inviscid wall design of the supersonic sections is provided. A subsonic diffuser may, or may not, be located downstream of the throat section. The two approaches for supersonic compression are analytically evaluated using three criteria: the efficiency of the pressure recovery process, the length of the supersonic portion of the diffuser, and the maximum value of the adverse pressure gradient along the wall of the supersonic section. Avoidance of a limit line and shock free inlet flow can be viewed as additional criteria. Both approaches have comparable efficiencies that significantly exceed current practice; the recovered pressure is close to its isentropic maximum value. A small maximum value for the pressure gradient is important for inhibiting turbulent boundary-layer separation. In turn, this is essential for the physical realization of the approach. 
         [0074]    Having thus described the present invention by reference to certain of its preferred embodiments, it is noted that the embodiments disclosed are illustrative rather than limiting in nature and that a wide range of variations, modifications, changes and substitutions are contemplated in the foregoing disclosure and, in some instances, some features of the present invention may be employed without a corresponding use of the other features. Many such variations and modifications may be considered obvious and desirable by those skilled in the art based upon a review of the foregoing description of preferred embodiments. Accordingly, it is appropriate that the appended claims be construed broadly and in a manner consistent with the scope of the invention.