Search NASA⌕ Search

SEARCH · Search NASA

Results for “production simulation codes”

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 127 records · Page 7

Development of a Ne gas target for {sup 22}Na production by proton irradiation

The article presents the design and development of a neon gas target for the production of {sup 22}Na using a proton beam from the room temperature cyclotron in Variable Energy Cyclotron Centre, Kolkata. The target design is made to handle a beam power of 85 W (17 MeV, 5 μA). The design is based on simulation using the computer code FLUKA for the beam dump and CFD-CFX for target cooling. The target has been successfully used for the production of {sup 22}Na in a 6 day long 17 MeV, 5 μA proton irradiation run.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Development of Critical Experiments to Provide Validation Data for Multiphysics Coupling Methods (Final Report)

Computational simulation has become one of the most important tools in the modern toolbox for nuclear reactor design. Reactor simulations, like the reactors themselves, involve multiple physical processes including neutronics, fluid flow and heat transfer (thermal-hydraulics), thermo-mechanics, fuel behavior, chemistry, and balance of plant. Each process can affect the behavior of the others, and advanced computational tools have been developed that couple separate physics models together to simulate the feedback effects. The objective of this work is to generate experimental data to support the validation of multi-physics reactor codes. It specifically targets the SHARP toolset from the NEAMS Reactor Product Line, but it is intended to be useful in any set of coded that couples thermal-hydraulics with neutronics simulations. The project takes advantage of the inherent flexibility of the RPI Reactor Critical Facility (RCF), a low power open-pool research reactor. The standard configuration of the RCF is used to explore the effects of moderator temperature on neutronic behavior and reactivity. Because the low power of the reactor generates no appreciable temperature change in the moderator, the moderator temperature is artificially changed with two 18 kW electric immersion heaters. The temperature feedback effects can be observed by recording the change in the reactor power as the reactor goes from a supercritical state, to critical, to a subcritical state. To expand the temperature range and temperature rate of change that can be investigated, an experimental apparatus consisting of a heated water loop was designed and installed in the RCF. In this setup, water from a heated reservoir tank is pumped through a pipe placed through the center of the reactor core. The water in the reservoir can be heated to around 70 °C and pumped through the system. The pump has a variable flow rate and electronic control, so a number of scenarios can be developed. A slug of water can be pushed through the system to generate a reasonably fast change in reactivity. The pump speed can be increased or decreased over the course of a measurement, or even oscillated at different periods, to observe the response of the reactor. To complement the experimental data generated by the RCF, a set of computational models have been developed. They are intended to illustrate the usefulness of the experimental data generated in validating coupled multi-physics simulations. A model previously constructed in MCNP has been used in conjunction with safety analysis at the RCF for many years and was modified to support the design of the heated water loop experiment. This model was translated into Serpent 2, a similar Monte Carlo code, which provided the added advantage of being able to produce collapsed cross-section data files for use in the deterministic neutron transport codes that are part of the SHARP toolset. A complete reconstruction of the RCF model, both with and without the heated fluid loop experiment, were necessary for use with the SHARP codes, PROTEUS (neutronics) and Nek5000 (thermal hydraulics), as these require an unstructured mesh format to operate. The computational requirements for fully coupled three dimensional simulations are significant, and challenges to full implementation and comparison with RCF experiments remain. A wide variety of experiments have been conducted for eventual benchmarking and validation of multi-physics reactor codes. A complete database of detailed time- and temperature-dependent measurements have been generated and stored. Several examples of the experimental results are presented here, with refinement and analysis of experiments ongoing. While these experiments still only address a finite range of operating conditions, in particular being limited to water in the sub-boiling regime, they represent a unique (to date) look at the feedback between temperature and neutronic behavior suited for validation of the coupling routines in multi-physics codes. It is hoped that these will be the foundation of future benchmarks that will provide confidence in the tools that will be used to develop the next generation of nuclear reactor systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulation of Gas Production from Multilayered Hydrate-Bearing Media with Fully Coupled Flow, Thermal, Chemical and Geomechanical Processes Using TOUGH+Millstone. Part 2: Geomechanical Formulation and Numerical Coupling

The TOUGH+Millstone simulator has been developed for the analysis of coupled flow, thermal and geomechanical processes associated with the formation and/or dissociation of CH 4 hydrates in geological media. It is composed of two constituent codes: (a) a significantly enhanced version of the TOUGH+HYDRATE simulator, v2.0, that accounts for all known flow, physical, thermodynamic and chemical processes associated with the behavior of hydrate-bearing systems undergoing changes and includes the most recent advances in the description of the system properties, coupled seamlessly with (b) Millstone v1.0, a new code that addresses the conceptual, computational and mathematical shortcomings of earlier codes used to describe the geomechanical response of these systems. The capabilities of the TOUGH+Millstone code are demonstrated in the simulation and analysis of the system flow, thermal and geomechanical behavior during gas production from a realistic complex offshore hydrate deposit. In the second part of this series, we describe the Millstone geomechanical simulator. The hydrate-dependent, rate-based poromechanical formulation is presented and solved using a finite element discretization. A novel multimesh coupling scheme is introduced, wherein interpolators are automatically built to transfer data between the finite difference discretization of TOUGH+ and the finite element discretization of Millstone. We provide verification examples against analytic solutions for poroelasticity and a simplified demonstration problem for mechanically induced phase change in a hydrate sediment.

58 GEOSCIENCES↗

SimApp 1.8.x

System for converting DESIRE (Direct Executing Simulation in Real Time) simulation models to executable code models (C# / C++). These converted models are used as the basis for 2 products, one of which is the focus of this review: an engineering design interface for Fissile Solution Systems. ModelConverter.exe (C#) converts DESIRE model inputs to a C++ plug-in subclass. SimApp.exe (C# / C++) has a C++ dll that is a numerical integration engine that the Converted models plug into. The engineering design interface is a C# GUI that calls on the C++ dll to execute a particular converted model. The GUI allows the user to specify the length of a simulation, run the simulation, view data plots and data tables. Stability plots may also be produced and manipulated. The engineering design interface will be released as source code and executable; the model converter source is also intended for release.

Determan, JohnC.↗

SimApp 1.8.x

System for converting DESIRE (Direct Executing Simulation in Real Time) simulation models to executable code models (C# / C++). These converted models are used as the basis for 2 products, one of which is the focus of this review: an engineering design interface for Fissile Solution Systems. ModelConverter.exe (C#) converts DESIRE model inputs to a C++ plug-in subclass. SimApp.exe (C# / C++) has a C++ dll that is a numerical integration engine that the Converted models plug into. The engineering design interface is a C# GUI that calls on the C++ dll to execute a particular converted model. The GUI allows the user to specify the length of a simulation, run the simulation, view data plots and data tables. Stability plots may also be produced and manipulated. The engineering design interface will be released as source code and executable; the model converter source is also intended for release.

Determan, John C.↗

Time-Dependent Boundary Modeling to Inform Design of SPARC Diagnostic and Actuators

The purpose of this work was to inform the design of diagnostic and actuator systems for a new experimental facility for Commonwealth Fusion Systems (CFS), called SPARC [Creely2020], using time-dependent plasma boundary simulations with the Scrape Off Layer Plasma Simulator (SOLPS) code [Wiesen2015]. The SPARC tokamak is a device that aims to demonstrate fusion energy production and high-field plasma operating scenarios. The developed technologies will be incorporated into a net-electricity pilot that will demonstrate approximately 200MW of electric power. A major difference from other private industry fusion companies is the use of a standard aspect ratio, but at a high magnetic field strength, which provides a well-established physics basis and a modest extrapolation to reactor operation. Time-dependent simulations were used to inform the design of main-chamber plasma facing components, which are protected by feedback control of the heat and particle loads during high power operation. The simulations provide the timescales and magnitudes of plasma and neutral particle response useful for determining the position and operation of diagnostics and actuators that will enable control. The plasma response and phase-space diagrams also provide input into advanced model-based control schemes, which are envisioned for later SPARC operation. Model-based control can reduce the risk of component failure, which would require expensive in-situ repairs and the associated delays, by predicting the system response to actuators over a short time-horizon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Activation of lead coolant in Westinghouse lead fast reactor

The Westinghouse Lead-cooled Fast Reactor (LFR) is a medium-sized, passively safe, scalable reactor being developed by Westinghouse in collaboration with domestic and international organizations with the primary goals of reducing front-end capital cost and generating flexible and cost-competitive electricity while satisfying the highest standards in terms of safety and sustainability. An important aspect of new plant design is an assessment of activation level of its components so that an appropriate nuclear waste management plan can be established. As lead coolant is a distinguishing design feature of the LFR, its activation due to irradiation and the implications on establishing an appropriate waste disposal approach are important for engagement of project stakeholders, customers and nuclear safety regulators. This paper provides a brief description of the Westinghouse LFR design, a background on lead activation phenomenology, and then describes the methodology used for lead activation analysis performed with the computer code Serpent 2. It was determined that impurities naturally occurring in lead have the biggest contribution to lead activation. Corrosion and erosion products from the components immersed in lead can also be important contributors to lead activity and their impact can be limited by appropriate selection of materials and operating conditions. Polonium isotopes, resulting from bismuth generated from lead irradiation and its impurities, bring a notable contribution to the inhalation and ingestion doses during the plant operation and in the first years after the plant shutdown. However, it is not by far as important as for lead-bismuth coolant.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Comparative study of MacCormack and TVD MacCormack schemes for three-dimensional separation at wing/body junctions in supersonic flows

A robust, discontinuity-resolving TVD MacCormack scheme containing no dependent parameters requiring adjustment is presently used to investigate the 3D separation of wing/body junction flows at supersonic speeds. Many production codes employing MacCormack schemes can be adapted to use this method. A numerical simulation of laminar supersonic junction flow is found to yield improved separation location predictions, as well as the axial velocity profiles in the separated flow region.

Lakshmanan, Balakrishnan↗

Ares I Flight Control System Overview

This paper describes the control challenges posed by the Ares I vehicle, the flight control system design and performance analyses used to test and verify the design. The major challenges in developing the control system are structural dynamics, dynamic effects from the powerful first stage booster, aerodynamics, first stage separation and large uncertainties in the dynamic models for all these. Classical control techniques were employed using innovative methods for structural mode filter design and an anti-drift feature to compensate for translational and rotational disturbances. This design was coded into an integrated vehicle flight simulation and tested by Monte Carlo methods. The product of this effort is a linear, robust controller design that is easy to implement, verify and test.

Hall, Charles↗

A Preliminary Study on the Toxic Combustion Products Testing of Polymers Used in High-Pressure Oxygen Systems

One likely cause of polymer ignition in a high-pressure oxygen system is adiabatic-compression heating of polymers caused by pneumatic impact. Oxidative _ pyrolysis or combustion of polymers in a high-pressure oxygen system could generate toxic gases. This paper reports the preliminary results of toxic combustion product testing of selected polymers in a pneumatic-impact test system. Five polymers commonly used in high-pressure oxygen systems, Nylon 6/6, polychlorotrifluoroethylene (CTFE), polytetrafluoroethylene (PTFE), fluoroelastomer (Viton(TradeMark) A), and nitrile rubber (Buna N), were tested in a pneumatic-impact test system at 2500- or 3500-psia oxygen pressure. The polymers were ignited and burned, then combustion products were collected in a stainless-steel sample bottle and analyzed by GC/MS/IRD, GC/FID, and GC/Methanizer/FID. The results of adiabatic-compression tests show that combustion of hydrocarbon polymers, nitrogen-containing polymers, and halogenated polymers in high-pressure oxygen systems are relatively complete. Toxicity of the combustion product gas is presumably much lower than the combustion product gas generated from ambient-pressure oxygen (or air) environments. The NASA-Lewis equilibrium code was used to determine the composition of combustion product gas generated from a simulated, adiabatic-compression test of nine polymers. The results are presented and discussed.

Hshieh, Fu-Yu↗

Comparison between Hayabusa 2 Spectral Measurements and Simulations

This paper compares recent Hayabusa 2 spectral measurements with state-of-the-art shock-layer radiation simulations resulting from the LAURA/HARA code suite. These simulations include coupled ablation, which accounts for the injection of ablation products into the flowfield, and coupled radiation, which accounts for radiative energy loss in the flowfield. To enable the coupled ablation simulations, a best-estimate model is developed for Hayabusa’s carbon-phenolic ablator, based on the limited available published information. The comparison between the simulations and measurements focuses on two atomic nitrogen lines and two atomic oxygen lines, as well as the CN Violet band system. For the atomic lines, the measurements and simulations agree within 25% over most of the trajectory. This excellent agreement is unprecedented for observed radiation measurements of atomic lines, where previous Stardust and Hayabusa 1 comparisons were significantly worse. The improved agreement for these Hayabusa 2 comparisons is both the result of improved measurement quality and enhanced flowfield/radiation modeling. To provide a link between these observed radiation measurements and the radiative heating to a vehicle surface, the recently developed flowfield property binning approach is used to identify flowfield properties that provide the dominant emission contribution to the measured spectrum. These identified flowfield properties are then shown to match those for the surface radiative heating to the currently developed Mars Sample Return (MSR) Earth Entry System (EES) at a specific trajectory point and surface location. This matching of flowfield properties indicates that this EES radiative heating and Hayabusa 2 observed radiation simulations are equivalent radiation problems. Therefore, the excellent agreement between Hayabusa 2 measurements and simulations may be leveraged to inform the radiation heating margin for the MSR EES.

Christopher O Johnston↗

Comparison between Hayabusa 2 Spectral Measurements and Simulations

Recent Hayabusa 2 spectral measurements are compared with state-of-the-art shock-layer radiation simulations resulting from the LAURA/HARA code suite. These simulations include coupled ablation, which accounts for the injection of ablation products into the flowfield, and coupled radiation, which accounts for radiative energy loss in the flowfield. To enable the cou- pled ablation simulations, a best-estimate model is developed for Hayabusa’s carbon-phenolic ablator based on available published information. The comparison between the simulations and measurements focuses on three atomic nitrogen lines, two atomic oxygen lines, one atomic hydrogen line, and the CN Violet band system. Parametric uncertainties are evaluated for these features using a recently developed capability in LAURA/HARA. These uncertainties capture the impact of flowfield and radiation modeling uncertainties on the observed radiation simulations, and include the impact of coupled radiation and ablation product emission. The resulting parametric uncertainties range from 50 to over 100% for the atomic lines, and 80% for CN Violet. Uncertainties in flowfield kinetics and non-Boltzmann modeling provide the dominant contribution to these uncertainties. For the atomic nitrogen and oxygen lines, com- parisons between the nominal simulations and measurements show agreement well within these simulation uncertainty bounds, with the measurements within 20% of the nominal simulation over most of the trajectory, and within 10% at peak emission. This excellent agreement is un- precedented for observed radiation measurements of atomic lines, where previous Stardust and Hayabusa 1 comparisons were significantly worse. The improved agreement for these Hayabusa 2 comparisons is both the result of improved measurement quality and enhanced flowfield and radiation modeling. For the ablation product emission from the 656 nm atomic hydrogen line, the agreement between measurements and simulations is within ±20% for all trajectory points except one outlier. This good comparison is surprising considering this species depends on the ablation rate and is shown to have strong wake contribution, which both contribute to the large ≈±70% parametric uncertainty. For the ablation product emission from CN Violet, the simulations over-predict the measurements by over a factor of two early in the trajectory, with the comparison improving later in the trajectory. This disagreement for CN Violet is consistent with the ≈±80% parametric uncertainty bounds evaluated for the simulations. To utilize the excellent agreement between the atomic line simulations and measurements in the assessment of the radiative heating margin for a future flight vehicle, the flowfield property binning (FPB) approach is applied. Using the Mars Sample Return (MSR) Earth Entry System (EES) as an example, the combination of the excellent Hayabusa comparisons and the FPB analysis justify nearly a 10% reduction in the radiative heating margin from an assumed baseline margin of 30%. This direct quantitative link between observed radiation measurements and the radia- tive heating margin of future flight vehicles demonstrates the value of these observed radiation measurements.

Shock-Layer Radiative Heating↗

Analysis and Simulation of a Blue Energy Cycle

The mixing process of fresh water and seawater releases a significant amount of energy and is a potential source of renewable energy. The so called ‘blue energy’ or salinity-gradient energy can be harvested by a device consisting of carbon electrodes immersed in an electrolyte solution, based on the principle of capacitive double layer expansion (CDLE). In this study, we have investigated the feasibility of energy production based on the CDLE principle. Experiments and computer simulations were used to study the process. Mesoporous carbon materials, synthesized at the Oak Ridge National Laboratory, were used as electrode materials in the experiments. Neutron imaging of the blue energy cycle was conducted with cylindrical mesoporous carbon electrodes and 0.5 M lithium chloride as the electrolyte solution. For experiments conducted at 0.6 V and 0.9 V applied potential, a voltage increase of 0.061 V and 0.054 V was observed, respectively. From sequences of neutron images obtained for each step of the blue energy cycle, information on the direction and magnitude of lithium ion transport was obtained. A computer code was developed to simulate the process. Experimental data and computer simulations allowed us to predict energy production.

30 DIRECT ENERGY CONVERSION↗

Computer simulation of auroral kilometric radiation

An investigation is carried out of the linear amplification and nonlinear saturation of electromagnetic waves associated with two components of electrons consisting of a cold Maxwellian background and a weakly relativistic population of electrons possessing a loss-cone distribution. The goal is to understand the mechanism responsible for the production of auroral kilometric radiation (AKR). A relativistic 1-2/2 dimensional electromagnetic particle simulation code is used in studying an initial value problem modeling the AKR source region. The simulation makes it possible to follow the growth of all modes of radiation, including the ordinary mode (O-mode), slow extraordinary mode (Z-mode), whistler mode and fast extraordinary mode (X-mode), past the point of saturation.

Wagner, J. S.↗

Using the Diurnal Variability in GeoNEX TOA Reflectances for Earth Monitoring

Observations from the third-generation geostationary satellite instruments (GOES 16/17 ABI, Himawari 8/9 AHI, and etc.) have spatial resolution and spectral band configurations comparable to flagship LEO sensors (e.g., MODIS/VIIRS). More importantly, these data are acquired at very high temporal resolution, faithfully recording the variations of the full disk of Earth at every 5-10 minutes. They thus provide unique information about Earth’s atmosphere and surface. In order to explore the unique information content of geostationary data, this study systematically analyzes the diurnal variability in the GeoNEX L1G TOA reflectance products and compares them to simulated results by state-of-the-art radiative transfer codes. Our results show that • The smoothness of the TOA reflectance diurnal cycle provides a convenient and reliable way to identify stable atmospheric conditions and filter out passing clouds/shadows. • The diurnal variability of the blue band (0.47µm) reflectance is regulated mainly by atmospheric optical conditions over a majority of land cover types. As such, the diurnal variability of the blue band data allows us to retrieve AOD without invoking the use of spectral band ratios (SRC) as in previous algorithms. • In comparison, the diurnal variability of the short-wave infrared band (2.2µm) BRFs is mainly regulated by surface reflectance and the sun-target-satellite geometry. This information allows us to test and, if suitable, retrieve surface BRDF parameters. • Spectral band ratios, especially those between the 2.2µm and 0.47µm bands, are not constant but vary by locations and sun-target-satellite geometries. Our analysis clearly demonstrates that the information provided in high-frequent geostationary observations is unique and complementary to LEO sensors. Therefore, a synergy of GEO and LEO (and other) sensors has the great potential to improve existing remote sensing models and algorithms for better Earth monitoring.

Diurnal Variability↗

Experimental validation of a nuclear forensics methodology for source reactor-type discrimination of chemically separated plutonium

An experimental validation of a nuclear forensics methodology for the source reactor-type discrimination of separated weapons-useable plutonium is presented. The methodology uses measured values of intra-element isotope ratios of plutonium and fission product contaminants. MCNP radiation transport codes were used for various reactor core modeling and fuel burnup simulations. A reactor-dependent library of intra-element isotope ratio values as a function of burnup and time since irradiation was created from the simulation results. The experimental validation of the methodology was achieved by performing two low-burnup experimental irradiations, resulting in distinct fuel samples containing sub-milligram quantities of weapons-useable plutonium. The irradiated samples were subjected to gamma and mass spectrometry to measure several intra-element isotope ratios. For each reactor in the library, a maximum likelihood calculation was utilized to compare the measured and simulated intra-element isotope ratio values, producing a likelihood value which is proportional to the probability of observing the measured ratio values, given a particular reactor in the library. The measured intra-element isotope ratio values of both irradiated samples and its comparison with the simulation predictions using maximum likelihood analyses are presented. Finally, the analyses validate the nuclear forensics methodology developed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Spatial calibration and synthetic diagnostic of a multi-energy hard x-ray camera at WEST tokamak

WEST (tungsten environment in steady-state tokamak) is starting operation for the first time with a water-cooled full tungsten divertor, enabling long pulse operation. Heating is provided by radiofrequency systems, including lower hybrid current drive (LHCD). In this context, a compact multi-energy hard x-ray camera has been installed for energy and space-resolved measurements of the electron temperature, the fast electron tail density produced by LHCD and runaway electrons, and the beam–target emission of tungsten at the target due to fast electron losses interacting with the divertor plates. The diagnostic is a pinhole camera based on a 2D pixel array detector (Pilatus 3 CdTe CMOS Hybrid-Pixel detector produced by DECTRIS). The novelty of this diagnostic technique is the detector’s capability of adjusting the threshold energy at pixel level. This innovation provides great flexibility in the energy configuration, allowing simultaneous space and energy-resolved x-ray measurements. This contribution details two important steps in the preparation of the diagnostic operation. First, the in-vessel spatial calibration that was carried out with a radioactive source. Second, the synthetic diagnostic is obtained by the suite of codes ALOHA/C3PO/LUKE/R5-X2, which simulates LH wave propagation and absorption, as well as the fast electron bremsstrahlung production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental studies of Micro- and Nano-grained UO 2 : Grain Growth Behavior, Sufrace Morphology, and Fracture Toughness

This activity is supported by the US Nuclear Energy Advanced Modeling and Simulation (NEAMS) Fuels Product Line (FPL) and aims at providing experimental data for the validation of the mesoscale simulation code MARMOT. MARMOT is a mesoscale multiphysics code that predicts the coevolution of microstructure and properties within reactor fuel during its lifetime in the reactor. It is an important component of the Moose-Bison-Marmot (MBM) code suite that has been developed by Idaho National Laboratory (INL) to enable next generation fuel performance modeling capability as part of the NEAMS Program FPL. In order to ensure the accuracy of the microstructure-based materials models being developed within the MARMOT code, extensive validation efforts must be carried out. In this report, we summarize the experimental efforts in FY16 including the following important experiments: (1) in-situ grain growth measurement of nano-grained UO 2 ; (2) investigation of surface morphology in micrograined UO 2 ; (3) Nano-indentation experiments on nano- and micro-grained UO 2 . The highlight of this year is: we have successfully demonstrated our capability to in-situ measure grain size development while maintaining the stoichiometry of nano-grained UO 2 materials; the experiment is, for the first time, using synchrotron X-ray diffraction to in-situ measure grain growth behavior of UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗