Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave propagation”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth-order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Incommensurate spin density wave and magnetocaloric effect in the metallic triangular lattice HoAl 2 Ge 2

Here, we report the magnetic structure and the magnetocaloric effect (MCE) of the ternary compound HoAl 2 Ge 2 with a trigonal CaAl 2 Si 2 -type crystal structure. A neutron powder diffraction experiment reveals that HoAl 2 Ge 2 exhibits an incommensurate spin density wave (SDW) with a propagation vector k=(0.23,0,0.06). The special arrangement of magnetic moments in HoAl 2 Ge 2 induces interesting physical phenomena and large magnetocaloric effects. The rise in resistivity at low temperatures indicates the effect of the SDW state in the electronic transport. The maximum magnetic-entropy change is –16.1J/kg K under a magnetic field change of 0–70 kOe for an isotropic HoAl 2 Ge 2 powder and it increases to –17.9J/kg K for a single crystal when the magnetic field (H) is applied parallel to the ab plane. A large rotating magnetic-entropy change of –5.1J/kg K for H=20 kOe in a HoAl 2 Ge 2 single crystal is obtained, which is closely associated to the magnetic anisotropy of the SDW order and its response to the external magnetic field. We discuss the large MCE in terms of the field-induced metamagnetic transition from the incommensurate SDW order to the ferromagnetic order. Our study establishes the triangular lattice R Al 2 Ge 2 (R=rare-earth elements) as a unique family of compounds to explore the existence of the incommensurate spin density waves and the correlated physical properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The Observed Effects of Cold Pools on Convection Triggering and Organization during DYNAMO/AMIE

Cold pools that form from existing convection can generate new clouds nearby, a process that has been suggested to contribute to cloud organization. However, some disagreement exists on the role of cold pools in cloud organization over tropical oceans, and much remains to be understood on how this role depends on the state of large-scale disturbances such as the Madden-Julian Oscillation (MJO). Here, this study addresses this question by examining the intraseasonal variability in the properties of cold pools and their effects on convection triggering using observations. The unique set of surface meteorology data and ground-based radars deployed during the DYNAMO/AMIE field campaign are used to identify cold pools and the associated spatiotemporal evolution of rainfall. The results suggest that cold pools enhance rainfall at their expected speed of 3–8 m s –1 , which includes the Doppler-shifting effect by the background wind. The cold pools also enhance rainfall mainly within a 20 km radius, especially in the direction of boundary-layer vertical wind shear. However, gravity waves appear to contribute predominantly to the propagation of convection over a wider range of distances due to their faster propagation speeds. The interface of convection triggering and propagation at different speeds of cold pools and gravity waves seems to contribute to cloud organization. The effectiveness of cold pools and gravity waves on cloud organization also depends highly on the MJO due to the combined effects of changes in environmental humidity and boundary to lower-tropospheric wind.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Non-equilibrium effective field theory and second sound

We investigate the phenomenon of second sound in various states of matter from the perspective of non-equilibrium effective field theory (EFT). In particular, for each state of matter considered, we find that at least two (though sometimes multiple) qualitatively different EFTs exist at finite temperature such that there is always at least one EFT with a propagating second-sound wave and at least one with no such second-sound wave. To aid in the construction of these EFTs, we use the method of cosets developed for non-equilibrium systems. It turns out that the difference between the EFTs with and without second-sound modes can be understood as arising from different choices of a new kind of inverse Higgs constraint. Finally, we demonstrate that it is possible to bypass the need for new inverse Higgs constraints by formulating EFTs on a new kind of manifold that is like the usual fluid worldvolume, but with reduced gauge symmetries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Role of α → ε → α phase transformation on the spall behavior of iron at atomic scales

Shock compression of iron microstructures above a threshold stress results in a α (BCC) → ε (HCP) transformation and the propagation of the release wave brings the metal back to the α phase following the ε→α transformation. Predicting failure behavior under shock loading conditions (spallation) relies on understanding the evolution of defects in the microstructure as it undergoes the α→ε→α phase transformation. This study uses molecular dynamics (MD) simulations to investigate the role of defect evolution during α→ε→α phase transformation on the spall strength values of single-crystal (sc) Fe microstructures. The MD simulations aim to characterize the ε phase fraction formed during shock compression and the defects during shock release for variations in loading orientations and shock stresses. The simulations are carried out for loading along the [100], [110], [111], and [112] orientations and for impact velocities ranging from 600 m/s to 1 km/s. The ε phase fractions during compression and defects (dislocations, twins) characterized during spall failure show an orientation dependence that affects the spall strength values. Here, the lowest value for spall strength is observed for the <110> loading orientation that shows a high density of twinning at the spall plane, whereas the highest value is observed for the <100> orientation and is associated with a α (BCC) → γ (FCC) transformation at the spall plane. The correlations of the spall strength values with the strain rates and with the ε phase fractions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The physics of desiccation cracks 1: Ductile fracturing and dependence on relative humidity

Ductile deformation is ubiquitously found in the shrinkage of geomaterials. The existence of ductility requires elastoplastic mechanics when analyzing the structure deformation under external stress. Here, in this work, we explore the physics of ductile fracturing based on the results from desiccation experiments and triaxial tests. By using the digital image correlation (DIC) method to generate the strain maps of samples undergoing desiccation cracking under different relative humidities, we obtain results showing the previously postulated Cnoidal Wave ductile failure patterns propagating under atmospheric condition-controlled crack velocities. We then correlate these observations with rate-dependent plasticity models calibrated through triaxial tests undergoing several unloading–reloading cycles and velocity stepping. This work demonstrates the necessity to consider ductility in soil cracking, indicating that the formation of crack patterns in soil desiccation is a slow and predictable process.

58 GEOSCIENCES↗
DOE OSTI2025

Quantum Simulation of Molecular Dynamics Processes─A Benchmark Study Using a Classical Simulator and Present-Day Quantum Hardware

Here, we explore how the fundamental problems in quantum molecular dynamics can be modeled using classical simulators (emulators) of quantum computers and the actual quantum hardware available to us today. The list of problems we tackle includes propagation of a free wave packet, vibration of a harmonic oscillator, and tunneling through a barrier. Each of these problems starts with the initial wave packet setup. Although Qiskit provides a general method for initializing wave functions, in most cases it generates deep quantum circuits. While these circuits perform well on noiseless simulators, they suffer from excessive noise on quantum hardware. To overcome this issue, we designed a shallower quantum circuit for preparing a Gaussian-like initial wave packet, which improves the performance of real hardware. Next, quantum circuits are implemented to apply the kinetic and potential energy operators for the evolution of a wave function over time. The results of our modeling on classical emulators of quantum hardware agree perfectly with the results obtained using the traditional (classical) methods. This serves as a benchmark and demonstrates that the quantum algorithms and Qiskit codes we developed are accurate. However, the results obtained on the actual quantum hardware available today, such as IBM’s superconducting qubits and IonQ’s trapped ions, indicate large discrepancies due to hardware limitations. This work highlights both the potential and challenges of using quantum computers to solve fundamental quantum molecular dynamics problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

An In-situ and Direct Confirmation of Super-Planckian Thermal Radiation Emitted From a Metallic Photonic-Crystal at Optical Wavelengths

Planck’s law predicts the distribution of radiation energy, color and intensity, emitted from a hot object at thermal equilibrium. The Law also sets the upper limit of radiation intensity, the blackbody limit. Recent experiments reveal that micro-structured tungsten can exhibit significant deviation from the blackbody spectrum. However, whether thermal radiation with weak non-equilibrium pumping can exceed the blackbody limit in the far field remains un-answered experimentally. Here, we compare thermal radiation from a micro-cavity/tungsten photonic crystal (W-PC) and a blackbody, which are both measured from the same sample and also in-situ. We show that thermal radiation can exceed the blackbody limit by >8 times at λ=1.7 μm resonant wavelength in the far-field. Our observation is consistent with a recent calculation by Wang and John performed for a 2D W-PC filament. This finding is attributed to non-equilibrium excitation of localized surface plasmon resonances coupled to nonlinear oscillators and the propagation of the electromagnetic waves through non-linear Bloch waves of the W-PC structure. This discovery could help create super-intense narrow band thermal light sources and even an infrared emitter with a laser-like input-output characteristic.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Giant enhancement of nonreciprocity in gyrotropic heterostructures

Abstract Nonreciprocity is a highly desirable feature in photonic media since it allows for control over the traveling electromagnetic waves, in a way that goes far beyond ordinary filtering. One of the most conventional ways to achieve nonreciprocity is via employing gyrotropic materials; however, their time-reversal-symmetry-breaking effects are very weak and, hence, large, bulky setups combined with very strong magnetic biases are required for technologically useful devices. In this work, artificial heterostructures are introduced to enhance the effective nonreciprocal behavior by reducing the contribution of the diagonal susceptibilities in the collective response; in this way, the off-diagonal ones, that are responsible for nonreciprocity, seem bigger. In particular, alternating gyrotropic and metallic or plasmonic films make an epsilon-near-zero (ENZ) effective-medium by averaging the diagonal permittivities of opposite sign, representing the consecutive layers. The homogenization process leaves unaltered the nonzero off-diagonal permittivities of the original gyrotropic substance, which become dominant and ignite strong nonreciprocal response. Realistic material examples that could be implemented experimentally in the mid-infrared spectrum are provided while the robustness of the enhanced nonreciprocity in the presence of actual media losses is discussed and bandwidth limitations due to the unavoidable frequency dispersion are elaborated. The proposed concept can be extensively utilized in designing optical devices that serve a wide range of applications from signal isolation and wave circulation to unidirectional propagation and asymmetric power amplification.

42 ENGINEERING↗
DOE OSTI2021

Understanding the phase transformation mechanisms that affect the dynamic response of Fe-based microstructures at the atomic scales

Large-scale molecular dynamics (MD) simulations were carried out to investigate the shock-induced evolution of microstructure in Fe-based systems comprising single-crystal and layered Cu/Fe alloys with a distribution of interfaces. The shock compression of pure single-crystal Fe oriented along [110] above a threshold pressure results in a BCC (α)→HCP (ε) phase transformation behavior that generates a distribution of ε phase variants in the phase transformed region of the microstructure behind the shock front. The propagation of the release wave through a phase transformed ε phase causes a reverse ε→α phase transformation and renders a distribution of twins for the [110] oriented Fe that serve as void nucleation sites during spall failure. The simulations reveal that the α→ε→α transformation-induced twinning for shock loading along the [110] direction is due to a dominant ε phase variant formed during compression that rotates on the arrival of the release wave followed by a reverse phase transformation to twins in the α phase. The modifications in the evolution of the ε phase variants and twins in Fe behavior are also studied for Cu–Fe layered microstructures due to the shock wave interactions with the Cu/Fe interfaces using a newly constructed Cu–Fe alloy potential. Here, the MD simulations suggest that interfaces affect the observed variants during shock compression and, hence, distributions of twins during shock release that affects the void nucleation stresses in the Fe phase of Cu/Fe microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A resistive MHD model and simulation on plasma flow evolution in the presence of resonant magnetic perturbation in a tokamak

Nonaxisymmetric magnetic fields such as the intrinsic error field and the externally applied resonant magnetic perturbation (RMP) in a tokamak are known to influence the plasma momentum transport and flow evolution through plasma response, which itself strongly depends on the plasma flow as well. The nonlinear interaction between plasma response and flow has been previously modeled in the conventional error field theory with the “no-slip” condition, which has been recently extended to allow the “free-slip” condition. In this work, we further target this specific process and numerically simulate the nonlinear plasma response and flow evolution in the presence of a single-helicity RMP in a circular-shaped model tokamak configuration, based on the full resistive MHD model in the initial-value code NIMROD. Time evolution of the parallel (to k) flow or “slip frequency” profile and its asymptotic steady state obtained from the NIMROD simulations are compared with both conventional and extended nonlinear response theories. Here, k is the wave vector of the propagating island. Good agreement with the extended theory with free-slip condition has been achieved for the parallel flow profile evolution in response to RMP in all resistive regimes, whereas the difference from the conventional theory with the no-slip condition tends to diminish as the plasma resistivity approaches zero.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Mesoscale simulation of granular materials under weak shock compaction–pore size distribution effects

This research established a systematic method to generate various pore-size distributions (PSDs) and studied the effect of PSDs on the shock compaction response of granular materials using two-dimensional mesoscale simulations under identical porosity. Simulations utilized various PSDs for three particle shapes (circle, ellipse, and square). Contacting particle configurations using three PSDs, characterized by spatially uniform distributed pores to heterogeneous distributed pores, and non-contacting particle configurations under a single case of PSD were tested. The PSD of generated particle sets was characterized using coordination number, mean diameter, and bimodality coefficient as statistical metrics. Mesoscale simulations showed that regardless of the conditions of pore distributions, shock compaction of granular materials consistently demonstrates a precursor, shock compaction front, and end. However, the shock compaction velocity of contacting particles was dependent on the PSDs despite the constant initial porosity. The compaction velocity was faster in particle configurations with relatively uniform pore distributions than in heterogeneous pore distributions, which our study demonstrated can be attributed to particle rearrangement during compaction. Circular-shaped particles had high sensitivity in shock compaction response to the various PSDs. Furthermore, a contacting particle configuration tended to propagate the shock compaction wave relatively faster than particles that were in a non-contact configuration. This study established the relative importance of considering PSD as a metric over the coordination number in studies of the shock compaction response of granular materials. Further, insights are provided on the evolving shock substructure to characterize the shock compaction response of granular materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Two types of superconducting pairs in stripe-ordered La 2 - x Ba x CuO 4 ( x = 1 / 8 ) : Evidence from resistivity measurements

Recent angle-resolved c-axis resistivity measurements of the stripe-ordered La 2-x Ba x CuO 4 (LBCO) with x = 1/8 revealed an unexpected dependence on the direction of the in-plane magnetic field. We argue that these and other available data for the c-axis transport point to the existence of superconducting pairs of two different types in the x = 1/8 LBCO below the stripe ordering temperature. The pairs of one type carry finite momentum and are confined to the Cu-O planes; the pairs of other type (probably the conventional d wave with zero momentum) propagate along narrow conducting channels traversing the sample in the c-axis direction. The evidence for this comes from the observed exponential temperature dependence of the c-axis resistivity ρ c (T ) which we attribute to the thermally excited slips of the superconducting phase and flux flows. We present a simple theory to fit the observed π/2-periodic dependence of ρ c on the direction of the in-plane magnetic field and the other data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Plasmon-enhanced optical nonlinearity in graphene nanomeshes

Using the density-matrix formalism, here we show that graphene nanomeshes (GNMs)—graphene sheets patterned with antidots—have large plasmon-enhanced nonlinear optical response. GNMs can be designed to behave as quasi-one-dimensional plasmonic crystals in which plasmons with large propagation lengths are efficiently excited. The associated third-order Kerr and third-harmonic-generation susceptibility can be as high as 10 –7 and 10 –9 m 2 V –2 , respectively, over the mid-to-near-infrared frequency range. Furthermore, carrier-density tuning in GNMs can flip the propagation direction of plasmonic waves and enables bidirectional switching of optical signals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Phonon-Driven Femtosecond Dynamics of Excitons in Crystalline Pentacene from First Principles

Nonradiative exciton relaxation processes are critical for energy transduction and transport in optoelectronic materials, but how these processes are connected to the underlying crystal structure and the associated electron, exciton, and phonon band structures, as well as the interactions of all these particles, is challenging to understand. Here, in this work, we present a first-principles study of exciton-phonon relaxation pathways in pentacene, a paradigmatic molecular crystal and optoelectronic semiconductor. We compute the momentum- and band-resolved exciton-phonon interactions, and use them to analyze key scattering channels. We find that both exciton intraband scattering and interband scattering to parity-forbidden dark states occur on the same ~100 fs timescale as a direct consequence of the longitudinal-transverse splitting of the bright exciton band. Consequently, exciton-phonon scattering exists as a dominant nonradiative relaxation channel in pentacene. We further show how the propagation of an exciton wave packet is connected with crystal anisotropy, which gives rise to the longitudinal-transverse exciton splitting and concomitant anisotropic exciton and phonon dispersions. Our results provide a framework for understanding the role of exciton-phonon interactions in exciton nonradiative lifetimes in molecular crystals and beyond.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗