Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation framework”

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 469 records · Page 26

System Modeling of Frost/Defrost Cycles in Heat Pumps

Accurately modeling frost/defrost cycles is crucial for guiding heat exchanger design, appropriately sizing the compressor, and quantifying seasonal energy consumption for heating. While previous modelling studies of frosting cycle have predominantly focused on the component level of fin-and-tube coils, there has been limited development of models for system defrosting cycles. This study enhances the existing knowledge by improving a detailed, segment-to-segment heat exchanger model to simulate frost/defrost cycles for both fin-and-tube coils and microchannel heat exchangers. Furthermore, it integrates this component-level model into a steady-state/quasi-steady-state system modeling framework to simulate complete frost/defrost processes, i.e. the ORNL Heat Pump Design Model (HPDM). This comprehensive approach will unveil differences in coil types, and their interactions with the outdoor fan. It will also predict energy consumption patterns and time intervals associated with these processes.

Shen, Bo↗

CHARGE-MAP: An integrated framework to study the multicriteria EV charging infrastructure expansion problem

The widespread adoption of electric vehicles (EVs) in recent years has necessitated the development of effective charging infrastructures. However, charging infrastructure expansion is a multifaceted problem that requires careful consideration of the existing infrastructure, spatiotemporal distribution of charging demands, power-grid capacity, and budget constraints. Here, to approach this complex problem, we present CHARGE-MAP, a data-driven simulation-optimization framework, focused on ensuring meaningful charging experience for individual EV owners. CHARGE-MAP integrates three modules: an agent-based simulation module that estimates spatiotemporal distribution of charging demands by modeling EV adopter mobility and charging behavior; an optimization module that determines optimal new charging station/charger locations and capacities, while minimizing expected detour distances and wait-times with a limited number of new stations; and a power module that determines how to connect the stations to the power grid while maintaining its stability. Using the state of Virginia (consisting of 95 counties and 38 independent cities) as a case study, our results show that CHARGE-MAP can meet the demand of ~198,600 predicted EVs with 1,305 new public charging stations and 2,164 new chargers. It reduces average detour distances for charging by 66% and wait-times at stations by 72% compared to the existing infrastructure. Furthermore, transformer capacity requirement analysis reveals that only 1.8% of residential transformers require upgrades, while over 80% of commercial charging locations can be supported with modest transformer infrastructure (25 to 50 kVA). This indicates that targeted investments can facilitate cost-effective EV integration. Consequently, CHARGE-MAP provides policymakers and urban planners with crucial data-driven insights for effective EV charging infrastructure expansion. Sign up for PNAS alerts.

charging infrastructure↗

A mortar thermomechanical contact computational framework for nuclear fuel performance simulation

Nuclear fuel performance simulations involve the modeling of complex physical phenomena, ranging from fission gas release to fuel swelling and other temperature-induced effects. For light-water reactors (LWRs), swelling of the fuel and the pressure it imposes on the clad when they come into contact causes permanent clad deformation. Accurately characterizing the fuel-cladding interaction, which involves multiple physics, is essential to accurately simulate the fuel/cladding system. Thermomechanical modeling of this problem using a variationally consistent enforcement (e.g., a mortar approach) has been shown to improve the quality of results and facilitate convergence. Here, we present a general multiphysics computational framework for solving nuclear fuel problems using a mortar approach in BISON, a nuclear fuel performance code. In this study analyses show that using the mortar approach, which enables variationally consistent constraint enforcement, improves the quality of results as compared to the more commonly used node-on-face enforcement for representative LWR nuclear fuel simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Possible routes to a strong science case for a joint CLEoPATRA-Roman microlensing survey

The Contemporaneous LEnsing Parallax and Autonomous TRansient Array (CLEoPATRA) is a mis sion concept for a space-based microlensing survey, proposed to enhance science for the Roman Space Telescope (RST) microlensing survey. From CLEoPATRA’s planned vantage point in geosynchronous equatorial orbit (GEO), it will observe microlensing events simultaneously and from a different point of view to RST at L2. This joint L2-GEO observation baseline gives a projected angular separation of order the size of the characteristic angular lensing scale of Free-Floating-Planets (FFPs) in the Galaxy, making it powerful for constraining FFP properties. We find that a strong science case will demonstrate which parts of the FFP parameter space a joint RST-CLEoPATRA survey will give lead ing lens constraints compared with competing ground- and space- based RST observatory companions (e.g., PRIME, Euclid). Both the simulations and framework that connect telescope capabilities to physical constraints already exist for RST FFPs and can be used to demonstrate the power of joint RST-CLEoPATRA survey, with the addition of a simulated CLEoPATRA event blending distribution. Population-level FFP constraints could also be calculated by adapting ongoing black hole microlensing research at Lawrence Livermore National Laboratory. We also find that considering CLEoPATRA as a third space-based observatory (in addition to RST and Euclid/PRIME), which will break microlensing degeneracies leading to unambiguous FFP constraints, is worth investigating as part of a science case. Finally, we highlight possible auxiliary CLEoPATRA science cases including relieving tension between the RST Core Community Surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Logistics simulation of a remediation effort for a hypothetical radiological contamination scenario

To mitigate the effects following a large-scale nuclear or radiological material release in an urban environment and to expedite recovery, the Integrated Wash-Aid Treatment Emergency Reuse System (IWATERS) was developed. IWATERS consists of three operations: washing contaminated surfaces with an ionic wash solution, collecting, and treating the contaminated wash solution on-site to remove contaminants, and reusing the treated solution throughout operations to preserve the clean water resource. This study develops a framework to simulate the logistics of IWATERS deployment, thereby gaining an understanding of the timeline for decontamination operations. For this purpose, the Analysis of Mobility Platform and GoldSim were leveraged for a hypothetical contamination scenario covering 65,200 m 2 of an urban center. The framework reveals that remediation progress is limited by several resources, notably the availability of vermiculite, a reactive clay that is required to treat the contaminated wash solution. Further, this study also presents how the simulation approach can be used to characterize alternatives to reduce the influence of limited resources on operational progress. Overall, this work lays the foundation for evaluating different decontamination methods through detailed logistics simulation, i.e., by refining simulation assumptions and expanding the range of scenarios the simulation can depict.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Application of community data to surface complexation modeling framework development: Iron oxide protolysis

This study presents a comprehensive community data-driven surface complexation modeling framework for simulating potentiometric titration of mineral surfaces. Compiled community data for ferrihydrite, goethite, hematite, and magnetite are fit to produce representative protolysis constants that can reproduce potentiometric titration data collected from multiple literature sources. Using this framework, the impact of surface complexation model type and surface site density (SSD) on the fit quality and protolysis constants can be readily evaluated. For example, the non-electrostatic model yielded a poor data fit compared to diffuse double layer model and constant capacitance models due to the absence of known surface charge effects. Regardless of the choice of iron oxide mineral, pK a1 decreased with increasing SSD while the opposite tendency was observed for pK a2 . This newly developed framework demonstrates a method to reconcile community data-wide potentiometric titration data using Findable, Accessible, Interoperable, Reusable data principles to produce mineral protolysis constants that improve robustness of surface complexation models for applications in metal sorption and reactive transport modeling. The framework is readily expandable (as community data increase) and extensible (as the number of minerals increase). The framework provides a path forward for developing self-consistent, comprehensive, and updateable surface complexation databases for surface complexation and reactive transport modeling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An integrated modeling framework with open architecture for phase field simulation of multi-component alloys

An integrated modeling framework (PanPhaseField) has been developed, which enables a direct and fast coupling between CALPHAD calculations and large-scale phase field simulations for multi-component alloys. Further, it adopts an open architecture allowing for integration of user-defined phase field models in a plug-and-play manner by taking full advantage of the user-friendly graphical interface of Pandat software. The developed modeling platform becomes an enabling tool that can be used to simulate the evolution of spatially varying microstructures of industrial complex alloys for various engineering applications.

36 MATERIALS SCIENCE↗

Null Space Monte Carlo Evaluation of the Plateau to River Model

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model provides a computational framework to simulate the F&T of contaminants in groundwater associated with the 200-PO-1, 200-UP-1, 200-BP-5, and 200-ZP-1 Groundwater Operable Units (OUs) in the Hanford Site Central Plateau. In addition, the model includes adjacent areas and facilities (e.g., the State Approved Land Disposal Site). Intended and anticipated uses of the model include calculating water levels, hydraulic gradients, and groundwater flows throughout the model domain (encompassing the 200 West and 200 East Areas) for use in subsequent F&T calculations for contaminants of concern and developing scale-appropriate, telescopic-mesh refinement models for detailed evaluation of areas within the model domain where required. The overall objective of the modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within the OU boundaries. The objective for the model development phase is to create a common modeling platform that can be used for investigations of the Central Plateau groundwater OUs and areas downgradient toward the Columbia River. The P2R Model calibration to historical data observed at the Hanford Site is documented in CP-57037, Model Package Report for the Plateau to River Model Version 8.3. The purpose of this environmental calculation is to describe a null space Monte Carlo (NSMC) evaluation was conducted with the historic calibration of the P2R model. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. Use of the NSMC approach results in 100 groundwater flow models that are variants of the calibrated P2R Model. These variant models can be used to evaluate uncertainty in model predictions made by the calibrated P2R Model for other analyses. A secondary purpose of the environmental calculation is to establish these variant models for use with other applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An uncertainty visualization framework for large-scale cardiovascular flow simulations: A case study on aortic stenosis

We present a generalizable uncertainty quantification (UQ) and visualization framework for lattice Boltzmann method simulations of high Reynolds number vascular flows, demonstrated on a patient-specific stenosed aorta. The framework combines EasyVVUQ for parameter sampling with large-eddy simulation turbulence modeling in HemeLB, and executes ensembles on the Frontier exascale supercomputer. Spatially resolved metrics, including entropy and isosurface-crossing probability, are used to map uncertainty in pressure and wall shear stress fields directly onto vascular geometries. Two sources of model variability are examined: inlet peak velocity and the Smagorinsky constant. Inlet velocity variation produces high uncertainty downstream of the stenosis where turbulence develops, while upstream regions remain stable. Smagorinsky constant variation has little effect on the large-scale pressure field but increases WSS uncertainty in localized high-shear regions. In both cases, the stenotic throat manifests low entropy, indicative of robust identification of elevated WSS. By linking quantitative UQ measures to three-dimensional anatomy, the framework improves interpretability over conventional 1D UQ plots and supports clinically relevant decision-making, with broad applicability to vascular flow problems requiring both accuracy and spatial insight.

Hemodynamics↗

Molecular Simulation of Functionalized Covalent Organic Framework Membranes for Inorganic Salt Separation

Covalent organic frameworks (COFs) enable molecular-level design of nanochannels for selective separation in pressure-driven membrane processes. Through variation of building blocks and, subsequently, the pore structure and chemistry, membrane performance can be tailored. This study employs nonequilibrium molecular dynamics simulations to theoretically demonstrate the tunable selectivity of COF membranes through a bottom-up functionalization approach. Water and salt transport are evaluated for six β-ketoenamine-linked COFs with varying multilayer thicknesses. For the thinnest multilayer (0.64–0.72 nm), all COFs exhibit low sodium sulfate (Na 2 SO 4 ) rejection (55–66%). However, 20 stacked sheets (6.4–7.2 nm) provide 67–98% rejection, with the sulfonated COF providing the highest Na 2 SO 4 rejection. Analysis of time-resolved ion density profiles reveals that solute rejection is primarily governed by interfacial exclusion arising from pore size and functional group chemistry. Although increasing salt rejection compromises water permeance, the permeance of all COF membranes is at least two orders of magnitude greater than that of a commercially available nanofiltration membrane. Overall, this work guides the rational design of COF membranes for aqueous salt separation.

Nanofiltration↗

Error-correcting codes for fermionic quantum simulation

Utilizing the framework of \mathbb{Z}_2 ℤ 2 lattice gauge theories in the context of Pauli stabilizer codes, we present methodologies for simulating fermions via qubit systems on a two-dimensional square lattice. We investigate the symplectic automorphisms of the Pauli module over the Laurent polynomial ring. This enables us to systematically increase the code distances of stabilizer codes while fixing the rate between encoded logical fermions and physical qubits. We identify a family of stabilizer codes suitable for fermion simulation, achieving code distances of d=2,3,4,5,6,7, allowing correction of any \lfloor \frac{d-1}{2} \rfloor ⌊ d − 1 2 ⌋ -qubit error. In contrast to the traditional code concatenation approach, our method can increase the code distances without decreasing the (fermionic) code rate. In particular, we explicitly show all stabilizers and logical operators for codes with code distances of d=3,4,5. We provide syndromes for all Pauli errors and invent a syndrome-matching algorithm to compute code distances numerically.

Physics↗

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation↗

Simulating long-range coherence of atoms and photons in quantum computers

Lasers and Bose-Einstein condensates (BECs) exhibit macroscopic quantum coherence in seemingly unrelated ways. Lasers possess a well-defined global phase and are characterized by large fluctuations in the number of photons. In BECs of atoms, instead, the number of particles is conserved and the global phase is undefined. Here, we present a unified framework to simulate lasers and BECs states in gate-based quantum computers, by mapping bosonic particles to qubit excitations. Our approach relies on a scalable circuit that measures the total number of particles without destroying long-range coherence. We introduce complementary probes to measure the global and relative phase coherence of a quantum state, and demonstrate their functionality on a Rigetti quantum computer. Our work shows that particle-number conservation enhances long-range phase coherence, highlighting a mechanism used by superfluids and superconductors to gain phase stiffness.

Torre, Emanuele Dalla↗

Rethinking materials simulations: Blending direct numerical simulations with neural operators

Abstract Materials simulations based on direct numerical solvers are accurate but computationally expensive for predicting materials evolution across length- and time-scales, due to the complexity of the underlying evolution equations, the nature of multiscale spatiotemporal interactions, and the need to reach long-time integration. We develop a method that blends direct numerical solvers with neural operators to accelerate such simulations. This methodology is based on the integration of a community numerical solver with a U-Net neural operator, enhanced by a temporal-conditioning mechanism to enable accurate extrapolation and efficient time-to-solution predictions of the dynamics. We demonstrate the effectiveness of this hybrid framework on simulations of microstructure evolution via the phase-field method. Such simulations exhibit high spatial gradients and the co-evolution of different material phases with simultaneous slow and fast materials dynamics. We establish accurate extrapolation of the coupled solver with large speed-up compared to DNS depending on the hybrid strategy utilized. This methodology is generalizable to a broad range of materials simulations, from solid mechanics to fluid dynamics, geophysics, climate, and more.

36 MATERIALS SCIENCE↗

Validation of seismic velocity models in southern California with full-waveform simulations

SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

58 GEOSCIENCES↗

Simulating Crystallization in a Colloidal System Using State Predictive Information Bottleneck Based Enhanced Sampling

Here, we investigate crystal nucleation in supersaturated colloid suspensions using enhanced molecular dynamics simulations augmented with machine learning techniques. The simulations reveal that crystallization in the model colloidal system studied here, with particles interacting through a repulsive screened Coulomb Yukawa potential, proceeds from vapor to dense liquid droplet to crystalline phases across multiple high barriers. Employing a one-dimensional reaction coordinate derived from the State Predictive Information Bottleneck framework, our simulations capture back-and-forth phase transitions across multiple barriers effectively in biased metadynamics simulations. We obtain relative free energy differences between different phases and also quantify the roles of different molecular level features in driving the phase changes.

Chemistry↗

Connected Traffic Signal Coordination Optimization Framework through Network-Wide Adaptive Linear Quadratic Regulator–Based Control Strategy

Traffic congestion in metropolitan areas causes several significant challenges, such as longer travel times, decreased productivity, increased fuel consumption and vehicle emissions, and even severe injuries during crashes. Traffic signal control is a management approach to reduce traffic congestion and allocate the appropriate right of way for safety and mobility efficiency, both in temporal and spatial domains. Here, this study proposes a network-wide adaptive signal control coordination optimization framework based on the linear quadratic regulator algorithm. The traffic flow conditions driven by signal control inputs are formulated based on their network-wide state-space representation. After modeling traffic control regulation constraints, an adaptive linear quadratic regulator algorithm is designed to maximize the network-wide total throughput under the current conditions. Optimal signal control split time durations for multiple intersections in the network are derived by solving the algebraic Riccati equation. Furthermore, the recursive least square parameter estimation method is employed to quantify dynamic traffic condition changes. To verify the effectiveness of this proposed signal control framework, both simulation and real-world experimental tests are conducted for multiple intersections in downtown Chattanooga, Tennessee, United States. In preparation for real-world experimental tests, pipelines for real-time data processing implementation and historical traffic flow data analysis are conducted. The test results demonstrate that the proposed control framework achieves a decrease in travel time by up to 19.4%, total time spent (TTS) by up to 11.9%, and relative queue balance (RQB) by up to 15.6%. The research findings indicate that the proposed signal control framework can be generalized to handle large scale signal control optimization network-wide.

97 MATHEMATICS AND COMPUTING↗

Evaluate the System-Level Impact of Connected and Automated Vehicles Coupled with Shared Mobility: An Agent-based Simulation Approach

With the rapid growth of information and communication technologies, Connected and Automated Vehicles (CAVs) are deemed to be disruptive with the potential to significantly improve overall transportation system efficiency, however, may bring Vehicle Miles Traveled (VMT) increase or other issues. Further, shared mobility systems are another disruptive force that is reshaping our travel patterns. To quantify the combined impact of CAV and shared mobility on travel behavior, traffic performance and energy efficiency, we develop a mesoscopic simulation-based framework for mobility and energy efficiency evaluation considering the disruptive transportation technologies. Under this framework, we develop novel models for energy intensity and modal activity, and evaluated a variety of energy scenarios for different combinations of CAV applications, various levels of automation, roadway characteristics, and traffic conditions, while also varying different vehicle types and fuel/powertrain technologies. Based applying this modeling suite to a calibrated BEAM simulation network in Riverside California, it was found that cooperative automated driving in general will improve mobility, but automated vehicles, even when deployed in a shared autonomous fleet, will likely bring an increase of VMT (up to 36%) due to mode shifts and deadheading. Ride-hailing vehicles typically have better energy efficiency and a higher share of electric vehicles, which helps offset the negative impact from VMT increases when estimating the system-level energy consumption. In general, simulation results show a 6% increase in energy consumption for the scenarios with an increasing shift to ride-hailing modes.

42 ENGINEERING↗