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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Antiphase boundary migration as a diffusion mechanism in a P3 sodium layered oxide

The electrode materials of Na-ion batteries must exhibit fast ion transport, a property that is often negatively impacted when mobile cations order. Many layered Na (and K) intercalation compounds adopt the P3 crystal structure, which tends to stabilize cation-vacancy orderings over a wide range of Na concentrations. The ordered phases have been predicted to consist of domains that are periodically separated by antiphase boundaries (APBs) to accommodate variations in Na concentration. The mechanisms and rates with which Na diffuses within such ordered phases of P3 are not understood. We have conducted a first-principles study of Na diffusion in the P3 structure, using Na x CoO 2 as a model system. We find that most simple nearest-neighbor hops in Na-vacancy ordered phases in P3 result in unstable configurations and are therefore forbidden. We identify an alternative mechanism of Na transport that is mediated by APB migration via the formation and expansion of kinks along domain boundaries. The limiting kinetic barriers for this mechanism are quite low (0.03 to 0.30 eV). Vacancy defects are not required, and furthermore have formation energies that are higher than those of the APB kink defects. Furthermore, our results suggest that APB migration is a fundamental diffusion mechanism in the P3 structure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗
DOE OSTI2023

First-Order Individual Gas Sensors as Next Generation Reliable Analytical Instruments

It is conventionally expected that the performance of existing gas sensors may degrade in the field compared to laboratory conditions because (i) a sensor may lose its accuracy in the presence of chemical interferences and (ii) variations of ambient conditions over time may induce sensor-response fluctuations (i.e., drift). Breaking this status quo in poor sensor performance requires understanding the origins of design principles of existing sensors and bringing new principles to sensor designs. Existing gas sensors are single-output (e.g., resistance, electrical current, light intensity, etc.) sensors, also known as zero-order sensors (Karl Booksh and Bruce R. Kowalski, Analytical Chemistry, DOI: 10.1021/ac00087a718). Any zero-order sensor is undesirably affected by variable chemical background and sensor drift that cannot be distinguished from the response to an analyte. To address these limitations, we are developing multivariable gas sensors with independent responses, which are first-order analytical instruments. Here, we demonstrate self-correction against drift in two types of first-order gas sensors that operate in different portions of the electromagnetic spectrum. Our radiofrequency sensors utilize dielectric excitation of semiconducting metal oxide materials on the shoulder of their dielectric relaxation peak and achieve self-correction of the baseline drift by operation at several frequencies. Our photonic sensors utilize nanostructured sensing materials inspired by Morpho butterflies and achieve self-correction of the baseline drift by operation at several wavelengths. These principles of self-correction for drift effects in first-order sensors open opportunities for diverse emerging monitoring applications that cannot afford frequent periodic maintenance that is typical of traditional analytical instruments.

Potyrailo, Radislav A.↗
DOE OSTI2025

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Learning Grain-Boundary Segregation: From First Principles to Polycrystals

The segregation of solute atoms at grain boundaries (GBs) can strongly impact the structural and functional properties of polycrystals. Yet, due to the limited availability of simulation tools to study polycrystals at the atomistic scale (i.e., interatomic potentials), there is a minimal understanding of the variation of solute segregation tendencies across the very complex space of GB microenvironments and the large range of alloys in which it can occur. Here, we develop an algorithmic framework that can directly learn the full spectrum of segregation energies for a metal solute atom in a metal polycrystal from ab initio methods, bypassing the need for alloy interatomic potentials. This framework offers a pathway to a comprehensive catalog of GB solute segregation with quantum accuracy, for the entire alloy space. As an initial demonstration in this pursuit, we build an extensive GB segregation database for aluminum-based alloys across the periodic table, including dozens of alloys for which there are substantially no prior data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Time-dependent relaxed magnetohydrodynamics: Inclusion of cross helicity constraint using phase-space action

In this work, a phase-space version of the ideal magnetohydrodynamic (MHD) Lagrangian is derived from first principles and shown to give a relabeling transformation when a cross-helicity constraint is added in Hamilton's Action Principle. A new formulation of time-dependent “relaxed” magnetohydrodynamics is derived using microscopic conservation of mass and macroscopic constraints on total magnetic helicity, cross helicity, and entropy under variations of density, pressure, fluid velocity, and magnetic vector potential. This gives Euler–Lagrange equations consistent with previous work on both ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether E+u×B=0 should be a constraint is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Striped Jets in Post–Neutron Star Merger Systems

Abstract Models invoking magnetic reconnection as the particle acceleration mechanism within relativistic jets often adopt a gradual energy dissipation profile within the jet. However, such a profile has yet to be reproduced in first-principles simulations. Here we perform a suite of 3D general relativistic magnetohydrodynamic simulations of post–neutron star merger disks with an initially purely toroidal magnetic field. We explore the variations in both the microphysics (e.g., nuclear recombination, neutrino emission) and system parameters (e.g, disk mass). In all of our simulations, we find the formation of magnetically striped jets. The stripes result from the reversals in the poloidal magnetic flux polarity generated in the accretion disk. The simulations display large variations in the distributions of stripe duration, τ , and power, 〈 P Φ 〉. We find that more massive disks produce more powerful stripes, the most powerful of which reaches 〈 P Φ 〉 ∼ 10 49 erg s −1 at τ ∼ 20 ms. The power and variability that result from the magnetic reconnection of the stripes agree with those inferred in short-duration gamma-ray bursts. We find that the dissipation profile of the cumulative energy is roughly a power law in both radial distance, z , and τ , with a slope in the range of ∼1.7–3; more massive disks display larger slopes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Representation Fusion Framework for Decoupling Diagnostic Information in Multimodal Learning

Modern medicine increasingly relies on multimodal data, ranging from clinical notes to imaging and genomics, to guide diagnosis and treatment. However, integrating these heterogeneous data sources in a principled and interpretable manner remains a major challenge. We present MODES (Multi-mOdal Disentangled Embedding Space), a representation fusion framework that explicitly separates shared and modality-specific factors of variation, offering a structured latent space for multimodal information that improves both prediction and interpretability. By leveraging pre-trained unimodal foundation models, MODES mitigates the dependency on extensive paired datasets, crucial in data-scarce clinical settings. We introduce a masking strategy that optimizes representation dimensionality by eliminating low-information dimensions, to achieve compact, information-rich representations. Our framework demonstrates superior performance in predicting diagnoses and phenotypes compared to unimodal and conventional fusion models. MODES also enables robust diagnostic inference in missing data scenarios, offering an opportunity toward interpretable and efficient multimodal diagnostics in personalized healthcare.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

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

Efficacious Qubit Mappings for Quantum Simulations of the 12 C Rotational Band

Solving atomic nuclei from first principles places enormous demands on computational resources, which grow exponentially with increasing number of particles and the size of the space they occupy. We present first quantum simulations based on the variational quantum eigensolver for the low-lying structure of the 12C nucleus that provide acceptable bound-state energies even in the presence of noise. Here, we achieve this by taking advantage of two critical developments. First, we utilize an almost perfect symmetry of atomic nuclei that, in a complete symmetry-adapted basis, drastically reduces the size of the model space. Second, we use the efficacious Gray encoding, for which it has been recently shown that it is resource efficient, especially when coupled with a near band-diagonal structure of the nuclear Hamiltonian.

Encoding↗
DOE OSTI2020

PyDA: A Hands-On Introduction to Dynamical Data Assimilation with Python

Dynamic data assimilation offers a suite of algorithms that merge measurement data with numerical simulations to predict accurate state trajectories. Meteorological centers rely heavily on data assimilation to achieve trustworthy weather forecast. With the advance in measurement systems, as well as the reduction in sensor prices, data assimilation (DA) techniques are applicable to various fields, other than meteorology. However, beginners usually face hardships digesting the core ideas from the available sophisticated resources requiring a steep learning curve. In this tutorial, we lay out the mathematical principles behind DA with easy-to-follow Python module implementations so that this group of newcomers can quickly feel the essence of DA algorithms. We explore a series of common variational, and sequential techniques, and highlight major differences and potential extensions. We demonstrate the presented approaches using an array of fluid flow applications with varying levels of complexity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Collisional simulations of the modulator section in coherent electron cooling

The first section of any coherent electron cooling (CeC) system is the modulator, where the density of the electron beam is modulated by the copropagating ion beam. This density modulation is a result of Coulomb collisions between the individual particles of the two beams. The pairwise, stochastic part of the interactions impacts the overall performance of the CeC process. We present the first simulations of the density modulations of the electron beams from a collisional picture of the dynamics, considering the proof-of-principle CeC experiments at Brookhaven National Laboratory. These simulations were performed using PHAD, which is the first efficient, large-scale collisional numerical method in beam physics that we have previously developed and benchmarked. Realistic beam distributions and external fields have been optimized to provide strong modulation signals necessary for variations of coherent electron cooling systems. Cooling performance limits and potential collisionless simulation pitfalls are pointed out. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Device response principles and the impact on energy resolution of epitaxial quantum dot scintillators with monolithic photodetector integration

Abstract Epitaxial quantum dot (QD) scintillator crystals with picosecond-scale timing and high light yield have been created for medical imaging, high energy physics and national security applications. Monolithic photodetector (PD) integration enables the sensing of photons generated within the waveguiding crystal and allows a wide range of scintillator-photodetector coupling geometries. Until recently, these doubly novel devices have suffered from complex, high variance responses to monoenergetic sources which significantly reduces their precision and accuracy. The principles governing the overall device response have now been discerned and embodied by an expression derived within a geometrical optics framework which considers optical properties, surface roughness and photodetector coupling geometry. Response variation due to these factors was sufficiently reduced to obtain material-related energy resolution values of 2.4% with alpha particles. These findings place energy resolution alongside luminescence timescale, photon yield, and radiation hardness as outstanding properties of these engineered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Theoretical investigation of charge density wave instability in CuS 2

The existence of a charge density wave (CDW) in transition-metal dichalcogenide (TMDC) CuS 2 has remained undetermined since its first experimental synthesis nearly 50 years ago. Despite conflicting experimental literature regarding its low-temperature structure, there exists no theoretical investigation of the phonon properties and lattice stability of this material. By studying the first-principles electronic structure and phonon properties of CuS 2 at various electronic temperatures, we identify temperature-sensitive soft phonon modes which unveil a previously unreported Kohn anomaly at approximately 100 K. Variation of the electronic temperature shows the presence of two distinct phases, characterized at low electronic temperature by a 2 × 2 × 2 periodic charge modulation associated with the motion of the S 2 dimers. We find this is driven by a slight orbital occupation imbalance of the copper d and sulfur p orbitals, reminiscent of the Jahn-Teller effect in finite systems. Investigation of the Fermi surface presents a potential Fermi surface nesting vector related to the location of the Kohn anomaly and observed band splittings in the unfolded band structure. The combination of these results suggests a strong possibility of CDW order in CuS 2 . Further study of CuS 2 in monolayer form finds no evidence of a CDW phase, as the identified bulk periodic distortions cannot be realized in two dimensions. This behavior sets this material apart from other transition-metal dichalcogenide materials, which exhibit a charge density wave phase down to the two-dimensional limit. As CDW in TMDC materials is considered to compete with superconductivity, the lack of a CDW in monolayer CuS 2 suggests the possibility of enhanced superconductivity relative to other transition-metal dichalcogenides. Altogether, our work identifies CuS 2 as a previously unrealized candidate to study the interplay of superconductivity, CDW order, and dimensionality.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evidence for Ga clusters in β-Ga 2 O 3 from Raman spectroscopy and density functional theory

Monoclinic gallium oxide (β-Ga 2 O 3 ) single crystals have a Raman mode at ∼250 cm −1 that is strongly correlated with free-electron density. Prior work attributed this peak to an electronic excitation of a shallow donor impurity band. However, heavily n-type thin films grown by metalorganic chemical vapor deposition or molecular beam epitaxy do not have the peak. In the present work, an alternate model is proposed: the 250 cm −1 Raman peak arises from Ga clusters, defined as two or more Ga atoms that form Ga–Ga bonds. Raman mapping reveals variations in the frequency that are consistent with a distribution of cluster sizes. The intensity of the peak decreases as the temperature is raised, attributed to melting of the Ga clusters. First-principles calculations indicate that the 250 cm −1 mode is due to Ga–Ga bond-stretching vibrations. As the Fermi energy is raised, the formation of Ga–Ga dimers becomes energetically favorable, explaining the correlation between n-type conductivity and the appearance of the Raman peak.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Proximity-enhanced valley Zeeman splitting at the WS 2 /graphene interface

The valley Zeeman physics of excitons in monolayer transition metal dichalcogenides provides valuable insight into the spin and orbital degrees of freedom inherent to these materials. Being atomically-thin materials, these degrees of freedom can be influenced by the presence of adjacent layers, due to proximity interactions that arise from wave function overlap across the 2D interface. Here, we report 60 T magnetoreflection spectroscopy of the A- and B- excitons in monolayer WS 2 , systematically encapsulated in monolayer graphene. While the observed variations of the valley Zeeman effect for the A- exciton are qualitatively in accord with expectations from the bandgap reduction and modification of the exciton binding energy due to the graphene-induced dielectric screening, the valley Zeeman effect for the B- exciton behaves markedly different. We investigate prototypical WS 2 /graphene stacks employing first-principles calculations and find that the lower conduction band of WS 2 at the <?CDATA $K/K^{^{\prime}}$?> K / K ′ valleys (the <:tex-math><?CDATA $\mathrm{CB}^-$?> C B − band) is strongly influenced by the graphene layer on the orbital level. Specifically, our detailed microscopic analysis reveals that the conduction band at the Q point of WS 2 mediates the coupling between <?CDATA $\mathrm{CB}^-$?> C B − and graphene due to resonant energy conditions and strong coupling to the Dirac cone. This leads to variations in the valley Zeeman physics of the B- exciton, consistent with the experimental observations. Our results therefore expand the consequences of proximity effects in multilayer semiconductor stacks, showing that wave function hybridization can be a multi-step energetically resonant process, with different bands mediating the interlayer interactions. Such effects can be further exploited to resonantly engineer the spin-valley degrees of freedom in van der Waals and moiré heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Precipitation of dopants on acceptor-doped LaMnO3±δ revealed by defect chemistry from first principles

Perovskite oxides degrade at elevated temperatures while precipitating dopant-rich particles on the surface. A knowledge-based improvement of surface stability requires a fundamental and quantitative understanding of the dopant precipitation mechanism on these materials. We propose that dopant precipitation is a consequence of the variation of dopant solubility between calcination and operating conditions in solid oxide fuel cells (SOFCs) and electrolyzer cells (SOECs). To study dopant precipitation, we use 20% (D = Ca, Sr, Ba)-doped LaMnO3+δ (LDM20) as a model system. We employ a defect model taking input from density functional theory calculations. The defect model considers the equilibration of LDM20 with a reservoir consisting of dopant oxide (DO), peroxide (DO2), and O2 in the gas phase. The equilibrated non-stoichiometry of the A-site and B-site as a function of temperature, T, and oxygen partial pressure, p(O2), reveals three regimes for LDM20: A-site deficient (oxidizing conditions), A-site rich (atmospheric conditions), and near-stoichiometric (reducing conditions). Assuming an initial A/B non-stoichiometry, we compute the dopant precipitation boundaries in a p–T phase diagram. Our model predicts precipitation both under reducing (DO) and under highly oxidizing conditions (DO2). We found precipitation under anodic, SOEC conditions to be promoted by large dopant size, while under cathodic, SOFC conditions precipitation is promoted by initial A-site excess. The main driving forces for precipitation are oxygen uptake by the condensed phase under oxidizing conditions and oxygen release assisted by B-site vacancies under reducing conditions. Possible strategies for mitigating dopant precipitation under in electrolytic and fuel cell conditions are discussed.

Hess, Franziska (ORCID:000000024460643X)↗
DOE OSTI2023

Performance analysis of functionally graded multifunctional piezoelectric energy harvesting microgyroscopes

Functionally graded materials (FGMs) are composite materials with varying material properties in one or more directions. These materials possess distinct characteristics compared to their constituent components. The ability to control material distribution and compositions in FGMs offers improved constraints over the natural frequencies of the system, making them highly suitable for energy harvesting applications. In this study, we focus on a piezoelectric FGM composed of Platinum (Pt) and Lead Titanate Zirconate (PZT). By utilizing FGMs for energy harvesting, we simplify the system from a multilayer structure to a single-layer system, thereby increasing power density by reducing volume. In this work, a power law distribution is used to model the material variation throughout the thickness of the single-layered beam. The governing equations of motion and boundary conditions for FG energy harvesting microgyroscope are derived using Euler-Bernoulli beam theory, the constitutive piezoelectric principle, and the extended Hamilton's principle. To simulate the nonlinear motion of the FG energy harvester, we employ the differential quadrature method (DQM). Subsequently, an investigation is carried out to examine the effects of various factors including material distribution, platinum percentage, base rotation, DC voltage, and electrical load resistance on the energy harvesting microgyroscope with functionally graded materials. Here, the findings suggest that functionally graded energy harvesting gyroscopes can be adjusted to obtain effective energy harvesting and sensing capabilities by carefully selecting the material distribution, input voltage, and electrical load resistance.

42 ENGINEERING↗