Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation fidelity”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Verification and validation of detonation-shock-dynamics relations for explosives described by general equation of state and chemical reaction models

Detonation shock dynamics is a powerful method to model the behaviour of High Explosives (HE). However in order to use this method, the underlying relationship between the local radius of curvature and the detonation speed must be known. Previous work has developed methods to calculate this effect using simple, single-step Arrhenius and polytropic gas, models for the chemical reaction and the equation of state, respectively. In recent years, more complex models for both reaction rates and equations of state have been developed which show better agreement with experimental data than these simple models, especially when considering condensed phase explosives.. This work presents the governing equations for solving these problems in a way that is generalised to use arbitrary equations of state as well as reaction models which may have more than a single step and multiple product species. This implementation is verified against exact solutions, demonstrating that the equations were implemented properly. The verified algorithm is then validated against experimental data and high fidelity simulations, showing that it is able to make accurate predictions in a regime where the underlying assumptions of the governing equations are valid. Importantly, this approach has many applications: from creating equivalent detonation shock dynamics models for existing reactive burn calibrations for HE; to developing new functional forms and calibrations of reactive burn models for condensed phase high explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Modeling laser-wakefield accelerators using the time-averaged ponderomotive approximation in a Lorentz boosted frame

Future, high-fidelity simulations of multi-GeV-class laser Wakefield accelerators (LWFAs) will need to model the propagation of high-intensity laser drivers over meter-scale plasmas with high spatial and temporal resolutions, thus requiring high amounts of computational resources. Various techniques have been devised over the years to reduce the computational cost of such simulations, including the time-averaged ponderomotive approximation, and the use of the Lorentz boosted frame technique. In this paper we discuss the combination of these two computational techniques, highlighting the resulting significant reduction in the computational cost of LWFA simulations and the limitations of this approach. The combination of the two techniques can potentially become essential for the modeling of a multi-TeV, LWFA-based collider.

Laser Wakefield Acceleration↗
DOE OSTI2026

DESI DR1 Ly α 1D power spectrum: Validation of estimators

The Data Release 1 (DR1) of the Dark Energy Spectroscopic Instrument (DESI) is the largest sample to date for small-scale Lyα forest cosmology, accessed through its one-dimensional power spectrum (P 1D ). The Lyα forest P 1D is extracted from quasar spectra that are highly inhomogeneous (both in wavelength and between quasars) in noise properties due to intrinsic properties of the quasar, atmospheric and astrophysical contamination, and also sensitive to low-level details of the spectral extraction pipeline. We employ two estimators in DR1 analysis to measure P 1D : the optimal estimator and the fast Fourier transform (FFT) estimator. To ensure robustness of our DR1 measurements, we validate these two power spectrum and covariance matrix estimation methodologies against the challenging aspects of the data. First, using a set of 20 synthetic 1D realizations of DR1, we derive the masking bias corrections needed for the FFT estimator and the continuum fitting bias needed for both estimators. We demonstrate that both estimators, including their covariances, are unbiased with these corrections using the Kolmogorov-Smirnov test. Second, we substantially extend our previous suite of CCD image simulations to include 675,000 quasars, allowing us to accurately quantify the pipeline's performance. This set of simulations reveals biases at the highest k values, corresponding to a resolution error of a few percent. We base the resolution systematics error budget of DR1 P 1D on these values, but do not derive corrections from them since the simulation fidelity is insufficient for precise corrections.

Lyman alpha forest↗
DOE OSTI2025

SICAS, a new integrated divertor-SOL-pedestal-core framework for self-consistent modeling of ion and impurity transport in fusion devices

We present a new framework for core-edge integration studies named SICAS (SOLPS-ITER coupled to ASTRA-STRAHL) which enables high fidelity simulations of the core, edge, and divertor regions encompassing the transport of ions as well as the impurities through the entire plasma domain. SICAS handles the exchanging of the particle and power fluxes as well as transport coefficients to ensure consistency through the codes. An overlapping region is defined from the inside of the separatrix of the plasma to the SOLPS-ITER core boundary allowing for matching profiles and fluxes between the two codes for a self-consistent approach. The results presented here demonstrated the flexibility of SICAS to simulate different configurations, scenarios, divertor geometries, and plasma species with good agreement with DIII-D experimental data. This tool opens new possibilities in integrated modeling of fusion devices integrating all relevant phenomena in the core and the divertor plasmas. These capabilities are required for the interpretation of current experiments as well as the design of new devices.

core-edge integration↗
DOE OSTI2020

Toward flow control: An assessment of the curled wake model in the FLORIS framework

In this work, a new controls-oriented wake model is modified and compared to an analytical Gaussian wake model, high-fidelity simulation data, and experimental wind tunnel campaign. This model, called the curled wake model, captures a wake phenomenon that occurs behind yawed turbines, modeled as a collection of vortices shed from the rotor plane. Through turbine simulations, these vortices are shown to have a significant impact on the prediction of the wake steering's performance. Overall, the results support the concept of secondary steering, or a yawed turbine's ability to deflect the wake of a downstream turbine, and suggest that future turbine wake studies and yaw optimizations should include the curled wake phenomenon.

17 WIND ENERGY↗
DOE OSTI2024

Validation of new and existing methods for time-domain simulations of turbulence and loads

We seek to obtain a second-by-second match between the simulated and measured structural loads of a utility-scale wind turbine. To obtain the one-to-one load simulations, we start with the furthest upstream component of the modeling chain: the turbulent inflow. We consider new and existing methods to generate constrained-turbulence flow fields. The new method is based on large-eddy simulations (LES) and machine learning (ML). The existing methods include Kaimal-based TurbSim and the superstatistical wind field model. The inflow measurements used to constrain these simulations are obtained with a nacelle-mounted scanning lidar. We compare the flow fields for the different inflow simulation approaches and validate their associated load predictions against measurements collected in the Rotor Aero-dynamics, Aeroelastics, and Wake (RAAW) field campaign. We find that the rotor-position control developed for this study is key in enabling the time match between measurements and simulations. When this control approach is used, the load simulation performance tracks with the inflow simulation fidelity, with LES+ML yielding errors ≤ 4% for the damage-equivalent loads of flapwise bending moment, and tower fore-aft bending moments.

17 WIND ENERGY↗
DOE OSTI2024

Plasma-based particle sources

High-brightness particle beams generated by advanced accelerator concepts have the potential to become an essential part of future accelerator technology. In particular, high-gradient accelerators can generate and rapidly accelerate particle beams to relativistic energies. The rapid acceleration and strong confining fields can minimize irreversible detrimental effects to the beam brightness that occur at low beam energies, such as emittance growth or pulse elongation caused by space charge forces. Due to the high accelerating gradients, these novel accelerators are also significantly more compact than conventional technology. Advanced accelerators can be extremely variable and are capable of generating particle beams with vastly different properties using the same driver and setup with only modest changes to the interaction parameters. So far, efforts have mainly been focused on the generation of electron beams, but there are concepts to extend the sources to generate spin-polarized electron beams or positron beams. The beam parameters of these particle sources are largely determined by the injection and subsequent acceleration processes. Although, over the last decade there has been significant progress, the sources are still lacking a sufficiently high 6-dimensional (D) phase-space density that includes small transverse emittance, small energy spread and high charge, and operation at high repetition rate. This is required for future particle colliders with a sufficiently high luminosity or for more near-term applications, such as enabling the operation of free-electron lasers (FELs) in the X-ray regime. Major research and development efforts are required to address these limitations in order to realize these approaches for a front-end injector for a future collider or next-generation light sources. In particular, this includes methods to control and manipulate the phase-space and spin degrees-of-freedom of ultrashort plasma-based electron bunches with high accuracy, and methods that increase efficiency and repetition rate. These efforts also include the development of high-resolution diagnostics, such as full 6D phase-space measurements, beam polarimetry and high-fidelity simulation tools. A further increase in beam luminosity can be achieve through emittance damping. Emittance cooling via the emission of synchrotron radiation using current technology requires kilometer-scale damping rings. For future colliders, the damping rings might be replaced by a substantially more compact plasma-based approach. Here, plasma wigglers with significantly stronger magnetic fields are used instead of permanent-magnet based wigglers to achieve similar damping performance but over a two orders of magnitude reduced length.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Dynamical structure factors of warm dense matter from time-dependent orbital-free and mixed-stochastic-deterministic density functional theory

Abstract We present the first calculations of the inelastic part of the dynamical structure factor (DSF) for warm dense matter (WDM) using time-dependent orbital-free density functional theory (TD-OF-DFT) and mixed-stochastic-deterministic (mixed) Kohn Sham TD-DFT (KS TD-DFT). WDM is an intermediate phase of matter found in planetary cores and laser-driven experiments, where the accurate calculation of the DSF is critical for interpreting x-ray Thomson scattering measurements. Traditional TD-DFT methods, while highly accurate, are computationally expensive, motivating the exploration of TD-OF-DFT and mixed TD-KS-DFT as more efficient alternatives. We applied these methods to experimentally measured WDM systems, including solid-density aluminum and beryllium, compressed beryllium, and carbon–hydrogen mixtures. Our results show that TD-OF-DFT requires a dynamical kinetic energy potential in order to qualitatively capture the plasmon response. Additionally, it struggles with capturing bound electron contributions. In contrast, mixed TD-KS-DFT offers greater accuracy in distinguishing bound and free electron effects, aligning well with experimental data, though at a higher computational cost. This study highlights the trade-offs between computational efficiency and accuracy, demonstrating that TD-OF-DFT remains a valuable tool for rapid scans of parameter space, while mixed TD-KS-DFT should be preferred for high-fidelity simulations. Our findings provide insight into the future development of DFT methods for WDM and suggest potential improvements for TD-OF-DFT.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Modeling flow and transport in fracture networks using graphs

Fractures form the main pathways for flow in the subsurface within low-permeability rock. For this reason, accurately predicting flow and transport in fractured systems is vital for improving the performance of subsurface applications. Fracture sizes in these systems can range from millimeters to kilometers. Although modeling flow and transport using the discrete fracture network (DFN) approach is known to be more accurate due to incorporation of the detailed fracture network structure over continuum-based methods, capturing the flow and transport in such a wide range of scales is still computationally intractable. Furthermore, if one has to quantify uncertainty, hundreds of realizations of these DFN models have to be run. To reduce the computational burden, we solve flow and transport on a graph representation of a DFN. We study the accuracy of the graph approach by comparing breakthrough times and tracer particle statistical data between the graph-based and the high-fidelity DFN approaches, for fracture networks with varying number of fractures and degree of heterogeneity. Due to our recent developments in capabilities to perform DFN high-fidelity simulations on fracture networks with large number of fractures, we are in a unique position to perform such a comparison. We show that the graph approach shows a consistent bias with up to an order of magnitude slower breakthrough when compared to the DFN approach. We show that this is due to graph algorithm's underprediction of the pressure gradients across intersections on a given fracture, leading to slower tracer particle speeds between intersections and longer travel times. We present a bias correction methodology to the graph algorithm that reduces the discrepancy between the DFN and graph predictions. We show that with this bias correction, the graph algorithm predictions significantly improve and the results are very accurate. In conclusion, the good accuracy and the low computational cost, with O(10 4 ) times lower times than the DFN, makes the graph algorithm an ideal technique to incorporate in uncertainty quantification methods.

58 GEOSCIENCES↗
DOE OSTI2020

Leveraging Multiple Connected Traffic Light Signals in an Energy-Efficient Speed Planner

Connecting automated vehicles to traffic lights can lead to significant energy savings by enabling them to pass through intersections in an energy-efficient way without unnecessary stops. A cellular-based communication system connecting multiple traffic lights can help realize the full potential of energy-efficient driving at intersections. Thus, we propose a hierarchical speed planner that can leverage information from multiple connected traffic lights. The proposed speed planner consists of two modules: a green window selector and a reference trajectory generator. The green window selector, based on Dijkstra’s algorithm, finds a series of "green windows" for connected traffic lights that builds an energy-optimal path for vehicles to follow. The reference trajectory generator finds optimal entering times, based on the selected green window at each intersection, and then computes reference trajectories. Deriving and using analytical optimal entering speeds as a function of entering times allows us to guarantee the computational simplicity suitable for real-time implementation. We also demonstrate how to balance energy and traffic flow perspectives in the reference trajectory generator. Lastly, a high-fidelity simulation framework is used to evaluate the proposed speed planner and quantify the extent to which it can save energy in various real-world urban route scenarios.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Blackstart Capability and Survivability of Wind Turbines with Fully Rated Converters

We report the blackstart capability of wind turbines is critical for the recovery of wind-dominant power systems from blackouts. If wind power plants are not able to restore a power system, the incorporation of wind resources into electric grids could be limited by blackstart capability. This paper tackles this problem by setting forth a grid-forming controller for Type 4 wind turbines which employ fully rated power converters. To this end, we engineer two-axis anti-windup proportional-integral regulators that serve to ride through restoration disturbances. We also design an active protection system to prevent wind turbines from stalling during recovery from blackouts. We demonstrate in a high-fidelity simulation environment that Type 4 wind turbines with the developed controllers and protection subsystem can: (i) re-energize a notional wind-dominant grid; (ii) ride through asymmetrical faults; and (iii) survive low wind speed events.

17 WIND ENERGY↗
DOE OSTI2021

Encoding Frequency Constraints in Preventive Unit Commitment Using Deep Learning With Region-of-Interest Active Sampling

With the increasing penetration of renewable energy, frequency response and its security are of significant concerns for reliable power system operations. Frequency-constrained unit commitment (FCUC) is proposed to address this challenge. Despite existing efforts in modeling frequency characteristics in unit commitment (UC), current strategies can only handle oversimplified low-order frequency response models and do not consider wide-range operating conditions. This paper presents a generic data-driven framework for FCUC under high renewable penetration. Here, deep neural networks (DNNs) are trained to predict the frequency response using real data or high-fidelity simulation data. Next, the DNN is reformulated as a set of mixed-integer linear constraints to be incorporated into the ordinary UC formulation. In the data generation phase, all possible power injections are considered, and a region-of-interest active sampling is proposed to include power injection samples with frequency nadirs closer to the UFLC threshold, which enhances the accuracy of frequency constraints in FCUC. The proposed FCUC is investigated on the IEEE 39-bus system. Then, a full-order dynamic model simulation using PSS/E verifies the effectiveness of FCUC in frequency-secure generator commitments.

42 ENGINEERING↗
DOE OSTI2025

Critical Clearing Time for Droop-Controlled Grid-Forming Inverters With Circular Current Limiting and Virtual Impedance Current Limiting

The critical clearing time is a crucial parameter for assessing the large-signal stability of grid-forming (GFM) inverters. Variations in current limiters can significantly impact an inverter's transient behavior during large disturbances. Here, this letter introduces a method based on the power-angle relationship to evaluate the fault recovery of a droop-controlled GFM inverter under a three-phase bolted fault scenario, considering both the circular current limiter and the virtual impedance method. High-fidelity simulations validate the accuracy of the proposed formulas in estimating the critical clearing time of the GFM inverter.

42 ENGINEERING↗
DOE OSTI2022

Blackstart and Fault Ride-Through Capability of DFIG-Based Wind Turbines

Grid-forming wind turbines must be able to re-energize transmission branches, motor loads, unbalanced loads, and power electronics loads following a blackout. Further, they must ride through faults that can arise during restoration. To tackle these challenges, this paper sets forth a grid-forming control strategy for wind turbines using doubly-fed induction generators (DFIGs). In particular, a special set of feedforward control loops are engineered to mitigate the impacts of restoration transients as well as symmetrical and asymmetrical faults on the regulation of DFIG rotor currents and stator voltages. This paper also leverages two-axis anti-windup regulators to bound two-axis voltage and current commands within ratings. As a result, the controller performance is demonstrated via a high-fidelity simulation model of a notional wind-dominant grid.

17 WIND ENERGY↗
DOE OSTI2024

Existence of a Nonzero Worst‐Case ACH for Short‐Term Exposure in Ventilated Indoor Spaces

A well‐ventilated room is essential to reduce the risk of airborne transmission. As such, the scientific community sets minimum limits on ventilation with the idea that increased ventilation reduces pathogen concentration and thus reduces the risk of transmission. In contrast, the upper limit on ventilation is usually determined by human comfort and the need to reduce energy consumption. While average pathogen concentration decreases with increased ventilation, local concentration depends on multiple factors and may not follow the same trend, especially within short exposure times over large separation distances. Here, we show through experiments and high‐fidelity simulations the existence of a worst‐case ventilation where local pathogen concentration increases near the receiving host. This occurs during the type of meetings that were recommended during the pandemic (and in some cases solely authorized) with reduced occupancy adhering to social distancing and short exposure times below 20 minutes. We maintain that for cases of high occupancy and long exposure time, increased ventilation remains necessary.

Construction & Building Technology↗
DOE OSTI2016

TuckerCompressMPI v. 1.0

As parallel computing trends towards the exascale, scientific data produced by high-fidelity simulations are growing increasingly massive. For instance, a simulation on a three-dimensional spatial grid with 512 points per dimension that tracks 64 variables per grid point for 128 time steps yields 8 TB of data. By viewing the data as a dense five-way tensor, we can compute a Tucker decomposition to find inherent low-dimensional multilinear structure, achieving compression ratios of up to 10000 on real-world data sets with negligible loss in accuracy. So that we can operate on such massive data, we present the first-ever distributed-memory parallel implementation for the Tucker decomposition, whose key computations correspond to parallel linear algebra operations, albeit with nonstandard data layouts. Our approach specifies a data distribution for tensors that avoids any tensor data redistribution, either locally or in parallel. This software provides a method for compressing large-scale multiway data.

Austin, Woody↗
DOE OSTI2024

High-Fidelity Numerical Wave Tank Verification & Validation Study: Wave Generation Through Paddle Motion: Preprint

This paper presents a numerical benchmark study of wave propagation due to a paddle motion using different high-fidelity numerical models, which are capable of replicating the nearly actual physical wave tank testing. A full time series of the measured wave generation paddle motion which was used to generate wave propagation in the physical wave tank will be utilized in each of the models contributed by IEA OES Task 10's participants, which includes both computational fluid dynamics (CFD) and smooth hydrodynamic particle (SPH). The high-fidelity simulations of the physical wave testcase will allow for the evaluation of the initial transient effects from wave ramp-up and its evolution in the wave tank over time for two representative regular waves with varying levels of nonlinearity. A couple of interesting metrics like the predicted wave surface elevation at select wave probes, wave period, and phase-shift in time will be assessed to evaluate the relative accuracy of numerical models versus experimental data within specified time intervals. These models will serve as a guide for modelers in the wave energy community and provide a base case to allow further and more detailed numerical modeling of the fixed Kramer Sphere Cases under wave excitation force wave tank testing.

HYDRO ENERGY,TIDAL AND WAVE POWER↗