Core-collapse supernovae simulations
Explore the source record for details and available documents.
Massively-parallel simulations, astrophysics, perf↗
SEARCH · Engineering Papers
Results for “core 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.
Simulation scaling studies of reactor core two-phase flow using direct numerical simulation
Tremendous growth in supercomputing power in recent years has resulted in the emergence of high-resolution flow analysis methods as an advanced research tool to evaluate single and two-phase flow behavior. In particular, unstructured mesh-based methods have been applied to analyze flows in complex reactor core geometries, including those of light water reactors (LWR). The finite-element based code, PHASTA, is utilized to perform large-scale simulations of two-phase bubbly flows in LWR geometries. Given the large computational cost of direct numerical simulation (DNS) coupled with interface tracking methods (ITM), typical domains encompass a portion of a single subchannel. In the presented research, the state-of-the-art analysis of turbulent two-phase flows in complex LWR subchannel geometries are demonstrated at both prototypical reactor parameters as well as scaled low pressure conditions. Three different cases are studied, a high-pressure simulation in prototypical reactor subchannel geometry, a low-pressure case in prototypical geometry and a final low-pressure case in a geometry scaled up to conserve the ratio between the bubble size and the domain pitch. Utilizing advanced statistical processing tools, these simulation conditions are compared to shed light on the relevancy of two-phase flow characteristics given the significant differences between LWR and low-pressure conditions. These findings can lead to the generation of useful guiding principles when researchers need to scale the two-phase flow behavior captured at low pressure and temperature conditions to those at reactor operating conditions.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
A GLOBAL TURBULENCE MODEL FOR NEUTRINO-DRIVEN CONVECTION IN CORE-COLLAPSE SUPERNOVAE
Simulations of core-collapse supernovae (CCSNe) result in successful explosions once the neutrino luminosity exceeds a critical curve, and recent simulations indicate that turbulence further enables explosion by reducing this critical neutrino luminosity. We propose a theoretical framework to derive this result and take the first steps by deriving the governing mean-field equations. Using Reynolds decomposition, we decompose flow variables into background and turbulent flows and derive self-consistent averaged equations for their evolution. As basic requirements for the CCSN problem, these equations naturally incorporate steady-state accretion, neutrino heating and cooling, non-zero entropy gradients, and turbulence terms associated with buoyant driving, redistribution, and dissipation. Furthermore, analysis of two-dimensional (2D) CCSN simulations validate these Reynolds-averaged equations, and we show that the physics of turbulence entirely accounts for the differences between 1D and 2D CCSN simulations. As a prelude to deriving the reduction in the critical luminosity, we identify the turbulent terms that most influence the conditions for explosion. Generically, turbulence equations require closure models, but these closure models depend upon the macroscopic properties of the flow. To derive a closure model that is appropriate for CCSNe, we cull the literature for relevant closure models and compare each with 2D simulations. These models employ local closure approximations and fail to reproduce the global properties of neutrino-driven turbulence. Motivated by the generic failure of these local models, we propose an original model for turbulence which incorporates global properties of the flow. This global model accurately reproduces the turbulence profiles and evolution of 2D CCSN simulations.
79 ASTRONOMY AND ASTROPHYSICS↗
Impact of the dynamical core on the direct simulation of tropical cyclones in a high-resolution global model
Our paper examines the impact of the dynamical core on the simulation of tropical cyclone (TC) frequency, distribution, and intensity. The dynamical core, the central fluid flow component of any general circulation model (GCM), is often overlooked in the analysis of a model's ability to simulate TCs compared to the impact of more commonly documented components (e.g., physical parameterizations). The Community Atmosphere Model version 5 is configured with multiple dynamics packages. This analysis demonstrates that the dynamical core has a significant impact on storm intensity and frequency, even in the presence of similar large-scale environments. In particular, the spectral element core produces stronger TCs and more hurricanes than the finite-volume core using very similar parameterization packages despite the latter having a slightly more favorable TC environment. Furthermore, these results suggest that more detailed investigations into the impact of the GCM dynamical core on TC climatology are needed to fully understand these uncertainties. Key Points The impact of the GCM dynamical core is often overlooked in TC assessments The CAM5 dynamical core has a significant impact on TC frequency and intensity A larger effort is needed to better understand this uncertainty
58 GEOSCIENCES↗
Simulation of SPND Beta-Current in AP1000 core
The Self-Powered Neutron Detectors (SPNDs) are widely used in Generation-III commercial PWR due to its small size and external-power free character. It can measure the in-core neutron flux density immediately to obtain the power profile in the core. To better understand the actual state of the simulated core, a SPND beta-current simulation module has been established in the NECP-Bamboo code. Firstly, neutron simulation is carried out in the lattice code Bamboo-Lattice together with nuclide depletion simulation to obtain the depleted composition of the entire lattice including the SPND, and the few-group parameters including density and cross section of emitter material. Secondly, the depleted lattice can then be piped into a Monte-Carlo code for neutron-photon-electron simulation to obtain the electron escape efficiency from the emitter of the SPND to the corresponding collector. Thirdly, the electron escape efficiency can be parameterized into the function of the lattice state parameters similarly to the common few-group constants. Finally, the beta-current was calculated by the core code Bamboo-Core simultaneously with the simulation of the core operation process. The beta-current of the Vanadium SPND in the AP1000 core was simulated and analyzed in this paper. Comparison between simulation and measurement of the beta-current distribution in the specific core state is shown in this paper to verify the effectiveness of the method and the code. The radial normalized current distribution is in good agreement with the measurement. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Enhancements to the Image Analysis Tool for Core Punch Experiments and Simulations (vs. 2014)
A previous paper (Hogden & Unal, 2012, Image Analysis Tool for Core Punch Experiments and Simulations) described an image processing computer program developed at Los Alamos National Laboratory. This program has proven useful so developement has been continued. In this paper we describe enhacements to the program as of 2014.
97 MATHEMATICS AND COMPUTING↗
Localized reactive flow in carbonate rocks: Core-flood experiments and network simulations
We conducted four core-flood experiments on samples of a micritic, reef limestone from Abu Dhabi under conditions of constant flow rate. The pore fluid was water in equilibrium with CO 2 , which, because of its lowered pH, is chemically reactive with the limestone. Flow rates were between 0.03 and 0.1 mL/min. The difference between up and downstream pore pressures dropped to final values <<1 MPa over periods of 3–18 h. Scanning electron microscope and microtomography imaging of the starting material showed that the limestone is mostly calcite and lacks connected macroporosity and that the prevailing pores are few microns large. During each experiment, a wormhole formed by localized dissolution, an observation consistent with the decreases in pressure head between the up and downstream reservoirs. Moreover, we numerically modeled the changes in permeability during the experiments. We devised a network approach that separated the pore space into competing subnetworks of pipes. Thus, the problem was framed as a competition of flow of the reactive fluid among the adversary subnetworks. The precondition for localization within certain time is that the leading subnetwork rapidly becomes more transmissible than its competitors. This novel model successfully simulated features of the shape of the wormhole as it grew from few to about 100 µm, matched the pressure history patterns, and yielded the correct order of magnitude of the breakthrough time. Finally, we systematically studied the impact of changing the statistical parameters of the subnetworks. Larger mean radius and spatial correlation of the leading subnetwork led to faster localization.
58 GEOSCIENCES↗
Anisotropic formation of magnetized cores in turbulent clouds
In giant molecular clouds (GMCs), shocks driven by converging turbulent flows create high-density, strongly magnetized regions that are locally sheetlike. In previous work, we showed that within these layers, dense filaments and embedded self-gravitating cores form by gathering material along the magnetic field lines. Here, we extend the parameter space of our three-dimensional, turbulent MHD core formation simulations. We confirm the anisotropic core formation model we previously proposed and quantify the dependence of median core properties on the pre-shock inflow velocity and upstream magnetic field strength. Our results suggest that bound core properties are set by the total dynamic pressure (dominated by large-scale turbulence) and thermal sound speed c{sub s} in GMCs, independent of magnetic field strength. For models with a Mach number between 5 and 20, the median core masses and radii are comparable to the critical Bonnor–Ebert mass and radius defined using the dynamic pressure for P{sub ext}. Our results correspond to M{sub core}=1.2c{sub s}{sup 4}(G{sup 3}ρ{sub 0}v{sub 0}{sup 2}){sup −1/2} and R{sub core}=0.34c{sub s}{sup 2}(Gρ{sub 0}v{sub 0}{sup 2}){sup −1/2} for ρ{sub 0} and v{sub 0}, the large-scale mean density and velocity. For our parameter range, the median M{sub core}∼0.1−1M{sub ⊙}, but a very high pressure cloud could have lower characteristic core mass. We find cores and filaments form simultaneously, and filament column densities are a factor of ∼2 greater than the surrounding cloud when cores first collapse. We also show that cores identified in our simulations have physical properties comparable to those observed in the Perseus cloud. Superthermal cores in our models are generally also magnetically supercritical, suggesting that the same may be true in observed clouds.
79 ASTRONOMY AND ASTROPHYSICS↗
Fully compressible simulations of waves and core convection in main-sequence stars
Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.
79 ASTRONOMY AND ASTROPHYSICS↗
Simulation of in-core neutron noise measurements for axial void profile reconstruction in boiling water reactors
A possibility to reconstruct the axial void profile from the simulated in-core neutron noise which is caused by density fluctuations in a Boiling Water Reactor (BWR) heated channel is considered. For this purpose, a self-contained model of the two-phase flow regime is constructed which has quantitatively and qualitatively similar properties to those observed in real BWRs. The model is subsequently used to simulate the signals of neutron detectors induced by the corresponding perturbations in the flow density. The bubbles are generated randomly in both space and time using Monte-Carlo techniques. The axial distribution of the bubble production is chosen such that the mean axial void fraction and void velocity follow the actual values of BWRs. The induced neutron noise signals are calculated and then processed by the standard signal analysis methods such as Auto-Power Spectral Density (APSD) and Cross-Power Spectral Density (CPSD). Two methods for axial void and velocity profiles reconstruction are discussed: the first one is based on the change of the break frequency of the neutron auto-power spectrum with axial core elevation, while the second refers to the estimation of transit times of propagating steam fluctuations between different axial detector positions. This paper summarizes the principles of the model and presents a numerical testing of the qualitative applicability to estimate the required parameters for the reconstruction of the void fraction profile from the neutron noise measurements. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Sample-morphology effects on x-ray photoelectron peak intensities. III. Simulated spectra of model core–shell nanoparticles
The National Institute of Standards and Technology database for the simulation of electron spectra for surface analysis has been used to simulate Cu 2p photoelectron spectra for four types of spherical copper–gold nanoparticles (NPs). These simulations were made to extend the work of Tougaard [J. Vac. Sci. Technol. A 14, 1415 (1996)] and of Powell et al. [J. Vac. Sci. Technol. A 31, 021402 (2013)] who performed similar simulations for four types of planar copper–gold films. The Cu 2p spectra for the NPs were compared and contrasted with analogous results for the planar films and the effects of elastic scattering were investigated. The new simulations were made for a monolayer of three types of Cu/Au core–shell NPs on a Si substrate: (1) an Au shell of variable thickness on a Cu core with diameters of 0.5, 1.0, 2.0, 5.0, and 10.0 nm; (2) a Cu shell of variable thickness on an Au core with diameters of 0.5, 1.0, 2.0, 5.0, and 10.0 nm; and (3) an Au shell of variable thickness on a 1 nm Cu shell on an Au core with diameters of 0.5, 1.0, 2.0, 5.0, and 10.0 nm. For these three morphologies, the outer-shell thickness was varied until the Cu 2p{sub 3/2} peak intensity was the same (within 2%) as that found in our previous work with planar Cu/Au morphologies. The authors also performed similar simulations for a monolayer of spherical NPs consisting of a CuAu{sub x} alloy (also on a Si substrate) with diameters of 0.5, 1.0, 2.0, 5.0, and 10.0 nm. In the latter simulations, the relative Au concentration (x) was varied to give the same Cu 2p{sub 3/2} peak intensity (within 2%) as that found previously. For each morphology, the authors performed simulations with elastic scattering switched on and off. The authors found that elastic-scattering effects were generally strong for the Cu-core/Au-shell and weak for the Au-core/Cu-shell NPs; intermediate elastic-scattering effects were found for the Au-core/Cu-shell/Au-shell NPs. The shell thicknesses required to give the selected Cu 2p{sub 3/2} peak intensity for the three types of core–shell NPs were less than the corresponding film thicknesses of planar samples since Cu 2p photoelectrons can be detected from the sides and, for the smaller NPs, bottoms of the NPs. Elastic-scattering effects were also observed on the Au atomic fractions found for the CuAu{sub x} NP alloys with different diameters.
36 MATERIALS SCIENCE↗
COOL CORE CLUSTERS FROM COSMOLOGICAL SIMULATIONS
We present results obtained from a set of cosmological hydrodynamic simulations of galaxy clusters, aimed at comparing predictions with observational data on the diversity between cool-core (CC) and non-cool-core (NCC) clusters. Our simulations include the effects of stellar and active galactic nucleus (AGN) feedback and are based on an improved version of the smoothed particle hydrodynamics code GADGET-3, which ameliorates gas mixing and better captures gas-dynamical instabilities by including a suitable artificial thermal diffusion. In this Letter, we focus our analysis on the entropy profiles, the primary diagnostic we used to classify the degree of cool-coreness of clusters, and the iron profiles. In keeping with observations, our simulated clusters display a variety of behaviors in entropy profiles: they range from steadily decreasing profiles at small radii, characteristic of CC systems, to nearly flat core isentropic profiles, characteristic of NCC systems. Using observational criteria to distinguish between the two classes of objects, we find that they occur in similar proportions in both simulations and observations. Furthermore, we also find that simulated CC clusters have profiles of iron abundance that are steeper than those of NCC clusters, which is also in agreement with observational results. We show that the capability of our simulations to generate a realistic CC structure in the cluster population is due to AGN feedback and artificial thermal diffusion: their combined action allows us to naturally distribute the energy extracted from super-massive black holes and to compensate for the radiative losses of low-entropy gas with short cooling time residing in the cluster core.
79 ASTRONOMY AND ASTROPHYSICS↗
PROGENITOR-DEPENDENT EXPLOSION DYNAMICS IN SELF-CONSISTENT, AXISYMMETRIC SIMULATIONS OF NEUTRINO-DRIVEN CORE-COLLAPSE SUPERNOVAE
We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11–28 M {sub ⊙}, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si–O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si–O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection timescales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.
79 ASTRONOMY AND ASTROPHYSICS↗
Micro Reactor Simulation Using the PROTEUS Suite in FY19
Two micro reactor cores were simulated using the MOC solver of PROTEUS: MegaPower for a heat pipe micro reactor and Holos-Quad for a high-temperature gas-cooled micro reactor. Verification of PROTEUS standalone was performed by simulating 2D and 3D core benchmark problems developed based on MegaPower and Holos-Quad. Simulation results indicated that PROTEUS successfully modelled 2D and 3D cores of the MegaPower and Holos-Quad micro reactors and PROTEUS eigenvalues agreed well with the SERPENT Monte Carlo solutions. For a heat pipe cooled micro reactor, multiphysics simulations were performed using PROTEUS / FLUENT / ANLHTP. The coupling system of the three codes were developed using Python-based drivers, which exchange the data (power, temperature, and heat rate) required for the coupling and control the boundary conditions of FLUENT and ANLHTP. The coupled system was qualitatively verified using a 3D unit assembly problem composed of six fuel rods and seven heat pipes, which was developed based on the MegaPower compositions. The steady-state problem was solved using FLUENT / ANLHTP with given asymmetric fuel powers, demonstrating the reasonable convergence of temperatures of 7 heat pipes. Transient simulations were also performed by making one out of 7 heat pipes fail. For the transient, the coupled system of PROTEUS / FLUENT / ANLHTP demonstrated reasonable changes of total power and heat pipe temperatures with time, accounting for temperature feedback effects. In this simulation, no effect of structural change was accounted for, which will be performed using ANSYS-mechanical in future along with further verification with larger or whole-core 3D problems.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Rhapsody-G simulations I: the cool cores, hot gas and stellar content of massive galaxy clusters
We present the rhapsody-g suite of cosmological hydrodynamic zoom simulations of 10 massive galaxy clusters at the M vir ~10 15 M ⊙ scale. These simulations include cooling and subresolution models for star formation and stellar and supermassive black hole feedback. The sample is selected to capture the whole gamut of assembly histories that produce clusters of similar final mass. We present an overview of the successes and shortcomings of such simulations in reproducing both the stellar properties of galaxies as well as properties of the hot plasma in clusters. In our simulations, a long-lived cool-core/non-cool-core dichotomy arises naturally, and the emergence of non-cool cores is related to low angular momentum major mergers. Nevertheless, the cool-core clusters exhibit a low central entropy compared to observations, which cannot be alleviated by thermal active galactic nuclei feedback. For cluster scaling relations, we find that the simulations match well the M 500 –Y 500 scaling of Planck Sunyaev–Zeldovich clusters but deviate somewhat from the observed X-ray luminosity and temperature scaling relations in the sense of being slightly too bright and too cool at fixed mass, respectively. Stars are produced at an efficiency consistent with abundance-matching constraints and central galaxies have star formation rates consistent with recent observations. In conclusion, while our simulations thus match various key properties remarkably well, we conclude that the shortcomings strongly suggest an important role for non-thermal processes (through feedback or otherwise) or thermal conduction in shaping the intracluster medium.
79 ASTRONOMY AND ASTROPHYSICS↗
Kernel fusion in atomistic spin dynamics simulations on Nvidia GPUs using tensor core
In atomistic spin dynamics simulations, the time cost of constructing the space- and time-displaced pair correlation function in real space increases quadratically as the number of spins N, leading to significant computational effort. The GEMM subroutine can be adopted to accelerate the calculation of the dynamical spin-spin correlation function, but the computational cost of simulating large spin systems (>40000 spins) on CPUs remains expensive. In this work, we perform the simulation on the graphics processing unit (GPU), a hardware solution widely used as an accelerator for scientific computing and deep learning. Here we show that GPUs can accelerate the simulation up to 25-fold compared to multi-core CPUs when using the GEMM subroutine on both. To hide memory latency, we fuse the element-wise operation into the GEMM kernel using CUTLASS that can improve the performance by 26% ~ 33% compared to implementation based on cuBLAS. Furthermore, we perform the on-the-fly calculation in the epilogue of the GEMM subroutine to avoid saving intermediate results on global memory, which makes the large-scale atomistic spin dynamics simulation feasible and affordable.
97 MATHEMATICS AND COMPUTING↗
Study on efficiency of time computation in x-ray imaging simulation base on Monte Carlo algorithm using graphics processing unit
Monte Carlo (MC) is one of the powerful techniques for simulation in x-ray imaging. MC method can simulate the radiation transport within matter with high accuracy and provides a natural way to simulate radiation transport in complex systems. One of the codes based on MC algorithm that are widely used for radiographic images simulation is MC-GPU, a codes developed by Andrea Basal. This study was aimed to investigate the time computation of x-ray imaging simulation in GPU (Graphics Processing Unit) compared to a standard CPU (Central Processing Unit). Furthermore, the effect of physical parameters to the quality of radiographic images and the comparison of image quality resulted from simulation in the GPU and CPU are evaluated in this paper. The simulations were run in CPU which was simulated in serial condition, and in two GPU with 384 cores and 2304 cores. In simulation using GPU, each cores calculates one photon, so, a large number of photon were calculated simultaneously. Results show that the time simulations on GPU were significantly accelerated compared to CPU. The simulations on the 2304 core of GPU were performed about 64 -114 times faster than on CPU, while the simulation on the 384 core of GPU were performed about 20 – 31 times faster than in a single core of CPU. Another result shows that optimum quality of images from the simulation was gained at the history start from 10{sup 8} and the energy from 60 Kev to 90 Kev. Analyzed by statistical approach, the quality of GPU and CPU images are relatively the same.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A GPU‐Based Ocean Dynamical Core for Routine Mesoscale‐Resolving Climate Simulations
Abstract We describe an ocean hydrostatic dynamical core implemented in Oceananigans optimized for Graphical Processing Unit (GPU) architectures. On 64 A100 GPUs, equivalent to 16 computational nodes in current state‐of‐the‐art supercomputers, our dynamical core can simulate a decade of near‐global ocean dynamics per wall‐clock day at an 8‐km horizontal resolution; a resolution adequate to resolve the ocean's mesoscale eddy field. Such efficiency, achieved with relatively modest hardware resources, suggests that climate simulations on GPUs can incorporate fully eddy‐resolving ocean models. This removes a major source of systematic bias in current IPCC coupled model projections, the parameterization of ocean eddies, and represents a major advance in climate modeling. We discuss the computational strategies, focusing on GPU‐specific optimization and numerical implementation details that enable such high performance.
Silvestri, Simone [Massachusetts Institute of Tech↗