Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational principles”

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 127 records · Page 7

Integrated machine learning-molecular dynamics framework for electrolyte property prediction

Electrochemical stability windows determine the operating range of battery electrolytes, yet accurate prediction remains challenging because stability emerges from statistical ensembles of local solvation environments rather than single ground-state molecular structures. Traditional density functional theory calculations on energy-minimized clusters cannot capture the thermal variations in local coordination environments and geometries that govern decomposition, while SMILES-based machine learning methods lack explicit representation of three-dimensional solvation structure and ion pairing. Here, we introduce a structure-aware machine learning framework that predicts frontier orbital energies (HOMO and LUMO) directly from molecular dynamics-sampled solvation configurations, achieving sub-0.6 eV accuracy at computational costs 3–4 orders of magnitude lower than first-principles methods. Across twelve representative battery electrolytes, we demonstrate that solvent-separated and contact ion pairs exhibit strong size- and local chemistry dependent electronic stability, with variations in coordination shifts of HOMO or LUMO level by 2–3 eV, and that extended solvation structure and partially desolvated environment further modulate stability by up to 3 eV. By encoding the statistical nature of electrochemical failure through ensemble sampling of explicit solvation geometries, our approach enables high-throughput screening and rational design of next-generation battery electrolytes with mechanistic understanding of structure–property relationships.

Energy - Storage↗

Modeling Short-Range Order in Disordered Rocksalt Cathodes by Pair Distribution Function Analysis

Pair distribution function (PDF) analysis is a powerful technique for the characterization of short-range order (SRO) in disordered materials. Accurate interpretation of experimental PDF data is critically reliant on the development of structural models that can account for local variations in site occupancies and bond lengths. To this end, we outline an approach to model SRO using first-principles calculations based on the cluster-expansion formalism. These methods are validated on neutron scattering data from two disordered rocksalt oxyfluorides, Li 1.3 Mn 0.4 Ti 0.3 O 1.7 F 0.3 and Li 1.3 Mn 0.4 Nb 0.2 Ti 0.1 O 1.7 F 0.3 . For each composition, we demonstrate that an average structure without any SRO fails to reproduce several key features in the experimental PDF. To pinpoint the origin of the suspected SRO in these materials, configurational and displacive effects were separately investigated using two disparate models. Special quasi-random structures were relaxed using density functional theory to account for local changes in bond lengths while maintaining a near-random ionic configuration. This leads to slightly improved accuracy but still misrepresents asymmetry in the first few peaks of the PDF. Monte Carlo simulations were performed to model configurational SRO on a fixed lattice, which by itself is shown to have a minimal influence on the PDF. Instead, we find that it is the bond length relaxations within environments created by SRO which controls the details of the PDF, thereby highlighting the subtle but important coupling between configurational and displacive SRO in disordered materials.

36 MATERIALS SCIENCE↗

Proximal Galerkin: A Structure-Preserving Finite Element Method for Pointwise Bound Constraints

The proximal Galerkin finite element method is a high-order, low iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of pointwise bound constraints in infinite-dimensional function spaces. This paper introduces the proximal Galerkin method and applies it to solve free boundary problems, enforce discrete maximum principles, and develop a scalable, mesh-independent algorithm for optimal design with pointwise bound constraints. This paper also introduces the latent variable proximal point (LVPP) algorithm, from which the proximal Galerkin method derives. When analyzing the classical obstacle problem, we discover that the underlying variational inequality can be replaced by a sequence of second-order partial differential equations (PDEs) that are readily discretized and solved with, e.g., the proximal Galerkin method. Throughout this work, we arrive at several contributions that may be of independent interest. These include (1) a semilinear PDE we refer to as the entropic Poisson equation; (2) an algebraic/geometric connection between high-order positivity-preserving discretizations and certain infinite-dimensional Lie groups; and (3) a gradient-based, bound-preserving algorithm for two-field, density-based topology optimization. The complete proximal Galerkin methodology combines ideas from nonlinear programming, functional analysis, tropical algebra, and differential geometry and can potentially lead to new synergies among these areas as well as within variational and numerical analysis. Open-source implementations of our methods accompany this work to facilitate reproduction and broader adoption.

97 MATHEMATICS AND COMPUTING↗

Toward extracting $$\gamma $$ from $$B\rightarrow DK$$ without binning

Abstract $$B^\pm \rightarrow DK^\pm $$ B ± → D K ± transitions are known to provide theoretically clean information about the CKM angle $$\gamma $$ γ , with the most precise available methods exploiting the cascade decay of the neutral D into CP self-conjugate states. Such analyses currently require binning in the D decay Dalitz plot, while a recently proposed method replaces this binning with the truncation of a Fourier series expansion. In this paper, we present a proof of principle of a novel alternative to these two methods, in which no approximations at the level of the data representation are required. In particular, our new strategy makes no assumptions about the amplitude and strong phase variation over the Dalitz plot. This comes at the cost of a degree of ambiguity in the choice of test statistic quantifying the compatibility of the data with a given value of $$\gamma $$ γ , with improved choices of test statistic yielding higher sensitivity. While our current proof-of-principle implementation does not demonstrate optimal sensitivity to $$\gamma $$ γ , its conceptually novel approach opens the door to new strategies for $$\gamma $$ γ extraction. More studies are required to see if these can be competitive with the existing methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Long-read sequencing transcriptome quantification with lr-kallisto

RNA abundance quantification has become routine and affordable thanks to high-throughput “short-read” technologies that provide accurate molecule counts at the gene level. Similarly accurate and affordable quantification of definitive full-length, transcript isoforms has remained a stubborn challenge, despite its obvious biological significance across a wide range of problems. “Long-read” sequencing platforms now produce data-types that can, in principle, drive routine definitive isoform quantification. However some particulars of contemporary long-read datatypes, together with isoform complexity and genetic variation, present bioinformatic challenges. We show here, using ONT data, that fast and accurate quantification of long-read data is possible and that it is improved by exome capture. To perform quantifications we developed lr-kallisto, which adapts the kallisto bulk and single-cell RNA-seq quantification methods for long-read technologies.

Loving, Rebekah K. (ORCID:0000000187250376)↗

Compositional Variation Tests on DuraMelter 100 with LAW Sub-Envelope C2 Feed in Support of the LAW Pilot Melter (Final Report)

One of the principal objectives of the DM100 melter tests was to demonstrate that the LAW Sub-Envelope C2 feeds with ± 15% variations in the amounts of waste simulant are suitable for large-scale tests on the LAW Pilot melter. The other principle objective for the DM100 tests was to determine the relative effects of sugar versus waste organics on glass redox state, NOx emissions, sulfur emissions, sulfur retention in the glass, secondary sulfate layer formation, and glass-melt foaming. Additional objectives included measurement of the feed rheological properties and collection of melter emissions data. Screening tests are conducted at the DM100 scale in order to identify and correct potential problems before testing at the larger scale LAW Pilot melter. Testing goals were met in that the results of the DM100 melter test resulted in selection of LAW Sub-Envelope C2 feed suitable for Pilot Melter testing.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Constructing nested coordinates inside strongly shaped toroids using an action principle

A new approach for constructing polar-like boundary-conforming coordinates inside a toroid with strongly shaped cross-sections is presented. A coordinate mapping is obtained through a variational approach, which involves identifying extremal points of a proposed action in the mapping space from [0,2π] 2 ×[0,1] to a toroidal domain in $\mathbb{R}$ 3 . This approach employs an action built on the squared Jacobian and radial length. Extensive testing is conducted on general toroidal boundaries using a global Fourier–Zernike basis via action minimisation. The results demonstrate successful coordinate construction capable of accurately describing strongly shaped toroidal domains. The coordinate construction is successfully applied to the computation of three-dimensional magnetohydrodynamic equilibria in the GVEC code where the use of traditional coordinate construction by interpolation from the boundary failed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

AtomAI framework for deep learning analysis of image and spectroscopy data in electron and scanning probe microscopy

Over the past several decades, electron and scanning probe microscopes have become critical components of condensed matter physics, materials science and chemistry research. At the same time, the infrastructure for establishing a connection between microscopy observations and materials behaviour over a broader parameter space is lacking. In this work, we introduce AtomAI, an open-source software package bridging instrument-specific Python libraries, deep learning and simulation tools into a single ecosystem. AtomAI allows direct applications of deep neural networks for atomic and mesoscopic image segmentation converting image and spectroscopy data into class-based local descriptors for downstream tasks such as statistical and graph analysis. For atomically resolved imaging data, the output is types and positions of atomic species, with an option for subsequent refinement. AtomAI further allows the implementation of a broad range of image and spectrum analysis functions, including invariant variational autoencoders for disentangling structural factors of variation and im2spec type of encoder–decoder models for mapping structure–property relationships. Finally, our framework allows seamless connection to the first principles modelling with a Python interface on the inferred atomic positions.

36 MATERIALS SCIENCE↗

Mapping the three-dimensional fermiology of the triangular lattice magnet EuAg 4 Sb 2

In this paper, we report the temperature-field phase diagram as well as present a comprehensive study of the electronic structure and three-dimensional fermiology of the triangular-lattice magnet EuAg 4 ⁢Sb 2 , utilizing quantum oscillation measurements, angle-resolved photoemission spectroscopy, and first-principles calculations. The complex magnetic phase diagram of EuAg 4 ⁢Sb 2 highlights many transitions through nontrivial AFM states. Shubnikov-de Haas and de Haas-van Alphen oscillations were observed in the polarized ferromagnetic state of EuAg 4 ⁢Sb 2 , revealing three pairs of distinct spin-split frequency branches with small effective masses. A comparison of the angle-dependent oscillation data with first-principles calculations in the ferromagnetic state and angle-resolved photoemission spectra shows good agreement, identifying tubular hole pockets and hourglass-shaped hole pockets at the Brillouin zone center, as well as diamond-shaped electron pockets at the zone boundary. As the temperature increases, the frequency branches of the tiny hourglass pockets evolve into a more cylindrical shape, while the larger pockets remain unchanged. This highlights that variations in exchange splitting, driven by changes in the magnetic moment, primarily impact the small Fermi pockets without significantly altering the overall band structure. As a result, this is consistent with first-principles calculations, which show minimal changes near the Fermi level across ferromagnetic and simple antiferromagnetic states or under varying on-site Coulomb repulsion.

36 MATERIALS SCIENCE↗

Biocrust Mosses and Cyanobacteria Exhibit Distinct Carbon Uptake Responses to Variations in Precipitation Amount and Frequency

ABSTRACT Dryland organisms exhibit varied responses to changes in precipitation, including event size, frequency, and soil moisture duration, influencing carbon uptake and reserve management strategies. This principle, central to the pulse‐reserve paradigm, has not been thoroughly evaluated in biological soil crusts (biocrusts), essential primary producers on dryland surfaces. We conducted two experiments to investigate carbon uptake in biocrusts under different precipitation regimes. In the first, we applied a gradient of watering amounts to biocrusts dominated by moss or cyanobacteria, hypothesising distinct pulse‐response strategies. The second experiment extended watering treatments over three months, varying pulse size and frequency. Our results revealed distinct carbon uptake patterns: moss crusts exhibited increased CO 2 uptake with larger, less frequent watering events, whereas cyanobacteria crusts maintained similar carbon uptake across all event sizes. These findings suggest divergent pulse‐response strategies across biocrust types, with implications for modelling dryland carbon dynamics and informing land management under changing precipitation regimes.

Young, Kristina E. [USDA‐ARS Jornada Experimental ↗

Leveraging neural control variates for enhanced precision in lattice field theory

Results obtained with stochastic methods have an inherent uncertainty due to the finite number of samples that can be achieved in practice. In lattice QCD this problem is particularly salient in some observables like, for instance, observables involving one or more baryons and it is the main problem preventing the calculation of nuclear forces from first principles. The method of control variables has been used extensively in statistics and it amounts to computing the expectation value of the difference between the observable of interest and another observable whose average is known to be zero but is correlated with the observable of interest. Recently, control variates methods emerged as a promising solution in the context of lattice field theories. In our current study, instead of relying on an educated guess to determine the control variate, we utilize a neural network to parametrize this function. Using 1 + 1 dimensional scalar field theory as a testbed, we demonstrate that this neural network approach yields substantial improvements. Notably, our findings indicate that the neural network ansatz is particularly effective in the strong coupling regime. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Nanoscale lattice dynamics in hexagonal boron nitride moiré superlattices

Twisted two-dimensional van der Waals (vdW) heterostructures have unlocked a new means for manipulating the properties of quantum materials. The resulting mesoscopic moiré superlattices are accessible to a wide variety of scanning probes. To date, spatially-resolved techniques have prioritized electronic structure visualization, with lattice response experiments only in their infancy. Here, we therefore investigate lattice dynamics in twisted layers of hexagonal boron nitride (hBN), formed by a minute twist angle between two hBN monolayers assembled on a graphite substrate. Nano-infrared (nano-IR) spectroscopy reveals systematic variations of the in-plane optical phonon frequencies amongst the triangular domains and domain walls in the hBN moiré superlattices. Our first-principles calculations unveil a local and stacking-dependent interaction with the underlying graphite, prompting symmetry-breaking between the otherwise identical neighboring moiré domains of twisted hBN.

36 MATERIALS SCIENCE↗

Editorial: Linking Hydrological and Biogeochemical Processes in Riparian Corridors

The riparian corridor is a key component of the critical zone and an essential component of watershed systems. According to Merriam-Webster Dictionary, the word riparian is derived from the Latin word riparius, meaning “existing alongside a river.” Riparian corridors typically extend from a few meters to hundreds of meters adjacent to a river and are marked by rich biodiversity, vegetation, and intense biogeochemical activity. They act as integrators of watershed processes and constitute the primary pathways for the subsurface geochemical exports from the watershed. Although riparian corridors comprise only 2–10% of a watershed's area, as much as 90–98% of biogeochemical processing in watersheds occurs in this region, thereby affecting the subsurface geochemical exports and downstream river water quality. Indeed, the riparian corridor is a good example of the Pareto principle. This outsized contribution occurs at the interface between aquatic (river) and terrestrial (land) environments, where interactions between hydrologic and biogeochemical processes are intensified. Variations in the river corridor over time can thus also have outsize impacts. Therefore, it is important to understand the hydrological and biogeochemical linkages in riparian corridors to determine water availability and quality for sustainable management.

54 ENVIRONMENTAL SCIENCES↗

Development of a health monitoring framework: Application to a supercritical pulverized coal-fired boiler

In this study, this work details the development of a physics-based equipment health monitoring framework for a supercritical boiler, using first-principles models to estimate the remaining useful life (RUL) of its components. The framework accounts for fatigue and creep life consumption, generating spatio-temporal variations in mechanical and thermal stress. Analysis of the stress profile throughout the boiler highlights the finishing superheater inlet steam header as a vulnerable location susceptible to damage from cycling operation. The framework also yields quantified uncertainty in the RUL projection for specific locations, accounting for uncertainties in material properties and boiler operation. Results indicate that operational uncertainties (e.g., seasonal variation and operational strategy) and material properties (e.g., rupture time coefficients, Young’s modulus, yield strength, and coefficient of thermal expansion) significantly impact the RUL of the finishing superheater inlet steam header. Additionally, case studies demonstrate the use of the health monitoring framework as a predictive tool for operational planning under uncertainty, including scenarios with and without updates on the operation of the boiler.

20 FOSSIL-FUELED POWER PLANTS↗

Ecological theory applied to environmental metabolomes reveals compositional divergence despite conserved molecular properties

Stream and river systems transport and process substantial amounts of dissolved organic matter (DOM) from terrestrial and aquatic sources to the ocean, with global biogeochemical implications. However, the underlying mechanisms affecting the spatiotemporal organization of DOM composition are under-investigated. To understand the principles governing DOM composition, we leverage the recently proposed synthesis of metacommunity ecology and metabolomics, termed ‘meta-metabolome ecology.’ Applying this novel approach to a freshwater ecosystem, we demonstrated that despite similar molecular properties across metabolomes, metabolite identity significantly diverged due to environmental filtering and variations in putative biochemical transformations. We refer to this phenomenon as ‘thermodynamic redundancy,’ which is analogous to the ecological concept of functional redundancy. We suggest that under thermodynamic redundancy, divergent metabolomes can support equivalent biogeochemical function just as divergent ecological communities can support equivalent ecosystem function. As these analyses are performed in additional ecosystems, potentially generalizable principles, like thermodynamic redundancy, can be revealed and provide insight into DOM dynamics.

Danczak, Robert E.↗

Lagrangian formulation of nuclear–electronic orbital Ehrenfest dynamics with real-time TDDFT for extended periodic systems

Here, we present a Lagrangian-based implementation of Ehrenfest dynamics with nuclear–electronic orbital (NEO) theory and real-time time-dependent density functional theory for extended periodic systems. In addition to a quantum dynamical treatment of electrons and selected protons, this approach allows for the classical movement of all other nuclei to be taken into account in simulations of condensed matter systems. Furthermore, we introduce a Lagrangian formulation for the traveling proton basis approach and propose new schemes to enhance its application for extended periodic systems. Validation and proof-of-principle applications are performed on electronically excited proton transfer in the o-hydroxybenzaldehyde molecule with explicit solvating water molecules. These simulations demonstrate the importance of solvation dynamics and a quantum treatment of transferring protons. This work broadens the applicability of the NEO Ehrenfest dynamics approach for studying complex heterogeneous systems in the condensed phase.

Calculus of variations↗

Organic Electrochemical Transistor Channel Materials: Copolymerization Versus Physical Mixing of Glycolated and Alkoxylated Polymers

Organic electrochemical transistors (OECTs) feature a polymer channel capable of conducting both ions and electronic charges. The choice of the channel material is critical for OECT performance. Many efforts have focused on improving performance via the chemical tunability of conjugated polymers – through backbone, side chain, and molar mass engineering – leading to useful design principles for accumulation-mode OECT materials. However, tuning the chemical structure of conjugated polymers often requires time-consuming optimization of the synthesis route. Meanwhile, variations in molar mass, dispersity, structural defects, and metal content present challenges when attempting to analyze the detailed effects of structural modifications, as multiple performance-determining factors are often (unintentionally) changed at the same time. Therefore, this study explores blended channel materials obtained by physically mixing glycolated and alkoxylated polymers in different ratios, and compares their OECT performance with the corresponding statistical copolymers. It is shown that mixing two well-performing materials creates blends that enable rational tuning of the transistor properties without compromising on performance. Thus, channels based on blends of alkoxylated and glycolated polymers hold promise for OECT technology with tailored response, as only two materials are needed to achieve any desired side chain ratio, simplifying the optimization of OECT characteristics.

copolymerization↗

Variational Self-Consistent Theory for Beam-Loaded Cavities

A new variational theory is presented for beam loading in microwave cavities. The beam{field interaction is formulated as a dynamical interaction whose stationarity according to Hamilton’s principle will naturally lead to steady-state solutions that indicate how a cavity’s resonant frequency, Q and optimal coupling coefficient will detune as a result of the beam loading. A driven cavity Lagrangian is derived from first principles, including the effects of cavity wall losses, input power and beam interaction. The general formulation is applied to a typical klystron input cavity to predict the appropriate detuning parameters required to maximize the gain (or modulation depth) in the average Lorentz factor boost, (Δ γ). Numerical examples are presented, showing agreement with the general detuning trends previously observed in the literature. The developed formulation carries several advantages for beam-loaded cavity structures. It provides a self-consistent model for the dynamical (nonlinear) beam{field interaction, a procedure for maximizing gain under beam-loading conditions, and a useful set of parameters to guide cavity-shape optimization during the design of beam-loaded systems. Enhanced clarity of the physical picture underlying the problem seems to be gained using this approach, allowing straightforward inclusion or exclusion of different field configurations in the calculation and expressing the final results in terms of measurable quantities. Two field configurations are discussed for the klystron input cavity, using finite magnetic confinement or no confinement at all. Formulating the problem in a language that is directly accessible to the powerful techniques found in Hamiltonian dynamics and canonical transformations may potentially carry an additional advantage in terms of analytical computational gains, under suitable conditions

43 PARTICLE ACCELERATORS↗