Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “advanced modelling and simulation”

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 811 records · Page 45

Early Exploration of a Flexible Framework for Efficient Quantum Linear Solvers in Power Systems

The rapid integration of renewable energy resources presents formidable challenges in managing power grids. While advanced computing and machine learning techniques offer some solutions for accelerating grid modeling and simulation, there remain complex problems that classical computers cannot effectively address. Quantum computing, a promising technology, has the potential to fundamentally transform how we manage power systems, especially in scenarios with a higher proportion of renewable energy sources. One critical aspect is solving linear systems of equations, crucial for power system applications like power flow analysis, for which the Harrow-Hassidim-Lloyd (HHL) algorithm is a well-known quantum solution. However, HHL quantum circuits often exhibit excessive depth, making them impractical for current Noisy-Intermediate-Scale-Quantum (NISQ) devices. In this paper, we introduce a versatile framework, powered by NWQSim, that bridges the gap between power system applications and quantum linear solvers available in Qiskit. This framework empowers researchers to efficiently explore power system applications using quantum linear solvers. Through innovative gate fusion strategies, reduced circuit depth, and GPU acceleration, our simulator significantly enhances resource efficiency. Power flow case studies have demonstrated up to a eight-fold speedup compared to Qiskit Aer, all while maintaining comparable levels of accuracy.

quantum computing, Harrow-Hassidim-Lloyd, high-per↗

Integrated End-to-end Performance Prediction and Diagnosis for Extreme Scientific Workflows (IPPD) (Final Report)

This report details the accomplishments from the ASCR funded project “Integrated End-to-end Performance Prediction and Diagnosis for Extreme Scientific Workflows” under the award numbers FWP-66406 and DE-SC0012630, with a focus on the UC San Diego (Award No. DE-SC0012630) part of the accomplishments. We refer to the project as IPPD. The main activities of IPPD were centered on the development and integration of provenance information to capture empirically workflow information and identify the sources of bottlenecks as well as variability, and modeling and simulation to provide insights and predictively explore multiple scenarios for development of advanced techniques for optimization of resources and workflow execution. After a Phase I of the project, a Phase II effort focused on three major aspects: a) observe how data is generated, distributed, and used; b) analyze how data is (repeatedly) consumed with a focus both on repeated patterns and anomalies; and c) explore how to optimize data motion. The project leveraged and extended our existing tools with new research and demonstrated our work on the Belle II workflow suite as well as on workflows from NSLS-II.

97 MATHEMATICS AND COMPUTING↗

FY 2021 Isolated Grids and Grid-Connected Turbine Reference Systems; Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL)

For individuals, businesses, and communities focused on building resilient electrical grid infrastructure, wind energy can provide an affordable, accessible, and compatible distributed energy resource option that also enhances the capabilities of local grid operations. The Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL) project is a multi-year distributed wind research effort, driven through a partnership between four Department of Energy National Laboratories and industry to develop and improve the planning, design, and operation of wind-centered microgrids to complement solar, energy storage, and other distributed energy resources for grid-tied and isolated operation (U.S. Department of Energy, 2021). This report documents the application of methods developed through the initial three years of the MIRACL project to two real-world distributed wind reference systems. Specifically, the methods demonstrated in this report include 1) a market valuation framework to comprehensively value the services distributed wind can provide and 2) a resilience framework that enables stakeholders to characterize distribution system resilience and compare grid investment decisions from a resilience perspective. Additional methods mentioned in this report include distributed hybrid system design methods for grid resilience, advanced control co-simulation platforms, and power hardware-in-the-loop (PHIL) models. Preliminary results from these additional methods are presented in this report and will be demonstrated and/or applied to the reference systems in the coming year. The purpose of applying these methods to reference systems is to drive technology transfer of the theories, methodologies, and technologies developed under the MIRACL project and increase the number of referenceable case studies available to stakeholders interested in additional value-added capabilities of wind systems beyond bulk energy supply (i.e. kilowatt-hours).

17 WIND ENERGY↗

Acquire and Test IDAES Framework Components

The cooperation under the Grid Modernization Initiative and the Applied Energy Tri-Laboratory umbrellas requires mutual understanding of capabilities and a common effort to identify gaps existing in the available computational modeling and analysis tools. In this work, we report acquisition and testing of the National Energy Technology Laboratory’s Institute for the Design of Advanced Energy Systems (IDAES) framework, particularly including both electric power systems production cost model simulator Prescient as well as system optimization IDAES itself. Further, we consider the methods in which developments in Prescient and IDAES for grid-energy simulation and optimization may be included in the Framework for Optimization of ResourCes and Economics (FORCE) ecosystem that has been under collaborative development by Idaho National Laboratory, Argonne National Laboratory, and Oak Ridge National Laboratory. We propose two integration methods to introduce Prescient and IDAES into existing FORCE automated workflows. One integration involves using Prescient as a market reduced-order model and IDAES models similar to how existing Modelica models are used. The other integration involves wrapping the Prescient-IDAES dispatch optimization process with the stochastic capacity optimization workflow in FORCE.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The Compelling Need for a Mid-Scale Stellarator Facility

In the pursuit of the goal of commercial fusion as an abundant and safe source of energy, the stellarator is a leading concept with compelling attractiveness and demonstrated performance. A new mid-size stellarator is needed to retire risks and innovate towards a high performance, economically attractive, stellarator Fusion Pilot Plant. In this presentation we, as a community of US researchers from Universities, National Laboratories, and Private Industry, involved in studying the stellarator concept, lay out the programmatic and technical motivation for a new and modern mid-size stellarator research facility. A new mid-scale stellarator is needed to realize the potential predicted by a solid body of theory and simulation along with advances in computational tools for optimization and non-linear turbulence modeling. Notably, it is possible to combine the advantages of the stellarator (steady state, no current drive, no disruptions) with the good confinement regularly achieved in tokamaks. The top priorities for experimental work, and the motivation for a mid-scale stellarator experiment are: turbulence control, non-resonant divertor, MHD stability at large beta, confinement of fast particles, and coil simplification. A new mid-size quasi-symmetric stellarator, built as a user facility, would complement existing research at Wendelstein 7-X and LHD and strongly augment private industry. It would provide a program of innovative research, concept validation, theoretical advancement, and workforce development. Growing support and interest for stellarators by the fusion community and private industry affirms this rationale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Techno-Economic Modelling of Tidal Energy Converter Arrays in the Tacoma Narrows

Hydrokinetic tidal energy converter (TEC) technology is yet to become cost competitive with other renewable energy sources. Understanding the interaction between energy production and the costs incurred harvesting that energy may unlock the economic potential of this technology. Although hydrodynamic simulation of TEC arrays has matured over time, including demonstration of how small and large arrays affect the resource, integration of cost modelling is often limited. The advanced ocean energy array techno-economic modelling tool ‘DTOcean’ enables designers to calculate and improve the levelised cost of energy (LCOE) of an array through parametric simulation of the energy extraction, design of the electrical network, moorings and foundations, and simulation of the installation and lifetime operations and maintenance of the array. This work presents a verification of DTOcean’s ability to simulate the techno-economic performance of TEC arrays by reproducing the hypothetical RM1 reference model, a semi-analytical model of a TEC array based in the Tacoma Narrows of Washington state, U.S.A. It is demonstrated that DTOcean can produce a reasonable estimate to the LCOE predicted by the reference model, giving (in Euro cents per kiloWatt hour) 36.69 ¢/kWh against the reference model’s 34.612 ¢/kWh for 10 TECs, while for 50 TECs, DTOcean calculated 20.34 ¢/kWh compared to 17.34 ¢/kWh for the reference model.

Topper, Mathew B. R. (ORCID:0000000347324347)↗

Evaluating Short-warning Mitigation via Intentional Robust Disruption of a Hypothetical Impact of Asteroid 2023 NT1

We investigate various short-warning mitigation scenarios via fragmentation for a hypothetical impact of asteroid 2023 NT1, a near-Earth object (NEO) that was discovered on 2023 July 15, two days after its closest approach to Earth on July 13. The asteroid passed by Earth within ∼0.25 lunar distances, with a closest approach of ∼1 × 10 5 km and a velocity of 11.27 km s −1 . Its size remains largely uncertain, with an estimated diameter range of 26–58 m and a most probable estimate of 34 m (JPL Sentry, 2023 September 15; weighted by the NEO size frequency distribution). If 2023 NT1 had collided with Earth, it could have caused significant local damage. Assuming a spherical asteroid with a diameter of 34 m, uniform density of 2.6 g cm −3 , and impact velocity of 15.59 km s −1 , a collision would have yielded an estimated impact energy of ∼1.5 Mt, approximately 3 times the energy of the Chelyabinsk airburst in 2013. We analyze the effectiveness of mitigation via intentional robust disruption for objects similar to 2023 NT1. We utilize Pulverize It (PI), a NASA Innovative Advanced Concepts study of planetary defense via fragmentation, to model potential mitigation scenarios through simulations of hypervelocity asteroid disruption and atmospheric ground effects in the case of a terminal defense mode. Simulations suggest that PI is an effective multimodal approach for planetary defense that can operate in extremely short interdiction modes, in addition to long interdiction timescales with extended warning. Our simulations support the proposition that threats like 2023 NT1 can be effectively mitigated with intercepts of 1 day (or less) prior to impact, yielding minimal to no ground damage.

Asteroids↗

Cloud-Feedback Model Intercomparison Project: Tier 2 Simulations (Final Report)

The University of Miami (Subcontractor)’s Research Scientist James Benedict (with oversight by PI Amy Clement) was tasked with producing global climate model simulations as part of the Cloud-Feedback Model Intercomparison Project (CFMIP) Tier 2 protocol, managing the output of these simulations, providing assistance (as requested) to staff at Lawrence Livermore National Lab (LLNL) regarding model output and setup, and meeting with LLNL scientists and team members to coordinate work on the project. The simulations represent an important contribution to the CFMIP data repository and will advance understanding of a wide range of critical cloud, circulation, and precipitation responses to climate change. The Subcontractor completed all proposed simulations, formatted the model output to be compliant with CFMIP protocols, and published the model data to the CFMIP repository. When requested, the Subcontractor provided assistance to LLNL for model output and configuration queries. Benedict and/or Clement also met in person or virtually 1-2 times per year with project scientists to coordinate work and review results.

54 ENVIRONMENTAL SCIENCES↗

Towards Learning-Based Architectures for Sensor Impact Evaluation in Building Controls

Advanced control algorithms for building systems are known to have significant potential in reducing energy consumption while optimizing thermal comfort. The success of such algorithms is critically contingent on several different types of sensor systems, which are in turn, used for continuous monitoring, identification and estimation of several important building states, such as temperatures, humidity, air quality, power consumption and occupancy status. Nonidealities in any of these sensors can lead to significant performance degradation of the control functionalities, and may lead to unwanted sub-optimal building operation. In this paper, we provide a simulation example with a high-fidelity building model, for a particular use-case of advanced optimization-based control in buildings, i.e., occupancy-based controls. We show how imperfections in occupancy sensing can offset performance. Subsequently, we discuss a novel learning-based architecture to efficiently evaluate the impact of sensor nonidealities for building systems, in context of advanced control algorithms.

Bhattacharya, Saptarshi↗

Outdoor Test Bed Performance of a Power Line Sensor Using a Real-Time Event Simulator

This report summarizes the design and application of an outdoor power line sensor testbed (OPLST) with a real-time simulator and power meter to compare a potential transformer (PT) and current transformer (CT) versus an advanced line-post sensor. The OPLST was created to validate advanced medium-voltage (20/34.5 kV) outdoor power line sensors used in electrical distribution systems. Electrical utilities have installed metering and relay-protection transformers like PTs and CTs for several decades. The PTs/CTs are iron core measurement transformers based on the electromagnetic induction principle and provide reliable data in normal grid operation. However, new outdoor power line sensors (OPLS) using other technologies like voltage divider, Rogowski coil and optical principles have become available and may have favorable performance and costs compared to PTs/CTs. Therefore, the importance of testing these technologies with the PT/CT, to compare the measured phase voltage/current at different power grid scenarios is crucial to understand the performance of these new OPLS. For this study a G&W Electric Model CVS-36-O power line sensor was chosen as the OPLS to test with voltage/current signals. An OPAL-RT Technologies Model OP4510 real-time simulator and Schweitzer Engineering Laboratories Model SEL-735 power meter were installed with the 20/34.5 kV OPLST to compare the measured transient events collected from an advanced OPLS and the PT/CT. This system is installed at the Distributed Energy Communications and Control (DECC) lab, Oak Ridge National Laboratory (ORNL). The simulator generated different power grid scenarios (electrical faults, capacitor bank operation, service restoration, etc.), and its analog-output signals were connected to the voltage/current amplifiers that feed the 20/34.5 kV aerial cable loop through the PT/CT devices. Additional PT/CT devices were also wired with the medium voltage aerial cable loop to measure the phase current/voltage signals and server as references. After each test, common format for transient data exchange (COMTRADE) files were collected from the SEL-735 power meter and used to compare the performance of the OPLS with the PT/CT. The behavior of analog signals, harmonic components, total harmonic distortion and crest factors were assessed and found favorable for the G&W sensor as compared to the reference PT/CT.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of point cloud and image-based models for calorimeter fast simulation

Score based generative models are a new class of generative models that have been shown to accurately generate high dimensional calorimeter datasets. Recent advances in generative models have used images with 3D voxels to represent and model complex calorimeter showers. Point clouds, however, are likely a more natural representation of calorimeter showers, particularly in calorimeters with high granularity. Furthermore, point clouds preserve all of the information of the original simulation, more naturally deal with sparse datasets, and can be implemented with more compact models and data files. In this work, two state-of-the-art score based models are trained on the same set of calorimeter simulation and directly compared.

47 OTHER INSTRUMENTATION↗

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

Li, Longlei [Cornell Univ., Ithaca, NY (United Sta↗

ParFlow Sand Tank: A tool for groundwater exploration

The ParFlow Sand Tank model is an open source application designed to allow users to interactively simulate and visualize groundwater movement through the subsurface. The app is designed for both research and education; teaching hydrogeology concepts and making it easy explore and run sophisticated groundwater simulations. Our goal is to support increased accessibility and usability of research grade hydrology tools for research and teaching. The Sand Tank application simulates groundwater and surface water fluxes as well as contaminant transport in real time using the integrated physical hydrology model ParFlow (Kollet & Maxwell, 2006; Maxwell & Miller, 2005; Osei-Kuffuor et al., 2014) and the particle tracking code EcoSlim (Maxwell et al., 2019). ParFlow is a numerical hydrology model that simulates spatially distributed groundwater and surface water flow. It is a well established research tool with more than 90 publications documenting its development use to advance our understanding of groundwater dynamics and groundwater surface water interactions from the hillslope to the continental scale e.g. (Condon et al., 2020; Condon & Maxwell, 2019; Maxwell & Condon, 2016). It is designed for efficient parallel computation and has been run on many platforms spanning from laptops to supercomputers. However, one of the challenges of ParFlow is that it requires significant training and hydrologic expertise to develop simulations. The Sand Tank application makes this model accessible to anyone for education and exploration. Our application uses ParFlow for its simulation backend and ParaView for the data loading and processing. The communication infrastructure relies on the ParaViewWeb framework. We use model templates deployed in Docker images to setup the Sand Tank framework. Users can build the application locally or interact with it through our web deployment. When interacting with a template users can interactively change model parameters like subsurface processes or pump/inject water into the subsurface and watch the system respond to their changes in real time as the simulation runs. Additionally, our template setup will allow more advanced users to build custom templates of increasing complexity for both research and educational purposes.

54 ENVIRONMENTAL SCIENCES↗

Sensitivity of Fine‐Resolution Urban Heat Island Simulations to Soil Moisture Parameterization

ABSTRACT Urban areas experience the impact of natural disasters, such as heatwaves and flash floods, disparately in different neighbourhoods across a city. The demand for precise urban hydrometeorological and hydroclimatological modelling to examine this disparity, and the interacting challenges posed by climate change and urbanisation, has thus surged. The Weather Research and Forecasting (WRF) model has served such operational and research purposes for decades. Recent advancements in WRF, including enhanced numerical schemes and sophisticated urban atmospheric‐hydrological parameterizations, have empowered the simulation of urban geophysical processes at high resolution (~1 km), but even this resolution misses significant urban microclimate variability. This study applies the large‐eddy simulations (LES) mode within WRF, coupled with single‐layer urban canopy models (SLUCM), to enable even finer‐scale modelling (150 m) of the Urban Heat Island (UHI) effect in the Baltimore metropolitan area. We run nine scenarios to evaluate various methods of initializing soil moisture and various spinup lead times, and to assess the impact of WRF's Mosaic approach in depicting subgrid‐scale processes. We evaluate the scenarios by comparing the WRF simulated land surface temperature (LST) against Landsat LST and the WRF simulated hourly 2‐m air temperatures (AT) with observations from eight weather stations across the domain. Results underscore the paramount influence of the lead spinup time on the spatiotemporal distribution of simulated soil moisture, consequently shaping WRF's efficacy in predicting the UHI. Furthermore, interpolating soil moisture‐related parameters from the parent for child domain initialization yields a notable reduction in mean and root‐mean‐squared errors. This improvement was particularly evident in simulations with the longest spinup time, affirming the importance of carefully designing the initialization of soil moisture for improved urban temperature predictions.

Talebpour, Mahdad↗

Analysis of permeability change in dissolving rough fractures using depth-averaged flow and reactive transport models

Many geophysical problems and engineering applications involve fluid flow through evolving faults and fractures. These situations are typical in CO 2 injection scenarios where supercritical CO 2 or CO 2 –saturated fluids are able to dissolve the solid matrix. Dissolution or precipitation within a fracture could lead to self-reinforcing or self-limiting behavior due to feedbacks from changes in permeability (k). There is limited predictive understanding of this phenomenon. Here, we investigated the case of fracture dissolution using a two-dimensional depth-averaged reactive transport model. The model simulated an acidic fluid dissolving a three-dimensional rough fracture. We observed different dissolution patterns across a wide range of geometric properties and different Peclet and Damkohler numbers (960 cases total). Despite this variation, a persistent linear increase in k with time (t) roughly manifested once the dissolution front advanced through the fracture. The time rate of change of k (k/t) is roughly constant for all simulations. The modeled rate is consistent with a theoretical value based on the parallel plates model and with experimentally derived values. Finally, the results confirm the common presence of linear time-dependence of fracture k due to self-reinforcement.

58 GEOSCIENCES↗