Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “moving shock”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Magnetic Fields in the Termination Shock, Heliosheath and Solar Wind

This paper will discuss recent magnetic field observations of the termination shock (TS), heliosheath (HSH), and supersonic solar wind (SW) made by Voyager 1 (V1) and Voyager 2 (V2). Voyager 2 crossed the TS at least 5 times at 83.7 AU and 27.5deg South heliographic latitude during 2007/242.14 to 244.80. The structure of the shock, observed on three of the crossings, was highly variable, ranging from a simple quasi-perpendicular supercritical shock to a complex profile suggesting reformation. Two merged interaction regions were observed in the SW upstream of the TS, and the magnetic field increased to relatively high values during the month prior to crossing the TS. Voyager 2 observations in the HSH will also be discussed. Voyager 1 observed a small radial gradient of the magnetic field strength B in the HSH. It is expected that B should increase from the TS to heliopause, but the available predictions of the profile are not adequate to compare with the observations. Voyager 1 has entered a unipolar region, as the heliospheric current sheet moved below the latitude of V1.

Burlaga, L. F.↗
DOE OSTI2020

Parallel Implicit Hydrodynamics for High Explosive Burn Calculations

High explosives in hostile environments will require calculational capabilities that model processes, which evolve on timescales from minutes to nanoseconds. Eventually the HE will begin to move metal. To handle this temporal evolution an implicit hydrodynamics coupled to the chemical release of the HE energy is required. In addition, the use of chemical kinetics to model the transition from, the initially, slow heating of a confined high explosive through to deflagration and on to detonation requires many computational zones to model high explosive engineering systems. This requirement means that a fully parallel implicit hydrodynamics is essential. In this paper we present the calculation of a nonlinear matrix equation for the advanced particle pressure that has been made parallel and implemented in our AMR code, BABBO. This new parallel implicit hydrodynamics has been applied to a cookoff problem, as well as, one and two dimensional shock problems. Results are presented and discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2004

Global Dayside Ionospheric Uplift and Enhancement Associated with Interplanetary Electric Fields

The interplanetary shock/electric field event of 5-6 November 2001 is analyzed using ACE interplanetary data. The consequential ionospheric effects are studied using GPS receiver data from the CHAMP and SAC-C satellites and altimeter data from the TOPEX/ Poseidon satellite. Data from ~100 ground-based GPS receivers as well as Brazilian Digisonde and Pacific sector magnetometer data are also used. The dawn-to-dusk interplanetary electric field was initially ~33 mV/m just after the forward shock (IMF BZ = -48 nT) and later reached a peak value of ~54 mV/m 1 hour and 40 min later (BZ = -78 nT). The electric field was ~45 mV/m (BZ = -65 nT) 2 hours after the shock. This electric field generated a magnetic storm of intensity DST = -275 nT. The dayside satellite GPS receiver data plus ground-based GPS data indicate that the entire equatorial and midlatitude (up to +/-50(deg) magnetic latitude (MLAT)) dayside ionosphere was uplifted, significantly increasing the electron content (and densities) at altitudes greater than 430 km (CHAMP orbital altitude). This uplift peaked ~2 1/2 hours after the shock passage. The effect of the uplift on the ionospheric total electron content (TEC) lasted for 4 to 5 hours. Our hypothesis is that the interplanetary electric field ''promptly penetrated'' to the ionosphere, and the dayside plasma was convected (by E x B) to higher altitudes. Plasma upward transport/convergence led to a ~55-60% increase in equatorial ionospheric TEC to values above ~430 km (at 1930 LT). This transport/convergence plus photoionization of atmospheric neutrals at lower altitudes caused a 21% TEC increase in equatorial ionospheric TEC at ~1400 LT (from ground-based measurements). During the intense electric field interval, there was a sharp plasma ''shoulder'' detected at midlatitudes by the GPS receiver and altimeter satellites. This shoulder moves equatorward from -54(deg) to -37(deg) MLAT during the development of the main phase of the magnetic storm. We presume this to be an ionospheric signature of the plasmapause and its motion. The total TEC increase of this shoulder is ~80%. Part of this increase may be due to a "superfountain effect." The dayside ionospheric TEC above ~430 km decreased to values ~45% lower than quiet day values 7 to 9 hours after the beginning of the electric field event. The total equatorial ionospheric TEC decrease was ~16%. This decrease occurred both at midlatitudes and at the equator. We presume that thermospheric winds and neutral composition changes produced by the storm-time Joule heating, disturbance dynamo electric fields, and electric fields at auroral and subauroral latitudes are responsible for these decreases.

ionosphere↗
NASA NTRS2005

STEREO WAVES Capabilities for Studying Initiation and Early-time Dynamics of Solar Eruptions

In 2006, NASA will launch the twin STEREO spacecraft from Kennedy Space Center into a pair of heliocentric orbits near 1 AU such that the spacecraft will move away from Earth (ahead and behind) at about 22 degrees per year. The purposes of the STEREO Mission are to understand the causes and mechanisms of coronal mass ejection (CME) initiation and to follow the propagation of CMEs through the heliosphere. Additionally, STEREO will study the mechanisms and sites of solar energetic particle (SEP) acceleration and determine 3-D time-dependent traces of the magnetic topology, temperature, density and velocity of the solar wind between the sun and Earth. To accomplish these goals, each STEREO spacecraft will be equipped with set of optical and particles and fields instruments including the STEREO WAVES (SWAVES) instrument which will use radio waves to track the location of CME-driven shocks (via type I1 bursts) and the 3-D topology of open field lines along which energetic particles flow (via the ubiquitous type I11 bursts). Type 11 bursts very often commence with a series of special type 111 bursts (called SA or type 111-L bursts) that likely coincide with CME liftoff time, thus SWAVES should be able to determine this time to within 15 sec. It is also known that the occurrence of SEP events is usually accompanied by type I1 radio bursts at decametric wavelengths as well as strong type III bursts at all wavelengths. SWAVES will be able to determine the initiation of these bursts to within 15 sec, and from the simultaneous measurements from the two spacecraft, should be able to triangulate their source locations. The utility of radio observations and the capabilities of SWAVES will be illustrated by showing a number of examples using the similar Wind WAVES instrument in combination with SOH0 coronagraph and RHESSI high energy X-ray/gamma ray observations.

Kaiser, M. L.↗
NASA NTRS1989

The solar wind - Advances in our knowledge through two solar cycles

As the Pioneer and Voyager spacecraft have moved outward they have gradually unfolded a view of distant regions of the heliosphere. Information on the solar wind velocity, density and temperature as a function of distance out to more than 40 AU has been gathered. Meanwhile the description of the solar wind has evolved. Long-standing questions on the sources of the wind causing geomagnetic activity were clarified by the discovery of coronal holes and coronal mass ejections. The propagation of the resultant solar wind disturbances through the heliosphere has been studied using both observations and models. Plasma physical processes have been studied. This review focuses on the development of the concepts that have been used to describe the solar wind in the three dimensional heliosphere over the last two solar cycles. Collisionless shocks, transient disturbances in space, disturbance propagation and the distant solar wind are discussed.

Feynman, Joan↗
DOE OSTI2023

CO Laser Absorption Measurements During Syngas Combustion at High Pressure

Abstract Syngas is a desirable fuel for combustion in the Allam-Fetvedt cycle, which involves combustion under supercritical-CO2 conditions. While some work has been conducted in collecting ignition delay times (IDT) at the extreme pressures required by these systems, significant model deficiencies remain. Additionally, considerable barriers in terms of nonideal gas dynamic effects have been shown for these experiments in shock tubes. Further investigation into the fundamental combustion kinetics of H2/CO/CO2 mixtures is required. Time-resolved speciation measurements for target species have been shown to better aid in improving the understanding of underlying chemical kinetics than global ignition delay time measurements. Therefore, laser absorption measurements of CO were measured behind reflected shock waves during combustion of syngas at 5 and 10 bar and temperatures between 1080 and 2100 K. The mixtures investigated utilized H2-to-CO ratios of 1:1 and 1:4, respectively, each at stoichiometric conditions, allowing for discussions of the effect of initial fuel composition. A ratio of fuel to CO2 of 1:2 was also utilized to represent commercially available syngas. The mixtures were diluted in helium and argon (20% He, 76.5% Ar) to minimize thermal effects and to expedite CO thermal relaxation during the experiment. The resulting CO time histories were then compared to modern chemical kinetics mechanisms, and disagreement is seen for this system, which is assumed to be fairly well known. This study elucidates particular chemistry that needs improvement in moving toward a better understanding of syngas combustion at elevated pressures.

Engineering↗
NASA NTRS2012

Study of Convective Flow Effects in Endwall Casing Treatments in Transonic Compressor Rotors

The unsteady convective flow effects in a transonic compressor rotor with a circumferential-groove casing treatment are investigated in this paper. Experimental results show that the circumferential-groove casing treatment increases the compressor stall margin by almost 50% for the current transonic compressor rotor. Steady flow simulation of the current casing treatment, however, yields only a 15% gain in stall margin. The flow field at near-stall operation is highly unsteady due to several self-induced flow phenomena. These include shock oscillation, vortex shedding at the trailing edge, and interaction between the passage shock and the tip clearance vortex. The primary focus of the current investigation is to assess the effects of flow unsteadiness and unsteady flow convection on the circumferential-groove casing treatment. Unsteady Reynolds-averaged Navier-Stokes (URANS) and Large Eddy Simulation (LES) techniques were applied in addition to steady Reynolds-averaged Navier-Stokes (RANS) to simulate the flow field at near-stall operation and to determine changes in stall margin. The current investigation reveals that unsteady flow effects are as important as steady flow effects on the performance of the circumferential grooves casing treatment in extending the stall margin of the current transonic compressor rotor. The primary unsteady flow mechanism is unsteady flow injection from the grooves into the main flow near the casing. Flows moving into and out of the grooves are caused due to local pressure difference near the grooves. As the pressure field becomes transient due to self-induced flow oscillation, flow injection from the grooves also becomes unsteady. The unsteady flow simulation shows that this unsteady flow injection from the grooves is substantial and contributes significantly to extending the compressor stall margin. Unsteady flows into and out of the grooves have as large a role as steady flows in the circumferential grooves. While the circumferential-groove casing treatment seems to be a steady flow device, unsteady flow effects should be included to accurately assess its performance as the flow is transient at near-stall operation.

Hah, Chunill↗
DOE OSTI2020

A reconstructed discontinuous Galerkin method for compressible flows in Lagrangian formulation

In this work, a high-order accurate reconstructed discontinuous Galerkin (rDG) method is developed for solving two-dimensional hydrodynamic problems in cell-centered updated Lagrangian formulation. This method is the Lagrangian limit of the unsplit rDG-ALE formulation, and is obtained by assuming the equality of the grid velocity to the fluid velocity only at cell boundaries. The conservative variables and the Taylor basis defined on the time-dependent moving mesh, provide the piece-wise polynomial expansion in the updated Lagrangian formulation. A multi-directional nodal Riemann solver is implemented for computing the grid velocity at the vertices and the numerical flux at the cell boundaries. A characteristic limiting procedure is extended from the primitive variable version to the conservative variable version, and its performance is compared with the limiter on physical variables. A number of benchmark test cases are conducted to assess the accuracy, robustness, and non-oscillatory property of the DG(P0), DG(P1) and rDG(P1P2) methods. The numerical experiments demonstrate that the developed rDG method is able to attain the designed order of accuracy and the characteristic limiting procedure outperforms the limiter on physical variables in terms of the monotonicity and symmetry preservation for shock problems.

42 ENGINEERING↗
DOE OSTI2025

A Bayesian approach to time-domain photonic Doppler velocimetry analysis

Photonic Doppler velocimetry (PDV) is an established technique for measuring the velocities of fast-moving surfaces in high-energy-density experiments. In the standard approach to PDV analysis, the short-time Fourier transform (STFT) is used to generate a spectrogram from which the velocity history of the target is inferred. The user chooses the form, duration, and separation of the window function. Here, in this study, we present a Bayesian approach to infer the velocity directly from the PDV oscilloscope trace, without using the spectrogram for analysis. This is clearly a difficult inference problem due to the highly periodic nature of the data, but we find that with carefully chosen prior distributions for the model parameters, we can accurately recover the injected velocity from synthetic data. We validate this method using PDV data collected at the STAR two-stage light gas gun at Sandia National Laboratories, recovering shock-front velocities in quartz that are consistent with those inferred using the STFT-based approach and are interpolated across regions of low signal-to-noise data. Although this method does not rely on the same user choices as the STFT, we caution that it can be prone to misspecification if the chosen model is not sufficient to capture the velocity behavior. Analysis using posterior predictive checks can be used to establish whether a better model is required, although more complex models come with additional computational cost, often taking more than several hours to converge when sampling the Bayesian posterior. We, therefore, recommend it be viewed as a complementary method to that of the STFT-based approach.

Allison, James R. [First Light Fusion Ltd., Yarnto↗
NASA NTRS2001

Jet Stability and the Generation of Superluminal and Stationary Components

We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet results in the formation of multiple conical shocks behind the main perturbation. These trailing components can be easily distinguished because they appear to be released from the primary superluminal component instead of being ejected from the core. Their oblique nature should also result in distinct polarization properties. Those appearing closer to the core show small apparent motions and a very slow secular decrease in brightness and could be identified as stationary components. Those appearing farther downstream are weaker and can reach superluminal apparent motions. The existence of these trailing components indicates that not all observed components necessarily represent major perturbations at the jet inlet; rather, multiple emission components can be generated by a single disturbance in the jet. While the superluminal component associated with the primary perturbation exhibits a rather stable pattern speed, trailing components have velocities that increase with distance from the core but move at less than the jet speed. The trailing components exhibit motion and structure consistent with the triggering of pinch modes by the superluminal component. The increase in velocity of the trailing components is an indirect consequence of the acceleration of the expanding fluid, which is assumed to be relativistically hot; if observed, such accelerations would therefore favor an electron-positron (as opposed to proton rest mass) dominated jet.

Agudo, Ivan↗
NASA NTRS1987

Finite element models of wire rope for vibration analysis

The usefulness of wire rope in shock and vibration isolation is briefly reviewed and its modeling for the purpose of vibration analysis is addressed. A model of a nominally straight segment of wire rope is described in which the rope structure is represented by a maiden, or central, strand of wire with one (or more) strand(s) wrapped around it in a helix (helices). The individual strands are modeled using finite elements and MSC NASTRAN. Small linear segments of each wire are modeled mathematically by dividing them lengthwise into triangular prisms representing each prism by a solid NASTRAN element. To model pretensioning and allow for extraction of internal force information from the NASTRAN model, the wound strands are connected to the maiden strand and each other using spring (scalar elastic) elements. Mode shapes for a length of wire rope with one and fixed to a moving base and the other attached to a point mass, are presented. The use of the NASTRAN derived mode shapes to approximate internal normal forces in equations of motion for vibration analyses is considered.

Cochran, J. E., Jr.↗
NASA NTRS2014

Spectroscopic Measurement Techniques for Aerospace Flows

The conditions that characterize aerospace flows are so varied, that a single diagnostic technique is not sufficient for its measurement. Fluid dynamists use knowledge of similarity to help categorize and focus on different flow conditions. For example, the Reynolds number represents the ratio of inertial to viscous forces in a flow. When the velocity scales, length scales, and gas density are large and the magnitude of the molecular viscosity is low, the Reynolds number becomes large. This corresponds to large scale vehicles (e.g Airbus A380), fast moving objects (e.g. artillery projectiles), vehicles in dense fluids (e.g. submarine in water), or flows with low dynamic viscosity (e.g. skydiver in air). In each of these cases, the inertial forces dominate viscous forces, and unsteady turbulent fluctuations in the flow variables are observed. In contrast, flows with small length scales (e.g. dispersion of micro-particles in a solid rocket nozzle), slow moving objects (e.g. micro aerial vehicles), flows with low density gases (e.g. atmospheric re-entry), or fluids with a large magnitude of viscosity (e.g. engine coolant flow), all have low Reynolds numbers. In these cases, viscous forces become very important and often the flows can be steady and laminar. The Mach number, which is the ratio of the velocity to the speed of sound in the medium, also helps to differentiate types of flows. At very low Mach numbers, acoustic waves travel much faster than the object, and the flow can be assumed to be incompressible (e.g. Cessna 172 aircraft). As the object speed approaches the speed of sound, the gas density can become variable (e.g. flow over wing of Learjet 85). When the object speed is higher than the speed of sound (Ma > 1), the presences of shock waves and other gas dynamic features can become important to the vehicle performance (e.g. SR-71 Blackbird). In the hypersonic flow regime (Ma > 5), large changes in temperature begin to affect flow properties, causing real-gas effects to occur (e.g. X-43 Scramjet). At even higher Mach numbers, chemistry and nonequilibrium effects come into play (e.g. Startdust re-entry capsule), further complicating the measurement. These limits can be predicted by calculating the ratio of chemical and thermal relaxation time to the flow time scales. Other non-dimensional numbers can be used to further differentiate types of aerospace flows.

Danehy, Paul M.↗
NASA NTRS2007

Dog-Bone Horns for Piezoelectric Ultrasonic/Sonic Actuators

A shape reminiscent of a dog bone has been found to be superior to other shapes for mechanical-amplification horns that are components of piezoelectrically driven actuators used in a series of related devices denoted generally as ultrasonic/sonic drill/corers (USDCs). The first of these devices was reported in Ultrasonic/Sonic Drill/Corers With Integrated Sensors (NPO-20856), NASA Tech Briefs, Vol. 25, No. 1 (January 2001), page 38. The dog-bone shape was conceived especially for use in a more recent device in the series, denoted an ultrasonic/ sonic gopher, that was described in Ultrasonic/Sonic Mechanisms for Drilling and Coring (NPO-30291), NASA Tech Briefs, Vol. 27, No. 9 (September 2003), page 65. The figure shows an example of a dog-bone-shaped horn and other components of an ultrasonic gopher. Prerequisite to a meaningful description of this development is an unavoidably lengthy recapitulation of the principle of operation of a USDC and, more specifically, of the ultrasonic/sonic gopher as described previously in NASA Tech Briefs. The ultrasonic actuator includes a stack of piezoelectric rings, the horn, a metal backing, and a bolt that connects the aforementioned parts and provides compressive pre-strain to the piezoelectric stack to prevent breakage of the rings during extension. The stack of piezoelectric rings is excited at the resonance frequency of the overall ultrasonic actuator. Through mechanical amplification by the horn, the displacement in the ultrasonic vibration reaches tens of microns at the tip of the horn. The horn hammers an object that is denoted the free mass because it is free to move longitudinally over a limited distance between hard stops: The free mass bounces back and forth between the ultrasonic horn and a tool bit (a drill bit or a corer). Because the longitudinal speed of the free mass is smaller than the longitudinal speed of vibration of the tip of the horn, contact between the free mass and the horn tip usually occurs at a phase of the vibration favorable to transfer of momentum from the horn to the free mass. Therefore, the free mass picks up momentum and is accelerated back to the tool bit. Upon impact of the free mass on the tool bit, momentum is transferred to the tool bit. The impacts of the free mass on the tool bit repeat at a sonic frequency that can range from tens of hertz to about 1 kHz. The shock waves caused by the impacts of the free mass propagate to the interface between the tool bit and the medium (typically, rock, ice, or other brittle material) to be drilled or cored. The medium becomes fractured when its ultimate strain is exceeded at the medium/tool-bit interface. This concludes the description of the principle of operation.

Sherrit, Stewart↗
NASA NTRS1994

Magnetic shaping of planetary nebulae and other stellar wind bubbles

As in the case of the solar wind, the magnetic field in the wind from a magnetized, rotating star becomes increasingly toroidal with distance from the star. The strength of the magnetic field can be characterized by sigma, the ratio of toroidal magnetic energy density to kinetic energy density in the equatorial plane of the wind. A fast wind shocks against the external medium and creates a bubble whose volume is dominated by shocked gas. The toroidal magnetic field increases in the shocked bubble and can dominate the thermal pressure. Because of the low velocities in the bubble, hydrostatic equilibrium is a good approximation and allows the calculation of the thermal and magnetic pressure in the bubble, as in the model of Begelman and Li (1992) for the Crab Nebula. The structure, which is axisymmetric and extended in the polar direction, depends on two parameters: sigma nu(sub w)/w(sub 0), where nu(sub w) is the wind velocity and w(sub 0) is the shell velocity in the polar direction, and lambda = nu(sub a)/w(sub 0), where nu(sub a) is the velocity of the slow wind. For small values of lambda, there is a cusp in the shell in the equatorial plane, i.e., there is an equatorial ring. For larger values of lambda, the maximum of the surface density moves away from the equator i.e., a double ring structure. Our models should apply to planetary nebulae, if their central stars are sufficiently magnetized; the calculated shapes do resemble the observed shapes of planetaries. In all cases, our model predicts that X-ray emission from the bubble is concentrated toward the polar axis. Finally, we briefly discuss the asymmetry of the Crab Nebula and 3C 58.

Chevalier, Roger A.↗
NASA NTRS1991

J-85 jet engine noise measured in the ONERA S1 wind tunnel and extrapolated to far field

Noise from a J-85 turbojet with a conical, convergent nozzle was measured in simulated flight in the ONERA S1 Wind Tunnel. Data are presented for several flight speeds up to 130 m/sec and for radiation angles of 40 to 160 degrees relative to the upstream direction. The jet was operated with subsonic and sonic exhaust speeds. A moving microphone on a 2 m sideline was used to survey the radiated sound field in the acoustically treated, closed test section. The data were extrapolated to a 122 m sideline by means of a multiple-sideline source-location method, which was used to identify the acoustic source regions, directivity patterns, and near field effects. The source-location method is described along with its advantages and disadvantages. Results indicate that the effects of simulated flight on J-85 noise are significant. At the maximum forward speed of 130 m/sec, the peak overall sound levels in the aft quadrant were attentuated approximately 10 dB relative to sound levels of the engine operated statically. As expected, the simulated flight and static data tended to merge in the forward quadrant as the radiation angle approached 40 degrees. There is evidence that internal engine or shock noise was important in the forward quadrant. The data are compared with published predictions for flight effects on pure jet noise and internal engine noise. A new empirical prediction is presented that relates the variation of internally generated engine noise or broadband shock noise to forward speed. Measured near field noise extrapolated to far field agrees reasonably well with data from similar engines tested statically outdoors, in flyover, in a wind tunnel, and on the Bertin Aerotrain. Anomalies in the results for the forward quadrant and for angles above 140 degrees are discussed. The multiple-sideline method proved to be cumbersome in this application, and it did not resolve all of the uncertainties associated with measurements of jet noise close to the jet. The simulation was complicated by wind-tunnel background noise and the propagation of low frequency sound around the circuit.

Soderman, Paul T.↗
NASA NTRS1994

Head-on parallel blade-vortex interaction

An experimental and computational study was carried out to investigate the parallel head-on blade-vortex interaction (BVI) and its noise generation mechanism. A shock tube, with an enlarged test section, was used to generate a compressible starting vortex which interacted with a target airfoil. The dual-pulsed holographic interferometry (DPHI) technique and airfoil surface pressure measurements were employed to obtain quantitative flow data during the BVI. A thin-layer Navier-Stokes code (BV12D), with a high-order upwind-biased scheme and a multizonal grid, was also used to simulate numerically the phenomena occurring in the head-on BVI. The detailed structure of a convecting vortex was studied through independent measurements of density and pressure distributions across the vortex center. Results indicate that, in a strong head-on BVI, the opposite pressure peaks are generated on both sides of the leading edge as the vortex approaches. Then, as soon as the vortex passes by the leading edge, the high-pressure peak suddenly moves toward the low-peak-reducing in magnitude as it moves--simultaneously giving rise to the initial sound wave. In both experiment and computation, it is shown that the viscous effect plays a significant role in head-on BVIs.

Lee, Soogab↗
NASA NTRS1947

Wing plan forms for high-speed flight

It is pointed out that, in the case of an airfoil of infinite aspect ratio moving at an angle of sideslip, the pressure distribution is determined solely by that component of the motion in a direction normal to the leading edge. It follows that the attachment of plane waves to the airfoil at near-sonic or supersonic speeds (Ackeret theory) may be avoided and the pressure drag may be reduced by the use of plan forms in which the angle of sweepback is greater than the Mach angle. The analysis indicates that for aerodynamic efficiency, wings designed for flight at supersonic speeds should be swept back at an angle greater than the Mach angle, and the angle of sweepback should be such that the component of velocity normal to the leading edge is less than the critical speed of the airfoil sections. This principle may also be applied to wings designed for subsonic speeds near the speed of sound, for which the induced velocities resulting from the thickness might otherwise be sufficiently great to cause shock waves.

Jones, Robert J↗
DOE OSTI2020

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗