Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional modeling”

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 559 records · Page 31

A valuation framework for customers impacted by extreme temperature-related outages

Extreme temperature outages can lead to not just economic losses but also various non-energy impacts (NEI), such as increased mortality rates, property damage, and reduced productivity, due to significant degradation of indoor operating conditions caused by service disruptions. However, existing resilience assessment approaches lack specificity for extreme temperature conditions. They often overlook temperature-related mortality and neglect the customer characteristics and grid response in the calculation, despite the significant influence of these factors on NEI-related economic losses. This paper aims to address these gaps by introducing a comprehensive framework to estimate the impact of resilience enhancement not only on the direct economic losses incurred by customers but also on potential NEI, including mortality and the value of statistical life during extreme temperature-related outages. The proposed resilience valuation integrates customer characteristics and grid response variables based on a scalable grid simulation environment. This study adopts a holistic approach to quantify customer-oriented economic impacts, utilizing probabilistic loss scenarios that incorporate health-related factors and damage/loss models as a function of exposure for valuation. The proposed methodology is demonstrated through comparative resilient outage planning, using grid response models emulating a Texas weather zone during the 2021 winter storm Uri. The case study results show that enhanced outage planning with hardened infrastructure can improve the system resilience and thereby reduce the relative risk of mortality by 16% and save the total costs related to non-energy impacts by 74%. In conclusion, these findings underscore the efficacy of the framework by assessing the financial implications of each case, providing valuable insights for decision-makers and stakeholders involved in extreme-weather related resilience planning for risk management and mitigation strategies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Coevolution of the Ess1-CTD axis in polar fungi suggests a role for phase separation in cold tolerance

Most of the world’s biodiversity lives in cold (-2° to 4°C) and hypersaline environments. To understand how cells adapt to such conditions, we isolated two key components of the transcription machinery from fungal species that live in extreme polar environments: the Ess1 prolyl isomerase and its target, the carboxy-terminal domain (CTD) of RNA polymerase II. Polar Ess1 enzymes are conserved and functional in the model yeast, Saccharomyces cerevisiae. By contrast, polar CTDs diverge from the consensus (YSPTSPS) 26 and are not fully functional in S. cerevisiae. These CTDs retain the critical Ess1 Ser-Pro target motifs, but substitutions at Y1, T4, and S7 profoundly affected their ability to undergo phase separation in vitro and localize in vivo. We propose that environmentally tuned phase separation by the CTD and other intrinsically disordered regions plays an adaptive role in cold tolerance by concentrating enzymes and substrates to overcome energetic barriers to metabolic activity.

59 BASIC BIOLOGICAL SCIENCES↗

hippynn Python Package

hippynn is a python package for defining, training, and applying neural networks to atomistic systems. In particular, it focuses on Hierarchical Interacting Particle Neural Networks (HIP-NNs), a deep learning architecture for atomistic systems. HIP-NNs take input data describing the properties of atomistic systems (often generated using ab-initio quantum mechanics) to learn fast and accurate models. A trained HIP-NN can predict potential energy surfaces, atomic charges, and more. hippynn uses PyTorch for portable high-performance code, including both CPU and GPU support. hippynn allows for extensive customization, including user-defined models and loss functions, to facilitate future research into extensions of HIP-NN as well as other atomistic deep learning models.

Lubbers, Nicolas↗

Cluster model of C 12 in the density functional theory framework

We employ the constrained density functional theory to investigate cluster phenomena for the 12 C nucleus. The proton and neutron densities are generated from the placement of three 4 He nuclei (α particles) geometrically. These densities are then used in a density constrained Hartree-Fock calculation that produces an antisymmetrized state with the same densities through energy minimization. In the calculations no a priori analytic form for the single-particle states is assumed and the full energy density functional is utilized. The geometrical scan of the energy landscape provides the ground state of 12 C as an equilateral triangular configuration of three αs with molecular bond like structures. The use of the nucleon localization function provides further insight to these configurations. One can conclude that these configurations are a hybrid between a pure mean-field and a pure α particle condensate. Furthermore, this development could facilitate density functional theory based fusion calculations with a more realistic 12 C ground state.

6 ≤ A ≤ 19↗

Temperature dependence of nuclear quantum effects on liquid water via artificial neural network model based on SCAN meta-GGA functional

We investigate the temperature dependence of nuclear quantum effects (NQEs) on structural and dynamic properties of liquid water by training a neural network force field using first-principles molecular dynamics (FPMD) based on the strongly constrained and appropriately normed meta-generalized gradient approximation exchange-correlation approximation. The FPMD simulation based on density functional theory has become a powerful computational approach for studying a wide range of condensed phase systems. However, its large computational cost makes it difficult to incorporate NQEs in the simulation and investigate temperature dependence of various properties. To circumvent this difficulty, we use an artificial neural network model and employ the thermostatted ring polymer MD approach for studying the temperature dependence of NQEs on various properties. The NQEs generally bring the radial distribution functions closer to the experimental measurements. Translational diffusivity and rotational dynamics of water molecules are both slowed down by the NQEs. The competing inter-molecular and intra-molecular quantum effects on hydrogen bonds, as discussed by Habershon, Markland, and Manolopoulos [J. Chem. Phys. 131(2), 024501 (2019)], can explain the observed temperature dependence of the NQEs on the dynamical properties in our simulation.

Yao, Yi (ORCID:0000000299576979)↗

Exploring the Structural, Dynamic, and Functional Properties of Metal–Organic Frameworks through Molecular Modeling

This review spotlights the role of atomic-level modeling in research on metal-organic frameworks (MOFs), especially the key methodologies of density functional theory (DFT), Monte Carlo (MC) simulations, and molecular dynamics (MD) simulations. The discussion focuses on how periodic and cluster-based DFT calculations can provide novel insights into MOF properties, with a focus on predicting structural transformations, understanding thermodynamic properties and catalysis, and providing information or properties that are fed into classical simulations such as force field parameters or partial charges. Classical simulation methods, highlighting force field selection, databases of MOFs for high-throughput screening, and the synergistic nature of MC and MD simulations, are described. By predicting equilibrium thermodynamic and dynamic properties, these methods offer a wide perspective on MOF behavior and mechanisms. Additionally, the incorporation of machine learning (ML) techniques into quantum and classical simulations is discussed. These methods can enhance accuracy, expedite simulation setup, reduce computational costs, as well as predict key parameters, optimize geometries, and estimate MOF stability. By charting the growth and promise of computational research in the MOF field, the aim is to provide insights and recommendations to facilitate the incorporation of computational modeling more broadly into MOF research.

36 MATERIALS SCIENCE↗

Modeling Adsorption in Silica Pores via Minkowski Functionals and Molecular Electrostatic Moments

Capillary condensation phenomena are important in various technological and environmental processes. Using molecular simulations, we study the confined phase behavior of fluids relevant to carbon sequestration and shale gas production. As a first step toward translating information from the molecular to the pore scale, we express the thermodynamic potential and excess adsorption of methane, nitrogen, carbon dioxide, and water in terms of the pore’s geometric properties via Minkowski functionals. This mathematical reconstruction agrees very well with molecular simulations data. Our results show that the fluid molecular electrostatic moments are positively correlated with the number of adsorption layers in the pore. Moreover, stronger electrostatic moments lead to adsorption at lower pressures. These findings can be applied to improve pore-scale thermodynamic and transport models.

42 ENGINEERING↗

Piecewise linear approximation with minimum number of linear segments and minimum error: A fast approach to tighten and warm start the hierarchical mixed integer formulation

In several areas of economics and engineering, it is often necessary to fit discrete data points or approximate nonlinear functions with continuous functions. Piecewise linear (PWL) functions are a convenient way to achieve this. PWL functions can be modeled in mathematical problems using only linear and integer variables. Moreover, there is a computational benefit in using PWL functions that have the least possible number of segments. This work proposes a novel hierarchical mixed integer linear programming (MILP) formulation that identifies a continuous PWL approximation with minimum number of linear segments for a given target maximum error. The proposed MILP formulation also identifies the solution with the least maximum error among the solutions with minimum number of segments. Then, this work proposes a fast iterative algorithm that identifies non necessarily continuous PWL approximations by solving O(S log N) linear programming (LP) problems, where N is the number of data points and S is the minimum number of segments in the non necessarily continuous case. This work demonstrates that tight bounds for the MILP problem can be derived from these approximations. Next, a fast algorithm is introduced to transform a non necessarily continuous PWL approximation into a continuous one. Finally, the tight bounds and the continuous PWL approximations are used to tighten and warm start the MILP problem. The tightened formulation is shown in experimental results to be more efficient, especially for large data sets, with a solution time that is up to two orders of magnitude less than the existing literature.

97 MATHEMATICS AND COMPUTING↗

Output-weighted sampling for multi-armed bandits with extreme payoffs

We present a new type of acquisition function for online decision-making in multi-armed and contextual bandit problems with extreme payoffs. Specifically, we model the payoff function as a Gaussian process and formulate a novel type of upper confidence bound acquisition function that guides exploration towards the bandits that are deemed most relevant according to the variability of the observed rewards. This is achieved by computing a tractable likelihood ratio that quantifies the importance of the output relative to the inputs and essentially acts as an attention mechanism that promotes exploration of extreme rewards. Our formulation is supported by asymptotic zero-regret guarantees, and its performance is demonstrated across several synthetic benchmarks, as well as two realistic examples involving noisy sensor network data. Finally, we provide a JAX library for efficient bandit optimization using Gaussian processes.

Science & Technology - Other Topics↗

An Image-Plane Approach to Gravitational Lens Modeling of Interferometric Data

Strong gravitational lensing acts as a cosmic telescope, enabling the study of the high-redshift universe. Astronomical interferometers, such as the Atacama Large Millimeter/submillimeter Array (ALMA), have provided high-resolution images of strongly lensed sources at millimeter and submillimeter wavelengths. To model the mass and light distributions of lensing and source galaxies from strongly lensed images, strong lens modeling for interferometric observations is conventionally performed in the visibility space, which is computationally expensive. In this paper, we implement an image-plane lens modeling methodology for interferometric dirty images by accounting for noise correlations. We show that the image-plane likelihood function produces accurate model values when tested on simulated ALMA observations with an ensemble of noise realizations. We also apply our technique to ALMA observations of two sources selected from the South Pole Telescope survey, comparing our results with previous visibility-based models. Our model results are consistent with previous models for both parametric and pixelated source-plane reconstructions. We implement this methodology for interferometric lens modeling in the open-source software package lenstronomy.

Zhang, Nan [Illinois U., Urbana (main)] (ORCID:000↗

Accessing the Single-Particle Structure of the Pygmy Dipole Resonance in 208 Pb

New experimental data on the neutron single-particle character of the Pygmy Dipole Resonance (PDR) in Pb 208 are presented. They were obtained from ( d , p ) and resonant proton scattering experiments performed at the Q3D spectrograph of the Maier-Leibnitz Laboratory in Garching, Germany. The new data are compared to the large suite of complementary, experimental data available for Pb 208 and establish ( d , p ) as an additional, valuable, experimental probe to study the PDR and its collectivity. Furthermore, besides the single-particle character of the states, different features of the strength distributions are discussed and compared to large-scale shell model (LSSM) and energy-density functional plus quasiparticle-phonon model theoretical approaches to elucidate the microscopic structure of the PDR in Pb 208 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On the impact of the galaxy window function on cosmological parameter estimation

One important source of systematics in galaxy redshift surveys comes from the estimation of the galaxy window function. Up until now, the impact of the uncertainty in estimating the galaxy window function on parameter inference has not been properly studied. In this paper, we show that the uncertainty and the bias in estimating the galaxy window function will be salient for ongoing and next-generation galaxy surveys using a simulation-based approach. With a specific case study of cross-correlating emission-line galaxies from the DESI Legacy Imaging Surveys and the Planck cosmic microwave background lensing map, we show that neural network-based regression approaches to modelling the window function are superior in comparison to linear regression-based models. Here, we additionally show that the definition of the galaxy overdensity estimator can impact the overall signal-to-noise of observed power spectra. Finally, we show that the additive biases coming from the window functions can significantly bias the modes of the inferred parameters and also degrade their precision. Thus, a careful understanding of the window functions will be essential to conduct cosmological experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Capturing the Coverage Dependence of Aromatics’ Adsorption through Mean-Field Models

To capture the dominant interactions (surface-mediated and through-space steric) in catalytic hydrodeoxygenation systems, coverage-dependent mean-field models of aromatic adsorption are developed on Pt(111) and Ru(0001). We derive three key insights from this work: (1) we can universally apply mean-field models to capture the coverage-dependent behavior of oxygenated aromatics on transition metal surfaces, (2) we can deconvolute surface-mediated and throughspace steric interactions from the mean-field model, and (3) we can develop relatively accurate models that predict the adsorption energy of aromatics on transition metal surfaces for the full coverage range using the work function at lowest modeled coverage. In conclusion, our approach enables the rapid prediction of coverage-dependent behavior of oxygenated aromatics on transition metal surfaces, reducing the computational cost associated with these studies by an order of magnitude.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Parametric matrix models

We present a general class of machine learning algorithms called parametric matrix models. In contrast with most existing machine learning models that imitate the biology of neurons, parametric matrix models use matrix equations that emulate physical systems. Similar to how physics problems are usually solved, parametric matrix models learn the governing equations that lead to the desired outputs. Parametric matrix models can be efficiently trained from empirical data, and the equations may use algebraic, differential, or integral relations. While originally designed for scientific computing, we prove that parametric matrix models are universal function approximators that can be applied to general machine learning problems. After introducing the underlying theory, we apply parametric matrix models to a series of different challenges that show their performance for a wide range of problems. For all the challenges tested here, parametric matrix models produce accurate results within an efficient and interpretable computational framework that allows for input feature extrapolation.

Computational science↗

Reviewing MACCS Capabilities for Assessing Tritium Releases to the Environment

Tritium has a unique physical and chemical behavior which causes it to be highly mobile in the environment. As it behaves similarly to hydrogen in the environment, it may also be readily incorporated into the water cycle and other biological processes. These factors and other environmental transformations may also cause the oxidation of an elemental tritium release, resulting in a multiple order of magnitude increase in dose coefficient and radiotoxicity. While source term development and understanding for advanced reactors are still underway, tritium may be a radionuclide of interest. It is thus important to understand how tritium moves through the environment and how the MACCS accident consequence code handles acute tritium releases in an accident scenario. Additionally, existing tritium models may have functionalities that could inform updates to MACCS to handle tritium. In this report tritium transport is reviewed and existing tritium models are summarized in view of potential updates to MACCS.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Spatial Correlations of the Poisson Model for Radiation Transport

Characterizing the relationship between bulk physical properties and mixing in randomly heterogeneous media is a central challenge across many areas of science and engineering. A benchmark model for such studies is the Poisson model, a random tessellation of space by a Poisson process of hyperplanes. In radiation transport studies, the lack of exact expressions for the Poisson model’s spatial multipoint functions has led to approximate methods being used, introducing unquantified sources of error. Here, we recently introduced an exact solution for the Poisson model’s multipoint functions and closely related conditional probability functions (CPFs), providing a new opportunity to understand and reduce these sources of error. In this paper, we enable a more rigorous investigation of radiation transport in stochastic media by applying the recently introduced exact solution for the Poisson model’s CPFs. This paper consists of three main contributions. First, we introduce a unified framework for CPFs of the Poisson model, encompassing the recently introduced exact CPFs as well as the previously introduced atomic mix, nearest-neighbor, and combination CPFs. This framework also includes existing pruning techniques for the approximate CPFs, such as angular exclusion, as well as a novel form of angular exclusion suitable for the exact CPFs. Second, we use the exact CPFs to characterize the spatial regions where each approximate three-point CPF is most accurate, thereby explaining the observed hierarchy of accuracy among the approximate models. Finally, we evaluate material transmittance, reflectance, and flux in a three-dimensional test problem using conditional point sampling, demonstrating the relationship between CPF accuracy and transport simulation accuracy.

Poisson model↗

Negative components in the representation $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ of heavy-ion fusion cross sections

Negative components are observed in the representation $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ in heavy-ion fusions related to recent studies of noticeable oscillations in the high-energy region of the fusion excitation functions. A modified multiGaussian model has been developed, which includes contributions for both positive and negative components in the $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ representation. The model reproduces the heavy-ion fusion excitation functions very well, including oscillations. The property of these negative components in the $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ spectrum and the reaction mechanism of fusion at energies above the Coulomb barrier are investigated. Finally, the model also reproduces well light heavy-ion fusion reactions between 12 C and 16 O at high energy; results imply that the compound nucleus channel effect may have an important contribution at the high energy region for all fusion systems.

low & intermediate energy heavy-ion reactions↗

Accurate simulation of surfaces and interfaces of ten FCC metals and steel using Lennard–Jones potentials

Abstract The earlier integration of validated Lennard–Jones (LJ) potentials for 8 fcc metals into materials and biomolecular force fields has advanced multiple research fields, for example, metal–electrolyte interfaces, recognition of biomolecules, colloidal assembly of metal nanostructures, alloys, and catalysis. Here we introduce 12-6 and 9-6 LJ parameters for classical all-atom simulations of 10 further fcc metals (Ac, Ca (α), Ce (γ), Es (β), Fe (γ), Ir, Rh, Sr (α), Th (α), Yb (β)) and stainless steel. The parameters reproduce lattice constants, surface energies, water interfacial energies, and interactions with (bio)organic molecules in 0.1 to 5% agreement with experiment, as well as qualitative mechanical properties under standard conditions. Deviations are reduced up to a factor of one hundred in comparison to earlier Lennard–Jones parameters, embedded atom models, and density functional theory. We also explain a quantitative correlation between atomization energies from experiments and surface energies that supports parameter development. The models are computationally very efficient and applicable to an exponential space of alloys. Compatibility with a wide range of force fields such as the Interface force field (IFF), AMBER, CHARMM, COMPASS, CVFF, DREIDING, OPLS-AA, and PCFF enables reliable simulations of nanostructures up to millions of atoms and microsecond time scales. User-friendly model building and input generation are available in the CHARMM-GUI Nanomaterial Modeler. As a limitation, deviations in mechanical properties vary and are comparable to DFT methods. We discuss the incorporation of reactivity and features of the electronic structure to expand the range of applications and further increase the accuracy.

Kanhaiya, Krishan (ORCID:0000000236222655)↗