Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Wave functions”

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 1,711 records · Page 95

Comparison of interband related optical transitions and excitons in ZnGeN 2 and GaN

The optical dielectric function of ZnGeN 2 is calculated from the interband transitions using the energy bands calculated in the quasiparticle self-consistent (QS)GW^ method using two different levels of approximation: the independent particle approximation (IPA) and the Bethe-Salpeter Equation (BSE) approach. The first allows us to relate peaks in ε 2 to specific bands and k-points but does not include electron-hole interaction effects. The second includes electron-hole interaction or excitonic effects. The corresponding changes in the shape of ε 2 (ω) are found to be similar to those in GaN. The screened Coulomb interaction W^ is here calculated including electron-hole interactions in the polarization function and gives a band structure already going beyond the random phase approximation. The static dielectric constants including only electronic screening, commonly referred to as ε∞ , were calculated separately by extrapolating the wave vector dependent macroscopic dielectric function in the long-wavelength limit. Below the quasiparticle gap, we find three bound excitons optically active for different polarization. The convergence of these bound excitons with respect to the density of the k-mesh used in the BSE is studied and found to require a fine mesh. It is also found that these bound excitons originate from only the lowest conduction band and the top three valence bands. To incorporate the lattice screening, we include a scaling factor (ε∞/ε 0 ) 2 , which allows us to obtain exciton binding energies of the correct order of magnitude similar to those in GaN. Here, the excitons are related to each of the three fold split valence bands and the splittings of the latter are also studied as function of strain. Finally, a relation between the anisotropic effective masses and the valence band splitting is pointed out and explained.

36 MATERIALS SCIENCE↗

Three-fluid solar wind model with Alfven waves

We present a study of a three-fluid solar wind model. with continuity, momentum and separate energy equations for protons. alpha particles and electrons. Allowing separate coronal heat sources for all three species, we study the flow properties of the solar wind as a function of heat input, Alfven wave energy input, and alpha particle abundance.

Esser, Ruth↗

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. So, we describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

97 MATHEMATICS AND COMPUTING↗

Observation of Chorus Waves by the Van Allen Probes: Dependence on Solar Wind Parameters and Scale Size

Highly energetic electrons in the Earths Van Allen radiation belts can cause serious damage to spacecraft electronic systems and affect the atmospheric composition if they precipitate into the upper atmosphere. Whistler mode chorus waves have attracted significant attention in recent decades for their crucial role in the acceleration and loss of energetic electrons that ultimately change the dynamics of the radiation belts. The distribution of these waves in the inner magnetosphere is commonly presented as a function of geomagnetic activity. However, geomagnetic indices are nonspecific parameters that are compiled from imperfectly covered ground based measurements. The present study uses wave data from the two Van Allen Probes to present the distribution of lower band chorus waves not only as functions of single geomagnetic index and solar wind parameters but also as functions of combined parameters. Also the current study takes advantage of the unique equatorial orbit of the Van Allen Probes to estimate the average scale size of chorus wave packets, during close separations between the two spacecraft, as a function of radial distance, magnetic latitude, and geomagnetic activity, respectively. Results show that the average scale size of chorus wave packets is approximately 13002300 km. The results also show that the inclusion of combined parameters can provide better representation of the chorus wave distributions in the inner magnetosphere and therefore can further improve our knowledge of the acceleration and loss of radiation belt electrons.

Aryan, Homayon↗

Revisiting Competing Paths in Electrochemical CO 2 Reduction on Copper via Embedded Correlated Wavefunction Theory

We re-evaluate two key steps in the mechanism of CO 2 reduction on copper at a higher level of theory capable of correcting inherent errors in density functional theory (DFT) approximations, namely, embedded correlated wavefunction (ECW) theory. Here, we consider the CO reduction step on Cu(111), which is critical to understanding reaction selectivity. Here, we optimize embedding potentials at the periodic plane-wave DFT level using density functional embedding theory (DFET). All possible adsorption sites (adsites) for each adsorbate then are screened with ECW theory at the catalytically active site to refine the local electronic structure. Unsurprisingly, DFT and ECW theory predict different adsite preferences, largely because of DFT’s inability to properly situate the CO 2π* level. Differing preferred adsites suggest that different reaction pathways could emerge from DFT versus ECW theory. Starting from these preferred ECW theory adsites, we then obtain reaction pathways at the plane-wave DFT level using the climbing-image nudged elastic band method to determine minimum energy paths. Thereafter, we perform ECW calculations at the catalytically active site to correct the energetics at each interpolated structure (image) along the reaction pathways. Via this approach, we confirm that the first step in CO reduction via hydrogen transfer on Cu(111) is to form hydroxymethylidyne (*COH) instead of formyl (*CHO). Although the prediction to preferentially form *COH is consistent with that of DFT, the two theories predict quite different structural and mechanistic behaviors, suggesting that verification is needed for other parts of the mechanism of CO 2 reduction, which is the subject of ongoing work.

30 DIRECT ENERGY CONVERSION↗

Resonance frequency shift of an acoustic chamber containing a rigid sphere

The resonance frequency shift of an acoustic rectangular chamber due to the presence of a rigid sphere has been measured for l = 1,2 modes as a function of sphere size and position. The frequency shift is the results of volume exclusion and wave scattering. An analytical Green's function calculation was used to explain the data, providing excellent agreement between the measured and the calculated values. Also reported are similar measurements for a thin disk and the ratio of second harmonic to fundamental pressure as a function of sphere position. The measurement shows that the sphere reduces the first harmonic content, with sharply peaked suppression minima at specific sphere positions.

Leung, E.↗

Radiation Belt Electron Dynamics Driven by Large-Amplitude Whistlers

Acceleration of radiation belt electrons driven by oblique large-amplitude whistler waves is studied. We show analytically and numerically that this is a stochastic process; the intensity of which depends on the wave power modified by Bessel functions. The type of this dependence is determined by the character of the nonlinear interaction due to coupling between action and phase. The results show that physically significant quantities have a relatively weak dependence on the wave power.

whistlers↗

Multi-dimensional MHD simple waves

In this paper we consider a formalism for multi-dimensional simple MHD waves using ideas developed by Boillat. For simple wave solutions one assumes that all the physical variables (the density rho, gas pressure p, fluid velocity V, gas entropy S, and magnetic induction B in the MHD case) depend on a single phase function phi(r,t). The simple wave solution ansatz and the MHD equations then require that the phase function has the form phi = r x n(phi) - lambda(phi)t, where = n(phi) = Delta phi / (absolute value of Delta phi) is the wave normal and lambda(phi) = omega/k = -phi t / (absolute value of Delta phi) is the normal speed of the wave front. The formalism allows for more general simple waves than that usually dealt with in which n(phi) is a constant unit vector that does not vary along the wave front. The formalism has implications for shock formation for multi-dimensional waves.

Webb, G. M.↗

Magnetic structure of the low beta, quasi-perpendicular shock

The structure of the low-beta quasi-perpendicular shock is examined in view of ISEE 1 and 2 magnetic field measurements. An analysis of shock overshoots indicates that the strength of the overshoots of low-beta, quasi-perpendicular shocks increases as the ratio of the Mach number to the first critical Mach number increases. Wave analysis indicates that the power of the downstream waves also increases as a function of this ratio of criticality. The thickness of the shock is a factor of 1-2 times greater than a precursor wavelength, countering the conjecture that the shock is the last amplified cycle of the precursor wave.

Farris, M. H.↗

Surface Instability of Liquid Propellant under Vertical Oscillatory Forcing

Fluid motion in a fuel tank produced during thrust oscillations can circulate sub-cooled hydrogen near the liquid-vapor interface resulting in increased condensation and ullage pressure collapse. The first objective of this study is to validate the capabilities of a Computational Fluid Dynamics (CFD) tool, CFD-ACE+, in modeling the fundamental interface transition physics occurring at the propellant surface. The second objective is to use the tool to assess the effects of thrust oscillations on surface dynamics. Our technical approach is to first verify the CFD code against known theoretical solutions, and then validate against existing experiments for small scale tanks and a range of transition regimes. A 2D axisymmetric, multi-phase model of gases, liquids, and solids is used to verify that CFD-ACE+ is capable of modeling fluid-structure interaction and system resonance in a typical thrust oscillation environment. Then, the 3D mode is studied with an assumed oscillatory body force to simulate the thrust oscillating effect. The study showed that CFD modeling can capture all of the transition physics from solid body motion to standing surface wave and to droplet ejection from liquid-gas interface. Unlike the analytical solutions established during the 1960 s, CFD modeling is not limited to the small amplitude regime. It can extend solutions to the nonlinear regime to determine the amplitude of surface waves after the onset of instability. The present simulation also demonstrated consistent trends from numerical experiments through variation of physical properties from low viscous fluid to high viscous fluids, and through variation of geometry and input forcing functions. A comparison of surface wave patterns under various forcing frequencies and amplitudes showed good agreement with experimental observations. It is concluded that thrust oscillations can cause droplet formation at the interface, which results in increased surface area and enhanced heat transfer between the liquid and gas phases as the ejected droplets travel well into the warmer gas region.

Yang, H. Q.↗

$S$ Hmax orientation in the Alpine region from observations of stress-induced anisotropy of nonlinear elasticity

The orientation of $S$ Hmax is commonly estimated from in situ borehole breakouts and earthquake focal mechanisms. Borehole measurements are expensive, and therefore sparse, and earthquake measurements can only be made in regions with many well-characterized earthquakes. Here, we derive the stress-field orientation using stress-induced anisotropy in nonlinear elasticity. In this method, we measure the strain derivative of velocity as a function of azimuth. We use a natural pump-probe (NPP) approach which consists of measuring elastic wave speed using empirical Green’s functions (probe) at different points of the earth tidal strain cycle (pump). The approach is validated using a larger data set in the Northern Alpine Foreland region where the orientation of maximum horizontal compressive stress is known from borehole breakouts and drilling-induced fractures. The technique resolves NNW-SSW to N-S directed $S$ Hmax which is in good agreement with conventional methods and the recent crustal stress model. We confirm that the NPP method can be applied to dense large-scale seismic arrays. The technique is then applied to the Southern Alps to understand the contemporary stress pattern associated with the ongoing deformation due to counterclockwise rotation of the Adriatic plate with respect to the European plate. Our results explain why the two major faults in Northeastern Italy, the Giudicarie Fault and the Periadriatic Line (Pustertal–Gailtal Fault) are currently inactive, while the currently acting stress field allows faults in Slovenia to deform actively. We have demonstrated that the pump-probe method has the potential to fill in the measurement gap left by conventional approaches, both in terms of regional coverage and in depth.

58 GEOSCIENCES↗

Comment on 'Generation of broadband noise in the magnetotail by the beam acoustic instability' by P.B. Dusenbery

It is shown here that of the three plasma models for the generation of broadband noise in the magnetotail proposed by Dusenbery and Lyons (1985), two result in the same types of instabilities that are also excited in the Grabbe and Eastman (GE, 1984) model, while the third model introduces other modes not present in the GE model. For the plasma parameters given in the GE model, the ion/ion acoustic instability is a nonresonant one, even though the phase velocity of the excited waves falls within the distribution function of the core ions. The nonlinear saturation mechanism of the ion/ion acoustic instability is the trapping of both the core and the beam ions. In a reply, Dusenbery addresses several ongoing controversies which have resulted from studies of wave particle interactions in the plasma sheet boundary layer.

Omidi, N.↗

On the generation of magnetohydrodynamic waves in a stratified and magnetized fluid. II - Magnetohydrodynamic energy fluxes for late-type stars

Magnetohydrodynamic (MHD) wave energy fluxes for late-type stars are calculated, using previously obtained formulae for the source functions for the generation of MHD waves in a stratified, but otherwise uniform, turbulent atmosphere; the magnetic fields in the wave generation region are assumed to be homogeneous. In contradiction to previous results, it is shown that in this uniform magnetic field case there is no significant increase in the efficiency of MHD wave generation, at least within the theory's limits of applicability. The major results are that the MHD energy fluxes calculated for late-type stars are less than those obtained for compressible modes in the magnetic field-free case, and that these MHD energy fluxes do not vary enough for a given spectral type to explain the observed range of UV and X-ray fluxes from such stars. It is therefore concluded that MHD waves in stellar atmospheres with homogeneous magnetic fields in the wave generation region cannot explain the observed stellar coronal emissions; if such MHD waves are responsible for a significant component of stellar coronal heating, then nonuniform fields within the generation region must be appealed to.

Musielak, Z. E.↗

​​Hydrogen Transportation and Distributed Energy Systems Seismic Risk Assessment for Cascadia Subduction Zone Airport Facilities​

Portland International Airport in Oregon is exploring operating a fleet of 28 fuel cell electric buses to support airport operations and possibly provide backup power during outages. In this report, we evaluate the risk of hydrogen deployment at an airport considering the potential for seismic activity. We present a seismic risk assessment for pressurized piping for a generic hydrogen refueling station at the airport. Only pressurized piping was considered in the risk analysis because of its vulnerability and availability of fragility information unlike the other components. Seismic capacity of other components should be included as more information becomes available. To characterize the seismic hazard, we considered the initiating event frequency of a Cascadia Subduction Zone earthquake event based on data from the literature. Seismic stressors that translate to the site from an earthquake event are expressed as peak ground acceleration. Site disturbance is characterized as a function of soil conditions, represented by different shear wave velocity values, which affect how seismic waves propagate and impact structures. Pipe fragility curves as a function of seismic capacity were used to correlate ground acceleration and failure probability. These fragility characteristics were combined with site-specific soil conditions to calculate the probability of pipe rupture during a seismic event. It should be noted that multiple simultaneous failures due to an earthquake event as a common cause are not considered. Some epistemic uncertainties such as self-exciting shocks, aftershocks, and the time variant nature of the ground motion are not considered either.

08 HYDROGEN↗

From linear mechanics to nonlinear mechanics

Consideration is given to the techniques used in telecommunication where a nonlinear system (the modulator) results in a linear transposition of a signal. It is then shown that a similar method permits linearization of electromechanical devices or nonlinear mechanical devices. A sweep function plays the same role as the carrier wave in radio-electricity. The linearizations of certain nonlinear functionals are presented.

RESEARCH TECHNIQUE, MATHEMATICS↗

The Effects of Wave Escape on Fast Magnetosonic Wave Turbulence in Solar Flares

One of the leading models for electron acceleration in solar flares is stochastic acceleration by weakly turbulent fast magnetosonic waves ("fast waves"). In this model, large-scale flows triggered by magnetic reconnection excite large-wavelength fast waves, and fast-wave energy then cascades from large wavelengths to small wavelengths. Electron acceleration by large-wavelength fast-waves is weak, and so the model relies on the small-wavelength waves produced by the turbulent cascade. In order for the model to work, the energy cascade time for large-wavelength fast waves must be shorter than the time required for the waves to propagate out of the solar-flare acceleration region. To investigate the effects of wave escape, we solve the wave kinetic equation for fast waves in weak turbulence theory, supplemented with a homogeneous wave-loss term.We find that the amplitude of large-wavelength fast waves must exceed a minimum threshold in order for a significant fraction of the wave energy to cascade to small wavelengths before the waves leave the acceleration region.We evaluate this threshold as a function of the dominant wavelength of the fast waves that are initially excited by reconnection outflows.

Pongkitiwanichakul, Peera↗