Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Notes”

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 1,765 records · Page 98

False metals, real insulators, and degenerate gapped metals

This paper deals with a significant family of compounds predicted by simplistic electronic structure theory to be metals but are, in fact, insulators. This false metallic state has been traditionally attributed in the literature to reflect the absence of proper treatment of electron-electron correlation (“Mott insulators”) whereas, in fact, even mean-field like density functional theory describes the insulating phase correctly if the restrictions posed on the simplistic theory are avoided. Such unwarranted restrictions included different forms of disallowing symmetry breaking described in this article. As the science and technology of conductors have transitioned from studying simple elemental metals such as Al or Cu to compound conductors such as binary or ternary oxides and pnictides, a special class of degenerate but gapped metals has been noticed. Their presumed electronic configurations show the Fermi level inside the conduction band or valence band, yet there is an “internal band gap” between the principal band edges. The significance of this electronic configuration is that it might be unstable toward the formation of states inside the internal band gap when the formation of such states costs less energy than the energy gained by transferring carriers from the conduction band to these lower energy acceptor states, changing the original (false) metal to an insulator. The analogous process also exists for degenerate but gapped metals with the Fermi level inside the valence band, where the energy gain is defined by transfer of electrons from the donor level to the unoccupied part of the valence band. We focus here on the fact that numerous electronic structure methodologies have overlooked some physical factors that could stabilize the insulating alternative, predicting instead false metals that do not really exist (note that this is in general not a physical phase transition, but a correction of a previous error in theory that led to a false prediction of a metal). Such errors include: (i) ignoring spin symmetry breaking, such as disallowing magnetic spin ordering in CuBi2O4 or disallowing the formation of polymorphous spin networks in paramagnetic LaTiO 3 and YTiO 3 ; (ii) ignoring structural symmetry breaking, e.g., not enabling energy-lowering bond disproportionation (Li-doped TiO 2 , SrBiO 3 , and rare-earth nickelates), or not exploring pseudo-Jahn–Teller-like distortions in LaMnO 3 , or disallowing spontaneous formation of ordered vacancy compounds in Ba 4 As 3 and Ag 3 Al 22 O 34 ; and (iii) ignoring spin–orbit coupling forcing false metallic states in CaIrO 3 and Sr 2 IrO 4 . The distinction between false metals vs real insulators is important because (a) predicting theoretically that a given compound is metal even though it is found to be an insulator often creates the temptation to invoke high order novel physical effects (such as correlation in d-electron Mott insulators) to explain what was in effect caused by a more mundane artifact in a lower-level mean-field band theory, (b) recent prediction of exotic physical effects such as topological semimetals were unfortunately based on the above compounds that were misconstrued by theory to be metal, but are now recognized to be stable insulators not hosting exotic effects, and (c) practical technological applications based on stable degenerate but gapped metals such as transparent conductors or electrides for catalysis must rely on the systematically correct and reliable theoretical classification of metals vs insulators.

36 MATERIALS SCIENCE↗

A high efficiency and low vibration liquid nitrogen cooling system for silicon crystal based x-ray optics

A compact low-cost cryocooling system has been designed, constructed, and tested at SLAC National Accelerator Laboratory. The cooling power is provided by natural convection and phase change of the liquid nitrogen. The initial application was to cool silicon crystal optics to the sub-100 K range. A silicon crystal of dimension (width × depth × height) 50 × 50 × 30 mm3 has been used with an electric heater on the top surface in this prototyping test. This system can effectively provide more than 80 W of cooling power to the optics with a consumption of liquid nitrogen less than 2.1 l/h. The vibration of the silicon crystal was monitored during the tests with added electric heater power on the crystal. We note the vibration of the silicon crystal due to liquid nitrogen boiling is negligible.

47 OTHER INSTRUMENTATION↗

Marcus–Hush–Chidsey kinetics at electrode–electrolyte interfaces

Electrochemical kinetics at electrode–electrolyte interfaces limit the performance of devices including fuel cells and batteries. While the importance of moving beyond Butler–Volmer kinetics and incorporating the effect of electronic density of states of the electrode has been recognized, a unified framework that incorporates these aspects directly into electrochemical performance models is still lacking. In this work, we explicitly account for the density functional theory-calculated density of states numerically in calculating electrochemical reaction rates for a variety of electrode–electrolyte interfaces. We first show the utility of this for two cases related to Li metal electrodeposition and stripping on a Li surface and a Cu surface (anode-free configuration). The deviation in reaction rates is minor for cases with flat densities of states such as Li, but is significant for Cu due to nondispersive d-bands creating large variation. Finally, we consider a semiconducting case of a solid-electrolyte interphase consisting of LiF and Li2CO3 and note the importance of the Fermi level at the interface pinned by the redox reaction occurring there. We identify the asymmetry in reaction rates as a function of discharge/charge naturally within this approach.

Kurchin, Rachel (ORCID:0000000221474809)↗

Cleavable magnetic materials from van der Waals layered transition metal halides and chalcogenides

The study of exfoliated 2D magnetic materials is a vibrant and rapidly progressing field and is impacting many areas of condensed matter research including fundamental magnetism, spintronics and optoelectronics, and topological spin and electronic systems. The availability of promising candidate materials has enabled much of the progress in this field. Here, I present my perspective on the development of cleavable magnetic materials with a focus on transition metal halides and chalcogenides and an emphasis on crystal structure and magnetic order. To give proper context for these discussions, brief and incomplete surveys of recent work are included, focusing on specific aspects that I find most useful for guiding work on emerging materials and motivating expansion into other compounds and material families. Several instances of structural changes that can differentiate behaviors of bulk and ultrathin specimens are noted. Probing and understanding potential structural differences present a challenge but also an opportunity for material and device development, if they can be predicted and controlled. It is clear that careful investigations of structure, layer stacking, and defects in materials, how they may relate to the crystal chemistry, and how they may be different in crystals and mono- or few-layer specimens provide invaluable context for understanding the behavior of van der Waals layered magnetic materials.

36 MATERIALS SCIENCE↗

A fluid-structure interaction solver for investigating torsional galloping in solar-tracking photovoltaic panel arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. Furthermore, these simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects. These results are supplemented by a discussion of stability trends noted during testing and possible implications when considering multi-row array interactions.

14 SOLAR ENERGY↗

Synchrotron studies of functional interfaces and the state of the art: A perspective

The dramatic improvements in synchrotron light sources across the world imply tremendous opportunities for interface science. In this work, we describe a variety of current scattering and spectroscopic techniques with an eye toward how these will evolve, particularly with the advent of diffraction-limited sources. We also note the importance of in situ and operando methods for both understanding interface creation and interface evolution in different environments. As interfaces are often where critical processes occur, whether in regard to energy/chemical/materials/quantum information science, understanding the utility of these x-ray techniques will be of considerable importance to the cultivation of next-generation scientists.

36 MATERIALS SCIENCE↗

Machine learning methods for probabilistic locked-mode predictors in tokamak plasmas

A rotating tokamak plasma can interact resonantly with the external helical magnetic perturbations, also known as error fields. This can lead to locking and then to disruptions. We leverage machine learning (ML) methods to predict the locking events. We use a coupled third-order nonlinear ordinary differential equation model to represent the interaction of the magnetic perturbation and the plasma rotation with the error field. This model is sufficient to describe qualitatively the locking and unlocking bifurcations. Here, we explore using ML algorithms with the simulation data and experimental data, focusing on the methods that can be used with sparse datasets. These methods lead to the possibility of the avoidance of locking in real-time operations. We describe the operational space in terms of two control parameters: the magnitude of the error field and the rotation frequency associated with the momentum source that maintains the plasma rotation. The outcomes are quan- tified by order parameters that completely characterize the state, whether locked or unlocked. We use unsupervised ML methods to classify locked/unlocked states and note the usefulness of a certain normalization of the order parameters. Three supervised ML classifiers are used in suite to estimate the probability of locking in the region of control parameter space with hysteresis, i.e., the set of control parameters for which both locked and unlocked states can exist. The results show that a neural network gives the best estimate of the locking probability. An analogy of the present locking model with the van der Waals equation of state is also provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Molecular dynamics simulations of inelastic x-ray scattering from shocked copper

By taking the spatial and temporal Fourier transforms of the coordinates of the atoms in molecular dynamics simulations conducted using an embedded-atom-method potential, we calculate the inelastic scattering of x rays from copper single crystals shocked along [001] to pressures of up to 70 GPa. Above the Hugoniot elastic limit, we find that the copious stacking faults generated at the shock front introduce strong quasi-elastic scattering (QES) that competes with the inelastic scattering signal, which remains discernible within the first Brillouin zone; for specific directions in reciprocal space outside the first zone, the QES dominates the inelastic signal overwhelmingly. The synthetic scattering spectra we generate from our Fourier transforms suggest that energy resolutions of order 10 meV would be required to distinguish inelastic from quasi-elastic scattering within the first Brillouin zone of shock-loaded copper. We further note that high-resolution inelastic scattering also affords the possibility of directly measuring particle velocities via the Doppler shift. These simulations are of relevance to future planned inelastic scattering experiments at x-ray Free Electron Laser facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Characterization of thermal transport and evolution of Au plasma in ICF experiments by Thomson scattering

Here, this paper demonstrates the capability of optical Thomson scattering (OTS) to measure thermal transport, and provides support to radiation hydrodynamic and kinetic simulations of electron thermal transport and plasma evolution. OTS theory and plasma simulations are applied to the interpretation of experimental measurements of laser-produced plasma from spherical gold targets on the OMEGA laser facility. The dynamical form factor, S($\vec{k}$, ω), of electron density fluctuations that is used in the fitting of Thomson scattering spectra includes ion–ion collisions and effects of non-Maxwellian distribution functions. OTS measurements and their interpretation are consistent with the nonlocal transport model in radiation hydrodynamic simulations as well as with kinetic simulations in the second half of the probe pulse duration. In particular, the reversal of heat transport during cooling is observed to be consistent with simulations, while some discrepancies are noted during the initial heating of the Au targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Einstein–Rosen waves and the Geroch group

Under the action of the Geroch group, the Minkowski metric can be transformed into any vacuum metric with two commuting Killing vectors. In principle, this reduces the problem of deriving vacuum metrics with two commuting Killing vectors to pure algebra. In this short note, we use these facts to give a purely algebraic derivation of the Einstein-Rosen metric, which describes a cylindrical gravitational wave. Furthermore, our derivation has a straightforward extension to gravitational pulse waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Defect tolerance in halide perovskites: A first-principles perspective

In recent years, the impressive photovoltaic performance of halide perovskites has been commonly attributed to their defect tolerance. This attribution is seemingly intuitive and has been widely promoted in the field, though it has not been rigorously assessed. In this Perspective, we critically discuss the proposition of defect tolerance in halide perovskites based on first-principles calculations. We show that halide perovskites actually do suffer from defect-assisted nonradiative recombination, i.e., they are not defect tolerant. The nonradiative recombination rates in halide perovskites are comparable to or even greater than those in more conventional semiconductors. We note that to obtain accurate defect properties in halide perovskites, the level of theory and computational details are highly important, which was previously not sufficiently recognized. A distinctive feature of halide perovskites is that they can be grown with moderate defect densities using low-cost deposition techniques. But, similar to the case of conventional semiconductors, defect engineering is still key to improving the efficiency of perovskite solar cells.

14 SOLAR ENERGY↗

Progress toward fusion energy breakeven and gain as measured against the Lawson criterion

The Lawson criterion is a key concept in the pursuit of fusion energy, relating the fuel density n, pulse duration τ or energy confinement time τE, and fuel temperature T to the energy gain Q of a fusion plasma. The purpose of this paper is to explain and review the Lawson criterion and to provide a compilation of achieved parameters for a broad range of historical and contemporary fusion experiments. Although this paper focuses on the Lawson criterion, it is only one of many equally important factors in assessing the progress and ultimate likelihood of any fusion concept becoming a commercially viable fusion-energy system. Only experimentally measured or inferred values of n, τ or τE, and T that have been published in the peer-reviewed literature are included in this paper, unless noted otherwise. For extracting these parameters, we discuss methodologies that are necessarily specific to different fusion approaches (including magnetic, inertial, and magneto-inertial fusion). This paper is intended to serve as a reference for fusion researchers and a tutorial for all others interested in fusion energy.

Wurzel, Samuel E. (ORCID:0000000260548161)↗

Toward a machine-guided approach to energetic material discovery

In this article, we trained a machine learning (ML) model to connect microstructural details of an energetic material formulation to its performance for the purpose of guiding the discovery of new explosive formulations. Our hypothesis was that the algorithm would robustly learn the training data and produce an accurate surrogate model. Specifically, the algorithm learned the relationship between details of the void size distribution (VSD), initiating shock pressure, and the energetic material performance. We used realistic constraints on the VSD and a range of cases were ingested by a physically informed reactive flow model working within a hydrodynamic solver running on high-performance computing resources. The ML algorithm produced a surrogate model that accurately predicted known test points around the parameter space. In addition to the utility of the model and the process used for its development, we noted interesting comparisons between what we, the authors—subject matter experts, would heuristically conclude from the training data and the surrogate model predictions. We detected nuanced details that were missed by the surrogate model; however, these details are not important to an energetic material formulator. We concluded that the algorithm did indeed robustly learn the training data and produce an accurate surrogate model. We further concluded that the surrogate model is a powerful tool to guide the formulator in the absence of subject matter experts and limited-access computing resources.

42 ENGINEERING↗

Understanding the chemical bonding of ground and excited states of HfO and HfB with correlated wavefunction theory and density functional approximations

Knowledge of the chemical bonding of HfO and HfB ground and low-lying electronic states provides essential insights into a range of catalysts and materials that contain Hf–O or Hf–B moieties. Here, we carry out high-level multi-reference configuration interaction theory and coupled cluster quantum chemical calculations on these systems. We compute full potential energy curves, excitation energies, ionization energies, electronic configurations, and spectroscopic parameters with large quadruple-ζ and quintuple-ζ quality correlation consistent basis sets. We also investigate equilibrium chemical bonding patterns and effects of correlating core electrons on property predictions. Differences in the ground state electron configuration of HfB(X 4 Σ - ) and HfO(X 1 Σ + ) lead to a significantly stronger bond in HfO than HfB, as judged by both dissociation energies and equilibrium bond distances. We extend our analysis to the chemical bonding patterns of the isovalent HfX (X = O, S, Se, Te, and Po) series and observe similar trends. We also note a linear trend between the decreasing value of the dissociation energy (D e ) from HfO to HfPo and the singlet–triplet energy gap (ΔE S–T ) of the molecule. Finally, we compare these benchmark results to those obtained using density functional theory (DFT) with 23 exchange–correlation functionals spanning multiple rungs of “Jacob’s ladder.” When comparing DFT errors to coupled cluster reference values on dissociation energies, excitation energies, and ionization energies of HfB and HfO, we observe semi-local generalized gradient approximations to significantly outperform more complex and high-cost functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The role of breakup and coalescence in fine-scale bubble-induced turbulence. I. Dynamics

We study the effect of bubble breakup and coalescence on fine-scale dynamics of bubbly turbulent flows using direct numerical simulations. We perform two different simulations of dilute bubbly flows of void fraction 0.5%: one with bubbles breaking up and coalescing and the other without these physical processes. The volume of the fluid method is used for simulating bubbles undergoing breakup and coalescence while the bubbles are treated as rigid spheres in the immersed boundary method simulation. The energy spectrum in both types of simulation, consistent with previous studies, exhibits a -3 slope. We follow a single infinitesimal fluid element as it evolves to understand velocity gradient dynamics using conditional mean trajectories. We note finite-time divergence when the fluid element evolves under the action of inertial and pressure dynamics. The inertial, pressure, and viscous velocity gradient dynamics, when considered individually, produce the same results for bubble-induced turbulence (BIT) as with the classical homogeneous isotropic turbulence (HIT). Yet when the overall velocity gradient dynamics is considered, BIT results in non-cyclic trajectories moving toward stable node and unstable saddle while classical HIT shows cyclic behavior in their trajectories that move toward the origin. Interestingly, both the volume of fluid and immersed boundary simulations produce similar results. Therefore, there are two main takeaways from this research. First, new velocity gradient models are needed for BIT as their velocity gradients behavior is entirely different from the HIT. Second, we can neglect the bubble topology, breakup, and coalescence while studying or modeling the fine-scale dynamics of BIT.

42 ENGINEERING↗

Enforcing detailed balance in the Borgnakke–Larsen redistribution method with temperature dependent relaxation models

For decades, it has been observed that the commonly used Borgnakke–Larsen method for energy redistribution in Direct Simulation Monte Carlo codes fails to satisfy the principle of detailed balance when coupled to a wide variety of temperature dependent relaxation models, while seemingly satisfying detailed balance when coupled to others. Many attempts have been made to remedy the issue, yet much ambiguity remains, and no consensus appears in the literature regarding the root cause of the intermittent compatibility of the Borgnakke–Larsen method with temperature dependent relaxation models. This paper alleviates that ambiguity by presenting a rigorous theoretical derivation of the Borgnakke–Larsen method's requirement for satisfying detailed balance. Specifically, it is shown that the Borgnakke–Larsen method maintains detailed balance if and only if the probability of internal-energy exchange during a collision depends only on collision invariants (e.g., total energy). Here, the consequences of this result are explored in the context of several published definitions of relaxation temperature, including translational, total, and cell-averaged temperatures. Of particular note, it is shown that cell-averaged temperatures, which have been widely discussed in the literature as a way to ensure equilibrium is reached, also fail in a similar, although less dramatic, fashion when the aforementioned relationship is not enforced. The developed theory can be used when implementing existing or new relaxation models and will ensure that detailed balance is satisfied.

74 ATOMIC AND MOLECULAR PHYSICS↗

The role of breakup and coalescence in fine-scale bubble-induced turbulence. II. Kinematics

This second part of our research explores the kinematic aspect of fine-scale bubble-induced turbulence (BIT) to (i) present the effect of bubble breakup and coalescence and (ii) compare it against the universal kinematic fine-scale turbulence characteristics reported in the literature. To this end, we simulate a dilute bubbly system of 0.5% void fraction using two distinct numerical simulations. In the volume-of-fluid (VoF) simulation, bubbles undergo breakup and coalescence. In the immersed boundary method (IBM) simulation, however, they act as rigid spheres. We also perform a simulation of classical homogeneous isotropic turbulence (HIT). The first important outcome of this study is that BIT is radically different from HIT in terms of its kinematic fine-scale characteristics. In the vorticity-dominating regions, BIT exhibits a weak vortex stretching. This weak vortex stretching is due to (a) the intermediate strain-rate eigenvalues skewed weakly to positive and (b) the extensive strain-rate eigenvector aligning perpendicular to the vorticity vector. The BIT has, on average, not only a weak enstrophy production but also a weak strain production in strain-dominating regions. The weak strain production is due to (a) the presence of vortex stretching in highly strained fluid elements and (b) the absolute magnitude of compressive strain-rate eigenvalue being as close to the extensive strain-rate eigenvalue. Thus, none of the kinematic fine-scale HIT characteristics is noted for BIT. The second important conclusion is that bubble breakup and coalescence play little to no influence on the kinematics of fine-scale BIT as VoF and IBM simulations produce similar results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluation of a new DIII-D Doppler backscattering system for higher wavenumber measurement and signal enhancement

The high density fluctuation poloidal wavenumber, k θ (k θ > 8 cm –1 , k θ ρ s > 5, ρs is the ion gyro radius using the ion sound velocity), measurement capability of a new Doppler backscattering (DBS) system at the DIII-D tokamak has been experimentally evaluated. In DBS, wavenumber (k) matching becomes more important at higher wavenumbers, owing to the exponential dependence of the measured signal loss factor on wave vector mismatch. Wave vector matching allows for the Bragg scattering condition to be satisfied, which minimizes the signal loss at higher k’s. In the previous DBS system, without toroidal wave vector matching, the measured DBS signal-to-noise ratio at higher k θ (>8 cm –1 ) is substantially reduced, making it difficult to measure higher k θ turbulence. The new DBS system has been optimized to access higher wavenumber, k θ ≤ 20 cm –1 , density turbulence measurement. The optimization hardware addresses fluctuation wave vector matching using toroidal steering of the launch mirror to produce a backscattered signal with improved intensity. The probe’s sensitivity to high-k density fluctuations has been increased by approximately an order of magnitude compared to the old system that has been in use at DIII-D. Note that typical measurement locations are above or below the tokamak midplane on the low field side with normalized radial ranges of 0.5–1.0. As a result, the new DBS probe system with the toroidal matching of fluctuation wave vectors is thought to be critical to understanding high-k turbulent transport in fusion-relevant research at DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗