SEARCH · Search NASA

Results for “wave speed”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Preconditioning for the Navier-Stokes equations with finite-rate chemistry

The extension of Van Leer's preconditioning procedure to generalized finite-rate chemistry is discussed. Application to viscous flow is begun with the proper preconditioning matrix for the one-dimensional Navier-Stokes equations. Eigenvalue stiffness is resolved and convergence-rate acceleration is demonstrated over the entire Mach-number range from nearly stagnant flow to hypersonic. Specific benefits are realized at the low and transonic flow speeds typical of complete propulsion-system simulations. The extended preconditioning matrix necessarily accounts for both thermal and chemical nonequilibrium. Numerical analysis reveals the possible theoretical improvements from using a preconditioner for all Mach number regimes. Numerical results confirm the expectations from the numerical analysis. Representative test cases include flows with previously troublesome embedded high-condition-number areas. Van Leer, Lee, and Roe recently developed an optimal, analytic preconditioning technique to reduce eigenvalue stiffness over the full Mach-number range. By multiplying the flux-balance residual with the preconditioning matrix, the acoustic wave speeds are scaled so that all waves propagate at the same rate, an essential property to eliminate inherent eigenvalue stiffness. This session discusses a synthesis of the thermochemical nonequilibrium flux-splitting developed by Grossman and Cinnella and the characteristic wave preconditioning of Van Leer into a powerful tool for implicitly solving two and three-dimensional flows with generalized finite-rate chemistry. For finite-rate chemistry, the state vector of unknowns is variable in length. Therefore, the preconditioning matrix extended to generalized finite-rate chemistry must accommodate a flexible system of moving waves. Fortunately, no new kind of wave appears in the system. The only existing waves are entropy and vorticity waves, which move with the fluid, and acoustic waves, which propagate in Mach number dependent directions. The nonequilibrium vibrational energies and species densities in the unknown state vector act strictly as convective waves. The essential concept for extending the preconditioning to generalized chemistry models is determining the differential variables which symmetrize the flux Jacobians. The extension is then straight-forward. This algorithm research effort will be released in a future version of the production level computational code coined the General Aerodynamic Simulation Program (GASP), developed by Walters, Slack, and McGrory.

Godfrey, Andrew G.↗

NASA NTRS1972

Sky-wave backscatter - A means for observing our environment at great distances.

During the last five years, much progress has been made in the understanding of sky-wave backscatter. An explanation of the various interacting phenomena is presented, as is a review of the current state of knowledge reflecting recent advances in observational methods and analytic techniques. New narrow-beam antennas, coupled with signal modulations that permit fine resolution in time delay, are beginning to yield information concerning the character of the scatterers, which now can be separately discerned. These narrow beams also permit study of polarization fading from small regions, and this shows promise as a means for learning the distant sea state. Doppler shifts of a fraction of a hertz on signals of tens of megahertz are separable, permitting isolation of sea returns from ground returns by virtue of the Doppler effect resulting from sea-wave speed; this also suggests a potential sea-monitoring principle. Despite these advances, there is little practical application of sky-wave backscatter as a means of environmental monitoring. This lack is attributed to the large remaining gaps in our understanding of the echoes and our inability to interpret the forms of data that can be acquired with equipment of reasonable cost.

Croft, T. A.↗

NASA NTRS1987

Computation of viscous blast wave solutions with an upwind finite volume method

A fully conservative, viscous, implicit, upwind, finite-volume scheme for the thin-layer Navier-Stokes equations is described with application to blast wave flow fields. In this scheme, shocks are captured without the oscillations typical of central differencing techniques and wave speeds are accurately predicted. The finite volume philosophy ensures conservation and since boundary conditions are also treated conservatively, accurate reflections of waves from surfaces are assured. Viscous terms in the governing equations are treated in a manner consistent with the finite volume philosophy, resulting in very accurate prediction of boundary layer quantities. Numerical results are presented for four viscous problems: a steady boundary layer, a shock-induced boundary layer, a blast wave/cylinder interaction and a blast wave/supersonic missile interaction. Comparisons of the results with an established boundary layer code, similarity solution, and experimental data show excellent agreement.

Molvik, Gregory A.↗

NASA NTRS1991

Computation of viscous blast wave flowfields

A method to determine unsteady solutions of the Navier-Stokes equations was developed and applied. The structural finite-volume, approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the interaction of blast-waves with stationary targets. The inviscid flux is evaluated using MacCormack's modified Steger-Warming flux or Roe flux difference splittings with total variation diminishing limiters, while the viscous flux is computed using central differences. The use of implicit boundary conditions in conjunction with a telescoping in time and space method permitted solutions to this strongly unsteady class of problems. Comparisons of numerical, analytical, and experimental results were made in two and three dimensions. These comparisons revealed accurate wave speed resolution with nonoscillatory discontinuity capturing. The purpose of this effort was to address the three-dimensional, viscous blast-wave problem. Test cases were undertaken to reveal these methods' weaknesses in three regimes: (1) viscous-dominated flow; (2) complex unsteady flow; and (3) three-dimensional flow. Comparisons of these computations to analytic and experimental results provided initial validation of the resultant code. Addition details on the numerical method and on the validation can be found in the appendix. Presently, the code is capable of single zone computations with selection of any permutation of solid wall or flow-through boundaries.

Atwood, Christopher A.↗

DOE OSTI2021

Individual Wave Detection and Tracking within a Rotating Detonation Engine through Computer Vision Object Detection applied to High-Speed Images

Known for their simplistic design and continuous detonation, rotating detonation engines (RDEs) constitute a majority of current pressure gain combustion (PGC) research efforts. Experimental RDE operation times have been continuously extended through the use of rig cooling techniques. As the window of observable behavior is expanded, and as the technology matures toward eventual integration within gas turbines, monitoring techniques must evolve to better match industrial diagnostics. High-speed image analysis techniques prove useful to capture and evaluate the unsteady detonation behavior within the RDE. Traditional image analysis techniques, however, require extensive processing times which prohibit simultaneous monitoring. To better address this problem, a computer vision object detection methodology is proposed to quickly detect individual detonation waves within a single down-axis image. Detonation waves are detected in individual images by the implemented computer vision method You Only Look Once (YOLO) object detection network. In order to detect detonation waves, the network must first be trained using RDE images of interest, for which each required phase of network development is outlined. Detection of waves is improved through proper treatment of the collected image set, variation of Intersection over Union (IoU) and confidence thresholding, and through a parametric study of annotation dimensions. Each detected wave is described by its location and rotational direction, and locations are tracked to calculate wave velocity across each frame, leading to a timestep resolution of 20 µs. Wave velocities are also calculated through a series of frames, leading to a suitable average velocity estimation using as few as 10 frames. Uncertainty analysis accounting for variation in camera framerate, pixel width and annotation centroid locations estimates a total uncertainty of ±4.3% for velocity calculations, using the smallest annotation boxes. This method offers great reductions in processing times, as a step toward real-time monitoring of detonation waves within an RDE. Improving on previous studies, this technique is impartial to wave modes not included in the original training set and calculates wave velocities independent of high-speed pressure data. The ability to isolate waves within predicted bounding boxes will likely facilitate analysis of pixel intensity variation as an estimation of wave strength in future work.

Johnson, Kristyn↗

DOE OSTI2022

Detonation Wave Profile (update of LA-UR-15-29498-rev1)

The Zel’dovich-von Neumann-Doering (ZND) profile of a detonation wave is derived. Two basic material properties are required: i. An equation of state (EOS) for a partly burned explosive; P(V, e, λ). ii. A burn rate for the reaction progress variable; d/dtλ = R(V, e, λ). For a steady planar detonation wave the reactive flow PDEs can be reduced to ODEs. The detonation wave profile can be determined from an ODE plus algebraic equations for points on the partly burned detonation loci with a specified wave speed. Furthermore, for the CJ detonation speed, the end of the reaction zone is sonic. A solution to the reactive flow equations can be constructed with a rarefaction wave following the detonation wave profile. This corresponds to an underdriven detonation wave, and the rarefaction is know as a Taylor wave.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

NASA NTRS1986

Analysis of the Giacobini-Zinner bow wave

The cometary bow wave of P/Giacobini-Zinner has been analyzed using the complete set of ICE field and particle observations to determine if it is a shock. Changes in the magnetic field and plasma flow velocities from upstream to downstream have been analyzed to determine the direction of the normal and the propagation velocity of the bow wave. The velocity has then been compared with the fast magnetosonic wave speed upstream to derive the Mach number and establish whether it is supersonic, i.e., a shock, or subsonic, i.e., a large amplitude wave. The various measurements have also been compared with values derived from a Rankine-Hugoniot analysis. The results indicate that, inbound, the bow wave is a shock with M = 1.5. Outbound, a subsonic Mach number is obtained, however, arguments are presented that the bow wave is also likely to be a shock at this location.

Smith, E. J.↗

NASA NTRS1986

Analysis of the Giacobini-Zinner bow wave

The cometary bow wave of P/Giacobini-Zinner was analyzed using the ICE field and particle observations to determine if it is a shock. Changes in the magnetic field and plasma flow velocities from upstream to downstream were analyzed to determine the direction of the normal and the propagation velocity of the bow wave. The velocity was compared with the fast magnetosonic wave speed upstream to derive the Mach number and establish whether it is supersonic, i.e., a shock, or subsonic, i.e., a large amplitude wave. The measurements were compared with values derived from a Rankine-Hugoniot analysis. The results indicate that inbound the bow wave is a shock with M = 1.5. Outbound, a subsonic Mach number is obtained. However, arguments that the bow wave is also likely to be a shock at this location are presented.

Smith, E. J.↗

NASA NTRS1989

A procedure for computing surface wave trajectories on an inhomogeneous surface

Equations are derived for computing surface waves on smooth surfaces, including surfaces with a nonuniform wave speed. The prior literature dealt primarily with the theoretical development with little consideration given to computational methods, and examples were limited to waves on surfaces of simple analytic description, such as cones, spheres, and cylinders. The computational procedure presented is a relatively general method. Sample calculations illustrate the procedure for a class of practical shapes of the type that include aerodynamic and hydrodynamic surfaces. Equations are also included for computing the spreading of rays into a surrounding medium that will support waves.

Barger, Raymond L.↗

NASA NTRS1978

Analytical study of laser-supported combustion waves in hydrogen

Laser supported combustion (LSC) waves are an important ingredient in the fluid mechanics of CW laser propulsion using a hydrogen propellant and 10.6 micron lasers. Therefore, a computer model has been constructed to solve the one-dimensional energy equation with constant pressure and area. Physical processes considered include convection, conduction, absorption of laser energy, radiation energy loss, and accurate properties of equilibrium hydrogen. Calculations for 1, 3, 10 and 30 atm were made for intensities of 10 to the 4th to 10 to the 6th W/sq cm, which gave temperature profiles, wave speed, etc. To pursue the propulsion application, a second computer model was developed to describe the acceleration of the gas emerging from the LSC wave into a variable-pressure, converging streamtube, still including all the above-mentioned physical processes. The results show very high temperatures in LSC waves which absorb all the laser energy, and high radiative losses.

Kemp, N. H.↗

NASA NTRS1992

Navier-Stokes simulations of unsteady transonic flow phenomena

Numerical simulations of two classes of unsteady flows are obtained via the Navier-Stokes equations: a blast-wave/target interaction problem class and a transonic cavity flow problem class. The method developed for the viscous blast-wave/target interaction problem assumes a laminar, perfect gas implemented in a structured finite-volume framework. The approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the full viscous terms are computed by central differencing. Comparisons of unsteady numerical, analytical, and experimental results are made in two- and three-dimensions for Couette flows, a starting shock-tunnel, and a shock-tube blockage study. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing of the predominantly inviscid flows. Viscous effects were increasingly significant at large post-interaction times. While the blast-wave/target interaction problem benefits from high-resolution methods applied to the Euler terms, the transonic cavity flow problem requires the use of an efficient scheme implemented in a geometrically flexible overset mesh environment. Hence, the Reynolds averaged Navier-Stokes equations implemented in a diagonal form are applied to the cavity flow class of problems. Comparisons between numerical and experimental results are made in two-dimensions for free shear layers and both rectangular and quieted cavities, and in three-dimensions for Stratospheric Observatory For Infrared Astronomy (SOFIA) geometries. The acoustic behavior of the rectangular and three-dimensional cavity flows compare well with experiment in terms of frequency, magnitude, and quieting trends. However, there is a more rapid decrease in computed acoustic energy with frequency than observed experimentally owing to numerical dissipation. In addition, optical phase distortion due to the time-varying density field is modelled using geometrical constructs. The computed optical distortion trends compare with the experimentally inferred result, but underpredicts the fluctuating phase difference magnitude.

Atwood, C. A.↗

DOE OSTI2021

Error estimates of finite difference methods for the Dirac equation in the massless and nonrelativistic regime

We present four frequently used finite difference methods and establish the error bounds for the discretization of the Dirac equation in the massless and nonrelativistic regime, involving a small dimensionless parameter 0 < ε &NestedLessLess; 1 inversely proportional to the speed of light. In the massless and nonrelativistic regime, the solution exhibits rapid motion in space and is highly oscillatory in time. Specifically, the wavelength of the propagating waves in time is at O(ε), while in space, it is at O(1) with the wave speed at O(ε -1 ). We adopt one leap-frog, two semi-implicit, and one conservative Crank-Nicolson finite difference methods to numerically discretize the Dirac equation in one dimension and establish rigorously the error estimates which depend explicitly on the time step τ, mesh size h, and the small parameter ε. The error bounds indicate that, to obtain the “correct” numerical solution in the massless and nonrelativistic regime, i.e., 0 < ε &NestedLessLess; 1, all these finite difference methods share the same ε-scalability as time step τ = O(ε 3/2 ) and mesh size h = O(ε 1/2 ). A large number of numerical results are reported to verify the error estimates.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1977

Nonlinear stress waves in a perfectly flexible string

This paper discusses nonlinear stress-wave propagation in a perfectly flexible string obeying a quasilinear (rate-dependent) constitutive equation. Wave speeds and compatibility relations valid along various families of characteristics were determined. It was shown that the compatibility relations associated with the transverse as well as the longitudinal waves readily yield a physical interpretation when they are expressed in suitable variables and in vector form. Coding based on the present information was completed for the machine solution of a class of mixed initial- and boundary-value problems of practical interest. Computer simulation of the stress-wave interaction in the 40-foot lanyard in the Arcas 'Rocoz' system during deployment was carried out using a stress-strain relation for nylon at the strain rate of 30/second. A method for estimating the maximum tension and strain in a string during the initial loading phase is proposed.

Fan, D.-N.↗

NASA NTRS1996

Nonstationary Gravity Wave Forcing of the Stratospheric Zonal Mean Wind

The role of gravity wave forcing in the zonal mean circulation of the stratosphere is discussed. Starting from some very simple assumptions about the momentum flux spectrum of nonstationary (non-zero phase speed) waves at forcing levels in the troposphere, a linear model is used to calculate wave propagation through climatological zonal mean winds at solstice seasons. As the wave amplitudes exceed their stable limits, a saturation criterion is imposed to account for nonlinear wave breakdown effects, and the resulting vertical gradient in the wave momentum flux is then used to estimate the mean flow forcing per unit mass. Evidence from global, assimilated data sets are used to constrain these forcing estimates. The results suggest the gravity-wave-driven force is accelerative (has the same sign as the mean wind) throughout most of the stratosphere above 20 km. The sense of the gravity wave forcing in the stratosphere is thus opposite to that in the mesosphere, where gravity wave drag is widely believed to play a principal role in decelerating the mesospheric jets. The forcing estimates are further compared to existing gravity wave parameterizations for the same climatological zonal mean conditions. Substantial disagreement is evident in the stratosphere, and we discuss the reasons for the disagreement. The results suggest limits on typical gravity wave amplitudes near source levels in the troposphere at solstice seasons. The gravity wave forcing in the stratosphere appears to have a substantial effect on lower stratospheric temperatures during southern hemisphere summer and thus may be relevant to climate.

Alexander, M. J.↗

NASA NTRS1992

Film condensation in a horizontal rectangular duct

Condensation heat transfer in an annular flow regime with and without interfacial waves was experimentally investigated. The study included measurements of heat transfer rate with condensation of vapor flowing inside a horizontal rectangular duct and experiments on the initiation of interfacial waves in condensation, and adiabatic air-liquid flow. An analytical model for the condensation was developed to predict condensate film thickness and heat transfer coefficients. Some conclusions drawn from the study are that the condensate film thickness was very thin (less than 0.6 mm). The average heat transfer coefficient increased with increasing the inlet vapor velocity. The local heat transfer coefficient decreased with the axial distance of the condensing surface, with the largest change at the leading edge of the test section. The interfacial shear stress, which consisted of the momentum shear stress and the adiabatic shear stress, appeared to have a significant effect on the heat transfer coefficients. In the experiment, the condensate flow along the condensing surface experienced a smooth flow, a two-dimensional wavy flow, and a three-dimensional wavy flow. In the condensation experiment, the local wave length decreased with the axial distance of the condensing surface and the average wave length decreased with increasing inlet vapor velocity, while the wave speed increased with increasing vapor velocity. The heat transfer measurements are reliable. And, the ultrasonic technique was effective for measuring the condensate film thickness when the surface was smooth or had waves of small amplitude.

Lu, Qing↗

NASA NTRS1991

Transport of gyration-dominated space plasmas of thermal origin. I - Generalized transport equations

Grad's 20-moment set of transport equations has been examined in the limit of strong external magnetic fields. This simplified set of equations describes the transport of mass and parallel momentum as well as the transport of parallel and perpendicular energy and heat flow in the magnetic field direction. The effect of collisions was calculated assuming a modified relaxation model. Wave speeds and normal modes of the simplified set of equations were examined for an ion and electron gas. It was found that four of the ten normal modes are electron thermal-heat waves which approximately decouple from the six ion waves in the system. When low-frequency waves are considered (slow-wave approximation), this allows the electron energy and heat flow equations to be solved separately from the ion equations and in a time-independent fashion. When this was done, it was found that, under certain conditions, these equations predict an infinite electron perpendicular temperature, Te-perpendicular, in the collisionless regime. This occurs whenever Te-perpendicular is greater than the parallel temperature, at any point along collisionless and diverging magnetic field lines. When applied appropriately, the significantly simplified set of equations derived here are well suited for application to a large variety of problems in planetary ionospheres and magnetospheres.

Gombosi, Tamas I.↗

DOE OSTI2021

Shock melting and the hcp-bcc phase boundary of Mg under dynamic loading

The high-temperature and high-pressure response of magnesium has been investigated through shock-release experiments performed up to shock melting. The longitudinal and bulk sound speeds of Mg are reported along the Hugoniot from 25 to 56 GPa and used to determine the elastic properties and Grüneisen parameter. The hexagonal close-packed (hcp)–body-centered cubic (bcc) phase transition is marked by a reduction in the determined shear wave speed. Thermal softening is observed to begin between 40 and 44 GPa, with incipient melt at 55.5 GPa, in close agreement with previous diffraction measurements under laser shock loading. Examination of the release profiles showed that two different responses were observed, depending on the peak stress. When plotted in a stress-energy phase diagram, the two responses are observed to form separate lines that intersect the static hcp-bcc phase transition and incipient melt. The results indicate that the hcp-bcc transition occurs on the Hugoniot at 28.4 GPa and place the hcp-bcc-liquid triple point ~20 GPa. In conclusion, a multiphase equation of state is developed which places the melt boundary below the previously reported static measurements.

36 MATERIALS SCIENCE↗

DOE OSTI2021

A Shallow-Water Model for Convective Self-Aggregation

Randomly distributed convective storms can self-aggregate in the absence of large-scale forcings. Here we present a 1D shallow-water model to study the convective self-aggregation. This model simulates the dynamics of the planetary boundary layer and represents convection as a triggered process. Once triggered, convection lasts for finite time and occupies finite length. Here we show that the model can successfully simulate self-aggregation, and that the results are robust to a wide range of parameter values. In the simulations, convection excites gravity waves. The gravity waves then form a standing wave pattern, separating the domain into convectively active and inactive regions. We analyze the available potential energy (APE) budget and show that convection generates APE, providing energy for self-aggregation. By performing dimensional analysis, we develop a scaling theory for the size of convective aggregation, which is set by the gravity wave speed, damping time scale, and number density of convective storms. This paper provides a simple modeling framework to further study convective self-aggregation.

54 ENVIRONMENTAL SCIENCES↗