Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

27 AlNMR chemical shift of Al(OH)$^{-}_{4}$ calculated from first principles: Assessment of error cancellation in chemically distinct reference and target systems

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O)$^{3+}_{6}$, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)$^{-}_{4}$ requires an understanding of how these errors cancel, in order to define the limits of accurately predicting 27 Al chemical shielding in Al(OH)$^{-}_{4}$. In this work, we estimate the absolute shielding of the 27 Al nucleus in Al(OH)$^{-}_{4}$ at the coupled cluster level (515.1± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab-initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O)$^{3+}_{6}$ and Al(OH)$^{-}_{4}$ ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the 27 Al chemical shift of Al(OH)$^{-}_{4}$ from first principles.

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

From predicting to learning dissipation from pair correlations of active liquids

Active systems, which are driven out of equilibrium by local non-conservative forces, can adopt unique behaviors and configurations. Towards designing such materials, an important challenge is to precisely connect the static structure of active systems to the dissipation of energy induced by the local driving. Here, we use tools from liquid-state theories and machine learning to take on these challenges. We first demonstrate analytically for an isotropic active matter system that dissipation and pair correlations are closely related when driving forces behave like an active temperature. We then extend a nonequilibrium mean-field framework for predicting these pair correlations which, unlike most existing approaches, is applicable even for strongly interacting particles and far from equilibrium, to predict dissipation in these systems. Based on this theory, we reveal analytically a robust relation between dissipation and structure which holds even as the system approaches a nonequilibrium phase transition. Finally, we construct a neural network which maps static configurations of particles to their dissipation rate without any prior knowledge of the underlying dynamics. Furthermore, our results open novel perspectives on the interplay between dissipation and organization out-of-equilibrium.

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

Hidden singularities in spontaneously polarized molecular solids

Films of dipolar molecules formed by physical vapor deposition are, in general, spontaneously polarized, with implications ranging from electron transfer in molecular optoelectronic devices to the properties of astrochemical ices in the interstellar medium. Polarization arises from dipole orientation, which should intuitively decrease with increasing deposition temperature, T. However, it is experimentally found that minimum or maximum values in polarization vs T may be observed for cis-methyl formate, 1-propanol, and ammonia. A continuous analytic form of polarization vs T is developed, which has the property that it is not differentiable at all T. The minima and maxima in polarization vs T are marked by singularities in the differential of this analytic form. In conclusion, this exotic behavior is presently unique to films of dipolar species and has not been reported, for example, in the related magnetic phases of spin glasses.

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

Surface coverage dynamics for reversible dissociative adsorption on finite linear lattices

Dissociative adsorption onto a surface introduces dynamic correlations between neighboring sites not found in non-dissociative absorption. We study surface coverage dynamics where reversible dissociative adsorption of dimers occurs on a finite linear lattice. We derive analytic expressions for the equilibrium surface coverage as a function of the number of reactive sites, N, and the ratio of the adsorption and desorption rates. Using these results, we characterize the finite size effect on the equilibrium surface coverage. For comparable N’s, the finite size effect is significantly larger when N is even than when N is odd. Moreover, as N increases, the size effect decays more slowly in the even case than in the odd case. The finite-size effect becomes significant when adsorption and desorption rates are considerably different. These finite-size effects are related to the number of accessible configurations in a finite system where the odd-even dependence arises from the limited number of accessible configurations in the even case. We confirm our analytical results with kinetic Monte Carlo simulations. We also analyze the surface-diffusion case where adsorbed atoms can hop into neighboring sites. As expected, the odd-even dependence disappears because more configurations are accessible in the even case due to surface diffusion.

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

A Gaussian field approach to the solvation of spherical ions in electrolyte solutions

In this paper, the electrostatic response of an electrolyte solution to a spherical ion is studied with a Gaussian field theory. In order to capture the ionic correlation effect in concentrated solutions, the bulk dielectric response function is described by a two-Yukawa response function. The modified response function of the solution is solved analytically in the spherical geometry, from which the induced charge density and the electrostatic energy are also derived analytically. Comparisons with results for small ions in electrolyte solutions from the hyper-netted chain theory demonstrate the validity of the Gaussian field theory.

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

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
DOE OSTI2025

Laser-spectroscopy testbed for impurity monitoring in liquid metal-cooled fast reactors

A significant challenge in sodium-cooled fast reactors is controlling impurities, in particular oxygen impurities, within the sodium coolant, as they can accelerate corrosion and indicate leaks. Optical methods offer the potential to rapidly detect small concentrations of both gaseous and metal impurities that accelerate corrosion, plug coolant channels, and lead to increased activation of isotopes in the coolant. We present the design and performance of an apparatus designed to enable the application of multiple optical analytical techniques, such as laser-induced breakdown spectroscopy, to detect elemental impurities in the sodium melt with high sensitivity. Here, we experimentally demonstrate the detection of characteristic sodium and oxygen spectral lines in liquid sodium, which sets the stage for the optimization of its analytical sensitivity. A robust sensor of this type integrated with the sodium cooling loop has the potential to significantly improve the safety and operational efficiency of generation IV nuclear reactors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Ion Exchange of Selected Group II Metals and Lead by Crystalline Silicotitanate and Competition for Cs Exchange Sites

A series of batch contact tests were conducted to evaluate the exchange behavior of Ba, Ca, Pb, and Sr onto crystalline silicotitanate (CST) in support of an expedited Cs removal and pretreatment system at the Hanford site. Binary Na/M 2+ and ternary Na/Cs/M 2+ isotherms were generated to understand selectivity, capacity, and competitive impact of each analyte on Cs uptake from a simple 1 M NaOH/4.6 M NaNO 3 simulant. Analyte loading from a 0.1 M NaOH/5.5 M NaNO 3 simulant was assessed to determine the effect of hydroxide concentration on binary Na/M 2+ isotherms. Finally, results from binary and ternary isotherms indicated that group II metals, and Pb do not impact CST performance toward CST at concentrations expected in Hanford tank waste supernate.

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

The force explosion condition is consistent with spherically symmetric CCSN explosions

ABSTRACT One of the major challenges in core-collapse supernova (CCSN) theory is to predict which stars explode and which collapse to black holes. The analytic force explosion condition (FEC) shows promise in predicting which stars explode in that the FEC is consistent with CCSN simulations that use the light-bulb approximation for neutrino heating and cooling. In this follow-up manuscript, we take the next step and show that the FEC is consistent with the explosion condition when using actual neutrino transport in gr1d simulations. Since most 1D simulations do not explode, to facilitate this test, we enhance the heating efficiency within the gain region. To compare the analytic FEC and radiation-hydrodynamic simulations, this manuscript also presents a practical translation of the physical parameters. For example: we replace the neutrino power deposited in the gain region, Lντg, with the net neutrino heating in the gain region; rather than assuming that $\dot{M}$ is the same everywhere, we calculate $\dot{M}$ within the gain region; and we use the neutrino opacity at the gain radius. With small, yet practical modifications, we show that the FEC predicts the explosion conditions in spherically symmetric CCSN simulations that use neutrino transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

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

Diaboloidal mirrors: algebraic solution and surface shape approximations

A new type of optical element that can focus a cylindrical wave to a point focus (or vice versa) is analytically described. Such waves are, for example, produced in a beamline where light is collimated in one direction and then doubly focused by a single optic. A classical example in X-ray optics is the collimated two-crystal monochromator, with toroidal mirror refocusing. The element here replaces the toroid, and in such a system provides completely aberration free, point-to-point imaging of rays from the on-axis source point. We present an analytic solution for the mirror shape in its laboratory coordinate system with zero slope at the centre, and approximate solutions, based on bending an oblique circular cone and a bent right circular cylinder, that may facilitate fabrication and metrology.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Use of carbonyldiimidazole as a derivatization agent for the detection of pinacolyl alcohol, a forensic marker for Soman, by EI-GC–MS and LC-HRMS in official OPCW proficiency test matrices

Pinacolyl alcohol (PA), a key forensic marker for the nerve agent Soman (GD), is a particularly difficult analyte to detect by various analytical methods. In this work, we have explored the reaction between PA and 1,1'-carbonyldiimidazole (CDI) to yield pinacolyl 1H-imidazole-1-carboxylate (PIC), a product that can be conveniently detected by gas chromatography–mass spectrometry (GC–MS) and liquid chromatography-high-resolution mass spectrometry (LC-HRMS). Regarding its GC–MS profile, this new carbamate derivative of PA possesses favorable chromatographic features such as a sharp peak and a longer retention time (RT = 16.62 min) relative to PA (broad peak and short retention time, RT = 4.1 min). Additionally, the derivative can also be detected by LC-HRMS, providing an avenue for the analysis of this chemical using this technique where PA is virtually undetectable unless present in large concentrations. From a forensic science standpoint, detection of this low molecular weight alcohol signals the past or latent presence of the nerve agent Soman (GD) in a given matrix (i.e., environmental or biological). The efficiency of the protocol was tested separately in the analysis and detection of PA by EI-GC–MS and LC-HRMS when present at a 10 μg/mL in a soil matrix featured in the 44th PT and in a glycerol-rich liquid matrix featured in the 48th Official Organization for the Prohibition of Chemical Weapons (OPCW) Proficiency Test when present at a 5 μg/mL concentration. In both scenarios, PA was successfully transformed into PIC, establishing the protocol as an additional tool for the analysis of this unnatural and unique nerve agent marker by GC–MS and LC-HRMS.

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

Gas‐Phase Ion‐Molecule Interactions in a Collision Reaction Cell with ICP‐MS / MS : Investigations with CO 2 as the Reaction Gas

Carbon dioxide (CO 2 ) was used as a reaction gas to investigate the gas‐phase ion‐molecule interactions using the Agilent 8900 ICP‐MS/MS. A solution containing forty‐five elements representative of the periodic table was used to supply the ions to react with CO 2 in the collision/reaction cell (CRC). The only significant product ions formed were monoxides. The general reactivity was shown to be consistent with density functional theory (DFT)‐predicted reaction enthalpies, such that all predicted exothermic reactions produced product ions at levels of at least 1% of the unreacted ion. Most endothermic reactions observed had sufficient kinetic energy in excess of the reaction enthalpies. Our results suggest that reaction enthalpy is a reasonable predictor of reactivity with CO 2 on the timescales of the interactions in non‐thermal ICP‐MS/MS systems. The ease and rapidity of data collection with the ICP‐MS/MS and DFT calculations using the NWChem suite has value given the scarcity of thermochemical data of CO 2 reactions in the literature. These studies are especially useful for the identification of targeted reaction chemistries to be leveraged for analytical method development, such as for the inline separation of isobaric interferences from analytes of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multilevel Convergence Analysis of Multigrid-Reduction-in-Time

This study presents a multilevel convergence framework for multigrid-reduction-in-time (MGRIT) as a generalization of previous two-grid estimates. The framework provides a priori upper bounds on the convergence of MGRIT V- and F-cycles, with different relaxation schemes, by deriving the respective residual and error propagation operators. The residual and error operators are functions of the time-stepping operator, analyzed directly and bounded in the norm, both numerically and analytically. We present various upper bounds of different computational cost and varying sharpness. These upper bounds are complemented by proposing analytic formulae for the approximate convergence factor of V-cycle algorithms that take the number of fine grid time points, the temporal coarsening factors, and the eigenvalues of the time-stepping operator as parameters. The paper concludes with supporting numerical investigations of parabolic (anisotropic diffusion) and hyperbolic (wave equation) model problems. We assess the sharpness of the bounds and the quality of the approximate convergence factors. Observations from these numerical investigations demonstrate the value of the proposed multilevel convergence framework for estimating MGRIT convergence a priori and for the design of a convergent algorithm. We further highlight that observations in the literature are captured by the theory, including that two-level Parareal and multilevel MGRIT with F-relaxation do not yield scalable algorithms and the benefit of a stronger relaxation scheme. An important observation is that with increasing numbers of levels MGRIT convergence deteriorates for the hyperbolic model problem, while constant convergence factors can be achieved for the diffusion equation. The theory also indicates that L-stable Runge--Kutta schemes are more amendable to multilevel parallel-in-time integration with MGRIT than A-stable Runge--Kutta schemes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Comparing Experimentally-Measured Sand’s Times with Concentrated Solution Theory Predictions in a Polymer Electrolyte

We compare the electrochemically measured Sand’s time, the time required for the cell potential to diverge when the applied current density exceeds the limiting current, with theoretical predictions for a 0.47 M poly(ethylene oxide) (5 kg mol −1 )/LiTFSI electrolyte. The theoretical predictions are made using concentrated solution theory which accounts for both concentration polarization and polymer motion, using independently measured parameters that depend on concentration, c : conductivity ( κ ), salt diffusion coefficient ( D ), cationic transference number with respect to the solvent velocity ( t + 0 ), thermodynamic factor 1 + dln f ± dln c , and partial molar volume of the salt ( V ̅ ); f ± is the mean molar activity coefficient of the salt. We find quantitative agreement between experimental data and theoretical predictions. We derive a generalized analytical expression for Sand’s time for electrolytes based on dilute solution theory. This expression correctly predicts the divergence of the Sand’s time at the limiting current, in agreement with experimental data and concentrated solution theory predictions. When the applied current is large compared to the limiting current, the analytical expression approaches the standard expression for Sand’s time used in the literature.

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

Data for A Hybrid Biophysical-Machine Learning Framework for Diurnal Surface Energy Flux Estimation Using Proximal Sensing

Thermal infrared-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy (LE) and sensible heat (H) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal infrared data sets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of a ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R2 = 0.81–0.94) and H (R2 = 0.46–0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical-machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.

AI/ML↗
DOE OSTI2026

Forensic characterization of surrogate nuclear explosion debris: radiochemical and spectroscopic strategies for method validation

Surrogate nuclear explosion debris (SNED) has emerged as a critical platform for advancing post-detonation nuclear forensic analysis in the absence of readily accessible historic materials. SNED enables controlled investigation and validation of analytical methodologies used to interrogate the chemical, isotopic, radiological, and microstructural signatures preserved in nuclear explosion debris. This review presents an integrated assessment of destructive and non-destructive analytical techniques commonly employed within decision-driven nuclear forensic workflows. Each technique is discussed individually while highlighting how it contributes to different stages of post-detonation analysis. Core methods – including gamma and alpha spectrometry, ICP-MS, TIMS, SIMS, SEM-EDS, XRF, LIBS, vibrational spectroscopy, and X-ray absorption spectroscopy – are critically evaluated with respect to forensic maturity, information content, and matrix limitations. Emphasis is placed on the role of SNED in benchmarking multi-modal workflows and identifying gaps in reproducing heterogeneity, fractionation, and radiation-driven evolution relevant to forensic attribution.

X-ray spectroscopic methods↗
DOE OSTI2020

Sequim Campus Master Plan

Based in Richland, Washington, Pacific Northwest National Laboratory (PNNL) is one of ten U.S. Department of Energy (DOE) Office of Science (SC) national laboratories. PNNL is DOE’s premier chemistry, earth sciences, and data analytics laboratory, delivering vital mission impacts in energy resiliency and national security. PNNL has 19 core capabilities, each a powerful combination of world-class staff, state-of-the- art equipment, and mission-ready facilities. These capabilities represent a collective set of skills and a body of world-leading scientific and engineering work that provides exceptional value and mission delivery to DOE, the National Nuclear Security Administration (NNSA), the U.S. Department of Homeland Security (DHS), and other federal agencies and industry through strategic partnership projects on a not-to-interfere basis. Currently operated by Battelle Memorial Institute, PNNL occupies more than 2.3 million gross square feet (GSF) of buildings and employs approximately 4,700 staff members. PNNL has implemented a rolling 10-year campus strategy, initiated in fiscal year 2013, designed to deliver the foundation for the multidisciplinary science and engineering expertise, equipment, instrumentation, facilities, and infrastructure required to continue providing exceptional value and mission delivery. This strategy is incorporated into PNNL’s annual Laboratory Plan, which serves as DOE-SC’s equivalent of the Five- Year Site Plan required by DOE Order 430.1C Chg. 1, Real Property Asset Management. The PNNL Sequim Campus is located in Clallam County, just outside of Sequim, Washington. PNNL Sequim builds upon a rich history of research related to marine and coastal resources, environmental chemistry, water resources modeling, ecotoxicology, biotechnology, and national security. Facilities at the PNNL Sequim Campus were designed to take advantage of two unique locational assets: the pristine coastal ocean water in Sequim Bay and the region’s exceptionally clean air. Research facilities on the waterfront are supplied with up to 200 gallons per minute (gpm) of seawater and/or groundwater that enables research on aquatic systems ranging from fully seawater to fresh water, with an onsite water treatment system that removes all biological and hazardous chemical residues prior to discharge. Multiple large seawater tanks outside the waterfront laboratories, along with a boat ramp, pier, and floating dock, enable lab-to-field research. The uplands facility takes advantage of the amazing air quality and is built and operated to enable ultratrace chemical/radiological analytical chemistry. This combination of capabilities is not available at any other marine laboratory in the nation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗