Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model exploration”

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 649 records · Page 36

PSA 2025 Presentation: "Modeling and Sensitivity Analysis of a Generation IV Pebble Bed Reactor Using MELCOR 2.2"

Accompanying the advancement of reactor technologies is the need for computational modeling and simulation to predict their behavior under normal operating conditions and accident scenarios. New Generation IV reactor designs which employ non-conventional fuel have a particular need for modeling the behavior and release of radionuclides and other material from the fuel. In this work, MELCOR version 2.2, a system-level safety and accident scenario code developed by Sandia National Laboratories, was used to model a 200-MWth pebble bed modular reactor and calculate the inventories of circulating and deposited graphite, metal dust, and elemental components released from the fuel elements. A base case modeling the reactor under standard operating conditions was calculated using MELCOR and the inventories were extrapolated to 30 years of operation time using a logarithmic regression fit. A sensitivity analysis was also performed in which several key parameters for the base case model were modified to explore the effect of these changes on the inventories calculated by MELCOR. A set of transient scenario simulations for a depressurized loss of forced cooling (DLOFC) accident were also performed. The results of the sensitivity analysis and transient simulations are reported and discussed in relation to the modeling techniques used for this study.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Kinetics Modeling and Reactor Design Study of Glucose-to-Terpenes Cell-Free Conversion

Cell-free systems offer many advantages over traditional biological conversion by eliminating biological growth constraints. It also offers easy manipulation and finetuning of the reaction conditions for each individual enzyme. The conversion of cellulosic glucose to Limonene, a terpene, is a promising pathway for producing fuels and chemicals. Recent advances in developing cell-free systems focuses on bench scale optimization of terpene yield and to demonstrate its feasibility towards commercialization [1,2]. There is significant knowledge gap regarding reaction kinetics of these cell-free systems to further study how it will perform at larger scale. We present here, our studies on reaction kinetics and reactor design implications of cell-free glucose to Limonene conversion to facilitate the further development and commercialization of this process. We developed a novel kinetic model based on the metabolic-network structure of the cell-free system with multi-substrate reversible Michaelis-Menten rate law. To estimate kinetic parameters for this system of rate equations, we employed Bayesian optimization to perform global search with the assistance of gaussian processes to balance exploration and exploitation. The model parameters estimated showed good results compared with experimental data. The estimated parameters were used to perform sensitivity analysis. We found that Hexokinase is one of the most critical enzymes that affect the conversion of the glucose. We also observed that abundance of co-factors is also critical to the conversion of glucose to limonene. We investigated packed bed reactors with enzymes immobilized on the surface of particles to convert glucose stream into Limonene for larger scale production. The reactor design such as particle size, enzyme loading, and flow rate are found to be critical for improving yields. [1] Dudley, Q.M., Nash, C.J. and Jewett, M.C., 2019. Synthetic Biology, 4(1), p.ysz003. [2] Korman, T.P., Opgenorth, P.H. and Bowie, J.U., 2017. Nature communications, 8(1), p.15526.

09 BIOMASS FUELS↗

Model-Independent Searches for New Physics in Multi-Body Invariant Masses

Model-independent searches for physics beyond the Standard Model typically focus on invariant masses of two objects (jets, leptons or photons). In this study, we explore opportunities for similar model-agnostic searches in multi-body invariant masses. In particular, we focus on the situations in which new physics can be observed in a model-independent way in three and four-body invariant masses of jets and leptons. Such searches may have good prospects in finding new physics in the situations when two-body invariant masses, which have been extensively explored at collider experiments in the past, cannot provide sufficient signatures for experimental observations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Co-Active Subspace Methods for the Joint Analysis of Adjacent Computer Models

Active subspace (AS) methods are a valuable tool for understanding the relationship between the inputs and outputs of a Physics simulation. In this article, an elegant generalization of the traditional ASM is developed to assess the co-activity of two computer models. This generalization, which we refer to as a Co-Active Subspace (Co-AS) Method, allows for the joint analysis of two or more computer models allowing for thorough exploration of the alignment (or non-alignment) of the respective gradient spaces. We define co-active directions, co-sensitivity indices, and a scalar “concordance” metric (and complementary “discordance” pseudo-metric) and we demonstrate that these are powerful tools for understanding the behavior of a class of computer models, especially when used to supplement traditional AS analysis. Details for efficient estimation of the Co-AS and an accompanying R package (concordance) are provided. Practical application is demonstrated through analyzing a set of simulated rate stick experiments for PBX 9501, a high explosive, offering insights into complex model dynamics.

97 MATHEMATICS AND COMPUTING↗

Statistical analysis and degradation pathway modeling of photovoltaic minimodules with varied packaging strategies

Degradation pathway models constructed using network structural equation modeling (netSEM) are used to study degradation modes and pathways active in photovoltaic (PV) system variants in exposure conditions of high humidity and temperature. This data-driven modeling technique enables the exploration of simultaneous pairwise and multiple regression relationships between variables in which several degradation modes are active in specific variants and exposure conditions. Durable and degrading variants are identified from the netSEM degradation mechanisms and pathways, along with potential ways to mitigate these pathways. A combination of domain knowledge and netSEM modeling shows that corrosion is the primary cause of the power loss in these glass/backsheet PV minimodules. We show successful implementation of netSEM to elucidate the relationships between variables in PV systems and predict a specific service lifetime. The results from pairwise relationships and multiple regression show consistency. This work presents a greater opportunity to be expanded to other materials systems.

electrical measurements↗

Small-Signal Stability Impacts of Load and Network Dynamics on Grid-Forming Inverters

This paper presents several stability analyses for grid-forming inverters and synchronous generators considering the dynamics of transmission lines and different load models. Load models are usually of secondary importance compared to generation source models, but as the results show, they play a crucial role in stability studies with the introduction of inverter-based resources. Given inverter control time scales, the implications of considering or neglecting electromagnetic transients of the network are very relevant in the stability assessments. In this paper, we perform eigenvalue analyses for inverter-based resources and synchronous machines connected to a load and explore the effects of multiples models under different network representations. We explore maximum loadability of inverter-based resources and synchronous machines, while analyzing the effects of load and network dynamic models on small-signal stability. The results show that the network representation plays a fundamental role in the stability of the system of different load models. The resulting stability regions are significantly different depending on the source and load model considered.

grid forming inverters↗

Machine-learning enabled thermodynamic model for the design of new rare-earth compounds

We employ a descriptor based machine-learning approach to assess the effect of chemical alloying on formation-enthalpy of rare-earth intermetallics. Application of machine-learning approaches in rare-earth intermetallic design have been sparse due to limited availability of reliable datasets. In this work, we developed an ‘in-house’ rare-earth database with more than 600 + compounds, each entry was populated with formation enthalpy and related atomic features using high-throughput density-functional theory (DFT). The SISSO (sure independence screening and sparsifying operator) based machine-learning method with meaningful atomic features was used for training and testing the formation enthalpies of rare earth compounds. The complex lattice function coupled with the machine-learning model was used to explore the effect of transition metal alloying on the energy stability of Ce based cubic Laves phases (MgCu 2 type). The SISSO predictions show good agreement with high-fidelity DFT calculations and X-ray powder diffraction measurements. Our study provides quantitative guidance for compositional considerations within a machine-learning model and discovering new metastable materials. The electronic-structure of Ce-Fe-Cu based compound was also analyzed to get an in-depth understanding of the electronic origin of phase stability. The interpretable analytical models in combination with density-functional theory and experiments provide a fast and reliable design guide for discovering technologically useful materials.

36 MATERIALS SCIENCE↗

Reduced-order modeling for efficient cross section library development in high-temperature gas reactor pebble-bed depletion analysis

Accurate modeling of running-in and equilibrium conditions in pebble-bed reactors (PBRs) requires precise microscopic multigroup neutron cross sections. In Griffin, deterministic neutronics calculations rely on multivariate interpolation over large cross section libraries, resulting in significant memory usage and performance bottlenecks. This work, together with a companion paper on Griffin integration, explores reduced-order models (ROMs) to replace interpolation with lightweight surrogates. Several ROM techniques are benchmarked, with deep neural networks (DNNs) demonstrating superior memory efficiency, scalability, and predictive accuracy. A total of 295 DNNs were trained to build a comprehensive isotope library, integrated into Griffin through a custom LibTorch interface for depletion analysis. Initial results demonstrate that DNN-based ROMs drastically reduce memory demands while preserving accuracy, enabling finer tabulations and additional state variables without overhead. In conclusion, the framework also supports online cross section generation and real-time DNN updates through transfer learning, improving fidelity by capturing self-shielding and evolving nuclide compositions during burnup.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Feasibility of fusion plasma burn control via real-time, sub-divertor neutral gas isotopic and compositional analysis

The ability to provide fusion burn control without requiring physical access through the first wall and fuel breeding blankets, would be vital for any future, magnetically confined fusion power reactor. A multi-sensor, fusion fuel cycle exhaust, neutral gas analysis system on JET, capable of delivering real time data, and accessing only the sub-divertor region, provides an excellent example of such capability. Optimized for and operated during the deuterium–tritium experimental campaigns 2 and 3 (DTE2, DTE3), it is proving valuable for planning to explore fusion reactor burn control in ITER with a comparable diagnostic system called the Diagnostic Residual Gas Analyzer (DRGA). This paper aims to show feasibility of developing model-based controllers for ITER and next generation, reactor-relevant devices, by building both on the empirical experience in JET-DTE2, and on the already emerging experience on developing such models specifically for ITER. The paper begins with a specific use-case from JET-DTE2, pertaining to the observed sensitivity of the fusion neutron yield on the concentration of isotopic helium-3 ( 3 He), with data from one of the high-performance DT shots exhibited with emphasis on the 3 He measurement via the sub-divertor. Then, a first model is developed and then explored with simulations that aim to discover how well the controllers in the model react to either insufficient levels of 3 He or excessive levels of 3 He. The simulations then explore potential impact from a delay in the measurement (or the response) that would be comparable to the ∼1 s, conductance limited response for the ITER DRGA system, currently in its final design. The simulations show that control is feasible, and that its effectiveness is not significantly impacted by such delay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Programmable exploration of magnetic states in Lieb-kagome interpolated lattices

We investigate a hybrid modeling framework in which a quantum annealer is used to simulate magnetic interactions in molecular qubit lattices inspired by experimentally realizable systems. Using phthalocyanine assemblies as a structurally constrained prototype, we model a continuous deformation from a Lieb to a kagome lattice, revealing frustration-driven disorder and magnetic field-induced reordering in the spin structure. Here, the goal is to show how a quantum annealer can operate as a physically instantiated, programmable platform to emulate experimentally relevant lattice deformations and produce observables in a manner analogous to an experimental measurement, enabling the characterization of magnetic arrangements beyond the reach of current molecular architectures. This surrogate modeling approach offers a pathway to explore and iteratively design tunable magnetic states in synthetic materials. The synthetic design, structural characterization, and quantum simulation framework established here defines a modular and scalable paradigm for probing the limits of engineered matter across chemistry, condensed matter, and quantum information science.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Quasilinear modelling of collisional trapped electron modes

The simplification of collision operators is necessary for quasilinear turbulence modelling used with integrated modelling frameworks, such as the gyrokinetic code QuaLiKiz. The treatment of collisions greatly impacts the accuracy of trapped electron mode (TEM) modelling, which is necessary to predict the electron heat flux and the balance between inward and outward particle fluxes. In particular, accurate particle flux predictions are necessary to successfully model density peaking in the tokamak core. We explored two ways of improving collisional TEM model reduction for tokamak core plasmas. First, we carried out linear GENE simulations to study the complex interplay between collisions and trapped electrons. We then used these simulations to define an effective trapped fraction to characterize the collisional TEM based on two key parameters, the local inverse-aspect ratio $\epsilon$ and the collisionality $\nu ^\ast$. One aspect missing from analytical TEM research is that the collisional frequency and the bounce-transit frequency are both velocity dependent; this effective trapped fraction takes both into account. In doing so, we determined that two parameters are not enough to model the collisional TEM, as an additional third free parameter was necessary. We determined that this model, as currently formulated, is not suitable for integrated modelling purposes. Second, we directly improved QuaLiKiz's Krook operator, which relies on two free parameters. We determined that these parameters required adjustments against higher-fidelity collisional models. In order to improve density profile predictions when paired with integrated models, we refined the Krook operator by using GENE simulations as a higher-fidelity point of comparison. We then demonstrate strong improvement of density peaking predictions of QuaLiKiz within the integrated modelling framework JETTO.

fusion plasma↗

Quantitative modeling of carcinogenesis induced by single beams or mixtures of space radiations using targeted and non-targeted effects

Abstract Ionizing radiations encountered by astronauts on deep space missions produce biological damage by two main mechanisms: (1) Targeted effects (TE) due to direct traversals of cells by ionizing tracks. (2) Non-targeted effects (NTE) caused by release of signals from directly hit cells. The combination of these mechanisms generates non-linear dose response shapes, which need to be modeled quantitatively to predict health risks from space exploration. Here we used a TE + NTE model to analyze data on APC (1638N/+) mouse tumorigenesis induced by space-relevant doses of protons, 4 He, 12 C, 16 O, 28 Si or 56 Fe ions, or γ rays. A customized weighted Negative Binomial distribution was used to describe the radiation type- and dose-dependent data variability. This approach allowed detailed quantification of dose–response shapes, NTE- and TE-related model parameters, and radiation quality metrics (relative biological effectiveness, RBE, and radiation effects ratio, RER, relative to γ rays) for each radiation type. Based on the modeled responses for each radiation type, we predicted the tumor yield for a Mars-mission-relevant mixture of these radiations, using the recently-developed incremental effect additivity (IEA) synergy theory. The proposed modeling approach can enhance current knowledge about quantification of space radiation quality effects, dose response shapes, and ultimately the health risks for astronauts.

60 APPLIED LIFE SCIENCES↗

A More La Niña–Like Response to Radiative Forcing after Flux Adjustment in CESM2

Abstract In response to greenhouse gas forcing, most coupled global climate models project the tropical Pacific SST trend toward an “El Niño–like” state, with a reduced zonal SST gradient and a weakened Walker circulation. However, observations over the last five decades reveal a trend toward a more “La Niña–like” state with a strengthening zonal SST gradient. Recent research indicates that the identified trend differences are unlikely to be entirely due to internal variability and probably result, at least in part, from systematic model biases. In this study, Community Earth System Model, version 2 (CESM2), is used to explore how mean-state biases within the model may influence its forced response to radiative forcing in the tropical Pacific. The results show that using flux adjustment to reduce the mean-state bias in CESM2 over the tropical regions results in a more La Niña–like trend pattern in the tropical Pacific, with a strengthening of the tropical Pacific zonal SST gradient and a relatively enhanced Walker circulation, as hypothesized to occur if the ocean thermostat mechanism is stronger than the atmospheric mechanisms which by themselves would weaken the Walker circulation. We also find that the historical strengthening of the tropical Pacific zonal gradient is transient but persists into the near term in a high-emissions future warming scenario. These results suggest the potential of flux adjustment as a method for developing alternative projections that represent a wider range of possible future tropical Pacific warming scenarios, especially for a better understanding of regional patterns of climate risk in the near term.

Zhuo, Jing-Yi [Lamont-Doherty Earth Observatory, C↗

Toward accurate prediction of partial-penetration laser weld performance informed by three-dimensional characterization – Part II: μCT based finite element simulations

The mechanical behavior of partial-penetration laser welds exhibits significant variability in engineering quantities such as strength and apparent ductility. Understanding the root cause of this variability is important when using such welds in engineering designs. In Part II of this work, we develop finite element simulations with geometry derived from micro-computed tomography (μCT) scans of partial-penetration 304L stainless steel laser welds that were analyzed in Part I. We use these models to study the effects of the welds’ small-scale geometry, including porosity and weld depth variability, on the structural performance metrics of weld ductility and strength under quasi-static tensile loading. We show that this small-scale geometry is the primary cause of the observed variability for these mechanical response quantities. Additionally, we explore the sensitivity of model results to the conversion of the μCT data to discretized model geometry using different segmentation algorithms, and to the effect of small-scale geometry simplifications for pore shape and weld root texture. The modeling approach outlined and results of this work may be applicable to other material systems with small-scale geometric features and defects, such as additively manufactured materials.

36 MATERIALS SCIENCE↗

LigninGraphs: lignin structure determination with multiscale graph modeling

Lignin is an aromatic biopolymer found in ubiquitous sources of woody biomass. Designing and optimizing lignin valorization processes requires a fundamental understanding of lignin structures. Experimental characterization techniques, such as 2D-heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectra, could elucidate the global properties of the polymer molecules. Computer models could extend the resolution of experiments by representing structures at the molecular and atomistic scales. We introduce a graph-based multiscale modeling framework for lignin structure generation and visualization. The framework employs accelerated rejection-free polymerization and hierarchical Metropolis Monte Carlo optimization algorithms. We obtain structure libraries for various lignin feedstocks based on literature and new experimental NMR data for poplar wood, pinewood, and herbaceous lignin. The framework could guide researchers towards feasible lignin structures, efficient space exploration, and future kinetics modeling. Its software implementation in Python, LigninGraphs, is open-source and available on GitHub.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fast Approximate Methods for Modified Gravity Cosmological Simulations

The accelerated expansion of the Universe poses a major theoretical puzzle. Although the assumption of a non-zero cosmological constant provides a minimal extension of general relativity that is consistent with observational data, many theories of modified gravity have been suggested as possible alternatives. Predictions of structure formation for these models in the fully nonlinear regime are very expensive and it is difficult, if not impossible, to explore such a huge space of models and parameters using high-resolution N-body simulations. Even in the mildly nonlinear regime, perturbative methods can become extremely complex. We explore whether simplified dynamical approximations, applicable for a certain set of cosmological probes, can be used to investigate models of modified gravity with acceptable accuracy in the latter instance. For the case of chameleon gravity, we find that these methods can indeed be used to explore the region around the baryon acoustic oscillation scale, k ~ 0:1 hMpc-1 but not much further.

cosmology↗

Quantum model learning agent: characterisation of quantum systems through machine learning

Accurate models of real quantum systems are important for investigating their behaviour, yet are difficult to distil empirically. Here, we report an algorithm—the quantum model learning agent (QMLA)—to reverse engineer Hamiltonian descriptions of a target system. We test the performance of QMLA on a number of simulated experiments, demonstrating several mechanisms for the design of candidate Hamiltonian models and simultaneously entertaining numerous hypotheses about the nature of the physical interactions governing the system under study. QMLA is shown to identify the true model in the majority of instances, when provided with limited a priori information, and control of the experimental setup. Our protocol can explore Ising, Heisenberg and Hubbard families of models in parallel, reliably identifying the family which best describes the system dynamics. We demonstrate QMLA operating on large model spaces by incorporating a genetic algorithm to formulate new hypothetical models. The selection of models whose features propagate to the next generation is based upon an objective function inspired by the Elo rating scheme, typically used to rate competitors in games such as chess and football. In all instances, our protocol finds models that exhibit F 1 score ≥ 0.88 when compared with the true model, and it precisely identifies the true model in 72% of cases, whilst exploring a space of over 250 000 potential models. By testing which interactions actually occur in the target system, QMLA is a viable tool for both the exploration of fundamental physics and the characterisation and calibration of quantum devices.

97 MATHEMATICS AND COMPUTING↗

Automated electrosynthesis reaction mining with multimodal large language models (MLLMs)

Leveraging the chemical data available in legacy formats such as publications and patents is a significant challenge for the community. Automated reaction mining offers a promising solution to unleash this knowledge into a learnable digital form and therefore help expedite materials and reaction discovery. However, existing reaction mining toolkits are limited to single input modalities (text or images) and cannot effectively integrate heterogeneous data that is scattered across text, tables, and figures. In this work, we go beyond single input modalities and explore multimodal large language models (MLLMs) for the analysis of diverse data inputs for automated electrosynthesis reaction mining. We compiled a test dataset of 65 articles (MERMES-T24 set) and employed it to benchmark five prominent MLLMs against two critical tasks: (i) reaction diagram parsing and (ii) resolving cross-modality data interdependencies. The frontrunner MLLM achieved ≥96% accuracy in both tasks, with the strategic integration of single-shot visual prompts and image pre-processing techniques. We integrate this capability into a toolkit named MERMES (multimodal reaction mining pipeline for electrosynthesis). Our toolkit functions as an end-to-end MLLM-powered pipeline that integrates article retrieval, information extraction and multimodal analysis for streamlining and automating knowledge extraction. This work lays the groundwork for the increased utilization of MLLMs to accelerate the digitization of chemistry knowledge for data-driven research.

Leong, Shi Xuan↗