Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Engineering evaluation of the EQSIM simulated ground-motion database: The San Francisco Bay Area region

Ground-motion simulations for infrequent earthquake scenarios are gaining increasing interest in the engineering community for their potential to inform performance-based structural design and assessment, particularly in regions where major earthquakes are expected, but recordings from consistent historical events are not available. However, the absence of empirical data makes the evaluation of such scenarios and the interpretation of the simulation results very challenging. In this context, this paper presents the evaluation of the first EQSIM (v.1.1.0) ground-motion simulations generated for the San Francisco Bay Area (SFBA) region. The current database, which is at the first stages of development, includes eight realizations of a Hayward Fault M w 7 event. The evaluation of the simulated motions is first conducted on the average ground-motion amplitudes through comparison against the NGA-W2 GMPEs and a selected population of real records. A methodology for selecting a population of real records consistent with the simulated event is proposed. The objective of the proposed evaluation is twofold: (1) build confidence in the realistic character of the simulated motions for their use in engineering applications and (2) offer a critical physics-based interpretation of the simulation results that can help improve key features of the simulation models. To further investigate the implications of using the simulated ground motions for site-specific structural assessments, the ASCE7-16 approach is employed for the analysis of two hazard-consistent near-field sites in the SFBA. Results are discussed. Furthermore, this study offers a critical review of the EQSIM (v.1.1.0) SFBA ground-motion simulations and suggestions for improving the earthquake rupture models.

58 GEOSCIENCES↗
DOE OSTI2017

Simulation and visualization of energy-related occupant behavior in office buildings

In current building performance simulation programs, occupant presence and interactions with building systems are over-simplified and less indicative of real world scenarios, contributing to the discrepancies between simulated and actual energy use in buildings. Simulation results are normally presented using various types of charts. However, using those charts, it is difficult to visualize and communicate the importance of occupants’ behavior to building energy performance. This study introduced a new approach to simulating and visualizing energy-related occupant behavior in office buildings. First, the Occupancy Simulator was used to simulate the occupant presence and movement and generate occupant schedules for each space as well as for each occupant. Then an occupant behavior functional mockup unit (obFMU) was used to model occupant behavior and analyze their impact on building energy use through co-simulation with EnergyPlus. Finally, an agent-based model built upon AnyLogic was applied to visualize the simulation results of the occupant movement and interactions with building systems, as well as the related energy performance. A case study using a small office building in Miami, FL was presented to demonstrate the process and application of the Occupancy Simulator, the obFMU and EnergyPlus, and the AnyLogic module in simulation and visualization of energy-related occupant behaviors in office buildings. Furthermore, the presented approach provides a new detailed and visual way for policy makers, architects, engineers and building operators to better understand occupant energy behavior and their impact on energy use in buildings, which can improve the design and operation of low energy buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

A Fast Algorithm for Scanning Transmission Electron Microscopy Imaging and 4D-STEM Diffraction Simulations

Scanning transmission electron microscopy (STEM) is an extremely versatile method for studying materials on the atomic scale. Many STEM experiments are supported or validated with electron scattering simulations. However, using the conventional multislice algorithm to perform these simulations can require extremely large calculation times, particularly for experiments with millions of probe positions as each probe position must be simulated independently. Recently, the plane-wave reciprocal-space interpolated scattering matrix (PRISM) algorithm was developed to reduce calculation times for large STEM simulations. Here, we introduce a new method for STEM simulation: partitioning of the STEM probe into “beamlets,” given by a natural neighbor interpolation of the parent beams. This idea is compatible with PRISM simulations and can lead to even larger improvements in simulation time, as well requiring significantly less computer random access memory (RAM). We have performed various simulations to demonstrate the advantages and disadvantages of partitioned PRISM STEM simulations. We find that this new algorithm is particularly useful for 4D-STEM simulations of large fields of view. We also provide a reference implementation of the multislice, PRISM, and partitioned PRISM algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Gyrokinetic GENE simulations of DIII-D near-edge L-mode plasmas

We present gyrokinetic simulations with the GENE code addressing the near-edge region of an L-mode plasma in the DIII-D tokamak. At radial position ρ = 0.80, simulations with the ion temperature gradient (ITG) increased by 40% above the nominal value give electron and ion heat fluxes that are in simultaneous agreement with the experiment. This gradient increase is consistent with the combined statistical and systematic uncertainty σ of the charge exchange recombination spectroscopy measurements at the 1.6σ level. Multiscale simulations are carried out with a realistic mass ratio and geometry for the first time in the near-edge. These multiscale simulations suggest that the highly unstable ion temperature gradient (ITG) modes of the flux-matched ion-scale simulations suppress electron-scale transport, such that ion-scale simulations are sufficient at this location. At radial position ρ = 0.90, nonlinear simulations show a hybrid state of ITG and trapped electron modes, which was not expected from linear simulations. The nonlinear simulations reproduce the total experimental heat flux with the inclusion of E × B shear effects and an increase in the electron temperature gradient by ∼23%. This gradient increase is compatible with the combined statistical and systematic uncertainty of the Thomson scattering data at the 1.3σ level. These results are consistent with previous findings that gyrokinetic simulations are able to reproduce the experimental heat fluxes by varying input parameters close to their experimental uncertainties, pushing the validation frontier closer to the edge region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

TU-H-207A-02: Relative Importance of the Various Factors Influencing the Accuracy of Monte Carlo Simulated CT Dose Index

Purpose: Monte Carlo simulation is a frequently used technique for assessing patient dose in CT. The accuracy of a Monte Carlo program is often validated using the standard CT dose index (CTDI) phantoms by comparing simulated and measured CTDI{sub 100}. To achieve good agreement, many input parameters in the simulation (e.g., energy spectrum and effective beam width) need to be determined. However, not all the parameters have equal importance. Our aim was to assess the relative importance of the various factors that influence the accuracy of simulated CTDI{sub 100}. Methods: A Monte Carlo program previously validated for a clinical CT system was used to simulate CTDI{sub 100}. For the standard CTDI phantoms (32 and 16 cm in diameter), CTDI{sub 100} values from central and four peripheral locations at 70 and 120 kVp were first simulated using a set of reference input parameter values (treated as the truth). To emulate the situation in which the input parameter values used by the researcher may deviate from the truth, additional simulations were performed in which intentional errors were introduced into the input parameters, the effects of which on simulated CTDI{sub 100} were analyzed. Results: At 38.4-mm collimation, errors in effective beam width up to 5.0 mm showed negligible effects on simulated CTDI{sub 100} (<1.0%). Likewise, errors in acrylic density of up to 0.01 g/cm{sup 3} resulted in small CTDI{sub 100} errors (<2.5%). In contrast, errors in spectral HVL produced more significant effects: slight deviations (±0.2 mm Al) produced errors up to 4.4%, whereas more extreme deviations (±1.4 mm Al) produced errors as high as 25.9%. Lastly, ignoring the CT table introduced errors up to 13.9%. Conclusion: Monte Carlo simulated CTDI{sub 100} is insensitive to errors in effective beam width and acrylic density. However, they are sensitive to errors in spectral HVL. To obtain accurate results, the CT table should not be ignored. This work was supported by a Faculty Research and Development Award from Cleveland State University.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Mesoscale to Microscale Simulations over Complex Terrain with the Immersed Boundary Method in the Weather Research and Forecasting Model

Improvements to the Weather Research and Forecasting (WRF) Model are made to enable multiscale simulations over highly complex terrain with dynamically downscaled boundary conditions from the mesoscale to the microscale. Over steep terrain, the WRF Model develops numerical errors that are due to grid deformation of the terrain-following coordinates. An alternative coordinate system, the immersed boundary method (IBM), has been implemented into WRF, allowing for simulations over highly complex terrain; however, the new coordinate system precluded nesting within mesoscale simulations using WRF’s native terrain-following coordinates. Here, the immersed boundary method and WRF’s grid-nesting framework are modified to seamlessly work together. This improved framework for the first time allows for large-eddy simulation over complex (urban) terrain with IBM to be nested within a typical mesoscale WRF simulation. Simulations of the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma, are performed using a multiscale five-domain nested configuration, spanning horizontal grid resolutions from 6 km to 2 m. These are compared with microscale-only simulations with idealized lateral boundary conditions and with observations of wind speed/direction and SF 6 concentrations from a controlled release from intensive observation period 3. The multiscale simulation, which is configured independent of local observations, shows similar model skill predicting wind speed/direction and improved skill predicting SF 6 concentrations when compared with the idealized simulations, which require use of observations to set mean flow conditions. Use of this improved multiscale framework shows promise for enabling large-eddy simulation over highly complex terrain with dynamically downscaled boundary conditions from mesoscale models.

42 ENGINEERING↗
DOE OSTI2018

Large-eddy Simulations of Spray Variability Effects on Flow Variability in a Direct-injection Spark-ignition Engine Under Non-combusting Operating Conditions

Large-eddy Simulations (LES) have been carried out to investigate spray variability and its effect on cycle-to-cycle flow variability in a direct-injection spark-ignition (DISI) engine under non-reacting conditions. Initial simulations were performed of an injector in a constant volume spray chamber to validate the simulation spray set-up. Comparisons showed good agreement in global spray measures such as the penetration. Local mixing data and shot-to-shot variability were also compared using Rayleigh-scattering images and probability contours. The simulations were found to reasonably match the local mixing data and shot-to-shot variability using a random-seed perturbation methodology. After validation, the same spray set-up with only minor changes was used to simulate the same injector in an optically accessible DISI engine. Particle Image Velocimetry (PIV) measurements were used to quantify the flow velocity in a horizontal plane intersecting the spark plug gap. The engine was operated in a skip-fired operating mode and comparisons focused on cycles that included fuel injection, but no spark event and therefore no combustion. 105 total LES engine cycles were simulated using a parallel cycle simulation approach and 3 different perturbation methods in an attempt to isolate the effects of shot-to-shot spray variability and the initial turbulent flow field as well as their interaction effects on overall engine CCVs. The experimental mean and standard deviations were reasonably well matched by the simulations, though quantitative comparisons near the injection event during the intake stroke were difficult due to the high uncertainty in the PIV measurements at these crank angles. The 3 simulation perturbation methods resulted in very similar results, though further analysis found the current parallel cycle approach may be limiting the ability of the simulations to isolate the spray and flow effects.

Cycle-to-cycle Variability↗
DOE OSTI2016

Single-Bubble Flow Boiling Phenomena - Interface Tracking Simulation

Boiling, as one of the most efficient heat transfer mechanisms, is widely used in various engineering systems. Better understanding and modeling of this process remains a major challenge in multiphase flow research. The distribution of vapor in boiling system affects the heat transfer rate and may cause unfavorable conditions, such as heater burn-out. A number of flow regimes have been observed experimentally. In general, there are three regimes under boiling conditions: nucleate boiling, transition boiling, and film boiling. The nucleate boiling regime is categorized into partial boiling regime and fully developed nucleate boiling, according to behavior of bubble dynamics and heat transfer characteristics. In the first region, as the wall temperature increases, the fraction of the wall surface subject to nucleate boiling increases until bubble formation occupies the entire heated surface. Heat transfer in this region is a complex mixture of single-phase forced convection and nucleate boiling. As bubble density increases, bubbles turn to coalesce and form insulating vapor patches, which impedes liquid flow back to the surface. The heat transfer rate under this condition reduces dramatically. This point is called Critical Heat Flux (CHF), which leads to boiling crisis. The behavior of boiling crisis is dependent on local fluid conditions due to bubble growth and departure. Experiments are the ultimate tool for boiling model development and verification, but they still lack some capabilities needed for better understanding of the local mechanism of this phenomena. It's difficult to capture many important details at the pressure/temperature conditions of interest using experimental techniques, such as full three dimensional measurements of velocity and temperature distribution and estimates of heat flux partitioning during the boiling process. The interface tracking simulation (ITS) is one of the promising approaches to describe heat transfer of boiling phenomena and their underlying mechanisms. The advances of high performance computing (HPC) in recent years made it possible to apply DNS to a wide variety of adiabatic flows. Some work has been done to implement evaporation / condensation models as well. Lee and Nydahl performed a numerical simulation of bubble growth and departure assuming that the bubble has the shape of a truncated sphere and a micro-layer exists during bubble growth. A 2-D cylindrical bubbles simulation has been done by Son with a simplified model to predict the heat flux in a thin liquid micro-layer. A direct simulation of 2D film boiling is conducted by Juric and Tryggvason. The film boiling, below critical pressure with a modified Level Set Method, is studied by Son and Dhir. Ose and Kunugi have developed a boiling and condensation model on subcooled boiling phenomena. Fuchs performed numerical simulations of bubble growth and departure considering wall conduction effects. Fully transient heat and fluid flow is modeled with a free surface of a rising bubble and a periodic calculation of the whole cycling of a growing, detaching, and rising bubble. Mukherjee and Kandikar have done a numerical simulation of an evaporating meniscus on a moving heated surface with level set method. The heat transfer associated with the advancing and receding contact angle has been studied. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate single-bubble growth scenarios. For a single bubble simulation in infinite liquid, the analytical solution exists and can be compared with the numerical results for verification purposes. The presented capabilities will be utilized in the large scale DNS of boiling flows in the future. The finite-element based code, PHASTA is used for the simulation. It is a parallel, hierarchic, higher-order accurate, adaptive, stabilized, transient analysis flow solver, which has been shown to be an effective tool for a wide range of single and two-phase flow problems. It has been shown to be highly scalable on top supercomputers ({approx}85% scaling on up to 3 million mesh partitions using a 92 billion element mesh . PHASTA has been shown to reliably predict the details of adiabatic bubbly flows as well as single-phase flows with heat transfer. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate a single-bubble growth scenario. For a single bubble growth in superheated water, the analytical solution exists and can be compared with the numerical results for verification purposes. The second-order accurate spatial and time discretization is used for the flow solver in the presented work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Exploring the impacts of physics and resolution on aqua-planet simulations from a nonhydrostatic global variable-resolution modeling framework

Advances in computing resources are gradually moving regional and global numerical forecasting simulations towards sub-10 km resolution, but global high resolution climate simulations remain a challenge. The non-hydrostatic Model for Prediction Across Scales (MPAS) provides a global framework to achieve very high resolution using regional mesh refinement. Previous studies using the hydrostatic version of MPAS (H-MPAS) with the physics parameterizations of Community Atmosphere Model version 4 (CAM4) found notable resolution dependent behaviors. This study revisits the resolution sensitivity using the non-hydrostatic version of MPAS (NH-MPAS) with both CAM4 and CAM5 physics. A series of aqua-planet simulations at global quasi-uniform resolutions ranging from 240 km to 30 km and global variable resolution simulations with a regional mesh refinement of 30 km resolution over the tropics are analyzed, with a primary focus on the distinct characteristics of NH-MPAS in simulating precipitation, clouds, and large-scale circulation features compared to H-MPAS-CAM4. The resolution sensitivity of total precipitation and column integrated moisture in NH-MPAS is smaller than that in H-MPAS-CAM4. This contributes importantly to the reduced resolution sensitivity of large-scale circulation features such as the inter-tropical convergence zone and Hadley circulation in NH-MPAS compared to H-MPAS. In addition, NH-MPAS shows almost no resolution sensitivity in the simulated westerly jet, in contrast to the obvious poleward shift in H-MPAS with increasing resolution, which is partly explained by differences in the hyperdiffusion coefficients used in the two models that influence wave activity. With the reduced resolution sensitivity, simulations in the refined region of the NH-MPAS global variable resolution configuration exhibit zonally symmetric features that are more comparable to the quasi-uniform high-resolution simulations than those from H-MPAS that displays zonal asymmetry in simulations inside the refined region. Overall, NH-MPAS with CAM5 physics shows less resolution sensitivity compared to CAM4. These results provide a reference for future studies to further explore the use of NH-MPAS for high-resolution climate simulations in idealized and realistic configurations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Biomolecular Solvation Structure Revealed by Molecular Dynamics Simulations

To compare ordered water positions from experiment with those from molecular dynamics (MD) simulations, a number of MD models of water structure in crystalline endoglucanase were calculated. The starting MD model was derived from a joint X-ray and neutron diffraction crystal structure, enabling the use of experimentally assigned protonation states. Simulations were performed in the crystalline state, using a periodic 2 × 2 × 2 supercell with explicit solvent. Water X-ray and neutron scattering density maps were computed from MD trajectories using standard macromolecular crystallography methods. In one set of simulations, harmonic restraints were applied to bias the protein structure toward the crystal structure. For these simulations, the recall of crystallographic waters using strong peaks in the MD water electron density was very good, and there also was substantial visual agreement between the boomerang-like wings of the neutron scattering density and the crystalline water hydrogen positions. An unrestrained simulation also was performed. For this simulation, the recall of crystallographic waters was much lower. For both restrained and unrestrained simulations, the strongest water density peaks were associated with crystallographic waters. Overall, the results demonstrate that it is now possible to recover crystallographic water structure using restrained MD simulations but that it is not yet reasonable to expect unrestrained MD simulations to do the same. Further development and generalization of MD water models for force-field development, macromolecular crystallography, and medicinal chemistry applications is now warranted. In particular, the combination of room-temperature crystallography, neutron diffraction, and crystalline MD simulations promises to substantially advance modeling of biomolecular solvation.

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

A Baseline for Global Weather and Climate Simulations at 1 km Resolution

In an attempt to advance the understanding of the Earth's weather and climate by representing deep convection explicitly, we present a global, four-month simulation (November 2018 to February 2019) with ECMWF's hydrostatic Integrated Forecasting System (IFS) at an average grid spacing of 1.4 km. The impact of explicitly simulating deep convection on the atmospheric circulation and its variability is assessed by comparing the 1.4 km simulation to the equivalent well-tested and calibrated global simulations at 9 km grid spacing with and without parametrized deep convection. The explicit simulation of deep convection at 1.4 km results in a realistic large-scale circulation, better representation of convective storm activity, and stronger convective gravity wave activity when compared to the 9 km simulation with parametrized deep convection. Comparison of the 1.4 km simulation to the 9 km simulation without parametrized deep convection shows that switching off deep convection parametrization at a too coarse resolution (i.e., 9 km) generates too strong convective gravity waves. Based on the limited statistics available, improvements to the Madden-Julian Oscillation or tropical precipitation are not observed at 1.4 km, suggesting that other Earth system model components and/or their interaction are important for an accurate representation of these processes and may well need adjusting at deep convection resolving resolutions. Overall, the good agreement of the 1.4 km simulation with the 9 km simulation with parametrized deep convection is remarkable, despite one of the most fundamental parametrizations being turned off at 1.4 km resolution and despite no adjustments being made to the remaining parametrizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Representation of Leaf-to-Canopy Radiative Transfer Processes Improves Simulation of Far-Red Solar-Induced Chlorophyll Fluorescence in the Community Land Model Version 5

Recent advances in satellite observations of solar-induced chlorophyll fluorescence (SIF) provide a new opportunity to evaluate and constrain the simulation of terrestrial gross primary productivity (GPP). Accurate representation of the processes driving SIF emission and the radiative transfer of SIF to remote sensing sensors is an essential prerequisite for the evaluation and data assimilation. Recently, SIF simulations have been incorporated into several land surface models, but the scaling of SIF from leaf-level to canopy level is usually not well-represented. In this work, we incorporate the simulation of far-red SIF observed at nadir into the Community Land Model version 5 (CLM5). An efficient and accurate method based on escape probability is developed to scale SIF from leaf-level to top-of-canopy while taking clumping and the radiative transfer processes into account. SIF simulated by CLM5 and a canopy-level model agreed well at sites except one in needle leaf forest (R2>0.91, root-mean-square error < 0.19W m -2 sr -1 um -1 ), and captured the day-to-day variation of tower-measured SIF at temperate forest sites (R 2 >0.68). At the global scale, simulated SIF generally captured the spatial and seasonal pat37 terns of satellite-observed SIF (R2 > 0.76 except for tropical forest). Factors including the fluorescence emission model, clumping, bidirectional effect, and canopy properties (leaf optical properties and leaf area index) had considerable impacts on SIF simulation, and the discrepancies between simulated and observed SIF varied with plant functional type. By improving the representation of radiative transfer for SIF simulation, our model allows better comparisons between simulated and observed SIF towards constraining and evaluating GPP simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Large‐Eddy Simulations of a Convection Cloud Chamber: Sensitivity to Bin Microphysics and Advection

Abstract Bin microphysics schemes are useful tools for cloud simulations and are often considered to provide a benchmark for model intercomparison. However, they may experience issues with numerical diffusion, which are not well quantified, and the transport of hydrometeors depends on the choice of advection scheme, which can also change cloud simulation results. Here, an atmospheric large‐eddy simulation model is adapted to simulate a statistically steady‐state cloud in a convection cloud chamber under well‐constrained conditions. Two bin microphysics schemes, a spectral bin method and the method of moments, as well as several advection methods for the transport of the microphysical variables are employed for model intercomparison. Results show that different combinations of microphysics and advection schemes can lead to considerable differences in simulated cloud properties, such as cloud droplet number concentration. We find that simulations using the advection scheme that suffers more from numerical diffusion tends to have a smaller droplet number concentration and liquid water content, while simulation with the microphysics scheme that suffers more from numerical diffusion tends to have a broader size distribution and thus larger mean droplet sizes. Sensitivities of simulations to bin resolution, spatial resolution, and temporal resolution are also tested. We find that refining the microphysical bin resolution leads to a broader cloud droplet size distribution due to the advection of hydrometeors. Our results provide insight for using different advection and microphysics schemes in cloud chamber simulations, which might also help understand the uncertainties of the schemes used in atmospheric cloud simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

CMS Full Simulation for Run 3

We report the status of the CMS full simulation for Run 3. During the long shutdown of the LHC a significant update has been introduced to the CMS code for simulation. The CMS geometry description is reviewed. Several important modifications were needed. CMS detector description software is migrated to the DD4Hep community developed tool. We will report on our experience obtained during the process of this migration. Geant4 10.7 is the CMS choice for Run 3 simulation productions. We will discuss arguments for this choice, the strategy of adaptation of a new Geant4 version, and will report on the physics performance of the CMS simulation. A special Geant4 Physics List configuration FTFP_BERT_EMM will be described, which provides a compromise between simulation accuracy and CPU performance. A significant fraction of time for simulation of CMS events is spent on tracking of charged particles in a magnetic field. In the CMS simulation a dynamic choice of Geant4 parameters for tracking in field is implemented. A new method is introduced into simulation of electromagnetic components of hadronic showers in the electromagnetic calorimeter of CMS. For low-energy electrons and positrons a parametrization of GFlash type is applied. Results of tests of this method will be discussed. In summary, we expect about 25% speedup of the CMS simulation production for Run 3 compared to the Run 2 simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

A new method for analyzing and visualizing plasma simulations using a phase-space tessellation

We apply a novel phase-space interpolation technique referred to as the simplex-in-cell (SIC) method to analyze two- and three-dimensional particle-in-cell (PIC) simulations of electromagnetic plasmas. SIC relies on a discretization of the initial phase-space distribution function into simplices, which allows an approximation to the full, continuously defined distribution function to be constructed at any later time in the simulation. This allows densities, currents, and even full momentum distribution functions to be measured at any point in the simulation domain without averaging over control volumes. The SIC approach applies to any PIC simulation for which a tessellation of the initial particle distribution can be constructed. In this study, we use outputs from standard PIC simulations of the Weibel instability and compare physical quantities such as charge and current densities calculated in postprocessing using SIC and standard particle deposits. Using 2D simulations with 1–65 536 particles-per-cell, we find that SIC eliminates discrete particle noise and in some cases can reach a given noise level using ~1000 times fewer simulation particles than with standard particle deposition schemes. In regions of low density, such as between current filaments, SIC is able to capture small amplitude features even with fewer particles than gridpoints due to the deformable nature of the SIC volume elements. Here, by calculating momentum distributions, we show how SIC can capture low density tails in the spectrum using far fewer particles than are necessary for standard particle deposits. We calculate the charge density on spatial grids of increasing resolution to demonstrate the ability of SIC to reveal fine-scale details that are not accessible with standard particle deposits. Finally, we show how SIC can be extended to 3D and give an example of its use to calculate the charge density from 3D PIC simulations of the Weibel instability. These results motivate the future implementation of SIC directly in the simulation force calculation for a novel low-noise electromagnetic plasma simulation method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

InSituNet: Deep Image Synthesis for Parameter Space Exploration of Ensemble Simulations

We propose InSituNet, a deep learning based surrogate model to support parameter space exploration for ensemble simulations that are visualized in situ. In situ visualization, generating visualizations at simulation time, is becoming prevalent in handling large-scale simulations because of the I/O and storage constraints. However, in situ visualization approaches limit the flexibility of post-hoc exploration because the raw simulation data are no longer available. Although multiple image-based approaches have been proposed to mitigate this limitation, those approaches lack the ability to explore the simulation parameters. Our approach allows flexible exploration of parameter space for large-scale ensemble simulations by taking advantage of the recent advances in deep learning. Specifically, we design InSituNet as a convolutional regression model to learn the mapping from the simulation and visualization parameters to the visualization results. With the trained model, users can generate new images for different simulation parameters under various visualization settings, which enables in-depth analysis of the underlying ensemble simulations. We demonstrate the effectiveness of InSituNet in combustion, cosmology, and ocean simulations through quantitative and qualitative evaluations.

deep learning↗
DOE OSTI2025

Reweighting simulated events using machine-learning techniques in the CMS experiment

Data analyses in particle physics rely on an accurate simulation of particle collisions and a detailed simulation of detector effects to extract physics knowledge from the recorded data. Event generators together with a GEANT -based simulation of the detectors are used to produce large samples of simulated events for analysis by the LHC experiments. These simulations come at a high computational cost, where the detector simulation and reconstruction algorithms have the largest CPU demands. This article describes how machine-learning (ML) techniques are used to reweight simulated samples obtained with a given set of parameters to samples with different parameters or samples obtained from entirely different simulation programs. The ML reweighting method avoids the need for simulating the detector response multiple times by incorporating the relevant information in a single sample through event weights. Results are presented for reweighting to model variations and higher-order calculations in simulated top quark pair production at the LHC. This ML-based reweighting is an important element of the future computing model of the CMS experiment and will facilitate precision measurements at the High-Luminosity LHC.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

SuperNeuro: A Fast and Scalable Simulator for Neuromorphic Computing

In many neuromorphic workflows, simulators play a vital role for important tasks such as training spiking neural networks, running neuroscience simulations, and designing, implementing, and testing neuromorphic algorithms. Currently available simulators cater to either neuroscience workflows (e.g., NEST and Brian2) or deep learning workflows (e.g., BindsNET). Problematically, the neuroscience-based simulators are slow and not very scalable, and the deep learning-based simulators do not support certain functionalities that are typical of neuromorphic workloads (e.g., synaptic delay). In this paper, we address this gap in the literature and present SuperNeuro, which is a fast and scalable simulator for neuromorphic computing capable of both homogeneous and heterogeneous simulations as well as GPU acceleration. We also present preliminary results that compare SuperNeuro to widely used neuromorphic simulators such as NEST, Brian2, and BindsNET in terms of computation times. We demonstrate that SuperNeuro can be approximately 10×--300× faster than some of the other simulators for small sparse networks. On large sparse and large dense networks, SuperNeuro can be approximately 2.2×--3.4× faster than the other simulators, respectively.

Date, Prasanna↗