SEARCH · Engineering Papers

Results for “wave fields”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Effects of Turbulence and Temperature Fluctuations on Knock Development in an Ethanol/Air Mixture

The effects of turbulence on knock development and intensity for a thermally inhomogeneous stoichiometric ethanol/air mixture at a representative end-gas autoignition condition in internal combustion engines are investigated using direct numerical simulations with a skeletal reaction mechanism. Two- and three-dimensional simulations are performed by varying the most energetic length scale of temperature, $l_T$ , and its relative ratio with the most energetic length scale of turbulence, $l_T/l_e$ , together with two different levels of the turbulent velocity fluctuation, $u'$ . It is found that $l_T$ / $l_e$ and the ratio of ignition delay time to eddy-turnover time, $\tau _{ig}/\tau _t$ , are the key parameters that control the detonation development. An increase in either $l_T$ or $l_e$ enhances the detonation propensity by allowing a longer run-up distance for the detonation development. The characteristic length scale of the temperature field, $l_T$ , is significantly modified by high turbulence intensity achieved by a large $l_e$ and $u'$ . The intense turbulence mixing effectively distributes the initial temperature field to broader scales to support the developing detonation waves, thereby increasing the likelihood of the detonation formation. On the contrary, high turbulence intensity with a short mixing time scale, achieved by a small $l_e$ and a large $u'$ , reduces the super-knock intensity attributed to the finer broken-up structures of detonation waves. Either $\tau _{ig}/\tau _t$ less than unity or $l_e = l_T$ even with a large $u'$ is found to have no significant effect on super-knock mitigation. Finally, high turbulent intensity may induce high-pressure spikes comparable to the von Neumann spike. Overall, increased temperature and pressure by combustion heating, noticeably after the peak of heat release rate, significantly enhance the collision and interaction of multiple emerging autoignition fronts near the ending combustion process, resulting in localized high-pressure spikes.

42 ENGINEERING↗

DOE OSTI2022

Topological d -wave superconductivity in two dimensions

Despite intensive searches for topological superconductors, the realization of topological superconductivity remains under debate. Previous proposals for the topological s-wave, p-wave, and chiral d-wave superconductivity have both advantages and disadvantages. In this review, we discuss two-dimensional topological superconductivity based on the non-chiral d-wave superconductors. It is shown that the noncentrosymmetric d-wave superconductors become topological superconductors under an infinitesimal Zeeman field without fine-tuning of parameters. Additionally, Floquet engineering for introducing the Zeeman field in a controllable way is also proposed. When the two-dimensional noncentrosymmetric superconductors are stacked to recover the global inversion symmetry, the field-induced parity transition may occur, and the high-field odd-parity superconducting state realizes various topological phases depending on the stacking structures. Two-dimensional heterostructures of strongly correlated electron systems, which have been developed by recent experiments, are proposed as a platform of the high-temperature topological superconductivity and the interplay of topology and strong correlations in superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Unsteady aerodynamics of large-scale floating offshore wind turbines in surge motion

Unsteady aerodynamic loads significantly influence the design and wake flow field of floating offshore wind turbines, especially due to wave- or vibration-induced tower top motions triggering various unsteady phenomena. Recent studies show that increasing turbine sizes amplify unsteady aerodynamic effects, as their impact typically grows with rotor diameter. This work combines recent findings from experiments and simulations on model-scale FOWT aerodynamics with new numerical analyses of large-scale rotors, providing a comprehensive understanding of unsteady phenomena occurrence and impact. Numerical analyses of the IEA 15-MW and 22-MW rotors undergoing surge motions characterise the combined influence of motion-induced unsteady phenomena on rotor thrust. Results indicate that unsteady effects can reduce thrust force variations by up to 40% at realistic surge periods. These findings contrast with prior model-scale rotor investigations, attributed to the specialised design of the model rotors. Comparisons between numerical methods - the dynamic blade element momentum method in OpenFAST and free vortex wake modules in panMARE and OpenFAST - reveal persistent differences in thrust predictions under both idealised and realistic conditions. This highlights that classical blade element momentum approaches require enhancement to accurately capture unsteady loads on large-scale floating offshore wind turbines at low wind speeds.

16 TIDAL AND WAVE POWER↗

DOE OSTI2021

Multiple fermion scattering in the weakly coupled spin-chain compound YbAlO 3

The Heisenberg antiferromagnetic spin-1/2 chain, originally introduced almost a century ago, is one of the best studied models in quantum mechanics due to its exact solution, but nevertheless it continues to present new discoveries. Its low-energy physics is described by the Tomonaga-Luttinger liquid of spinless fermions, similar to the conduction electrons in one-dimensional metals. In this work we investigate the Heisenberg spin-chain compound YbAlO3 and show that the weak interchain coupling causes Umklapp scattering between the left- and right-moving fermions and stabilizes an incommensurate spin-density wave order at q = 2k F under finite magnetic fields. These Umklapp processes open a route to multiple coherent scattering of fermions, which results in the formation of satellites at integer multiples of the incommensurate fundamental wavevector Q = nq. Our work provides surprising and profound insight into bandstructure control for emergent fermions in quantum materials, and shows how neutron diffraction can be applied to investigate the phenomenon of coherent multiple scattering in metals through the proxy of quantum magnetic systems.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Intrinsic nature of chiral charge order in the kagome superconductor Rb V 3 Sb 5

Superconductors with kagome lattices have been identified for over 40 years, with a superconducting transition temperature T c up to 7 K. Recently, certain kagome superconductors have been found to exhibit an exotic charge order, which intertwines with superconductivity and persists to a temperature being one order of magnitude higher than T c . Here, we use scanning tunneling microscopy to study the charge order in kagome superconductor Rb V 3 Sb 5 . We observe both a 2×2 chiral charge order and nematic surface superlattices (predominantly 1×4). We find that the 2×2 charge order exhibits intrinsic chirality with magnetic field tunability. Defects can scatter electrons to introduce standing waves, which couple with the charge order to cause extrinsic effects. While the chiral charge order resembles that discovered in KV 3 Sb 5 , it further interacts with the nematic surface superlattices that are absent in KV 3 Sb 5 but exist in CsV 3 Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Elastic scattering of 3 He + 4 He with the SONIK scattering chamber

Measurements of the elastic scattering cross section of 3 He and 4 He are important in order to improve constraints on theoretical models of 4 He ⁢( 3 He ,𝛾)⁢ 7 Be, a key reaction in Big Bang nucleosynthesis and solar neutrino production. The astrophysical 𝑆 factor for this reaction is a significant source of uncertainty in the standard-solar-model prediction of the 7 Be and 8 B solar neutrino fluxes. The elastic scattering measurements reported in the literature do not extend to low energies and lack proper uncertainty quantification. A new measurement of the 4 He ⁢( 3 He, 3 He)⁢ 4 He reaction has been made at center-of-mass energies 𝐸c.m.=0.38−3.13 MeV using the Scattering of Nuclei in Inverse Kinematics (SONIK) scattering chamber: a windowless, extended gas target surrounded by an array of 30 collimated silicon charged particle detectors situated at TRIUMF. This is the first elastic scattering measurement of 3 He + 4 He made below 500 keV and it has greater angular range and better precision than previous measurements. The elastic scattering data were analyzed using both 𝑅-matrix and halo effective field theory frameworks, and values of the 𝑠-wave scattering length and effective range were extracted. Here, the resulting improvement in knowledge of the 𝑠-wave effective-range function at low energies reduces the overall uncertainty in 𝑆 34 at solar energies.

Elastic scattering reactions↗

DOE OSTI2026

Codimension-two spiral spin liquid in the effective honeycomb-lattice compound Cs 3 ⁢Fe 2 ⁢Cl 9

A codimension-two spiral spin liquid is a correlated paramagnetic state with one-dimensional ground state degeneracy hosted within a three-dimensional lattice. Here, in this work, via neutron scattering experiments and numerical simulations, we establish the existence of a codimension-two spiral spin liquid in the effective honeycomb-lattice compound Cs 3 ⁢Fe 2⁢ Cl 9 , which demonstrates an alternate path to spiral spin liquids by overcoming the long-standing impediment of weak further-neighbor interactions. In the long-range ordered regime, competing spiral and spin density wave orders emerge as a function of applied magnetic field, among which a possible order-by-disorder transition is identified.

Monte Carlo methods↗

DOE OSTI2022

First demonstration of tuning between the Kitaev and Ising limits in a honeycomb lattice

Recent observations of novel spin-orbit coupled states have generated interest in 4d/5d transition metal systems. A prime example is the J eff = $\frac{1}{2}$ state in iridate materials and α-RuCl 3 that drives Kitaev interactions. Here, by tuning the competition between spin-orbit interaction (λ SOC ) and trigonal crystal field (Δ T ), we restructure the spin-orbital wave functions into a previously unobserved μ = $\frac{1}{2}$ state that drives Ising interactions. This is done via a topochemical reaction that converts Li 2 RhO 3 to Ag 3 LiRh 2 O 6 . Using perturbation theory, we present an explicit expression for the μ = $\frac{1}{2}$ state in the limit Δ T ≫ λ SOC realized in Ag 3 LiRh 2 O 6 , different from the conventional J eff = $\frac{1}{2}$ state in the limit λ SOC ≫ Δ T realized in Li 2 RhO 3 . The change of ground state is followed by a marked change of magnetism from a 6 K spin-glass in Li 2 RhO 3 to a 94 K antiferromagnet in Ag 3 LiRh 2 O 6 .

36 MATERIALS SCIENCE↗

DOE OSTI2025

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

DOE OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

DOE OSTI2020

Lorentz Symmetry Group, Retardation, Intergalactic Mass Depletion and Mechanisms Leading to Galactic Rotation Curves

The general theory of relativity (GR) is symmetric under smooth coordinate transformations, also known as diffeomorphisms. The general coordinate transformation group has a linear subgroup denoted as the Lorentz group of symmetry, which is also maintained in the weak field approximation to GR. The dominant operator in the weak field equation of GR is thus the d’Alembert (wave) operator, which has a retarded potential solution. Galaxies are huge physical systems with dimensions of many tens of thousands of light years. Thus, any change at the galactic center will be noticed at the rim only tens of thousands of years later. Those retardation effects are neglected in the present day galactic modelling used to calculate rotational velocities of matter in the rims of the galaxy and surrounding gas. The significant differences between the predictions of Newtonian instantaneous action at a distance and observed velocities are usually explained by either assuming dark matter or by modifying the laws of gravity (MOND). In this paper, we will show that, by taking general relativity seriously without neglecting retardation effects, one can explain the radial velocities of galactic matter in the M33 galaxy without postulating dark matter. It should be stressed that the current approach does not require that velocities v are high; in fact, the vast majority of galactic bodies (stars, gas) are substantially subluminal—in other words, the ratio of $\frac{v}{c} \ll 1$. Typical velocities in galaxies are 100 km/s, which makes this ratio 0.001 or smaller. However, one should consider the fact that every gravitational system, even if it is made of subluminal bodies, has a retardation distance, beyond which the retardation effect cannot be neglected. Every natural system, such as stars and galaxies and even galactic clusters, exchanges mass with its environment, for example, the sun loses mass through solar wind and galaxies accrete gas from the intergalactic medium. This means that all natural gravitational systems have a finite retardation distance. The question is thus quantitative: how large is the retardation distance? For the M33 galaxy, the velocity curve indicates that the retardation effects cannot be neglected beyond a certain distance, which was calculated to be roughly 14,000 light years; similar analysis for other galaxies of different types has shown similar results. We demonstrate, using a detailed model, that this does not require a high velocity of gas or stars in or out of the galaxy and is perfectly consistent with the current observational knowledge of galactic and extra galactic material content and dynamics

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Digital Q-Meter for continuous-wave NMR

A method to perform continuous-wave NMR measurements of nuclear magnetization at high magnetic fields, above 2.5 T, without analog down-mixing is described. An FPGA controls a digital clock pulse which is used to stimulate a resonant circuit and provide a reference signal. An algorithm determines the real portion of a resonant circuit signal near the Larmor frequency of the species of interest using only two measurements of the waveform per cycle. The FPGA automatically alters a variable capacitance to tune the resonant circuit to the Larmor frequency.

Maxwell, James↗

DOE OSTI2023

Comparative investigation of iterative solutions of the all-frequency stable formulation and its vector-potential-only variation

An all-frequency stable formulation has been proposed and applied to solve low-frequency and multiscale electromagnetic problems to avoid the low-frequency breakdown catastrophe of the traditional vector wave equation. Based on the potential representation of the fields, this formulation involves vector and scalar potentials, as well as an auxiliary potential to enforce the gauge condition. In this paper, such a formulation is first simplified to involve vector potential only. Their iterative solutions are then sought by using an incomplete LU preconditioner and an iterative solver. The iterative performance of both formulations are investigated in a comparative study.

Mekonnen, Minyechil↗

DOE OSTI2022

Effective Model for Fractional Topological Corner Modes in Quasicrystals

High-order topological insulators (HOTIs), as generalized from topological crystalline insulators, are characterized with lower-dimensional metallic boundary states protected by spatial symmetries of a crystal, whose theoretical framework based on band inversion at special k points cannot be readily extended to quasicrystals because quasicrystals contain rotational symmetries that are not compatible with crystals, and momentum is no longer a good quantum number. In this work, we develop a low-energy effective model underlying HOTI states in 2D quasicrystals for all possible rotational symmetries. By implementing a novel Fourier transform developed recently for quasicrystals and approximating the long-wavelength behavior by their large-scale average, we construct an effective k · p Hamiltonian to capture the band inversion at the center of a pseudo-Brillouin zone. We show that an in-plane Zeeman field can induce mass kinks at the intersection of adjacent edges of a 2D quasicrystal topological insulators and generate corner modes (CMs) with fractional charge, protected by rotational symmetries. Our model predictions are confirmed by numerical tight-binding calculations. Furthermore, when the quasicrystal is proximitized by an s-wave superconductor, Majorana CMs can also be created by tuning the field strength and chemical potential. Our work affords a generic approach to studying the low-energy physics of quasicrystals, in association with topological excitations and fractional statistics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Magnetoelastic distortion of multiferroic BiFeO 3 in the canted antiferromagnetic state

Using THz spectroscopy, here we show that the spin-wave spectrum of multiferroic BiFeO 3 in its high-field canted antiferromagnetic state is well described by a spin model that violates rhombohedral symmetry. We demonstrate that the monoclinic distortion of the canted antiferromagnetic state is induced by the single-ion magnetoelastic coupling between the lattice and the two nearly antiparallel spins. The revised spin model for BiFeO 3 contains two new single-ion anisotropy terms that violate rhombohedral symmetry and depend on the direction of the magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Edge supercurrent reveals competition between condensates in a Weyl superconductor

In topological materials, the edge states are readily distinguished from the bulk states. The situation where a topological semimetal becomes superconducting so that Cooper pairs occupy both the bulk and the edge states is not well understood. In particular, we do not know if we can force their pairing symmetries to be different. Here we show that, when supercurrent is injected into the superconducting Weyl semimetal MoTe 2 from Nb contacts, the invasive s-wave pairing potential from Nb is incompatible with the intrinsic Cooper pair condensate in MoTe 2 . This incompatibility leads to strong stochasticity in the switching current and an unusual anti-hysteretic behaviour in the current–voltage loops. There is also an asymmetry in the edge oscillations where, as the magnetic field crosses zero, the phase noise switches from one with a noisy spectrum to one that is noise free. Furthermore, using the noise spectrum as a guide, we track the anomalous features to field-induced switching of the device gap function between s-wave symmetry and the unconventional symmetry intrinsic to MoTe 2 . We infer that the behaviour of the gap function along the edges is different from that in the bulk.

Superconducting properties and materials↗

DOE OSTI2024

Slow Electron Holes in the Earth's Magnetosheath

Here, we present a statistical analysis of electrostatic solitary waves observed aboard Magnetospheric Multiscale spacecraft in the Earth's magnetosheath. Applying single-spacecraft interferometry to several hundred solitary waves collected in about 2-minute interval, we show that almost all of them have the electrostatic potential of positive polarity and propagate quasi-parallel to the local magnetic field with plasma frame velocities of the order of 100 km/s. The solitary waves have typical parallel half-widths from 10 to 100 m that is between 1 and 10 Debye lengths and typical amplitudes of the electrostatic potential from 10 to 200 mV that is between 0.01% and 1% of local electron temperature. The solitary waves are associated with quasi-Maxwellian ion velocity distribution functions, and their plasma frame velocities are comparable with ion thermal speed and well below electron thermal speed. We argue that the solitary waves of positive polarity are slow electron holes and estimate the time scale of their acceleration, which occurs due to interaction with ions, to be of the order of one second. The observation of slow electron holes indicates that their lifetime was shorter than the acceleration time scale. We argue that multi-spacecraft interferometry applied previously to these solitary waves is not applicable because of their too-short spatial scales. The source of the slow electron holes and the role in electron-ion energy exchange remain to be established.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Analytical model and dynamical phase-field simulation of terahertz transmission across ferroelectrics

Here, we theoretically investigate the steady-state transmission of continuous terahertz (THz) wave across a freestanding ferroelectric slab. Based on the Landau-Ginzburg-Devonshire theory of ferroelectrics and the coupled equations of motion for polarization and electromagnetic (EM) waves, we derive the analytical expressions of the frequency- and thickness-dependent dielectric susceptibility and transmission coefficient at the thin slab limit in the harmonic excitation regime. When the slab thickness is much smaller than the THz wavelength in the ferroelectric, the analytical predictions agree well with the numerical simulations from a dynamical phase-field model that incorporates the coupled dynamics of strain, polarization, and EM wave in multiphase systems. At larger thicknesses, the transmission is mainly determined by the frequency-dependent attenuation of THz waves in the ferroelectric and the formation of a standing polarization/THz wave. Our results advance the understanding of the interaction between THz wave and ferroelectrics and suggest the potential of exploiting ferroelectrics to achieve low-heat-dissipation, nonvolatile voltage modulation of THz transmission for high-data-rate wireless communication.

36 MATERIALS SCIENCE↗