Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large system 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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A geo-visual analysis for exploring the socioeconomic benefits of the heating electrification using geothermal energy

In parallel to population growth and climate change, the rapid pace of urbanization worldwide has led to an enormous increase in energy demand and costs in urban areas. The subsequent energy burden has become an increasing concern for many households in the U.S. Previous studies have revealed that geothermal resources can effectively lower the electricity demand and carbon emissions in large cities. In this paper, we focus on the socioeconomic impacts of geothermal energy on urban systems by presenting an interactive visual analytics dashboard. The dashboard allows urban planners to spatially examine geothermal energy's practical benefits on energy affordability, urban livability, and resilience across the U.S. We compiled a list of socioeconomic metrics by integrating the simulation results from multiple geothermal and building models with multi-domain urban datasets (socioeconomic, demographic, and electricity utility). These metrics are created to characterize the benefits of the heating electrification of buildings using Geothermal Heat Pumps (GHP) for lowering the energy burden of middle-and low-income households nationwide. The visual dashboard employs a combination of multivariate, glyph-based, and geospatial visualization to reveal the variability and patterns in our metrics. We present a pilot study to demonstrate the GHPs' potential as a renewable and affordable solution for increasing the economic and energy grid resilience in U.S cities.

Xu, Haowen↗
DOE OSTI2025

RAG for FLAG: AI Assistance for a Physics Code

Artificial intelligence (AI) has quickly become an important tool in scientific research, where significant efforts are underway to develop tools that will expedite the research process. One area of particular impact is scientific software, which can be particularly complex, and therefore time consuming to learn and use effectively. AI assistants are increasingly helping to streamline the process by performing tasks such as interactively answering user questions or suggesting solutions. Los Alamos National Laboratory (LANL) develops several advanced scientific codes, such as FLAG, which can be used to run multiphysics simulations. With this study, our goal was to develop an AI assistant for FLAG that could help make the process of understanding the software and running physics simulations more efficient. To develop an AI assistant for FLAG, we used a method called retrieval-augmented generation (RAG), which is a technique that uses information from relevant data sources to enhance the accuracy of large language models (LLMs). We used the FLAG user manual and other FLAG documentation as the knowledge base for the RAG system. When a user provides a query, RAG retrieves relevant sections from the knowledge base in response, then uses those excerpts to generate grounded and contextually rich answers. We found that our AI assistant was able to provide context aware answers and source references to user queries. To evaluate performance, we developed a set of 40 benchmark questions and compared the accuracy of the responses to those of two standard LLMs without retrieval. Our AI assistant significantly outperformed the standard LLMs at answering FLAG-related questions, with an 82.5% accuracy rate, compared to 47.5% for both of the standard LLMs. This has the potential to make the process of learning and using FLAG much easier, especially for new users. Ultimately, it supports LANL’s broader mission by empowering scientists and engineers to focus more on discovery and analysis rather than on navigating complex software systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Global fluid turbulence simulations in the scrape-off layer of a stellarator island divertor

Isothermal fluid turbulence simulations have been performed in the edge and scrape-off layer (SOL) of an analytic stellarator configuration with an island divertor, thereby providing numerical insight into edge turbulence in regions around islands in a stellarator. The steady-state transport follows the a curvature drive that is inverse to the major radius ( $1/R$ ) toward the outboard side, but large fluctuations are present throughout the island divertor region, with the average wavelength of similar size to the island width. The system exhibits a prominent $m=2$ , $n=5$ mode, where m is the poloidal mode number and n is the toroidal mode number, although other modes are present. The amplitude and radial extent of the density fluctuations are similar throughout the edge and SOL, but can decrease near island O-points. The fluctuations exhibit a predominantly positive skewness on the outboard midplane, indicating blob-like perturbations for the transport into the outer SOL. It is determined that a point on the separatrix is generally more correlated with regions outside of the SOL than a nearby reference point which does not lie on the separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Wind Energy Forecasting with the Weather Research and Forecasting Model

This was a collaborative effort between Lawrence Livermore National Security, LLC as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Siemens Energy, Inc. (Siemens) to develop a wind resource forecasting tool. LLNL was to develop an independent high-resolution mesoscale modeling capability forecasting tool that could be implemented in conjunction with existing wind farm control and monitoring software to provide forecasting of wind resources using local observations of winds and temperature. Research with LLNL’s state-of-the-art large-eddy simulation meteorological prediction model, based on the community WRF model and innovative turbulence parameterizations, would improve that model’s applicability to large wind farms offshore and in complex terrain. The modeling capability would include uncertainty quantification. Finally, the application of the modeling tool and existing global climate change predictions would enable the delineation of the likely effects of climate change on wind resources. Siemens was to provide high time resolution hub-height wind speed and other meteorological data streams, including temperature profiles from wind farms, for LLNL to incorporate into the modeling system, to validate and tune this forecasting model for their locations of interest. These data streams would also be used for longer-term studies of correlations of wind resources to climate oscillations to indicate how long-term climate change trends may affect the available wind resource. Siemens would also provide information and observations of turbine wakes for incorporation into the modeling tool. By implementing state-of-the-art turbulence parameterizations into a simulation model and/or ensembles of simulation models, and by integrating real-time hub height wind speed and other meteorological datastreams from wind farms into that model or ensemble of models, LLNL would develop a forecasting tool that could be implemented by Siemens as an add-on to existing wind farm control and monitoring software to provide owners with useful resource forecasting. The desired outcome was that the accuracy level of the output would be sufficient to substantiate power output commitments. The final deliverable for this work would consist of a document outlining the algorithms and software tools that could be integrated into Siemens Wind Park Supervisor.

17 WIND ENERGY↗
DOE OSTI2022

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↗
DOE OSTI2022

Numerical and experimental investigation of the geometry dependent layer-wise evolution of temperature during laser powder bed fusion of Ti–6Al–4V

Abstract Laser powder bed fusion (L-PBF) is currently the additive manufacturing process with the widest industrial use for metal parts. Yet some hurdles persist on the way to a widespread industrial serial production, with reproducibility of the process and the resulting part properties being a major concern. As the geometry changes, so do the local boundary conditions for heat dissipation. Consequently, the use of global, geometry-independent processing parameters, which are today’s state of the art, may result in varying part properties or even defects. This paper presents a numerical simulation as a method to predict the geometry-dependent temperature evolution during the build. For demonstration, an overhang structure with varying angles towards the build platform was manufactured using Ti–6Al–4V. A calibrated infrared camera was integrated into a commercial L-PBF system to measure the temperature evolution over time for a total build height of 10 mm, and the results are used for validation of the simulation. It is shown that the simulation is capable of predicting the temperature between layers. The deviations between simulation and measurement remain in single digit range for smaller overhang structures (90°, 60° and 45°). For large overhang structures (30°), the simulation tends to over-predict the temperatures up to 15 °C. Experiments with varying process parameters showed the feasibility of energy reduction as compensation of the heat accumulation produced by overhang structures.

Li, Gefei (ORCID:0000000174837952)↗
DOE OSTI2022

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗
DOE OSTI2022

Automatic Generation of Event Trees and Fault Trees: A Model-Based Approach

In the past few decades, the increasing complexity of modern engineering systems has been driven by the integration of a large number of components whose operations may involve many disciplines (e.g., thermal hydraulics, plant operations, cybersecurity). Most computational tools used by industry and regulators for system safety and reliability assessments are still based on the traditional fault tree (FT) and event tree (ET) approach, which may not be able to capture complex interactions among system constituents. The use of simulation tools has widely increased in the past few decades to improve the fidelity of the reliability and safety analyses. However, the direct use of simulation tools as part of dynamic probabilistic risk assessment (DPRA) methods is not getting traction since (1) modeling the whole system under consideration with DPRA methods may be computationally expensive and unnecessary, and (2) the manual integration of DPRA models into existing state-of-practice probabilistic risk assessment models (i.e., based on FTs and ETs) can be time consuming and prone to errors. Here, in this paper we propose a procedure to overcome this limitation by presenting several algorithms designed to automatically construct subsystem ETs and FTs from DPRA methods for integration into an existing ET/FT system model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Neural network potential for Zr–Rh system by machine learning

Zr–Rh metallic glass has enabled its many applications in vehicle parts, sports equipment and so on due to its outstanding performance in mechanical property, but the knowledge of the microstructure determining the superb mechanical property remains yet insufficient. Here, we develop a deep neural network potential of Zr–Rh system by using machine learning, which breaks the dilemma between the accuracy and efficiency in molecular dynamics simulations, and greatly improves the simulation scale in both space and time. Here, the results show that the structural features obtained from the neural network method are in good agreement with the cases in ab initio molecular dynamics simulations. Furthermore, we build a large model of 5400 atoms to explore the influences of simulated size and cooling rate on the melt-quenching process of Zr 77 Rh 23 . Our study lays a foundation for exploring the complex structures in amorphous Zr 77 Rh 23 , which is of great significance for the design and practical application.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamics of qudit gates and effects of spectator modes on optimal control pulses

Qudit gates for high-dimensional quantum computing can be synthesized with high precision using numerical quantum optimal control techniques. Large circuits are broken down into modules and the tailored pulses for each module can be used as primitives for a qudit compiler. Application of the pulses of each module in the presence of extra modes may decrease their effectiveness due to crosstalk. In this paper, we address this problem by simulating qudit dynamics for circuit quantum electrodynamics systems. As a test case, we take pulses for single-qudit swap gates optimized in isolation and then apply them in the presence of spectator modes each of which are in Fock states. We provide an experimentally relevant scaling formula that can be used as a bound on the fidelity decay. Our results show that frequency shift from spectator mode populations has to be ≲ 0.1 % of the qudit's nonlinearity in order for high-fidelity single-qudit gates to be useful in the presence of occupied spectator modes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Adaptive Path-Following Control for Ground Vehicles Using a Switching Non-Quadratic Lyapunov Function

The application of adaptive control techniques in the development of control systems for intelligent vehicles, especially for ground vehicle path-following controllers, has gained popularity due to their ability to handle large-scale parametric uncertainties. However, the use of a standard quadratic Lyapunov function in existing adaptive control-based path-following controllers can lead to poor transient performance, such as slow convergence and/or large overshoot. To address this limitation, this study proposes the use of a switching non-quadratic Lyapunov function to design a model reference adaptive path-following controller that aims to provide superior transient performance. The stability and signal convergence of the closed-loop system are demonstrated through a Lyapunov-like analysis. Through dSPACE ASM simulation, the effectiveness of the proposed controller is illustrated, which confirms improved tracking performance over a baseline solution.

Zhou, Xingyu↗
DOE OSTI2020

Multi-resolution, Multi-scale Modeling for Scalable Macroalgae Production

Macroalgae harvested along coastlines and in the open ocean has traditionally been used as food for human consumption, animal feed, and fertilizer. The U.S. ARPA-E MARINER Program estimates that the nation has suitable conditions and geography to produce at least 500 million dry metric tons of macroalgae per year, yielding approximately 2.7 quadrillion BTUs of energy in the form of liquid fuel, which is roughly 10% of the nation’s annual transportation energy demand. Adverse environmental effects of nutrient overload and ocean acidification may also be reduced by large-scale macroalgae cultivation in many coastal ocean regions. However, the successful deployment of large-scale marine macroalgae farms for fuel production depends on ambient hydrodynamic conditions and nutrient availability, as well as their interactions with macroalgae farm structures. Pacific Northwest National Laboratory led an ARPA-E MARINER project to develop a set of numerical modeling tools capable of simulating ocean hydrodynamic and biogeochemical processes, macroalgae trajectories for free-floating systems, macroalgae growth and biomass yields, and hydrodynamic load on macroalgae canopies and farm structures using a multi-resolution and multi-scale approach. This set of modeling tools provides a suite of information essential for system design, optimal project siting, risk analysis, and management of macroalgae production systems in the ocean. Better clarity can also help macroalgae system developers reduce deployment costs, operational risk, and potential impacts on the local marine environment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Quantitative Analysis of the Effects of Natural and Human Factors on a Hydrological System in Zhangweinan Canal Basin

Drought is a water deficit state caused by large-scale climate change that cannot be avoided by a water resource management system. Water scarcity refers to the unsustainable utilization of water resources over a long time, which is the result of water management policies. However, water deficits caused by drought and water scarcity often occur simultaneously and are indistinguishable. This study proposes a model-based simulation framework that can quantitatively distinguish natural factors (drought) from human factors (water scarcity) in a hydrological system. The simulation was applied to the Zhangweinan Canal Basin, based on the runoff sequences from 1950 to 2004. The results show that the runoff curve number, soil depth, soil available water, soil evaporation compensation coefficient, base runoff α coefficient and the maximum canopy interception have the highest sensitivity to runoff, and that the calibrated and validated SWAT model can effectively simulate the runoff process in the Zhangweinan Canal Basin and similar areas. Abrupt changes in human activities in 1975 and water scarcity led to the disappearance of the summer peak runoff period in both wet and dry years. Human factors are the main reason for the change in the hydrological system in the study area; the runoff loss caused by human factors is four times that caused by natural factors according to the proposed variable threshold. This study proposes a model-based simulation framework that can help water resource managers to distinguish the effects of drought and water scarcity in water-stressed areas and adjust management accordingly.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Modeling the Interaction Between Energy Efficiency and Demand Response on Regional Grid Scales: Preprint

With increasing penetration of intermittent renewable generation at grid and distributed scales, flexible building loads can provide significant system value and support the evolving needs of the grid. The growing value of load flexibility may complicate the traditional separation between energy efficiency (EE) and demand response (DR). EE measures may compete in some cases with a building’s DR capabilities but complement one another in other cases. EE can also increase or decrease the need for DR at the system level and change the availability of DR to meet system needs. In this study we present a bottom-up approach to modeling interactive effects between EE and DR in buildings within two regions of the US electricity grid. From a library of building simulation models for different buildings and climates, we synthesize system-level demand profiles and the impacts of potential future EE portfolios. Coupling the underlying building models with a database of DR-enabling technologies, we then compute the quantity of DR that can be delivered in each scenario. The results show that EE and DR interactions are largely driven by the timing of EE savings that are measure-specific and the coincidence with system peak demand that is region-specific. We also find that perspective of the impacts matters – for instance that some EE measures reduce the system need for DR but also reduce the DR potential. Our results imply that utility EE and DR programs developed without considering interactive effects may lead to increased grid-management challenges over the long term.

buildings↗
DOE OSTI2022

Coastal Acoustic Buoy for Offshore Wind: Project Synthesis

The population of North Atlantic right whales is critically endangered and their habitat overlaps with offshore windfarm leases. It is therefore imperative that effective mitigation strategies be used to avoid impacts on right whales during the construction of offshore windfarms. The Department of Energy issued FOA Number DE-FOA-0001924 to encourage the development of technology that could monitor large exclusion zones for right whales in order to mitigate potential impact of construction noise on right whales. This report summarizes the past two years of the development and evaluation of the Coastal Acoustic Buoy for Offshore Wind (CABOW) project which aimed to develop technology to monitor large exclusion zones for North Atlantic right whales. Over the course of the project SMRU Consulting have implemented a rigorous design process including comparison of different approaches (e.g., single sensor vs multiple sensors), as well as consideration of placement and timing of acoustic monitoring. We have evaluated critical components of the CABOW system including, reliability, detection range, and bearing accuracy in areas adjacent to offshore windfarm leases in Maryland by conducting 3,536 playbacks of simulated right whale upcalls. The maximum call detection range was 7.5 km when noise was 99 dB re 1µPa rms (50-225 Hz), but this reduced to < 1 km when ambient noise levels were high. Our detection probability in the field was measured as a function of range as well as the source-to-noise level ratio allowing us to build a model to predict the probability of detection under various scenarios (sample size: 3,536 calls x 5 buoys = 17,680). The median bearing error was -0.25° but this is likely an underestimate of error due to experimental design. Using the published recall and precision of the two detectors we implemented in the CABOW system, we estimate that at a recall of 80%, our precision was > 80%, within the range of what we were aiming for in this project. To estimate our exclusion zone false negative and false positive rates, we built a simulation model using the empirical data from our field trial. We modelled three to nine CABOW units placed on the 10 km exclusion zone and estimated our false negative rate to be 1% or less (which was our project goal) and our false positive rate to be between seven and nine percent, slightly above our goal of 5%. However, we also modelled an equivalent PAM system that does not have bearing capabilities and found the false positive rates for that system to be six to eight times higher than the CABOW rate. This higher false positive rate of PAM systems with low spatial information could have significant cost repercussions for offshore wind developers by adding work shutdowns or delays without providing additional protection for right whales. The model we built allows us to explore the placement of PAM systems under various scenarios and will thus help facilitate planning of PAM mitigation systems to meet NOAA Incidental Harassment Authorizations for specific windfarms. We achieved an average system uptime of 98.3%, just below our goal of 99%. The issues that caused these short losses of data have been identified and fixed. Right whale detections and audio clips were typically transferred via radio from the buoys to the base station in two to four seconds. We therefore believe we have developed a highly robust real-time PAM system. Based on the above, we feel we have achieved the stated funding goal of developing a cost-effective and robust real-time PAM system that enables the monitoring of large exclusion zones for right whales during the constructions of offshore windfarms. This should lead to decreased costs and risks for the offshore wind sector while providing robust mitigation for right whales. It is important to state that PAM mitigation will need to be implemented with other mitigation strategies (e.g., visual observers) to provide a complete mitigation strategy to ensure that any effects on right whales from offshore windfarm construction is minimized.

17 WIND ENERGY↗
DOE OSTI2025

Floquet Engineering of Interactions and Entanglement in Periodically Driven Rydberg Chains

Neutral atom arrays driven into Rydberg states constitute a promising approach for realizing programmable quantum systems. Enabled by strong interactions associated with Rydberg blockade, they allow for simulation of complex spin models and quantum dynamics. We introduce a new Floquet engineering technique for systems in the blockade regime that provides control over novel forms of interactions and entanglement dynamics in such systems. Our approach is based on time-dependent control of Rydberg laser detuning and leverages perturbations around periodic many-body trajectories as resources for operator spreading. These time-evolved operators are utilized as a basis for engineering interactions in the effective Hamiltonian describing the stroboscopic evolution. As an example, we show how our method can be used to engineer strong spin exchange, consistent with the blockade, in a one-dimensional chain, enabling the exploration of gapless Luttinger liquid phases. In addition, we demonstrate that combining gapless excitations with Rydberg blockade can lead to dynamic generation of large-scale multipartite entanglement. Experimental feasibility and possible generalizations are discussed.

Floquet systems↗
DOE OSTI2025

Electrochemical Characterization of Dissolved Oxides in Molten FLiNaK

Molten fluoride salts have long been the subject of investigation given their application in a variety of industrial processes. Oxides are a common exogeneous impurity within these systems, and can lead to many operational and regulatory concerns. In recent studies, the electro-oxidation of dissolved metallic oxides has been attributed to a single step oxygen evolution reaction despite voltametric evidence suggesting otherwise. Here, we first use square wave voltammetry to confirm that the oxidation of dissolved oxides in fluorides proceeds via a two-electron transfer reaction and then demonstrate that the presence of a peroxide-mediated redox reaction better explains the electrochemical data. Numerical simulations were then used to support the determination of improved diffusion coefficients to enable electroanalytical measurements of the oxide concentration across at a variety of temperatures. A linear Arrhenius relationship was observed when the experimental data was corrected for ohmic resistance effects. Concentration measurements using the corrected square wave voltammetry data demonstrated high accuracy across a large range of concentrations, while the uncorrected data showed plateauing and high errors. In total, this work serves to rectify previously misconstrued electrochemical data of oxides in molten salts and demonstrate how accurate, in situ concentration measurements can be achieved in real-world systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗