SEARCH · Engineering Papers

Results for “correspondence”

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 685 records · Page 38

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Chronostratigraphy of talus flatirons and piedmont alluvium along the Book Cliffs, Utah – Testing models of dryland escarpment evolution

Research on the evolution of dryland escarpments and what drives and controls their erosional retreat have been limited by a lack of chronologically constrained records. Prior investigations of escarpments in the southwestern U.S. and globally garner two endmember conceptual models -- one focused on bottom-up baselevel drivers and autogenic variations in piedmont and toeslope erosion processes, the other focused on the role of climate as a top-down forcing mechanism on sediment production and transport from cliffs. In this work, we test these conceptual models along a section of the Book Cliffs in central Utah, where climate has varied strongly and baselevel has fallen significantly over the late Quaternary. Four generations of talus flatirons and piedmont terraces are preserved in the study area we date them by optically stimulated luminescence and 10Be-exposure techniques. OSL depositional ages cluster into four generations of deposits with mean ages and mean standard errors of 117.8 ± 6.7 ka, 83.4 ± 3.3 ka, 52.4 ± 1.8 ka, and 5.6 ± 0.3 ka. 10Be-exposure results from talus boulders protruding from flatirons confirm that the youngest generation of talus is Holocene, with recent rock-fall boulders still being received on the landform top. Older flatiron generations produce incoherent, minimum exposure ages reflecting an increasing dominance of in situ boulder weathering and a steady-state denudation rate of ~45 mm/ky after tens of thousands of years. Based on the chronostratgraphy and correlation of talus and piedmont deposits, we interpret climate as the primary control on escarpment evolution. Depositional ages do not correspond to glacial epochs as historically presumed, but instead to episodes of intermediate and unstable or variable Pleistocene climate during Marine Isotope Stage-3 and 5 a/b. Results support the significance of high climate variability at millennial timescales in driving both sediment production by mass wasting from the cliff as well as sediment storage and reworking along the piedmont. In contrast, wetter climates of glacial epochs apparently caused greater erosion and sediment transport to trunk drainages. Our chronology also reveals a systematic decrease in deposit age with increasing distance up the piedmont from the regional baselevel of the Price River, consistent with transient fluvial incision and terrace abandonment. The Book Cliffs record notably does not correlate with similarly constrained records in Spain or the contrasting record from the Negev, yet they all share the feature of enhanced talus deposition and cliff retreat corresponding to climate instability rather than in response to steady glacial-climate conditions. Thus, we outline a model of escarpment evolution and retreat governed largely by top-down climate-forcing of sediment production and transport during climate variability, with local changes in baselevel and erosional processes on the piedmont being a modifying, secondary factor.

58 GEOSCIENCES↗

DOE OSTI2022

What lies beneath: Vertical temperature heterogeneity in a Mediterranean woodland savanna

As the availability of satellite and airborne thermal infrared remote sensing (TIR-RS) data increases and their spatial, temporal, and spectral resolutions improve, researchers are finding diverse applications for TIR-RS measurements. TIR-RS is now commonly applied in regional- and continental-scale analyses, such as those focused on fire and surface energy balance. However, its application lags in plant physiology and ecology, for which a finer-scale understanding of plant canopy temperatures would be useful to elucidate plant water dynamics, for example. In particular, while methods to disaggregate TIR-RS pixels in horizontal space have advanced, possible vertical stratification of plant canopy temperature and its implications for understanding the correspondence between TIR-RS and finer-scale, field-based thermal measurements (e.g. made with a thermal camera) remain unexplored. Here, we use data from a thermal camera deployed concurrently with the recent ECOSTRESS mission to quantify vertical temperature gradients within tree canopies and temperatures of over- vs. under-story plants in a Mediterranean woodland savanna. We then leverage diverse ancillary data to maximize the geometric comparability of ECOSTRESS and thermal camera measurements, in order to assess the extent to which the two forms of thermal measurements correspond. Specifically, we ask: (1) What are the patterns of intra-canopy and over- vs. under-story vertical temperature in a Mediterranean woodland savanna?, and (2) How can vertically-resolved, but spatially-limited field-based temperature measurements be reconciled with spatially-extensive, but surface-only, temperature measurements of a space-borne remote sensor? Further, we found consistent patterns of vertical thermal heterogeneity both within tree canopies and between ecosystem over- and under-stories. The daytime difference between the top and bottom thirds of blue oak canopies was, on average, 0.48 ° C – and sometimes several times larger. Notably, canopy tops are cooler, likely associated with the under-story grass reaching daytime temperatures often exceeding over-story temperatures by 10° C. Given the consistency of the intra-canopy temperature gradients, we expected the ECOSTRESS sensor would be in better agreement with camera measurements of canopy tops than bulk canopies or canopy bottoms. However, within-canopy gradients were overwhelmed by other sources of disagreement between the measurements, in part associated with upscaling camera measurements across space. Overall, thermal camera and ECOSTRESS measurements were largely in agreement at night (pixel RMSE = 1.1°C), but they were more divergent during times of low (but >0 W/m2) and high incoming solar radiation (daytime pixel RMSE = 3.5°C).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Influence of the air-cooled supercritical carbon dioxide recompression cycle design-point on molten salt power tower techno-economic performance

Here, this paper investigates the techno-economic performance of a molten salt power tower integrated with an air-cooled supercritical carbon dioxide recompression cycle. Prior research into the cycle design point showed trade-offs between efficiency, cost, and temperature difference across the heat input, resulting in trade-offs between the design-point generation of the molten salt power tower system and the cost of sensible-heat thermal energy storage. Prior research into cycle off-design showed that colder ambient temperature designs have higher efficiencies on cold days but lower efficiency and capacity degradation on hot days. We present a framework to evaluate how cycle design, cost, and corresponding off-design performance affect the techno-economic performance of the molten salt power tower system. We use a steady-state off-design cycle model assuming fixed shaft speeds, inventory control, and air-cooler fan power control and integrate the cycle model with the molten salt power tower model from the National Renewable Energy Laboratory's System Advisor Model. First, we present a surrogate model that uses a heuristic approach to correct the interpolation when the cycle reaches capacity limitations at ambient temperatures hotter than design. Next, we evaluate how the recompression cycle design-point ambient temperature, air-cooler approach temperature, air-cooler fan power, and recuperator conductance influence cycle and system design metrics. Finally, we complete annual simulations using an hourly "duck curve" pricing schedule and present annual performance metrics to evaluate the cumulative effects of the cycle design parameters and corresponding off-design performance on molten salt power tower techno-economics.

14 SOLAR ENERGY↗

DOE OSTI2020

Artificial linear brush abrasion of coatings for photovoltaic module first-surfaces

Natural soiling and the subsequent necessary cleaning of photovoltaic (PV) modules result in abrasion damage to the cover glass. The durability of the front glass has important economic consequences, including determining the use of anti-reflective and/or anti-soiling coatings as well as the method and frequency of operational maintenance (cleaning). Artificial linear brush abrasion using Nylon 6/12 bristles was therefore examined to explore the durability of representative PV first-surfaces, i.e., the surface of a module incident to direct solar radiation. Specimens examined include silane surface functionalized-, roughened (etched)-, porous silica-coated-, fluoropolymer-coated-, and ceramic (TiO 2 or ZrO 2 /SiO 2 /ZrO 2 /SiO 2 )-coated-glass, which are compared to monolithic-poly(methyl methacrylate) and -glass coupons. Characterization methods used in this study include: optical microscopy, ultraviolet–visible–near-infrared (UV-VIS-NIR) spectroscopy, sessile drop goniometry, white-light interferometry, atomic force microscopy (AFM), and depth-profiling X-ray photoelectron spectroscopy (XPS). The corresponding characteristics examined include: surface morphology, transmittance (i.e., optical performance), surface energy (water contact angle), surface roughness, scratch width and depth, and chemical composition, respectively. The study here was performed to determine coating failure modes; identify characterization methods that can detect nascent failures; compare the durability of popular contemporary coating materials; identify their corresponding damage characteristics; and compare slurry and dry-dust abrasion. This study will also aid in developing an abrasion standard for the PV industry.

14 SOLAR ENERGY↗

DOE OSTI2023

Topological band gap in intercalated epitaxial graphene

The study of functional manipulation of graphene is a critical subject, both for fundamental research and practical applications. In this study, we present that the intercalation of 5d transition metals into epitaxial graphene on SiC is a promising strategy for the realization of topologically nontrivial phases with a finite band gap in graphene. Employing first-principles calculations, grounded in density functional theory, we demonstrate that Re- and Ta-intercalated graphene evolve into two-dimensional topological insulators. These exhibit linear Dirac cones and quadratic bands with topological band gaps, respectively. The emergence of these topological states is attributed to the strong spin–orbit coupling strength of the intercalants. We show that the corresponding topological edge states persist within the finite bulk band gap, aligning with the bulk-boundary correspondence. Furthermore, we explore the spin splitting of the band structure, brought about by the inversion symmetry breaking and the spin–orbit coupling. Finally, our study underscores that the intercalation of graphene is an effective and a feasible approach for manipulating the band gap and the topological nature of graphene. Such intercalated graphene systems hold potential utility for spintronics and low-dimensional quantum device applications.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Implications of stop-and-go traffic on training learning-based car-following control

Learning-based car-following control (LCC) of connected and autonomous vehicles (CAVs) is gaining significant attention with the advancement of computing power and data accessibility. While the flexibility and large model capacity of model-free architecture enable LCC to potentially outperform the model-based car-following (CF) model in improving traffic efficiency and mitigating congestion, the generalizability of LCC for traffic conditions different from the training environment/dataset is not well-understood. Herein, this study seeks to explore the impact of stop-and-go traffic in the training dataset on the generalizability of LCC. It uses the characteristics of lead vehicle trajectories to describe stop-and-go traffic, and links the theory of identifiability (i.e., obtaining a unique parameter estimation result using sensor measurements) to the generalizability of behavior cloning (BC) and policy-based deep reinforcement learning (DRL). Correspondingly, the study shows theoretically that: (i) stop-and-go traffic can enable the property of identifiability and enhance the control performance of BC-based LCC in different traffic conditions; (ii) stop-and-go traffic is not necessary for DRL-based LCC to generalize to different traffic conditions; (iii) DRL-based LCC trained with only constant-speed lead vehicle trajectories (not sufficient to ensure identifiability) can be generalized to different traffic conditions; and (iv) stop-and-go traffic increases variance in the training dataset, which improves the convergence of parameter estimation while negatively impacting the convergence of DRL to the optimal control policy. Numerical experiments validate the above findings, illustrating that BC-based LCC entails comprehensive training datasets for generalizing to different traffic conditions, while DRL-based LCC can achieve generalization with simple free-flow traffic training environments. This further suggests DRL as a more promising and cost-effective LCC approach to reduce operational costs, mitigate traffic congestion, and enhance safety and mobility, which can accelerate the deployment and acceptance of CAVs.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2020

Evaluation of historical CMIP6 model simulations of extreme precipitation over contiguous US regions

Simulated historical precipitation is evaluated for Coupled Model Intercomparison Project Phase 6 (CMIP6) models using precipitation indices defined by the Expert Team on Climate Change Detection and Indices. The model indices are evaluated against corresponding indices from the CPC unified gauge-based analyses of precipitation over seven geographical regions across the contiguous US (CONUS). The regions assessed match those in recent US National Climate Assessment Reports. To estimate observational uncertainty, precipitation indices for three other observational datasets (HadEx2, Livneh and PRISM) are evaluated against the CPC analyses. Both the moderate and extreme mean precipitation intensities are overestimated over the western CONUS and underestimated in the areas of the Central Great Plains (CGP) in most CMIP6 models tested. Most CMIP6 models overestimate the mean and variability of wet spell durations and underestimate the mean and variability of dry spell durations across the CONUS. Biases in interannual variability of most of the indices have similar patterns to those in corresponding mean biases. The median and interquartile model spreads in CMIP6 model biases are clearly smaller than those in CMIP5 model biases for wet spell durations. Multimodel medians of CMIP6 (CMIP6-MMM) and CMIP5 (CMIP5-MMM) have similar biases in climatology and variability but biases tend to be smaller in CMIP6-MMM. Depending on the index, extreme precipitation is slightly better in parts of the eastern half of the CONUS in CMIP6-MMM, otherwise, the biases in climatology and variability are similar to CMIP5-MMM. CMIP6-MMM performs better than individual models and even observational datasets in some cases. Differences between observational datasets for most indices are comparable to the CMIP6 interquartile model spread. The better-performing observational and model datasets are different in different parts of the CONUS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Heterogeneity of neutrophils in arterial hypertension

Highlights: • Normal-density and low-density neutrophils from patients with hypertension influx more Na{sup +} compared to normotensive individuals. • Normal-density neutrophils from patients with hypertension generate low levels of reactive oxygen species. • Normal-density and low-density neutrophils from patients with hypertension exhibit lower levels of apoptosis. • Normal-density neutrophils from patients with hypertension maintain an activated and degranulated phenotype. • Low-density neutrophils from patients with hypertension show atypical mitochondrial morphology. Cellular heterogeneity and diversity are recognized to contribute to the functions of neutrophils under homeostatic and pathological conditions. We previously suggested that the chronic inflammatory responses associated with hypertension (HTN) are related to the participation of different subpopulations of neutrophils. Two populations of neutrophils can be obtained by density gradient centrifugation: normal-density neutrophils (NDN) and low-density neutrophils (LDN). However, the lack of standardized functional protocols has limited phenotypic characterization and functional comparisons of LDN and NDN. Based on their capability to incorporate Na{sup +}, maturity and activation stage, we characterized NDN and LDN in blood samples from ten patients with HTN and ten healthy individuals (HI) using flow cytometry. We compared the levels of reactive oxygen species (ROS), generation of neutrophil extracellular traps (NETs) and levels of apoptosis in NDN and LDN. In general, the NDN and LDN subpopulations from patients with HTN exhibited higher levels of sodium influx and ROS, and lower levels of apoptosis than the corresponding NDN and LDN subsets from HI. Transmission electron microscopy revealed NDN and LDN from patients with HTN exhibited alterations to mitochondrial morphology and fewer cytoplasmic granules than the corresponding HI subpopulations. Our results indicate both the NDN and LDN subpopulations enhance the effects of inflammation that contribute to the pathophysiology of HTN. Further detailed studies are required to characterize the events during ontogeny of the myeloid lineage that result in the diverse phenotypic characteristics of each subpopulation of LDN and NDN.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Molecular basis of functional exchangeability between ezrin and other actin-membrane associated proteins during cytokinesis

Highlights: • Anillin or supervillin depletion affects localization of different domains of ezrin. • Actin binding domains are compatible between ezrin and anillin for furrow ingression. • All ezrin and anillin domains are incompatible for polar cortical control. The mechanism that mediates the interaction between the contractile ring and the plasma membrane during cytokinesis remains elusive. We previously found that ERM (Ezrin/Radixin/Moesin) proteins, which usually mediate cellular pole contraction, become over-accumulated at the cell equator and support furrow ingression upon the loss of other actin-membrane associated proteins, anillin and supervillin. In this study, we addressed the molecular basis of the exchangeability between ezrin and other actin-membrane associated proteins in mediating cortical contraction during cytokinesis. We found that depletion of anillin and supervillin caused over-accumulation of the membrane-associated FERM domain and actin-binding C-terminal domain (C-term) of ezrin at the cleavage furrow, respectively. This finding suggests that ezrin differentially shares its binding sites with these proteins on the actin cytoskeleton or inner membrane surface. Using chimeric mutants, we found that ezrin C-term, but not the FERM domain, can substitute for the corresponding anillin domains in cytokinesis and cell proliferation. On the other hand, either the membrane-associated or the actin/myosin-binding domains of anillin could not substitute for the corresponding ezrin domains in controlling cortical blebbing at the cell poles. Our results highlight specific designs of actin- or membrane-associated moieties of different actin-membrane associated proteins with limited exchangeability, which enables them to support diverse cortical activities on the shared actin-membrane interface during cytokinesis.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Stability of expanding accretion shocks for an arbitrary equation of state

We present a theoretical stability analysis for an expanding accretion shock that does not involve a rarefaction wave behind it. The dispersion equation that determines the eigenvalues of the problem and the explicit formulae for the corresponding eigenfunction profiles are presented for an arbitrary equation of state and finite-strength shocks. For spherically and cylindrically expanding steady shock waves, we demonstrate the possibility of instability in a literal sense, a power-law growth of shock-front perturbations with time, in the range of $h_c< h<1+2 {\mathcal {M}}_2$ , where $h$ is the D'yakov-Kontorovich parameter, $h_c$ is its critical value corresponding to the onset of the instability and ${\mathcal {M}}_2$ is the downstream Mach number. Shock divergence is a stabilizing factor and, therefore, instability is found for high angular mode numbers. As the parameter $h$ increases from $h_c$ to $1+2 {\mathcal {M}}_2$ , the instability power index grows from zero to infinity. This result contrasts with the classic theory applicable to planar isolated shocks, which predicts spontaneous acoustic emission associated with constant-amplitude oscillations of the perturbed shock in the range $h_c< h<1+2 {\mathcal {M}}_2$ . Examples are given for three different equations of state: ideal gas, van der Waals gas and three-terms constitutive equation for simple metals.

Huete, César (ORCID:0000000232278520)↗

DOE OSTI2020

Theory of the tertiary instability and the Dimits shift within a scalar model

The Dimits shift is the shift between the threshold of the drift-wave primary instability and the actual onset of turbulent transport in a magnetized plasma. It is generally attributed to the suppression of turbulence by zonal flows, but developing a more detailed understanding calls for consideration of specific reduced models. The modified Terry–Horton system has been proposed by St-Onge ( J. Plasma Phys. , vol. 83, 2017, 905830504) as a minimal model capturing the Dimits shift. Here, we use this model to develop an analytic theory of the Dimits shift and a related theory of the tertiary instability of zonal flows. We show that tertiary modes are localized near extrema of the zonal velocity $U(x)$ , where $x$ is the radial coordinate. By approximating $U(x)$ with a parabola, we derive the tertiary-instability growth rate using two different methods and show that the tertiary instability is essentially the primary drift-wave instability modified by the local $U'' \doteq {\rm d}^2 U/{\rm d} x^2 $ . Then, depending on $U''$ , the tertiary instability can be suppressed or unleashed. The former corresponds to the case when zonal flows are strong enough to suppress turbulence (Dimits regime), while the latter corresponds to the case when zonal flows are unstable and turbulence develops. This understanding is different from the traditional paradigm that turbulence is controlled by the flow shear $| {\rm d} U / {\rm d} x |$ . Our analytic predictions are in agreement with direct numerical simulations of the modified Terry–Horton system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗

DOE OSTI2024

Nonlinear solution of classical three-wave interaction via finite-dimensional quantum model

The quantum three-wave interaction, the lowest-order nonlinear interaction in plasma physics, describes energy–momentum transfer between three resonant waves in the quantum regime. We describe how it may also act as a finite-degree-of-freedom approximation to the classical three-wave interaction in certain circumstances. By promoting the field variables to operators, we quantize the classical system, show that the quantum system has more free parameters than the classical system and explain how these parameters may be selected to optimize either initial or long-term correspondence. We then numerically compare the long-time quantum–classical correspondence far from the fixed point dynamics. We discuss the Poincaré recurrence of the system and the mitigation of quantum scrambling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Geometrical optics without singularities: using the ray time as the coordinate space

Geometrical optics (GO) is widely used for reduced modelling of waves in plasmas, but it fails near reflection points, where it predicts a spurious singularity of the wave amplitude. We show how to avoid this singularity by adopting a different representation of the wave equation. Instead of the physical coordinate 𝑥 and the wavevector 𝑘, we use the ray time 𝜏 as the new canonical coordinate and the ray energy ℎ as the associated canonical momentum. To derive the envelope equation in the 𝜏-representation, we construct the Weyl symbol calculus on the (𝜏,ℎ) space and show that the corresponding Weyl symbols are related to their (𝑥,𝑘) counterparts by the Airy transform. This allows us to express the coefficients in the envelope equation through the known properties of the original dispersion operator. When necessary, solutions of this equation can be mapped to the 𝑥-space using a generalised metaplectic transform. However, the field per se might not even be needed in practice. Instead, knowing the corresponding Wigner function usually suffices for linear and quasilinear calculations. As a Weyl symbol itself, the Wigner function can be mapped analytically, using the aforementioned Airy transform. We show that the standard Airy patterns that form in regions where conventional GO fails are successfully reproduced within metaplectic GO (MGO) simply by remapping the field from the 𝜏-space to the 𝑥-space. An extension to mode-converting waves is also presented. This formulation, which we call generalised MGO, can be particularly useful, for example, for reduced modelling of the O–X conversion in inhomogeneous plasma near the critical density, an effect that is important for fusion applications and also occurs in the ionosphere. Overall, MGO can replace GO for any practical purposes, because it better handles cutoffs and is similar otherwise.

plasma waves↗

DOE OSTI2020

Learning and meta-learning of stochastic advection–diffusion–reaction systems from sparse measurements

Physics-informed neural networks (PINNs) were recently proposed in [18] as an alternative way to solve partial differential equations (PDEs). A neural network (NN) represents the solution, while a PDE-induced NN is coupled to the solution NN, and all differential operators are treated using automatic differentiation. Here, we first employ the standard PINN and a stochastic version, sPINN, to solve forward and inverse problems governed by a non-linear advection–diffusion–reaction (ADR) equation, assuming we have some sparse measurements of the concentration field at random or pre-selected locations. Subsequently, we attempt to optimise the hyper-parameters of sPINN by using the Bayesian optimisation method (meta-learning) and compare the results with the empirically selected hyper-parameters of sPINN. In particular, for the first part in solving the inverse deterministic ADR, we assume that we only have a few high-fidelity measurements, whereas the rest of the data is of lower fidelity. Hence, the PINN is trained using a composite multi-fidelity network, first introduced in [12], that learns the correlations between the multi-fidelity data and predicts the unknown values of diffusivity, transport velocity and two reaction constants as well as the concentration field. For the stochastic ADR, we employ a Karhunen–Loève (KL) expansion to represent the stochastic diffusivity, and arbitrary polynomial chaos (aPC) to represent the stochastic solution. Correspondingly, we design multiple NNs to represent the mean of the solution and learn each aPC mode separately, whereas we employ a separate NN to represent the mean of diffusivity and another NN to learn all modes of the KL expansion. For the inverse problem, in addition to stochastic diffusivity and concentration fields, we also aim to obtain the (unknown) deterministic values of transport velocity and reaction constants. The available data correspond to 7spatial points for the diffusivity and 20 space–time points for the solution, both sampled 2000 times. We obtain good accuracy for the deterministic parameters of the order of 1–2% and excellent accuracy for the mean and variance of the stochastic fields, better than three digits of accuracy. In the second part, we consider the previous stochastic inverse problem, and we use Bayesian optimisation to find five hyper-parameters of sPINN, namely the width, depth and learning rate of two NNs for learning the modes. Here, we obtain much deeper and wider optimal NNs compared to the manual tuning, leading to even better accuracy, i.e., errors less than 1% for the deterministic values, and about an order of magnitude less for the stochastic fields.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

NASICON Na 3 V 2 (PO 4 ) 3 Enables Quasi-Two-Stage Na + and Zn 2+ Intercalation for Multivalent Zinc Batteries

Identifying positive electrode materials capable of reversible multivalent electrochemistry in electrolytes containing divalent ions such as Mg 2+ , Ca 2+ , and Zn 2+ at high operating potentials remains an ongoing challenge in “beyond lithium-ion” research. In this paper, we explore the Zn 2+ charge-storage mechanism of a vanadium-based Na + superionic conductor (NASICON), Na 3 V 2 (PO 4 ) 3 . By using X-ray synchrotron techniques to unravel potential-dependent structure–property relationships, we ascribe the reversible electrochemical behavior of Na 3 V 2 (PO 4 ) 3 to a quasi-two-stage intercalation process that involves both Na+ and Zn 2+ . Initial charging of Na 3 V 2 (PO 4 ) 3 leads to a Na + -extracted phase corresponding to Na 3 V 2 (PO 4 ) 3 , whereas subsequent discharge results predominantly in Na + intercalation followed by Zn 2+ intercalation. Operando X-ray diffraction of Na 3 V 2 (PO 4 ) 3 was used to study the phase changes associated with the first charge/discharge process, and ex situ measurements were used to precisely link the changes in the crystal structure to a quasi-two-stage intercalation of Na + and Zn 2+ . The corresponding changes in the V-oxidation state, V-O coordination, and the presence of Zn 2+ were confirmed by X-ray absorption spectroscopy. The results of this work present a comprehensive understanding of the charge-storage properties for a well-established NASICON structure that confers both the high capacity (~100 mA h g -1 ) and high potential (1.35 and 1.1 V vs Zn/Zn 2+ ).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Localized High Concentration Electrolytes for High Voltage Lithium-Metal Batteries: Correlation between the Electrolyte Composition and Its Reductive/Oxidative Stability

We demonstrate a first-principles screening methodology as an effective tool to explore electrolyte formulations for the new generation of high energy density rechargeable batteries. We study the liquid structure and electronic properties in dilute electrolytes, high concentration electrolytes (HCE), and localized high concentration electrolytes (LHCE), with focus on electrolyte formulations based on lithium bis(fluorosulfonyl)imide (LiFSI), dimethyl carbonate (DMC), and bis(2,2,2-trifluoroethyl) ether (BTFE) as a diluent. We describe the solvation complexes in the dilute electrolyte and explore structural changes triggered by the increase in lithium salt concentration for HCEs and the diluent effects in LHCEs. In HCE formulations, there is a 4-fold coordination environment of lithium-ions as in the dilute electrolyte, but the number of lithium-ion interactions with O atoms from FSI- anions dominates. In these solutions, the ability of the FSI- anions to interact with multiple lithium-ions allows complex 3D network formation and influences the reductive/oxidative behavior of the electrolyte. Interestingly, in LHCEs, the BTFE diluent molecules do not change the 3D solution structure when diluting the HCE formulation from 5.49 to 3.83 M. However, there is a composition threshold where the structural and electronic behavior may change. We show that diluting the HCE electrolyte with BTFE down to 1.77 M breaks the three-dimensional solution structure into an island-like solvation complex. We relate these structural changes to the electronic properties of the electrolytes finding a causal relationship between the reductive/oxidative behavior and the lithium-oxygen interaction mechanisms in the solvated complexes. The coordination with lithium-ions lowers the electrolyte LUMO and HOMO levels: the higher is the number of interactions with lithium-ions, the more likely the solvent molecule, FSI- anion, or diluent molecule is to be reduced and the less likely it is to become oxidized. The evolution of the solvated ion structure in HCE and LHCE suggests a close connection to a corresponding change in the lithium-ion transport mechanisms for these electrolytes.

anions, molecules, oxygen, solvents, electrolytes↗