Search NASA⌕ Search

Engineering topics

Feng, B.

Publications and source records attributed to Feng, B..

18 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Fast Reactor Physics Model Verification Studies using ARC and PyARC Workflows

PyARC was recently developed at Argonne National Laboratory to automate many of the tasks required in the ARC (Argonne Reactor Computation) fast reactor simulation workflow, from input file generation, code execution, data transfer between ARC codes, and output postprocessing. PyARC will likely be the path forward to train new users of the ARC codes with the goal of wide adoption by the national laboratories, academia, and industry. In particular, for the ANL-JAEA collaboration under the Civil Nuclear Working Group (CNWG) project agreement NE-01, PyARC will be used to model the Joyo and EBR-II reactors for comparisons with measured data and calculated results from JAEA (Task 3: Fast Reactor Fuel and Core). As an additional avenue for verification and validation, this report investigates the use of PyARC towards a variety of existing ARC-based reactor models, in order to understand its efficacy in replicating the behavior of base ARC codes and better understand any limitations within modeling realistic fast reactor problems. To this end, PyARC was used to model the Joyo MKI, RBEC Benchmark-M, PRISM Mod-B, and EBR-II Run 138B cores, and its results were compared to those from existing ARC-based models. It was found that for hexagonal-based geometries PyARC was able to replicate the behavior of ARC codes to within 10 pcm for small reactor cores, and ~150pcm difference in eigenvalue for larger cores. These discrepancies are attributed primarily to differences in local mesh refinement options between ARC and PyARC, which currently cannot be resolved with PyARC’s latest version (1.6.0). In some of these cases, PyARC was used to model steady-state problems with initial core compositions originating from a prior REBUS depletion calculation. While PyARC was not designed to support such steady-state calculations, workarounds were applied to replicate the behavior of ARC-based calculations as closely as possible. Thus, these results demonstrate the wide extent to which they can be applied to fast reactor problems while still providing immense benefit to the user in terms of automating and standardizing common routines within the fast reactor analysis workflow. This study concluded that PyARC will be suitable for modeling the steady-state conditions of the EBR-II and Joyo fast reactors as part of the CNWG project agreement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

The Virtual Test Bed (VTB) repository: a library of multiphysics reference reactor models using NEAMS tools

With the next generation of nuclear reactors under development, modeling and simulation (MS) tools are being developed by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in order to support their design, licensing, and future operation. Mirroring the physical test beds currently under construction (i.e., EBR-II and ZPPR), the Virtual Test Bed (VTB) was launched by the National Reactor Innovation Center (NRIC) in collaboration with NEAMS to support the advanced reactor community. This collaborative effort, which involves multiple teams at both Idaho National Laboratory and Argonne National Laboratory aims to use NEAMS tools to model a wide range of reactor designs. Those models are automatically tested to ensure their continued functionality as the tools are further developed. Examples are extensively documented, each acting as a tutorial for applying the relevant NEAMS tools to that reactor design. Currently, five advanced reactor types (with a total of eight specific design variants) are simulated by a variety of different models. These models range from steady-state, core multiphysics simulations to integrated plant analysis during loss-of flow transients. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Application of sensitivity analysis in DYMOND/Dakota to fuel cycle transition scenarios

The ability to perform sensitivity analysis has been enabled for the nuclear fuel cycle simulator DYMOND through its coupling with the design and analysis toolkit Dakota. To test and demonstrate these new capabilities, a transition scenario and multi-parameter study were devised. The transition scenario represents a partial transition from the US nuclear fleet to a closed fuel cycle with small modular LWRs and fast reactors fueled by reprocessed used nuclear fuel. Four uncertain parameters in this transition were studied – start date of reprocessing, total reprocessing capacity, the nuclear energy demand growth, and the rate at which the fast reactors are deployed – with respect to their impact on four response metrics. The responses – total natural uranium consumed, maximum annual enrichment capacity required, total disposed mass, and total cost of the nuclear fuel cycle – were chosen based on measures known to be of interest in transition scenarios and to be significantly impacted by the varying parameters. Furthermore, analysis of this study was performed both from the direct sampling and through surrogate models developed in Dakota to calculate the global sensitivity measures Sobol’ indices. This example application of this new capability showed that the most consequential parameter to most metrics was the share of new build capacity that is fast reactors. However, for the cost metric, the scaling factor of the energy demand growth was significant and had synergistic behavior with the fast reactor new build share.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Detailed analyses of a TRISO-fueled microreactor: Modeling of a Micro-Reactor System using NEAMS Tools

Our analysis provides preliminary confirmation of the feasibility of applying the NEAMS codes for multiphysics simulations to model load following transients and accidental scenarios. Detailed demonstration was performed on a simple unit-cell model to perform transient simulations, while verifying accuracy of the multiphysics solution obtained with the NEAMS codes. Results on the full core focused on showcasing the multiphysics workflow rather than detailed verification of the results, at least for this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Sensitivity and Uncertainty Quantification of Transition Scenario Simulations

This report documents the first collective attempt at developing and applying capabilities to quantify uncertainties, assess parametric sensitivities, and optimize multiple parameters and metrics in fuel cycle simulations generated by the SA&I Campaign. To do this, external codes that were designed to perform sensitivity analysis and uncertainty quantification (SA&UQ) needed to be coupled to the SA&I Campaign’s nuclear fuel cycle simulators (NFCS). In FY20, two approaches were pursued: 1) coupling Cyclus to an ORNL-internal code called MOT (Metaheuristic Optimization Tool) and 2) coupling DYMOND to the opensource SA&UQ tool kit Dakota. The primary objective of having these NFCS/SA&UQ coupled capabilities is to better inform DOE-NE and other stakeholders on the results generated from the NFCS. For a given set of fuel cycle strategies, policies, and technology assumptions that make up a fuel cycle scenario, these NFCS have traditionally been used by the SA&I Campaign to provide quantitative answers in terms of year-by-year mass flows, infrastructure requirements, costs, etc. With these newly developed coupled capabilities, the SA&I Campaign can now efficiently simulate hundreds or thousands of these scenarios, sample large ranges of parameters and assumptions, and use the unique features of the SA&UQ tools to process the data. This enables providing answers with known and propagated uncertainties, determining the sensitivity of important metrics to different parameters and assumptions, quantifying how much fuel cycle and technology parameters impact each other, and producing optimized fuel cycle strategies for single and multiple variables. To demonstrate these new capabilities, the Cyclus/MOT was used to model several scenarios ranging from simple fleet retirements to transitions to advanced reactors. Specifically, for a transition scenario from LWRs to SFRs and advanced LWRs, uncertainty quantification, sensitivity analysis, and optimization studies were applied to cases involving single and multiple parameter (input) and single and multiple metric (output) variations. In addition, a similar transition scenario was modeled to demonstrate how to optimize the reprocessing capacity parameter to minimize two performance metrics while taking into account uncertainties from two other parameters. Lastly, a depletion module based on SCALE/ORIGEN was added in Cyclus to simulate the third scenario that was designed to quantify the impact of the modeling assumption that all LWR used nuclear fuel have the same burnup. The newly developed DYMOND/Dakota capability was also applied to a transition scenario from the existing fleet to small modular reactors and fast reactors. This particular scenario involves not only explicit isotopic depletion via ORIGEN-2, but also includes multirecycling and utilizing the criticality search feature to determine the fresh fuel composition of recycled fuel, a feature unique to the DYMOND NFCS. A large database of simulations were run with 4 main parameters that were sampled: start date of reprocessing, reprocessing capacity, energy demand growth rate, and advanced reactor share of the fleet. The 4 main metrics were uranium consumption, enrichment requirements, waste generation, and levelized cost of electricity using data from the Cost Basis Report. The demonstrated SA&UQ results include those that inform on how to choose parameters to avoid “failed” scenarios, Sobol’ indices that inform on the importance of various parameters individually and synergistically, and Analysis of Variance (ANOVA) studies that decompose parameter ranges into groups and informs on whether variations are statistically significant.

Feng, B.↗
DOE OSTI2020

Multi-Physics Simulations of Molten Chloride Fast Reactor using Nek5000 and PROTEUS-NODAL

This report documents the FY20 work to demonstrate a specific multi-physics capability within NEAMS tools to model molten chloride fast reactors (MCFR). Specifically, the CFD code Nek5000 and the neutronics code PROTEUS-NODAL were coupled off-line to calculate the delayed neutron precursor distribution within the flowing fuel salts in an MCFR core. An MCFR benchmark problem and a 3D test problem based on a conceptual MSR design were successfully solved through the coupled simulation of Nek5000 and PROTEUS-NODAL. In the coupled simulations for these test problems, Nek5000 provided high-quality flow field information, and PROTEUS-NODAL performed the neutron diffusion calculation along with the precursor drift model. The impact of precursor drift on the core status was also investigated. Additionally, the precursor drift model implemented in PROTEUS-NODAL was verified by making use of the reference precursor distribution generated from Nek5000. Results of these test problems suggest that the multi-physics capability of Nek5000 and PROTEUS-NODAL was properly developed for tracking delayed neutron precursor drift.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Molten salt reactor core simulation with PROTEUS

The Nuclear Energy Advance Modeling and Simulation Program (NEAMS) of the US Department of Energy (DOE) aims to develop advanced simulation tools to aid the development of advanced nuclear technologies. The high-fidelity PROTEUS code was developed in the NEAMS program for reactor physics analysis. This paper tests the applicability of the PROTEUS code to simulate a Molten Salt Reactor (MSR). The OpenMC code based on the Monte Carlo Method was employed to generate the cross sections and to provide stochastic solutions for verification. The PROTEUS code was successfully employed to simulate the steady-state performance of an MSR represented by the ThorCon core design. The 2D lattice, 2D core, and 3D full core calculations performed by PROTEUS show good agreement with that predicted by OpenMC. Furthermore, for the nodal option within PROTEUS, the precursor drift model was also tested, and this effect resulted in a reduction of the k-eff by 300 pcm.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Simulation of Thermal Spectrum MSR using PROTEUS

In recent years, there has been renewed interest in Molten Salter Reactors (MSRs) for their potential advantages compared to reactors that rely on solid fuel. Given the current interest for MSRs, it is a good opportunity to demonstrate that the advanced tools developed in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program are useful for modeling MSRs. In addition, one objective of the program is to develop modeling and simulation capabilities to predict the performance and safety of all types of advanced nuclear reactor. For this study, PROTEUS [1], the primary reactor physics tool within the program, was assessed and used to model one of the preliminary core designs provided by ThorCon, one of the leading MSR companies from industry. The focus of this first step in the modeling and simulation work is to predict the core design’s steady state performance.

Fei, T.↗
DOE OSTI2019

Coupled nonlinear elasticity, plastic slip, twinning, and phase transformation in single crystal titanium for plate impact loading

High purity single crystal titanium (Ti) under shockwave loading in plate impact experiment is modelled and simulated, with help of the finite element. A thermodynamically consistent system of equations is formulated in the frame of large deformation to consider material anisotropy, rate dependence, and multi-physics including nonlinear elasticity, dislocation based plastic slip, deformation twinning, and phase transformation. A novel kinematics is proposed to consider phase transformation in the α twin variants and also to consider the dislocation based slip in all of components of parent material, primary twins and the ω phase. The thermodynamically consistent driving forces for plastic slip, twinning, and phase transformation are developed based upon the second law of thermodynamics. For nonlinear elasticity, the stiffness matrix is dependent on the volume fractions of all components, and the volumetric part of Cauchy stress is obtained from the nonlinear equation of state. Dislocation based plastic slip in multi-variant and multiphase heterogeneous materials is developed and interactions of dislocations among all slip modes are taken into account. A mechanism for dislocation density evolution during twinning and phase transformation is proposed. The shock loading along the [0001] and [10$\bar{1}1$ directions of single crystal high purity Ti is investigated computationally. Very good correspondence between simulation and experiment is obtained, which includes pole figures, volume fractions of components, free surface velocity, peak pressure, and phase transformation pressure. Multiple experimental phenomena are interpreted based upon the progression of dislocation slip, deformation twinning, and phase transformation. In compression with the [10$\bar{1}1$ crystal, a higher volume fraction in the primary twins but a lower secondary twin volume fraction in the [0001] crystal in experiment was observed. Finally, the higher propensity for phase transformation occurs in the [0001] crystal is reproduced. In addition to material texture, distributions of temperature, stresses, and plastic strains dependent on the impact loading directions are revealed.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Coupled elasticity, plastic slip, and twinning in single crystal titanium loaded by split-Hopkinson pressure bar

Here coupled elasticity, plastic slip, and twinning in high-purity single crystal titanium loaded by split-Hopkinson pressure bar (SHPB) are investigated, in the framework of large deformation and by utilizing the finite element method (FEM). A thermodynamically consistent system of equations for combined plastic slip and twinning is formulated. Novel kinematics is proposed to develop the driving forces for slip and twinning processes and to consider plastic slip in twins. Dislocation based crystal plasticity is proposed, with an emphasis on the interactions among all slip modes in multi-variant/multiphase heterogeneous materials. A mechanism for dislocation density evolution during twinning is proposed. The impact loading along the [0001] and [1$0\bar{1}$1] directions of single crystal high purity titanium is investigated. The evolution of stress-strain field, dislocation density, and volume fraction for each variant in the sample during the entire loading process are obtained and discussed in detail. Results in simulations show that for the [0001] specimen due to the existence of six symmetry primary twin variants with respect to the loading direction and fifteen different slip systems in each variant, the fields of stress and dislocation density are very homogenous in the sample. The simulation results are compared with experimental data for the [0001] and [1$0\bar{1}$1] specimens. The experimental results are explained and interpreted for both specimens, which include the pole figures, the stress-strain curves, and volume fraction of variants. The simulation results provide an important insight into mechanical responses of high-purity single crystal titanium under high rate loading, and this paper advances the mesoscale modeling and simulations in coupled plastic slip and twinning.

36 MATERIALS SCIENCE↗
NASA NTRS1980

Lensless matched spatial filter correlator experiments

By combining a matched spatial filter and holographic optical element on a single plate, a compact optical frequency plane correlator results. Experimental pattern recognition correlation data and a theoretical analysis with experimental verification of the system's shift invariance are provided.

Shen, M.↗
NASA NTRS1980

HOE/lensless matched spatial filter wavelength-scaling correlator

A scaling correlator optical pattern recognition system is described in which a lensless matched spatial filter (with the second Fourier transform lens and the matched spatial filter recorded on the same plate) is used with a first Fourier transform lens that is also an holographic optical element. The matched spatial filter is recorded at one wavelength and correlation is obtained at a second wavelength. Experimental demonstration and output correlation SNR data are reported, together with a comparison of the system's noise level using conventional optics and holographic elements.

Shen, M.↗
NASA NTRS1980

Laser diode lensless MSF-HOE correlator

Description of an optical pattern recognition system employing a laser-diode lensless MSF correlator based on the use of holographic optical elements (HOE). The system has size, cost, and weight advantages over conventional optical processors, and it is believed to be the first unification of MSF correlation, laser diode, and HOE technologies.

Caimi, F.↗