SEARCH · Engineering Papers
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.
Weakly nonlinear surface waves are examined in the context of the beam model for jetlike radio sources. By introducing a finite scale length, viz. the beam radius, geometrical dispersion can act to balance nonlinear wave growth and thereby produce solitons, localized wave packets of stable waveform. A method for obtaining a soliton equation from the MHD equations is presented and then applied to radio knots in jets.
Spheroidal-type seismogram study of wave generation, propagation, dispersion and attenuation in realistic spherical Earth models, especially Gutenberg-Bullen A-prime model
Available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, anelastic dispersion, sphericity, and gravity. The resulting models, for average Earth, average ocean, and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low-veocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studied. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggestive. For most of the Earth SHSV in the vicinity of the low-velocity zone. Near the East Pacific Rise, however, SVSH at depth, consistent with ascending flow. Anisotropy is as important as temperature in causing radial and lateral variations in velocity. The models have a high velocity nearly isotropic layer at the top of the mantle that thickens with age. This layer defines the LID, or seismic lithosphere. In the Pacific, the LID thickens with age to a maximum thickness of about 50 km. This thickness is comparable to the thickness of the elastic lithosphere. The LID thickness is thinner than derived using isotropic or pseudo-isotropic procedures A new model for Average Earth is obtained which includes a thin LID. This model extends the fit of a P.R.E.M. type model to shorter period surface waves.
Long period Rayleigh wave and Love wave dispersion data, particularly for oceanic areas, were not simultaneously satisfied by an isotropic structure. Available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, an elastic dispersion, sphericity, and gravity. The resulting models, for the average Earth, average ocean and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low velocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studies. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggested.
Long period Rayleigh wave and Love wave dispersion data, particularly for oceanic areas, were not simultaneously satisfied by an isotropic structure. Available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, an elastic dispersion, sphericity, and gravity. The resulting models, for the average earth, average ocean and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low velocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studies. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggested.
NASA grant NAG 5 2866, funded by the Astrophysics Theory Program, enabled the study the Kelvin-Helmholtz instability in protostellar jets. In collaboration with co-investigator Philip Hardee, the PI derived the analytic dispersion relation for the instability in including a cooling term in the energy equation which was modeled as one of two different power laws. Numerical solutions to this dispersion relation over a wide range of perturbation frequencies, and for a variety of parameter values characterizing the jet (such as Mach number, and density ratio) were found It was found that the growth rates and wavelengths associated with unstable roots of the dispersion relation in cooling jets are significantly different than those associated with adiabatic jets, which have been studied previously. In collaboration with graduate student Jianjun Xu (funded as a research associate under this grant), hydrodynamical simulations were used to follow the growth of the instability into the nonlinear regime. It was found that asymmetric surface waves lead to large amplitude, sinusoidal distortions of the jet, and ultimately to disruption Asymmetric body waves, on the other hand, result in the formation of shocks in the jet beam in the nonlinear regime. In cooling jets, these shocks lead to the formation of dense knots and filaments of gas within the jet. For sufficiently high perturbation frequencies, however, the jet cannot respond and it remains symmetric. Applying these results to observed systems, such as the Herbig-Haro jets HH34, HH111 and HH47 which have been observed with the Hubble Space Telescope, we predicted that some of the asymmetric structures observed in these systems could be attributed to the K-H modes, but that perturbations on timescales associated with the inner disk (about 1 year) would be too rapid to cause disruption. Moreover, it was discovered that weak shock 'spurs' in the ambient gas produced by ripples in the jet surface due to nonlinear, modes of surface and/or body waves could accelerate the ambient gas to low velocity. This latter effect represents a new mechanism by which supersonic jets can accelerate low velocity outflows.
SUMMARY The spectral element method is currently the method of choice for computing accurate synthetic seismic wavefields in realistic 3-D earth models at the global scale. However, it requires significantly more computational time, compared to normal mode-based approximate methods. Source stacking, whereby multiple earthquake sources are aligned on their origin time and simultaneously triggered, can reduce the computational costs by several orders of magnitude. We present the results of synthetic tests performed on a realistic radially anisotropic 3-D model, slightly modified from model SEMUCB-WM1 with three component synthetic waveform ‘data’ for a duration of 10 000 s, and filtered at periods longer than 60 s, for a set of 273 events and 515 stations. We consider two definitions of the misfit function, one based on the stacked records at individual stations and another based on station-pair cross-correlations of the stacked records. The inverse step is performed using a Gauss–Newton approach where the gradient and Hessian are computed using normal mode perturbation theory. We investigate the retrieval of radially anisotropic long wavelength structure in the upper mantle in the depth range 100–800 km, after fixing the crust and uppermost mantle structure constrained by fundamental mode Love and Rayleigh wave dispersion data. The results show good performance using both definitions of the misfit function, even in the presence of realistic noise, with degraded amplitudes of lateral variations in the anisotropic parameter ξ. Interestingly, we show that we can retrieve the long wavelength structure in the upper mantle, when considering one or the other of three portions of the cross-correlation time series, corresponding to where we expect the energy from surface wave overtone, fundamental mode or a mixture of the two to be dominant, respectively. We also considered the issue of missing data, by randomly removing a successively larger proportion of the available synthetic data. We replace the missing data by synthetics computed in the current 3-D model using normal mode perturbation theory. The inversion results degrade with the proportion of missing data, especially for ξ, and we find that a data availability of 45 per cent or more leads to acceptable results. We also present a strategy for grouping events and stations to minimize the number of missing data in each group. This leads to an increased number of computations but can be significantly more efficient than conventional single-event-at-a-time inversion. We apply the grouping strategy to a real picking scenario, and show promising resolution capability despite the use of fewer waveforms and uneven ray path distribution. Source stacking approach can be used to rapidly obtain a starting 3-D model for more conventional full-waveform inversion at higher resolution, and to investigate assumptions made in the inversion, such as trade-offs between isotropic, anisotropic or anelastic structure, different model parametrizations or how crustal structure is accounted for.
In this work, we derive from first principles the dispersion relation for electromagnetic radiation that is spontaneously emitted at a 2D curved plasma–vacuum interface in the absence of both background and external electric and magnetic fields. This plasma geometry is motivated by plasmas generated through hypervelocity impact projectiles—an area that is of importance to satellite safety in the aerospace community, and may be of importance for the developing study of projectile impact-driven inertial confinement fusion. We find spontaneous radiation occurs when thermal surface waves propagate along the edge of a circular plasma, with a frequency lower than the plasma frequency and spatial decay proportional to 1/r. This is distinct from more well-known mechanisms of radiation, in which the radiation frequency is greater than or equal to the plasma frequency. This theory is validated against 2D, electromagnetic particle-in-cell simulations and shows reasonable agreement. Finally, we discuss the effects of plasma geometry and attribute the radiation generation mechanism to the acceleration of charged surface waves as they traverse azimuthally along a curved path.
We demonstrate that a phonon stopband can be synthesized from an aperiodic structure comprising a discrete set of phononic filter stages. Each element of the set has a dispersion relation that defines a complete bandgap when calculated under a Bloch boundary condition. Hence, the effective stopband width in an aperiodic phononic filter (PnF) may readily exceed that of a phononic crystal with a single lattice constant or a coherence scale. With simulations of multi-moded phononic waveguides, we discuss the effects of finite geometry and mode-converting junctions on the phonon transmission in PnFs. The principles described may be utilized to form a wide stopband in acoustic and surface wave media. Relative to the quantum of thermal conductance for a uniform mesoscopic beam, a PnF with a stopband covering 1.6-10.4 GHz is estimated to reduce the thermal conductance by an order of magnitude at 75 mK.
In liquid rocket propulsion, the knowledge and the understanding of liquid-gas interfacial phenomena are very important. This is important for predicting the onset of cavitation occurring in swirl injection elements used in STME, as well as atomization processes in shear-induced injectors (co-axial) and impinging injector elements. From the fact that all the physical processes including droplet size distribution, droplet dispersion, mixing and combustion are controlled by atomization processes, it is expected that the successful incorporation of the volume of fraction (VOF) will greatly enhance the analytical capability of predicting spray combustion processes in liquid-fueled engines. In this paper, a methodology is developed to define and track interfaces between two fluids in non-orthogonal, body-fitted grids using a single fractional volume of fluid (VOF) variable to describe the distribution of the liquid phase in a gas-liquid flow field. This method was implemented in a mature CFD code MAST (Multiphase All-Speed Transient) utilizing the general PISO-C algorithm. For the preliminary study on the analysis of spray combustion and tracking of the interface between two phases, we will report on the progress of the simulation of the instability on the liquid column; the surface wave instability and the droplet breakup from the liquid surface.
Experimental data on the saturation and decay of a previously described induced scattering instability are presented. The instability is driven by static background field asymmetries and initially produces a single growing mode (ω 1 , k 1 ). As the instability develops, several harmonically related (ω n =nω 1 , k n =nk 1 ) eigenmodes of the plasma column grow up, resulting in a spatially localized moving potential perturbation with δϕ/T ≈ 20%. After a period of rapid growth, the mode amplitudes saturate and vary slowly for roughly 10 5 ω 1 −1 . This quasi-equilibrium ends with an unusual change in the eigenmode phases, wherein δϕ changes from positive to negative. Plasma heating during this period increases the damping of all modes and rapidly quenches the instability. This instability is of particular interest because the coupling with the field asymmetry produces an enhanced level of radial transport.
Hydromagnetic surface waves propagating on finitely conducting fluid using alternating magnetic field, including electromagnetic skin depth effects in computing theoretical wave dispersion
Both the theoretical and experimental investigations into current-wave interactions are discussed. The following three problems were studied: (1) the dispersive relation of a random gravity-capillary wave field; (2) the changes of the statistical properties of surface waves under the influence of currents; and (3) the interaction of capillary-gravity with the nonuniform currents. Wave current interaction was measured and the feasibility of using such measurements for remote sensing of surface currents was considered. A laser probe was developed to measure the surface statistics, and the possibility of using current-wave interaction as a means of current measurement was demonstrated.
The application of surface acoustic wave (SAW) signal processing devices in the space shuttle was explored. In order to demonstrate the functions which a SAW device might perform, a breadboard pulse compression filter (PCF) module was assembled. The PCF permits a pulse radar to operate with a large duty cycle and low peak power, a regime favorable to the use of solid state RF sources. The transducer design, strong coupling compensation, circuit model analysis, fabrication limitations, and performance evaluation of a PCF are described. The nominal value of the compression ratio is 100:1 with 10-MHz bandwidth centered at 60 MHz and 10-microsecond dispersive delay. The PCF incorporates dispersive interdigital transducers and a piezoelectric lithium niobate substrate.
Anisotropic waves in composites are considered, taking into account wave speeds, wave surfaces, flexural waves in orthotropic plates, surface waves, edge waves in plates, and waves in coupled composite plates. Aspects of dispersion in composites are discussed, giving attention to pulse propagation and dispersion, dispersion in rods and plates, dispersion in a layered composite, combined material and structural dispersion, continuum theories for composites, and variational methods for periodic composites. The characteristics of attenuation and scattering processes are examined and a description is given of shock waves and impact problems in composites. A number of experiments are also reported.
Differences in experimental results for surface gravity wave phase speeds obtained by cross spectra and slope-height spectra are compared. It is shown that, for a simple two-dimensional model containing both dispersive and non-dispersive waves, the observed differences can be explained readily. The slope-height technique measures a weighted average of the wavenumber squared; for frequencies sufficiently large compared with the dominant frequency, the computed phase-speed variation with frequency is representative of the wave having the larger wavenumber - the dispersive wave. For the cross-spectral method, it is shown that the small-probe-separation assumption usually employed is not valid for frequencies much larger than the dominant frequency, except at the singular point where both waves have exactly equal spectral densities and the phase function represents an average of the two modes. For all other cases, the phase function approaches that of the wave having the larger spectral density, and essentially ignores the presence of the other wave, even if both modes are relatively close in their contributions to the total spectral density.
The solar atmosphere and solar wind are magnetically structured. The structuring can include tangential discontinuities, which can support surface waves. Such waves can be dispersive. This means that dispersion and nonlinearity can balance in such a way that solitary waves (or solitons) can result. This general point is illustrated by a two-dimensional nonlinear analysis which explicitly demonstrates the presence of long-wavelength solitary waves propagating on tangential discontinuities. If the waves are only weakly nonlinear, then they obey the Korteweg-de Vries equation and are true solitons.
The results of a spatial stability analysis of a two-dimensional slab jet, in which optically thin radiative cooling is dynamically important, are presented. We study both magnetized and unmagnetized jets at external Mach numbers of 5 and 20. We model the cooling rate by using two different cooling curves: one appropriate to interstellar gas, and the other to photoionized gas of reduced metallicity. Thus, our results will be applicable to both protostellar (Herbig-Haro) jets and optical jets from active galactic nuclei. We present analytical solutions to the dispersion relations in useful limits and solve the dispersion relations numerically over a broad range of perturbation frequencies. We find that the growth rates and wavelengths of the unstable Kelvin-Helmholtz (K-H) modes are significantly different from the adiabatic limit, and that the form of the cooling function strongly affects the results. In particular, if the cooling curve is a steep function of temperature in the neighborhood of the equilibrium state, then the growth of K-H modes is reduced relative to the adiabatic jet. On the other hand, if the cooling curve is a shallow function of temperature, then the growth of K-H modes can be enhanced relative to the adiabatic jet by the increase in cooling relative to heating in overdense regions. Inclusion of a dynamically important magnetic field does not strongly modify the important differences between an adiabatic jet and a cooling jet, provided the jet is highly supermagnetosonic and not magnetic pressure-dominated. In the latter case, the unstable modes behave more like the transmagnetosonic magnetic pressure-dominated adiabatic limit. We also plot fluid displacement surfaces associated with the various waves in a cooling jet in order to predict the structures that might arise in the nonlinear regime. This analysis predicts that low-frequency surface waves and the lowest order body modes will be the most effective at producing observable features in the jet.