SEARCH · Engineering Papers

Results for “3D temperature modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Ozone Measurements and a 3D Chemical Transport Model

We have used our three-dimensional chemical transport model (CTM) to calculate the expected reponse of stratospheric composition over the past 30 years to forcing by chlorine and bromine compounds, solar ultraviolet, and volcanic aerosols. The CTM uses off-line winds and temperatures fiom a 50-year run of the finite volume general circulation model (FVGCM). We compare the total column ozone and the ozone profile fiom the CTM output to a variety of data sources. These include a merged total ozone data set from TOMS and SBUV using the new version 8 algorithm. Total ozone fiom the CTM are compared to ground-station measurements of total ozone at specific locations. Ozone profiles are compared to satellite meausrements fiom SBUV, SAGE, and HALOE. Profiles are also compared to ozonesondes over several locations. The results of the comparisons are quantified by using a time-series statistical analysis to determine trends, solar cycle, and volcanic reponse in both the model and in the data. Initial results indicate that the model responds to forcings in a way that is similar to the observed atmospheric response. The model does seem to be more sensitive to the chlorine and bromine perturbation ihan is the data. Further details and comparisons wiii be discussed.

Stolarski, Richard S.↗

DOE OSTI2023

RELAP5-3D Solutions to Exercise 1 of the OECD-NEA HTTF Benchmark

The high-temperature gas-cooled reactor (HTGR) is an advanced reactor concept that has received considerable attention over the last 60 years. HTGRs boast high coolant outlet temperatures, passive safety, and large margins to fuel failure. Recent years have seen increasing interest in HTGRs, whether large-scale HTGRs for electricity or process heat applications or as microreactors serving remote communities. HTGRs have considerable potential to fulfill these applications, but modeling and simulation tools typically used for reactor safety analysis have not been validated for HTGR modeling. The High-Temperature Test Facility (HTTF) was constructed at Oregon State University (OSU) to provide validation data for HTGR modeling. Recently, Idaho National Laboratory (INL), OSU, Argonne National Laboratory, Canadian Nuclear Laboratories, and the University of Tennessee Knoxville have collaborated to develop an HTGR thermal hydraulics benchmark based on HTTF experiments. This summary presents results from a portion of that benchmark using RELAP5-3D

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Effect of Inlet Throttling on Thermohydraulic Instability in a Large Scale Waterbased RCCS: A System-Level Analysis with RELAP5-3D

This paper presents results from system -level modeling of a water -based reactor cavity cooling system using RELAP5-3D. The computational model is benchmarked with experimental data from a half -scale RCCS test facility at Argonne National Laboratory. The model prediction is first compared with a two-phase oscillatory baseline experimental case where mixed accuracy is obtained. The model shows reasonable prediction of mass flow rate, pressure, and temperature but significant overprediction of void fraction. The model prediction is then compared with a fault case where the inlet of the risers is gradually reduced using a throttling valve. As the valve is closed, the model is able to predict some major flow phenomena observed in the experiment such as the dampening of oscillations, the reintroduction of oscillations, as well as boiling, flashing, and geysering in the risers. However, the timeline of these events are not well captured by the model. The model is also used to investigate the evolution of flow regime in the chimney. This work highlights that the semiempirical constitutive relations used in RELAP-3D could have a strong influence on the accuracy of the model in two-phase oscillatory flows.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

A Study of Propane Combustion in a Spark-Ignited Cooperative Fuel Research (CFR) Engine

Liquefied petroleum gas (LPG), whose primary composition is propane, is a promising candidate for heavy-duty vehicle applications as a diesel fuel alternative due to its CO 2 reduction potential and high knock resistance. To realize diesel-like efficiencies, spark-ignited LPG engines are proposed to operate near knock-limit over a wide range of operating conditions, which necessitates an investigation of fuel-engine interactions that leads to end-gas autoignition with propane combustion. This work presents both experimental and numerical studies of stoichiometric propane combustion in a sparkignited (SI) cooperative fuel research (CFR) engine. Engine experiments are initially conducted at different compression ratio (CR) values, and the effects of CR on engine combustion are characterized. A three-pressure analysis (TPA) model based on the two-zone combustion concept is developed in GT-Power and validated using test results to estimate in-cylinder wall temperatures, residual gas fraction, etc. This model is further utilized to examine end-gas chemistry by enabling the SI turbulent flame combustion and unburned gas chemical kinetics modules. Finally, a three-dimensional (3D) computational fluid dynamic (CFD) model of the CFR engine is developed in CONVERGE, where the G-equation and SAGE detailed chemical kinetics models are implemented for combustion modeling. Here, a 153 species reduced chemical kinetics mechanism derived from the detailed NUIGMech1.1 mechanism based on the ignition delay and laminar flame speed (LFS) studies is used to generate an LFS lookup table and to describe end-gas autoignition chemistry. Multi-cycle Reynolds-averaged Navier-Stokes (RANS) simulations are then performed for the tested CRs, and the numerical model is shown to be capable of predicting the propane combustion characteristics, particularly the end-gas autoignition behavior.

20 FOSSIL-FUELED POWER PLANTS↗

NASA NTRS1999

A Lagrangian Simulation of Subsonic Aircraft Exhaust Emissions

To estimate the effect of subsonic and supersonic aircraft exhaust on the stratospheric concentration of NO(y), we employ a trajectory model initialized with air parcels based on the standard release scenarios. The supersonic exhaust simulations are in good agreement with 2D and 3D model results and show a perturbation of about 1-2 ppbv of NO(y) in the stratosphere. The subsonic simulations show that subsonic emissions are almost entirely trapped below the 380 K potential temperature surface. Our subsonic results contradict results from most other models, which show exhaust products penetrating above 380 K, as summarized. The disagreement can likely be attributed to an excessive vertical diffusion in most models of the strong vertical gradient in NO(y) that forms at the boundary between the emission zone and the stratosphere above 380 K. Our results suggest that previous assessments of the impact of subsonic exhaust emission on the stratospheric region above 380 K should be considered to be an upper bound.

Schoeberl, M. R.↗

NASA NTRS2006

Modeling the Observed QBO and Inter-Annual Variations of the Diurnal Tide in the Mesosphere

In the current version of the Numerical Spectral Model (NSM), the Quasi-biennial Oscillation (QBO) is generated primarily by small-scale gravity waves (GW) from Hines' Doppler Spread Parameterization (DSP). The model does not have topography, and the planetary waves are solely generated by instabilities. We discuss a 3D modeling study that describes the QBO extending from the stratosphere into the upper mesosphere, where the oscillation produces significant inter-annual variations in the diurnal tide. The numerical results are compared with temperature measurements from the SABER (TIMED) and MLS (UARS) instruments obtained by Huang et al. (2006). With a GW source that peaks at the Equator and is taken to be isotropic and independent of season, the NSM generates a QBO with variable periods around 26 months and zonal wind amplitudes of almost 25 m/s at 30 km. As reported earlier, the NSM reproduces the observed equinoctial maxima in the diurnal tide at altitudes around 95 km. The modeled QBO modulates the tide such that the seasonal amplitude maxima can vary from one year to another by as much as 30%. To shed light on the underlying mechanisms, the relative importance of the advection terms are discussed, and they are shown to be important in the stratosphere. At altitudes above 80 km, however, the QBO-related inter-annual variations of the tide are generated primarily by GW momentum deposition. In qualitative agreement with the SABER measurements, the model generates distinct zonal-mean QBO temperature variations in the stratosphere and mesosphere. In the stratosphere, the computed amplitudes are not much smaller than those observed, and the rate of downward propagation at the Equator is reproduced. The modeled temperature amplitudes in the mesosphere, however, are much smaller than those observed. The observed and computed temperature variations of the QBO peak at the Equator but extend with phase reversals to high latitudes, in contrast to the zonal winds that are confined to equatorial latitudes. Hemispherical asymmetries also appear in both the model results and the observations. The temperature amplitudes outside the equatorial region however tend to occur at lower latitudes in the model results. While there is qualitative agreement between the TIMED measurements and the model prediction, there are some areas of significant disagreement that require us to reexamine the present version of the NSM. The numerical results critically depend on the chosen parameters that determine the wave forcing, and there are a number of avenues to improve the performance of the model that had not been tuned to fit the observations. The GW spectrum and its latitude dependence in the troposphere are not well known, and numerical experiments are discussed that describe the related model response. While it appears that eastward propagating Kelvin waves and westward propagating Rossby gravity waves are not the primary source to generate the QBO, the GW forcing can seed the oscillation and act as a catalyst to enhance effectiveness of these planetary waves.

Mayr, Hans G.↗

NASA NTRS2015

LEWICE3D/GlennHT Particle Analysis of the Honeywell Al502 Low Pressure Compressor

A flow and ice particle trajectory analysis was performed for the booster of the Honeywell AL502 engine. The analysis focused on two closely related conditions one of which produced a rollback and another which did not rollback during testing in the Propulsion Systems Lab at NASA Glenn Research Center. The flow analysis was generated using the NASA Glenn GlennHT flow solver and the particle analysis was generated using the NASA Glenn LEWICE3D v3.56 ice accretion software. The flow and particle analysis used a 3D steady flow, mixing plane approach to model the transport of flow and particles through the engine. The inflow conditions for the rollback case were: airspeed, 145 ms; static pressure, 33,373 Pa; static temperature, 253.3 K. The inflow conditions for the non-roll-back case were: airspeed, 153 ms; static pressure, 34,252 Pa; static temperature, 260.1 K. Both cases were subjected to an ice particle cloud with a median volume diameter of 24 microns, an ice water content of 2.0 gm3 and a relative humidity of 100 percent. The most significant difference between the rollback and non-rollback conditions was the inflow static temperature which was 6.8 K higher for the non-rollback case.

Ice Accretion Codes↗

NASA NTRS2019

Joint 3D-Wind Retrievals with Stereoscopic Views from MODIS and GOES

Atmospheric motion vectors (AMVs), derived by tracking patterns, represent the winds in a layer characteristic of the pattern. AMV height (or pressure), important for applications in atmospheric research and operational meteorology, is usually assigned using observed IR brightness temperatures with a modeled atmosphere and can be inaccurate. Stereoscopic tracking provides a direct geometric height measurement of the pattern that an AMV represents. We extend our previous work with multi-angle imaging spectroradiometer (MISR) and GOES to moderate resolution imaging spectroradiometer (MODIS) and the GOES-R series advanced baseline imager (ABI). MISR is a unique satellite instrument for stereoscopy with nine angular views along track, but its images have a narrow (380 km) swath and no thermal IR channels. MODIS provides a much wider (2330 km) swath and eight thermal IR channels that pair well with all but two ABI channels, o e ring a rich set of potential applications. Given the similarities between MODIS and VIIRS, our methods should also yield similar performance with VIIRS. Our methods, as enabled by advanced sensors like MODIS and ABI, require high-accuracy geographic registration in both systems but no synchronization of observations. AMVs are retrieved jointly with their heights from the disparities between triplets of ABI scenes and the paired MODIS granule. We validate our retrievals against MISR-GOES retrievals, operational GOES wind products, and by tracking clear-sky terrain. We demonstrate that the 3D-wind algorithm can produce high-quality AMV and height measurements for applications from the planetary boundary layer (PBL) to the upper troposphere, including cold-air outbreaks, wildfire smoke plumes, and hurricanes.

Carr, James L.↗

NASA NTRS2008

Multiscale Modeling of Ablation and Pyrolysis in PICA-Like materials

During atmospheric entry of planetary probes, the thermal protection system (TIPS) of the probe is exposed to high temperatures under low pressures. In these conditions, carbonous fibrous TIPS materials may undergo oxidation leading to mass loss and wall recession called ablation. This work aims to improve the understanding of material/environment interactions through a study of the coupling between oxygen transport in the Knudsen regime, heterogeneous oxidation of carbon, and surface recession. A 3D Random Walk Monte Carlo simulation tool is used for this study. The fibrous architecture of a model material, consisting of high porosity random array of carbon fibers, is numerically represented on a 3D Cartesian grid. Mass transport in the Knudsen regime from the boundary layer to the surface, and inside this porous material is simulated by random walk. A reaction probability is used to simulate the heterogeneous oxidation reaction. The surface recession of the fibers is followed by front tracking using a simplified marching cube approach. The output data of the simulations are ablation velocity and dynamic evolution of the material porosity. A parametric study is carried out to analyze the material behavior as a function of Knudsen number for the porous media (length of the mean free path compared to the mean pore diameter) and the intrinsic reactivity of the carbon fibers. The model is applied to Stardust mission reentry conditions and explains the unexpected behavior of the TIPS material that underwent mass loss in volume.

Lachaud, Jean↗

DOE OSTI2023

Insights on Simulating Summer Warming of the Great Lakes: Understanding the Behavior of a Newly Developed Coupled Lake-Atmosphere Modeling System

The Laurentian Great Lakes are the world's largest freshwater system and regulate the climate of the Great Lakes region, which has been increasingly experiencing climatic, hydrological, and ecological changes. An accurate mechanistic representation of the Great Lakes thermal structure in Regional Climate Models (RCMs) is paramount to studying the climate of this region. Currently, RCMs have primarily represented the Great Lakes through coupled one-dimensional (1D) column lake models; this approach works well for small inland lakes but is unable to resolve the realistic hydrodynamics of the Great Lakes and leads to inaccurate representations of lake surface temperature (LST) that influence regional climate and weather patterns. This work overcomes this limitation by developing a fully two-way coupled modeling system using the Weather Research and Forecasting model and a three-dimensional (3D) hydrodynamic model. The coupled model system resolves the interactive physical processes between the atmosphere, lake, and surrounding watersheds; and validated against a range of observational data. The model is then used to investigate the potential impacts of lake-atmosphere coupling on the simulated summer LST of Lake Superior. By evaluating the difference between our two-way coupled modeling system and our observation-driven modeling system, we find that coupled-lake atmosphere dynamics can lead to a higher LST during June-September through higher net surface heat flux entering the lake in June and July and a lower net surface heat flux entering the lake in August and September. The unstratified water in June distributes the entering surface heat flux throughout the water column leading to a minor LST increase, while the stratified waters of July create a conducive thermal structure for the water surface to warm rapidly under the higher incoming surface heat flux. This research provides insight into the coupled modeling system behavior, which is critical for enhancing our predictive understanding of the Great Lakes climate system.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Volumetric recombination in EMC3-EIRENE: Implementation and first application to the pre-fusion power operation phase in ITER

An upgrade to the 3D plasma boundary model EMC3-EIRENE for detached divertor plasmas is introduced and evaluated. Stabilization of the iterative approximation of a self-consistent plasma–neutral gas solution is required at low divertor plasma temperatures, and this has been achieved by linearization of the energy loss term associated with the electron–neutral gas interaction. An implicit method for volumetric recombination has been implemented, which improves agreement with SOLPS-ITER results for the pre-fusion power operation phase in ITER. Resonant magnetic perturbations for control of edge localized modes are found to result in a partially detached divertor plasma with a significantly broader reduction of particle loads in the traditional strike zone, which is attributed to volumetric recombination. The additional non-axisymmetric far scrape-off layer strike points, on the other hand, remain attached at relatively high temperature so that volumetric recombination is not effective there.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2010

Computational Fluid Dynamics (CFD) Simulations of Jet Mixing in Tanks of Different Scales

For long-duration in-space storage of cryogenic propellants, an axial jet mixer is one concept for controlling tank pressure and reducing thermal stratification. Extensive ground-test data from the 1960s to the present exist for tank diameters of 10 ft or less. The design of axial jet mixers for tanks on the order of 30 ft diameter, such as those planned for the Ares V Earth Departure Stage (EDS) LH2 tank, will require scaling of available experimental data from much smaller tanks, as well designing for microgravity effects. This study will assess the ability for Computational Fluid Dynamics (CFD) to handle a change of scale of this magnitude by performing simulations of existing ground-based axial jet mixing experiments at two tank sizes differing by a factor of ten. Simulations of several axial jet configurations for an Ares V scale EDS LH2 tank during low Earth orbit (LEO) coast are evaluated and selected results are also presented. Data from jet mixing experiments performed in the 1960s by General Dynamics with water at two tank sizes (1 and 10 ft diameter) are used to evaluate CFD accuracy. Jet nozzle diameters ranged from 0.032 to 0.25 in. for the 1 ft diameter tank experiments and from 0.625 to 0.875 in. for the 10 ft diameter tank experiments. Thermally stratified layers were created in both tanks prior to turning on the jet mixer. Jet mixer efficiency was determined by monitoring the temperatures on thermocouple rakes in the tanks to time when the stratified layer was mixed out. Dye was frequently injected into the stratified tank and its penetration recorded. There were no velocities or turbulence quantities available in the experimental data. A commercially available, time accurate, multi-dimensional CFD code with free surface tracking (FLOW-3D from Flow Science, Inc.) is used for the simulations presented. Comparisons are made between computed temperatures at various axial locations in the tank at different times and those observed experimentally. The affect of various modeling parameters on the agreement obtained are assessed.

Breisacher, Kevin↗

DOE OSTI2023

Bentonite buffer under high temperature: laboratory experiments and coupled process modeling

Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.

Zheng, Liange↗

NASA NTRS2012

Axisymmetric Plume Simulations with NASA's DSMC Analysis Code

A comparison of axisymmetric Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) results to analytic and Computational Fluid Dynamics (CFD) solutions in the near continuum regime and to 3D DAC solutions in the rarefied regime for expansion plumes into a vacuum is performed to investigate the validity of the newest DAC axisymmetric implementation. This new implementation, based on the standard DSMC axisymmetric approach where the representative molecules are allowed to move in all three dimensions but are rotated back to the plane of symmetry by the end of the move step, has been fully integrated into the 3D-based DAC code and therefore retains all of DAC s features, such as being able to compute flow over complex geometries and to model chemistry. Axisymmetric DAC results for a spherically symmetric isentropic expansion are in very good agreement with a source flow analytic solution in the continuum regime and show departure from equilibrium downstream of the estimated breakdown location. Axisymmetric density contours also compare favorably against CFD results for the R1E thruster while temperature contours depart from equilibrium very rapidly away from the estimated breakdown surface. Finally, axisymmetric and 3D DAC results are in very good agreement over the entire plume region and, as expected, this new axisymmetric implementation shows a significant reduction in computer resources required to achieve accurate simulations for this problem over the 3D simulations.

Stewart, B. D.↗

DOE OSTI2021

A measure of active interfaces in supported catalysts for high-temperature reactions

Formulating knowledge of structure-function relationships in heterogeneous catalysis is central to the design of efficient catalysts; yet, the elucidation of dominant reaction sites has remained as a challenge. Here, in this paper, we present a methodology that can be used to visualize metal-gas and metal-oxide-gas interfaces in three dimensions and to quantify their catalytic activity levels. As a case study, CH 4 oxidation occurring in a Pt/CeO 2 system is chosen. By employing thermally robust Pt@CeO 2 model catalysts with size-tunable and monodisperse cores, and gas-permeable shells, we reconstruct a series of structures in 3D via electron tomography and match the information to activity data and theoretical calculations. This strategy reveals that two different interfaces catalyze the CH 4 oxidation and that their contribution to the overall rate changes with the Pt size, temperature, and gas atmosphere. Our results provide an analytic platform on which to explore reaction pathways and mechanisms applicable to multiple reactions and materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Modeling Nondestructive Defect Detection in Additively Manufactured Metallic Structures for Nuclear Applications

The future of quickly, economically produced metallic nuclear reactor parts with minimal supply-chain dependence lies in Laser Powder Bed Fusion (LPBF) Additive Manufacturing (AM): a 3D printing method involving laser melting and net shaping stainless steel and Inconel metallic powder into a solid structure. However, intrinsic features in LPBF frequently leads to the formation of materials defects, such as pores, within 3D printed structures. As safe long-term use in energy applications requires knowledge of all relevant defects before deployment in a reactor, we must develop methods for nondestructive detection of these defects. We are investigating Pulsed Thermal Tomography (PTT), which is a non-contact nondestructive imaging method scalable to arbitrary structure size. Thermal tomography (TT) is a computational method for 3D spatial reconstruction of material thermal effusivity from flash or pulsed thermography temperature data cube. Thermography data cube consists of 2D surface temperature measurements at different times. The objective of the present work is to investigate limits on defect detection in AM metallic structures with PTT. To this effect, we modeled PTT with COMSOL heat transfer computer simulations. We developed a layered media COMSOL simulation consisting of a Stainless Steel 316 (SS316) plate with an internal layer of un-sintered SS316 powder. Thermophysical properties of the powder layer were modeled with equivalent volume mixing model. To account for partial sintering at the boundary of the defect, the transition between solid and powder layers was modeled as a Gaussian. Using data from COMSOL simulations, we reconstructed depth-dependent thermal effusivity, which allowed defect visibility estimation. A series of parametric studies determined that at 1mm depth, 50µm is the smallest detectable defect. In addition, classification of the defects which can lead to early fatigue of the metallic structure in a reactor is briefly discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Artificial correlation heating in PIC simulations

The Particle-in-Cell (PIC) method, a cornerstone in plasma modeling, is widely employed for its ability to simulate kinetic phenomena in device-scale domains. Part of what makes this possible is that computational macroparticles represent many physical particles. It converges under certain constraints, including a grid spacing that resolves the Debye length and a time step small enough to respect the Courant–Friedrichs–Lewy condition and plasma frequency stability limit. Here, we introduce a new constraint necessary to avoid Artificial Correlation Heating (ACH). This requires that the macroparticle coupling strength be smaller than one, Γ ω < 1, where Γ ω ≡ Γω 2/3 ⁠, Γ = Z 2 e 2 /(4πε 0 ak B T) is the physical coupling strength, and w is the macroparticle weight. This is particularly relevant to 3D simulations of dense plasmas, which are becoming common with modern computing power. If this condition is violated, the finite macroparticle weight artificially enhances the coupling strength and causes the plasma to heat until the macroparticle coupling strength is near unity, depending on the grid resolution. A comprehensive model of ACH is developed that incorporates electron density, temperature, macroparticle weight, and grid resolution. It is then tested using PIC simulations, delineating the boundaries of the method's applicability and offering a predictive framework for ACH. Moreover, the research explores a runaway heating process induced by ACH in the presence of ionization, which can lead to numerical instability. A conclusion of this study is that the onset of ACH can impose a more stringent constraint on the macroparticle weight and average number of macroparticles per cell than what is typically expected, particularly in 3D simulations of dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Hybrid AI-ML and FE-based Digital Twin Predictive Modeling Framework for a PWR Coolant System Components: Updates on Multi-Time-Series-3D-Location Dependent Usages Factor Prediction

The long-term operation (LTO) of nuclear power plants (NPP) beyond their original design life of 40 years can lead to more material damage associated with cyclic fatigue under thermal-mechanical loading cycles and associated long-term exposure of reactor material to the deleterious reactor-coolant environments. However, under this LTO condition, the reactor components can still safely operate but may require more frequent Nondestructive Evaluation (NDE) of reactor components. Requiring frequent NDE inspections may lead to frequent NPP shutdowns which can lead to power outages and additional NDE inspection cost-related economic loss. The economic loss can be minimized by reducing uncertainty in life estimation of safety-critical pressure boundary components and by implementing a more digital approach such as using upcoming digital-twin (DT) technology for predicting the structural states (e.g., time and location dependent inside/outside thickness temperature, stress, strain, plastic deformation, etc.) and associated fatigue life of a component in real time. The DT framework is based on limited experimental data, Artificial-intelligence (AI)-Machine-Learning (ML) and multiphysics-computationalmechanics such as finite element- (FE) based models. Given the real-time thermal-hydraulic process measurements from several existing plant sensors, the overall goal of the DT framework is to predict the cumulative usages factors or equivalent fatigue lives in real time and at any random 3D location of the components. This includes inaccessible locations such as inside the thickness or location of a component. This prediction can be at thousands to millions of 3D point clouds or locations like conventional FE-based models, but without running an FE model in real time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗