Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core simulation”

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 1,189 records · Page 66

Light Water Reactor Sustainability Program Support and Modeling for the Boiling Water Reactor Station Black Out Case Study Using RELAP and RAVEN

The existing fleet of nuclear power plants is in the process of extending its lifetime and increasing the power generated. In order to evaluate the impact of these two factors on the safety of the plant, the Risk Informed Safety Margin Characterization (RISMC) project aims to provide insight to decision makers through a series of simulations of the plant dynamics for different initial conditions (e.g., probabilistic analysis and uncertainty quantification). This report focuses, in particular, on the impact of power uprate on the safety of a boiled water reactor system. The case study considered is a loss of off-site power followed by the loss of diesel generators, i.e., a station black out (SBO) event. Analysis is performed by using a thermo-hydraulic code, i.e. RELAP-5, and a stochastic analysis tool currently under development at INL, i.e. RAVEN. Starting from the event tree models contained in SAPHIRE, we built the input file for RELAP-5 that models in great detail system dynamics under SBO conditions. We also interfaced RAVEN with RELAP-5 so that it would be possible to run multiple RELAP-5 simulation runs by changing specific keywords of the input file. We both employed classical statistical tools, i.e. Monte-Carlo, and more advanced machine learning based algorithms to perform uncertainty quantification in order to quantify changes in system performance and limitations as a consequence of power uprate. We also employed advanced data analysis and visualization tools that helped us to correlate simulation outcome such as maximum core temperature with a set of input uncertain parameters. Results obtained gave a detailed overview of the issues associated to power uprate for a SBO accident scenario. We were able to quantify how timing of safety related events were impacted by a higher reactor core power. Such insights can provide useful material to the decision makers to perform risk-infomed safety margins management.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Quasilinear Carbon Transport In An Impurity Hole Plasma In LHD

Comprehensive electrostatic gyrokinetic linear stability calculations for ion-scale microinstabilities in an LHD plasma with an ion-ITB and carbon "impurity hole" are used to make quasilinear estimates of particle flux to explore whether microturbulence can explain the observed outward carbon fluxes that flow "up" the impurity density gradient. The ion temperature is not stationary in the ion-ITB phase of the simulated discharge, during which the core carbon density decreases continuously. To fully sample these varying conditions the calculations are carried out at three radial locations and four times. The plasma parameter inputs are based on experimentally measured profiles of electron and ion temperature, as well as electron and carbon density. The spectroscopic line-average ratio of hydrogen and helium densities is used to set the density of these species. Three ion species (H,He,C) and the electrons are treated kinetically, including collisions. Electron instability drive does enhance the growth rate significantly, but the most unstable modes have characteristics of ion temperature gradient (ITG) modes in all cases. As the carbon density gradient is scanned between the measured value and zero, the quasilinear carbon flux is invariably inward when the carbon density profile is hollow, so turbulent transport due to the instabilities considered here does not explain the observed outward flux of impurities in impurity hole plasmas. The stiffness of the quasilinear ion heat flux is found to be 1.7-2.3, which is lower than several estimates in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pulsed Magnetic Welding for Advanced Core and Cladding Steel

To investigate a solid-state joining method, pulsed magnetic welding (PMW), for welding the advanced core and cladding steels to be used in Generation IV systems, with a specific application for fuel pin end-plug welding. As another alternative solid state welding technique, pulsed magnetic welding (PMW) has not been extensively explored on the advanced steels. The resultant weld can be free from microstructure defects (pores, non-metallic inclusions, segregation of alloying elements). More specifically, the following objectives are to be achieved: 1. To design a suitable welding apparatus fixture, and optimize welding parameters for repeatable and acceptable joining of the fuel pin end-plug. The welding will be evaluated using tensile tests for lap joint weldments and helium leak tests for the fuel pin end-plug; 2 Investigate the microstructural and mechanical properties changes in PMW weldments of proposed advanced core and cladding alloys; 3. Simulate the irradiation effects on the PWM weldments using ion irradiation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Enhancements to the Subplane Capability in MPACT

In recent years, much of the development effort for the Virtual Environment for Reactor Applications (VERA) has been devoted to improving performance and targeting industry-class computers. This includes reducing both runtime (improving efficiency) and core-count requirements to enable simulation on smaller machines. In MPACT, in particular, many aspects of the neutronic iteration strategy have been evaluated and improved, such as the subgroup and method of characteristics (MOC) runtimes, and many efforts have focused on improving the coarse mesh finite difference (CMFD) runtime, which has become a dominant component of the runtime as other aspects have been improved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Summary Report for U.S. and Japan Bilateral Collaborations (FY2019)

Within the frame work of the Civil Nuclear Energy R&D Working Group (CNWG), researchers from Argonne National Laboratory (ANL), Japan Atomic Energy Agency (JAEA) and Central Research Institute of Electric Power Industry (CRIEPI) have been participating in the R&D collaborations in three tasks that are related to the fast reactor technologies, which include metal fuel sodium-cooled fast reactor (SFR) core bowing analysis, modeling and simulation of Monju Turbine Trip Benchmark, and collaborative improvement of the SAS4A code.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

AI-Enabled Discovery and Physics-Based Optimization of Energy Efficient Processing Strategies for Advanced Turbine Alloys (Final Technical Report)

In this project, the multi-organizational team of academic and industrial researchers from the University of Kentucky an aerospace and energy generation OEM partner has leveraged novel Digital Process Twin (DPT) models of process/structure interactions (i.e., process-induced surface integrity) to advance a paradigm of fully integrated computational materials engineering (ICME). Using efficient process models as the core of a digital process simulator for a reinforcement learning algorithm, the team has integrated industrial data and metrics of structure/performance/energy relationships and manufacturing-related energy metrics to optimize dynamic processing parameters for significantly improved life-cycle energy efficiency of advanced γ-TiAl low-pressure turbine (LPT) alloys, as indicated by a set of design relevant parameters (e.g., residual stresses and scrap rate). The key objective and anticipated outcome of the project was at least a 10% reduction in life-cycle embodied energy for a recently developed, γ-TiAl low-pressure turbine (LPT) alloy and nickel-based superalloy Inconel 718, through the adoption of the proposed AI-enabled process optimization approach. The final project outcomes significantly exceeded this original target, realizing manufacturing-related energy efficiency improvements of more than 130% for TiAl and up to 80% for Inconel 718. Rather than following the prevailing and highly inefficient empirical paradigm, the proposed study demonstrated the feasibility of adopting a digital, physics-based process design and optimization paradigm. The recurring need for manual intervention, rework, reinspection causes significant production bottlenecks and unnecessary expense associated with delivering the requisite component quality. The OEM partner, and turbine industry in general, expect to reap significant cost and resource savings if an AI-optimized set of parameters can be applied to specific machining operations. The technical scope of the proposed project involved the paving of a realistic path towards model-based and AI-enabled Integrated Computational Materials Engineering (ICME), and away from inefficient empirical process optimization and legacy manufacturing practices, which are no longer able to efficiently process novel high-performance turbine alloy materials. The project team will address the fundamental knowledge gap that currently exists within the ICME paradigm with respect to the process/structure/performance/energy impacts of finishing processes. While significant resources have been devoted to the ‘early stages’ of manufacturing, such as alloy design, primary and secondary processing, finishing processes have not been adequately integrated within ICME. To provide an actionable path towards model-based finishing process design (e.g., machining, burnishing, grinding, polishing), we will employ a novel AI-enabled process optimization paradigm, based on a computationally efficient, physics-based process simulator. Through limited experimental work to calibrate and validate our process simulator model via an advanced in-situ characterization technique and process optimization via reinforcement learning, the project will seek to demonstrate a viable alternative to the inefficient ‘legacy’ processing strategies, empirical testing and broad scope machining learning approaches, all of which fail to adequately consider complex process physics. The project team has identified an intermetallic γ-TiAl LPT alloy, which is currently being used as part of the OEM partner’s advanced gas turbine designs. This particular alloy poses significant manufacturing challenges during finishing operations, which limit the degree to which the current turbine design can be manufactured in an energy- and cost-efficient manner. Empirical testing and numerical modeling efforts to optimize processing parameters for γ-TiAl have not been able to resolve these manufacturing challenges, so the proposed physics-based AI-enabled optimization technology would offer a truly novel and transformative capability. The multi-organizational team of academic and industry experts from the UKY and the OEM partner will work together closely to demonstrate the analytical and experimental critical function and characteristic proof of concept of this novel approach.

20 FOSSIL-FUELED POWER PLANTS↗

New Formation Models for the Kepler-36 System

Formation of the planets in the Kepler-36 system is modeled here by detailed numerical simulations according to the core-nucleated accretion scenario. The standard model is updated to include the dissolution of accreting rocky planetesimals in the gaseous envelope of the planet, leading to substantial enrichment of the envelope mass in heavy elements and a non-uniform composition with depth. For Kepler-36 c, models involving in situ formation and models involving orbital migration are considered. The results are compared with standard formation models. The calculations include the formation (accretion) phase as well as the subsequent cooling phase, up to the age of Kepler-36 (7 Gyr). During the latter phase, mass loss induced by stellar XUV radiation is included. In all cases, the results fit the measured mass, 7.84 M ⊕ , and radius, 3.68 R ⊕ , of Kepler-36 c. Two parameters are varied to obtain these fits: the disk solid surface density at the formation location and the "efficiency" factor in the XUV mass-loss rate. The updated models are hotter and therefore less dense in the silicate portion of the planet and in the overlying layers of H/He, as compared with standard models. The lower densities mean that only about half as much H/He is needed to be accreted to fit the present-day mass and radius constraints. For Kepler-36 b, an updated in situ calculation shows that the entire H/He envelope is lost, early in the cooling phase, in agreement with observation.

79 ASTRONOMY AND ASTROPHYSICS↗

Grand-design Spiral Arms in a Young Forming Circumstellar Disk

We study formation and long-term evolution of a circumstellar disk in a collapsing molecular cloud core using a resistive magnetohydrodynamic simulation. While the formed circumstellar disk is initially small, it grows as accretion continues, and its radius becomes as large as 200 au toward the end of the Class-I phase. A pair of grand-design spiral arms form due to gravitational instability in the disk, and they transfer angular momentum in the highly resistive disk. Although the spiral arms disappear in a few rotations as expected in a classical theory, new spiral arms form recurrently as the disk, soon becoming unstable again by gas accretion. Such recurrent spiral arms persist throughout the Class-0 and I phases. We then perform synthetic observations and compare our model with a recent high-resolution observation of a young stellar object Elias 2–27, whose circumstellar disk has grand-design spiral arms. We find good agreement between our theoretical model and the observation. Our model suggests that the grand-design spiral arms around Elias 2–27 are consistent with material arms formed by gravitational instability. If such spiral arms commonly exist in young circumstellar disks, it implies that young circumstellar disks are considerably massive and gravitational instability is the key process of angular momentum transport.

79 ASTRONOMY AND ASTROPHYSICS↗

Formation and Evolution Simulations of Saturn, Including Composition Gradients and Helium Immiscibility

In this paper, the formation of Saturn is modeled by detailed numerical simulations according to the core-nucleated accretion scenario. Previous models are enhanced to include the dissolution of accreting planetesimals, composed of water ice, rock, and iron, in the gaseous envelope of the planet, leading to a nonuniform composition with depth. The immiscibility of helium in metallic hydrogen layers is also considered. The calculations start at a mass of 0.5 Earth masses and are extended to the present day. At 4.57 Gyr, the model, proceeding outward, has the following structure: (i) a central core composed of 100% heavy elements and molecules, (ii) a region with a decreasing heavy-element mass fraction, down to a value of 0.1, (iii) a layer of uniform composition with the helium mass fraction Y enhanced over the primordial value, (iv) a helium rain region with a gradient in Y, (v) an outer convective, adiabatic region with uniform composition in which Y is reduced from the primordial value, and (vi) the very outer layers where cloud condensation of the heavy elements occurs. Models of the distribution of heavy elements as a function of radius are compared with those derived to fit the observations of the Cassini mission, with rough qualitative agreement. The helium mass fraction in Saturn’s outer layers is estimated to be around 20%. Models are found that provide good agreement with Saturn’s intrinsic luminosity and radius.

79 ASTRONOMY AND ASTROPHYSICS↗

N2O and CO production by electric discharge - Atmospheric implications

Enhanced levels of N2O and CO were measured in tropospheric air samples exposed to a 17,500-J laboratory discharge. These enhanced levels correspond to an N2O production rate of about 4 trillion molecules/J and a CO production rate of about 10 to the 14th molecules/J. The CO measurements suggest that the primary region of chemical production in the discharge is the shocked air surrounding the lightning channel, as opposed to the slower-cooling inner core. Additional experiments in a simulated Venus atmosphere (CO2 - 95%, N2 - 5%, at one atmosphere) indicate an enhancement of CO from less than 0.1 ppm prior to the laboratory discharge to more than 2000 ppm after the discharge. Comparison with theoretical calculations appears to confirm the ability of a shock-wave/thermochemical model to predict the rate of production of trace species by an electrical discharge.

Levine, J. S.↗

NTP system simulation and detailed nuclear engine modeling

The topics are presented in viewgraph form and include the following: nuclear thermal propulsion (NTP) & detailed nuclear engine modeling; modeling and engineering simulation of nuclear thermal rocket systems; nuclear thermal rocket simulation system; INSPI-NTVR core axial flow profiles; INSPI-NTRV core axial flow profiles; specific impulse vs. chamber pressure; turbine pressure ratio vs. chamber pressure; NERVA core axial flow profiles; P&W XNR2000 core axial flow profiles; pump pressure rise vs. chamber pressure; streamline of jet-induced flow in cylindrical chamber; flow pattern of a jet-induced flow in a chamber; and radiative heat transfer models.

Anghaie, Samim↗

Ionospheric Simulation System for Satellite Observations and Global Assimilative Modeling Experiments (ISOGAME)

ISOGAME is designed and developed to assess quantitatively the impact of new observation systems on the capability of imaging and modeling the ionosphere. With ISOGAME, one can perform observation system simulation experiments (OSSEs). A typical OSSE using ISOGAME would involve: (1) simulating various ionospheric conditions on global scales; (2) simulating ionospheric measurements made from a constellation of low-Earth-orbiters (LEOs), particularly Global Navigation Satellite System (GNSS) radio occultation data, and from ground-based global GNSS networks; (3) conducting ionospheric data assimilation experiments with the Global Assimilative Ionospheric Model (GAIM); and (4) analyzing modeling results with visualization tools. ISOGAME can provide quantitative assessment of the accuracy of assimilative modeling with the interested observation system. Other observation systems besides those based on GNSS are also possible to analyze. The system is composed of a suite of software that combines the GAIM, including a 4D first-principles ionospheric model and data assimilation modules, an Internal Reference Ionosphere (IRI) model that has been developed by international ionospheric research communities, observation simulator, visualization software, and orbit design, simulation, and optimization software. The core GAIM model used in ISOGAME is based on the GAIM++ code (written in C++) that includes a new high-fidelity geomagnetic field representation (multi-dipole). New visualization tools and analysis algorithms for the OSSEs are now part of ISOGAME.

Pi, Xiaoqing↗

New Formation Models for the Kepler-36 System

Formation of the planets in the Kepler-36 system is modeled by detailed numerical simulations according to the core-nucleated accretion scenario. The standard model is updated to include the dissolution of accreting rocky planetesimals in the gaseous envelope of the planet, leading to substantial enrichment of the envelope mass in heavy elements and a non-uniform composition with depth. For Kepler-36 c, models involving in situ formation and models involving orbital migration are considered. The results are compared with standard formation models. The calculations include the formation (accretion) phase, as well as the subsequent cooling phase, up to the age of Kepler-36 (7 Gyr). During the latter phase, mass loss induced by stellar XUV radiation is included. In all cases, the results fit the measured mass, 7.84 M, and radius, 3.68 R, of Kepler-36 c. Two parameters are varied to obtain these fits: the disk solid surface density at the formation location, and the "efficiency" factor in the XUV mass loss rate. The updated models are hotter and therefore less dense in the silicate portion of the planet and in the overlying layers of H/He, as compared with standard models. The lower densities mean that only about half as much H/He is needed to be accreted to fit the present-day mass and radius constraints. For Kepler-36 b, an updated in situ calculation shows that the entire H/He envelope is lost, early in the cooling phase, in agreement with observation.

planets and satellites:formation↗

Design, Instrumentation, and Data Analysis for the SLS Core Stage Green Run Test Series

The Space Launch System (SLS) Core Stage (CS) Thrust Vector Control (TVC) system is comprised of eight mechanical feedback Shuttle heritage Type III TVC actuators and four RS-25 engines, each attached to a Shuttle heritage gimbal block/bearing. The actuators are powered by a Shuttle-derived hydraulic Core Auxiliary Power Unit (CAPU), and integrated with an all-new Core Stage thrust structure. The actuators are interfaced to the SLS Vehicle Management (VM) software via an all-new TVC Actuator Control (TAC) avionics subsystem. Despite the significant test and flight experience of the Shuttle hardware, the SLS Green Run ambient and hot fire test activities revealed a number of new findings associated with the dynamic response of the TVC integrated system. Test responses suggested that the TVC system did not meet its performance specifications and its step and frequency responses exhibited unexpected departures from prior lab tests and modeled behavior. This paper is the fifth installment in a seven-paper series surveying the design, engineering, test validation, and flight performance of the Core Stage Thrust Vector Control system. In this paper, the design of the TVC analyses conducted during the Core Stage Green Run test series are discussed in detail. Throughout the course of the test activities, the SLS flight control team worked diligently with the Core Stage contractor to revise test command profiles and ensure sufficient instrumentation was available to collect data. Post-test analysis combined the Green Run modal, ambient, and hot fire test data, MSFC 2- axis Core Stage TVC Inertial Load Simulator (ILS) data, Hardware-In-the-Loop (HWIL) Systems Integration Lab (SIL) results, and actuator Acceptance Testing Procedure (ATP) responses. These data were used to characterize the response, validate critical math models of the TVC subsystem, and isolate the probable cause of the unexpected responses. Through comprehensive analysis of the available test data sources, the integrated team identified the dominant contributors to the observed response and developed test-correlated rationale for vehicle flight control system performance, ultimately leading to a confident posture for the Artemis I mission.

John H. Wall↗

Efficiency of static core turn-off in a system-on-a-chip with variation

A processor-implemented method for improving efficiency of a static core turn-off in a multi-core processor with variation, the method comprising: conducting via a simulation a turn-off analysis of the multi-core processor at the multi-core processor's design stage, wherein the turn-off analysis of the multi-core processor at the multi-core processor's design stage includes a first output corresponding to a first multi-core processor core to turn off; conducting a turn-off analysis of the multi-core processor at the multi-core processor's testing stage, wherein the turn-off analysis of the multi-core processor at the multi-core processor's testing stage includes a second output corresponding to a second multi-core processor core to turn off; comparing the first output and the second output to determine if the first output is referring to the same core to turn off as the second output; outputting a third output corresponding to the first multi-core processor core if the first output and the second output are both referring to the same core to turn off.

42 ENGINEERING↗

Bellerophon

The Bellerophon software system supports CHIMERA, a production-level HPC application that simulates the evolution of core-collapse supernovae. Bellerophon enables CHIMERA's geographically dispersed team of collaborators to perform job monitoring and real-time data analysis from multiple supercomputing resources, including platforms at OLCF, NERSC, and NICS. Its multi-tier architecture provides an encapsulated, end-to-end software solution that enables the CHIMERA team to quickly and easily access highly customizable animated and static views of results from anywhere in the world via a cross-platform desktop application.

Lingerfelt, Eric J↗

Solution of the Hot Zero Power BEAVRS Using a Pin-Cell Homogenized Model

The validation of a reactor physics code system is an essential part of the code development process. It allows for the evaluation of the code's ability to predict and model the complex physical processes that take place in a nuclear reactor under various operational conditions. Idaho National Laboratory (INL) has been developing a set of reactor physics analysis tools within the MOOSE framework. MAMMOTH is the main reactor physics application, which integrates other MOOSE applications to solve complex multiphysics problems. Results from this study demonstrate that the MOOSE based reactor physics tools, currently under development at INL, have the flexibility to model with high accuracy the HZP condition for a standard PWR. The DRAGON5-Rattlesnake full core pin-cell homogenized simulation produces eigenvalues and detector signals that are comparable to those generated with OpenMC and nTRACER. However, many improvements are necessary to achieve higher fidelity results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computational Study of Droplet Dynamics in Dispersed Flow Film Boiling Regime in PWR Geometries

Dispersed flow film boiling (DFFB) regime is exhibited during low-flooding rate conditions following the postulated Loss-of-Coolant Accident (LOCA) in nuclear reactor flow channels. LOCA is a type of design-basis accident scenario, which is crucial for reactor design and licensing. Accurate models to predict the heat transfer in DFFB regime is critical for the reactor to maintain safety limits. The nominal operating linear heat rate is also set based on the maximum linear heat rate determined for a LOCA transient. It is understood that the behavior of entrained droplets in DFFB regime greatly affects the heat transfer by evaporating and desuperheating the surrounding steam. Furthermore, interaction with flow structures also aids in increasing the evaporation rate of entrained droplets. In order to quantify this enhancement in heat transfer, the size of droplets downstream of flow structures (Spacer Grids and Mixing Vanes) is of crucial importance. Thermal hydraulic analysis codes like COBRA-TF, CATHARE-3 and adopt a three-field equation model that allows them to model the droplet field separately. However, to ensure the accuracy of these field equations, various correlations and physical models are required. This paper proposes a method for characterizing the mean Sauter Diameter of droplets across spacer-grids and mixing vanes in PWR sub-channels using PHASTA - a high fidelity flow solver with multiphase capabilities using the level-set method. The PHASTA code is a finite element based computational fluid dynamics code designed for solving compressible and incompressible turbulent flows. It has been demonstrated to be massively scalable (up to 786,432 cores), which is crucial for simulating large realistic domains. The code incorporates the level-set method for interface tracking of the gas-liquid interface.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗