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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Molecular interactions from the density functional theory for chemical reactivity: Interaction chemical potential, hardness, and reactivity principles

In the first paper of this series, the authors derived an expression for the interaction energy between two reagents in terms of the chemical reactivity indicators that can be derived from density functional perturbation theory. While negative interaction energies can explain reactivity, reactivity is often more simply explained using the “|dμ| big is good” rule or the maximum hardness principle. Expressions for the change in chemical potential (μ) and hardness when two reagents interact are derived. A partial justification for the maximum hardness principle is that the terms that appear in the interaction energy expression often reappear in the expression for the interaction hardness, but with opposite sign.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Larger than uncorrelated vacancy diffusion contributions in chemically disordered crystalline materials

Arising from a variational approach to compute diffusion coefficients, we introduce “superkinetic kinosons”, which are contributions to the diffusion flux of mobile defects and atomic species by single jump mechanisms that can exceed their uncorrelated values. In vacancy-mediated diffusion in crystalline random and high entropy materials, these contributions primarily arise from intermittent jumps of slow moving atoms that remove the vacancy out of correlation traps and enable a quantification of such effects on the diffusion coefficients. They can be significant contributors to diffusion, even when faster moving atoms have infinite percolation networks. Furthermore, their existence and importance provide fundamental insights about underlying aspects of vacancy diffusion such as formation of localized correlation traps and deviation of vacancy diffusion from percolation behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2021

SARS-CoV-2 RBD antibodies that maximize breadth and resistance to escape

An ideal therapeutic anti-SARS-CoV-2 antibody would resist viral escape 1,2,3 , have activity against diverse sarbecoviruses 4,5,6,7 , and be highly protective through viral neutralization 8,9,10,11 and effector functions 12,13 . Understanding how these properties relate to each other and vary across epitopes would aid the development of therapeutic antibodies and guide vaccine design. Here we comprehensively characterize escape, breadth and potency across a panel of SARS-CoV-2 antibodies targeting the receptor-binding domain (RBD). Despite a trade-off between in vitro neutralization potency and breadth of sarbecovirus binding, we identify neutralizing antibodies with exceptional sarbecovirus breadth and a corresponding resistance to SARS-CoV-2 escape. One of these antibodies, S2H97, binds with high affinity across all sarbecovirus clades to a cryptic epitope and prophylactically protects hamsters from viral challenge. Antibodies that target the angiotensin-converting enzyme 2 (ACE2) receptor-binding motif (RBM) typically have poor breadth and are readily escaped by mutations despite high neutralization potency. Nevertheless, we also characterize a potent RBM antibody (S2E12 8 ) with breadth across sarbecoviruses related to SARS-CoV-2 and a high barrier to viral escape. These data highlight principles underlying variation in escape, breadth and potency among antibodies that target the RBD, and identify epitopes and features to prioritize for therapeutic development against the current and potential future pandemics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Model reduction methods for nuclear emulators

The field of model order reduction (MOR) is growing in importance due to its ability to extract the key insights from complex simulations while discarding computationally burdensome and superfluous information. We provide an overview of MOR methods for the creation of fast & accurate emulators of memory- and compute-intensive nuclear systems, focusing on eigen-emulators and variational emulators. As an example, we describe how 'eigenvector continuation' is a special case of a much more general and well-studied MOR formalism for parameterized systems. We continue with an introduction to the Ritz and Galerkin projection methods that underpin many such emulators, while pointing to the relevant MOR theory and its successful applications along the way. Here, we believe that this guide will open the door to broader applications in nuclear physics and facilitate communication with practitioners in other fields.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Thermodynamically consistent semi-compressible fluids: a variational perspective

This paper presents (Lagrangian) variational formulations for single and multicomponent semi-compressible fluids with both reversible (entropy-conserving) and irreversible (entropy-generating) processes. Semi-compressible fluids are useful in describing low-Mach dynamics, since they are soundproof. These models find wide use in many areas of fluid dynamics, including both geophysical and astrophysical fluid dynamics. Specifically, the Boussinesq, anelastic and pseudoincompressible equations are developed through a unified treatment valid for arbitrary Riemannian manifolds, thermodynamic potentials and geopotentials. By design, these formulations obey the 1st and 2nd laws of thermodynamics, ensuring their thermodynamic consistency. This general approach extends and unifies existing work, and helps clarify the thermodynamics of semi-compressible fluids. To further this goal, evolution equations are presented for a wide range of thermodynamic variables: entropy density s, specific entropy η, buoyancy b, temperature T, potential temperature θ and a generic entropic variable χ; along with a general definition of buoyancy valid for all three semicompressible models and arbitrary geopotentials. Finally, the elliptic equation for the pressure perturbation (the Lagrange multiplier that enforces semi-compressibility) is developed for all three equation sets in the case of reversible dynamics, and for the Boussinesq/anelastic equations in the case of irreversible dynamics; and some discussion is given of the difficulty in formulating it for the pseudoincompressible equations with irreversible dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

BUQEYE guide to projection-based emulators in nuclear physics

The BUQEYE collaboration (Bayesian Uncertainty Quantification: Errors in Your effective field theory) presents a pedagogical introduction to projection-based, reduced-order emulators for applications in low-energy nuclear physics. The term emulator refers here to a fast surrogate model capable of reliably approximating high-fidelity models. As the general tools employed by these emulators are not yet well-known in the nuclear physics community, we discuss variational and Galerkin projection methods, emphasize the benefits of offline-online decompositions, and explore how these concepts lead to emulators for bound and scattering systems that enable fast and accurate calculations using many different model parameter sets. We also point to future extensions and applications of these emulators for nuclear physics, guided by the mature field of model (order) reduction. All examples discussed here and more are available as interactive, open-source Python code so that practitioners can readily adapt projection-based emulators for their own work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Tailoring the Selective Oxidation of Hydroxyl-Containing Compounds via Precisely Tuning the Hydrogen-Bond Strength of Catalyst H-Bond Acceptors

The unique performance of the enzyme is mainly achieved via weak interactions between the “outer coordination sphere” and the substrate. Inspired by this process, we developed 3D encapsulated-structure catalysts with hydrogen-bond engineering on the shell, which mimics the “outer coordination sphere” of an enzyme. Various hydrogen bond acceptors (C=O, S=O, and N–O groups) are imparted in the shell. Concentration-dependent 1H NMR, inverse-phase gas Chromatography (IGC) measurements, and DFT calculations underscore that the hydrogen bond strength between the acceptor groups and alcohol follows the order of C=O < S=O < N–O. The hydroxyl compound oxidation rate vs the hydrogen bond strength follows a volcano behavior, reminiscent of Sabatier’s principle. The performance variation among catalysts is attributed to the adsorption strength of the substrate. The proposed bioinspired design principle expands the scope of encapsulated catalysts, enabling fine regulation of catalytic activity through precise microenvironment control via weak interactions with substrates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Tetraquarks made of sufficiently unequal-mass heavy quarks are bound in QCD

Tetraquarks, bound states composed of two quarks and two antiquarks, have been the subject of intense study but are challenging to understand from first principles. We apply variational and Green’s function Monte Carlo methods to compute tetraquark ground-state energies in potential nonrelativistic QCD using a wide range of color and spatial wave functions. We find no evidence for bound tetraquarks composed of equal-mass quarks and antiquarks. Conversely, we find clear evidence for the existence of bound tetraquarks for sufficiently unequal quark/antiquark mass ratios at all overall mass scales where our effective theory results are applicable. We predict the critical mass ratios for bound state formation and study tetraquark bound states’ spatial and color structure at leading order and next-to-leading order in potential nonrelativistic QCD. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Stability analysis of the Eulerian–Lagrangian finite volume methods for nonlinear hyperbolic equations in one space dimension

In this paper, we construct a novel Eulerian–Lagrangian finite volume (ELFV) method for nonlinear scalar hyperbolic equations in one space dimension. It is well known that the exact solutions to such problems may contain shocks though the initial conditions are smooth, and direct numerical methods may suffer from restricted time step sizes. To relieve the restriction, we propose an ELFV method, where the space-time domain was separated by the partition lines originated from the cell interfaces whose slopes are obtained following the Rakine–Hugoniot junmp condition. Unfortunately, to avoid the intersection of the partition lines, the time step sizes are still limited. To fix this gap, we detect effective troubled cells (ETCs) and carefully design the influence region of each ETC, within which the partitioned space-time regions are merged together to form a new one. Then with the new partition of the space-time domain, we theoretically prove that the proposed first-order scheme with Euler forward time discretization is total-variation-diminishing and maximum-principle-preserving with at least twice larger time step constraints than the classical first order Eulerian method for Burgers’ equation. Numerical experiments verify the optimality of the designed time step sizes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Metal-to-insulator transition in oxide semimetals by anion doping

Dirac semimetals exhibiting nontrivial topological characteristics stand as exemplary parent compounds, offering great promise for the realization of novel electronic states. These materials also serve as a versatile platform for systematically investigating topological quantum phase transitions. In this study, we present compelling evidence of profound structural and transport phase shifts in a recently uncovered oxide Dirac semimetal, SrNbO3, achieved through effective in situ anion doping. Notably, a remarkable increase in resistivity of more than three orders of magnitude at room temperature is observed upon nitrogen-doping. The extent of electronic modulation in SrNbO3 is strongly correlated with the misfit strain, underscoring its phase instability to both chemical doping and crystallographic symmetry variations. Using first-principles calculations, we discern that elevating the level of nitrogen doping induces an upward shift in the conductive bands of SrNbO3-?N?. Consequently, a transition from a metallic state to an insulating state becomes apparent as the nitrogen concentration reaches a threshold of 1/3. This investigation sheds light on the potential of anion engineering in oxide Dirac semimetals, offering pathways for manipulating their physical properties. These insights hold promise for future applications that harness these materials for tailored functionalities.

Hong, Haitao↗
DOE OSTI2021

OptimiSM v.0.0.1

SAND2021-15062 O OptimiSM is a library for solving solid mechanics problems and for rapid prototyping of computational solid mechanics research ideas. It was originally created to facilitate research and development of algorithms for contact mechanics. It is a realization of a philosophy that problems can be represented with variational least-action principles, even when there is history-dependent material behavior and contact. This opens the door to using powerful methods from non-convex optimization for solutions. The code was successful in its original purpose, but along the way, the authors realized that it could be of general use to the computational solid mechanics research community as an R&D platform. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Talamini, Brandon↗
DOE OSTI2025

Robustness of quasi-symmetry along parametric boundary variation

Quasi-symmetry is a guiding principle to modern stellarator optimization for improved plasma confinement. However, the robustness of optimized configurations, which can be crucial for maintaining performance under diverse engineering constraints and practical limitations, has received relatively little attention. Here, in this work, we present various case studies on this robustness, by investigating variations in 1/ v neoclassical transport when plasma configurations are smoothly altered across distinct optimized targets. These targets, optimized from different families—quasi-axisymmetric, quasi-helical, and quasi-isodynamic—are approximately matched in major radius as part of a flexible stellarator design. Our study shows that an optimized target does not always represent a local minimum in transport and that the robustness of a local minimum when present can vary significantly. Furthermore, there are configurations which belong to no established families but have transport levels as low as those of optimized targets. These results highlight the importance of conducting extended searches with key parametric variations around optimized configurations, to ensure its robustness as well as flexibility if desired.

neoclassical transport↗
DOE OSTI2022

Visualizing the interplay of Dirac mass gap and magnetism at nanoscale in intrinsic magnetic topological insulators

In intrinsic magnetic topological insulators, Dirac surface-state gaps are prerequisites for quantum anomalous Hall and axion insulating states. Unambiguous experimental identification of these gaps has proved to be a challenge, however. Here, we use molecular beam epitaxy to grow intrinsic MnBi 2 Te 4 thin films. Using scanning tunneling microscopy/spectroscopy, we directly visualize the Dirac mass gap and its disappearance below and above the magnetic order temperature. We further reveal the interplay of Dirac mass gaps and local magnetic defects. We find that, in high defect regions, the Dirac mass gap collapses. Ab initio and coupled Dirac cone model calculations provide insight into the microscopic origin of the correlation between defect density and spatial gap variations. Finally, this work provides unambiguous identification of the Dirac mass gap in MnBi 2 Te 4 and, by revealing the microscopic origin of its gap variation, establishes a material design principle for realizing exotic states in intrinsic magnetic topological insulators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Mirror dark sector solution of the Hubble tension with time-varying fine-structure constant

Here we explore a model introduced by Cyr-Racine et al. that resolves the Hubble tension by invoking a “mirror world” dark sector with energy density a fixed fraction of the “ordinary” sector of Λ+ cold dark matter. Although it reconciles cosmic microwave background and large-scale structure observations with local measurements of the Hubble constant, H 0 , the model requires a value of the primordial helium mass fraction, Y p = 0.170 ± 0.025, that is discrepant with observations and with the predictions of big bang nucleosynthesis. We consider a variant of the model with standard helium mass fraction but with the value of the electromagnetic fine-structure constant, α, slightly different during photon decoupling from its present value. If α at that epoch is lower than its current value by Δα ≃ -2 x 10 -5 , then we can achieve the same Hubble tension resolution as in the work by Cyr-Racine et al. but with consistent helium abundance. As an example of such time evolution of α, we consider a toy model of an ultralight scalar field, with mass m < 4 x 10 -29 eV, coupled to electromagnetism, which evolves after photon decoupling at redshift z ≃ 10 3 , and that with appropriate coupling appears to be consistent with late-time constraints on α variation and the weak equivalence principle.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Evaluation of Leak Detection Technologies for Low Global Warming Potential (GWP), Flammable Refrigerants

Current commercial refrigeration systems use refrigerants with global warming potential (GWP) values ranging from 1250 to 4000. The emergence of low GWP alternatives (GWP <150) is expected to significantly reduce direct emissions in this sector, playing a crucial role in the ongoing electrification and decarbonization initiatives. However, many of these low GWP alternatives pose a flammability risk, necessitating robust sensing solutions to ensure the reliable and safe operation of the equipment. This paper examines various sensing mechanisms suitable for potential applications in systems that employ flammable refrigerants, specifically those designated as A2L class. It provides a summary of A2L refrigerants and their properties, followed by a comprehensive review of sensor classes, covering their working principles, features, advantages, and limitations. Additionally, the article delves into key performance characteristics such as accuracy, selectivity, sensitivity, dynamic characteristics, and durability, among other properties. The article discusses areas for improvement and suggests corresponding approaches for potential sensors in facilitating the successful adoption of flammable refrigerants. Finally, this paper presents the latest findings from experimental evaluation of 5 different sensing principles in detecting the composition variation as a result of various operational conditions. Reliability and sensitivity of the sensor in responding to shifts in true composition and the resultant LFL value is also discussed.

Reshniak, Viktor↗
DOE OSTI2025

HydroGEN STCH node support: FY25 Q1 quarterly report – LLNL's input

For reducing the cost of H 2 production by enhancing the efficiency and cyclability of water splitting reactions, we aim to understand the relationship between microstructure and surface reaction performance on perovskite oxides. In this quarter, mesoscale simulations were carried with various 3D porous microstructures to characterize the microstructure-dependent H 2 production. During the water splitting process, incoming gas including steam (H 2 O) passes through porous regions, which can induce local pressure variations according to the Bernoulli’s principle. Accordingly, we developed the mesoscale simulation method to characterize the reaction kinetics with inhomogeneous local pressures within porous microstructures.

08 HYDROGEN↗
DOE OSTI2025

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Insight into ideal shear strength of Ni-based dilute alloys using first-principles calculations and correlational analysis

Here the present work examines the effect of alloying elements (denoted X) on the ideal shear strength for 26 dilute Ni-based alloys, Ni11X, as determined by first-principles calculations of pure alias shear deformations. The variations in ideal shear strength are quantitatively explored with correlational analysis techniques, showing the importance of atomic properties such as size and electronegativity. The shear moduli of the alloys are affirmed to show a strong linear relationship with their ideal shear strengths, while the shear moduli of the individual alloying elements were not indicative of alloy shear strength. Through combination with available ideal shear strength data on Mg alloys, a potential application of the Ni alloy data is demonstrated in the search for a set of atomic features suitable for machine learning applications to mechanical properties. As another illustration, the calculated Ni ideal shear strengths play a key role in a predictive multiscale framework for deformation behavior of single crystal alloys at large strains, as shown by simulated stress–strain curves.

36 MATERIALS SCIENCE↗