Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer 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 883 records · Page 49

An assessment of total depletion uncertainties for a GBC-32 dry storage cask with PLUS7 fuels using Monte Carlo sampling and two-way ANOVA F-test

It is important to accurately evaluate depletion uncertainties that exist due to the limitations of a computer code on accuracy in a nuclear depletion calculation for criticality safety analysis. In this work, the bias and bias uncertainty in the effective neutron multiplication factor (k{sub eff}) for the generic 32 PWR-assembly burnup credit (GBC-32) cask system were assessed using the Monte Carlo uncertainty sampling method for three different burnups and the isotopic biases and uncertainties obtained with the calculated and measured isotopic concentrations. The PLUS7 used nuclear fuel assemblies discharged from the Hanbit Nuclear Power Plant Unit 3 were assumed to be loaded into the GBC-32 cask. In addition, the statistical F-test in the two-way analysis of variance (ANOVA) was suggested to efficiently estimate the statistically converged total depletion uncertainty in k{sub eff}, with a significant reduction of computing efforts in the Monte Carlo criticality calculations, and it was applied to the GBC-32 cask system after a sensitivity analysis was performed to select the important nuclides having significant contributions to the total depletion uncertainty. With the statistical F-tests in the two-way ANOVA, 100 criticality calculations were needed to get a statistically converged total depletion uncertainty for 10 and 30 GWD/MTU, while 120 criticality calculations were needed for 50 GWD/MTU. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Investigations into validation of plutonium solutions for criticality safety analysis

Plutonium solution chemistry is among the most complex of all the elements in the periodic table. Under acidic conditions up to four different oxidation states can be present in significant amounts in a single solution. From the perspective of nuclear criticality safety, accurate information regarding the concentration of plutonium and other constituents in solution must be well-known in order to correctly model plutonium solution systems for analysis. This requires knowledge of the solution system along with the ability to make accurate assumptions and predictions of the chemical and physical properties of the solution system. If neutron transport methods are used to assess subcriticality in a nuclear criticality safety analysis, the models must be compared with benchmark experiments which are similar and the suitability of the computer code system must be quantified. This paper describes some investigations undertaken to support the validation of plutonium solution systems for nuclear criticality safety analysis and how the Whisper sensitivity/uncertainty (S/U) tool can provide information to support validation of such systems. Case study results for various solution systems such as aqueous plutonium nitrate, aqueous plutonium chloride and hypothetical metal-water mixture are presented. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Flow Field and Pressure Measurements in a 61-Pin Wire-Wrapped Bundle

Texas A and M University is conducting isothermal flow tests in a wire-wrapped, 61-pin hexagonal fuel bundle to support the research on advanced nuclear fuel development sponsored by the Department of Energy (DOE). The purpose of these tests is to produce high spatial- and temporal-resolution measurements of the flow and pressure measurements at different locations of the test bundle. The data resulting from these tests will be used to benchmark numerical simulations of the experiments using both industry-standard Reynolds-averaged Navier-Stokes (RANS) based Computation Fluid Dynamics (CFD) methods and a highly resolved Large Eddy Simulation (LES) computer code. This article described the main features of the experimental facility, its capabilities in producing experimental data of pressure drop and velocity fields at different axial locations and within different exterior and interior sub-channels, using advanced flow visualization techniques. Examples of experimental data outputs and typical locations of the measurements that will included in the research scope are also shown in this article. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evaluation of the Inclusive Electron Scattering Observables in the Resonance Region from the Experimental Data

The method for evaluation of inclusive electron scattering observables in the resonance region from the experimental data is presented. The computer codes in the web page allow user to evaluate inclusive electron scattering, virtual photon-proton cross sections and the F{sub 1} and F{sub 2} structure functions at the invariant masses of the final hadrons 1.0 < W < 4.0 GeV and at photon virtualities 0.5 < Q {sup 2} < 7.0 GeV{sup 2}. The developed approach is of interest for analyses of the experimental data on exploration of the hadron structure in the experiments with electromagnetic probes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

RISMC Industry Application no. 1 (ECCS/LOCA) - Core Characterization Automation: Lattice Codes Interface for PHISICS/RELAP5-3D

The U.S. NRC is currently proposing rulemaking '10 CFR 50.46c' to revise the loss-of-coolant-accident (LOCA) and emergency core cooling system (ECCS) acceptance criteria to include the effects of higher burn-up on cladding performance. The key implications of this proposition are that the cladding performance cannot be evaluated in isolation anymore. Both cladding and ECCS performance need to be considered in a coupled manner. This also implies that models for cladding performance as well as LOCA methodologies will need to be updated for regulatory purposes. Furthermore, efforts need to be made to find the limiting case within the regulatory input space. The Risk-Informed Safety Margin Characterization (RISMC) Pathway, as part of the DOE Light Water Reactor Sustainability (LWRS) research and development program initialed several programs at the Idaho National Laboratory (INL) to support the industry in the transition to the new rule and to offer potential solutions to this LOCA/ECCS analysis. One of these programs is the Industry Application no. 1 (IA1). Its purpose is to provide to the plant owner/operator a vehicle to inform decisions to manage margins related to compliance with the proposed 10 CFR 50.46c rule. IA1's goal is to develop an Integrated Evaluation Model (IEM) to revisit how uncertainties are propagated across the physical disciplines and data involved, as well as well as how risks are evaluated in a LOCA under 10 CFR 50.46c. This IEM is called LOCA Tool-kit for the US (LOTUS) and it connects five major disciplines involved in LOCA analysis. The focus of LOTUS is to establish the automation interfaces among the five disciplines. LOTUS will utilize in a first step existing computer codes. The risk-informed margins management approach for IA1 will provide a means of quantifying the impact on the key LOCA analysis figures-of-merit like peak cladding temperature (PCT), equivalent cladding reacted (ECR), etc. of a change in LOCA analysis inputs. The information that the risk analysis and associated tools provide can then be used for decision-making and margin management. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Simplex-stochastic collocation method with improved scalability

The Simplex-Stochastic Collocation (SSC) method is a robust tool used to propagate uncertain input distributions through a computer code. However, it becomes prohibitively expensive for problems with dimensions higher than 5. The main purpose of this paper is to identify bottlenecks, and to improve upon this bad scalability. In order to do so, we propose an alternative interpolation stencil technique based upon the Set-Covering problem, and we integrate the SSC method in the High-Dimensional Model-Reduction framework. In addition, we address the issue of ill-conditioned sample matrices, and we present an analytical map to facilitate uniformly-distributed simplex sampling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Neutrino: A Better Understanding Through Astrophysics: Final Report

The final report for the award "The Neutrino: A Better Understanding Through Astrophysics" is given. The goals of the work were the following: to construct new theoretical approaches to the problem of neutrino propagation in media including where neutrino-neutrino interactions are important; to pioneer the use of new approaches, including super-scattering operators, for the evolution of neutrino thermal and statistical ensembles; to implement these new approaches in computer codes to study neutrino evolution in supernovae and other hot, dense environments; to increase the realism of simulated signals of a Galactic supernovae neutrino burst in current and future neutrino detectors; to study the simulated signals to determine the ability to extract information on the missing neutrino mixing parameters and the dynamics of the supernova explosion; and to study sterile neutrinos and non-standard interactions of neutrinos in supernovae and their effect upon the signal. Accomplishments made in these areas are described.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

RELAP5-3D Pre-test Prediction for High Temperature Test Facility Test PG-26

The RELAP5-3D computer code has been used to perform a pre-test prediction for Test PG-26 in the High Temperature Test Facility at Oregon State University. Test PG-26 will simulate a pressurized conduction cooldown transient following a trip of the primary coolant system gas circulator. The facility is initialized with a steady-state power of 700 kW and prototypic core inlet and outlet helium coolant temperatures of 259 and 687°C, respectively. Sensitivity calculations addressed the flow coastdown time, the reactor cavity cooling system flow rate and temperature, and decay power.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Particle Resuspension Simulation Capability to Substantiate DOE-HDBK-3010 Data

The data provided in the U.S. Department of Energy (DOE) Handbook, DOE-HDBK-3010 provides data such as airborne release fractions and respirable fractions which are widely used in safety analysis and facility design. These fractions are often used as conservative bounds, despite the fact that the original data may have been obtained at a different scale or in a different configuration than the end application. The use of computational techniques to augment these data can be used to extend laboratory-scale data applicability to facility-scale decisions. Our prior report described an approach using Sandia National Laboratories (SNL) computer codes to simulate particle releases from various accident scenarios. The relevant physics controlling the particle resuspension phenomena can be highly problem-dependent. They can include surface tension for particles or surfaces with liquid present, van der Waals forces for small particles, electrostatic forces for charged particles, and mechanical interactions with surface features larger than the particles sizes. The complexity of this phenomenon and recent developments have been described in a recent review by Henry and Minier. (authors)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Converged accelerated finite difference scheme for the multigroup neutron diffusion equation

Computer codes involving neutron transport theory for nuclear engineering applications always require verification to assess improvement. Generally, analytical and semi-analytical benchmarks are desirable, since they are capable of high precision solutions to provide accurate standards of comparison. However, these benchmarks often involve relatively simple problems, usually assuming a certain degree of abstract modeling. In the present work, we show how semi-analytical equivalent benchmarks can be numerically generated using convergence acceleration. Specifically, we investigate the error behavior of a 1D spatial finite difference scheme for the multigroup (MG) steady-state neutron diffusion equation in plane geometry. Since solutions depending on subsequent discretization can be envisioned as terms of an infinite sequence converging to the true solution, extrapolation methods can accelerate an iterative process to obtain the limit before numerical instability sets in. The obtained results have been compared to the analytical solution to the 1D multigroup diffusion equation when available, using FORTRAN as the computational language. Finally, a slowing down problem has been solved using a cascading source update, showing how a finite difference scheme performs for ultra-fine groups (104 groups) in a reasonable computational time using convergence acceleration. (authors)

97 MATHEMATICS AND COMPUTING↗

ITER Port Interspace Pressure Calculations

The ITER ('the way' in Latin) project in Cadarache, France, will operate the largest tokamak ever built to demonstrate the physical and technical feasibility of nuclear fusion as an energy source. The ITER vacuum vessel is equipped with 50 access ports. Each port has an opening in the bio-shield that communicates with a dedicated port cell. During tokamak operation, the bio-shield opening must be closed with a concrete plug to shield the radiation coming from the plasma. This port plug separates the port cell into a port interspace (between vacuum vessel closure lid and port plug) on the inner side and a port cell on the outer side. This paper presents calculations of pressures and temperatures in the ITER port interspace after a double-ended guillotine break (DEGB) of a pipe in the tokamak cooling water system (TCWS). It is assumed that the DEGB occurs during the worst possible conditions-during water baking operation, with water at 523 K (250 deg. C) and 4.4 MPa. These conditions are more severe than during normal tokamak operation with the water at 398 K (125 deg. C) and 2 MPa. Two computer codes are employed in these calculations: RELAP5-3D Version 4.2.1 to calculate blowdown releases from the pipe break, and MELCOR, Version 1.8.6 to calculate pressures and temperatures in the port interspace. A sensitivity study has been performed to optimize the maximum peak pressure in the port interspace volume of ITER. The optimum case was for a labyrinth area (FL131) of 0.5 m{sup 2} and a relief panel (FL135) kept open after the first cycle. This case keeps the pressure in the port interspace under 150 kPa and the temperatures under 385 K. (authors)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Synthesis of polycaprolactone-titanium oxide multilayer films by nanosecond laser pulses and electrospinning technique for better implant fabrication

In this study, a combination of electrospinning and laser texturing is introduced as a novel method for increasing the biocompatibility of metal implants. Besides having a rough laser treated surface, the implant benefits from the high porosity and better wettability of an electrospun fibrous structure, which is a more favorable environment for cell proliferation. Titanium samples were patterned using a nanosecond laser beam and were placed as collectors in an electrospinning machine. They were then soaked in simulated body fluid for four weeks. Energy Dispersive X-ray and X-Ray Diffraction results indicate significantly more hydroxyapatite formation on laser treated samples with nanoscale fibers deposited on their surface. This shows that having a laser treated surface underneath the fibrous layer can improve short-term biocompatibility even before degradation of fibers. The thermal conductivity of the electrospun layer, measured using a Hot Disk Transient Plane Source instrument and computer code, was shown to be considerably lower than that of titanium and very close to bone. The presence of this layer can therefore be beneficial in making the implant more compatible to a biological medium. In case of dental implants, it was shown that this layer can act as a thermal barrier while a hot beverage is consumed and it can decrease the temperature rise by about 60%, which avoids any possible damage to newly formed cells during the healing period.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Study on Radionuclides Selection for DCGLs Development for Site Release of KRR Units 1 and 2 site in Korea - 20256

The selection of Derived Concentration Guideline Levels (DCGLs) in nuclear facilities decommissioning is important because they serve as criteria for site characterization and final status surveys. DCGLs should be derived for KAERI TRIGA Mark II and III research reactor site clearance and prior to the work, a radionuclides list for DCGLs consideration was created. Radionuclides selection for DCGLs was performed according to the procedure proposed by the US NRC based on the characterization survey related to the decommission site. By evaluating the dose contribution of each radionuclide using D and D computer code, radionuclides with a dose contribution fraction of 0.1% or more for the entire radionuclide were selected. As a result, the selected radionuclides were beta radionuclides H-3, C-14 and Sr-90 and gamma radionuclides Co-60, Cs-137, Eu-152 and Eu-154. The selected radionuclides can be presented for consultation with domestic regulatory agencies in the future and will be used as the technical base of DCGLs for the final status survey. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

Coupling Exascale Multiphysics Applications: Methods and Lessons Learned

With the growing computational complexity of science and the complexity of new and emerging hardware, it is time to re-evaluate the traditional monolithic design of computational codes. One new paradigm is constructing larger scientific computational experiments from the coupling of multiple individual scientific applications, each targeting their own physics, characteristic lengths, and/or scales. We present a framework constructed by leveraging capabilities such as in-memory communications, workflow scheduling on HPC resources, and continuous performance monitoring. This code coupling capability is demonstrated by a fusion science scenario, where differences between the plasma at the edges and at the core of a device have different physical descriptions. This infrastructure not only enables the coupling of the physics components, but it also connects in situ or online analysis, compression, and visualization that accelerate the time between a run and the analysis of the science content. Results from runs on Titan and Cori are presented as a demonstration.

Choi, Jong Youl↗

NSRD-11: Computational Capability to Substantiate DOE-HDBK-3010 Data.

Safety basis analysts throughout the U.S. Department of Energy (DOE) complex rely heavily on the information provided in the DOE Handbook, DOE - HDBK - 3010, Airborne Release Fractions/Rates and Respirable Fractions for Nonreactor Nuclear Facilities, to determine radionuclide source terms. In calculating source terms, analysts tend to use the DOE Handbook's bounding values on airborne release fractions (ARFs) and respirable fractions (RFs) for various categories of insults (representing potential accident release categories). This is typically due to both time constraints and the avoidance of regulatory critique. Unfortunately, these bounding ARFs/RFs represent extremely conservative values. Moreover, they were derived from very limited small-scale bench/laboratory experiments and/or from engineered judgment. Thus, the basis for the data may not be representative of the actual unique accident conditions and configurations being evaluated. The goal of this research is to develop a more accurate and defensible method to determine bounding values for the DOE Handbook using state-of-art multi-physics-based computer codes. This enables us to better understand the fundamental physics and phenomena associated with the types of accidents in the handbook. In this year, this research included improvements of the high-fidelity codes to model particle resuspension and multi-component evaporation for fire scenarios. We also began to model ceramic fragmentation experiments, and to reanalyze the liquid fire and powder release experiments that were done last year. The results show that the added physics better describes the fragmentation phenomena. Thus, this work provides a low-cost method to establish physics-justified safety bounds by taking into account specific geometries and conditions that may not have been previously measured and/or are too costly to perform.

54 ENVIRONMENTAL SCIENCES↗

Multi-cycle boiling water reactor fuel cycle optimization

In this work a new computer code, BWROPT (Boiling Water Reactor Optimization), is presented. BWROPT uses the Parallel Simulated Annealing (PSA) algorithm to solve the out-of-core optimization problem coupled with an in-core optimization that determines the optimum fuel loading pattern. However it uses a Haling power profile for the depletion instead of optimizing the operating strategy. The result of this optimization is the optimum new fuel inventory and the core loading pattern for the first cycle considered in the optimization. Several changes were made to the optimization algorithm with respect to other nuclear fuel cycle optimization codes that use PSA. Instead of using constant sampling probabilities for the solution perturbation types throughout the optimization as is usually done in PSA optimizations the sampling probabilities are varied to get a better solution and/or decrease runtime. The new fuel types available for use can be sorted into an array based on any number of parameters so that each parameter can be incremented or decremented, which allows for more precise fuel type selection compared to random sampling. Also, the results are sorted by the new fuel inventory of the first cycle for ease of comparing alternative solutions. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evaluation of 3 numerical methods for propulsion integration studies on transonic transport configurations

An effort has been undertaken at the NASA Langley Research Center to assess the capabilities of available computational methods for use in propulsion integration design studies of transonic transport aircraft, particularly of pylon/nacelle combinations which exhibit essentially no interference drag. The three computer codes selected represent state-of-the-art computational methods for analyzing complex configurations at subsonic and transonic flight conditions. These are: EULER, a finitie volume solution of the Euler equation; VSAERO, a panel solution of the Laplace equation; and PPW, a finite difference solution of the small disturbance transonic equations. In general, all three codes have certain capabilities that allow them to be of some value in predicting the flows about transport configurations, but all have limitations. Until more accurate methods are available, careful application and interpretation of the results of these codes are needed.

Yaros, S. F.↗

Evaluation of three numerical methods for propulsion integration studies on transonic transport configurations

An effort has been undertaken at the NASA Langley Research Center to assess the capabilities of available computational methods for use in propulsion integration design studies of transonic transport aircraft, particularly of pylon/nacelle combinations which exhibit essentially no interference drag. The three computer codes selected represent state-of-the-art computational methods for analyzing complex configurations at subsonic and transonic flight conditions. These are: EULER, a finite volume solution of the Euler equation; VSAERO, a panel solution of the Laplace equation; and PPW, a finite difference solution of the small disturbance transonic equations. In general, all three codes have certain capabilities that allow them to be of some value in predicting the flows about transport configurations, but all have limitations. Until more accurate methods are available, careful application and interpretation of the results of these codes are needed.

Yaros, Steven F.↗