Engineering Papers⌕ Search

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.

At least 415 records · Page 23

Analysis of Lunar Highland Regolith Samples From Apollo 16 Drive Core 64001/2 and Lunar Regolith Simulants - an Expanding Comparative Database

This slide presentation reviews the work to analyze the lunar highland regolith samples that came from the Apollo 16 core sample 64001/2 and simulants of lunar regolith, and build a comparative database. The work is part of a larger effort to compile an internally consistent database on lunar regolith (Apollo Samples) and lunar regolith simulants. This is in support of a future lunar outpost. The work is to characterize existing lunar regolith and simulants in terms of particle type, particle size distribution, particle shape distribution, bulk density, and other compositional characteristics, and to evaluate the regolith simulants by the same properties in comparison to the Apollo sample lunar regolith.

Schrader, Christian M.↗

Simulation of the full-core pin-model by JMCT Monte Carlo neutron-photon transport code

Since the large numbers of cells over a million, the tallies over a hundred million and the particle histories over ten billion, the simulation of the full-core pin-by-pin model has become a real challenge for the computers and the computational methods. On the other hand, the basic memory of the model has exceeded the limit of a single CPU, so the spatial domain and data decomposition must be considered. JMCT (J Monte Carlo Transport code) has successful fulfilled the simulation of the full-core pin-by-pin model by the domain decomposition and the nested parallel computation. The k{sub eff} and flux of each cell are obtained. (authors)

97 MATHEMATICS AND COMPUTING↗

Simulating High Flux Isotope Reactor Core Thermal-Hydraulics via Interdimensional Model Coupling

A coupled interdimensional model is presented for the simulation of the thermal-hydraulic characteristics of the High Flux Isotope Reactor core at Oak Ridge National Laboratory. The model consists of two domains a solid involute fuel plate and the surrounding liquid coolant channel. The fuel plate is modeled explicitly in three-dimensions. The coolant channel is approximated as a twodimensional slice oriented perpendicular to the fuel plate s surface. The two dimensionally-inconsistent domains are linked to one another via interdimensional model coupling mechanisms. The coupled model is presented as a simplified alternative to a fully explicit, fully three-dimensional model. Involute geometries were constructed in SolidWorks. Derivations of the involute construction equations are presented. Geometries were then imported into COMSOL Multiphysics for simulation and modeling. Both models are described in detail so as to highlight their respective attributes in the 3D model, the pursuit of an accurate, reliable, and complete solution; in the coupled model, the intent to simplify the modeling domain as much as possible without affecting significant alterations to the solution. The coupled model was created with the goal of permitting larger portions of the reactor core to be modeled at once without a significant sacrifice to solution integrity. As such, particular care is given to validating incorporated model simplifications. To the greatest extent possible, the decrease in solution time as well as computational cost are quantified versus the effects such gains have on the solution quality. A variant of the coupled model which sufficiently balances these three solution characteristics is presented alongside the more comprehensive 3D model for comparison and validation.

COMSOL HFIR LEU↗

Large-scale Multiphysics Simulations of Small Modular Reactors Operating in Natural Circulation

Thanks to the advancements in high-performance computing, advanced modeling and simulation have become crucial in driving the development and deployment of next-generation nuclear reactors, such as small modular reactors (SMRs). SMRs offer the promise of cost-effective baseload electricity production and improved safety, while addressing some of the challenges associated with large reactor designs, such as high capital costs and extended construction timelines. As part of the Exascale Computing Project, the large-scale multiphysics simulation of an entire SMR primary system has been achieved by combining computational fluid dynamics and neutronics. In addition to the successful demonstration of full-core SMR simulations, the current study integrated the impact of natural circulation into the system. Natural circulation is the primary mechanism driving coolant circulation in SMRs. The mass flow rate in the core depends on the core power, and a numerical model has been developed to predict it. The pressure drop caused by the helical coil steam generator was also accounted for by developing a pressure drop correlation based on high-fidelity large eddy simulation results, further improving prediction accuracy. In conclusion, the results of the study demonstrate that the implemented natural circulation model is effective in predicting the responses of SMR full-core multiphysics simulations.

ECP↗

Accelerating Our Understanding of Supernova Explosion Mechanism via Simulations and Visualizations with GenASiS

Core-collapse supernovae are among the most powerful explosions in the Universe, releasing about 1053 erg of energy on timescales of a few tens of seconds. These explosion events are also responsible for the production and dissemination of most of the heavy elements, making life as we know it possible. Yet exactly how they work is still unresolved. One reason for this is the sheer complexity and cost of a self-consistent, multi-physics, and multi-dimensional core-collapse supernova simulation, which is impractical, and often impossible, even on the largest supercomputers we have available today. To advance our understanding we instead must often use simplified models, teasing out the most important ingredients for successful explosions, while helping us to interpret results from higher fidelity multi-physics models. In this paper we investigate the role of instabilities in the core-collapse supernova environment. We present here simulation and visualization results produced by our code GenASiS.

Budiardja, R. D.↗

Equation of state effects in the core collapse of a 20 – M ⊙ star

Uncertainties in our knowledge of the properties of dense matter near and above nuclear saturation density are among the main sources of variations in multimessenger signatures predicted for core-collapse supernovae (CCSNe) and the properties of neutron stars (NSs). We construct 97 new finite-temperature equations of state (EOSs) of dense matter that obey current experimental, observational, and theoretical constraints and discuss how systematic variations in the EOS parameters affect the properties of cold nonrotating NSs and the core collapse of a 20–M ⊙ progenitor star. The core collapse of the 20–M ⊙ progenitor star is simulated in spherical symmetry using the general-relativistic radiation-hydrodynamics code GR1D where neutrino interactions are computed for each EOS using the NuLib library. We conclude that the effective mass of nucleons at densities above nuclear saturation density is the largest source of uncertainty in the CCSN neutrino signal and dynamics even though it plays a subdominant role in most properties of cold NS matter. Meanwhile, changes in other observables affect the properties of cold NSs, while having little effect in CCSNe. To strengthen our conclusions, we perform six octant three-dimensional CCSN simulations varying the effective mass of nucleons at nuclear saturation density. Here, we conclude that neutrino heating and, thus, the likelihood of explosion is significantly increased for EOSs where the effective mass of nucleons at nuclear saturation density is large.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

RELAP5-3D Modeling of High Temperature Test Facility (HTTF)

HTTF at Oregon State University (OSU) Reference: General Atomics? modular high-temperature gas-cooled reactor Helium cooled, electrically heated Prismatic graphite blocks in the core and reflectors Alumina ceramic blocks are used to simulate the core and top and bottom reflectors One-fourth scale in length and diameter Most of the coolant channels in the core are full scale Lower pressure compared to the prototype reactor Over 500 instruments Designed primarily to investigate depressurized conduction cooldown (DCC) transients

97 MATHEMATICS AND COMPUTING↗

RELAP5-3D Modeling of the OECD-NEA HTTF Benchmark

HTTF at Oregon State University (OSU) Reference: General Atomics’ modular high-temperature gas-cooled reactor Helium cooled, electrically heated (2.2 MW) Prismatic graphite blocks in the core and reflectors Alumina ceramic blocks are used to simulate the core and top and bottom reflectors One-fourth scale in length and diameter Most of the coolant channels in the core are full scale Lower pressure compared to the prototype reactor (0.7 MPa) Over 500 instruments Designed primarily to investigate depressurized (DCC) and pressurized (PCC) conduction cooldown transients

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Exploring the effects of a nonhydrostatic dynamical core in high-resolution aquaplanet simulations

This study explores the impact of a nonhydrostatic dynamical core in high-resolution regional climate simulations using an aquaplanet framework. The Weather Research and Forecasting (WRF) model is used to conduct simulations with both hydrostatic (H) and nonhydrostatic (NH) solvers at horizontal grid spacings (?x) of 36, 12, and 4 km. The differences between the H and NH simulated precipitation (?P) are notable even at ?x = 12 km in the inter-tropical convergence zone and the transition region to the drier subtropics. At gray-zone grid spacing (12 km and 4 km) over the tropics, ?P is sensitive to whether a cumulus parameterization scheme is used or not. With an idealized Witch of Agnesi land mountain, differences in the precipitation and circulation (vertical velocity) between the H and NH simulations are significant even at ?x = 36 km in the tropics due largely to the strong feedbacks related to moist processes. The differences increase as the model grid spacing and mountain half width (a) are reduced, accompanied by a shift towards a more nonhydrostatic flow regime at a = 24 km. Latent heat release drastically enhances the differences between the NH and H simulations and extends the effect of nonhydrostatic dynamics to a broader region over the mountain and downstream over the tropics. Overall, differences exist between H and NH simulations even at resolutions between 12 and 36 km, but the differences are sensitive to the representations of moist physics and other features such as horizontal diffusion used in the WRF model.

nonhydrostatic dynamics, latent heating, nonhydros↗

On the Origin of Pulsar and Magnetar Magnetic Fields

Abstract In order to address the generation of neutron star magnetic fields, with particular focus on the dichotomy between magnetars and radio pulsars, we consider the properties of dynamos as inferred from other astrophysical systems. With sufficiently low (modified) Rossby number, convective dynamos are known to produce dipole-dominated fields whose strength scales with convective flux, and we argue that these expectations should apply to the convective protoneutron stars (PNSs) at the centers of core-collapse supernovae. We analyze a suite of three-dimensional simulations of core collapse, featuring a realistic equation of state and full neutrino transport, in this context. All our progenitor models, ranging from 9 M ⊙ to 25 M ⊙ , including one with initial rotation, have sufficiently vigorous PNS convection to generate dipole fields of order ∼10 15 Gauss, if the modified Rossby number resides in the critical range. Thus, the magnetar/radio pulsar dichotomy may arise naturally in part from the distribution of core rotation rates in massive stars.

79 ASTRONOMY AND ASTROPHYSICS↗

PWR core and spent fuel pool analysis using scale and nestle

The SCALE nuclear analysis code system [SCALE, 2011], developed and maintained at Oak Ridge National Laboratory (ORNL) is widely recognized as high quality software for analyzing nuclear systems. The SCALE code system is composed of several validated computer codes and methods with standard control sequences, such as the TRITON/NEWT lattice physics sequence, which supplies dependable and accurate analyses for industry, regulators, and academia. Although TRITON generates energy-collapsed and space-homogenized few group cross sections, SCALE does not include a full-core nodal neutron diffusion simulation module within. However, in the past few years, the open-source NESTLE core simulator [NESTLE, 2003], originally developed at North Carolina State Univ. (NCSU), has been updated and upgraded via collaboration between ORNL and the Univ. of Tennessee (UT), so it now has a growingly seamless coupling to the TRITON/NEWT lattice physics [Galloway, 2010]. This study presents the methodology used to couple lattice physics data between TRITON and NESTLE in order to perform a three-dimensional full-core analysis employing a 'real-life' Duke Energy PWR as the test bed. The focus for this step was to compare the key parameters of core reactivity and radial power distribution versus plant data. Following the core analysis, following a three cycle burn, a spent fuel pool analysis was done using information generated from NESTLE for the discharged bundles and was compared to Duke Energy spent fuel pool models. The KENO control module from SCALE was employed for this latter stage of the project. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Numerical study of core formation of asymmetrically driven cone-guided targets

Compression of a directly driven fast ignition cone-sphere target with a finite number of laser beams is numerically studied using a three-dimensional hydrodynamics code IMPACT-3D. The formation of a dense plasma core is simulated for 12-, 9-, 6-, and 4-beam configurations of the GEKKO XII laser. The complex 3D shapes of the cores are analyzed by elucidating synthetic 2D x-ray radiographic images in two orthogonal directions. Finally, the simulated x-ray images show significant differences in the core shape between the two viewing directions and rotation of the stagnating core axis in the top view for the axisymmetric 9- and 6-beam configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ON THE IMPACT OF THREE DIMENSIONS IN SIMULATIONS OF NEUTRINO-DRIVEN CORE-COLLAPSE SUPERNOVA EXPLOSIONS

We present one-dimensional (1D), two-dimensional (2D), and three-dimensional (3D) hydrodynamical simulations of core-collapse supernovae including a parameterized neutrino heating and cooling scheme in order to investigate the critical core neutrino luminosity (L{sub crit}) required for explosion. In contrast to some previous works, we find that 3D simulations explode later than 2D simulations, and that L{sub crit} at fixed mass accretion rate is somewhat higher in three dimensions than in two dimensions. We find, however, that in two dimensions L{sub crit} increases as the numerical resolution of the simulation increases. In contrast to some previous works, we argue that the average entropy of the gain region is in fact not a good indicator of explosion but is rather a reflection of the greater mass in the gain region in two dimensions. We compare our simulations to semi-analytic explosion criteria and examine the nature of the convective motions in two dimensions and three dimensions. We discuss the balance between neutrino-driven buoyancy and drag forces. In particular, we show that the drag force will be proportional to a buoyant plume's surface area while the buoyant force is proportional to a plume's volume and, therefore, plumes with greater volume-to-surface-area ratios will rise more quickly. We show that buoyant plumes in two dimensions are inherently larger, with greater volume-to-surface-area ratios, than plumes in three dimensions. In the scenario that the supernova shock expansion is dominated by neutrino-driven buoyancy, this balance between buoyancy and drag forces may explain why 3D simulations explode later than 2D simulations and why L{sub crit} increases with resolution. Finally, we provide a comparison of our results with other calculations in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗

A new parameter-free soft-core potential for silica and its application to simulation of silica anomalies

A core-softening of the effective interaction between oxygen atoms in water and silica systems and its role in developing anomalous thermodynamic, transport, and structural properties have been extensively debated. For silica, the progress with addressing these issues has been hampered by a lack of effective interaction models with explicit core-softening. In this work, we present an extension of a two-body soft-core interatomic force field for silica recently reported by us [S. Izvekov and B. M. Rice, J. Chem. Phys. 136(13), 134508 (2012)] to include three-body forces. Similar to two-body interaction terms, the three-body terms are derived using parameter-free force-matching of the interactions from ab initio MD simulations of liquid silica. The derived shape of the O–Si–O three-body potential term affirms the existence of repulsion softening between oxygen atoms at short separations. The new model shows a good performance in simulating liquid, amorphous, and crystalline silica. By comparing the soft-core model and a similar model with the soft-core suppressed, we demonstrate that the topology reorganization within the local tetrahedral network and the O–O core-softening are two competitive mechanisms responsible for anomalous thermodynamic and kinetic behaviors observed in liquid and amorphous silica. The studied anomalies include the temperature of density maximum locus and anomalous diffusivity in liquid silica, and irreversible densification of amorphous silica. We show that the O–O core-softened interaction enhances the observed anomalies primarily through two mechanisms: facilitating the defect driven structural rearrangements of the silica tetrahedral network and modifying the tetrahedral ordering induced interactions toward multiple characteristic scales, the feature which underlies the thermodynamic anomalies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiscale thermal hydraulic coupling methods for boiling water reactor simulation

During the last two years, the Virtual Environment for Reactor Applications (VERA) has been extended to simulate boiling water reactors (BWRs). The thermal hydraulic effects present in BWRs are far more complex than those in pressurized water reactors (PWRs). Therefore, the runtime is significantly increase and convergence behavior worse compared to PWRs. Most of the additional runtime is spent during the first few coupled iterations, when the power shape is still rapidly evolving, dramatically affecting the thermal hydraulics (TH). To alleviate the increased computational expense, a multiscale TH coupling approach was developed in VERA in which a highly efficient, simplified TH model is solved for several coupled iterations until the power shape is partially converged. The simplified solution is used to inform the assembly-wise flow distribution in the high-fidelity TH solver CTF, reducing the amount of work required to properly balance the pressure drop in each channel and reducing runtime for the couple calculation. Because the final TH calculations are performed with CTF, there is ultimately no impact on the accuracy of the converged solution. This paper presents the details of this multiscale TH coupling approach, along with results for a 4x4 array of GE-14 fuel bundles and whole-core coupled simulations of the Hatch core. These two cases show that the multiscale approach dramatically reduces the runtime of coupled BWR simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗