Search NASA⌕ Search

Engineering topics

Ghoos, Kristel

Publications and source records attributed to Ghoos, Kristel.

Neutronics Evaluation of the Preliminary Design of the Service Cell Door for the Second Target Station

The service cell door needs to provide sufficient shielding to achieve 0.25 mrem/h dose rates in the adjacent room while remote handling operations of waste materials are going on within the service cell for the Second Target Station (STS). This analysis provides a neutronics evaluation of the preliminary design of the service cell door to evaluate the required thickness and the effect of streaming gaps around the door. As radiation source, we use the decay gamma source of the moderator reflector assembly (MRA) of the STS with a decay time of either two days or six months. As the MRA contains a significant amount of beryllium, photonuclear reactions induce a neutron source, which dictates the required door thickness in certain situations. Four scenarios with different MRA positions and decay times have been analyzed. The required thickness differs drastically depending on the chosen scenario and on the material of the door. The streaming through the gaps around the door is relatively small with high dose rates limited to small sections of the geometry.

43 PARTICLE ACCELERATORS↗

Moderator Optimization for the Second Target Station Final Design

This report details the results from the optimization simulations performed for the final design of the Second Target Station (STS). This is a continuation of the analysis performed in 2022 for the preliminary design. To evaluate the impact of the design changes since 2022, the analysis was repeated with updated target and moderator models. Additionally, more degrees of freedom have been considered in the premoderator dimensions which provide a refinement compared to the previous optimization analysis.

43 PARTICLE ACCELERATORS↗

Overview of Shielding Analyses at Oak Ridge National Laboratory’s Spallation Neutron Source Second Target Station

The Neutronics Group of the Second Target Station project at Oak Ridge National Laboratory is responsible for all neutronics analyses related to design and construction. This paper provides an overview of four neutronics analyses performed for the Second Target Station project, which is representative of all the ongoing work but especially tasks related to shielding. These analyses extend from the proton accelerator through the target monolith and bunker to the end of a neutron beamline. Each example highlights the tools and methods the Neutronics Group uses for analysis. The primary tool is the Monte Carlo radiation transport code MCNP 6.2, but this is augmented by other codes that supplement the input of MCNP. The other codes specifically highlighted in this paper include Attila4MC, ADVANTG, and AARE.

Miller, Thomas↗

Neutronics Analysis of Heating and DPA of the Proton Beam Window and Core Vessel Nozzle Neutronics Heating Analysis for the Second Target Station

This reports details the analysis of the heating rates, dpa values and He production in the Proton Beam Window (PBW) assembly and the Core Vessel (CV) nozzle for the preliminary design for the Second Target Station (STS). This analysis focuses on the heating rates in the cooled beltline component of the CV. Neutronics calculations use a detailed unstructured mesh (UM) converted from a CAD model in SpaceClaim of the PBW and CV design of February 2024. The rest of the geometry is taken from the STS master model in constructive solid geometry (CSG). Heating is evaluated for four beam profiles, which are chosen as the most likely to produce the highest heating in the beltline component. The maximum energy deposition values in the beltline are 0.02-0.08 J/cc/pulse, depending on the beam profile. Detailed point cloud data of the full UM geometry is made available for mechanical analysis. The maximum dpa rate in the beltline steel is approximately 0.02 dpa/year, which corresponds to a lifetime of 3,000 years. These results are in good agreement with previous results for the CV. For the PBW itself, maximal heating, dpa, and He production rates are very similar to the previous analysis.

43 PARTICLE ACCELERATORS↗

Proton Beam Window In Situ Neutronics Analysis for the Preliminary Design for the Second Target Station

These analyses provide estimates for the residual dose rates for various geometric configurations related to proton beam window (PBW) removal and replacement maintenance events. This information will be used to plan a remote handling strategy for the replacement and maintenance of the PBW. The 5 and 100 mrem/h dose rate contour lines are of specific interest as they classify radiation areas (between 5 and 100 mrem/h) and high radiation areas (>100 mrem/h).

43 PARTICLE ACCELERATORS↗

Optimization of the Second Target Station cold source moderators using an automated workflow

The Second Target Station (STS) at the US Department of Energy’s Oak Ridge National Laboratory is designed to become the world’s highest peak-brightness spallation source of cold neutrons. Successful completion of the STS, which is currently in the preliminary design phase, will provide transformative new capabilities to examine novel materials for future technologies. At STS, neutrons will be generated by spallation reactions in a solid tungsten target. They will be moderated and thermalized in two cold (20 K) para-hydrogen moderators. Careful optimization of these moderators is essential to the project’s success. To find optimal moderator designs, an advanced optimization workflow integrates high-fidelity neutronics calculations using the Monte Carlo N-Particle (MCNP) transport code MCNP6.2 with state-of-the-art optimization algorithms in the Dakota optimization toolkit. For each design iteration, a parametrized solid CAD geometry is generated in Creo and automatically converted into an unstructured mesh geometry by Attila 4MC for the neutronics calculation with MCNP. Iterations repeat until optimal designs are found. Herein this paper presents the results of a sensitivity and optimization study for the cylindrical and tube moderators. Both moderators can be optimized for maximum peak brightness, maximum time-integrated brightness, or any combination between these extremes. Maximum peak brightness is achieved by using smaller optimal dimensions of the moderators, whereas maximum time-integrated brightness is achieved by using larger dimensions. A Pareto front details the designs that optimally balance both brightness metrics. The Pareto front can be found in only 40–110 iterations with 4–5 design parameters when using the efficient global and Pareto-set optimization algorithms in Dakota. Additionally, important engineering constraints can be taken into account, such as the coupling between the cylindrical moderator radius and aluminum vessel wall thicknesses required to ensure structural integrity of the vessels. This interaction has a significant impact on the resulting optimal designs. Our new, highly efficient, fully automated optimization workflow will be used to optimize additional STS components in the future and can be adopted for design and optimization studies at other experimental neutron and accelerator facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Shielding Thickness Recommendations for an Activated Target Segment Vertical Stave

This analysis examines the shielding requirements needed to plan maintenance and handling operations of the target segment assembly vertical stave. It provides the analyses and results that fulfill the requirements of Task Order TO-150. The vertical stave is a stainless-steel welded assembly that carries cooling water to and from the target block. In this report, we characterize the dose rates and shielding requirements for the vertical stave only, i.e., separated from the target block. Shielding and dose rates for the target block only have been calculated in response to TO-008.

61 RADIATION PROTECTION AND DOSIMETRY↗

Coupled Target-Beam-Moderator Optimization for the Second Target Station

This report describes the results for a coupled target-beam-moderator optimization analysis for the Second Target Station (STS) at ORNL's Spallation Neutron Source. This study is a continuation of the optimization analysis for the moderators in the preliminary design of STS performed in 2022. In the 2022 analysis the dimensions of the moderators are parameterized, while the target and the proton beam profile are kept constant. In this analysis the target height and the proton beam profile are added as parameters. This allows to study the coupled effects of changing target, moderator and beam dimensions. Similar to the 2022 analysis, this work is performed with an automated optimization workflow that uses the optimization toolbox DAKOTA, parameterized geometries in CREO and SpaceClaim, the unstructured mesh generation in Attila4MC, and the particle transport code MCNP6.2©. This workflow enables an efficient optimization using high-fidelity geometries. The main conclusions of this analysis are the following: • Coupled beam-target-moderator optimization provides a few additional percent performance gain over stand-alone moderator optimization. • The moderator performance is not very sensitive to the target height (between ≈60 and ≈80 mm) as long as the beam profile is chosen adequately. • The moderator performance is sensitive to the choice of beam spatial standard deviations, even when the footprint is kept constant. • The optimal moderator radius is the same for a beam footprint of 30 cm 2 , 62.5 cm 2 , and 90 cm 2 . Also the slope of the super-gaussian beam profile does not significantly impact the optimal moderator radius. • The optimal parameters and sensitivities are very similar to the 2022 optimization analysis. These results only indicate a a difference in the optimal radius of the cylindrical moderator, however, this has been corrected in the final design moderator optimization. The main purpose of this report is to document the simulations, results and lessons learned. The most impactful results are summarized in. We also note that the target geometry used in this work is not the final design.

43 PARTICLE ACCELERATORS↗

Activation Analysis in Preparation for a Tungsten Irradiation Experiment at LANSCE

To organize the safe handling of activated material, knowing the residual dose rates is crucial. Here, in this work, we present the pre-experiment activation analysis for an experiment in which tungsten blocks are irradiated by 800-MeV protons. In this analysis, we use the Monte Carlo N-Particle (MCNP) code for radiation transport, Attila4MC for unstructured mesh generation, and Activation in Accelerator Radiation Environments (AARE), including CINDER2008, for activation analysis. If the tungsten blocks must be removed within a day after the experiment, then exposure to personnel entering the room must be reduced. One exposure-reduction strategy is to add carbon steel shielding around the tungsten blocks, efficiently reducing the dose from the activated tungsten. However, the shielding becomes activated itself during irradiation: 56Mn is the dominant contributor for short decay times. The actual schedule at the time of the experiment allowed sufficient cool-off time for the tungsten in the room so that additional shielding was not necessary. A less rigorous comparison of the calculated values with the post-experiment measurements showed reasonable agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Second Target Station Project: LANSCE WNR Target 2 (Blue Room) Experiment 2022

The Second Target Station (STS) at SNS will address emerging scientific challenges by providing a source of intense, cold neutrons to instruments optimized for this source, using rotating, tungsten target blocks. The STS target will receive 1.3 GeV proton beam pulses from the SNS accelerator at a repetition rate of 15 Hz. The facility life is planned for 40 years, and each target assembly life is expected to be approximately 10 years. An accurate strain prediction is then critical for fatigue life assessment of STS target blocks because they will be subject to approximately 108 beam pulses per lifetime. As an R&D activity, the LANSCE WNR Target 2 (blue room) facility was used to test the strain response of prototypical target blocks to the thermal shock of a proton pulse. The blue room was well suited for a pulsed proton beam impact test of subscale STS target blocks; the 800 MeV proton energy is approximately 60% of the 1.3 GeV proton energy expected from the SNS accelerator to the STS. Both the LANSCE Proton Storage Ring (PSR) and SNS are short-pulse proton beam sources with nominal pulse widths of 250 ns and 661 ns, respectively, so the energy deposition in the target occurs in < 1 microsecond pulse duration.

43 PARTICLE ACCELERATORS↗

An unstructured mesh based neutronics optimization workflow

We have developed a fully automated workflow to optimize the neutronics performance of the Second Target Station (STS) at the Oak Ridge National Laboratory’s Spallation Neutron Source. The optimization workflow starts with the parametrized solid CAD engineering models and converts them into the unstructured mesh (UM) models for the neutronics calculations with MCNP6.2. Calculations are executed and their results are loaded into the Dakota optimization toolkit. Dakota analyzes the results and proposes new geometry parameters for the next design iteration. The cycle repeats until the optimal parameters are found. The automated CAD to MCNP conversion, the use of high-fidelity UM models, and the use of modern optimizer are the key elements that advance the entire optimization workflow in comparison with the original workflow. The original workflow was based on a simplified constructive solid geometry (CSG) modeling with MCNPX, mcnp_pstudy tool, and an in-house optimizer. Herein to demonstrate the new workflow, we present a case of neutronics optimization of the moderator–reflector assembly (MRA). Apart from the MRA, the workflow can optimize other major STS components, such as the spallation target, neutron beamlines, radiation shielding, and various accelerator components. Importantly, the new workflow opens the door to the advanced multi-physics multi-parameter optimization and has the potential for use in other nuclear physics and accelerator applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Numerical error analysis of SOLPS-ITER simulations of EAST

Abstract Plasma edge simulations with codes like SOLPS-ITER are widely employed to interpret fusion experiments. However, numerical errors appearing in such simulations are rarely investigated, despite their potential large impact on simulation results. These errors consist of the statistical error and the bias, both resulting from the finite number of employed EIRENE Monte Carlo particles and incomplete convergence, and the discretization error due to the finite resolution of the computational grids. In this contribution, the resulting numerical errors on simulations of pure deuterium and neon seeded H-mode EAST discharges are examined. The statistical error can be kept small compared to other numerical error contributions by averaging the plasma profiles. This allows investigating the bias and discretization errors using Richardson extrapolation. It is shown that grid refinement and the number of employed Monte Carlo particles have the largest influence on the result, in agreement with similar studies of an ITER deuterium case. For the first time, numerical error bars on the entire simulated target profiles are determined showing that the largest numerical error is 17.9%, mainly due to the plasma grid discretization. On top, also numerical errors on simulated neutral pressures are investigated in detail, for which the statistical error is dominant. The analysis demonstrates which setup is needed to keep numerical errors limited: the SOLPS-ITER averaging procedure should be employed including enough EIRENE particles, and the involved grids should be sufficiently refined to reduce discretization errors.

Boeyaert, Dieter (ORCID:0000000309208660)↗

Neutronics Analysis in Preparation of the Blue Room Experiment

The goal of the Blue Room experiment is to investigate the ‘thermal shock’ stress in a tungsten target provoked by proton pulse impact. Short, strong pulses of highly energetic protons cause a rapid energy deposition in tungsten, which leads to a sudden temperature rise. This thermal shock causes a stress response and leads to fatigue, which can limit the lifetime of the target. In the Blue Room experiment, tungsten blocks will be irradiated by pulses of 800 MeV protons with several intensities to gain understanding in the strain response of the tungsten block.

36 MATERIALS SCIENCE↗

Heating, DPA, and He production in the Proton Beam Window Assembly for the Second Target Station

This report details the analysis of the heating rates, dpa values, and He production in the Proton Beam Window (PBW), the assembly around the PBW, and the shielding above it. Neutronics calculations use a detailed Unstructured Mesh (UM) converted from the solid CAD model of the PBW design of June 2022. The rest of the geometry is taken from the CSG master model (as of June 2022). For this analysis, a proton beam profile with a footprint of ≈90 cm 2 has been used. This proton profile was created by a Ring to Second Target magnet array configuration using octupole magnets. Maximum values in the PBW (aluminum) are 5.4 J/cc/pulse heating rate, 1.1 dpa per year, and 660 appm He-production per year. These results within ~ 10% with the previous calculations from June 2020. Average heating and dpa rates have been calculated for each component in the PBW assembly and will be used to evaluate and design the cooling needs and make decisions on removable components for the PBW. It is desirable to optimize and simplify the PBW design by only replacing components that are at the end of their lifetime. Additional analysis for the energy deposition in the PBW is performed for two super-Gaussian beam profiles with a footprint of 60 cm 2 and 30 cm 2 respectively. Also, an estimate of the PBW lifetime is performed for two realistic quadrupole beam profiles.

43 PARTICLE ACCELERATORS↗

Optimization of the moderators in the STS preliminary design

This report details the results for an optimization of the dimensions of the moderators in the preliminary design of the Spallation Neutron Source Second Target Station (STS). This study uses the optimization algorithms of Dakota and an unstructured mesh model for the moderators in MCNP. More details on the unstructured mesh model and the automated mesh generation can be found in [3]. Parallel to this effort, the same moderator geometries have been optimized using a constructive solid geometry (CSG) MCNP model. More details on this model and its results can be found in [4]. Three optimal designs are selected for each moderator: one that is optimized for maximum peak brightness, one for maximum time-integrated brightness, and one for a combination of peak and time-integrated brightness. The backbone of the optimization work flow is provided by Dakota. For each set of design parameters requested by Dakota, a new solid geometry is automatically built in Creo and SpaceClaim, and subsequently exported to Attila4MC to generate an unstructured mesh geometry for MCNP. After the MCNP calculation is finished, the objective function (e.g., brightness metric) is returned to Dakota. After the new design has been evaluated, a result-file is written, and Dakota proposes the next set of design parameters to be evaluated. The loop continues until a specified convergence criterion has been met. The design parameters of the cylindrical (upper) moderator include the hydrogen radius, the premoderator thickness (top, bottom, radial), the beryllium radius and the horizontal position of the moderator. The crucial design choice is the hydrogen radius. A radius of 62 mm is shown to provide the maximum time-integrated brightness. The maximum peak brightness occurs with a radius of 40 mm. A combined (middle) design, which balances peak and time-integrated brightnesses, is obtained with a hydrogen radius of 50 mm. The premoderator thicknesses and the beryllium radius are slightly larger in the design optimized for time-integrated brightness than in the design optimized for peak brightness. The sensitivity to these two parameters is relatively small close to the optimal configurations. The hydrogen vessel and vacuum vessel wall thicknesses are dependent on the radius of the liquid hydrogen due to structural integrity requirements. The increased wall thicknesses for larger vessels significantly penalize the time-integrated brightness, with the maximum obtainable value reduced by more than 10% relative to earlier studies which used fixed vessel wall thicknesses. The impact of the variable wall thicknesses is much less for the peak brightness and combined brightness designs. The design parameters of the tube (lower) moderator selected for the optimization are the tube length, the annular premoderator thickness, the beryllium radius and the horizontal position of the moderator. The tube length is the crucial parameter and is chosen large (210 mm) and small (125 mm) in the designs optimized for time-integrated and peak brightness respectively. A combined optimal design has a tube length of 170 mm. The premoderator thickness and the beryllium radius are chosen larger in the design optimized for time-integrated brightness.

42 ENGINEERING↗