Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bluff body”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

Application of an Upwind High Resolution Finite-Differencing Scheme and Multigrid Method in Steady-State Incompressible Flow Simulations

The analysis and design of a submarine propulsor requires the ability to predict the characteristics of both laminar and turbulent flows to a higher degree of accuracy. This report presents results of certain benchmark computations based on an upwind, high-resolution, finite-differencing Navier-Stokes solver. The purpose of the computations is to evaluate the ability, the accuracy and the performance of the solver in the simulation of detailed features of viscous flows. Features of interest include flow separation and reattachment, surface pressure and skin friction distributions. Those features are particularly relevant to the propulsor analysis. Test cases with a wide range of Reynolds numbers are selected; therefore, the effects of the convective and the diffusive terms of the solver can be evaluated separately. Test cases include flows over bluff bodies, such as circular cylinders and spheres, at various low Reynolds numbers, flows over a flat plate with and without turbulence effects, and turbulent flows over axisymmetric bodies with and without propulsor effects. Finally, to enhance the iterative solution procedure, a full approximation scheme V-cycle multigrid method is implemented. Preliminary results indicate that the method significantly reduces the computational effort.

NAVIER STOKES EQUATIONS↗

Study of the Unsteady Flow Features on a Stalled Wing

The occurrence of large scale structures in the post stall flow over a rectangular wing at high angles of attack was investigated in a small-scale subsonic wind tunnel. Mean and time dependent measurements within the separated flow field suggest the existence of two distinct angle of attack regimes beyond wing stall. The shallow stall regime occurs over a narrow range of incidence angles (2-3 deg.) immediately following the inception of leading edge separation. In this regime, the principal mean flow structures, termed stall cells, are manifested as a distinct spanwise periodicity in the chordwise extent of the separated region on the model surface with possible lateral mobility not previously reported. Within the stall cells and on the wing surface, large amplitude pressure fluctuations occur with a frequency much lower than anticipated for bluff body shedding, and with minimum effect in the far wake. In the deep stall regime, stall cells are not observed and the separated region near the model is relatively free of large amplitude pressure disturbances.

Yon, Steven A.↗

Large-Vortex Capture by a Wing at Very High Angles of Attack

In generating the lift on a wing, the static stall is a severe barrier. As the angle of attack, alpha, increases to the stall angle, alpha(sub stall) the flow separation point on the upper surface of the wing moves to the leading edge, so that on a two-dimensional airfoil or a large-aspect-ratio wing, the lift abruptly drops to a very low level. Therefore, the first generation of aeronautical flow type, i.e., the attached steady flow, has been limited to alpha less than alpha(sub stall). Owing to the obvious importance in applications, therefore, a great effort has been made in the past two decades to enlarge the range of usable angles of attack by various flow controls for a large-aspect-ratio wing. Basically, relevant works fall into two categories. The first category is usually refereed to as separation control, which concentrates on partially separated flow at alpha less than alpha(sub stall). Since the first experimental study of Collins and Zelenevitz, there has been ample literature showing that a partially separated flow can be turned to almost fully attached by flow controls, so that the lift is recovered and the stall is delayed (for a recent work see Seifert et al.). It has been well established that, in this category, unsteady controls are much more effective than steady ones and can be realized at a very low power-input level (Wu et al.; Seifert et al.). The second and more ambitious category of relevant efforts is the post-stall lift enhancement. Its possibility roots at the existence of a second lift peak at a very high angle of attack. In fact, As alpha further increases from alpha(sub stall), the completely separated flow develops and gradually becomes a bluff-body flow. This flow gives a normal force to the airfoil with a lift component, which reaches a peak at a maximum utilizable angle of attack, alpha(sub m) approx.= 40 deg. This second peak is of the same level as the first lift peak at alpha(sub stall). Meanwhile, the drag is also quickly increased (e.g., Fage and Johansen ; Critzos et al.). Figure 1 shows a typical experimental lift and drag coefficients of NACA-0012 airfoil in this whole range of angle of attack. Obviously, without overcoming the lift crisis at alpha(sub stall) the second lift peak is completely useless. Thus, the ultimate goal of post-stall lift enhancement is to fill the lift valley after stall by flow controls, so that a wing and/or flap can work at the whole range of 0 deg less than alpha less than alpha(sub m). Relevant early experimental studies have been extensively reviewed by Wu et al., who concluded that, first, similar to the leading-edge vortex on a slender wing, the lift enhancement on a large-aspect-ratio wing should be the result of capturing a vortex on the upper surface of the wing; and, second, using steady controls cannot reach the goal, and one must rely on unsteady controls with low-level power input as well. Wu et al. also conjectured that the underlying physics of post-stall lift enhancement by unsteady controls consists of a chain of mechanisms: vortex layer instability - receptivity resonance - nonlinear streaming.

Wu, J. M.↗

Numerical Investigation of Flow Around Rectangular Cylinders with and Without Jets

The problem of flow past bluff bodies was studied extensively in the past. The problem of drag reduction is very important in many high speed flow applications. Considerable work has been done in this subject area in case of circular cylinders. The present study attempts to investigate the feasibility of drag reduction on a rectangular cylinder by flow injection by flow injection from the rear stagnation region. The physical problem is modeled as two-dimensional body and numerical analysis is carried out with and without trailing jets. A commercial code is used for this purpose. Unsteady computation is performed in case of rectangular cylinders with no trailing jets where as steady state computation is performed when jet is introduced. It is found that drag can be reduced by introducing jets with small intensity in rear stagnation region of the rectangular cylinders.

Tiwari, S. N .↗

Interaction of a Boundary Layer with a Turbulent Wake

The objective of this grant was to study the transition mechanisms on a flat-plate boundary layer interacting with the wake of a bluff body. This is a simplified configuration presented and designed to exemplify the phenomena that occur in multi-element airfoils, in which the wake of an upstream element impinges on a downstream one. Some experimental data is available for this configuration at various Reynolds numbers. The first task carried out was the implementation and validation of the immersed-boundary method. This was achieved by performing calculations of the flow over a cylinder at low and moderate Reynolds numbers. The low-Reynolds number results are discussed, which is enclosed as Appendix A. The high-Reynolds number results are presented in a paper in preparation for the Journal of Fluid Mechanics. We performed calculations of the wake-boundary-layer interaction at two Reynolds numbers, Re approximately equal to 385 and 1155. The first case is discussed and a comparison of the two calculations is reported. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. This is shown: long streaky structures appear in the boundary layer in correspondence of the three-dimensionalities in the rollers. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established. A third simulation was subsequently carried out at a higher Reynolds number, Re=3900. This calculation gave results similar to those of the Re=l155 case. Turbulence was established at fairly low Reynolds number, as a consequence of the high level of the free-stream perturbation. An instantaneous flow visualization for that case is shown. A detailed examination of flow statistics in the transitional and turbulent regions, including the evolution of the turbulent kinetic energy (TKE) budget and frequency spectra showed the formation and evolution of turbulent spots characteristic of the bypass transition mechanism. It was also observed that the turbulent eddies achieved an equilibrium, fully developed turbulent states first, as evidenced by the early agreement achieved by the terms in the TKE budget with those observed in turbulent flows. Once a turbulent Reynolds stress profile had been established, the velocity profile began to resemble a turbulent one, first in the inner region and later in the outer region of the wall layer. An extensive comparison of the three cases, including budgets, mean velocity and Reynolds stress profiles and flow visualization, is included. The results obtained are also presented.

Piomelli, Ugo↗

Field Effects of Buoyancy on a Premixed Turbulent Flame Studied by Particle Image Velocimetry

Typical laboratory flames for the scientific investigation of flame/turbulence interactions are prone to buoyancy effects. Buoyancy acts on these open flame systems and provides upstream feedbacks that control the global flame properties as well as local turbulence/flame interactions. Consequently the flame structures, stabilization limits, and turbulent reaction rates are directly or indirectly coupled with buoyancy. The objective of this study is to characterize the differences between premixed turbulent flames pointing upwards (1g), pointing downwards (-1g), and in microgravity (mg). The configuration is an inverted conical flame stabilized by a small cone-shaped bluff body that we call CLEAN Flames (Cone-Stabilized Lean Flames). We use two laser diagnostics to capture the velocity and scalar fields. Particle image velocimetry (PIV) measures the mean and root mean square velocities and planar imaging by the flame fronts method outlines the flame wrinkle topology. The results were obtained under typical conditions of small domestic heating systems such as water heaters, ovens, and furnaces. Significant differences between the 1g and -1g flames point to the need for including buoyancy contributions in theoretical and numerical calculations. In Earth gravity, there is a complex coupling of buoyancy with the turbulent flow and heat release in the flame. An investigation of buoyancy-free flames in microgravity will provide the key to discern gravity contributions. Data obtained in microgravity flames will provide the benchmark for interpreting and analyzing 1g and -1g flame results.

Cheng, Robert K.↗

High-Speed Schlieren Movies of Decelerators at Supersonic Speeds

Tests were conducted on several types of porous parachutes, a paraglider, and a simulated retrorocket. Mach numbers ranged from 1.8-3.0, porosity from 20-80 percent, and camera speeds from 1680-3000 feet per second (fps) in trials with porous parachutes. Trials of reefed parachutes were conducted at Mach number 2.0 and reefing of 12-33 percent at camera speeds of 600 fps. A flexible parachute with an inflatable ring in the periphery of the canopy was tested at Reynolds number 750,000 per foot, Mach number 2.85, porosity of 28 percent, and camera speed of 36oo fps. A vortex-ring parachute was tested at Mach number 2.2 and camera speed of 3000 fps. The paraglider, with a sweepback of 45 degrees at an angle of attack of 45 degrees was tested at Mach number 2.65, drag coefficient of 0.200, and lift coefficient of 0.278 at a camera speed of 600 fps. A cold air jet exhausting upstream from the center of a bluff body was used to simulate a retrorocket. The free-stream Mach number was 2.0, free-stream dynamic pressure was 620 lb/sq ft, jet-exit static pressure ratio was 10.9, and camera speed was 600 fps.

SUPERSONIC SPEED↗

Development of Hairpin Vortices in Turbulent Spots and End-Wall Transition

The end-stage phase of boundary layer transition is characterized by the development of hairpin-like vortices which evolve rapidly into patches of turbulent behavior. In general, the characteristics of the evolution form this hairpin stage to the turbulent stage is poorly understood, which has prompted the present experimental examination of hairpin vortex development and growth processes. Two topics of particular relevance to the workshop focus will be covered: 1) the growth of turbulent spots through the generatio and amalgamation of hairpin-like vortices, and 2) the development of hairpin vortices during transition in an end-wall junction flow. Brief summaries of these studies are described below. Using controlled generation of hairpin vortices by surface injection in a critical laminar boundary layer, detailed flow visualization studies have been done of the phases of growth of single hairpin vortices, from the initial hairgin generation, through the systematic generation of secondary hairpin-like flow structures, culminating in the evolution to a turbulent spot. The key to the growth process is strong vortex-surface interactions, which give rise to strong eruptive events adjacent to the surface, which results in the generation of subsequent hairpin vortex structures due to inviscid-viscuous interactions between the eruptive events and the free steam fluid. The general process of vortex-surface fluid interaction, coupled with subsequent interactions and amalgamation of the generated multiple hairpin-type vortices, is demonstrated as a physical mechanism for the growth and development of turbulent spots. When a boundary layer flow along a surface encounters a bluff body obstruction extending from the surface (such as cylinder or wing), the strong adverse pressure gradients generated by these types of flows result in the concentration of the impinging vorticity into a system of discrete vortices near the end-wall juncture of the obstruction, with the extensions of the vortices engirdling the obstruction to form "necklace" or "horseshoe" vortices. Recent hydrogen bubble and particle image visualization have shown that as Reynolds number is increased for a laminar approach flow, the flow will become critical, and a destabilization of the necklace vortices results in the development of an azimuthal waviness, or "kinks", in the vortices. These vortex kinks are accentuated by Biot-Savart effects, causing portions of a distorted necklace vortex to make a rapid approach to the surface, precipitating processes of localized, three-dimensional surface interactions. These interactions result in the rapid generation, focussing, and ejection of thin tongues of surface fluid, which rapidly roll-over and appear as hairpin vortices in the junction region. Subsequent amalgamation of these hairpin vortices with the necklace vortices produces a complex transitional-type flow. A presentation of key results from both these studies will be done, emphasizing both the ubiquity of such hairpin-type flow structures in manifold transitional-type flows, and the importance of vortex-surface interactions n the development of hairpin vortices.

Smith, Charles R.↗

Validation of US3D for Capsule Aerodynamics using 05-CA Wind Tunnel Test Data

Several comparisons of computational fluid dynamics to wind tunnel test data are shown for the purpose of code validation. The wind tunnel test, 05-CA, uses a 7.66% model of NASA's Multi-Purpose Crew Vehicle in the 11-foot test section of the Ames Unitary Plan Wind tunnel. A variety of freestream conditions over four Mach numbers and three angles of attack are considered. Test data comparisons include time-averaged integrated forces and moments, time-averaged static pressure ports on the surface, and Strouhal Number. The applicability of the US3D code to subsonic and transonic flow over a bluff body is assessed on a comprehensive data set. With close comparison, this work validates US3D for highly separated flows similar to those examined here.

Schwing, Alan↗

Hot-Film and Hot-Wire Anemometry for a Boundary Layer Active Flow Control Test

Unsteady active flow control (AFC) has been used experimentally for many years to minimize bluff-body drag. This technology could significantly improve performance of rotorcraft by cleaning up flow separation. It is important, then, that new actuator technologies be studied for application to future vehicles. A boundary layer wind tunnel was constructed with a 1ft-x-3ft test section and unsteady measurement instrumentation to study how AFC manipulates the boundary layer to overcome adverse pressure gradients and flow separation. This unsteady flow control research requires unsteady measurement methods. In order to measure the boundary layer characteristics, both hot-wire and hot-film Constant Temperature Anemometry is used. A hot-wire probe is mounted in the flow to measure velocity while a hot-film array lays on the test surface to measure skin friction. Hot-film sensors are connected to an anemometer, a Wheatstone bridge circuit with an output that corresponds to the dynamic flow response. From this output, the time varying flow field, turbulence, and flow reversal can be characterized. Tuning the anemometers requires a fan test on the hot-film sensors to adjust each output. This is a delicate process as several variables drastically affect the data, including control resistance, signal input, trim, and gain settings.

Control Test↗

Helicopter Fuselage Active Flow Control in the Presence of a Rotor

This work extends previous investigations of active flow control for helicopter fuselage drag and download reduction to include the effects of the rotor. The development of the new wind tunnel model equipped with fluidic oscillators is explained in terms of the previous test results. Large drag reductions greater than 20% in some cases were measured during powered testing without increasing, and in some cases decreasing download in forward flight. As confirmed by Particle Image Velocimetry (PIV), the optimum actuator configuration that provided a decrease in both drag and download appeared to create a virtual (fluidic) boat-tail fairing instead of attaching flow to the ramp surface. This idea of a fluidic fairing shifts the focus of 3D separation control behind bluff bodies from controlling/reattaching surface boundary layers to interacting with the wake flow.

Martin, Preston B↗

An Examination of Parameters Affecting Large Eddy Simulations of Flow Past a Square Cylinder

Separated flow over a bluff body is analyzed via large eddy simulations. The turbulent flow around a square cylinder features a variety of complex flow phenomena such as highly unsteady vortical structures, reverse flow in the near wall region, and wake turbulence. The formation of spanwise vortices is often times artificially suppressed in computations by either insufficient depth or a coarse spanwise resolution. As the resolution is refined and the domain extended, the artificial turbulent energy exchange between spanwise and streamwise turbulence is eliminated within the wake region. A parametric study is performed highlighting the effects of spanwise vortices where the spanwise computational domain's resolution and depth are varied. For Re=22,000, the mean and turbulent statistics computed from the numerical large eddy simulations (NLES) are in good agreement with experimental data. Von-Karman shedding is observed in the wake of the cylinder. Mesh independence is illustrated by comparing a mesh resolution of 2 million to 16 million. Sensitivities to time stepping were minimized and sampling frequency sensitivities were nonpresent. While increasing the spanwise depth and resolution can be costly, this practice was found to be necessary to eliminating the artificial turbulent energy exchange.

Mankbadi, M. R.↗

Hypersonic Flight and the Re-Entry Problem: The Twenty-First Wright Brothers Lecture

Up to this point the discussion of the problems of rocket vehicles has been confined to the effects of phenomena which have in the past been important ones for lower speed aircraft and will continue to be important for aircraft of all speeds. Now with considerable extension of both speed and altitude, other phenomena also become important. The nature of some problems will be altered, as a result, and new problems, of course, will be encountered. First, it is well to note that our interest in bluff bodies for ballistic vehicles in particular, and in rounded-nosed bodies generally, has changed our emphasis in aerodynamics. The detached bow waves which occur with such bodies at high supersonic speeds complicate the calculations of the flow-field characteristics. In the present period, much attention is being given to such studies. In addition, at the very high altitudes attained by most of the rocket-craft, the mean free path of air molecules can be of the same order, or long, compared to the dimensions of the vehicles. Thus, slip-flow and free-molecule-flow studies are of interest, particularly for satellite vehicles." The aerodynatnicist must deal with air having unfamiliar states and properties. Second, at hypersonic speeds where, for example, air is greatly decelerated, it may undergo considerable change in composition, the degree of change depending upon many factors. Dissociation of oxygen and nitrogen molecules can occur and, in addition, thermal ionization of many of the constituents. It is naturally to be expected that the convective heat transfer will, as a result, be altered from what it was for the "perfect" gas, and this has been the subject of much recent research effort. Moreover, the decelerated gas becomes capable of radiating energy and the radiative heat transfer must generally be considered for hypersonic vehicles, particularly for long-range ballistic rockets. It is not only the aerodynamic heating problems that are affected. The fact that at very high air temperature the gas becomes electrically conductive introduces new problems in radio wave transmission and reception. In addition, a conducting gas flow can, of course, be influenced by a magnetic field. The study of such flows, which has been termed "magneto gas dynamics," is still in too primitive a state to indicate how important a role it can play, but many interesting possibilities suggest themselves. Third, our experience with airplanes powered by air-breathing engines has naturally been restricted to the stratosphere, or lower. Our ignorance increases with altitude. For rockets, literally, "the sky's the limit," and it is not surprising that a great emphasis has now been placed on obtaining a more thorough understanding of the whole atmosphere. These studies are not aimed at an understanding of the chemical and physical characteristics alone, but also of the occurrence of high-energy particles, from meteors to cosmic rays, and the nature of the problems they will promote.

Wright Brother↗

Experimental Verification of Buffet Calculation Procedure Using Unsteady PSP

Typically a limited number of dynamic pressure sensors are employed to determine the unsteady aerodynamic forces on large, slender aerospace structures. The estimated forces are known to be very sensitive to the number of the dynamic pressure sensors and the details of the integration scheme. This report describes a robust calculation procedure, based on frequency-specific correlation lengths, that is found to produce good estimation of fluctuating forces from a few dynamic pressure sensors. The validation test was conducted on a flat panel, placed on the floor of a wind tunnel, and was subjected to vortex shedding from a rectangular bluff-body. The panel was coated with fast response Pressure Sensitive Paint (PSP), which allowed time-resolved measurements of unsteady pressure fluctuations on a dense grid of spatial points. The first part of the report describes the detail procedure used to analyze the high-speed, PSP camera images. The procedure includes steps to reduce contamination by electronic shot noise, correction for spatial non-uniformities, and lamp brightness variation, and finally conversion of fluctuating light intensity to fluctuating pressure. The latter involved applying calibration constants from a few dynamic pressure sensors placed at selective points on the plate. Excellent comparison in the spectra, coherence and phase, calculated via PSP and dynamic pressure sensors validated the PSP processing steps. The second part of the report describes the buffet validation process, for which the first step was to use pressure histories from all PSP points to determine the "true" force fluctuations. In the next step only a selected number of pixels were chosen as "virtual sensors" and a correlation-length based buffet calculation procedure was applied to determine "modeled" force fluctuations. By progressively decreasing the number of virtual sensors it was observed that the present calculation procedure was able to make a close estimate of the "true" unsteady forces only from four sensors. It is believed that the present work provides the first validation of the buffet calculation procedure which has been used for the development of many space vehicles.

Buffet↗

Flow Energy Piezoelectric Bimorph Nozzle Harvester

There is a need for a long-life power generation scheme that could be used downhole in an oil well to produce 1 Watt average power. There are a variety of existing or proposed energy harvesting schemes that could be used in this environment but each of these has its own limitations. The vibrating piezoelectric structure is in principle capable of operating for very long lifetimes (decades) thereby possibly overcoming a principle limitation of existing technology based on rotating turbo-machinery. In order to determine the feasibility of using piezoelectrics to produce suitable flow energy harvesting, we surveyed experimentally a variety of nozzle configurations that could be used to excite a vibrating piezoelectric structure in such a way as to enable conversion of flow energy into useful amounts of electrical power. These included reed structures, spring mass-structures, drag and lift bluff bodies and a variety of nozzles with varying flow profiles. Although not an exhaustive survey we identified a spline nozzle/piezoelectric bimorph system that experimentally produced up to 3.4 mW per bimorph. This paper will discuss these results and present our initial analyses of the device using dimensional analysis and constitutive electromechanical modeling. The analysis suggests that an order-of-magnitude improvement in power generation from the current design is possible.

transducers↗

Wavenumber-frequency Spectra of Pressure Fluctuations Measured via Fast Response Pressure Sensitive Paint

The recent advancement in fast-response Pressure-Sensitive Paint (PSP) allows time-resolved measurements of unsteady pressure fluctuations from a dense grid of spatial points on a wind tunnel model. This capability allows for direct calculations of the wavenumber-frequency (k-ω) spectrum of pressure fluctuations. Such data, useful for the vibro-acoustics analysis of aerospace vehicles, are difficult to obtain otherwise. For the present work, time histories of pressure fluctuations on a flat plate subjected to vortex shedding from a rectangular bluff-body were measured using PSP. The light intensity levels in the photographic images were then converted to instantaneous pressure histories by applying calibration constants, which were calculated from a few dynamic pressure sensors placed at selective points on the plate. Fourier transform of the time-histories from a large number of spatial points provided k-ω spectra for pressure fluctuations. The data provides first glimpse into the possibility of creating detailed forcing functions for vibro-acoustics analysis of aerospace vehicles, albeit for a limited frequency range.

Pressure Sensitive Paint↗

Computational Fluid Dynamic Optimization of an Experimental Rotating Detonation Rocket Engine Nozzle

A parametric optimization study is performed on the nozzle of a laboratory rotating detonation rocket engine (RDRE) using a three-dimensional computational fluid dynamic simulation. The primary optimization objective is maximum nozzle thrust. The basic nozzle configuration is a shrouded, truncated plug. The fluid in the RDRE chamber leading to the nozzle is choked at its exit so that its cyclic behavior is unaffected by any changes to the nozzle design. Optimization is performed for a single operating point. Parameters varied are the overall nozzle area expansion ratio and the fraction of the expansion area that is provided by the shroud. These two parameters indirectly affect the angle of the plug nozzle cone, and the bluff body area associated with its truncation. Nozzle thrust is evaluated as the difference between the thrust of the RDRE chamber-plus-nozzle combination and that of the chamber alone. The nozzle produces approximately 20% of the total engine thrust. The baseline nozzle is found to perform well, yielding 58.1% of the thrust calculated for a notional ideal RDRE nozzle which can instantaneously change shape to allow isentropic expansion of every fluid element. Optimization improves the performance, bringing the nozzle thrust to 70.0% of the notional ideal, and total engine thrust (chamber-plus-nozzle) to 94% of the ideal.

detonation↗

Computational Fluid Dynamic Optimization of an Experimental Rotating Detonation Rocket Engine Nozzle

A parametric optimization study is performed on the nozzle of a laboratory rotating detonation rocket engine (RDRE) using a three-dimensional computational fluid dynamic simulation. The primary optimization objective is maximum nozzle thrust. The basic nozzle configuration is a shrouded, truncated plug. The fluid in the RDRE chamber leading to the nozzle is choked at its exit so that its cyclic behavior is unaffected by any changes to the nozzle design. Optimization is performed for a single operating point. Parameters varied are the overall nozzle area expansion ratio and the fraction of the expansion area that is provided by the shroud. These two parameters indirectly affect the angle of the plug nozzle cone, and the bluff body area associated with its truncation. Nozzle thrust is evaluated as the difference between the thrust of the RDRE chamber-plus-nozzle combination and that of the chamber alone. The nozzle produces approximately 20% of the total engine thrust. The baseline nozzle is found to perform well, yielding 58.1% of the thrust calculated for a notional ideal RDRE nozzle which can instantaneously change shape to allow isentropic expansion of every fluid element. Optimization improves the performance, bringing the nozzle thrust to 70.0% of the notional ideal, and total engine thrust (chamber-plus-nozzle) to 94% of the ideal.

detonation↗