Time-dependent density functional theory with the orthogonal projector augmented wave method
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
No abstract available
The travel times of the seismic waves obtained for the Apollo-14 and -16 active seismic experiments and the Apollo-16 grenade launches are shown to be consistent with a powder-layer model of the shallow lunar crust. The velocity variation with depth determined from these data is: V(z) = approximately 110 z to the 1/6 power m/sec for z less than 10 meters and V(z) is nearly = to 250 m/sec for z greater than 10 meters. The velocity values found for the 10 meter depth are similar to those found by Kovach, et al. (1972). The z to the 1/6 power depth dependence for the velocity of the topmost layer is that predicted on the basis of a powder layer (Gangi, 1972). The Amplitude variation of the direct waves as a function of source-to-receiver separation, x, is A(x) = A(o)x to the -n power exp(-ax) where 1.5 n 2.2 and a is nearly = to 0.047 neper/m. Velocity-spectra analyses of the direct, surface-reflected, bottom-reflected and refracted waves give results that are consistent with the velocity model inferred from the traveltime data.
Massive coalescing binary systems are candidate sources of gravitational radiation in the millihertz frequency band accessible to spacecraft Doppler tracking experiments. This paper discusses signal processing and detection probability for waves from coalescing binaries in the regime where the signal frequency increases linearly with time, i.e., 'chirp' signals. Using known noise statistics, thresholds with given false alarm probabilities are established for one- and two-spacecraft experiments. Given the threshold, the detection probability is calculated as a function of gravitational wave amplitude for both one- and two-spacecraft experiments, assuming random polarization states and under various assumptions about wave directions. This allows quantitative statements about the detection efficiency of these experiments and the utility of coincidence experiments. In particular, coincidence probabilities for two-spacecraft experiments are insensitive to the angle between the directions to the two spacecraft, indicating that near-optical experiments can be done without constraints on spacecraft trajectories.
Surface fluctuation spectroscopy is a dynamic light scattering technique which involves the detection of thermal fluctuations (capillary waves) that cover every liquid interface. The advantages of this technique over conventional mechanical methods are discussed. The signal resolution and precision are improved by using a grating heterodyne technique coupled with a digital correlator. Detailed calculations showed that, for the measured power spectrum (the Fourier transform of the correlation function) to take the Lorentzian shape, it is necessary that certain optical conditions be met. These conditions are demonstrated with experimental data collected for ethanol at the ambient temperature as a function of the wave number. Once the optical parameters are adjusted correctly, the experimental results show excellent agreement with the literature values. It is also demonstrated that, using microcomputer technology, the entire system can be automated.
Lower transition region models with a balance between mechanical heating and radiative losses are expanded to include wave pressure effects. The models are used to study the simple damping length form of the heating function. The results are compared to the results obtained by Woods et al. (1990) for solutions in the lower transition region. The results suggest that a mixture of fast-mode and slow-mode waves may provide the appropriate heating mechanism in the lower transition region, with the decline in effective vertical wave speed caused by the refraction and eventual total reflection of the fast-mode wave resulting from the decreasing atmospheric density.
One of the outstanding problems of collisionless shock physics is the microphysical description of electron thermalization and associated wave processes. Observations at the Earth's bow shock have shown various features in the electron velocity distribution such as a loss-cone, temperature anisotropies, and beams before a thermalized, flat-topped distribution is reached downstream. It is generally believed that in addition to the macroscopic fields, which determine the overall mapping of the distribution function, electromagnetic whistler waves and lower hybrid waves are responsible for the pitch-angle scattering and thermalization of the electrons. However, a genuine physical explanation requires analysis of the evolution of the electron distribution in the self-consistent wave fields. Particle simulations are ideally suited to do this, and include processes in which both electrons and ions are involved, such as cross-field streaming instabilities. We show results from l-D and 2-D simulations using an implicit code, which allows time steps large compared to the inverse electron plasma frequency. It is shown that in order to capture the correct propagation characteristics of the waves, a large mass ratio between the protons and electrons is generally necessary. Detailed diagnostics of the waves as well as the evolution of the electron velocity distribution will be presented and will be discussed in the frame work of satellite observations.
High-time-resolution magnetic-field data for different separations of ISEE 1 and 2 in the earth's ion foreshock region are examined to study the coherence length of upstream ULF waves. Examining the correlation coefficients of the low-frequency waves as a function of separation distance shows that the correlation coefficient depends mainly on the separation distance of ISEE 1 and 2 transverse to the solar-wind flow. It drops to about 0.5 when the transverse separation is about 1 earth radius, a distance much larger than the proton thermal gyroradius in the solar wind. Thus the coherence length of the low-frequency waves is about one earth radius, which is of the order of the wavelength, and is consistent with that estimated from the bandwidth of the waves.
An acoustic analogy is proposed that directly includes nonlinear propagation effects. We examine the Lighthill acoustic analogy and replace the Green's function of the wave equation with numerical solutions of the generalized Burgers' equation. This is justified mathematically by using similar arguments that are the basis of the solution of the Lighthill acoustic analogy. This approach is superior to alternatives because propagation is accounted for directly from the source to the far-field observer instead of from an arbitrary intermediate point. Validation of a numerical solver for the generalized Burgers' equation is performed by comparing solutions with the Blackstock bridging function and measurement data. Most importantly, the mathematical relationship between the Navier- Stokes equations, the acoustic analogy that describes the source, and canonical nonlinear propagation equations is shown. Example predictions are presented for nonlinear propagation of jet mixing noise at the sideline angle
An acoustic analogy is proposed that directly includes nonlinear propagation effects. We examine the Lighthill acoustic analogy and replace the Green's function of the wave equation with numerical solutions of the generalized Burgers' equation. This is justified mathematically by using similar arguments that are the basis of the solution of the Lighthill acoustic analogy. This approach is superior to alternatives because propagation is accounted for directly from the source to the far-field observer instead of from an arbitrary intermediate point. Validation of a numerical solver for the generalized Burgers' equation is performed by comparing solutions with the Blackstock bridging function and measurement data. Most importantly, the mathematical relationship between the Navier-Stokes equations, the acoustic analogy that describes the source, and canonical nonlinear propagation equations is shown. Example predictions are presented for nonlinear propagation of jet mixing noise at the sideline angle.
A simplified method called the transfer function technique (TFT) was devised for evaluating the stress wave environment in a structure containing internal equipment. The TFT consists of following the initial in-plane stress wave that propagates through a structure subjected to a dynamic load and characterizing how the wave is altered as it is transmitted through intersections of structural members. As a basis for evaluating the TFT, impact experiments and detailed stress wave analyses were performed for structures with two or three, or more members. Transfer functions that relate the wave transmitted through an intersection to the incident wave were deduced from the predicted wave response. By sequentially applying these transfer functions to a structure with several intersections, it was found that the environment produced by the initial stress wave propagating through the structure can be approximated well. The TFT can be used as a design tool or as an analytical tool to determine whether a more detailed wave analysis is warranted.
Mariner 10 measurements at the Mercury bow shock provide examples where the magnetic field is approximately parallel or perpendicular to the bow shock normal. Upstream of a broad irregular parallel shock, left hand circularly polarized waves are observed which cut off very sharply at approximately 4 Hz. Upstream of a perpendicular shock, right hand circularly polarized waves are observed which persist up to the Nyquist frequency of 12 Ha. Determination of the wave propagation vector as a function of frequency helps conclusively identify the waves as whistler mode waves propagating from the shock. The magnetosheath downstream of the parallel shock is disturbed more than that downstream of the perpendicular shock particularly below 1 Hz. In the latter case regular left hand polarized waves observed slightly above the proton gyrofrequency are identified as ion cyclotron waves with wavelength approximately 300 km which are Doppler shifted up to their observed frequency.
We report on the scaling and units of the Rayleigh wave elastic response function A R (ω) that is output from the widely used Computer Programs in Seismology (CPS) to aid in the modeling of ground motion sourced by atmospheric explosions. The program uses mixed units (km, second, km/s, gm/cc) to keep A R (ω) near 10 0 and prevent any numerical underflow or overflow. We compare two models for the response of an elastic half space to the output from CPS. Our application inputs the recommended, mixed unit geological models to determine how researchers must scale this output to obtain physical units for A R (ω) that represents the amplitude scaling for the minimum group velocity (Airy phase) contribution to Rayleigh waves. We determine that a CPS user must scale the output for A R (ω) by 10 -12 (m/km) 2 (g/cc/m 3 /kg) to obtain MKS (meter, kg, second) units and then must multiply this result by the vertical component eigenfunction squared, that is evaluated at the free surface.
This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.
This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. [1995] that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.
A full-wave propagation model was developed that describes the propagation of gravity waves from the Earth's surface to the upper boundary, which can be placed anywhere between 150 and 500 km altitude. The model includes a realistic background atmosphere, and includes the effects of mean horizontal winds and their vertical shears, mean vertical temperature gradients, the eddy and molecular diffusion of heat and momentum, and the effects of ion-drag. This model solves five coupled second-order differential equations (continuity, momentum, and energy) in the vertical coordinate to derive the perturbation variables u', v', w' (horizontal and vertical velocity components), T' (temperature) and p' (pressure). The upper boundary can be automatically selected based on tests using the radiation condition at the upper boundary, wherein the height is increased until the wave is experiencing severe dissipation at the upper boundary, ensuring that substantial absorption occurs for any waves reflected from the upper boundary. The determination of wave amplitude is a key requirement of wave energetics. Therefore, the fullwave model has been applied to airglow observations in order to determine wave amplitudes as a function of altitude. This was accomplished by using the full-wave model output to drive a chemistry perturbation module that describes minor species perturbations and the resulting airglow perturbations. The full-wave output was multiplied by an altitude-independent factor such that the modeled and observed relative airglow intensity perturbations were equal. The effects of mean winds were included in these studies, and found to be the most important model input affecting the calculations (being more important than the choice of eddy diffusion profiles and chemical kinetic coefficients). In one study (Hickey et al., 1997a) these winds could not be well estimated from the measurements, whereas in the second study (Hickey et al.,1997b) the mean were well defined with a sodium wind-temperature lidar.
This is essentially a progress report on a theoretical investigation of the propagation of transient waves in a random medium. The emphasis in this study is on applications to sonic-boom propagation, particularly as regards the effect of atmospheric turbulence on the sonic-boom waveform. The analysis is general, however, and is applicable to other types of waves besides sonic-boom waves. The phenomenon of primary concern in this investigation is the fine structure of the wave. A figure is used to illustrate what is meant by finestructure.
Using the Fokker-Planck equation, an examination is made of the acceleration of charged particles by fluctuational electric fields caused by the propagation of intense radio emission in the long-wave range. Correlation functions of electric fields near the boundary of plasma instability were computed. The interaction of charged particles with these fields is examined. It is shown that due to the interaction of particles with fluctuational fields, solar cosmic rays can be accelerated to energies of E=100 to 1000 MeV.