Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated waves”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Magnetoelastic coupling and spin contributions to entropy and thermal transport in biferroic yttrium orthochromite

Direct engineering of material properties through exploitation of spin, phonon, and charge-coupled degrees of freedom is an active area of development in materials science. However, the relative contribution of the competing orders to controlling the desired behavior is challenging to decipher. In particular, the independent role of phonons, magnons, and electrons, quasiparticle coupling, and relative contributions to the phase transition free energy largely remain unexplored, especially for magnetic phase transitions. Here, we study the lattice and magnetic dynamics of biferroic yttrium orthochromite using Raman, infrared, and inelastic neutron spectroscopy techniques, supporting our experimental results with first-principles lattice dynamics and spin-wave simulations across the antiferromagnetic transition at T N ~ 138 K. Spectroscopy data and simulations together with the heat capacity (C p ) measurements, allow us to quantify individual entropic contributions from phonons (0.01 ± 0.01k B atom -1 ), dilational (0.03 ± 0.01k B atom -1 ), and magnons (0.11 ± 0.01k B atom -1 ) across TN. High-resolution phonon measurements conducted in a magnetic field show that anomalous T-dependence of phonon energies across TN originates from magnetoelastic coupling. Additionally, phonon scattering is primarily governed by the phonon–phonon coupling, with little contribution from magnon–phonon coupling, short-range spin correlations, or magnetostriction effects; a conclusion further supported by our thermal conductivity measurements conducted up to 14 T, and phenomenological modeling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

High-field side scrape-off layer density profile measurements and implications for high-field side LHCD coupling in DIII-D

The high-field side (HFS) scrape-off layer (SOL) is an often under-diagnosed region of tokamak plasmas. Situated in a region with favorable curvature, the HFS SOL has minimal turbulence-induced radial transport compared to the low-field side. Using profile reflectometry, the HFS SOL density profile is measured with high temporal resolution for a wide range of DIII-D plasma discharges for the first time. The magnetic configuration, particularly the location of the secondary separatrix (SS) largely determines the HFS SOL density profile. Additionally, density perturbations induced by edge-localized modes (ELMs) are observed in the HFS SOL, and the impact of ELMs on the SOL density profile is characterized for a range of magnetic configurations. It is found that ELM-induced HFS SOL perturbations are generally localized to within the SS and ELM-induced changes to the HFS wall density are minimized in near double-null configurations. The characterization of the HFS SOL density profile on DIII-D is crucial for the effective coupling of the planned HFS lower hybrid current drive (LHCD) launcher. Experimental HFS SOL measurements are used as input to full-wave simulations of LHCD coupling. Furthermore, this work simulates the effect of the magnetic configuration on LHCD coupling. It also simulates LHCD coupling in ELMs and high q min discharges. During such discharges, the LHCD coupling is predicted to be resilient to ELM-induced SOL density perturbations.

high-field side↗
DOE OSTI2023

Broadening of microwave heating beams in the DIII-D tokamak by edge turbulence

Abstract We have demonstrated for the first time that turbulent plasma density fluctuations in the edge of the DIII-D tokamak are responsible for substantial broadening of an injected microwave beam by successful quantitative comparison between experimental observations and first principles 2D full-wave simulations. The broadening of the beam has important implications for control of tokamak discharges through localized electron cyclotron deposition needed for eliminating magnetohydrodynamic instabilities. This new predictive capability is mandatory to design & operate present & future tokamaks in such a way that microwave heating schemes achieve their intended objectives.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Conditional diffusion machine-learning framework for mapping valence electron distribution from convergent beam electron diffraction

Quantitative convergent beam electron diffraction (CBED) enables determination of aspherical valence electron distributions through refinement of low-order structure factors, which are highly sensitive to chemical bonding and charge density variations. However, conventional quantitative CBED (QCBED) requires solving a highly nonlinear inverse problem with many coupled parameters, and computationally intensive dynamical diffraction calculations, making it time-consuming and difficult to apply to complex systems. More broadly, reconstructing charge density and orbital electron distribution from diffraction data has long been a central challenge in both x-ray and electron crystallography. Here, in this study, we introduce an artificial-intelligence (AI)-based framework that replaces traditional refinement with a data-driven inverse solver. Using a large synthetic CBED dataset generated by Bloch-wave simulations, we train a conditional diffusion model to directly infer crystal structural parameters and multipole density formalism parameters, and hence valence electron distributions, from CBED patterns alone. By learning from forward simulations across realistic parameter space, the model effectively solves the inverse problem. Compared with direct regression approaches, the diffusion-based framework provides posterior parameter distributions for rigorous uncertainty quantification while preserving quantitative fidelity and reducing analysis time by orders of magnitude. By eliminating the need for external single-crystal x-ray diffraction data and complex nonlinear refinement, this approach enables practical, high-throughput, and in situ quantitative CBED, enabling real-time mapping of valence electron distributions and their correlation with functional responses in quantum and energy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

A Survey on the Expanding Scope and Interdisciplinary Opportunities for Processing-in-Memory Techniques

Processing-in-Memory (PIM) is emerging as a practical path to overcome the limitations of traditional von Neumann architectures. At its core, PIM systems implement computing primitives such as logic operations and multiply-accumulate acceleration through compute-in-memory, near-memory processing, or hybrid designs. The role of memory cells varies widely across technologies, acting as inputs, outputs, or analog accumulators through bit-lines and sense amplifiers. This diversity creates trade-offs in precision, bandwidth, latency, and programmability, making it difficult to build a unified understanding on the progress of the field. In this survey, we organize recent advances of PIM into three areas. First, we discuss the progress on the architectural optimizations of PIM and its integration with both DRAM and emerging non-volatile memories. Second, we examine how PIM is being used to accelerate key computing domains, including generative AI workloads and high-performance kernels, along with new approaches. Third, we highlight the growing adoption of PIM in computational sciences, where it is being applied to solve interdisciplinary problems such as genome analysis, mRNA quantification, mass spectrometry, quantum circuit simulation, wave modeling, and secure computation. Finally, we synthesize the major challenges that continue to slow PIM adoption, including manufacturing constraints, power delivery, thermal reliability, data consistency, runtime and memory-management coordination, and the difficulty of building portable software abstractions without sacrificing commercial viability. This work provides an updated, structured perspective on PIM’s potential across computing and computational sciences and the barriers that must be solved for it to reach its full impact.

Asifuzzaman, Kazi [Oak Ridge National Laboratory (↗
DOE OSTI2021

A Novel Miniaturized MIMO Microstrip Patch Antenna Array with EBG- and Slot-Loading

A compact four-element multiple-input multiple-output (MIMO) antenna array based on the miniaturized mi-crostrip patch is presented. The reduction in array form factor is realized by partially populating the ground plane with elec-tromagnetic band gap (EBG) cells constructed from cross slots, in combination with reactively loading each patch with a slot. The resultant 2 × 2 array system demonstrates excellent inter-element isolation and can be easily fabricated using traditional milling and printing techniques, as it exploits a 2D layout that can be conveniently implemented on two sides of a single-layer substrate, without the need for any via structures. The operation of the array near the 2.5 GHz band is verified with full-wave simulations. The size of the new design is approximately 42% of that employing conventional patch elements.

Liao, DaHan↗
DOE OSTI2024

Conductive Spray Paint Fabricated Cylindrical Coplanar Waveguide for Noncontact Liquid-Level Sensing

A noncontact liquid-level sensor using a cylindrical coplanar waveguide (CCPW) is proposed. The concept is demonstrated using a prototype device made by directly spraying conductive paint on the outer surface of a glass cylinder. The position of the air–liquid interface is measured using the time-domain reflection (TDR) signal from the CCPW. The geometric parameters of the prototype sensor were optimized using a full-wave simulation tool. The experimental results show that the prototype CCPW-based liquid-level sensor has a resolution of 29.8μm over a range of 210 mm at the sweeping frequency bandwidth of 500 MHz–4 GHz. The tradeoff between the resolution and the measurement range was discussed. The proposed sensor provides a low-cost, easy-to-fabricate, and noncontact method for liquid-level measurement, which has potential applications in fuel storage, chemical production, and wastewater treatment.

42 ENGINEERING↗
DOE OSTI2021

Penetration Through Slots in Overmoded Cavities

A resonant cavity undergoes three distinct behaviors with increasing frequency: 1) fundamental modes, localized in frequency with well defined modal distribution; 2) undermoded region, where modes are still separated, but are sufficiently perturbed by small imperfections that their spectral positions (and distributions) are statistical in nature; and 3) overmoded region, where modes overlap, field distributions follow stochastic distributions, and the slot acts as if in free space. Understanding the penetration through slots in the overmoded region is of great interest, and is the focus of this article. Since full-wave solvers may not be able to provide a timely answer for very high frequencies due to a lack of memory and/or computation resources, we develop bounding methods to estimate worst-case average and maximum fields within the cavity. Finally, after discussing the bounding formulation, we compare its results to full-wave simulations at the first, second, and third resonance supported by the slot in the case of a cylindrical cavity. Note that the bounding formulation indicates that results are nearly independent of cavity shape: only the cavity volume, frequency, and cavity quality factor affect the overmoded region, making this formulation a powerful tool to assess electromagnetic interference and electromagnetic compatibility effects within cavities.

overmoded cavity↗
DOE OSTI2022

TRUST Nonlinear Dynamics (ND) Report: FY22

The objective of the Delivery Environments (DE) Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) work package, is to quantify and help increase confidence in specific areas of computational and experimental capabilities that are applicable to current and future delivery environments. More complete quantification of confidence in experimental and computational capabilities and the sufficient increase of confidence in those capabilities is critical to improving weapons engineering design, qualification, and assessment efforts that are relevant to the current and future stockpile. This work will use and provide feedback on analysis tools and experimental results databases for efficient and responsive engineering which are currently under development: engineering common model framework (ECMF), engineering quantification of margins and uncertainties (EQMU), and the test information management system (TIMS). Additionally, this work is being used to test a newly available framework in W-13, Weapons Analysis for Verified Engineering Simulation (WAVES).

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

6G Sparse Transmitarray

This powerpoint presents a sparse transmitarray design for 6G technology. The proposed approach leverages aperiodic arrays with large inter-element spacings, offering grating lobe-free operation, enhanced gain, and narrow beamwidths. The transmitarray topology, designed for sub-THz and THz frequencies, simplifies feed-network complexity, and minimizes losses. Numerical analysis and full wave simulation results show this approach is feasible.

Christian, Thomas Edward↗
DOE OSTI2020

Penetration through Slots in Cylindrical Cavities Operating at Fundamental Cavity Modes in the Presence of Electromagnetic Absorbers

Placing microwave absorbing materials into a high-quality factor resonant cavity may in general reduce the large interior electromagnetic fields excited under external illumination. In this study, we aim to combine two analytical models we previously developed: 1) an unmatched formulation for frequencies below the slot resonance to model shielding effectiveness versus frequency; and 2) a perturbation model approach to estimate the quality factor of cavities in the presence of absorbers. The resulting model realizes a toolkit with which design guidelines of the absorber's properties and location can be optimized over a frequency band. Analytic predictions of shielding effectiveness for three transverse magnetic modes for various locations of the absorber placed on the inside cavity wall show good agreement with both full-wave simulations and experiments, and validate the proposed model. This analysis opens new avenues for specialized ways to mitigate harmful fields within cavities.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Wildfire smoke triggers cirrus formation: lidar observations over the eastern Mediterranean

The number of intense wildfires may increase further in upcoming years as a consequence of climate change. It is therefore necessary to improve our knowledge about the role of smoke in the climate system, with emphasis on the impact of smoke particles on the evolution of clouds, precipitation, and cloud radiative properties. Presently, one key aspect of research is whether or not wildfire smoke particles can initiate cirrus formation. In this study, we present lidar observations over Limassol, Cyprus, from 27 October to 3 November 2020, when extended wildfire smoke fields crossed the Mediterranean Basin from Portugal to Cyprus. We found strong evidence that aged smoke (organic aerosol particles) originating from wildfires in North America triggered significant ice nucleation at temperatures from -47 to -53 °C and caused the formation of extended cirrus layers. The observations suggest that the ice crystals were nucleated just below the tropopause in the presence of smoke particles serving as ice-nucleating particles (INPs). The main part of the 2–3 km thick smoke layer was, however, in the lower stratosphere just above the tropopause. With actual radiosonde observations of temperature and relative humidity and lidar-derived smoke particle surface area concentrations used as starting values, gravity wave simulations show that the lofting of air by 100–200 m is sufficient to initiate significant ice nucleation on the smoke particles, leading to ice crystal number concentrations of 1–100 L -1 .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

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↗
DOE OSTI2021

Various facets of magnetic charge correlation: Micromagnetic and distorted-wave Born approximation simulations study

The emergent concept of the magnetic charge quasiparticle provides a new realm to study the evolution of magnetic properties in two-dimensional artificially frustrated magnets. Here we report on the exploration of magnetic phases due to various magnetic charge correlation using the complementary numerical techniques of micromagnetic and distorted-wave Born approximation simulations in artificial permalloy honeycomb lattice. The honeycomb element length varies between 10 nm and 100 nm, while the width and thickness are kept within the single domain limit. In addition to the charge ordered loop state, we observe disordered charge arrangement, characterized by the random distribution of ±Q charges, in single domain size honeycomb lattice. As the length of the honeycomb element increases, low multiplicity magnetic charges tend to form contiguous bands in thinner lattice. Thin honeycomb lattice with 100 nm element length exhibits a perfect spin ice pattern, which remains unaffected by the modest increase in the width of element size. We simulate scattering profiles under the pretext of distorted-wave Born approximation formalism for the micromagnetic phases. The results are expected to provide useful guidance in the experimental investigation of magnetic phases in an artificial honeycomb magnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Hybrid simulation of Alfvén wave parametric decay instability in a laboratory relevant plasma

Large-amplitude Alfvén waves are subject to parametric decays, which can have important consequences in space, astrophysical, and fusion plasmas. Although the Alfvén wave parametric decay instability (PDI) was predicted decades ago, observational evidence is limited, stimulating considerable interest in laboratory demonstration of the instability and associated numerical modeling. Here, we report an important step toward direct hybrid simulation of the Alfvén wave PDI in a laboratory plasma, using antenna-like wave injection of a circularly polarized wave and realistic wave-plasma parameters. Considering collisionless damping, we identify the threshold Alfvén wave amplitudes and frequencies required for triggering the instability. These threshold behaviors are corroborated by simple theoretical analysis, considering the balance between PDI growth and ion Landau damping and PDI development in a bounded plasma. Other effects not included in the present model such as finite transverse wave scale and ion–neutral collisions are briefly discussed. These hybrid simulations promise to be a useful tool for investigating laboratory Alfvén wave dynamics and may provide guidance for future laboratory demonstration of the PDI process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Recommendations on setup in simulating atmospheric gravity waves under conventionally neutral boundary layer conditions

Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.

17 WIND ENERGY↗
DOE OSTI2022

Direct Numerical Simulations of Turbulent Flow over Misaligned Traveling Waves

Direct numerical simulations of turbulent flow over prescribed traveling waves in a half-channel flow setup subject to a combined streamwise and spanwise pressure gradient are performed at a friction Reynolds number of 180. The simulations undertaken in this study consider both flow-aligned (including opposing and following) and misaligned waves as we attempt to quantify the effect of wave misalignment on turbulence statistics, the pressure drag force, and the wave attenuation rate. For the simulations, we consider three characteristic wave age values corresponding to slow-, intermediate-, and fast-moving waves. Wave misalignment is taken into account by applying a spanwise pressure gradient vertical to the traveling waves, which results in a three-dimensional turbulent flow field above the moving waves. Key flow quantities such as the mean velocity, velocity variances, and momentum fluxes are found to vary with the wave parameters, confirming the findings of previous studies.

direct numerical simulation↗
DOE OSTI2022

Quantum signal processing for simulating cold plasma waves

Numerical modeling of radio-frequency waves in plasma with sufficiently high spatial and temporal resolution remains challenging even with modern computers. However, such simulations can be sped up using quantum computers in the future. In this work, we propose how to do such modeling for cold plasma waves, in particular, for an X wave propagating in an inhomogeneous one-dimensional plasma. The wave system is represented in the form of a vector Schrödinger equation with a Hermitian Hamiltonian. Block encoding is used to represent the Hamiltonian through unitary operations that can be implemented on a quantum computer. To perform the modeling, we apply the so-called quantum signal processing algorithm and construct the corresponding circuit. Quantum simulations with this circuit are emulated on a classical computer, and the results show agreement with traditional classical calculations. We also discuss how our quantum circuit scales with the resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗