Search NASA⌕ Search

SEARCH · Search NASA

Results for “Integrated 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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Multi frame synchrotron radiography of pulsed power driven underwater single wire explosions

We present herein the first use of synchrotron-based phase contrast radiography to study pulsed-power driven high energy density physics experiments. Underwater electrical wire explosions have become of interest to the wider physics community due to their ability to study material properties at extreme conditions and efficiently couple stored electrical energy into intense shock waves in water. The latter can be shaped to provide convergent implosions, resulting in very high pressures (1-10 Mbar) produced on relatively small pulsed power facilities (100s of kA-MA). Multiple experiments have explored single-wire explosions in water, hoping to understand the underlying physics and better optimize this energy transfer process; however, diagnostics can be limited. Optical imaging diagnostics are usually obscured by the shock wave itself; and until now, diode-based X-ray radiography has been of relatively low resolution and rather a broad x-ray energy spectrum. Utilising phase contrast imaging capabilities of the ID19 beamline at the European Synchrotron Radiation Facility, we were able to image both the exploding wire and the shock wave. Probing radiation of 20-50 keV radiographed 200 μm tungsten and copper wires, in ~2-cm diameter water cylinders with resolutions of 8 μm and 32 μm. The wires were exploded by a ~30-kA, 500-ns compact pulser, and 128 radiographs, each with a 100-ps X-ray pulse exposure, spaced at 704 ns apart were taken in each experiment. Abel inversion was used to obtain the density profile of the wires, and the results are compared to two dimensional hydrodynamic and one dimensional magnetohydrodynamic simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A novel, magnetically driven convergent Richtmyer–Meshkov platform

Here, we introduce a novel experimental platform for the study of the Richtmyer–Meshkov instability in a cylindrically converging geometry using a magnetically driven cylindrical piston. Magnetically driven solid liner implosions are used to launch a shock into a liquid deuterium working fluid and, ultimately, into an on-axis rod with a pre-imposed perturbation. The shock front trajectory is tracked through the working fluid and up to the point of impacting the rod through the use of on axis photonic Doppler velocimetry. This configuration allows for precise characterization of the shock state as it impacts the perturbed rod interface. Monochromatic x-ray radiography is used to measure the post-shock interface evolution and rod density profile. The ALEGRA MHD model is used to simulate the dynamics of the experiment in one dimension. We show that late in time the perturbation growth becomes non-linear as evidenced by the observation of high-order harmonics, up to n = 5. Two dimensional simulations performed using a combination of the GORGON MHD code and the xRAGE radiation hydrodynamics code suggest that the late time non-linear growth is modified by convergence effects as the bubbles and spikes experience differences in the pressure of the background flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dynamics of implosion, bounce, and re-implosion of a shell directly driven by a time-shaped laser pulse

Gas-filled plastic spherical shells have been imploded using the 60-beam OMEGA laser. The same laser pulse was used as in a previous experiment, which demonstrated the dynamic formation of a shell by irradiating a homogeneous plastic foam sphere. In the present experiment, the shell initially implodes, bounces, and expands; a new shell then forms and undergoes a second implosion. Analysis of x-ray self-emission-streaked images and sequences of 2D x-ray self-emission framing camera images show shell evolution that is in reasonable agreement with 1D simulations. Shell integrity and symmetry are well maintained until the formation of the new shell. The final implosion, instead, is affected by substantial asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Dynamic Risk Framework for the Physical Security of Nuclear Power Plants

This paper describes ongoing work within the Light Water Reactor Sustainability pathway at Idaho National Laboratory (INL) to optimize the security and cost of nuclear power plants. It introduces the dynamic risk assessment tool developed at INL, Event Modeling Risk Assessment using Linked Diagrams (EMRALD). EMRALD is leveraged to optimize the security posture of a nuclear power plant by integrating force-on-force (FOF) simulations and operator mitigation actions, including dynamic and flexible coping strategies (FLEX). To illustrate the methodology, four attack scenarios are modeled in a commercially available FOF simulation tool using a hypothetical nuclear power plant facility. The simulation results provide valuable insights into possible attack outcomes, as well as the probabilistic risk of a core damage event given these outcomes. Safety mitigation procedures are modeled in EMRALD dependent on the attack outcomes by considering human operator uncertainties. The results demonstrate that the number of armed responders can be optimized, while still maintaining the same protection level as the initial security posture. The proposed modeling and simulation framework of integrating FLEX equipment with FOF models enables the nuclear power plants to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Advanced co-simulation framework for assessing the interplay between occupant behaviors and demand flexibility in commercial buildings

With buildings contributing significantly to electricity usage, enabling demand flexibility becomes a challenge, especially when accounting for occupant comfort. This study introduces an innovative co-simulation framework integrating multiple models: heating, ventilation, and air conditioning (HVAC) system, building zone load, indoor airflow, supervisory control, and occupant comfort and behavior. Uniquely, this framework allows for a comprehensive and dynamic analysis of building systems and occupant interactions in demand response events. Using this framework, we conducted a case study using a typical small office building model. Specifically, we focused on three areas: (1) the impact of indoor airflow modeling on energy use, occupant comfort, and behaviors forecasting, (2) the impact of occupant behaviors on demand flexibility, and (3) occupant comfort and behaviors under demand response events. Key performance indicators such as energy use, flexibility factor, durations of occupant discomfort and occupant behaviors were analyzed. Our findings indicated variations in energy usage and occupant comfort within demand flexibility events, marked by uncertainty boundaries, with variability in demand shedding up to 57.9%. Here, we concluded that this framework is suitable for analyzing typical commercial buildings and their HVAC systems in terms of demand flexibility potential under the impact of occupant behaviors.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2015

Near-Earth asteroid satellite spins under spin-orbit coupling

We develop a fourth-order numerical integrator to simulate the coupled spin and orbital motions of two rigid bodies having arbitrary mass distributions under the influence of their mutual gravitational potential. We simulate the dynamics of components in well-characterized binary and triple near-Earth asteroid systems and use surface of section plots to map the possible spin configurations of the satellites. For asynchronous satellites, the analysis reveals large regions of phase space where the spin state of the satellite is chaotic. For synchronous satellites, we show that libration amplitudes can reach detectable values even for moderately elongated shapes. The presence of chaotic regions in the phase space has important consequences for the evolution of binary asteroids. It may substantially increase spin synchronization timescales, explain the observed fraction of asychronous binaries, delay BYORP-type evolution, and extend the lifetime of binaries. The variations in spin rate due to large librations also affect the analysis and interpretation of light curve and radar observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

THE M31 VELOCITY VECTOR. III. FUTURE MILKY WAY M31-M33 ORBITAL EVOLUTION, MERGING, AND FATE OF THE SUN

We study the future orbital evolution and merging of the Milky Way (MW)-M31-M33 system, using a combination of collisionless N-body simulations and semi-analytic orbit integrations. Monte Carlo simulations are used to explore the consequences of varying all relevant initial phase-space and mass parameters within their observational uncertainties. The observed M31 transverse velocity from Papers I and II implies that the MW and M31 will merge t = 5.86{sup +1.61}{sub -0.72} Gyr from now. The first pericenter occurs at t = 3.87{sup +0.42}{sub -0.32} Gyr, at a pericenter distance of r = 31.0{sup +38.0}{sub -19.8} kpc. In 41% of Monte Carlo orbits, M31 makes a direct hit with the MW, defined here as a first-pericenter distance less than 25 kpc. For the M31-M33 system, the first-pericenter time and distance are t = 0.85{sup +0.18}{sub -0.13} Gyr and r = 80.8{sup +42.2}{sub -31.7} kpc. By the time M31 gets to its first pericenter with the MW, M33 is close to its second pericenter with M31. For the MW-M33 system, the first-pericenter time and distance are t = 3.70{sup +0.74}{sub -0.46} Gyr and r = 176.0{sup +239.0}{sub -136.9} kpc. The most likely outcome is for the MW and M31 to merge first, with M33 settling onto an orbit around them that may decay toward a merger later. However, there is a 9% probability that M33 makes a direct hit with the MW at its first pericenter, before M31 gets to or collides with the MW. Also, there is a 7% probability that M33 gets ejected from the Local Group, temporarily or permanently. The radial mass profile of the MW-M31 merger remnant is significantly more extended than the original profiles of either the MW or M31, and suggests that the merger remnant will resemble an elliptical galaxy. The Sun will most likely ({approx}85% probability) end up at a larger radius from the center of the MW-M31 merger remnant than its current distance from the MW center, possibly further than 50 kpc ({approx}10% probability). There is a {approx}20% probability that the Sun will at some time in the next 10 Gyr find itself moving through M33 (within 10 kpc), but while dynamically still bound to the MW-M31 merger remnant. The arrival and possible collision of M31 (and possibly M33) with the MW is the next major cosmic event affecting the environment of our Sun and solar system that can be predicted with some certainty.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

ON THE CLUSTER PHYSICS OF SUNYAEV-ZEL'DOVICH AND X-RAY SURVEYS. II. DECONSTRUCTING THE THERMAL SZ POWER SPECTRUM

Secondary anisotropies in the cosmic microwave background are a treasure-trove of cosmological information. Interpreting current experiments probing them are limited by theoretical uncertainties rather than by measurement errors. Here we focus on the secondary anisotropies resulting from the thermal Sunyaev-Zel'dovich (tSZ) effect; the amplitude of which depends critically on the average thermal pressure profile of galaxy groups and clusters. To this end, we use a suite of hydrodynamical TreePM-SPH simulations that include radiative cooling, star formation, supernova feedback, and energetic feedback from active galactic nuclei. We examine in detail how the pressure profile depends on cluster radius, mass, and redshift and provide an empirical fitting function. We employ three different approaches for calculating the tSZ power spectrum: an analytical approach that uses our pressure profile fit, a semianalytical method of pasting our pressure fit onto simulated clusters, and a direct numerical integration of our simulated volumes. We demonstrate that the detailed structure of the intracluster medium and cosmic web affect the tSZ power spectrum. In particular, the substructure and asphericity of clusters increase the tSZ power spectrum by 10%-20% at l {approx} 2000-8000, with most of the additional power being contributed by substructures. The contributions to the power spectrum from radii larger than R {sub 500} is {approx}20% at l = 3000, thus clusters interiors (r < R {sub 500}) dominate the power spectrum amplitude at these angular scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Framework for assessment of magnetic equilibrium controller performance on the MAST upgrade spherical tokamak

Here, in this work we present the assessment framework for magnetic equilibrium controllers on MAST Upgrade spherical tokamak (MAST-U) spherical tokamak. Such controllers are essential for the MAST-U since exhaust physics and core-edge integration studies require advanced divertor plasma configurations. The developed framework is based on the TokSys suite of plasma control codes, which was adapted and upgraded for MAST-U. However, extra capabilities were added on top of TokSys to support the development of new control algorithms, deployment of controllers to the plasma control system (PCS) and evaluation of their performance. The controller assessment was realized via closed-loop integrated control simulations with the actual MAST-U PCS and different physics-based plasma models. Since all components of the assessment chain were experimentally validated, these simulations provide qualified controllers applicable for direct use in the experiment. This resulted in the successful experimental demonstration of advanced plasma shape control on MAST-U with minimal on-machine development time. A similar methodology would be beneficial to other tokamaks, both existing and future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Divertor detachment and heat exhaust mitigation control in KSTAR with tungsten divertor

KSTAR has recently undergone an upgrade to use a new tungsten divertor to run experiments in ITER-relevant scenarios. Even with a high melting point of tungsten, it is important to control the heat flux impinging on tungsten divertor targets to minimize sputtering and contamination of the core plasma. Heat flux on the divertor is often controlled by increasing the degree of detachment of scrape-off layer plasma from the target plates. In this work, we have demonstrated successful divertor detachment and heat exhaust dissipation control experiments using two different methods. The first method uses attachment fraction as a control variable which is estimated using ion saturation current measurements from embedded Langmuir probes in the divertor. The second method uses a novel machine-learning-based surrogate model of 2D UEDGE simulation database, DivControlNN. We demonstrated running inference operation of DivControlNN in realtime to estimate heat flux at the divertor and use it as the control variable in a feedback loop with impurity gas flow. We present interesting insights from these experiments including a systematic approach to tuning controllers and discuss future improvements in the control infrastructure and control variables for future burning plasma experiments.

KSTAR tungsten divertor operations↗
DOE OSTI2018

Bringing global gyrokinetic turbulence simulations to the transport timescale using a multiscale approach

The vast separation dividing the characteristic times of energy confinement and turbulence in the core of toroidal plasmas makes first-principles prediction on long timescales extremely challenging. Here in this work, we report the demonstration of a multiple-timescale method that enables coupling global gyrokinetic simulations with a transport solver to calculate the evolution of the self-consistent temperature profile. This method, which exhibits resiliency to the intrinsic fluctuations arising in turbulence simulations, holds potential for integrating nonlocal gyrokinetic turbulence simulations into predictive, whole-device models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Self-consistent core-pedestal ITER scenario modeling

The purpose of these integrated ITER simulations is to identify dependencies that can impact the performance of ITER. The dependence of fusion power production, temperature and density pedestals, and core profiles on varying magnetic-q, edge density fueling strength, neoclassical transport, magnetic field strength, and alpha heating, are examined. It is noted that ion modes usually dominate the interior, while electron modes dominate the barrier. The electron mode without flow shear is substantially stronger. The q increase tends to cause more kinetic ballooning modes at the edge and the q decrease means that the peeling is getting worse, but the ballooning is getting better. The large edge q tends to provide hollow density while small edge q gives normal density profile. This is consistent with the observation that the density of particles is inversely proportional to the density of the plasma current. The rise in the source of the edge particle increases the density of the edge and the reaction rate close to the edge increases temperature fluxes, resulting in weak pedestal barriers. When neoclassical transport is turned off or the strength of B-field is increased, the pronounced edge barriers are identified. In addition, a distinct effect on the magnetohydrodynamic (MHD) modes is observed due to the increased B-field. The interesting aspect is that the barrier collapses about q = 2 as the stability of MHD has decreased to replace the barrier with a continuous slope but with roughly the same fusion power. This would actually be good for a reactor. The slope of the H-mode pedestal is found to be reduced due to the alpha heating. It is expected that the problem of edge-localized modes causing damage to the first wall will be eased. The effects of circular, elongated, and general geometry on plasma profiles are also compared and contrasted. The circular geometry has a hollow density profile and no edge or internal thermal transport barrier (ETB or ITB). However, both ETBs and ITBs can be found in the elongated geometry. In general geometry, there is a weaker ETB, comparable ITBs, but a higher edge particle transport barrier than in elongated geometry. Higher density near the edge, on the other hand, causes more wall erosion, so elongated geometry might be the best option.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development of a data-driven neural network model for electron thermal transport in NSTX

A data-driven electron thermal transport neural network (ETT-NN) model, trained on TRANSP interpretative analysis results of National Spherical Torus Experiment (NSTX), was developed to enable faster and more accurate ETT computation for spherical tokamaks (STs). The model incorporates both convolutional NNs and recurrent NNs, allowing it to simultaneously account for the spatial and temporal non-localities and multi-scale features of turbulent transport, which have been considered only in a limited manner in conventional models. The model was validated through interpretative analysis and predictive simulations using Tokamak Reactor Integrated Automated Suite for Simulation and Computation, demonstrating relatively high accuracy. Additionally, parameter scans were performed on test discharges known to exhibit specific turbulent modes, such as microtearing mode, trapped electron mode, kinetic ballooning mode, and electron temperature gradient mode. The scanning results revealed that the ETT-NN model exhibits the same trends as those observed in conventional gyrokinetic simulations or theories, while also capturing the global nature of turbulent transport, indicating that the data-driven model accurately reflects the underlying physical characteristics. Furthermore, due to the dimensionless nature of the model, we can feasibly expand its applicability by incorporating data from other devices and uncovering the characteristics of ETT in STs in the future.

NSTX↗
DOE OSTI2012

A theoretical study of a nano-opto-mechanical sensor using a photonic crystal-cantilever cavity

In this simulation study, integration of a nanocantilever inside a two-dimensional (2D) photonic crystal (PC) cavity resulted in a unique photonic crystal-cantilever cavity (PC 3 ), where the cantilever served as a tunable mechanical defect of the PC slab. Strong nano-opto-mechanical interactions between the cantilever and the defect-mode field inside the PC 3 gave rise to a high sensitivity of the resonance wavelength to surface stress-induced cantilever deflection. Mechanical and optical responses of the PC 3 to surface stress changes on the cantilever surface were studied by using a finite-element method (FEM) and a finite-difference time-domain (FDTD) method, respectively. Theoretical analysis revealed that the devised PC 3 sensor could resolve a conservative minimum surface stress at the level of ~0.8 mN m -1 , representing state-of-the-art cantilever sensor performance. Also, the PC 3 sensor design used an ultracompact structure with an on-chip optical length of only several microns, while a conventional reflected laser beam detection scheme requires a ~1 m long free-space optical path.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hybrid Energy Management with Real-Time Control of a High-Power EV Charging Site

Decarbonization of transportation systems is driving higher capacity energy storage and faster charging power requirements in electric vehicles (EVs). Given the potential advantages - such as increased efficiency, reduced inverter capacity, and less total cable mass - there is a demand in the industry for more DC distribution for high-power charging (HPC) hubs. However, the cost-effective, adaptive, and robust operation of the DC-coupled HPC hub necessitates a robust site energy management system (SEMS). Validating SEMS operation using a digital twin of an HPC hub in a real-time simulator (RTS) platform is crucial before field deployment. In this study, we propose a hybrid energy management site controller designed to achieve high-level, long-term operational objectives while managing low-level power sharing control between hub assets. We develop a centralized model predictive controller (MPC) to optimize hub operating points and use these points to update the droop parameters of the site energy storage system (ESS). This approach ensures the hub follows an optimal operating point while maintaining the flexibility to respond to load surges. We tested and verified our proposed approach both offline and on a Controller Hardware-in-the-loop (C-HIL) simulation platform integrated with a SEMS framework, demonstrating real-time site operation and validating a cost-effective and robust site controller.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Cyber-Physical Reconfiguration for Disaster Resilience of Power Distribution Systems

Cyber-physical distribution systems (CPDS) have emerged from the integration of information technology into distribution systems. While offering substantial benefits, this integration also introduces vulnerabilities. The interaction between cyber networks and distribution systems renders CPDS susceptible to disasters. To ensure critical load supply and system resilience, rapid post-disaster load restoration is required. The paper proposes a critical load restoration (CLR) framework in CPDS using a network reconfiguration approach that exploits the existing post-disaster resources to restore critical loads within the shortest possible time. Using graph theory, the cyber network and distribution system are integrated into a single digraph, minimizing the CLR complexity in CPDS. A cost metric is also defined to satisfy network-specific objectives and constraints. A heuristic is proposed to guide the load restoration process using the cost metric within the integrated digraph. Simulation results confirm the framework's superiority over existing literature, which either overlooks cyber components or prolongs restoration with additional resource deployment.

cyber-physical system↗
DOE OSTI2022

A framework for evaluating power electronic systems for grid integration

As power electronic systems (PESs) grow increasingly complex, integrating and testing these systems in the development stages is becoming more crucial and more challenging. New modes and resource integration strategies are in development and must be tested before pilot demonstrations and mass production. This paper proposes a framework for testing of PES technologies including early validation with integration of simulated resources. The framework is demonstrated to work in two example configurations as proof of principle.

Dean, Ben↗