Some techniques for stable hypersonic cfd simulations
Explore the source record for details and available documents.
Maeng, Jungyeoul Brad [Sandia National Laboratorie↗
SEARCH · Search NASA
Results for “simulation technique”
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.
Probabilistic Data Association for Orbital-Element Estimation Using Multistage Expectation–Maximization
Tracking space objects is important for managing space traffic and predicting collisions, but is difficult in part due to data association and orbit model uncertainty. Expectation–maximization (EM) is a commonly used tracking method that has not been widely considered for tracking space objects. The technique consists of iteratively computing data association probabilities with a set of current element estimates, and updating estimates of the elements by solving a nonlinear weighted least-squares regression problem where the weights are the data association probabilities. In this paper, we demonstrate the use of EM for probabilistic data association and orbital-element estimation by applying the technique to simulated data from two angles-only tracking scenarios. In both scenarios, EM provides correct data associations and accurate maximum likelihood estimates of orbital elements. One scenario considers tracking a single object in clutter and quantifies the improvement of the orbital-element estimates and data associations as the detection probability increases. However, standard application of EM requires knowing the number of objects or may fail when a large number of objects are present. To address these issues, this paper employs a multistage version of EM that is applicable when there are a large and possibly unknown number of objects.
97 MATHEMATICS AND COMPUTING↗
Scenario Creation and Power-Conditioning Strategies for Operating Power Grids with Two-Stage Stochastic Economic Dispatch: Preprint
A significant difficultly associated with the use of stochastic programming to solve optimal power flow problems on a 5-minute timescale is the quality of renewable energy scenarios input by the user. This is especially true when considering power systems with high penetrations of renewable energy, e.g. wind power. This paper introduces the use of stochastic programming to solve the DCOPF problem with scenarios drawn directly from high-fidelity data sets. Hence, the proposed method avoids the problem of lost physics by finding high-fidelity analogs that can describe future states of the system. Furthermore, this method can be simply extended to output multi-period scenarios to the stochastic program. We demonstrate the effectiveness of this technique on simulated dispatch operations of the RTS-GMLC over a week.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Scenario Creation and Power-Conditioning Strategies for Operating Power Grids with Two-Stage Stochastic Economic Dispatch
A significant difficulty associated with the use of stochastic programming to solve optimal power flow problems on a 5-minute timescale is the quality of renewable energy scenarios input by the user. This is especially true when considering power systems with high penetrations of renewable energy, e.g. wind power. This paper introduces the use of stochastic programming to solve the DC optimal power flow problem with scenarios drawn directly from high-fidelity data sets. Hence, the proposed method avoids the problem of lost physics by finding high-fidelity analogs that can describe future states of the system. Furthermore, this method can be simply extended to output multi-period scenarios to the stochastic program. We demonstrate the effectiveness of this technique by simulating dispatch operations on a synthetic test system over the course of a week.
analog forecast↗
Uniform Distribution Technique for Neutrino Beam Scan Simulation
In Fermilab's neutrino facilities such as the Neutrinos at the Main Injector (NuMI) and the upcoming Long Baseline Neutrino Facility (LBNF), a proton beam strikes high-power target, producing positively and negatively charged pions and kaons. There is a need for detailed simulations in order to capture all particle interactions and beam propagation from protons on target to short-lived mesons decaying into muons and neutrinos. The generation of individual beam simulations is a resource-intensive and time-consuming process. In this paper, we describe a method through which many simulation samples with high statistics can be generated to study the effects of beam scan across a target for given beam configurations.
43 PARTICLE ACCELERATORS↗
Data-driven Quasi-static Surrogate Models for Nuclear-powered Integrated Energy Systems
The integration of nuclear power into energy systems presents a promising avenue to address the growing global energy demands while minimizing greenhouse gas emissions. In this paper, we introduce a data-driven quasi-static surrogate model for nuclear-powered Integrated Energy Systems (IES) that comprises various components, including a small modular reactor (SMR), steam manifold, balance of plant (BOP), high-temperature steam electrolysis (HTSE), and district heating (DH) system. Traditional physics-based models for these components often entail significant computational resource and time consumption, necessitating the development of efficient surrogate models. The development of a complete surrogate model for the IES involves the creation of individual surrogate models for each component, leveraging machine learning techniques and simulated data. These isolated surrogate models are subsequently integrated, enabling a holistic view of the IES and reducing the computational burden associated with detailed physics-based simulations. This paper outlines the development process, validation, and the performance evaluation of the surrogate models. The exceptional performance, with low root-mean-squared errors and R-squared scores of at least 99.8% across all individual surrogate models, underscores their accuracy and practical applicability. These results demonstrate the potential of these models to expedite the analysis of nuclear-powered IES, offering insights that can shape future research and development efforts.
08 HYDROGEN↗
The DESI DR1 Peculiar Velocity Survey: Growth Rate Measurements from the Galaxy Power Spectrum
The large-scale structure of the Universe and its evolution encapsulate a wealth of cosmological information. A powerful means of unlocking this knowledge lies in measuring the auto-power spectrum and/or the cross-power spectrum of the galaxy density and momentum fields, followed by the estimation of cosmological parameters based on these spectrum measurements. In this study, we generalize the cross-power spectrum model to accommodate scenarios where the density and momentum fields are derived from distinct galaxy surveys. The growth rate of the large-scale structures of the Universe, commonly represented as $fσ_8$, is extracted by jointly fitting the monopole and quadrupole moments of the auto-density power spectrum, the monopole of the auto-momentum power spectrum, and the dipole of the cross-power spectrum. Our estimators, theoretical models and parameter-fitting framework have been tested using mocks, confirming their robustness and accuracy in retrieving the fiducial growth rate from simulation. These techniques are then applied to analyze the power spectrum of the DESI Bright Galaxy Survey and Peculiar Velocity Survey, and the fit result of the growth rate is $fσ_8=0.440^{+0.080}_{-0.096}$ at effective redshift $z_{\rm eff}=0.07$. By synthesizing the fitting outcomes from correlation functions, maximum likelihood estimation and power spectrum, yields a consensus value of $fσ_8(z_{\rm eff}=0.07) = 0.450 ^{+0.055}_{-0.055}$, and correspondingly we obtain $γ=0.580^{+0.110}_{-0.110}$, $Ω_\mathrm{m}=0.301^{+0.011}_{-0.011}$ and $σ_8=0.834^{+0.032}_{-0.032}$. The measured $fσ_8$ and $γ$ are consistent with the prediction of the $Λ$ Cold Dark Matter Model and General Relativity.
Qin, F. [Marseille, CPPM] (ORCID:0000000179507864)↗
Design and Modeling of the Charge Readout of a SiMOS Quantum Dot with a Single Electron Transistor and CryoCMOS
Single electron spin qubits trapped in SiMOS quantum dots are a promising technology for scaling to thousand- or million-qubit systems due to their compatibility with mature CMOS manufacturing processes. A readout system that combines a single electron transistor with a custom cryogenic CMOS amplification and digitization chain offers key advantages by avoiding the use of bulky RF components or room-temperature interconnects. We present design techniques and simulation results for an optimized qubit-SET cryoCMOS interface, culminating in the design of the QNDR1 ASIC, the first cryogenic readout ASIC designed under the Quandarum project, which targets the development of a many-channel spin qubit based detector for use in high energy physics.
Quinn, Adam [Fermilab] (ORCID:0009000605056371)↗
Geant4 Event Biasing and Fast Simulation
Geant4 offers advanced event biasing techniques to significantly accelerate simulations involving rare events. Various biasing methods, such as leading particle selection, cross-section biasing, radioactive decay enhancement, and bremsstrahlung splitting, enable efficient event sampling, though they require careful handling. Additionally, Geant4 provides a Fast Simulation Interface, allowing the replacement of standard processes in specific region and for selected particles, enabling faster execution or external code integration. Applications of fast simulation include electromagnetic shower modeling in calorimeters, machine learning inference, and offloading tasks to specialized hardware like GPUs, making Geant4 a powerful tool for computationally demanding simulations.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Variance reduction techniques for Monte Carlo neutron noise simulations
The small fluctuations of the neutron flux caused by small perturbations of the macroscopic cross-sections take the name of neutron noise. Advanced Monte Carlo methods have been recently proposed in order to solve the neutron noise equations in the frequency domain, which allows establishing reference solutions to validate faster but approximate deterministic solvers. Due to the presence of particles carrying two statistical weights (for the real and imaginary components of the noise field), both of which may be positive or negative, Monte-Carlo simulations of neutron noise pose distinct challenges in terms of variance reduction. In this work we investigate two variance-reduction techniques, namely branchless collisions and weight cancellation, and probe their effectiveness for a benchmark con- figuration concerning the noise field induced by a pin with oscillating cross sections in a fuel assembly.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
An Enhanced Virtual Crack Closure Technique for Stress Intensity Factor Calculation along Arbitrary Crack Fronts and the Application in Hydraulic Fracturing Simulation
The virtual crack closure technique (VCCT) is widely used for calculating energy release rates along crack fronts and modeling the propagation of cracks in solid materials. Although the VCCT formulation for smooth crack fronts has been sufficiently addressed in the literature, the application of VCCT to a nonsmoothed crack front with sharp corners warrants further investigation. In this study, we describe an enhanced VCCT to calculate energy release rates and stress intensity factors for cracks with arbitrary shapes in 3D domains discretized on structured grids. The formulations of the enhanced VCCT were developed and implemented into a multiphysics simulation environment capable of simulating crack propagation in the framework of linear elastic fracture mechanics. Comparisons with existing analytical/numerical solutions and other VCCT approaches were performed to verify the enhanced VCCT in terms of SIF calculation along nonsmoothed crack fronts. We then applied the enhanced VCCT to a hydraulically driven penny-shaped fracture problem to further demonstrate its capability to simulate nonsmoothed fracture propagation.
02 PETROLEUM↗
Magnetic Designs of New First Target Beamline Magnets for the ORNL SNS Upgrade
The Spallation Neutron Source (SNS) at ORNL currently is being upgraded from 1.0 GeV to 1.3 GeV. Several water-cooled magnets should be upgraded to transport 30% of higher beam energy. New chicane, injection/extraction septum, and Lambertson magnets were designed. Designing the magnets was a challenging task because the new magnets required good combined integrated field quality and needed to occupy the old magnets space but with about 20% greater integrated magnetic field. Additional strong requirements applied to the magnets fringe field do not disturb the circulating beam. The special field profiles had to be provided in foil areas between magnets. The analysis described here was based on OPERA3D simulations. A special technique was used for analyzing the integrated field harmonics. Initially, the particle track was simulated, and integrated field components were calculated along this track for the reference radius, which were used for the harmonics analysis. In addition, 3D field maps were provided for beam optics simulations. The final beamline analysis confirmed good beam transmission and low losses.
3d field maps↗
Scalable Data Reduction Techniques for Extreme-Scale Unstructured PDE Simulations
The present report summarizes activities performed over the 3-year (plus one year of NCE) duration of the award DE-SC0022283, “Scalable Data Reduction Techniques for Extreme-Scale Unstructured PDE Simulations,” funded by the DOE Office of Advanced Scientific Computing Research (ASCR).
Doostan, Alireza↗
Quantum simulation of boson-related Hamiltonians: techniques, effective Hamiltonian construction, and error analysis
Elementary quantum mechanics proposes that a closed physical system consistently evolves in a reversible manner. However, control and readout necessitate the coupling of the quantum system to the external environment, subjecting it to relaxation and decoherence. Consequently, system-environment interactions are indispensable for simulating physically significant theories. A broad spectrum of physical systems in condensed-matter and high-energy physics, vibrational spectroscopy, and circuit and cavity QED necessitates the incorporation of bosonic degrees of freedom, such as phonons, photons, and gluons, into optimized fermion algorithms for near-future quantum simulations. In particular, when a quantum system is surrounded by an external environment, its basic physics can usually be simplified to a spin or fermionic system interacting with bosonic modes. Nevertheless, troublesome factors such as the magnitude of the bosonic degrees of freedom typically complicate the direct quantum simulation of these interacting models, necessitating the consideration of a comprehensive plan. This strategy should specifically include a suitable fermion/boson-to-qubit mapping scheme to encode sufficiently large yet manageable bosonic modes, and a method for truncating and/or downfolding the Hamiltonian to the defined subspace for performing an approximate but highly accurate simulation, guided by rigorous error analysis. In this pedagogical tutorial review, we aim to provide such an exhaustive strategy, focusing on encoding and simulating certain bosonic-related model Hamiltonians, inclusive of their static properties and time evolutions. Specifically, we emphasize two aspects: (1) the discussion of recently developed quantum algorithms for these interacting models and the construction of effective Hamiltonians, and (2) a detailed analysis regarding a tightened error bound for truncating the bosonic modes for a class of fermion-boson interacting Hamiltonians.
bosonic Hamiltonian↗
Simulation of Creep Deformation and Failure in Graded AM Microstructures
This report describes modeling tools and techniques developed to simulate the long-term material performance of 316H stainless steel manufactured using Laser Powder Bed Fusion (LPBF). A physics-based Crystal Plasticity Finite Element model is used to simulate creep in microstructures and to study the roles of grain morphology, porosity, and texture. We describe our modeling methodology, including an orientation-mapping technique to capture the spatially varying crystallographic orientation that results from the build conditions. Our study of microstructural features shows that AM microstructures produced by LPBF tend to creep faster in the build direction, while texture and grain boundaries strengthen the transverse directions. However, when grain-boundary porosity and the consequent cavity growth are included in the model, the transverse directions begin to creep faster. In examining texture, the results indicate that spatially varying orientation arising from the build conditions increases anisotropy in the material, making it critical to account for orientation gradients in the material to accurately model its mechanical behavior. We also describe a material-model calibration campaign in which we calibrated the constitutive model specifically for LPBF 316H stainless steel at 725℃ for both solution-annealed and as-built conditions. Finally, these tools and techniques are used to model creep in microstructures representing different regions of an LPBF material with graded microstructure, owing to intentional variation in processing conditions. The creep simulation results show good agreement with experimental data across all three microstructures, with future work planned to study rupture in the material.
36 MATERIALS SCIENCE↗
Assessment of simulated and observed cavitation-induced erosion damage in Spallation Neutron Source target vessels
Cavitation-induced erosion damage in different Spallation Neutron Source (SNS) target designs are simulated using explicit finite element–based techniques and compared with observations of erosion in targets after operation. The efficacy of the previously developed method, called saturation time, was evaluated using erosion-damaged samples from new target designs. A new metric called maximum bubble size was implemented under the rationale that larger cavitation bubbles will collapse more intensely. The maximum cavitation bubble size over 1 ms of simulated time was calculated based on the Rayleigh–Plesset equation for each element integration point and presented as a contour map at the vessel surface for assessing with erosion observations. SNS targets are now operated with helium gas injection to reduce cavitation damage. A simulation method using a material model for the mixture of mercury and gas bubbles was recently developed and used to account for the effect of small gas bubbles on the structural response of the target vessel. Furthermore, this work compares the new method's results with observed cavitation damage. Maps of the calculated maximum bubble size for targets operated with and without gas injection were compared with photographs of erosion damage observed in SNS targets. The patterns in maximum bubble size maps correlated well with observations of erosion patterns in target vessels after service. Advantages and challenges of the maximum bubble size simulation technique are provided, and differences between results from the previous and the newly proposed metric are discussed.
Jiang, Hao↗
Simulating the Multi-epoch Direct Detection Technique to Isolate the Thermal Emission of the Non-transiting Hot Jupiter HD187123b
We report the 6.5σ detection of water from the hot Jupiter HD187123b with a Keplerian orbital velocity K {sub p} of 53 ± 13 km s{sup −1}. This high-confidence detection is made using a multi-epoch, high-resolution, cross-correlation technique, and corresponds to a planetary mass of 1.4{sub −0.3}{sup +0.5} M {sub J} and an orbital inclination of 21° ± 5°. The technique works by treating the planet/star system as a spectroscopic binary and obtaining high signal-to-noise, high-resolution observations at multiple points across the planet’s orbit to constrain the system’s binary dynamical motion. All together, seven epochs of Keck/NIRSPEC L-band observations were obtained, with five before the instrument upgrade and two after. Using high-resolution SCARLET planetary and PHOENIX stellar spectral models, we were able to drastically increase the confidence of the detection by running simulations that could reproduce, and thus remove, the nonrandom structured noise in the final likelihood space well. The ability to predict multi-epoch results will be extremely useful for furthering the technique. Here, we use these simulations to compare three different approaches to combining the cross correlations of high-resolution spectra and find that the Zucker log(L) approach is least affected by unwanted planet/star correlation for our HD187123 data set. Furthermore, we find that the same total signal-to-noise ratio (S/N) spread across an orbit in many, lower S/N epochs rather than fewer, higher S/N epochs could provide a more efficient detection. This work provides a necessary validation of multi-epoch simulations, which can be used to guide future observations and will be key to studying the atmospheres of farther separated, non-transiting exoplanets.
79 ASTRONOMY AND ASTROPHYSICS↗
Electromagnetic Transient Simulation of Large-Scale Inverter-Based Resources With High-Granularity
The power grid is undergoing a significant transformation with the rapid increase in inverter-based resources (IBRs), including large-scale photovoltaic (PV) plants. Ensuring reliable and resilient grid operation in this new paradigm necessitates high-granularity electromagnetic transient (EMT) modeling that accurately captures the behavior of individual inverters and their interactions within IBR plants. Central to this approach is the detailed representation of both the IBR plant’s collector system and the dynamics of individual inverters. To achieve this, a high-granularity EMT model of a large-scale PV plant has been developed using advanced simulation algorithms, including matrix splitting and the Schur complement. These proposed techniques significantly enhance simulation speed, numerical stability, and accuracy while improving the modularity and efficiency of the collector system’s representation. The effectiveness of the proposed methods is validated through simulations of a representative large-scale PV plant consisting of 125 individual PV inverters, 25 IBR unit transformers, and a 52-bus collector system.
Choi, Jongchan [Oak Ridge National Laboratory (ORN↗