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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Kinetic modeling of hot tail runaway electron generation during plasma disruptions using the JOREK code

The generation of runaway electrons (REs) during disruptions poses a significant challenge for the operation of tokamaks. The production of these high-energy electrons can cause substantial damage, particularly when the plasma current is high, making it a critical concern for ITER. For the high-temperature plasmas anticipated in ITER, the primary generation of REs may be dominated by the hot tail mechanism, which consists of the acceleration of hot electrons from the pre-disruption population which have not yet thermalized with the bulk following the rapid cooling of the plasma. To account for the significant 3D effects on RE production, a hot tail modeling framework has been developed within the non-linear 3D extended MHD code JOREK. This paper presents the structure of this framework, which is based on test electrons evolving in MHD fields. The verification of the method shows good agreement with the reference DREAM code for 0D test cases, as well as for axisymmetric simulations of 15 MA ITER H-mode disruption scenarios. Furthermore, a proof-of-principle application to a DIII-D case demonstrates the framework’s capability to capture for the first time the hot tail generation in 3D MHD simulations in realistic geometry. Preliminary results suggest that the production of REs is significantly reduced by stochastic losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Transport model comparison studies of intermediate-energy heavy-ion collisions

Transport models are the main method to obtain physics information on the nuclear equation of state and in-medium properties of particles from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. To this end, calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, which test separately selected ingredients of a transport code, and more realistic calculations of heavy-ion collisions. Over the years, six studies have been performed within this project. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. After a compact description of the underlying transport approaches, we review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known under the names of Boltzmann–Uehling–Uhlenbeck (BUU) and Quantum Molecular Dynamics (QMD) type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. Furthermore, our evaluation studies often indicate improved strategies in performing transport simulations and thus can provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Simulation of a Flexible Production System

The development of computer models for a flexible production system is considered, as a step toward automated design. The formulation of modeling algorithms for flexible production systems at the level of technological operations is outlined. The corresponding software for a computerized modeling system is described.

42 ENGINEERING↗
DOE OSTI2017

Atomic scale simulations for improved CRUD and fuel performance modeling

A more mechanistic description of fuel performance codes can be achieved by deriving models and parameters from atomistic scale simulations rather than fitting models empirically to experimental data. The same argument applies to modeling deposition of corrosion products on fuel rods (CRUD). Here are some results from publications in 2016 carried out using the CASL allocation at LANL.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Onset of Fluidization in MP-PIC Simulations using MFIX-Exa

Fluidized bed reactors are used across a variety of industries, including for energy processes like pyrolysis that result in low-cost energy products. Design and scale-up of fluidized beds is de-risked by modeling and simulation, utilizing tools like NETL’s MFIX-Exa High-Performance Computing (HPC) code for reacting multiphase flow. This report summarizes an investigation into the breadth of problems to which MFIX-Exa may be applied, specifically with regard to low fluid velocities and the onset of fluidization. A simple fluidization study is conducted both experimentally and numerically for particles of interest, then reactor simulations are compared to cold flow experiments for uniform distributor plates. Approaches for modeling bubble caps are also presented.

discrete particle method↗
NASA NTRS1996

Quiet High-Speed Fan

A calibration of the acoustic and aerodynamic prediction methods was performed and a baseline fan definition was established and evaluated to support the quiet high speed fan program. A computational fluid dynamic analysis of the NASA QF-12 Fan rotor, using the DAWES flow simulation program was performed to demonstrate and verify the causes of the relatively poor aerodynamic performance observed during the fan test. In addition, the rotor flowfield characteristics were qualitatively compared to the acoustic measurements to identify the key acoustic characteristics of the flow. The V072 turbofan source noise prediction code was used to generate noise predictions for the TFE731-60 fan at three operating conditions and compared to experimental data. V072 results were also used in the Acoustic Radiation Code to generate far field noise for the TFE731-60 nacelle at three speed points for the blade passage tone. A full 3-D viscous flow simulation of the current production TFE731-60 fan rotor was performed with the DAWES flow analysis program. The DAWES analysis was used to estimate the onset of multiple pure tone noise, based on predictions of inlet shock position as a function of the rotor tip speed. Finally, the TFE731-60 fan rotor wake structure predicted by the DAWES program was used to define a redesigned stator with the leading edge configured to minimize the acoustic effects of rotor wake / stator interaction, without appreciably degrading performance.

Lieber, Lysbeth↗
DOE OSTI2022

Modeling mesoscale fission gas behavior in UO2 by directly coupling the phase field method to spatially resolved cluster dynamics

Abstract Fission gas release within uranium dioxide nuclear fuel occurs as gas atoms diffuse through grains and arrive at grain boundary (GB) bubbles; these GB bubbles grow and interconnect with grain edge bubbles; and grain edge tunnels grow and connect to free surfaces. In this study, a hybrid multi-scale/multi-physics simulation approach is presented to investigate these mechanisms of fission gas release at the mesoscale. In this approach, fission gas production, diffusion, clustering to form intragranular bubbles, and re-solution within grains are included using spatially resolved cluster dynamics in the Xolotl code. GB migration and intergranular bubble growth and coalescence are included using the phase field method in the MARMOT code. This hybrid model couples Xolotl to MARMOT using the MultiApp and Transfer systems in the MOOSE framework, with Xolotl passing the arrival rate of gas atoms at GBs and intergranular bubble surfaces to MARMOT and MARMOT passing evolved GBs and bubble surface positions to Xolotl. The coupled approach performs well on the two-dimensional simulations performed in this work, producing similar results to the standard phase field model when Xolotl does not include fission gas clustering or re-solution. The hybrid model performs well computationally, with a negligible cost of coupling Xolotl and MARMOT and good parallel scalability. The hybrid model predicts that intragranular fission gas clustering and bubble formation results in up to 70% of the fission gas being trapped within grains, causing the increase in the intergranular bubble fraction to slow by a factor of six. Re-solution has a small impact on the fission gas behavior at 1800 K but it has a much larger impact at 1000 K, resulting in a twenty-times increase in the concentration of single gas atoms within grains. Due to the low diffusion rate, this increase in mobile gas atoms only results in a small acceleration in the growth of the intergranular bubble fraction. Finally, the hybrid model accounts for migrating GBs sweeping up gas atoms. This results in faster intergranular bubble growth with smaller initial grain sizes, since the additional GB migration results in more immobile gas clusters reaching GBs.

Kim, Dong-Uk↗
NASA NTRS2012

Submillimetric Spectroscopic Observations of Volatiles in Comet C-2004 Q2 (Machholz)

Submillimeter spectroscopic observations of comets provide an important tool for understanding their chemical composition and enable a taxonomic classification. Aims. We aim to determine the production rates of several parent- and product volatiles and the C-12/C-13 isotopic carbon ratio in the long-period comet C/2004 Q2 (Machholz), which is likely to originate from the Oort Cloud. Methods. The line emission from several molecules in the coma was measured with high signal-to-noise ratio in January 2005 at heliocentric distance of 1.2 AU by means of high-resolution spectroscopic observations using the Submillimeter Telescope (SMT) at the Arizona Radio Observatory (ARO). Results. We have obtained production rates of several volatiles (CH3OH, HCN, H(sup 13)CN, HNC, H2CO, CO, and CS) by comparing the observed and simulated line-integrated intensities. We calculated the synthetic profiles using a radiative transfer code that includes collisions between neutrals and electrons, and the effects of radiative pumping of the fundamental vibrational levels by solar infrared radiation. Furthermore, multiline observations of the CH3OH J = 7-6 series allow us to estimate the rotational temperature using the rotation diagram technique. We find that the CH3OH population distribution of the levels sampled by these lines can be described by a rotational temperature of 40 +/- 3 K. Derived mixing ratios relative to hydrogen cyanide are CO/CH3OH/H2CO/CS/HNC/HC-13N/HCN= 30.9/24.6/4.8/0.57/0.031/0.013/1 assuming a pointing offset of 8" due to the uncertain ephemeris at the time of the observations and the telescope pointing error. Conclusions. The measured relative molecular abundances in C/2004 Q2 (Machholz) are between low- to typical values of those obtained in Oort Cloud comets, suggesting that it has visited the inner solar system previously and undergone thermal processing. The HNC/HCN abundance ratio of approx 3.1% is comparable to that found in other comets, accounting for the dependence on the heliocentric distance, and could possibly be explained by ion-molecule chemical processes in the low-temperature atmosphere. From a tentative HC-13N detection, the measured value of 97 +/- 30 for the HC-12N/HC-13N isotopologue pair is consistent with a telluric value. The outgassing variability observed in the HCN production rates over a period of two hours is consistent with the rotation of the nucleus derived using different observational techniques.

chemical composition↗
DOE OSTI2016

Modeling and Depletion Simulations for a High Flux Isotope Reactor Cycle with a Representative Experiment Loading

The purpose of this report is to document a high-fidelity VESTA/MCNP High Flux Isotope Reactor (HFIR) core model that features a new, representative experiment loading. This model, which represents the current, high-enriched uranium fuel core, will serve as a reference for low-enriched uranium conversion studies, safety-basis calculations, and other research activities. A new experiment loading model was developed to better represent current, typical experiment loadings, in comparison to the experiment loading included in the model for Cycle 400 (operated in 2004). The new experiment loading model for the flux trap target region includes full length 252 Cf production targets, 75 Se production capsules, 63 Ni production capsules, a 188 W production capsule, and various materials irradiation targets. Fully loaded 238 Pu production targets are modeled in eleven vertical experiment facilities located in the beryllium reflector. Other changes compared to the Cycle 400 model are the high-fidelity modeling of the fuel element side plates and the material composition of the control elements. Results obtained from the depletion simulations with the new model are presented, with a focus on time-dependent isotopic composition of irradiated fuel and single cycle isotope production metrics.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

A Rule-Based Expert System Applied to Moisture Durability of Building Envelopes

The moisture durability of an envelope component such as a wall or roof is difficult to predict. Moisture durability depends on all the construction materials used, as well as the climate, orientation, air tightness, and indoor conditions. Modern building codes require more insulation and tighter construction but provide little guidance about how to ensure these energy-efficient assemblies remain moisture durable. Furthermore, as new products and materials are introduced, builders are increasingly uncertain about the long-term durability of their building envelope designs. Oak Ridge National Laboratory and the US Department of Energy's Building America Program are applying a rule-based expert system methodology in a web tool to help designers determine whether a given wall design is likely to be moisture durable and provide expert guidance on moisture risk management specific to a wall design and climate. The expert system is populated with knowledge from both expert judgment and probabilistic hygrothermal simulation results.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2015

Study of the characteristics of reconfigurable plasma antenna array

This paper presents a design and simulation of a reconfigurable array of plasma antenna. The plasma column is used as radiating elements instead of metal to create an antenna. The advantages of the plasma antenna over the conventional antenna are its possible to change the operating parameters, such as the working pressure, input power, radius of the discharge tube, resonant frequency, and length of the plasma column. In addition, plasma antenna can be reconfigurable with respect to shape, frequency and radiation parameters in a very short time. The plasma discharge tube was designed with a length of 200 mm, the radius of the plasma column was 2.5 mm and the coupling sleeve was connected to the SMA as the ground. This simulation was performed by using the simulation software Computer Simulation Technology (CST). The frequency is set in the range of 1 GHz to 10 GHz. The performance of the designed antenna was analyzed in term of return loss, gain and radiation pattern. For reconfigurable plasma antenna array, it shows that the gain is increase when the number of antenna element is increase. The combination of the discharge tube and metal rod as an antenna array has been done, and the result shows that an array with the plasma element can achieve higher gain.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1990

Analysis of generic scramjet external nozzle flowfields employing simulant gases

Experimental and computational results are presented for a generic scramjet nozzle/afterbody model. Mixtures of argon and Freon were employed as simulant scramjet exhaust gases because they correctly model the inviscid simulation parameters of actual scramjet combustion products, allowing test and analysis at significantly reduced temperatures. Air was also employed as an exhaust gas to provide experimental and computational comparisons with simulant gas effects. Several computational aerodynamic analysis codes, employing different levels of theoretical modeling, were applied to compute solutions about the scramjet geometry. The solutions are compared with measured surface static and flowfield pitot pressure data taken in a Mach 6 freestream.

Tatum, Kenneth E.↗
NASA NTRS2011

Lunar Solar Origins Exploration (LunaSOX)

The Moon offers a unique vantage point from which to investigate the Sun and its interaction via the solar wind magnetic fields, plasma, and energetic particles with the geospace system including the Moon itself. The lunar surface and exosphere provide in part a record of solar coronal plasma material input and resultant space weathering over billions of years. The structure and dynamics of solar wind interactions with the Moon provide an accessible near-Earth laboratory environment for study of general solar wind interactions with the vast multitude of airless asteroidal bodies of the inner solar system. Spacecraft in lunar orbit have the often simultaneous opportunity, except when in the Earth's magnetosphere, to make in-situ compositional measurements of the solar wind plasma and to carry out remote observations from the Moon of the solar corona, potentially enabled by lunar limb occultation of the solar disk. The LunaSOX project at NASA Goddard Space Flight Center is addressing these heliophysical science objectives from and of the Moon with support from NASA's Lunar Advanced Science and Exploration Research (LASER) program: (1) specify history of solar wind parameters at and sunward of the Moon through enhanced access (http://lunasox.gsfc.nasa.gov/) to legacy and operational mission data products from the Apollo era to the present, (2) model field and plasma interactions with the lunar surface, exosphere, and wake, as constrained by the available data, through hybrid kinetic code simulations, and (3) advance mission concepts for heliophysics from and of the Moon.

Cooper, John F.↗
DOE OSTI2015

Approach for the isotopic characterization of irradiated start-up sources for disposal at WIPP - 15674

The National Nuclear Security Administration, Global Threat Reduction Initiative (GTRI)/ Off-site Source Recovery Program (OSRP) requires an approved final disposition pathway for sealed irradiated startup sources with uncertain irradiation histories. To create this path a reliable method is required to determine yields of transmutation, fission and activation products contained within the source at the time of disposal. Because of the uncertain histories of these sources, a number of likely histories are modeled with the code ORIGEN in order to establish the bounds of the radioisotope yields. The ORIGEN simulations are confirmed with gamma spectroscopy, specifically examining the activities of several fission and transmutation isotopes. These calculations, coupled with other previously established characterization techniques for transuranic neutron sources, will be used to establish a waste profile that is in compliance with Waste Isolation Pilot Plant (WIPP) waste acceptance criteria. (authors)

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

Hardening of DOE UQ Codes for Ultra-Deep Proactive Geosteering in Oilfield Service

The primary objective of this project was to harden the DOE uncertainty quantification (UQ) code to solve the inverse problem for the geosteering service in the oilfield drilling and exploration. The primary focus is to verify the feasibility of using QUESO to solve geosteering inverse problems through statistical inversion. QUESO is a software library developed by QUEST, a DOE's Scientific Discovery through Advanced Computing (SciDAC) institute. In the process of drilling a borehole, geosteering is the act of adjusting the well trajectory on the fly to reach the geological targets or to keep the wellbore within desired formation. Geosteering inversion is to infer underground formation structure using measurements such as azimuthal resistivity logging data. The geosteering process relies on the real-time reconstruction of the subsurface earth model, and a robust and efficient inversion algorithm is always necessary to satisfy the job requirement. The new generation of azimuthal resistivity geosteering tool is capable of detecting multiple formation boundaries 20-100 ft away, and the corresponding geosteering inversion problem can be very challenging as many parameters need to be inverted within such a large depth of investigation. The success of the project will support US energy security by helping to advance subsurface modeling and characterization technologies in oilfield applications. An effective subsurface inversion algorithm and the implementation as a software package will promote the usage of ultra-deep azimuthal resistivity measurements, which can help well placement in complicated underground environments, reduce costs and risks of drilling, and increase oil and gas production. It will help medium and small size independent service companies and tool manufactures to enter the more advanced and lucrative directional drilling area. In addition, by hardening the QUESO code and building a transparent interface for easy manipulation of the QUESO functions, we provide a user-friendly tool for a larger scientific and industry communities that rely heavily on modeling and simulation tools for forecasting, parameter studies, design, and decision making. These are all activities that can strongly benefit from meaningful representations of modeling uncertainty supported by our product.

58 GEOSCIENCES↗
DOE OSTI2020

Advanced Simulation and Experiments of Strongly Coupled Geomechanics and Flow for Gas Hydrate Deposits: Validation and Field Application. Final Report

The objectives of the proposed research are (1) to investigate geomechanical responses induced by depressurization experimentally and numerically; (2) to enhance the current numerical simulation technology in order to simulate complex physically coupled processes by depressurization and (3) to perform in-depth numerical analyses of two selected potential production test sites: one based on the deposits observed at the Ulleung basin UBGH2-6 site; and the other based on well-characterized accumulations from the westend Prudhoe Bay. To these ends, the recipient will have the following specific objectives: 1). Information obtained from multi-scale experiments previously conducted at the recipient’s research partner (the Korean Institute of Geoscience and Mineral Resources (KIGAM)) that were designed to represent the most promising known Ulleung Basin gas hydrate deposit as drilled at site UBGH2-6 will be evaluated (Task 2). These findings will be further tested by new experimental studies at Lawrence Berkeley National Laboratory (LBNL) and Texas A&M (TAMU) (Task 3) that are designed capture complex coupled physical processes between flow and geomechanics, such as sand production, capillarity, and formation of secondary hydrates. The findings of Tasks 2 and 3 will be used to further improve numerical codes. 2) Develop (in Tasks 4 through 6) an advanced coupled geomechanics and non-isothermal flow simulator (T+MAM) to account for large deformation and strong capillarity. This new code will be validated using data from the literature, from previous work by the project team, and with the results of the proposed experimental studies. The developed simulator will be applied to both Ulleung Basin and Prudhoe Bay sites, effectively addressing complex geomechanical and petrophysical changes induced by depressurization (e.g., frost-heave, strong capillarity, cryo-suction, induced fracturing, and dynamic permeability).

02 PETROLEUM↗
NASA NTRS1981

Simulation of charge exchange plasma propagation near an ion thruster propelled spacecraft

A model describing the charge exchange plasma and its propagation is discussed, along with a computer code based on the model. The geometry of an idealized spacecraft having an ion thruster is outlined, with attention given to the assumptions used in modeling the ion beam. Also presented is the distribution function describing charge exchange production. The barometric equation is used in relating the variation in plasma potential to the variation in plasma density. The numerical methods and approximations employed in the calculations are discussed, and comparisons are made between the computer simulation and experimental data. An analytical solution of a simple configuration is also used in verifying the model.

Robinson, R. S.↗
NASA NTRS2000

Performance of the OVERFLOW-MLP and LAURA-MLP CFD Codes on the NASA Ames 512 CPU Origin System

The shared memory Multi-Level Parallelism (MLP) technique, developed last year at NASA Ames has been very successful in dramatically improving the performance of important NASA CFD codes. This new and very simple parallel programming technique was first inserted into the OVERFLOW production CFD code in FY 1998. The OVERFLOW-MLP code's parallel performance scaled linearly to 256 CPUs on the NASA Ames 256 CPU Origin 2000 system (steger). Overall performance exceeded 20.1 GFLOP/s, or about 4.5x the performance of a dedicated 16 CPU C90 system. All of this was achieved without any major modification to the original vector based code. The OVERFLOW-MLP code is now in production on the inhouse Origin systems as well as being used offsite at commercial aerospace companies. Partially as a result of this work, NASA Ames has purchased a new 512 CPU Origin 2000 system to further test the limits of parallel performance for NASA codes of interest. This paper presents the performance obtained from the latest optimization efforts on this machine for the LAURA-MLP and OVERFLOW-MLP codes. The Langley Aerothermodynamics Upwind Relaxation Algorithm (LAURA) code is a key simulation tool in the development of the next generation shuttle, interplanetary reentry vehicles, and nearly all "X" plane development. This code sustains about 4-5 GFLOP/s on a dedicated 16 CPU C90. At this rate, expected workloads would require over 100 C90 CPU years of computing over the next few calendar years. It is not feasible to expect that this would be affordable or available to the user community. Dramatic performance gains on cheaper systems are needed. This code is expected to be perhaps the largest consumer of NASA Ames compute cycles per run in the coming year.The OVERFLOW CFD code is extensively used in the government and commercial aerospace communities to evaluate new aircraft designs. It is one of the largest consumers of NASA supercomputing cycles and large simulations of highly resolved full aircraft are routinely undertaken. Typical large problems might require 100s of Cray C90 CPU hours to complete. The dramatic performance gains with the 256 CPU steger system are exciting. Obtaining results in hours instead of months is revolutionizing the way in which aircraft manufacturers are looking at future aircraft simulation work. Figure 2 below is a current state of the art plot of OVERFLOW-MLP performance on the 512 CPU Lomax system. As can be seen, the chart indicates that OVERFLOW-MLP continues to scale linearly with CPU count up to 512 CPUs on a large 35 million point full aircraft RANS simulation. At this point performance is such that a fully converged simulation of 2500 time steps is completed in less than 2 hours of elapsed time. Further work over the next few weeks will improve the performance of this code even further.The LAURA code has been converted to the MLP format as well. This code is currently being optimized for the 512 CPU system. Performance statistics indicate that the goal of 100 GFLOP/s will be achieved by year's end. This amounts to 20x the 16 CPU C90 result and strongly demonstrates the viability of the new parallel systems rapidly solving very large simulations in a production environment.

Taft, James R.↗