Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Criticality and axial offset searches based on the integrated neutron balance approach

Criticality and axial offset (AO) searches are key features in the simulation of plant operation procedures. Typically, the dedicated algorithms consist of either a Newton method or a secant method that features the successive computation of difference quotients to be used as derivatives. Within this context, computational robustness and efficiency are of paramount importance. This applies in particular to when the variations imposed during the search are spatially heterogeneous, such as in the case of tuning control rod positions for achieving target AO. In order to optimize this kind of search capability in AREVA NP's reactor code ARTEMIS in accordance with this design principle, a generalized approach has been implemented that harmonizes very well numerically with the overall multi-physics iterative solution process. As embedded in the latter, the new neutronic integral re-balancing approach is defined through periodic whole-core space-energy integrations of the principal terms in the core-integrated process-rate balance terms (i.e. neutron absorption, production and leakage). This procedure yields sequences of single zero-dimensional equations from which the chosen tuning parameters can be solved directly in dependence of the imposed (and systematically fixed) values for the response quantities k{sub eff} and/or AO. The converged result of the iteration sequence of successively obtained search parameter values is final in terms of being fully consistent with the entire set of multi-physics equations while enabling the accurate fulfillment of the target response value. The k etc and AO searches can be pursued simultaneously. Judging from the results of pursued verifications, the neutronic integral re-balancing approach fulfils the above-mentioned expectations convincingly. Specific verification examples are presented, such as the determination of the insertion depth of a critical control bank, a double search on target criticality and target AO by adjusting the boron concentration and the position of a control bank, and a critical power search during reactor cycle stretch out. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Cluster dynamics simulation of xenon diffusion during irradiation in UO 2

Diffusion of fission gas in UO 2 nuclear fuel impacts several important performance metrics, such as fission gas release, swelling, and thermal conductivity. Current empirical models of fission gas release have significant uncertainty, some of which derives from the bulk diffusion rate and its dependence on, for example, fuel chemistry and irradiation. In this work, we have applied the previously-developed Free Energy Cluster Dynamics (FECD) methodology in the code Centipede to calculate xenon cluster concentrations in UO 2 under intrinsic (high temperature) and irradiation-enhanced (intermediate temperature) conditions in order to develop a model of the xenon diffusion coefficient based on the atomic scale mechanisms responsible for transport. While the diffusion mechanism for xenon in UO 2 is adequately described by the Xe + U 2 O vacancy cluster for intrinsic conditions, a similar process is not capable of capturing measured in-pile fission gas diffusivity at intermediate temperatures. Therefore, a different diffusion mechanism must dominate under this regime. Using calculated atomistic data, we have shown that irradiation-enhanced diffusion at intermediate temperatures occurs via the larger Xe + U 4 O y vacancy clusters, which have lower migration barriers and increase in concentration by several orders of magnitude compared to intrinsic conditions. This mechanism is enabled by the increased uranium vacancy concentration under irradiation due to Frenkel pair production. In addition, the fast migration of uranium interstitials with two attached oxygen interstitials lowers the total uranium interstitial concentration through reactions with sinks. This allows the extended defects, such as Xe + U 4 O y vacancy clusters, to maintain high concentrations by limiting annihilation with attached vacancies. Furthermore, predictions using the Xe + U 4 O y diffusion mechanism are in good agreement with experiment, albeit with some differences in the Arrhenius slope, which we believe may be related to either experimental or model parameter uncertainty. Lastly, an analytical expression suitable for application in fuel performance simulations was derived to capture the predictions of the Centipede simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Structure of oxides prepared by decomposition of layered double Mg–Al and Ni–Al hydroxides

Abstracts: Thermal decomposition of Mg–Al and Ni–Al layered double hydroxides LDH at temperatures lower than 800 °C leads to the formation of oxides with different structures. Mg–Al oxide has a very defective structure and consists of octahedral layers as in periclase MgO and mixed octahedral–tetrahedral layers as in spinel MgAl{sub 2}O{sub 4}. Mixed Ni–Al oxide has a sandwich-like structure, consisting of a core with Al-doped NiO-like structure and some surface layers with spinel NiAl{sub 2}O{sub 4} structure epitaxial connected with the core. Suggested models were verified by simulation of X-ray diffraction patterns using DIFFaX code, as well as HRTEM, IR-, UV-spectroscopies, and XPS. - Graphical abstract: In the Mg–Al layered double hydroxide Al{sup 3+} ions migrate into interlayers during decomposition. The Mg–Al oxide represents sequence of octahedral and octahedral–tetrahedral spinel layers with vacancies. The Ni–Al oxide has a sandwich-like structure with NiO-like core and surface spinel layers as a result of migration of Al{sup 3+} ions on the surface. The models explain the presence and absence of “memory effect” for the Mg–Al and Ni–Al oxides, respectively. - Highlights: • We study products of Mg(Ni)–Al LDH decomposition by calcination at 500(400)–800 °C. • In Mg–Al/Ni–Al LDH Al ions migrate into interlayers/on the surface during decomposition. • Mg–Al oxide represents sequence of periclase- and spinel-like layers with vacancies. • Ni–Al oxide has a sandwich-like structure with NiO-like core and surface spinel layers. • The models explain the presence/absence of “memory effect” for Mg–Al/Ni–Al oxides.

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

NEML2 - THE NEW ENGINEERING MATERIAL MODEL LIBRARY, VERSION 2

SF-22-157 NEML2 is an extension of an earlier open source code developed at ANL (SF-xxxx) called NEML. These libraries provide tools for constructing and using models of engineering materials, with a focus on high temperature materials for nuclear reactor applications. The software include tools for setting up material models, calibrating models, storing models for later use in other simulations, and running simple "material point" simulations of the material behavior. The main difference between NEML and NEML2 is the new software has been complete restructured to run on GPUs and other accelerators in order to improve the performance of methods used to train constitutive models against data and to improve the performance of the final models in production environments, like finite element software. While over time we expect NEML2 to become more-or-less feature complete with NEML we are filing a new software disclosure because (1) the library has been complete restructured and (2) we want to host NEML2 in a separate repository.

MESSNER, MARK↗
NASA NTRS1985

An expanding vortex site for the r-process in rotating stellar collapse

The astrophysical site of r-process nucleosynthesis is investigated theoretically using two-dimensional expanding-vortex stellar-collapse simulations based on the Eulerian adaptive-mesh MHD code of Symbalisty (1984). The results are presented graphically, and it is found that the classical r-process can be explained as the result of the collapse of the highly rotating iron core of a 15-solar-mass star, in which the angular velocity along the rotation axis reaches a maximum of 1000 rad/s, corresponding to an angular momentum of 3.4 x 10 to the 48th erg/s for the inner 1.5 solar mass. The ejected jets are shown to yield about 0.0004 solar mass per supernova, sufficient to explain the observed abundances of r-process products.

Symbalisty, E. M. D.↗
NASA NTRS1988

Numerical simulation of hypersonic inlet flows with equilibrium or finite rate chemistry

An efficient numerical program incorporated with comprehensive high temperature gas property models has been developed to simulate hypersonic inlet flows. The computer program employs an implicit lower-upper time marching scheme to solve the two-dimensional Navier-Stokes equations with variable thermodynamic and transport properties. Both finite-rate and local-equilibrium approaches are adopted in the chemical reaction model for dissociation and ionization of the inlet air. In the finite rate approach, eleven species equations coupled with fluid dynamic equations are solved simultaneously. In the local-equilibrium approach, instead of solving species equations, an efficient chemical equilibrium package has been developed and incorporated into the flow code to obtain chemical compositions directly. Gas properties for the reaction products species are calculated by methods of statistical mechanics and fit to a polynomial form for C(p). In the present study, since the chemical reaction time is comparable to the flow residence time, the local-equilibrium model underpredicts the temperature in the shock layer. Significant differences of predicted chemical compositions in shock layer between finite rate and local-equilibrium approaches have been observed.

Yu, Sheng-Tao↗
NASA NTRS2011

Contributions of TetrUSS to Project Orion

The NASA Constellation program has relied heavily on Computational Fluid Dynamics simulations for generating aerodynamic databases and design loads. The Orion Project focuses on the Orion Crew Module and the Orion Launch Abort Vehicle. NASA TetrUSS codes (GridTool/VGRID/USM3D) have been applied in a supporting role to the Crew Exploration Vehicle Aerosciences Project for investigating various aerodynamic sensitivities and supplementing the aerodynamic database. This paper provides an overview of the contributions from the TetrUSS team to the Project Orion Crew Module and Launch Abort Vehicle aerodynamics, along with selected examples to highlight the challenges encountered along the way. A brief description of geometries and tasks will be discussed followed by a description of the flow solution process that produced production level computational solutions. Four tasks conducted by the USM3D team will be discussed to show how USM3D provided aerodynamic data for inclusion in the Orion aero-database, contributed data for the build-up of aerodynamic uncertainties for the aero-database, and provided insight into the flow features about the Crew Module and the Launch Abort Vehicle.

Mcmillin, Susan N.↗
DOE OSTI2025

Experimental characterization and potential energy savings of insulated cladding for U.S. residential buildings

This study evaluates a novel insulated cladding composed of foamed cement and vinyl siding that is designed as a retrofit product to enhance thermal performance in U.S. residential buildings. Experimental measurements showed that foamed cement exhibits a thermal conductivity of 0.0334–0.0365 W/m-K (R-value of 4.0–4.3 °F⋅ft 2 ⋅h/BTU), which is significantly higher than that of conventional fiber cement boards. Although mechanical testing confirmed the material’s brittleness and low strength, it also indicated its suitability for nonstructural insulation applications. Water vapor permeability testing demonstrated effective moisture resistance with the polymer coating. To evaluate the insulated cladding’s energy performance, we used the U.S. Department of Energy’s prototype single-family building model and simulated across 50 locations spanning 16 International Energy Conservation Code climate zones. The R-7 (R-value of 7 °F⋅ft 2 ⋅h/BTU per inch) retrofit cladding resulted in heating energy savings up to 275 therms (Fairbanks, AK) and cooling energy savings up to 1,547 kWh (Phoenix, AZ). The energy cost savings varied by region, but the highest annual cost savings ($\$$434/year) were observed in Santa Maria, CA, and the highest annual cost savings percentage (37 %) were noted in Monterey, CA. On a national scale, the insulated retrofit provided an average annual energy cost reduction of 16 %, highlighting the cladding system’s broad applicability and financial viability.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Dissipation and particle energization in moderate to low beta turbulent plasma via PIC simulations

Here, we simulate decaying turbulence in electron-positron pair plasmas using a fully-kinetic particle-in-cell (PIC) code. We run two simulations with moderate-to-low plasma β (the ratio of thermal pressure to magnetic pressure). The energy decay rate is found to be similar in both cases. The perpendicular wave-number spectrum of magnetic energy shows a slope between ${k}_{\perp }^{-1.3}$ and ${k}_{\perp }^{-1.1}$, where the perpendicular (⊥) and parallel (∥) directions are defined with respect to the magnetic field. The particle kinetic energy distribution function shows the formation of a non-thermal feature in the case of lower plasma β, with a slope close to E-1. The correlation between thin turbulent current sheets and Ohmic heating by the dot product of electric field (E) and current density (J) is investigated. By heating the parallel E∥ centerdot J∥ term dominates the perpendicular E⊥ centerdot J⊥ term. Regions of strong E∥ centerdot J∥ are spatially well-correlated with regions of intense current sheets, which also appear correlated with regions of strong E∥ in the low β simulation, suggesting an important role of magnetic reconnection in the dissipation of low β plasma turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Evaluating the effectiveness of program data features for guiding memory management

Recent trends have led to the adoption of larger and more complex memory systems, often with multiple tiers of memory performance within the same platform. To utilize complex memory systems efficiently, current data management strategies must be altered to map usage demands to the underlying hardware. Applications, as the generators of memory accesses, are well-suited to guide this process, but building and maintaining separate source code versions for different memory systems is not feasible in most cases. One potential solution is to employ automated program profiling and analysis to facilitate the production of application-based guidance. By attaching memory usage information to static or lightweight program features, compilers and runtime systems can generate fine-grained guidance without additional efforts from users or developers. Recent works have employed this approach with some success, but it is not clear which program features are most useful for guiding data management.This work evaluates the effectiveness of using different program data features to predict memory usage and guide memory management. It employs a custom set of simulation tools, based in the Intel® Pin framework, to collect and analyze the usage characteristics of application data associated with common program data features, such as allocation sites, types, and context. The results show that even relatively simple features, such as object size, are useful for selecting data with similar usage properties, but finer-grained features, such as the instructions that access a particular object, are often much more effective. Additionally, this work evaluates the performance of using different data features and different program inputs to guide data placement on a heterogeneous memory platform with a limited amount of high performance memory.

Effler, T. Chad↗
DOE OSTI2018

Bison Improvements for Robustness and Speed

Early in its development, the primary emphasis of efforts to improve the Bison fuel performance code were focused on initial development of new capabilities. As Bison has matured and become adopted by a wider set of users, it has become clear that a number of robustness issues needed to be addressed to ensure that the code can reliably produce a converged solution. This need became particularly evident when a recent effort by Westinghouse Electric Company to simulate all of the fuel rods in a reactor using Bison resulted in approximately 30% of the models failing to run to completion because of lack of convergence or other errors that were encountered. It is important to emphasize that there was no evidence that Bison generated results that were incorrect for the set of models exercised. In many cases, issues such as those encountered could be addressed by adjusting solver-related parameters. However, there is no question that for Bison to become widely adopted for production work, it needs to be inherently far more robust than it was. Because of this, a concerted effort was undertaken to improve Bison's robustness, or its ability to reliably obtain converged solutions under a wider variety of conditions. As a result of this work, a number of high-priority areas for improvement for modeling light water reactor (LWR) fuel rods have been identified and addressed. As a result of this work, the number of failing Bison analyses in the aforementioned set of simulations for a full core has been reduced to a fraction of a percent. In addition, the robustness of Bison models using the VERA core simulator has been significantly improved. This is important because all of the Bison models representing individual fuel rods must run to completion successfully for the full core model to run to completion. This report provides summaries of the identified sources of robustness issues in Bison and the development that was done to resolve these issues. It also summarizes the addition of new capabilities to address deficiencies in Bison's ability to include needed aspects of LWR fuel behavior, and improvements that were made to improve Bison's solution efficiency.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Simulation of xenon, uranium vacancy and interstitial diffusion and grain boundary segregation in UO 2

In light water reactor fuel, gaseous fission products segregate to grain boundaries, resulting in the nucleation and growth of large intergranular fission gas bubbles. Based on the mechanisms established from density functional theory (DFT) and empirical potential calculations 1, continuum models for diffusion of xenon (Xe), uranium (U) vacancies and U interstitials in UO 2 have been derived for both intrinsic conditions and under irradiation. Segregation of Xe to grain boundaries is described by combining the bulk diffusion model with a model for the interaction between Xe atoms and three different grain boundaries in UO 2 ( Σ5 tilt, Σ5 twist and a high angle random boundary),as derived from atomistic calculations. All models are implemented in the MARMOT phase field code, which is used to calculate effective Xe and U diffusivities as well as redistribution for a few simple microstructures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil" in the Journal of Geophysical Research-Biogeosciences. We simulated biogeochemical cycling in arctic soils using the PFLOTRAN geochemical model combined with measurements from previous NGEE Arctic incubations of polygonal permafrost soils in northern Alaska (Zheng et al., 2018). Simulated iron cycling, carbon dioxide production, and methane production were compared with incubation measurements and the parameterized model was then used to simulate coupled iron and carbon cycling over repeated oxic-anoxic cycles at different levels of carbon substrate availability and pH. The most recent data version (2.0) in the archive incorporates changes to the model and simulations as suggested by reviewers during the manuscript review process. These changes include an updated parameterization of the model; a new set of simulations omitting the iron cycle for direct evaluation of how iron cycle processes affect modeled outcomes; and a set of simulations testing different scenarios of carbon substrate availability in addition to scenarios of initial soil pH. This archive contains simulation code, model output, and analysis code for PFLOTRAN simulations. All scripts are python except the batch script for submitting multiprocessor jobs. Note that the model also requires compiled versions of the Alquimia interface and the NGEE Arctic fork of the PFLOTRAN geochemical simulator (see the README_INSTALL document for basic instructions). The Output directory contains eight data files in netCDF format generated by the model. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 10-year research effort (2012-2022) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1991

Hydrodynamical instabilities and mixing in SN 1987A - Two-dimensional simulations of the first 3 months

Results are presented from numerical simulations of the early stages of the explosion of SN 1987A. Using a two-dimensional cylindrical geometry version of a smooth particle hydrodynamics code, the explosion is followed for three months to investigate both the early hydrodynamical instabilities and the effect of the subsequent radioactive decay of Ni-56 and Co-56 with half-lives of 6.1 and 77.8 days, respectively. It is shown that the mixing induced by hydrodynamical instabilities occurring during the first few hours is substantially modified at later time by the radioactive decay of Ni-56 and Co-56. The inner cavity of the expanding supernova remnant fills up with nickel, its decay products thus forming a giant 'nickel bubble'. The peak velocity of the nickel increases by approximately 30 percent after the decays. While these results adequately model the core of the observed Fe line profiles, they fail to reproduce the high velocity wings of the spectra.

Herant, Marc↗
NASA NTRS2003

The Formation of Asteroid Satellites in Catastrophic Impacts: Results from Numerical Simulations

We have performed new simulations of the formation of asteroid satellites by collisions, using a combination of hydrodynamical and gravitational dynamical codes. This initial work shows that both small satellites and ejected, co-orbiting pairs are produced most favorably by moderate-energy collisions at more direct, rather than oblique, impact angles. Simulations so far seem to be able to produce systems qualitatively similar to known binaries. Asteroid satellites provide vital clues that can help us understand the physics of hypervelocity impacts, the dominant geologic process affecting large main belt asteroids. Moreover, models of satellite formation may provide constraints on the internal structures of asteroids beyond those possible from observations of satellite orbital properties alone. It is probable that most observed main-belt asteroid satellites are by-products of cratering and/or catastrophic disruption events. Several possible formation mechanisms related to collisions have been identified: (i) mutual capture following catastrophic disruption, (ii) rotational fission due to glancing impact and spin-up, and (iii) re-accretion in orbit of ejecta from large, non-catastrophic impacts. Here we present results from a systematic investigation directed toward mapping out the parameter space of the first and third of these three collisional mechanisms.

Durda, D. D.↗
DOE OSTI2021

An efficient numerical model for predicting residual stress and strain in parts manufactured by laser powder bed fusion

Abstract Computational modeling of additively manufactured structures plays an increasingly important role in product design and optimization. For laser powder bed fusion processes, the accurate modeling of stress and distortion requires large amount of computational cost due to very localized heat input and evolving complex geometries. The current study takes advantage of a graphics processing unit accelerated explicit finite element analysis code and approximated heat conduction analysis to predict the macroscopic thermo-mechanical behavior in laser selective melting. Adjacent layers and tracks were lumped to reduce the number of time steps and elements in the finite element model. The effects of track and layer grouping on prediction accuracy and solution efficiency are investigated to provide a guidance for a cost-effective simulation. Thin-wall builds from Inconel alloy 625 (IN625) powders were simulated by applying the developed modeling approach to get the detailed residual stress and distortion at a computational speed 50 times higher than conventional approach. Under repeated heating and cooling cycles, a high tensile stress was produced near surfaces of a build due to a larger shrinkage on surface than that in central area. It is also shown that horizontal stresses concentrate near the root and top layers of the IN625 build. The predicted residual elastic strain distribution was validated by the experimental measurement using x-ray synchrotron diffraction.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Three-dimensional calculation of supersonic reacting flows using an LU scheme

A new three-dimensional numerical program incorporated with comprehensive real gas property models has been developed to simulate supersonic reacting flows. The code employs an implicit finite volume, Lower-Upper (LU) time-marching method to solve the complete Navier-Stokes and species equations in a fully-coupled and very efficient manner. A chemistry model with nine species and eighteen reaction steps are adopted in the program to represent the chemical reaction of H2 and air. To demonstrate the capability of the program, flow fields of underexpanded hydrogen jets transversely injected into supersonic air stream inside the combustors of scramjets are calculated. Results clearly depict the flow characteristics, including the shock structure, separated flow regions around the injector, and the distribution of the combustion products.

Yu, Sheng-Tao↗
NASA NTRS1989

Three-dimensional calculation of supersonic reacting flows using an LU scheme

A new three-dimensional numerical program incorporated with comprehensive real gas property models was developed to simulate supersonic reacting flows. The code employs an implicit finite volume, Lower-Upper (LU) time-marching method to solve the complete Navier-Stokes and species equations in a fully-coupled and very efficient manner. A chemistry model with nine species and eighteen reaction steps are adopted in the program to represent the chemical reaction of H2 and air. To demonstrate the capability of the program, flow fields of underexpanded hydrogen jets transversely injected into supersonic air stream inside the combustors of scramjets are calculated. Results clearly depict the flow characteristics, including the shock structure, separated flow regions around the injector, and the distribution of the combustion products.

Yu, Sheng-Tao↗