Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple shocks”

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 109 records · Page 6

Sonic Booms from Aircraft in Maneuvers

Superboom measurements and calculated pressure patterns have been made for fighter aircraft in level accelerated flight and in turn maneuvers. A summary of the main findings of these superboom studies, qualitative comparisons with analysis, and a physical explanation of some of the observed phenomena are presented. From 2 to 6 booms were observed (compared with normally overserved 2 booms for steady flight), and pressure-buildup factors of from 2 to 4 were measured (depending on the type of maneuver and the location of the observer).The relative positions of the multiple shocks were predicted within 1,000 feet, and the location of the initial superboom impact was found to be predictable within about ±2 miles when the temperature profile and airplane position information were available.

SONIC BOOM↗

Experimental investigation of hypersonic buzz on a high cross-range shuttle configuration.

An experimental study confirming the existence of hypersonic buzz in the separated region of a space shuttle's downward deflected control surface is described. Wind tunnel tests in helium on a 1:100 scale model of a shuttle type configuration at M = 17.5 clearly demonstrated the occurrence of large scale pressure oscillations throughout a substantial range of angles of attack and control surface deflections. Measurements of surface pressure fluctuations along with high speed schlieren photography served to delineate the boundaries of the phenomenon and to identify the physical parameters necessary for the onset of a buzz condition. Prominent spectral energy peaks were also observed in the surface pressures that could be directly correlated with the instability. These conditions were found to be strongly dependent on the presence of multiple shock interactions and the position of these interactions relative to the location of a reattaching shear layer.

Goldman, R. L.↗

Wing analysis using a transonic potential flow computational method

The ability of the method to compute wing transonic performance was determined by comparing computed results with both experimental data and results computed by other theoretical procedures. Both pressure distributions and aerodynamic forces were evaluated. Comparisons indicated that the method is a significant improvement in transonic wing analysis capability. In particular, the computational method generally calculated the correct development of three-dimensional pressure distributions from subcritical to transonic conditions. Complicated, multiple shocked flows observed experimentally were reproduced computationally. The ability to identify the effects of design modifications was demonstrated both in terms of pressure distributions and shock drag characteristics.

Henne, P. A.↗

Shuttle system ascent aerodynamic and plume heating

The shuttle program provided a challenge to the aerothermodynamicist due to the complexity of the flow field around the vehicle during ascent, since the configuration causes multiple shock interactions between the elements. Wind tunnel tests provided data for the prediction of the ascent design heating environment which involves both plume and aerodynamic heating phenomena. The approach for the heating methodology based on ground test firings and the use of the wind tunnel data to formulate the math models is discussed.

Foster, L. D.↗

A multiple type-II burst associated with a coronal transient and its MHD simulation

A large coronal transient took place on May 8, 1981. The transient was related to an M7.7/2B flare and was associated with at least two coronal type-II bursts. The velocities of the type-II bursts were in the range 1100-1800 km/s, in excess of the transient velocity of 1000 km/s. Two dimensional positions of the type-II radio sources are available from both the Clark Lake and the Culgoora Radio Observatories. Two dimensional MHD simulations of the event are carried out, taking into account the observed velocity, position, and size of the type-II bursts. The multiple shocks observed during the event and their interaction are simulated, and results of the simulation are discussed.

Gergely, T. E.↗

Application of NCOREL to conical multi-finned and multi-faceted configurations

A more versatile analytic conformal mapping approach for grid generation is implemented in a full potential supersonic flow code (NCOREL). Configurations such as multi-finned bodies and wings with vertical tails cannot be treated using a single conformal mapping transformation for grid generation. Instead, a series of analytic conformal mappings are used in progression to generate grids that are capable of resolving the complex multiple shock flow fields that exist about these configurations. Aplications of the grid generation techniques are presented for a variety of cross sections along with their conical full potential flow solutions.

Siclari, M. J.↗

Numerical simulations of adiabatic axisymmetric accretion flow. I - A new mechanism for the formation of jets

Numerical simulations of the uniform axisymmetric flow past a gravitating sphere have been studied. It is found that the structure of the flow is extremely sensitive to the boundary condition at the surface of the gravitating object. For the case in which the boundary is totally absorbing, a steady state flow is reached. However, for a boundary which is not totally absorbing, steady state flows are not obtained. The morphology of the flow is also sensitive to the Mach number at infinity and to the ratio of the free-fall velocity at the surface of the gravitating object to the flow velocity at inifinity. A new mechanism for the formation of jets is identified in which a fraction of the accretion energy is tapped to drive an anisotropic supersonic outflow with collimation provided by a combination of the inertia of matter which surrounds the beam and the development of multiple shock structures.

Fryxell, B. A.↗

Schlieren studies of compressibility effects on dynamic stall of airfoils in transient pitching motion

Compressibility effects on the flowfield of an airfoil executing rapid transient pitching motion from 0 - 60 degrees over a wide range of Mach numbers and pitching rates were studied using a stroboscopic schlieren flow visualization technique. The studies have led to the first direct experimental documentation of multiple shocks on the airfoil upper surface flow for certain conditions. Also, at low Mach numbers, additional coherent vortical structures were found to be present along with the dynamic stall vortex, whereas at higher Mach numbers, the flow was dominated by a single vortex. The delineating Mach number for significant compressibility effects was 0.3 and the dynamic stall process was accelerated by increasing the Mach number above that value. Increasing the pitch rate monotonically delayed stall to angles of attack as large as 27 degrees.

Chandrasekhara, M. S.↗

Idealized tip-to-tail waverider model

The flow field of an idealized cone-derived waverider is axisymmetric. This forebody feature is preserved for the rest of the vehicle, including the inlet, cowl, combustor, and nozzle. There is thus an inviscid, tip-to-tail model in which both the external and internal flows are axisymmetric. The assumption of axial symmetry provides a major simplification for the analysis and allows for a systematic integration of the propulsion unit with the aerodynamics. The code is an initial formulation that provides only the most basic engineering data, such as lift, thrust, drag, and fuel consumption for a point-designed vehicle that may be cruising at a low hypersonic Mach number. The user may specify flight altitude and Mach number, a multiple shock configuration for the inlet, a few basic geometric parameter, H2 or CH4 as fuel, the fuel/air ratio, etc. A new design concept is used for the nozzle that avoids shock waves, minimizes the nozzle length, and may maximize its thrust. A general description of the model is provided with emphasis on the design of the nozzle and fins. Preliminary results are presented that compare the cruise flight range using H2 or CH4 as the fuel.

Emanuel, G.↗

Cratering mechanics on Venus - Pressure enhancement by the atmospheric 'ocean'

The impedance match technique and EOSs of equations of state (EOSs) of geologically relevant materials are used to investigate cratering mechanics on Venus, specifically, the coupling of impactor kinetic energy and momentum into the target surface. These EOSs are modified to account for multiple shocks. Peak impact pressures from both first reflection and later reverberations are determined. These are compared to values obtained using an atmosphereless model, and the differences between and implications for atmosphere-affected and atmosphereless impacts are discussed.

Brackett, Robert A.↗

When clusters collide - A numerical Hydro/N-body simulation of merging galaxy clusters

A 3D numerical simulation of two merging clusters of galaxies, using a hybrid Hydro/N-body code, is presented. The hydrodynamics of the code is solved by an Eulerian finite difference method. Initial results disclose that the X-ray emission of the dominant cluster becomes elongated and broadened; heating occurs at the core of the dominant cluster as a result of multiple shocks, and high velocity gas motions within the intracluster medium. It is predicted that clusters which have undergone recent mergers and do not have cooling flows will have high peculiar gas velocities and that the shocks and turbulence generated during the merger may power cluster-wide radio halos. Prolonged high-velocity gas motions through the dominant cluster core possibly play a major role in the formation and shaping of wide-angle tailed radio sources associated with central dominant galaxies. The N-body component of the simulation reveals the subcluster to be dispersed as it passes through the dominant cluster.

Roettiger, Kurt↗

Interferometric investigations of compressible dynamic stall over a transiently pitching airfoil

The dynamic stall flow field over NACA 0012 airfoil pitching transiently from 0 - 60 at a constant rate under compressible flow conditions has been studied using the real-time technique of point diffraction interferometry. This investigation using nonintrusive diagnostics provides a quantitative description of the overall flow field, including the finer details of dynamic stall vortex formation, growth and the concomitant changes in the pressure distribution. Analysis of several hundred interferograms obtained for a range of flow conditions shows that the peak leading edge suction pressure coefficient that stall is nearly constant for a given free stream Mach number at all nondimensional pitch rates. Also, this value is below that seen in steady flow at static stall for the same Mach number, indicating that dynamic effects significantly effect the separation behavior. Further, for a given Mach number, the dynamic stall vortex seems to form rapidly at nearly the same angle of attack for all pitch rates studied. As the vortex is shed, it induces an anti-clockwise trailing edge vortex, which grows in a manner similar to that of a starting vortex. The measured peak suction pressure coefficient drops as the free stream Mach number increases. For free stream Mach numbers above 0.4, small multiple shocks appear near the leading edge.

Chandrasekhara, M. S.↗

Compressibility effects on dynamic stall of airfoils undergoing rapid transient pitching motion

The research was carried out in the Compressible Dynamic Stall Facility, CDSF, at the Fluid Mechanics Laboratory (FML) of NASA Ames Research Center. The facility can produce realistic nondimensional pitch rates experienced by fighter aircraft, which on model scale could be as high as 3600/sec. Nonintrusive optical techniques were used for the measurements. The highlight of the effort was the development of a new real time interferometry method known as Point Diffraction Interferometry - PDI, for use in unsteady separated flows. This can yield instantaneous flow density information (and hence pressure distributions in isentropic flows) over the airfoil. A key finding is that the dynamic stall vortex forms just as the airfoil leading edge separation bubble opens-up. A major result is the observation and quantification of multiple shocks over the airfoil near the leading edge. A quantitative analysis of the PDI images shows that pitching airfoils produce larger suction peaks than steady airfoils at the same Mach number prior to stall. The peak suction level reached just before stall develops is the same at all unsteady rates and decreases with increase in Mach number. The suction is lost once the dynamic stall vortex or vortical structure begins to convect. Based on the knowledge gained from this preliminary analysis of the data, efforts to control dynamic stall were initiated. The focus of this work was to arrive at a dynamically changing leading edge shape that produces only 'acceptable' airfoil pressure distributions over a large angle of attack range.

Chandrasekhara, M. S.↗

Interferometric investigations of compressible dynamic stall over a transiently pitching airfoil

The compressible dynamic stall flowfield over a NACA 0012 airfoil transiently pitching from 0 to 60 deg at a constant rate under compressible flow conditions has been studied using real-time interferometry. A quantitative description of the overall flowfield, including the finer details of dynamic stall vortex formation, growth, and the concomitant changes in the airfoil pressure distribution, has been provided by analyzing the interferograms. For Mach numbers above 0.4, small multiple shocks appear near the leading edge and are present through the initial stages of dynamic stall. Dynamic stall was found to occur coincidentally with the bursting of the separation bubble over the airfoil. Compressibility was found to confine the dynamic stall vortical structure closer to the airfoil surface. The measurements show that the peak suction pressure coefficient drops with increasing freestream Mach number, and also it lags the steady flow values at any given angle of attack. As the dynamic stall vortex is shed, an anti-clockwise vortex is induced near the trailing edge, which actively interacts with the post-stall flow.

Chandrasekhara, M. S.↗

A Structured Grid Based Solution-Adaptive Technique for Complex Separated Flows

The objective of this work was to enhance the predictive capability of widely used computational fluid dynamic (CFD) codes through the use of solution adaptive gridding. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. In order to study the accuracy and efficiency improvements due to the grid adaptation, it is necessary to quantify grid size and distribution requirements as well as computational times of non-adapted solutions. Flow fields about launch vehicles of practical interest often involve supersonic freestream conditions at angle of attack exhibiting large scale separate vortical flow, vortex-vortex and vortex-surface interactions, separated shear layers and multiple shocks of different intensity. In this work, a weight function and an associated mesh redistribution procedure is presented which detects and resolves these features without user intervention. Particular emphasis has been placed upon accurate resolution of expansion regions and boundary layers. Flow past a wedge at Mach=2.0 is used to illustrate the enhanced detection capabilities of this newly developed weight function.

Thornburg, Hugh↗

Vapor Flow Patterns During a Start-Up Transient in Heat Pipes

The vapor flow patterns in heat pipes are examined during the start-up transient phase. The vapor core is modelled as a channel flow using a two dimensional compressible flow model. A nonlinear filtering technique is used as a post process to eliminate the non-physical oscillations of the flow variables. For high-input heat flux, multiple shock reflections are observed in the evaporation region. The reflections cause a reverse flow in the evaporation and circulations in the adiabatic region. Furthermore, each shock reflection causes a significant increase in the local pressure and a large pressure drop along the heat pipe.

Issacci, F.↗

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗

Estimation of Supersonic Stage Separation Aerodynamics of Winged-Body Launch Vehicles Using Response Surface Methods

A collection of statistical and mathematical techniques referred to as response surface methodology was used to estimate the longitudinal stage separation aerodynamic characteristics of a generic, bimese, winged multi-stage launch vehicle configuration using data obtained on small-scale models at supersonic speeds in the NASA Langley Research Center Unitary Plan Wind Tunnel. The simulated Mach 3 staging was dominated by multiple shock wave interactions between the orbiter and booster vehicles throughout the relative spatial locations of interest. This motivated a partitioning of the overall inference space into several contiguous regions within which the separation aerodynamics were presumed to be well-behaved and estimable using cuboidal and spherical central composite designs capable of fitting full second-order response functions. The primary goal was to approximate the underlying overall aerodynamic response surfaces of the booster vehicle in belly-to-belly proximity to the orbiter vehicle using relatively simple, lower-order polynomial functions that were piecewise-continuous across the full independent variable ranges of interest. The quality of fit and prediction capabilities of the empirical models were assessed in detail, and the issue of subspace boundary discontinuities was addressed. The potential benefits of augmenting the central composite designs to full third order using computer-generated D-optimality criteria were also evaluated. The usefulness of central composite designs, the subspace sizing, and the practicality of fitting low-order response functions over a partitioned inference space dominated by highly nonlinear and possibly discontinuous shock-induced aerodynamics are discussed.

Erickson, Gary E.↗