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 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

New reactor cavity cooling system having passive safety features using novel shape for HTGRs and VHTRs

A new, highly efficient reactor cavity cooling system (RCCS) with passive safety features without a requirement for electricity and mechanical drive is proposed for high temperature gas cooled reactors (HTGRs) and very high temperature reactors (VHTRs). The RCCS design consists of continuous closed regions; one is an ex-reactor pressure vessel (RPV) region and another is a cooling region having heat transfer area to ambient air assumed at 40 (°C). The RCCS uses a novel shape to efficiently remove the heat released from the RPV with radiation and natural convection. Employing the air as the working fluid and the ambient air as the ultimate heat sink, the new RCCS design strongly reduces the possibility of losing the heat sink for decay heat removal. Therefore, HTGRs and VHTRs adopting the new RCCS design can avoid core melting due to overheating the fuels. The simulation results from a commercial CFD code, STAR-CCM+, show that the temperature distribution of the RCCS is within the temperature limits of the structures, such as the maximum operating temperature of the RPV, 713.15 (K) = 440 (°C), and the heat released from the RPV could be removed safely, even during a loss of coolant accident (LOCA). Finally, when the RCCS can remove 600 (kW) of the rated nominal state even during LOCA, the safety review for building the HTTR could confirm that the temperature distribution of the HTTR is within the temperature limits of the structures to secure structures and fuels after the shutdown because the large heat capacity of the graphite core can absorb heat from the fuel in a short period. Therefore, the capacity of the new RCCS design would be sufficient for decay heat removal.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRSDate not supplied

A Computational Study of Plume Modeling For Space Launch System Abort Scenarios

This extended abstract details work on Space Launch System (SLS) core stage drag during abort scenarios that is currently being conducted. In the final paper, viscous computational fluid dynamic simulations will be used to study the axial force coefficient on the SLS core stage as a function of its distance to the upstream crew module and launch abort system, the freestream conditions, and the chemical model used to represent the launch abort system abort motor and attitude control motor plumes. In addition, the computational expense of each method will be studied in order to help determine the appropriate trade-off between fidelity and simulation turn around time.

CFD↗
DOE OSTI2016

Benchmark Progress towards TREAT M8CAL Core Loading to Support Validation of TREAT Operations

Benchmark experiment data is necessary to validate modeling and simulation activities to support Transient Reactor Experiment Test (TREAT) Facility restart and operations. Key measurements of interest include criticality, control rod worth, excess reactivity, and shutdown margin for varying core loading sizes. Benchmark evaluations are being developed according the guidelines provided in the International Handbook of Evaluated Reactor Physics Benchmark Experiments (IRPhEP Handbook) that can also support advanced modeling and simulation activities, new experiment design for advanced reactors and accident tolerant fuel, and potential low-enriched uranium conversion of the TREAT reactor. This summary discusses development activities supporting a benchmark evaluation of the TREAT M8CAL core loading, which will be submitted for peer-review and publication in the IRPhEP Handbook with the identifier TREAT-FUND-RESR-002. Detailed models of TREAT were developed and evaluated to prepare benchmark models of the minimum critical mass configuration. Calculations were performed using MCNP6.1 with ENDF/B-VII.1 nuclear data libraries. Discussion of results.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

INFLUENCE OF MAGNETOROTATIONAL INSTABILITY ON NEUTRINO HEATING: A NEW MECHANISM FOR WEAKLY MAGNETIZED CORE-COLLAPSE SUPERNOVAE

We investigated the impact of magnetorotational instability (MRI) on the dynamics of weakly magnetized, rapidly rotating core-collapse supernovae by conducting high-resolution axisymmetric MHD simulations with simplified neutrino transfer. We found that an initially sub-magnetar-class magnetic field is drastically amplified by MRI and substantially affects the dynamics thereafter. Although the magnetic pressure is not strong enough to eject matter, the amplified magnetic field efficiently transfers angular momentum from small to large radii and from higher to lower latitudes, which causes the expansion of the heating region due to the extra centrifugal force. This then enhances the efficiency of neutrino heating and eventually leads to neutrino-driven explosion. This is a new scenario of core-collapse supernovae that has never been demonstrated by past numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Interaction of external n = 1 magnetic fields with the sawtooth instability in low- q RFX-mod and DIII-D tokamaks

External n = 1 magnetic fields are applied in RFX-mod and DIII-D low safety factor Tokamak plasmas to investigate their interaction with the internal MHD dynamics and in particular with the sawtooth instability. In these experiments the applied magnetic fields cause a reduction of both the sawtooth amplitude and period, leading to an overall stabilizing effect on the oscillations. In RFX-mod sawteeth eventually disappear and are replaced by a stationary m = 1, n = 1 helical equilibrium without an increase in disruptivity. However toroidal rotation is significantly reduced in these plasmas, thus it is likely that the sawtooth mitigation in these experiments is due to the combination of the helically deformed core and the reduced rotation. The former effect is qualitatively well reproduced by nonlinear MHD simulations performed with the PIXIE3D code. The results obtained in these RFX-mod experiments motivated similar ones in DIII-D L-mode diverted Tokamak plasmas at low q 95 . These experiments succeeded in reproducing the sawtooth mitigation with the approach developed in RFX-mod. In DIII-D this effect is correlated with a clear increase of the n = 1 plasma response, that indicates an enhancement of the coupling to the marginally stable n = 1 external kink, as simulations with the linear MHD code IPEC suggest. A significant rotation braking in the plasma core is also observed in DIII-D. Finally, numerical calculations of the neoclassical toroidal viscosity (NTV) carried out with PENT identify this torque as a possible contributor for this effect.

3D magnetic fields↗
DOE OSTI2017

Development of High-Performance Cast Crankshafts. Final Technical Report

The objective of this project was to develop technologies that would enable the production of cast crankshafts that can replace high performance forged steel crankshafts. To achieve this, the Ultimate Tensile Strength (UTS) of the new material needs to be 850 MPa with a desired minimum Yield Strength (YS; 0.2% offset) of 615 MPa and at least 10% elongation. Perhaps more challenging, the cast material needs to be able to achieve sufficient local fatigue properties to satisfy the durability requirements in today’s high performance gasoline and diesel engine applications. The project team focused on the development of cast steel alloys for application in crankshafts to take advantage of the higher stiffness over other potential material choices. The material and process developed should be able to produce high-performance crankshafts at no more than 110% of the cost of current production cast units, perhaps the most difficult objective to achieve. To minimize costs, the primary alloy design strategy was to design compositions that can achieve the required properties with minimal alloying and post-casting heat treatments. An Integrated Computational Materials Engineering (ICME) based approach was utilized, rather than relying only on traditional trial-and-error methods, which has been proven to accelerate alloy development time. Prototype melt chemistries designed using ICME were cast as test specimens and characterized iteratively to develop an alloy design within a stage-gate process. Standard characterization and material testing was done to validate the alloy performance against design targets and provide feedback to material design and manufacturing process models. Finally, the project called for Caterpillar and General Motors (GM) to develop optimized crankshaft designs using the final material and manufacturing processing path developed. A multi-disciplinary effort was to integrate finite element analyses by engine designers and geometry-specific casting simulations with existing materials models to optimize crankshaft cost and performance. Prototype crankshafts of the final design were to be produced and validated using laboratory bench testing and on-engine durability testing. ICME process simulation tools were used to investigate a broad range of processing concepts. These concepts included casting orientation, various mold and core materials, and various filling and feeding strategies. Each crankshaft was first simulated without gating and risers, which is termed natural solidification. The natural solidification results were used as a baseline for strategy development of each concept. Casting process simulations and ICME tools were proven to be reasonable predictors of real world results. Potential alloys were developed that could meet the project material property goals with appropriate normalization and temper treatments. For the alloys considered, post-normalization temper treatments proved to be necessary to achieve the desired yield strengths and elongations and appropriate heat treatments were designed using ICME tools. The experimental data of all the alloys were analyzed in combination with ICME tools to establish chemistry-process-structure relations. Several GM small gas engine (SGE) crankshafts were successfully cast in sand molds using two different sprue, runner, gate, riser, chill designs. These crankshafts were cast in two different steel alloys developed during the project, but casting finishing (e.g. riser removal) remains a cost challenge. A long list of future work was left unfinished when this project was unexpectedly terminated.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Benchmark Development for TREAT Minimum Critical Mass Core Loading to Support Validation of TREAT Operations

Benchmark experiment data is necessary to validate modeling and simulation activities to support Transient Reactor Experiment Test (TREAT) Facility restart and operations. Key measurements of interest include criticality, control rod worth, excess reactivity, and shutdown margin for varying core loading sizes. Benchmark evaluations are being developed according the guidelines provided in the International Handbook of Evaluated Reactor Physics Benchmark Experiments (IRPhEP Handbook) that can also support advanced modeling and simulation activities, new experiment design for advanced reactors and accident tolerant fuel, and potential low-enriched uranium conversion of the TREAT reactor. This summary discusses development activities supporting a benchmark evaluation of the TREAT minimum critical mass core loading, which will be submitted for peer-review and publication in the IRPhEP Handbook with the identifier TREAT-FUND-RESR-001. Detailed models of TREAT were developed and evaluated to prepare benchmark models of the minimum critical mass configuration. Calculations were performed using MCNP6.1 with ENDF/B-VII.1 nuclear data libraries. Discussion of results.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Data for: Spatial access and resource limitations control carbon mineralization in soils

This dataset contains data and code used for the paper "Spatial access and resource limitations control carbon mineralization in soils", https://doi.org/10.1016/j.soilbio.2021.108427. Core-scale soil carbon fluxes are ultimately regulated by pore-scale dynamics of substrate availability and microbial access. These are constrained by physicochemical and biochemical phenomena (e.g. spatial access and hydrologic connectivity, physical occlusion, adsorption-desorption with mineral surfaces, nutrient and resource limitations). We conducted an experiment to determine how spatial access and resource limitations influence core-scale water-soluble soil organic matter (SOM) mineralization, and how these are regulated by antecedent moisture conditions. Intact soil cores were incubated at field-moist vs. drought conditions, after which they were saturated from above (to simulate precipitation) or below (to simulate groundwater recharge). Soluble carbon (acetate) and nitrogen (nitrate) forms were added to some cores during the rewetting process to alleviate potential nutrient limitations. Soil respiration was measured during the incubation, after which pore water was extracted from the saturated soils and analyzed for water soluble organic carbon concentrations and characterization. Our results showed that carbon (C) amendments increased the cumulative carbon dioxide (CO2) evolved from the soil cores, suggesting that the soils were C-limited. Drought and rewetting increased soil respiration, and there was a greater abundance of complex aromatic molecules in pore waters sampled from these soils. This newly available substrate appeared to alleviate nutrient limitations on respiration, because there were no further respiration increases with subsequent C and N amendments. We had hypothesized that respiration would be influenced by wetting direction, as simulated precipitation would mobilize C from the surface. However, as a main effect, this response was seen only in the C-amended soils, indicating that surface-C may not have been bioavailable. At the pore scale (pore water samples), drought and the C, N amendments caused a net loss of identified molecules when the soils were rewet from below, whereas wetting from above caused a net increase in identified molecules, suggesting that fresh inputs stimulated the C-and N-limited microbial populations present deeper in the soil profile. Our experiment highlights the complex and interactive role of antecedent moisture conditions, wetting direction, and resource limitations in driving core-scale C fluxes.This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/kaizadp/TES_spatial_access_2021 (DOI: 10.5281/zenodo.5522938).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

COUPLED MULTIPHYSICS PRIMARY LOOP SIMULATIONS OF THE Mk1-FHR IN THE VIRTUAL TEST BED

The Virtual Test Bed (VTB) repository hosts a wide range of challenge problems to showcase modeling and simulation capabilities to support advanced reactor demonstration. A coupled multiphysics model of the Mk1 pebble-bed fluoride-salt-cooled high-temperature reactor (PB-FHR) is presented here. The analysis leverages NEAMS tools (Griffin), SAM, Pronghorn, and the MOOSE heat conduction module) for core neutronics, thermal hydraulics of the core and the primary loop, and multiscale fuel performance simulations. The analysis was entirely created by coupling standalone simulations of the reactor previously available on the VTB. All input files and documentation developed for this example are available on the VTB website: \url{mooseframework.inl.gov/virtual_test_bed/}}. This model was featured in the NRIC Tech Talk presented in December 2021.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Radiative divertor detachment and impurity transport with nitrogen and neon seeding in KSTAR H-mode plasmas

Achieving core-edge compatible divertor detachment is a critical requirement for stable operation in future fusion devices. This study compares nitrogen and neon seeding in KSTAR H-mode plasmas with carbon walls, combining experiments and SOLPS-ITER modelling to evaluate their radiative dissipation and core-edge compatibility. Experimentally, N seeding achieved stronger divertor detachment, with larger reductions in target particle and heat fluxes, a higher divertor radiation fraction, and a lower core radiation fraction than Ne. In contrast, Ne seeding triggered a significant rise in core radiation followed by H–L back transitions, limiting the maximum total radiated power fraction to roughly half that of N. SOLPS-ITER simulations reproduced the experimental trends and revealed that the better core-edge compatibility of N arises from its higher divertor retention in addition to its higher cooling factor. The relative positions of the stagnation points of impurity poloidal velocity and ionization sources did not explain the different divertor compression. Instead, in the present modelling, the higher impurity parallel particle flux, resulting from the higher parallel impurity velocity, explains the stronger nitrogen impurity compression in the divertor region. The parallel temperature distribution with N was more favourable for achieving higher impurity parallel velocity than with Ne, because the impurity velocity is governed by modifications of the main ion flow due to friction and thermal forces, both of which strongly depend on the temperature. Ultimately, this behaviour is attributed to the strongly divertor-localized radiation of N. These results demonstrate that N is more effective than Ne in achieving radiative divertor detachment while maintaining low core contamination in KSTAR, consistent with observations in other present tokamaks.

KSTAR↗
DOE OSTI2013

Evaluation of Preindustrial to Present-day Black Carbon and its Albedo Forcing from Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP)

As a part of the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), we evaluate the historical black carbon (BC) aerosols simulated by 8 ACCMIP models against the observations including 12 ice core records, a long-term surface mass concentrations and recent Arctic BC snowpack measurements. We also estimate BC albedo forcing by performing additional simulations using the NCAR Community Land and Sea-Ice model 4 with prescribed meteorology from 1996-2000, which includes the SNICAR BC-snow model. We evaluated the vertical profile of BC snow concentrations from these offline simulations to using recent BC snowpack measurements. Despite using the same BC emissions, global BC burden differs by approximately a factor of 3 among models due to the differences in aerosol removal parameterizations and simulated meteorology among models; 34 Gg to 103 Gg in 1850 and 82 Gg to 315 Gg in 2000. However,models agree well on 2.5~3 times increase in the global BC burden from preindustrial to present-day, which matches with the 2.5 times increase in BC emissions. We find a large model diversity at both NH and SH high latitude regions for BC burden and at SH high latitude regions for deposition fluxes. The ACCMIP simulations match the observed BC mass concentrations well in Europe and North America except at Jungfrauch and Ispra. However, the models fail to capture the Arctic BC seasonality due tosevere underestimations during winter and spring. Compared to recent snowpack measurements, the simulated vertically resolved BC snow concentrations are, on average, within a factor of 2-3 of observations except for Greenland and Arctic Ocean. However, model and observation differ widely due to missing interannual variations in emissions and possibly due to the choice of the prescribed meteorology period (i.e., 1996-2000).

ACCMIP↗
DOE OSTI2023

On dysprosium utilisation in multi-main-phase Nd–Dy–Fe–B magnets with core–shell microstructures

The development of high-performance Nd–Dy–Fe–B magnets that minimise the consumption of the scarce rare earth (RE) element Dy remains a major global scientific and technological quest. Here, we designed an alloy microstructure comprising of a uniform Dy-lean core–Dy-rich shell in a series of multi-main-phase (MMP) Nd–Dy–Fe–B magnets. The resulting MMP Dy1 and Dy3 magnets with an overall Dy level of 1 and 3 wt.% possessed values of 0.48 and 0.29 T/wt.% of coercivity increment per unit weight percentage of the Dy addition, respectively. Most importantly, the resulting MMP Dy3 magnet exhibited a high coercivity (2.38 T), an excellent thermal stability of the coercivity (|β| = 0.531%/°C), a high squareness factor (> 95%), all with little diminishment in the remanent magnetisation (1.35 T) and maximum energy product (43.6 MGOe). These properties are superior to the currently available sintered Nd–Dy–Fe–B magnets which utilise higher levels of Dy of 5 wt.%. Via magnetic and multi-scale microstructural characterisation experiments and micromagnetic simulations, the formation of the Dy-lean core–Dy-rich shell microstructure is rationalised via solid-state-diffusion and solution reprecipitation during liquid-phase sintering. The Dy-lean core–Dy-rich shell microstructure and the non-ferromagnetic low-Fe RE-rich grain boundary phase led to the synergistic magnetic performance. This is significant in the context of the MMP Nd–Dy–Fe–B magnets being applied to large-scale production. The present work establishes a pathway for the more sustainable utilisation of Dy in permanent magnets via formation of a uniform core–shell microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Rotationally induced fragmentation in the prestellar core L1544

Recent observations indicate that there is no correlation between the level of turbulence and fragmentation in detected protostellar cores, suggesting that turbulence works mainly before gravitationally bound prestellar cores form and that their inner parts are likely to be velocity coherent. Based on this evidence, we simulate the collapse and fragmentation of an isolated, initially centrally condensed, uniformly rotating core of total mass M = 5.4 M {sub ☉}, using the smoothed particle hydrodynamics code GADGET-2 modified with the inclusion of sink particles, in order to compare the statistical properties of the resulting stellar ensembles with previous gravoturbulent fragmentation models. The initial conditions are intended to fit the observed properties of the prestellar core L1544. We find that for ratios of the rotational to the gravitational energy β ≥ 0.05, a massive disk is formed at the core center from which a central primary condenses after ∼50 kyr. Soon thereafter the disk fragments into secondary protostars, consistent with an intermediate mode of star formation in which groups of 10-100 stars form from a single core. The models predict peak accretion rates between ∼10{sup –5} and 10{sup –4} M {sub ☉} yr{sup –1} for all stars and reproduce many of the statistical properties predicted from gravoturbulent fragmentation, suggesting that on the small scales of low-mass, dense cores these are independent of whether the contracting gas is turbulent or purely rotating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SU-D-BRF-06: A Brachytherapy Simulator with Realistic Haptic Force Feedback and Real-Time Ultrasounds Image Simulation for Training and Teaching

Purpose: Surgical procedures require dexterity, expertise and repetition to reach optimal patient outcomes. However, efficient training opportunities are usually limited. This work presents a simulator system with realistic haptic force-feedback and full, real-time ultrasounds image simulation. Methods: The simulator is composed of a custom-made Linear-DELTA force-feedback robotic platform. The needle tip is mounted on a force gauge at the end effector of the robot, which responds to needle insertion by providing reaction forces. 3D geometry of the tissue is using a tetrahedral finite element mesh (FEM) mimicking tissue properties. As the needle is inserted/retracted, tissue deformation is computed using a mass-tensor nonlinear visco-elastic FEM. The real-time deformation is fed to the L-DELTA to take into account the force imparted to the needle, providing feedback to the end-user when crossing tissue boundaries or needle bending. Real-time 2D US image is also generated synchronously showing anatomy, needle insertion and tissue deformation. The simulator is running on an Intel I7 6- core CPU at 3.26 MHz. 3D tissue rendering and ultrasound display are performed on a Windows 7 computer; the FEM computation and L-DELTA control are executed on a similar PC using the Neutrino real-time OS. Both machines communicate through an Ethernet link. Results: The system runs at 500 Hz for a 8333-tetrahedron tissue mesh and a 100-node angular spring needle model. This frame rate ensures a relatively smooth displacement of the needle when pushed or retracted (±20 N in all directions at speeds of up to 2 m/s). Unlike commercially-available haptic platforms, the oblong workspace of the L-DELTA robot complies with that required for brachytherapy needle displacements of 0.1m by 0.1m by 0.25m. Conclusion: We have demonstrated a real-life, realistic brachytherapy simulator developed for prostate implants (LDR/HDR). The platform could be adapted to other sites or training for other types of needle-based procedures.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

CFD Analysis of Upper Plenum Flow for a Sodium-Cooled Small Modular Reactor

Upper plenum flow behavior is important for many operational and safety issues in sodium fast reactors. The Prototype Gen-IV Sodium Fast Reactor (PGSFR), a pool-type, 150 MWe output power design, was used as a reference case for a detailed characterization of upper plenum flow for normal operating conditions. Computational Fluid Dynamics (CFD) simulation was utilized with detailed geometric modeling of major structures. Core outlet conditions based on prior system-level calculations were mapped to approximate the outlet temperatures and flow rates for each core assembly. Core outlet flow was found to largely bypass the Upper Internal Structures (UIS). Flow curves over the shield and circulates within the pool before exiting the plenum. Cross-flows and temperatures were evaluated near the core outlet, leading to a proposed height for the core outlet thermocouples to ensure accurate assembly-specific temperature readings. A passive scalar was used to evaluate fluid residence time from core outlet to IHX inlet, which can be used to assess the applicability of various methods for monitoring fuel failure. Additionally, the gas entrainment likelihood was assessed based on the CFD simulation results. Based on the evaluation of velocity gradients and turbulent kinetic energies and the available gas entrainment criteria in the literature, it was concluded that significant gas entrainment is unlikely for the current PGSFR design.

computational fluid dynamics↗
DOE OSTI2014

DARK MATTER HEATING AND EARLY CORE FORMATION IN DWARF GALAXIES

We present more results from a fully cosmological ΛCDM simulation of a group of isolated dwarf galaxies that has been shown to reproduce the observed stellar mass and cold gas content, resolved star formation histories, and metallicities of dwarfs in the Local Volume. Here we investigate the energetics and timetable of the cusp-core transformation. As suggested by previous work, supernova-driven gas outflows remove dark matter (DM) cusps and create kiloparsec-size cores in all systems having a stellar mass M {sub *} > 10{sup 6} M {sub ☉}. The {sup D}M core mass removal efficiency{sup —}dark mass ejected per unit stellar mass—ranges today from a few to a dozen, and increases with decreasing host mass. Because dwarfs form the bulk of their stars prior to redshift 1 and the amount of work required for DM heating and core formation scales approximately as M{sub vir}{sup 5/3}, the unbinding of the DM cusp starts early and the formation of cored profiles is not as energetically onerous as previously claimed. DM particles in the cusp typically migrate to 2-3 core radii after absorbing a few percent of the energy released by supernovae. The present-day slopes of the inner DM mass profiles, Γ ≡ dlog M/dlog R ≅ 2.5-3, of the simulated ''Bashful'' and ''Doc'' dwarfs are similar to those measured in the luminous Fornax and Sculptor dwarf spheroidals. None of the simulated galaxies has a circular velocity profile exceeding 20 km s{sup –1} in the inner 1 kpc, implying that supernova feedback is key to solve the ''too-big-to-fail'' problem for Milky Way subhalos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A multi-group Monte Carlo core analysis method and its application in SCWR design

Complex geometry and spectrum have been the characteristics of many newly developed nuclear energy systems, so the suitability and precision of the traditional deterministic codes are doubtable while being applied to simulate these systems. On the contrary, the Monte Carlo method has the inherent advantages of dealing with complex geometry and spectrum. The main disadvantage of Monte Carlo method is that it takes long time to get reliable results, so the efficiency is too low for the ordinary core designs. A new Monte Carlo core analysis scheme is developed, aimed to increase the calculation efficiency. It is finished in two steps: Firstly, the assembly level simulation is performed by continuous energy Monte Carlo method, which is suitable for any geometry and spectrum configuration, and the assembly multi-group constants are tallied at the same time; Secondly, the core level calculation is performed by multi-group Monte Carlo method, using the assembly group constants generated in the first step. Compared with the heterogeneous Monte Carlo calculations of the whole core, this two-step scheme is more efficient, and the precision is acceptable for the preliminary analysis of novel nuclear systems. Using this core analysis scheme, a SCWR core was designed based on a new SCWR assembly design. The core output is about 1,100 MWe, and a cycle length of about 550 EFPDs can be achieved with 3-batch refueling pattern. The average and maximum discharge burn-up are about 53.5 and 60.9 MWD/kgU respectively. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Computational Benefit of GPU Optimization for the Atmospheric Chemistry Modeling

Abstract Global chemistry‐climate models are computationally burdened as the chemical mechanisms become more complex and realistic. Optimization for graphics processing units (GPU) may make longer global simulation with regional detail possible, but limited study has been done to explore the potential benefit for the atmospheric chemistry modeling. Hence, in this study, the second‐order Rosenbrock solver of the chemistry module of CAM4‐Chem is ported to the GPU to gauge potential speed‐up. We find that on the CPU, the fastest performance is achieved using the Intel compiler with a block interleaved memory layout. Different combinations of compiler and memory layout lead to ~11.02× difference in the computational time. In contrast, the GPU version performs the best when using a combination of fully interleaved memory layout with block size equal to the warp size, CUDA streams for independent kernels, and constant memory. Moreover, the most efficient data transfer between CPU and GPU is gained by allocating the memory contiguously during the data initialization on the GPU. Compared to one CPU core, the speed‐up of using one GPU alone reaches a factor of ~11.7× for the computation alone and ~3.82× when the data transfer between CPU and GPU is considered. Using one GPU alone is also generally faster than the multithreaded implementation for 16 CPU cores in a compute node and the single‐source solution (OpenACC). The best performance is achieved by the implementation of the hybrid CPU/GPU version, but rescheduling the workload among the CPU cores is required before the practical CAM4‐Chem simulation.

54 ENVIRONMENTAL SCIENCES↗