Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reach code”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Characterization of cartridge filters from the IEA-R1 Nuclear Reactor

The management of radioactive waste ensures safety to human health and the environment nowadays and for the future, without overwhelming the upcoming generations. The primary characterization of radioactive waste is one of the main steps in the management of radioactive waste. This step permits to choose the best treatment for the radioactive waste before forwarding it to its final disposal. The aim of the present work is the primary characterization of cartridge filters from the IEA-R1 nuclear reactor utilizing gamma-ray spectrometry, and the method of Monte Carlo for calibration. The IEA-R1 is located in the Nuclear and Energy Research Institute (IPEN - CNEN) in the city of Sao Paulo, Brazil. Cartridge filters are used for purification of the cooling water that is pumped through the core of the pool type nuclear research reactors. Once worn out, these filters are replaced and then become radioactive waste. Determination of the radioactive inventory is of paramount importance in the management of such radioactive waste, and one of the main methods for doing so is the gamma-ray spectrometry, which can identify and quantify high energy photon emitters. The technique chosen for the characterization of radioactive waste in the present work is the gamma-ray spectrometry with High purity Germanium (HPGe) detectors. From the energy identified in the experimental spectrum, three radioisotopes were identified in the cartridge filter: {sup 108m}Ag, {sup 110m}Ag, {sup 60}Co. For the estimated activity of the filter, the calibration in efficiency was made utilizing the MCNP4C code of the Monte Carlo method. Such method was chosen because there is no standard source available in the same geometry of the cartridge filter, therefore a simulation had to be developed in order to reach a calibration equation, necessary to estimate the activity of the radioactive waste. The results presented an activity value in the order of MBq for all radioisotopes. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

High Accuracy Transition Metal Effective Cores for the Many-Body Diffusion Monte Carlo Method

Practical applications of the real-space diffusion Monte Carlo (DMC) method require the removal of core electrons, where currently localization approximations of semilocal potentials are generally used in the projector. Accurate calculations of complex solids and large molecules demand minimizing the impact of approximated atomic cores. Prior works have shown that the errors from such approximations can be sizable in both finite and periodic systems. In this work, we show that a class of differential pseudopotentials, known as pseudo-Hamiltonians, can be constructed for the 3d transition metal atoms, entirely removing the need for any localization scheme in the DMC projector. As a proof of principle, we demonstrate the approach for the case of Co. In order to minimize errors in the pseudo-Hamiltonian at the many-body level, we generalize the recently proposed correlation-consistent pseudopotential generation scheme to successively close semilocal representations of the differential potentials. Our generation scheme successfully produces potentials tailored specifically for real space projector quantum Monte Carlo methods with low error at the many-body level, i.e., with many-body scattering properties very close to relativistic all-electron results. In particular, we show that the agreement with respect to atomic and molecular quantities reach chemical accuracy in many cases-on par with the most accurate semilocal pseudopotentials available. Further, our pseudo-Hamiltonian generation scheme utilizes standard quantum chemistry codes designed only to work with semilocal pseudopotentials, enabling straightforward generation of pseudo-Hamiltonians for additional elements in future works.

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

Summary of 2020 edge-plasma modeling for FNSFand plans for 2021

This year, following the study of liquid-wall options, our edge-plasma modeling has returned to analyzing issues associated with a solid tungsten-wall, Fusion Nuclear Science Facility (FNSF)-size device, primarily using the 2D UEDGE plasma/neutral transport code. The focus is to predict plasma/neutral heat and particle fluxes to various wall components and determine the concentration of injected and wall-sputtered impurities that reach the core plasma region. Our models have been advanced in a number of areas compared to the initial modeling done in the 2017 timeframe: this year we include each individual charge-state of impurity ions instead of assuming a fixed concentration of impurities relative to the deuterium-tritium (DT) fuel plasma, including the interaction between inner and outer divertor legs. Here edge fueling of the impurity and pumping in the private-flux region is incorporated. Increasingly detailed atomic-physics rate files for impurity radiation have been compared, as well as the role of molecules. The impact of spatial variation of the turbulence-driven radial plasma transport is considered to contribute to the understanding of the heat-flux-limit effects of choosing a double-null or single-null divertor configuration. Detailed plasma/neutral wall fluxes, densities, and temperatures have been provided to the ERO group for their more detailed studies of wall sputtering and near-surface impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Using Hydrodynamic Similarity as a Verification Method for Impact Cratering Simulations in the FLAG Hydrocode

Hydrodynamic codes (hydrocodes) are common tools for modeling hypervelocity impacts to provide insight into the physical phenomenon. Hydrocodes can simulate impacts from micrometer to kilometer spatial scales and reach impact velocities difficult to achieve in experimental settings. However, numerical models are approximations, and demonstrating that a numerical method is capable of providing physical results for these models is essential. In this work, we employ a hydrocode verification technique that leverages hydrodynamic similarity, a mathematical property of the conservation equations of fluid mechanics that form the basis for hydrocode models. Using the FLAG hydrocode, we simulate aluminum (Al) and basalt projectiles and targets at spatial scales spanning 7 orders of magnitude (hundreds of micrometers to kilometers). These materials were chosen because Al-6061 is a common material in spacecraft and satellites and basalt is a useful approximation of rocky astronomical bodies. Our results show that hydrodynamic similarity holds for each material model used and across spatial scales. We show that under certain conditions hydrodynamic similarity can apply in the presence of gravity and that similarity does not hold in the presence of strength models. We conclude that the FLAG hydrocode preserves important mathematical properties of fluid dynamics in hypervelocity impacts of Al-6061 and basalt.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Multi-cycle reload analysis of a long cycle gas-cooled fast modular reactor

There is currently significant interest in deploying HALEU-fueled fast reactors, including the General Atomics (GA) Fast Modular Reactor (FMR). Such reactors can achieve very long fuel cycles, but with multi-batch loading will take decades to reach equilibrium. This motivates design and analysis of both the initial core and multi-cycle reload, which is typically performed using fast-running, deterministic fast reactor codes such as the Argonne Reactor Computation (ARC) codes. In this paper, multicycle reload of the GA FMR is analyzed using the ARC codes. The GA FMR utilizes 19.75 % enriched fuel in a 16 year cycle with a three-batch strategy, with twice-burned fuel placed on the core periphery. The GA FMR has a softened neutron spectrum due to reflecting elements in the core, so the neutronic solution is first benchmarked against the OpenMC Monte Carlo code. Discrepancy on k eff is 400–600 pcm, likely due to the softened neutron spectrum, heterogeneous fuel assembly design and central reflector. However, the rms discrepancy on the assembly power distribution is only 0.6 %, despite the presence of the central reflector. A reload strategy is devised for the first three cycles of such a reactor, ultimately spanning the first 45–48 years of its operation. The fresh core uses 19.75 %, 19.25 % and 16.75 % enriched fuel in place of fresh, once-burned and twice-burned and is then subsequently refueled with only 19.75 % enriched fuel. The cycle length is varied over 3 cycles of operation to balance fuel utilization and reactor availability, specifically with use of an extended 18-year Cycle 1, followed by a shortened 11-year Cycle 2. Cycle 3 is close to the target 16-year length. Finally, placing twice burned assemblies next to the GA FMR central reflector can reduce power peaking by 3 %, at the expense of slightly reducing the cycle length.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Inhomogeneous energy injection in the 21-cm power spectrum: Sensitivity to dark matter decay

The 21-cm signal provides a novel avenue to measure the thermal state of the Universe during cosmic dawn and reionization (redshifts z ~ 5–30), and thus to probe energy injection from decaying or annihilating dark matter (DM). These DM processes are inherently inhomogeneous: both decay and annihilation are density-dependent, and furthermore, the fraction of injected energy that is deposited at each point depends on the gas ionization and density, leading to further anisotropies in absorption and propagation. In this work, we develop a new framework for modeling the impact of spatially inhomogeneous energy injection and deposition during cosmic dawn, accounting for ionization and baryon density dependence, as well as the attenuation of propagating photons. We showcase how this first completely inhomogeneous treatment affects the predicted 21-cm power spectrum in the presence of exotic sources of energy injection, and forecast the constraints that upcoming HERA measurements of the 21-cm power spectrum will set on DM decays to photons and to electron/positron pairs. These projected constraints considerably surpass those derived from CMB and Lyman-a measurements, and for decays to electron/positron pairs they exceed all existing constraints in the sub-GeV mass range, reaching lifetimes of ~10 28 s. Our analysis demonstrates the unprecedented sensitivity of 21-cm cosmology to exotic sources of energy injection during the cosmic dark ages. Our code, dm21cm, includes all these effects and is publicly available in an accompanying release.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Reduced fast electron transport in shock-heated plasma in multilayer targets due to self-generated magnetic fields

Fast electron transport has been studied in cold solid density CH, cold CH foam (200 mg / cm 3 ), and CH plasma (40 eV 30 mg/cm 3 ) targets—the latter created by shocking the CH foam with a 1.2 kJ long pulse laser and allowing it to expand. The fast electrons were produced using the OMEGA EP laser pulse (800 J, 8 ps) incident on a Au flat target. With the CH plasma, the fluence of fast electrons reaching a Cu foil at the far side of the transport was reduced significantly (25 × weaker peak K $α$ emission). Particle-in-cell simulations using the osiris code modeled fast electron transport in the unshocked foam and plasma cases assuming fixed ionization and including source generation, transport in Au and CH layers, Coulomb collisions, and refluxing. Simulations indicate two main mechanisms which alter electron energy transport through the target between the foam and plasma cases, both due to the magnetic field: a collimating field in the CH region, caused by the resistivity of the return current and more prevalent in the foam; and an insulating field at the Au-CH interface, present only with the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

ADG: automated generation and evaluation of many-body diagrams

The goal of the present paper is twofold. First, a novel expansion many-body method applicable to superfluid open-shell nuclei, the so-called Bogoliubov in-medium similarity renormalization group (BIMSRG) theory, is formulated. This generalization of standard single-reference IMSRG theory for closed-shell systems parallels the recent extensions of coupled cluster, self-consistent Green’s function or many-body perturbation theory. Within the realm of IMSRG theories, BIMSRG provides an interesting alternative to the already existing multi-reference IMSRG (MR-IMSRG) method applicable to open-shell nuclei. The algebraic equations for low-order approximations, i.e., BIMSRG(1) and BIMSRG(2), can be derived manually without much difficulty. However, such a methodology becomes already impractical and error prone for the derivation of the BIMSRG(3) equations, which are eventually needed to reach high accuracy. Based on a diagrammatic formulation of BIMSRG theory, the second objective of the present paper is thus to describe the third version (v3.0) of the code that automatically (1) generates all valid BIMSRG(n) diagrams and (2) evaluates their algebraic expressions in a matter of seconds. This is achieved in such a way that equations can easily be retrieved for both the flow equation and the Magnus expansion formulations of BIMSRG. Expanding on this work, the first future objective is to numerically implement BIMSRG(2) (eventually BIMSRG(3)) equations and perform ab initio calculations of mid-mass open-shell nuclei.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Investigation of Application Specific k eff Sensitivities to Existing and New Subcritical and Critical Experiments

Integral subcritical and critical experiments are essential in providing validation of computational codes which improve confidence in the predicted k eff results of specific applications. This work discusses the relevance of existing integral subcritical and critical experiments to applications far reaching across the nuclear complex. These applications include Travelling Wave Reactors (TWR), Kilopower Reactor Using Stirling Technology (KRUSTY), and plutonium reflected by tantalum. The primary goal of this work is to improve knowledge of energy-dependent cross-section sensitivities for specific applications, identify nuclear data gaps in existing experiments, enable an optimized design of new experiments to address these gaps, and design practical new experiments that could be executed at the National Criticality Experiments Research Center (NCERC) that have a high-value to these applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Extend an innovative HPC-Compatible Multiple Temporal-spatial Resolution Concurrent Finite Element Modeling Approach to Guide Laser Powder Bed Fusion Additive

Laser power bed fusing (PBF) additive manufacturing is a key enabling technology to manufacture highly complex and integrated automotive structures. However, the geometric complexity of PBF-AM technique also leads to highly non-uniform heating and cooling rate in the manufactured part, which may cause flaw formation and produce excessive and nonuniform residual stresses, which increase quality uncertainties and manufacture issues, leading to increases in cost and energy consumption in the form of rejected parts. In this research project, we developed an innovative Multi-Spatial-Temporal-Resolution Finite Element (MUST-FE) method and completed the corresponding high performance computation (HPC) platform-based in-house code, which enables high accuracy prediction of temperature and residual stress fields for component-scale PBF-AM manufacture in efficient computation time. The MUST-FE model is calibrated and validated with a “2D pad” AlSi10Mg experiments by matching the melt pool shape and dimension, and with a “XY-cross” AlSi10Mg experiment by matching the thermal distortion and residual stress. The innovative multi-resolution and concurrent modeling approach adopted in this code ensures accuracy and computational efficiency, which will enable energy-efficient and high-yield, low-cost manufacturing of optimized, qualifiable automotive structures and contribute towards reaching technical targets outlined in AMO’s Program Plan to develop additive manufacturing systems that deliver consistently reliable parts with predictable properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Integration of Ultrawide Bandgap Photoconducting Semiconductor Switches into High-Voltage, High-Current Pulser Systems

In this feasibility study we will evaluate the viability of integrating ultrawide bandgap (UWBG) photoconducting semiconductor switches (PCSS) in high-voltage pulser systems capable of driving low-impedance dynamic loads critical for next-generation accelerator and RF-source applications. The project consisted of developing modelling tools that would inform design of UWBG PCSS packages that could be mounted onto conventional circuit boards. Specifically, we are interested in evaluating the following through simulation software: (1) resistance to electrical breakdown (2) proximity effects due to high-frequency skin effect (3) transient behavior from switching (4) effective transmission of RF signals (5) thermal dynamics of the switch under operation (6) sensitivity to inhomogeneities in optical excitation of the material (7) cross-talk effects from parallel switching structures, and (8) susceptibility of system to damaged components. The project resulted in the development of a special-purpose circuit code that implements photoconductivity physics. This code was used to evaluate objectives 1-6 above, identifying a critical requirement associated with the need for fine control of dopant concentrations in the UWBG material to reach an optimal solution for voltage standoff and laser light responsivity requirements. In addition, we identified that the temperature dependence of the electrical conductivity of this material was typically not well-measured and developed a platform for measuring both the “dark” and “illuminated” conductivity on real physical switches in the lab to better inform our models.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

An algorithm for computing the 2D structure of fast rotating stars

Stars may be understood as self-gravitating masses of a compressible fluid whose radiative cooling is compensated by nuclear reactions or gravitational contraction. The understanding of their time evolution requires the use of detailed models that account for a complex microphysics including that of opacities, equation of state and nuclear reactions. The present stellar models are essentially one-dimensional, namely spherically symmetric. However, the interpretation of recent data like the surface abundances of elements or the distribution of internal rotation have reached the limits of validity of one-dimensional models because of their very simplified representation of large-scale fluid flows. In this article, we describe the ESTER code, which is the first code able to compute in a consistent way a two-dimensional model of a fast rotating star including its large-scale flows. Compared to classical 1D stellar evolution codes, many numerical innovations have been introduced to deal with this complex problem. First, the spectral discretization based on spherical harmonics and Chebyshev polynomials is used to represent the 2D axisymmetric fields. A nonlinear mapping maps the spheroidal star and allows a smooth spectral representation of the fields. The properties of Picard and Newton iterations for solving the nonlinear partial differential equations of the problem are discussed. It turns out that the Picard scheme is efficient on the computation of the simple polytropic stars, but Newton algorithm is unsurpassed when stellar models include complex microphysics. Finally, we discuss the numerical efficiency of our solver of Newton iterations. This linear solver combines the iterative Conjugate Gradient Squared algorithm together with an LU-factorization serving as a preconditioner of the Jacobian matrix.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

High energy conversion efficiency in laser-proton acceleration by controlling laser-energy deposition onto thin foil targets

An all-optical approach to laser-proton acceleration enhancement is investigated using the simplest of target designs to demonstrate application-relevant levels of energy conversion efficiency between laser and protons. Controlled deposition of laser energy, in the form of a double-pulse temporal envelope, is investigated in combination with thin foil targets in which recirculation of laser-accelerated electrons can lead to optimal conditions for coupling laser drive energy into the proton beam. This approach is shown to deliver a substantial enhancement in the coupling of laser energy to 5–30 MeV protons, compared to single pulse irradiation, reaching a record high 15% conversion efficiency with a temporal separation of 1 ps between the two pulses and a 5 μm-thick Au foil. A 1D simulation code is used to support and explain the origin of the observation of an optimum pulse separation of ∼1 ps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Advancing Artificial Intelligence with Liquid Argon Neutrino Experiments (Technical Report)

The grant allowed two main contributions: 1) The development of a first successful demonstration of the employment of Optimal Transport in liquid argon time projection chamber neutrino detectors. Optimal Transport, used in other contexts and specifically with LHC calorimetric data, was adapted to address a key particle identification challenge in LArTPCs: the separation of pi0 backgrounds from single-electrons produced in charged-current electron neutrino interactions. The work, leveraging ML methods such as k-nearest-neighbor (kNN) and support-vector-machine (SVM), showed an increase in background rejection of a factor of two or more. Work is now ongoing to incorporate this development in physics analyses for LArTPC experiments and more broadly expand the use of OT in LArTPC detectors including DUNE. This work was done in collaboration with the phenomenology group led by Nathaniel Craig at UCSB. 2) The deployment of NuGraph2, a graph neural network developed for LArTPC reconstruction, in the MicroBooNE experiment. NuGraph2 uses novel graph-neural-network methods on the rather simple LArTPC inputs of reconstructed hits, greatly simplifying the workflow compared to the use of waveform or signal-deconvolved wire ROIs. The network performed particle classification and was shown to address many challenging problems in LArTPC imaging including track-shower separation and the identification of protons and charged pions from primary muons. Our group collaborated with Giuseppe Cerati (FNAL scientist) who is one of the core developers of NuGraph2 to integrate this tool in MicroBooNE’s analysis framework. This consisted in tow key contributions: a) Studying performance on real data, which came with several months of iterations because the MC-trained version of the network was found to show significant bias that our group investigated and addressed. b) Integrating the output hit labeling of NuGraph2 into the existing particle tracking and shower reconstruction code. As a result of this work led by our team NuGraph2 is now enabling a suite of new analyses which benefit from enhanced capabilities and thus broader physics reach. The grant supported primarily the salary of UCSB graduate student Chuyue “Michaelia” Fang as well as partial summer salary support for PI Caratelli. Some funds were used for travel by Michaelia to ML related schools and conferences.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Physical Properties of Sub-galactic Clumps at 0.5 ≤ Z ≤ 1.5 in the UVUDF

We present an investigation of clumpy galaxies in the Hubble Ultra Deep Field at 0.5⩽z⩽1.5 in the rest-frame far-ultraviolet (FUV) using Hubble Space Telescope Wide Field Camera 3 broadband imaging in F225W, F275W, and F336W. An analysis of 1404 galaxies yields 209 galaxies that host 403 kpc scale clumps. These host galaxies appear to be typical star-forming galaxies, with an average of 2 clumps per galaxy and reaching a maximum of 8 clumps. We measure the photometry of the clumps and determine the mass, age, and star formation rates (SFR) using the spectral energy distribution fitting code FAST. We find that clumps make an average contribution of 19% to the total rest-frame FUV flux of their host galaxy. Individually, clumps contribute a median of 5% to the host galaxy SFR and an average of ∼4% to the host galaxy mass, with total clump contributions to the host galaxy stellar mass ranging widely from lower than 1% up to 93%. Clumps in the outskirts of galaxies are typically younger, with higher SFRs, than clumps in the inner regions. The results are consistent with clump migration theories in which clumps form through violent gravitational instabilities in gas-rich turbulent disks, eventually migrate toward the center of the galaxies, and coalesce into the bulge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Environmental Flow Requirements from FERC Licenses Across the US

Environmental flow requirements included in Federal Energy Regulatory Commission (FERC) hydropower licenses are important for balancing natural properties and benefits of river ecosystems (e.g., healthy species, recreation, water supply, flood control) supporting hydropower production. In some cases, environmental flow requirements may limit operational flexibility given current operational schemes and make a hydropower plant less able to provide power to the electric grid on demand. Hydropower plants may gain some flexibility as hydropower scheduling time periods are made to be more responsive to the short-term needs of an energy grid increasingly reliant on intermittent renewables. However, many flow requirements focus on the daily, monthly, or seasonal flow fluctuations which matches the time scale of most paradigms linking flow alterations to the health of river ecosystems. This dataset seeks to provide a greater understanding of how flexibility in environmental requirements can be leveraged to create positive outcomes for both the power system and the environment. It contains information on environmental flow requirements from the Protection, Mitigation, and Enhancement section of 50 randomly selected FERC licenses: 25 issued from 1998-2013 that were also included in the ORNL Mitigation Database (Schramm et al. 2015) and 25 issued from 2014-present. The information on environmental flow requirements was extracted from the PM&E section of 50 randomly selected FERC licenses: 25 issued from 1998-2013 that were also included in the ORNL Mitigation Database (Schramm et al. 2015) and 25 issued from 2014-present. The flow requirements were then categorized into flow augmentation categories based on whether the license stated a specific water management purpose for the given requirement called augmentation categories (i.e., fisheries or habitat, recreation or boating, industry, and general or unspecified; Table B). Requirements were also grouped into flow type categories (e.g., minimum flow rate, maximum flow rate, ramping rate). Additional information related to flow requirements such as the augmentation time-period and whether the flow rate was continuous (i.e., condition must be present at all-times) or instantaneous (i.e., condition present at a point in time) was also extracted from the licenses. Some licenses had specific flow requirements based on whether the project was in a wet, dry, or normal water year. If that information was presented in the license, it was also included in the data set. The location within the project was noted, hereafter, zone, in the dataset for flow requirements relating to specific areas of hydropower projects (Dam, Powerhouse, Bypass Reach). Maximum discharge capacities of hydropower facilities were also extracted from both the Existing Hydropower Assets (EHA) data set and the National Inventory of Dams (NID) databases. Each facility was coded with project identification codes from the EHA dataset to facilitate cross-referencing between datasets.

13 HYDRO ENERGY↗
DOE OSTI2021

Predictions of core plasma performance for the SPARC tokamak

SPARC is designed to be a high-field, medium-size tokamak aimed at achieving net energy gain with Ion Cyclotron Range-of-Frequencies (ICRF) as its primary auxiliary heating mechanism. Empirical predictions with conservative physics indicate that SPARC baseline plasmas would reach Q~11, well above its mission objective of Q>2. To build confidence that SPARC will be successful, physics-based integrated modeling has also been performed. The TRANSP code coupled with the theory-based TGLF turbulence model and EPED predictions for pedestal stability find that Q~9 is attainable in standard H-mode operation and confirms Q>2 operation is feasible even with adverse assumptions. In this analysis, ion cyclotron waves are simulated with the full wave TORIC code and alpha heating is modeled with the Monte-Carlo fast ion NUBEAM module. Detailed analysis of expected turbulence regimes with linear and nonlinear CGYRO simulations is also presented, demonstrating that profile predictions with the TGLF reduced model are in reasonable agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Development of the FHR advanced natural circulation analysis code and application to FHR safety analysis

The University of California, Berkeley (UCB) is performing thermal hydraulics safety analysis to develop the technical basis for design and licensing of fluoride-salt-cooled, high-temperature reactors (FHRs). FHR designs investigated by UCB use natural circulation for emergency, passive decay heat removal when normal decay heat removal systems fail. The FHR advanced natural circulation analysis (FANCY) code has been developed for assessment of passive decay heat removal capability and safety analysis of these innovative system designs. The FANCY code uses a one-dimensional, semi-implicit scheme to solve for pressure-linked mass, momentum and energy conservation equations. Graph theory is used to automatically generate a staggered mesh for complicated pipe network systems. Heat structure models have been implemented for three types of boundary conditions (Dirichlet, Neumann and Robin boundary conditions). Heat structures can be composed of several layers of different materials, and are used for simulation of heat structure temperature distribution and heat transfer rate. Control models are used to simulate sequences of events or trips of safety systems. A proportional-integral controller is also used to automatically make thermal hydraulic systems reach desired steady state conditions. A point kinetics model is used to model reactor kinetics behavior with temperature reactivity feedback. The underlying large sparse linear systems in these models are efficiently solved by using direct and iterative solvers provided by the SuperLU code on high performance machines. Input interfaces are designed to increase the flexibility of simulation for complicated thermal hydraulic systems. In conclusion, this paper mainly focuses on the methodology used to develop the FANCY code, and safety analysis of the Mark 1 pebble-bed FHR under development at UCB is performed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗