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.
Numerical solutions of a damped, nonlinear wave equation are presented. The equation describes the propagation of waves in a narrow thermocline or inversion which lose energy by exciting internal waves in the weakly stratified ambient environment. The results provide estimates for the persistence of finite-amplitude internal waves propagating in a thermoclinic waveguide.
The nonlinear parabolic equation describing the propagation of the electromagnetic wave in a semiconductor with the superlattice is analyzed. The possibility of the existence of the solitary waves is proved both for a small amplitude of the electrical field and the latter moderate values.
Laboratory experiments on important plasma physics issues of electrodynamic tethers were performed. These included current propagation, formation of wave wings, limits of current collection, nonlinear effects and instabilities, charging phenomena, and characteristics of transmission lines in plasmas. The experiments were conducted in a large afterglow plasma. The current system was established with a small electron-emitting hot cathode tethered to an electron-collecting anode, both movable across the magnetic field and energized by potential difference up to V approx.=100 T(sub e). The total current density in space and time was obtained from complete measurements of the perturbed magnetic field. The fast spacecraft motion was reproduced in the laboratory by moving the tethered electrodes in small increments, applying delayed current pulses, and reconstructing the net field by a linear superposition of locally emitted wavelets. With this technique, the small-amplitude dc current pattern is shown to form whistler wings at each electrode instead of the generally accepted Alfven wings. For the beam electrode, the whistler wing separates from the field-aligned beam which carries no net current. Large amplitude return currents to a stationary anode generate current-driven microinstabilities, parallel electric fields, ion depletions, current disruptions and time-varying electrode charging. At appropriately high potentials and neutral densities, excess neutrals are ionized near the anode. The anode sheath emits high-frequency electron transit-time oscillations at the sheath-plasma resonance. The beam generates Langmuir turbulence, ion sound turbulence, electron heating, space charge fields, and Hall currents. An insulated, perfectly conducting transmission line embedded in the plasma becomes lossy due to excitation of whistler waves and magnetic field diffusion effects. The implications of the laboratory observations on electrodynamic tethers in space are discussed.
Recent advances in computational fluid mechanics are discussed in reviews and reports. Sections are devoted to (1) the modeling of local phenomena and edge effects in solids, (2) stochastic modeling and simulation of fracture toughness, and (3) partitioning strategy and new finite elements. Particular attention is given to global and local finite-element/spectral-boundary-element techniques for failure analysis, simulations of microfracture in metal-matrix composites, fatigue analysis of cracked anisotropic plates under stochastic loading, mathematical modeling for the analysis of nonlinear aircraft dynamics, physical and mathematical modeling of wave propagation in the Ariane 5 VEB structure, partitioning based on symmetry transformations, an FEM approach to adaptive reliability assurance, and time-domain FEMs for the large rotational dynamics of multibody systems.
The Internal Propagation problems, Fan Noise problem, and Turbomachinery Noise problems are solved using the space-time conservation element and solution element (CE/SE) method. The problems in internal propagation problems address the propagation of sound waves through a nozzle. Both the nonlinear and linear quasi 1D Euler equations are solved. Numerical solutions are presented and compared with the analytical solution. The fan noise problem concerns the effect of the sweep angle on the acoustic field generated by the interaction of a convected gust with a cascade of 3D flat plates. A parallel version of the 3D CE/SE Euler solver is developed and employed to obtain numerical solutions for a family of swept flat plates. Numerical solutions for sweep angles of 0, 5, 10, and 15 deg are presented. The turbomachinery problems describe the interaction of a 2D vortical gust with a cascade of flat-plate airfoils with/without a downstream moving grid. The 2D nonlinear Euler Equations are solved and the converged numerical solutions are presented and compared with the corresponding analytical solution. All the comparisons demonstrate that the CE/SE method is capable of solving aeroacoustic problems with/without shock waves in a simple and efficient manner. Furthermore, the simple non-reflecting boundary condition used in the CE/SE method which is not based on the characteristic theory works very well in 1D, 2D and 3D problems.
Linear and nonlinear propagating characteristics of right-hand polarized, slow electromagnetic, magnetoplasma waves (whistlers) are discussed in terms of stability and dispersion. An analysis of the stability of whistlers propagating at an angle to the static magnetic field is presented. A new mechanism is derived for the onset of stimulated emissions, and modulational instability for nonlinear whistlers are discussed.
Recent thermal plasma data and a computer code by Horne (1989) are employed to evaluate the linear-path-integrated gain of waves propagating through the Io to predict the Jovian plasma-wave environment. Estimates of the nonlinear saturation amplitudes are utilized with the thermal plasma data from two frequency bands to study the convective growth of the ion-cyclotron (IC) waves. Strong cyclotron resonant damping is theorized to prevent wave propagation to the lower latitudes, and the thermal plasma and cyclotron resonant energetic ions are expected to further confine the IC waves. L-mode waves below the O(+) gyrofrequency in the equatorial region of the torus are shown to inhabit an unstable region. The IC waves probably achieve nonlinear amplitudes regardless of plasma properties due to the rapid amplification in this region. It is suggested that the Ulysses data complicate the identification of the waves because the magnetometer is not adequately sensitive and because of the low frequency of the plasma-wave detector.
Topological mechanical metamaterials have enabled new ways to control stress and deformation propagation. Exemplified by Maxwell lattices, they have been studied extensively using a linearized formalism. Herein, we study a two-dimensional topological Maxwell lattice by exploring its large deformation quasi-static response using geometric numerical simulations and experiments. We observe spatial nonlinear wave-like phenomena such as harmonic generation, localized domain switching, amplification-enhanced frequency conversion, and solitary waves. We further map our linearized, homogenized system to a non-Hermitian, nonreciprocal, one-dimensional wave equation, revealing an equivalence between the deformation fields of two-dimensional topological Maxwell lattices and nonlinear dynamical phenomena in one-dimensional active systems. Our study opens a regime for topological mechanical metamaterials and expands their application potential in areas including adaptive and smart materials and mechanical logic, wherein concepts from nonlinear dynamics may be used to create intricate, tailored spatial deformation and stress fields greatly transcending conventional elasticity.
A simple theory is presented which determines details of the breaking and saturation of a gravity wave as it propagates upward in the atmosphere. Breaking and saturation are here due to nonlinear advection analogous to the breaching of a surface wave and to the breaking of a planetary wave. Much simplification is obtained by the assumption that in a wave packet consisting of a primary wave and its harmonics, the primary wave remains dominant. This assumption, referred to a quasi-monochromatic approximation, is suggested by observations. Determined by this approximate theory are: a detailed picture of the waveform as it steepens and breaks; harmonics of the wave; the turbulence generation; and an underlying relationship between superadiabatic lapse rate and saturation by wave-wave interactions.
The effect of convective melt motion on the growth of morphological instabilities in crystal growth has been the focus of many studies in the past decade. While most of the efforts have been directed towards investigating the linear stability aspects, relatively little attention has been devoted to experimental and numerical studies. In a pure morphological case, when there is no flow, morphological changes in the solid-liquid interface are governed by heat conduction and solute distribution. Under the influence of a convective motion, both heat and solute are redistributed, thereby affecting the intrinsic morphological phenomenon. The overall effect of the convective motion could be either stabilizing or destabilizing. Recent investigations have predicted stabilization by a flow parallel to the interface. In the case of non-parallel flows, e.g., stagnation point flow, Brattkus and Davis have found a new flow-induced morphological instability that occurs at long wavelengths and also consists of waves propagating against the flow. Other studies have addressed the nonlinear aspects (Konstantinos and Brown, Wollkind and Segel)). In contrast to the earlier studies, our present investigation focuses on the effects of the potential flow fields typically encountered in Hele-Shaw cells. Such a Hele-Shaw cell can simulate a gravity-free environment in the sense that buoyancy-driven convection is largely suppressed, and hence negligible. Our interest lies both in analyzing the linear stability of the solidification process in the presence of potential flow fields, as well as in performing high-accuracy nonlinear simulations. Linear stability analysis can be performed for the flow configuration mentioned above. It is observed that a parallel potential flow is stabilizing and gives rise to waves traveling downstream. We have built a highly accurate numerical scheme which is validated at small amplitudes by comparing with the analytically predicted results for the pure morphological case. We have been able to observe nonlinear effects at larger times. Preliminary results for the case when flow is imposed also provide good validation at small amplitudes.
A simplified nonlinear numerical model for the development of incompressible magnetohydrodynamics in the presence of a strong magnetic field B|| and stratification, nicknamed 'Shell-Atm,' is presented. In planes orthogonal to the mean field, the nonlinear incompressible dynamics is replaced by two-dimensional shell models for the complex variables u and b, allowing one to reach large Reynolds numbers while at the same time carrying out sufficiently long integrations to obtain good statistics at moderate computational cost. The shell models of different planes are coupled by Alfve'n waves propagating along B||. The model may be applied to open or closed magnetic field configurations where the axial field dominates and the plasma pressure is low; here we apply it to the specific case of a magnetic loop of the solar corona heated by means of turbulence driven by photospheric motions, and we use statistics for its analysis. The Alfven waves interact nonlinearly and form turbulent spectra in the directions perpendicular and, through propagation, also parallel to the mean field. A heating function is obtained and shown to be intermittent; the average heating is consistent with values required for sustaining a hot corona and is proportional to the aspect ratio of the loop to the -1.5 power, and haracteristic properties of heating events are distributed as power laws. Crosscorrelations show a delay of dissipation compared with energy content.
Data from the LIMS instrument for January 1979 are used to provide further evidence for the often observed vacillation between the amplitudes of waves 1 and 2 in the stratosphere. The vacillation is shown to result primarily from nonlinear wave-wave interactions within the stratosphere. Two ways of interpreting nonlinearity are discussed. In the first, the basic state is defined to include large amplitude waves as well as the mean zonal wind. A forced wave propagates with respect to this asymmetric basic state, which can lead to changes in the conventional zonal wavenumber measured at one latitude. The other view of nonlinearity, interaction of wave with the zonal flow and with other wavenumbers are considered separately. Wave-wave interactions among waves 1, 2 and 3 are calculated. The derivation and computation of wave-wave interaction terms in the potential enstrophy balance are given. The observations indicate that enstrophy transfer among waves can be substantial even when the amplitude of one of the contributing waves is small. The computed enstrophy balance also demonstrates that wave-wave interactions can have a large effect on the interaction of waves with the mean flow.
SPE is a project in order to develop new, more physics-based, seismic models of explosions (see Nelson et al. 2013). One key component of this effort is numerical modeling enabled by modern state-of-art code/software. Accurate modeling of the shape and amplitude of seismic waves from their generation to their propagation to remote monitoring seismic stations is important to our ability to determine the origin and strength of the source from remote recording. In this project, the modeling is performed by coupling two codes HOSS and SPECFEM3D. HOSS models the dynamic nonlinear processes happening near the explosion. SPECFEM3D computes the propagation of seismic waves as they travel through 3D complex Earth models. SPECFEM3D was modified in order to be driven by a set of time-series calculated by HOSS in lieu of a seismic source. The goal of this summer project is the investigation of several questions pertinent to the establishment of a full end-to-end modeling capability from the high-rate strain regime area to remote distances where seismic station record seismic waves generated by explosions. The investigated questions are: (1) How to perform proper filtering? Direct-solution modeling only sustains a limited range of frequency depending on the grid size. As we go from one modeling domain to the other via coupling, the mesh size is getting coarser to allow modeling at large scale but also to account for the fact that high-frequency waves do not physically travel to large distances. So filtering of the time-series generated by the near-field hydrodynamic modeling is a current practice often employed but its effects on the modeled waveforms has to be investigated. (2) Quantitative assessment of the efficiency of attenuation to remove high-frequency content of the wavefield. This assessment will allow to create meshes with a grid size appropriate to the actual physics of wave propagation for a given explosion. (3) Checking that the coupling process and the two codes respect the cylindrical symmetry that is expected in the case of a pure explosion in a half-space. (4) The effect of the state-of-stress in the near-source area on the modeled seismic waveforms. These questions will be investigated through the modeling of SPE-4P, the fourth explosion of this series because it was designed to have little interaction with the geologic surrounding and the free surface of the Earth so that it is the most explosion-like experiment, with the most symmetries to be verified.
For electromagnetic wave propagation and scattering in an anisotropic random medium, the Dyson equation for the mean field and the Bethe-Salpeter equation for the correlation or the covariance of the field were derived. With the random permittivity expressed in a general anisotropic form, the bilocal and the nonlinear approximations are employed to solve the Dyson equation, and the ladder approximation to solve the Bethe-Salpeter equation. The mean dyadic Green's function for a two-layer anisotropic random medium with arbitrary three-dimensional correlation functions has been investigated with the zeroth-order solutions to the Dyson equation under the nonlinear approximation. The effective propagation constants are calculated for the four characteristic waves associated with the coherent vector fields, propagating in an anisotropic random-medium layer, which are the ordinary and extraordinary waves with upward- and downward-propagating vectors.
Whistlers nonlinear interaction during cold magnetoplasma propagation, considering energy exchanges between waves
Five relativistic plane nonlinear waves were investigated: circularly polarized waves and electrostatic plasma oscillations propagating parallel to the magnetic field, relativistic Alfven waves, linearly polarized transverse waves propagating in zero magnetic field, and the relativistic analog of the extraordinary mode propagating at an arbitrary angle to the magnetic field. When the ions are driven relativistic, they behave like electrons, and the assumption of an 'electron-positron' plasma leads to equations which have the form of a one-dimensional potential well. The solutions indicate that a large-amplitude superluminous wave determines the average plasma properties.
A test particle approach is used to compare gyroresonant pitch angle scattering of energetic electrons by coherent versus incoherent whistler mode waves, for the case in which the coherent wave amplitude is below the nonlinear phase trapping threshold. Wave packets of 400 ms duration propagating along the magnetic field at L = 4 within the plasmasphere are considered, and the wave-induced pitch angle scattering along the propagation path from one hemisphere to the other and the resulting precipitation flux are computed. An incoherent wave spectrum is simulated by random modulation of the wave frequency at intervals of 1 ms, thereby generating signals with nearly constant power spectral density over a bandwidth of 2 kHz centered at 5.5 kHz. The associated pitch angle scattering is compared with that of a monochromatic 5.5-kHz signal of 400 ms duration. Results of the test particle analysis are compared with those expected on the basis of a classical diffusion treatment, and an expression is derived for an effective “diffusion” coefficient for pitch angle scattering by coherent waves. The trajectory followed by a particle when interacting with incoherent waves essentially represents a random walk in velocity space, while for coherent waves the pitch angle of the particle varies in a well-defined manner. In spite of the fact that individual particle scatterings are typically larger for coherent waves, the peak precipitation fluxes induced by incoherent waves are found to be approximately the same as those for coherent waves having the same total power. This results from the fact that incoherent waves interact with particles over a wider range of energies. As a consequence, the energy spectrum and the temporal extent of transient precipitation pulses due to incoherent wave packets are broader than those for equivalent coherent ones.
Recent observations of binary, pre-main-sequence, solar-type stars provide evidence that such systems may coexist with circumstellar disks. The binary disk systems, besides being of general interest for the study of star formation, potentially provide useful tests of companion-disk interaction theories prominent in current hypotheses of planet formation. In this paper, we apply an asymptotic analysis of the nonlinear, resonant interaction of a stellar companion with a disk to understand the dependence of such interactions on the properties of the system: the binary mass ratio, the physical properties of the disk, and the effective dissipation (treated herein as viscosity). The method is based on a WKBJ approximation and exploits the conditions that the disk is thin and much less massive than the primary, but does not require that the companion-induced disturbance be small. Both isothermal and adiabatic responses are treated. Only circular orbit resonances are considered in this paper. It is demonstrated that the temperature of the disk as well as the relative mass of the companion affects the degree of nonlinearity, and that nonlinearity promotes high wave compression ratios, long wavelengths, and increased propagation distances. Nevertheless, the total torque exerted between the companion and the disk is well represented by linear theory. The amplitudes of density disturbances are reduced by viscosity and nonisothermality. Because resonant interactions are generally strong and capable of driving rapid evolution, one might expect observations of systems undergoing strong, resonant-driven evolution to be rare. In this connection, it is pointed out that the m = 1 resonance is distinguished by being anomalously weaker than the others and is therefore of observational interest. It is speculated that, in conditions of intrinsically small dissipation, the propagation of resonant-driven density waves is limited by the tendency of their wavelength to diminish with distance, and that the propagation distance (and therefore the region of the disk to which angular momentum is redistributed) is set by the distance at which the wavelength becomes comparable to the disk thickness.