SEARCH · Engineering Papers
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.
Forum papers are thought-provoking opinion pieces or essays founded in fact, sometimes containing speculation, on a civil engineering topic of general interest and relevance to the readership of the journal. The views expressed in this Forum article do not necessarily reflect the views of ASCE or the Editorial Board of the journal.
Nonlinear WKB is a multiscale technique for studying locally plane-wave solutions of nonlinear partial differential equations (PDEs). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, here we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether’s theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by S 1 . As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow.
Magnetic reconnection, topological changes in magnetic fields, is a fundamental process in magnetized plasmas. It is associated with energy release in regions of magnetic field annihilation, but this is only one facet of this process. Astrophysical fluid flows normally have very large Reynolds numbers and are expected to be turbulent, in agreement with observations. In strong turbulence, magnetic field lines constantly reconnect everywhere and on all scales, thus making magnetic reconnection an intrinsic part of the turbulent cascade. We note in particular that this is inconsistent with the usual practice of magnetic field lines as persistent dynamical elements. A number of theoretical, numerical, and observational studies starting with the paper done by Lazarian and Vishniac [Astrophys. J. 517, 700–718 (1999)] proposed that 3D turbulence makes magnetic reconnection fast and that magnetic reconnection and turbulence are intrinsically connected. In particular, we discuss the dramatic violation of the textbook concept of magnetic flux-freezing in the presence of turbulence. We demonstrate that in the presence of turbulence, the plasma effects are subdominant to turbulence as far as the magnetic reconnection is concerned. The latter fact justifies a magnetohydrodynamiclike treatment of magnetic reconnection on all scales much larger than the relevant plasma scales. We discuss the numerical and observational evidence supporting the turbulent reconnection model. In particular, we demonstrate that the tearing reconnection is suppressed in 3D, and unlike the 2D settings, 3D reconnection induces turbulence that makes magnetic reconnection independent of resistivity. We show that turbulent reconnection dramatically affects key astrophysical processes, e.g., star formation, turbulent dynamo, and acceleration of cosmic rays. We provide criticism of the concept of “reconnection-mediated turbulence” and explain why turbulent reconnection is very different from enhanced turbulent resistivity and hyper-resistivity and why the latter have fatal conceptual flaws.
The elementary basis of intelligence in organisms with a central nervous system includes neurons and synapses and their complex interconnections forming neural circuits. In non-neural organisms such as slime mold with gel-like media, viscosity modulation enables adaptation to changing environments. At a larger scale, collective intelligence emerges via social interactions and feedback in animal colonies. Learning and memory are therefore multi-scale features that evolve as a result of constant interactions with the environment. There is growing interest in emulating such features of intelligence in computing machines and autonomous systems. Materials that can respond to their environment in a manner similar to organisms (referred to as "organismic materials") therefore may be of interest as hardware components in artificial intelligence machines. In this brief review, we present a class of semiconductors called correlated oxides as candidates for learning machines. The term "correlated" refers to the fact that electrons in such lattices strongly interact and the ground state is not what is predicted by classical band theory. Such materials can undergo insulator-metal transitions at near ambient conditions under external stimuli such as thermal or electrical fields, strain, and chemical doping. Depending on the mechanism driving the transition, intermediate states can be metastable with different volatilities, and the time scales of phase change can be controlled over many orders of magnitude. The change in electronic properties can be sharp or gradual, leading to digital or analog behavior. These properties enable the realization of artificial neurons and synapses and emulate the associative and non-associative learning characteristics found in various organisms. We examine microscopic properties concerning electronic and structural transitions leading to collective behavior and theoretical treatments of the ground state and dynamical response, showcasing VO 2 as a model system. Next, we briefly review algorithms designed from the plasticity demonstrated by phase changing systems. Finally, we conclude the brief review with suggestions for future research toward realizing non-von Neumann machines.
This work presents a novel method of obtaining in situ strain measurements at high temperature by simultaneous digital image correlation (DIC), which provides the total strain on the specimen surface, and synchrotron x-ray diffraction (XRD), which provides lattice strains of crystalline materials. DIC at high temperature requires specialized techniques to overcome the effects of increased blackbody radiation that would otherwise overexpose the images. The technique presented herein is unique in that it can be used with a sample enclosed in an infrared heater, remotely and simultaneously with synchrotron XRD measurements. The heater included a window for camera access, and the light of the heater lamps is used as illumination. High-temperature paint is used to apply a random speckle pattern to the sample to allow the tracking of displacements and the calculation of the DIC strains. An inexpensive blue theatrical gel filter is used to block interfering visible and infrared light at high temperatures. This technique successfully produces properly exposed images at 870 °C and is expected to perform similarly at higher temperatures. The average strains measured by DIC were validated by an analytical calculation of the theoretical strain. Simultaneous DIC and XRD strain measurements of Inconel 718 (IN718) tensile test specimens were performed under thermal and mechanical loads and evaluated. Finally, this approach uses the fact that with DIC, the total strain is measured, including plastic strain, while with XRD, only elastic strain is captured. The observed differences were discussed with respect to the effective deformation mechanisms.
Nanoshell quantum dots (QDs) represent a novel class of colloidal semiconductor nanocrystals (NCs), which supports tunable optoelectronic properties over the extended range of particle sizes. Traditionally, the ability to control the bandgap of colloidal semiconductor NCs is limited to small-size nanostructures, where photoinduced charges are confined by Coulomb interactions. A notorious drawback of such a restricted size range concerns the fact that assemblies of smaller nanoparticles tend to exhibit a greater density of interfacial and surface defects. This presents a potential problem for device applications of semiconductor NCs where the charge transport across nanoparticle films is important, as in the case of solar cells, field-effect transistors, and photoelectrochemical devices. Furthermore, the morphology of nanoshell QDs addresses this issue by enabling the quantum-confinement in the shell layer, where two-dimensional excitons can exist, regardless of the total particle size. Such a geometry exhibits one of the lowest surface-to-volume ratios among existing QD architectures and, therefore, could potentially lead to improved charge-transport and multi-exciton characteristics. The expected benefits of the nanoshell architecture were recently demonstrated by a number of reports on the CdS bulk /CdSe nanoshell model system, showing an improved photoconductivity of solids and increased lifetime of multi-exciton populations. Along these lines, this perspective will summarize the recent work on CdS bulk /CdSe nanoshell colloids and discuss the possibility of employing other nanoshell semiconductor combinations in light-harvesting and lasing applications.
The spectral properties associated with laminar, anti-parallel reconnection are examined using a 2.5D kinetic particle in cell simulation. Both the reconnection rate and the energy spectrum exhibit three distinct phases: an initiation phase where the reconnection rate grows, a quasi-steady phase, and a declining phase where both the reconnection rate and the energy spectrum decrease. During the steady phase, the energy spectrum exhibits approximately a double power law behavior, with a slope near −5/3 at wave numbers smaller than the inverse ion inertial length and a slope steeper than −8/3 for larger wave numbers up to the inverse electron inertial length. This behavior is consistent with a Kolmogorov energy cascade and implies that laminar reconnection may fundamentally be an energy cascade process. Consistent with this idea is the fact that the reconnection rate exhibits a rough correlation with the energy spectrum at wave numbers near the inverse ion inertial length. The 2D spectrum is strongly anisotropic with most energy associated with the wave vector direction normal to the current sheet. Reconnection acts to isotropize the energy spectrum, reducing the Shebalin angle from an initial value of 70° to about 48° (nearly isotropic) by the end of the simulation. The distribution of energy over length scales is further analyzed by dividing the domain into spatial subregions and employing structure functions.
Recent reports of the superconductivity in hydrides of two different families (covalent lattice, as in SH3 and clathrate-type H-cages containing La and Y atoms, as in LaH10 and YH6) have revealed new families of high-Tc materials with Tc’s near room temperature values. These findings confirm earlier expectations that hydrides may have very high Tc’s due to the fact that light H atoms have very high vibrational frequencies, leading to high Tc values within the conventional Bardeen–Cooper–Schrieffer phonon mechanism of superconductivity. However, as is pointed out by Ashcroft, it is important to have the metallic hydrogen “alloyed” with the elements added to it. This concept of a metallic alloy containing a high concentration of metal-like hydrogen atoms has been instrumental in finding new high-Tc superhydrides. These new superhydride “room-temperature” superconductors are stabilized only at very high pressures above 100 GPa, making the experimental search for their superconducting properties very difficult. We will review the current experimental and theoretical results for LaH10−x and YH6−x superhydrides.
Multiferroic BiFeO$_3$ (BFO) shows several phonon modes at infrared (IR) to THz energies, which are expected to carry information on any sample property coupled to crystal lattice vibrations. While macroscopic IR studies of BFO are often limited by single-crystal size, scattering-type scanning near-field optical microscopy (s-SNOM) allows for IR thin film spectroscopy of nanoscopic probing volumes with negligible direct substrate contribution to the optical signal. In fact, polaritons such as phonon polaritons of BFO introduce a resonant tip–sample coupling in s-SNOM, leading to both stronger signals and enhanced sensitivity to local material properties. Here in this paper, we explore the near-field response of BFO thin films at three consecutive resonances (centered around 5 THz, 13 THz, and 16 THz), by combining s-SNOM with a free-electron laser. We study the dependence of these near-field resonances on both the wavelength and tip–sample distance. Enabled by the broad spectral range of the measurement, we probe phonon modes connected to the predominant motion of either the bismuth or oxygen ions. Therefore, we propose s-SNOM at multiple near-field resonances as a versatile and very sensitive tool for the simultaneous investigation of various sample properties.
In recent years, several magnetic Mott insulators with strong spin-orbit coupling were suggested to be proximate to the Kitaev quantum spin liquid, whose one of the most exciting features is the fractionalization of spin excitations into itinerant Majorana fermions and static Z 2 fluxes. Unfortunately, the ground states of these systems cannot be easily captured by experiment, remaining featureless to conventional local probes. Here, we propose to study the signatures of fractionalized excitations by exploiting their coupling to the lattice vibrations, dubbed magnetoelastic coupling, which arises from the fact that the interaction between spins depends on the relative distance between them. We argue that the magnetoelastic coupling can lead to the distinct modification of the phonon dynamics, which can be observed by measuring renormalized phonon spectrum, the sound attenuation, and the phonon Hall viscosity. This makes the phonon dynamics a promising tool for the characterization and identification of quantum spin liquid phases. In this work, we focus on the magnetoelastic effects in the three-dimensional Kitaev model realized on the hyperhoneycomb lattice. Here, the hyperhoneycomb Kitaev spin liquid is particularly interesting since the strong Kitaev interaction was observed in the Kitaev magnet β-Li 2 IrO 3 , for which the spin-orbit entangled J eff = 1/2 moments of iridium ions form precisely the hyperhoneycomb lattice.
Shear localization is characterized by the concentration of plastic deformation within narrow bands of material and constitutes the precursor to failure in materials subjected to high strain rates. These bands have been observed experimentally, but classical continuum material models have not been able to fully capture their development. The aim of this research is to investigate the role of microstructural differences in the dynamic shear response of wrought and additively manufactured (AM) 316L SS. Compact Forced Simple Shear specimens for Digital Image Correlation (CFSS-DIC) of wrought and additively manufactured 316L SS were tested on a Split Hopkinson Pressure Bar (SHPB) in conjunction with infrared thermometry measurements. Results indicate that the microstructures of wrought and AM 316L SS exhibit differences that affect the shear response of this materials. In conclusion, CFSS-DIC samples tested on the SHPB did not exhibit recrystallization inside the shear region, fact confirmed by thermometry measurements.
Cross beam energy transfer (CBET), invoked by setting a wavelength difference, Δλ, between inner and outer beam cones, can be used to increase the drive on the waist in indirectly driven inertial confinement fusion experiments at the National Ignition Facility (NIF). Historically, hot spot symmetry control in capsule implosions in high (≥0.9 mg/cm 3 4 He) gas fill Hohlraums was enabled by substantial CBET. However, these implosion designs suffered from inflight symmetry swings, high SRS backscatter on the inner cones, and significant hot electron generation posing a threat to DT fuel preheat. Subsequent experiments in larger, low (≤0.6 mg/cm 3 4 He) gas fill Hohlraums demonstrated round implosions by varying the inner cone fraction throughout the laser drive at Δλ = 0 Å while keeping backscatter and hot electron generation very low. To enable driving larger capsules at a given Hohlraum size, additional tools for implosion symmetry control are required. With this goal in mind, here we present a detailed experimental study of using CBET in low gas fill Hohlraums near NIF's current peak power capability. We find a ~2.5× higher sensitivity of the P 2 Legendre mode with respect to Δλ changes compared to that of high gas fill designs. We attribute this observation to the fact that backscatter remains very low and that CBET remains in a linear regime, as suggested by simulations. As a result, a much smaller Δλ of order 1 Å is sufficient for sustaining implosion symmetry while keeping laser-to- Hohlraum coupling high and hot electron generation very low. While this study used plastic ablator capsules, our findings can be generalized to other ablator materials and, hence, show great promise for using wavelength detuning as a strong lever for implosion symmetry control in future low gas fill designs that require smaller case to capsule ratios in order to increase the energy coupled to the capsule.
The physical processes taking place at the separatrix and scrape-off layer regions are crucial for the operation of tokamaks as they govern the interaction of hot plasma with the vessel walls. Numerical modeling of the edge with state-of-the-art codes attempts to elucidate the complex interactions between neoclassical drifts, turbulence, poloidal and parallel flows that control the physical set-up of the SOL region. Here, we present post-processing analysis of simulation results from the gyrokinetic code XGC1, comparing and contrasting edge turbulence characteristics from a simulation of the DIII-D tokamak against a simulation of the Alcator C-Mod tokamak. We find that the equilibrium ExB flux across the separatrix has a similar poloidal pattern in both discharges which can be explained by rB-drifts and trapped ion excursions. Furthermore, collisionality is noted to play a major role in the way that it prevents local charge accumulations from having more global effects in the C-Mod case. In both cases, turbulent electron heat flux is observed to be higher than the ion one and is possibly related to the need of electrons to maintain quasineutrality through the only channel available to them for exiting the confinement. By Fourier analysis, we identify turbulent frequencies and growth rates of the dominant mode in both simulations. In the case of C-Mod, these numbers point to the presence of a drift wave. In the DIII-D case, further linear simulations with the Gene code reveal a trapped electron mode. Furthermore, using a blob detection and tracking tool, we present the amplitude and size distributions of the blobs from both simulations. The amplitude distributions are in qualitative agreement with experimental observations while the size distributions are consistent with the fact that most of the blobs are not connecting to the divertor plates and suggest that they are generated by the shearing of the turbulent modes.
The structure of static MHD equilibria that admit continuous families of Euclidean symmetries is well understood. Such field configurations are governed by the classical Grad–Shafranov equation, which is a single elliptic partial differential equation in two space dimensions. By revealing a hidden symmetry, here we show that in fact all smooth solutions of the equilibrium equations with non-vanishing pressure gradients away from the magnetic axis satisfy a generalization of the Grad–Shafranov equation. In contrast to solutions of the classical Grad–Shafranov equation, solutions of the generalized equation are not automatically equilibria, but instead only satisfy force balance averaged over the one-parameter hidden symmetry. We then explain how the generalized Grad–Shafranov equation can be used to reformulate the problem of finding exact three-dimensional smooth solutions of the equilibrium equations as finding an optimal volume-preserving symmetry.
The Hartree–Fock problem provides the conceptual and mathematical underpinning of a large portion of quantum chemistry. As efforts in quantum technology aim to enhance computational chemistry algorithms, the Hartree–Fock method, central to many other numerical approaches, is a natural target for quantum enhanced algorithms. While quantum computers and quantum simulation offer many prospects for the future of modern chemistry, the non-deterministic polynomial-complete Hartree–Fock problem is not a likely candidate. Here, we highlight this fact from a number of perspectives including computational complexity, practical examples, and the full characterization of energy landscapes for simple systems.
Through an extensive data analysis of powder X-ray diffraction data obtained at pressures up to at least 78 GPa, we report the experimental equations of state for all lanthanide nitrides between NdN and LuN, excluding the radioactive Pm. By fitting the obtained unit cell volumes as a function of pressure with the third order Birch–Murnaghan equation of state, we find that the bulk modulus increases with an increasing lanthanide number from K 0 = 146(12) GPa for NdN to 182(7) GPa in EuN. Hereafter, the bulk modulus reaches a plateau for the rest of the series except for TmN which has a lower bulk modulus. We find that the first derivative of the bulk modulus is around 4 for all compounds except TbN, which displays a significantly different compression behavior. In addition, we find a B1 to B10 pressure-induced phase transition in NdN, SmN, EuN, and GdN at increasingly higher pressures. In fact, we observe that the onset pressure of the phase transition increases linearly with Ln atomic number. From TbN and onwards, we do not observe any sign of a B1 to B10 transition indicating that the transition pressure exceeds the studied pressure range. Therefore, we believe that, for the heavier lanthanides, the linear relationship between the onset pressure and the lanthanide number does not hold and even higher pressures are needed to observe the transition. This coherent study of the series of lanthanide nitrides offers a unique opportunity for benchmark studies of computational methods applied to compounds with 4f electrons.
A toroidal asymptotic matching model of the response of a tokamak plasma to a static resonant magnetic perturbation (RMP) is used to simulate the n = 3 RMP-induced edge-localized-mode-suppression windows in q 95 that are evident when the plasma current is slowly ramped in DIII-D discharge #145380. All quantities employed in the simulation are derived from experimental measurements, apart from the neutral particle data. Three cases are considered. In the first case, the natural frequencies of tearing modes resonant in the plasma are determined by the ion flows at the corresponding resonant surfaces, which is the prediction of nonlinear tearing mode theory. In the second case, the natural frequencies are determined by the local E x B velocities at the resonant surfaces. In the third case, the natural frequencies are determined by the electron flows at the resonant surfaces, which is the prediction of linear tearing mode theory. The second case gives the best agreement between the simulations and the experimental observations. The first and third cases only lead to partial agreement between the simulations and the observations. In the first case, the lack of complete agreement may be a consequence of using an inaccurate assumption for the neutral particle distribution in the pedestal. In the third case, the lack of complete agreement is probably due to the fact that the response of a tokamak plasma to an RMP is not accurately described by linear tearing mode theory.
Finite temperature auxiliary field-based quantum Monte Carlo methods, including determinant quantum Monte Carlo and Auxiliary Field Quantum Monte Carlo (AFQMC), have historically assumed pivotal roles in the investigation of the finite temperature phase diagrams of a wide variety of multidimensional lattice models and materials. Despite their utility, however, these techniques are typically formulated in the grand canonical ensemble, which makes them difficult to apply to condensates such as superfluids and difficult to benchmark against alternative methods that are formulated in the canonical ensemble. Working in the grand canonical ensemble is furthermore accompanied by the increased overhead associated with having to determine the chemical potentials that produce desired fillings. Given this backdrop, in this work, we present a new recursive approach for performing AFQMC simulations in the canonical ensemble that does not require knowledge of chemical potentials. To derive this approach, we exploit the convenient fact that AFQMC solves the many-body problem by decoupling many-body propagators into integrals over one-body problems to which non-interacting theories can be applied. We benchmark the accuracy of our technique on illustrative Bose and Fermi–Hubbard models and demonstrate that it can converge more quickly to the ground state than grand canonical AFQMC simulations. We believe that our novel use of HS-transformed operators to implement algorithms originally derived for non-interacting systems will motivate the development of a variety of other methods and anticipate that our technique will enable direct performance comparisons against other many-body approaches formulated in the canonical ensemble.