Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation analytics”

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 973 records · Page 54

Effect of Solvent on the Interaction of Lignin with a Zeolite Nanosheet in the Condensed Phase

Lignin is the essential building block of lignocellulosic biomass, an excellent renewable source of different aromatic monomers for the polymer and biofuel industry. The depolymerization of lignin into value-added chemicals and fuels through the catalytic process poses a significant challenge due to the complex structure of lignin. Understanding lignin’s conformational diversity and dynamics in the liquid phase is crucial to designing an effective depolymerization process. Here, we conducted all-atom molecular dynamics simulations to understand the conformation and dynamics of softwood lignin on the all-silica zeolite nanosheet based on the MFI topology in a binary mixture of water–methanol at three different molar compositions (0%, 50%, and 100% methanol). We observed that the methanol–surface interaction is stronger than the water–surface interaction, and methanol readily diffused into the MFI core. Lignin surface contacts decrease with increasing methanol composition due to higher solubility and dynamics. Lignin dynamics on the surface in neat water is an order of magnitude smaller than methanol. Here, we also found that lignin adopts a slightly extended conformation when it stays on the surface than in the bulk solution phase for the pure water case, whereas for pure methanol and the binary solution structures are statistically similar.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First principles reactive simulation for equation of state prediction

The high cost of density functional theory (DFT) has hitherto limited the ab initio prediction of the equation of state (EOS). In this article, we employ a combination of large scale computing, advanced simulation techniques, and smart data science strategies to provide an unprecedented ab initio performance analysis of the high explosive pentaerythritol tetranitrate (PETN). Comparison to both experiment and thermochemical predictions reveals important quantitative limitations of DFT for EOS prediction and thus the assessment of high explosives. In particular, we find that DFT predicts the energy of PETN detonation products to be systematically too high relative to the unreacted neat crystalline material, resulting in an underprediction of the detonation velocity, pressure, and temperature at the Chapman–Jouguet state. The energetic bias can be partially accounted for by high-level electronic structure calculations of the product molecules. Furthermore we demonstrate a modeling strategy for mapping chemical composition across a wide parameter space with limited numerical data, the results of which suggest additional molecular species to consider in thermochemical modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dissolution of a cylindrical disk in Hele-Shaw flow: a conformal-mapping approach

Here we apply conformal mapping to find the evolving shapes of a dissolving cylinder in a potential flow. Similar equations can be used to describe melting in a flowing liquid phase. Results are compared with microfluidic experiments and numerical simulations. Shapes predicted by conformal mapping agree almost perfectly with experimental observations, after a modest (20 %) rescaling of the time. Finite-volume simulations show that the differences with experiment are connected to the underlying assumptions of the analytical model: potential flow and diffusion-limited dissolution. Approximate solutions of the equations describing the evolution of the shape of the undissolved solid can be derived from a Laurent expansion of the mapping function from the unit circle. Asymptotic expressions for the evolution of the area of the disk and the shift in its centre of mass have been derived at low and high Péclet number. Analytic approximations to the leading-order Laurent coefficients provide additional insight into the mechanisms underlying pore-scale dissolution.

58 GEOSCIENCES↗

Transferable Water Potentials Using Equivariant Neural Networks

Machine learning interatomic potentials (MLIPs) have emerged as a technique that promises quantum theory accuracy for reduced cost. It has been proposed [J. Chem. Phys. 2023, 158, 084111] that MLIPs trained on solely liquid water data cannot accurately transfer to the vapor–liquid equilibrium while recovering the many-body decomposition (MBD) analysis of gas-phase water clusters. This suggests that MLIPs do not directly learn the physically correct interactions of water molecules, limiting transferability. In this work, we show that MLIPs using equivariant architecture and trained on 3200 liquid water structures reproduces liquid-phase water properties (e.g., density within 0.003 g/cm 3 between 230 and 365 K), vapor–liquid equilibrium properties up to 550 K, the MBD analysis of gas-phase water cluster up to six-body interactions, and the relative energy and the vibrational density of states of ice phases. We show that potentials developed using equivariant MLIPs allow transferability for arbitrary phases of water that remain stable in nanosecond long simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Wettability of ultra-small pores of carbon electrodes by size-asymmetric ionic fluids

Recently, we studied the phase behavior of ionic fluids under confinement using the classical density functional theory within the framework of the restricted primitive model. Additionally, the theoretical results indicate that narrowing the pore size may lead to a drastic reduction in the electric double layer capacitance, while increasing the surface electrical potential would improve the ionic accessibility of micropores. In this work, we extend the theoretical investigation to systems containing size-asymmetric electrolytes that may exhibit a vapor-liquid like phase transition in the bulk phase. The effects of pore size and surface electric potential on the phase diagram and microscopic structures of the confined electrolytes were studied over a broad range of parameters. We found that decreasing the pore size or increasing the surface potential could destabilize the liquid phase in micropores, and capillary evaporation could occur regardless of the size asymmetry between cations and anions. Compared to that in a symmetric ionic system, the vapor-liquid phase separation is more likely to take place as the size asymmetry becomes more pronounced. The phase transition would alter the “accessibility” of ions to micropores and lead to coexisting micropores with different surface charge densities as identified by Monte Carlo simulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laser polarization effects on crossed-beam energy transfer in inertial confinement fusion

Crossed-beam energy transfer (CBET) is a consequential process in inertial confinement fusion (ICF) experiments. CBET depends on the polarization state of the interacting lasers and can, in turn, alter their polarization. Here, we derive analytical expressions for two-beam CBET with linearly and circularly polarized light, and present simulations of CBET that include polarization effects for ICF conditions relevant to experiments on the National Ignition Facility. The results show that CBET mixes the polarization of the beams as they propagate through the target and leads to azimuthal power imbalances between beams within a given cone (i.e., at the same polar angle in the target chamber). These variations are shown to correlate with variations in measured backscatter between beams of a same quadruplet. Using circularly polarized, instead of linearly polarized, light produces similar overall levels of CBET within a cone of beams but tends to reduce the azimuthal variations. This is expected to slightly improve irradiation symmetry and reduce the risk of backscatter from stimulated Brillouin scattering. Circular polarization may therefore be beneficial for the next generation of laser fusion drivers.

Physics - Plasma physics↗

Scale-free, programmable design of morphable chain loops of kilobots and colloidal motors

Micron-scale robots require systems that can morph into arbitrary target configurations controlled by external agents such as heat, light, electricity, and chemical environment. Achieving this behavior using conventional approaches is challenging because the available materials at these scales are not programmable like their macroscopic counterparts. To overcome this challenge, we propose a design strategy to make a robotic machine that is both programmable and compatible with colloidal-scale physics. Our strategy uses motors in the form of active colloidal particles that constantly propel forward. We sequence these motors end-to-end in a closed chain forming a two-dimensional loop that folds under its mechanical constraints. We encode the target loop shape and its motion by regulating six design parameters, each scale-invariant and achievable at the colloidal scale. We demonstrate the plausibility of our design strategy using centimeter-scale robots called kilobots . We use Brownian dynamics simulation to explore the large design space beyond that possible with kilobots, and present an analytical theory to aid the design process. Multiple loops can also be fused together to achieve several complex shapes and robotic behaviors, demonstrated by folding a letter shape “M,” a dynamic gripper , and a dynamic pacman . The material-agnostic, scale-free, and programmable nature of our design enables building a variety of reconfigurable and autonomous robots at both colloidal scales and macroscales.

36 MATERIALS SCIENCE↗

Personalized Tucker Decomposition: Modeling Commonality and Peculiarity on Tensor Data

In this paper, we propose a personalized Tucker decomposition (perTucker) to address the limitations of traditional tensor decomposition methods in capturing heterogeneity across different datasets. perTucker decomposes tensor data into shared global components and personalized local components. We introduce an order orthogonality assumption and develop a proximal gradient regularized block coordinate descent algorithm guaranteed to converge to a stationary point. The unique and common representations learned by perTucker reveal intrinsic statistical patterns in data and provide valuable information for a wide range of downstream analytics, including anomaly detection, source classification, and clustering. We demonstrate perTucker’s effectiveness through a simulation study and two case studies on solar flare detection and tonnage signal classification.

14 SOLAR ENERGY↗

Validation of semi-analytical, semi-empirical covariance matrices for two-point correlation function for early DESI data

ABSTRACT We present an extended validation of semi-analytical, semi-empirical covariance matrices for the two-point correlation function (2PCF) on simulated catalogs representative of luminous red galaxies (LRGs) data collected during the initial 2 months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). We run the pipeline on multiple effective Zel’dovich (EZ) mock galaxy catalogs with the corresponding cuts applied and compare the results with the mock sample covariance to assess the accuracy and its fluctuations. We propose an extension of the previously developed formalism for catalogs processed with standard reconstruction algorithms. We consider methods for comparing covariance matrices in detail, highlighting their interpretation and statistical properties caused by sample variance, in particular, non-trivial expectation values of certain metrics even when the external covariance estimate is perfect. With improved mocks and validation techniques, we confirm a good agreement between our predictions and sample covariance. This allows one to generate covariance matrices for comparable data sets without the need to create numerous mock galaxy catalogs with matching clustering, only requiring 2PCF measurements from the data itself. The code used in this paper is publicly available at https://github.com/oliverphilcox/RascalC.

79 ASTRONOMY AND ASTROPHYSICS↗

Qubit control noise spectroscopy with optimal suppression of dephasing

Here, we extend quantum noise spectroscopy (QNS) of amplitude control noise to settings where dephasing noise or detuning errors make significant contributions to qubit dynamics. Previous approaches to characterize amplitude noise are limited by their vulnerability to low-frequency dephasing noise and static detuning errors, which can overwhelm the target control noise signal and introduce bias into estimates of the amplitude noise spectrum. To overcome this problem, we leverage optimal control to identify a family of amplitude control waveforms that optimally suppress low-frequency dephasing noise and detuning errors, while maintaining the spectral concentration in the amplitude filter essential for spectral estimation. The waveforms found via numerical optimization have surprisingly simple analytic forms, consisting of oscillating sine waves obeying particular amplitude and frequency constraints. In numerically simulated QNS experiments, these waveforms demonstrate superior robustness, enabling accurate estimation of the amplitude noise spectrum in regimes where existing approaches are biased by low-frequency dephasing noise and detuning errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Many-Body Level Statistics of Single-Particle Quantum Chaos

We consider a noninteracting many-fermion system populating levels of a unitary random matrix ensemble (equivalent to the q = 2 complex Sachdev-Ye-Kitaev model)—a generic model of single-particle quantum chaos. We study the corresponding many-particle level statistics by calculating the spectral form factor analytically using algebraic methods of random matrix theory, and match it with an exact numerical simulation. Despite the integrability of the theory, the many-body spectral rigidity is found to have a surprisingly rich landscape. In particular, we find a residual repulsion of distant many-body levels stemming from single-particle chaos, together with islands of level attraction. These results are encoded in an exponential ramp in the spectral form factor, which we show to be a universal feature of nonergodic many-fermion systems embedded in a chaotic medium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Theory of Defect-Induced Crystal Field Perturbations in Rare-Earth Magnets

We present a theory describing the single-ion anisotropy of rare-earth (RE) magnets in the presence of point defects. Taking the RE-lean 1 ∶ 12 magnet class as a prototype, we use first-principles calculations to show how the introduction of Ti substitutions into SmFe 12 perturbs the crystal field, generating new coefficients due to the lower symmetry of the RE environment. We then demonstrate that these perturbations can be described extremely efficiently using a screened point charge model. We provide analytical expressions for the anisotropy energy that can be straightforwardly implemented in atomistic spin dynamics simulations, meaning that such simulations can be carried out for an arbitrary arrangement of point defects. The significant crystal field perturbations calculated here demonstrate that a sample that is single phase from a structural point of view can nonetheless have a dramatically varying anisotropy profile at the atomistic level if there is compositional disorder, which may influence localized magnetic objects like domain walls or skyrmions. Published by the American Physical Society 2024

25 ENERGY STORAGE↗

Co-design of Advanced Architectures for Graph Analytics using Machine Learning

A graph is an excellent way of representing relationships among entities. We can use graph analytics to synthesize and analyze such relational data, and extract relevant features that are useful for various tasks such as machine learning. Considering the crucial role of graph analytics in various domains, it is important and timely to investigate the right hardware configurations that can achieve optimal performance for graph workloads on future high-performance computing systems. Design space exploration studies facilitate the selection of appropriate configurations (e.g. memory) to achieve a desired system performance. Recently, the approach of accelerating graph analytics using persistent non-volatile memory has gained a lot of attention. Traditional system simulators such as Gem5 and NVMain can be used to explore the design space of these advanced memory architectures for graph workloads. However, these simulators are slow in execution thus limiting the efficiency of design space exploration studies. To overcome this challenge, we proposed a machine learning based approach to co-design advanced memory architectures for graph workloads. We tested our approach with DRAM, non-volatile memory, and hybrid memory (DRAM+NVM) using a breadth first search benchmark algorithm. Our results showed the applicability of the proposed machine learning based approach to the co-design of the advanced memory architectures. In this paper, we provide recommendations on selecting advanced memory architectures to achieve desired performance for graph workloads. We also discuss the performances of different machine learning models that were considered in this study.

Kurte, Kuldeep↗

Electrifying Airport GSE: Monte Carlo Grid Impacts

Airports globally are shifting from ICE-powered to electric Ground Support Equipment (eGSE) to enhance efficiency, reduce operational costs, and improve operator health. Leveraging predictable routes, flat terrain, and low operational speeds, airports provide ideal conditions for electrification. This study evaluates freight GSE electrification at Dallas-Fort Worth International Airport (DFW), USA, using the Agile@ platform, which integrates three analytical methods: Freight Facility Model (FFM), Activity-Structure-Intensity-Fuel (ASIF), and Monte Carlo simulations. Results from 10,000 simulations indicate modest but critical increases in electricity demand and significant variability in GSE energy consumption. These insights emphasize the importance of data-driven scheduling, targeted maintenance, and strategic infrastructure planning. For high-uncertainty scenarios, airports are advised to deploy buffer energy storage systems (battery banks), implement demand-response charging strategies, schedule flexible workforce shifts, and prioritize proactive maintenance-particularly for equipment with higher operational uncertainty, such as tug tractors with trailers. Agile@ thus offers a robust, scalable, and data-driven framework to optimize long-term GSE planning and enhance reliability across diverse airport environments.

Bose, Ranjan [ORNL] (ORCID:0009000791026327)↗

Estimation of the Water Table Position in Unconfined Aquifers with MODFLOW 6

Abstract The numerical estimation of the position of the water table in unconfined aquifers is important for many practical applications. Its determination through observations or analytical methods is restricted to a few cases. Therefore, it is often estimated through numerical simulations, which may be affected by numerical artifacts and/or poor stability. We use MODFLOW to estimate the position of the water table for a seemingly simple example problem and demonstrate difficulties that can be faced when performing this kind of numerical simulation. We explain the causes for the numerical challenges that originate from the properties of the mathematical equations that must be solved. Based on the results of more than 600 steady‐state simulations, we show how the stability of the numerical solution can be affected by the values of physical parameters that define the problem (e.g., recharge rate, anisotropy ratio, and other parameters that control the numerical algorithm such as settings of the linear and nonlinear solution methods). Finally, we comment on some best practices to apply numerical simulations to estimate the water table position.

Geology↗

A population data-driven workflow for COVID-19 modeling and learning

CityCOVID is a detailed agent-based model that represents the behaviors and social interactions of 2.7 million residents of Chicago as they move between and colocate in 1.2 million distinct places, including households, schools, workplaces, and hospitals, as determined by individual hourly activity schedules and dynamic behaviors such as isolating because of symptom onset. Disease progression dynamics incorporated within each agent track transitions between possible COVID-19 disease states, based on heterogeneous agent attributes, exposure through colocation, and effects of protective behaviors of individuals on viral transmissibility. Throughout the COVID-19 epidemic, CityCOVID model outputs have been provided to city, county, and state stakeholders in response to evolving decision-making priorities, while incorporating emerging information on SARS-CoV-2 epidemiology. Here we demonstrate our efforts in integrating our high-performance epidemiological simulation model with large-scale machine learning to develop a generalizable, flexible, and performant analytical platform for planning and crisis response.

97 MATHEMATICS AND COMPUTING↗

Investigations of nonlinear polarization transfer between obliquely intersecting beams

Electromagnetic waves propagating through plasma can interact nonlinearly through a variety of different mechanisms. The excitation of a plasma beat wave (ions or electrons) can create a refractive index modulation that changes the dispersion of the interacting beams. Alternatively, high-intensity beams can enter the regime where relativistic nonlinearities influence the propagation dynamics. In recent studies [ Opt. Express 29 , 1162 ( 2021 ) OPEXFF 1094-4087 10.1364/OE.413064 ], it was proposed that two beams propagating along the same axis can exchange their polarization state due to nonlinear interaction. Here we present a numerical analysis of two laser beams intersecting in a nonlinear medium at varying angles. Polarization transfer is observed as predicted by analytical theory for a range of angles. For small angles, it is found that filamentation of the interacting beams becomes important. Analytical estimates of the filamentation threshold are presented, and good agreement is found with the simulation data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electrification of Pakistan's Transport System: Modeling Electric Vehicle Penetration and Energy Supply Chain Impacts

Initial analysis of the National Electric Vehicle Policy (NEVP) of Pakistan was undertaken by an integrated energy planning team from November 2019 through August 2020. This technical reference document seeks to present the background, inputs and assumptions, methodology, and results of the policy analysis in detail. The intended audience is the technical modeler, analyst, or reviewer who seeks to understand specifically the modeling effort within this initial analytical phase, with the end-goal of interpreting, re-producing, or modifying the simulations or extending the models to conduct follow-on studies subsequently. These studies may examine in greater detail the energy sector, impacts to emissions, charging station infrastructure needs, and overall benefits-cost trade of the policy, among other areas of interest.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗