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 235 records · Page 13

A Computational Study of the Thermodynamic Conditions Leading to Autoignition in Nanosecond Pulsed Discharges

Nanosecond pulsed discharges have attracted the attention of engine manufacturers due to the possibility of attaining distributed ignition sites that accelerate burn rates while resulting in very little electrode erosion. Multidimensional modeling tools currently capture the electrical structure of such discharges accurately, but resolving the chemical structure remains a challenging problem owing to the disparity of time-scales in streamer propagation (nanoseconds) and ignition phenomena (microseconds). The purpose of this study is to extend multidimensional results toward resolving the chemical structure in the wake of streamers (or the afterglow) by using a batch reactor model (BRM). This can afford the use of very detailed chemical kinetic information. The full nonequilibrium nature of the electrons is taken into account, along with fast gas heating, shock wave propagation, and thermal diffusion. The results shed light on ignition phenomena brought about by such discharges.

42 ENGINEERING↗

X-ray diffraction strains in laser-ablated aluminum, nickel, sodium, and invar: Pressures to 475 GPa

Dynamically compressed materials in longitudinal waves are described by two physical models: hydrostatic pressure, with equal, normal, principal stresses or material uniaxially strained in the wave propagation direction. These models are disparate, so experimental comparisons and evaluations are important. Polycrystalline material in a state of hydrostatic pressure will have no eccentricity of x-ray-diffracted Debye–Scherrer rings. A general three-dimensional solution of Bragg diffracted x rays based on principal crystallographic strains in the compression wave was found. The distortion of x-ray diffraction beams has been used for strain measurements; the analysis developed incorporates a strained reciprocal lattice and the incident x-ray beam. Strain-distorted Polanyi surfaces form an annulus of compression with an ellipsoid of revolution in reciprocal space, which is intersected by Ewald's sphere for Bragg diffraction. The in situ measurements for strain describe nanosecond diffraction evaluated using two planes, (hkl) and (h′k′l′), both in the same crystallographic phase. Diffraction from Al, Ni, Na, and Invar quantifies the compression axial strains in these materials: The compression axial ratios are 0.65, 1.05, 0.88, and 1.58 at pressures of 291, 402, 409, and 367 GPa for the respective materials. Crystal structure transformations with homogeneous pressurized stresses, mandating equal normal strains, should not be anticipated to agree with heterogeneous, uniaxially strained, and sheared crystalline phases. Measurements support in-plane strains increasing with pressure, p , in fcc and hcp aluminum as ε 22 = −0.21 p /TPa.

Alloys↗

Electromechanical response of piezoelectric thin films subjected to large strain

The use of piezoelectric thin sensors, when subjected to large static strain levels in ultrasonic applications is reported in this paper. A novel flipped beam bending method is proposed to induce static strain level nearing twice the incipient yield strain in the beam material (~0.39%). First, the effect of static large strain on ultrasonic guided wave propagation in the beam was studied using the velocity measurements obtained from a laser Doppler vibrometer. Next, the ultrasonic response of thin sensors is studied under compressive and tensile static strain (±0.72%) with the sensors symmetrically bonded on the lower and upper surfaces, and at the same axial location. The ultrasonic guided wave response to varying levels of static strain is recorded at different center frequencies (50–400 kHz) of the narrow banded ultrasonic excitation. As a case study, the variation of the ultrasonic guided wave response of polyvinylidene difluoride (PVDF) thin sensors with center frequencies at different static strain levels was analyzed and compared with that of lead zirconate titanate (PZT) thin sensors. PZT sensors show a higher variation of sensing performance with strain when compared to PVDF sensors. Apart from static strain, the possible effects on the change in performance of sensors, such as the change in the capacitance is discussed. The characterization method provides a reusable beam arrangement to study the ultrasonic sensor performance at large static strain levels. Characterization at large strain level allows the development of fail-safe sensors even under large strain and permanent plastic deformation of structure on which they are bonded, allowing accurate ultrasonic structural health monitoring even in the presence of damage progression. Piezopolymer sensor like PVDF are least susceptible to static strain due to large yield strain and low modulus and can serve as fail-safe sensors under large strain.

Rathod, Vivek↗

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗

Graphene-based analog of single-slit electron diffraction

In this study, we report on the experimental demonstration of single-slit diffraction exhibited by electrons propagating in encapsulated graphene with an effective de Broglie wavelength corresponding to their attributes as massless Dirac fermions. Nanometer-scale device designs were implemented to fabricate a single-slit followed by five detector paths. Predictive calculations were also utilized to readily understand the observations reported. These calculations required the modeling of wave propagation in ideal case scenarios of the reported device designs to more accurately describe the observed single-slit phenomenon. This experiment was performed at room temperature and 190 K, where data from the latter highlighted the exaggerated asymmetry between electrons and holes, recently ascribed to slightly different Fermi velocities near the K point. This observation and device concept may be used for building diffraction switches with versatile applicability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Coupling CTH to Linear Acoustic Propagation across an Air-Earth Interface

The interface between the Earth and the atmosphere forms a strong contrast in material properties. As such, numerical issues can arise when simulating an elastic wavefield across such a boundary when using a numerical simulation scheme. This is exacerbated when two different simulation codes are coupled straddling that interface. In this report we document how we implement the coupling of CTH, a nonlinear shock physics code, to a linearized elastic/acoustic wave propagation algorithm, axiElasti, across the air-earth interface. We first qualitatively verify that this stable coupling between the two algorithms produces expected results with no visible effects of the coupling interface. We then verify the coupling interface quantitatively by checking consistency with results from previous work and with coupled acoustic-elastic seismo-acoustic source inversions in three earth materials.

58 GEOSCIENCES↗

The Response of the Polar Vortex to Tropospheric Temperature Eddies in an Idealized General Circulation Model

Abstract A dry-core idealized general circulation model with a stratospheric polar vortex in the Northern Hemisphere is run with a combination of simplified topography and imposed tropospheric temperature perturbations, each located in the Northern Hemisphere with a zonal wavenumber of 1. The phase difference between the imposed temperature wave and the topography is varied to understand what effect this has on the occurrence of polar vortex displacements. Geometric moments are used to identify the centroid of the polar vortex for the purposes of classifying whether or not the polar vortex is displaced. Displacements of the polar vortex are a response to increased tropospheric wave activity. Compared to a model run with only topography, the likelihood of the polar vortex being displaced increases when the warm region is located west of the topography peak, and decreases when the cold region is west of the topography peak. This response from the polar vortex is due to the modulation of vertically propagating wave activity by the temperature forcing. When the southerly winds on the western side of the topographically forced anticyclone are collocated with warm- or cold-temperature forcing, the vertical wave activity flux in the troposphere becomes more positive or negative, respectively. This is in line with recent reanalysis studies that showed that anomalous warming west of the surface pressure high, in the climatological standing wave, precedes polar vortex disturbances.

54 ENVIRONMENTAL SCIENCES↗

Ultrasonic nondestructive evaluation of additively manufactured wear coatings

Here, this work pertains to the nondestructive evaluation of additively deposited coatings using ultrasound measurements. The specific objective was to evaluate using ultrasound surface wave measurements the quality of Stellite 21 coatings deposited on Inconel 718 substrates by the directed energy deposition (DED) additive manufacturing (AM) process. The surface wave speed and diffuse backscatter amplitude of the ultrasound waves were correlated with both cracking and warpage of the coating. This research is important, because, the integrity of DED-processed coatings currently requires a combination of destructive metallographic studies and X-ray computed tomography analysis, which are both expensive and time consuming. Instead of using a normal incidence configuration for the ultrasound measurements, the surface-wave approach allows the inspection of the near surface coating integrity. Three different frequencies were used to excite surface waves because each frequency has a different penetration depth. The signals obtained were used to quantify and compare the surface wave speed and surface wave-diffuse ultrasonic backscatter amplitude from five different DED-processed Stellite coated samples. Considering destructive metallurgical characterization as a reference, surface wave measurements were found to be effective for examining coating integrity. Further, the backscatter amplitude of the surface waves was correlated with the crack density and warpage. These flaws also change the stress state of the Stellite coating and consequently the surface wave propagation speed and its scattering behavior.

36 MATERIALS SCIENCE↗

Anisotropic magnon damping by zero-temperature quantum fluctuations in ferromagnetic CrGeTe 3

Spin and lattice are two fundamental degrees of freedom in a solid, and their fluctuations about the equilibrium values in a magnetic ordered crystalline lattice form quasiparticles termed magnons (spin waves) and phonons (lattice waves), respectively. In most materials with strong spin-lattice coupling (SLC), the interaction of spin and lattice induces energy gaps in the spin wave dispersion at the nominal intersections of magnon and phonon modes. Here we use neutron scattering to show that in the two-dimensional (2D) van der Waals honeycomb lattice ferromagnetic CrGeTe 3 , spin waves propagating within the 2D plane exhibit an anomalous dispersion, damping, and breakdown of quasiparticle conservation, while magnons along the c axis behave as expected for a local moment ferromagnet. These results indicate the presence of dynamical SLC arising from the zero-temperature quantum fluctuations in CrGeTe 3 , suggesting that the observed in-plane spin waves are mixed spin and lattice quasiparticles fundamentally different from pure magnons and phonons.

36 MATERIALS SCIENCE↗

Effect of Semiconducting Layer on Cable Insulation Damage Detection Using Surface Wave Reflectometry

This paper presents a study of the effect of the semiconducting layer present in medium and high voltage cables on surface wave propagation and its feasibility in detecting insulation damage on cables. Simulation studies were conducted using Ansys High Frequency Structure Simulator (HFSS). It was observed that attenuation increased with thicker semiconducting layer and higher conductivity of the layer. For instance, on a 24.4 mm diameter power cable, the electric field strength associated with a 1 GHz surface wave declines with distance and fades away after 30 m when a 2 mm thick semiconducting layer is present. More results as a function of semiconducting layer thickness, material conductivity and permittivity, distance of a damaged section, etc. will be presented.

reflectometry, unshielded power cable, semiconduct↗

Phason-Dominated Thermal Transport in Fresnoite

Fast-propagating waves in the phase of incommensurate structures, called phasons, have long been argued to enhance thermal transport. Although supersonic phason velocities have been observed, the lifetimes, from which mean free paths can be determined, have not been resolved. Using inelastic neutron scattering and thermal conductivity measurements, we establish that phasons in piezoelectric fresnoite make a major contribution to thermal conductivity by propagating with higher group velocities and longer mean free paths than phonons. In this work, the phason contribution to thermal conductivity is maximum near room temperature, where it is the single largest contributing degree of freedom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Large Propagation-Direction-Dependent Circularly Polarized Emission and Scattering Anisotropies of a Chiral Organic-Inorganic Semiconductor

Chiral materials are important tools for transducing circularly polarized light within many emerging opto-electronic and spin-based technologies. Here, we demonstrate that thin films of a bismuth iodide-based 0D chiral hybrid organic-inorganic semiconductor (HOIS) exhibit large anisotropy values in circularly polarized light emission (CPLE) that approach 50%, with mirror-image responses from front- and back-side measurements. A comprehensive analysis of light-wave propagation, absorption, emission, and scattering is constructed on the basis of a symmetry-derived exciton fine structure model, which accurately describes the direction and polarization dependence of the observed excitonic circular dichroism and CPLE, including contributions from both photoluminescence and resonant Raman scattering. Our analysis indicates that molecular chirality drives preferential film alignment with respect to the out-of-plane lattice vector direction, producing the observed anisotropies. This first demonstration in an HOIS system provides a unique route for enhancing polarization-dependent emission, and circularly polarized light transduction more broadly, in self-assembled HOIS.

14 SOLAR ENERGY↗

Connection between classical and quantum descriptions of spin waves using quantum circuits

A quantum computing circuit is presented that approximates a single spin wave quantum on a linear chain of spin 1/2 particles described by a Heisenberg Hamiltonian. The circuit is a product state where each qubit represents a spin. The spin wave motion is represented by opening the cone angle using Y rotations and then adding progressive Z rotations along the chain to represent wave propagation. We show analytically that this product state yields the correct dispersion relation in the limit of an unbounded chain. This observation is confirmed using both a simulator and various quantum processors. The use of the quantum computing paradigm in this case does not lead to a computational advantage, but rather leads to a novel conceptual connection between classical and quantum descriptions of spin waves, and may also be useful for characterizing the error in quantum processors.

magnons↗

Calculating RF current condensation with consistent ray-tracing and island heating

Finally, by exploiting the nonlinear amplification of the power deposition of RF waves, current condensation promises new pathways to the stabilization of magnetic islands. We present a numerical analysis of current condensation, coupling a geometrical optics treatment of wave propagation and damping to a thermal diffusion equation solver in the island. Taking into account the island geometry and relativistic damping, previous analytical theory can be made more precise and specific scenarios can be realistically predicted. With this more precise description, bifurcations and associated hysteresis effects could be obtained in an ITER-like scenario at realistic parameter values. Moreover, it is shown that dynamically varying the RF wave launching angles can lead to hysteresis and help to avoid the nonlinear shadowing effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Damage initiation and evolution in Al-Si layered microstructures under shock loading conditions at atomic scales

Designing materials with microstructural features like multiphase interfaces have shown significant promise for usage in the next generation defense and nuclear applications. The dynamic response of these interfaces in extreme environments is observed to vary with the deformability of individual phases, which could subsequently alter the favored damage nucleation sites related to spall failure. Here, we investigate the role of spacing of interfaces in a nanocrystalline layered Aluminum-Silicon system on the shock wave propagation behavior, and microstructural and defect evolution using classical molecular dynamics simulations. The simulations are carried out with different hypothetical Al/Si microstructures including variations in the distribution of Si grains as layers in a nanocrystalline Al matrix. The molecular dynamics simulations suggest that spall failure is preferentially initiated at Al/Si interfaces. In addition, for the same system volume and same concentration of Si, a higher number of layers is marked by an increase in shock wave velocity and reduced resistance to spall failure.

36 MATERIALS SCIENCE↗

Simulation of an inductively coupled plasma with a two-dimensional Darwin particle-in-cell code

A two-dimensional particle-in-cell code for the simulation of low-frequency electromagnetic processes in laboratory plasmas has been developed. The code uses the Darwin method omitting the electromagnetic wave propagation. The Darwin method separates the electric field into solenoidal and irrotational parts. The solenoidal electric field is calculated with a new algorithm based on the equation for the electric field vorticity. The system of linear equations in the new algorithm is readily solved using a standard iterative method. The irrotational electric field is the electrostatic field calculated with the direct implicit algorithm. The code is verified by reproducing the two-stream instability, electron electromagnetic waves, and shear Alfvén waves. The code is applied to simulate an inductively coupled plasma with the driving current flowing around the plasma region. In this simulation, a ring of dense plasma forms at the initial stage but then the density becomes maximal in the center and decays monotonically toward the walls. The skin effect is in the transitional mode between local and non-local, and the electron velocity distribution function is non-Maxwellian.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Collaborative Seismic Earth Model: Generation 2

Geological interpretations, earthquake source inversions and ground motion modeling, among other applications, require models that jointly resolve crustal and mantle structure. With the second generation of the Collaborative Seismic Earth Model (CSEM2), we present a global multi-resolution tomographic Earth model that serves this purpose. The model evolves through successive regional- and global-scale refinements. While the first generation aggregated regional models, with this study, we ensure consistency between all individual submodels, resulting in a model that accurately explains wave propagation across scales. Recent regional tomographic models were incorporated, comprising continental-scale inversions for Asia and Africa, as well as regional inversions for the Western US, Central Andes, Iran, and Southeast Asia. Across all regional refinements, over 793,000 source-receiver pairs contributed. Moreover, the long-wavelength Earth model (LOWE) introduces large-scale structures outside of pre-existing local refinements. A full-waveform inversion for global anisotropic P-and S-wave speed structure over a total of 194 iterations with a minimum period of 50 s on a large data set of 1 hr of waveform data from 2,423 earthquakes and over 6 million source-receiver pairs ensures that regional updates in the crust and uppermost mantle translate into updates of deeper, global-scale structure. To test the performance of CSEM2, we evaluate waveform fits between observed and synthetic seismograms at 50 s for an independent data set on the global scale, and on the regional scale for lower periods. We accurately simulate waveforms within and across regional refinements, maintaining the original resolution of the submodels embedded in the global framework.

58 GEOSCIENCES↗

A novel Doppler backscattering (DBS) system to simultaneously measure radio frequency plasma fluctuations and low frequency turbulence

A novel quadrature Doppler Backscattering (DBS) system has been developed and optimized for the E-band (60–90 GHz) frequency range using either O-mode or X-mode polarization in DIII-D plasmas. In general, DBS measures the amplitude of density fluctuations and their velocity in the lab frame. The system can simultaneously monitor both low-frequency turbulence (f < 10 MHz) and radiofrequency plasma density fluctuations over a selectable frequency range (20–500 MHz). Detection of high-frequency fluctuations has been demonstrated for low harmonics of the ion cyclotron frequency (e.g., 2f ci ~ 23 MHz) and externally driven high-frequency helicon waves (f = 476 MHz) using an adjustable frequency down conversion system. Importantly, this extends the application of DBS to a high-frequency spectral domain while maintaining important turbulence and flow measurement capabilities. This unique system has low phase noise, good temporal resolution (sub-millisecond), and excellent wavenumber coverage (k θ ~ 1–20 cm –1 and kr ≲ 30 cm –1 ). As a demonstration, localized internal DIII-D plasma measurements are presented from turbulence (f ≤ 5 MHz), Alfvenic waves (f ~ 6.5 MHz), ion cyclotron waves (f ≥ 20 MHz), as well as fluctuations around 476 MHz driven by an external high-power 476 MHz helicon wave antenna. In the future, helicon measurements will be used to validate GENRAY and AORSA modeling tools for prediction of helicon wave propagation, absorption, and current drive location for the newly installed helicon current drive system on DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗