Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled 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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Effects of Surface Turbulence Flux Parameterizations on the MJO: The Role of Ocean Surface Waves

This study investigates the sensitivity of the Madden–Julian oscillation (MJO) to changes to the bulk flux parameterization and the role of ocean surface waves in air–sea coupling using a fully coupled ocean–atmosphere–wave model. The atmospheric and ocean model components of the Energy Exascale Earth System Model (E3SM) are coupled to a spectral wave model, WAVEWATCH III (WW3). Two experiments with wind speed–dependent bulk algorithms (NCAR and COARE3.0a) and one experiment with wave-state-dependent flux (COR3.0a-WAV) were conducted. We modify COARE3.0a to include surface roughness calculated within WW3 and also account for the buffering effect of waves on the relative difference between air-side and ocean-side momentum flux. Differences in surface fluxes, primarily caused by discrepancies in drag coefficients, result in significant differences in MJO’s properties. While COARE3.0a has better convection–circulation coupling than NCAR, it exhibits anomalous MJO convection east of the date line. The wave-state-dependent flux (COR3.0-WAV) improves the MJO representation over the default COARE3.0 algorithm. Strong easterlies over the Pacific Ocean in COARE3.0a enhance the latent heat flux (LHFLX). This is responsible for the anomalous MJO propagation after the date line. In COR3.0a-WAV, waves reduce the anomalous easterlies, leading to a decrease in LHFLX and MJO dissipation after the date line. These findings highlight the role of surface fluxes in MJO simulation fidelity. Most importantly, we show that the proper treatment of wave-induced effects in bulk flux parameterization improves the simulation of coupled climate variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

A comparative study on magnetism in Zn-doped AlN and GaN from first-principles

First-principles calculations have been used to comparatively investigate electronic and magnetic properties of Zn-doped AlN and GaN. A total magnetic moment of 1.0 μ B {sub B} induced by Zn is found in AlN, but not in GaN. Analyses show that the origin of spontaneous polarization not only depend on the localized atomic orbitals of N and sufficient hole concentration, but also the relative intensity of the covalency of matrix. The relatively stronger covalent character of GaN with respect to AlN impedes forming local magnetic moment in GaN matrix. Our study offers a fresh sight of spontaneous spin polarization in d⁰ magnetism. The much stronger ferromagnetic coupling in c-plane of AlN means that it is feasible to realize long-range ferromagnetic order via monolayer delta-doping. This can apply to other wide band-gap semiconductors in wurtzite structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Integrated Commercial Carbon Capture and Storage (CCS) Prefeasibility Study at Rock Springs Uplift, Wyoming (Final Report)

To establish the potential of safe, long-term commercial-scale carbon capture & storage (CCS) in southwestern Wyoming, the University of Wyoming formed a CCS coordination team (CCT) to complete a Phase I CarbonSAFE pre-feasibility study adjacent to the region’s largest coal-fired power plant, the Jim Bridger Plant (JBP). The study’s primary objective was to evaluate the potential of storing 50+ million metric tons of anthropogenic carbon dioxide (CO 2 ) in stacked reservoirs over a period of 25 years. In addition, this integrated study evaluated related aspects of carbon capture utilization & storage (CCUS) at the study site, including: (1) the CO 2 character of post-combustion flue gas from JBP and the technical capacity for retrofitting proven commercial-scale CO 2 capture technology relative to CO 2 source character and technological constraints; (2) utilization of the existing CO 2 pipeline network in the immediate vicinity of JBP for enhanced oil recovery (CO 2 -EOR) and utilizing that network as part of a regional hub for CO 2 storage from other anthropogenic sources; (3) assessment of the challenges and benefits of meeting all of CarbonSAFE programmatic Phase I goals relative to Wyoming’s carbon management regulatory framework; and (4) assessment of the community and environmental factors that may impact CCS at JBP. With respect to storage, we analyzed the potential of stacked Mesozoic reservoirs at the study site, focusing on the Entrada and Nugget sandstones. The reservoir response to injection simulations varied greatly between formations, though they both were capable of CO 2 storage. The Entrada Sandstone is relative thin (~40’) and heterogeneous; storage using one injection and production well was limited to 1.2 Mt over 25 years at the study site. The Nugget Sandstone is much thicker (> 400’) and generally homogeneous; storage using one injection and production well approached 15.0 Mt over 25 years; this scenario suggests that the storage potential in the Nugget Sandstone is close to the estimated P90 calculation of 9.6 Mt/mi 2 as models indicate a storage potential of 8.3 Mt/mi 2 at saturation. Both reservoirs’ storage capacity benefited greatly from coupling an injection well with an offset production well; total storage capacity increased and reservoir pressure stayed below critical thresholds. This indicates that implementing a pressure management strategy at the study site would optimize storage and decrease risk. Seals associated with the targeted reservoirs were shown to be capable of retaining the potential storage volumes, and the site benefits from multiple seal redundancies (i.e. > 7,000’ of seal formations). Furthermore, assessment of different reservoir variables (e.g., fluid composition, pressure, etc.) indicates all reservoirs at the study site are isolated. Storage assessments from Mesozoic reservoirs were coupled with simulations from the site’s Paleozoic reservoirs (the Weber Sandstone and Madison Limestone), which were the focus of previous study. The fully integrated stacked reservoir pre-feasibility studies suggest that the Nugget and Madison formations alone could safely store 50 million metric tons of CO 2 within the study site. Furthermore, both formations are true saline reservoirs near JBP, making them ideal target reservoirs for CO 2 sequestration. Analysis of Wyoming’s laws and regulations confirm that the State has a favorable policy environment for CCS/CCUS. The project’s economic models estimated that the following revenues collectively are sufficient to finance the project over its lifespan if CO 2 utilization is coupled with CCS via the following strategies: (1) sales of CO 2 for EOR (approximately $69 million/year); (2) use of CO 2 tax credits such as amended §45Q and §48A (approximately $484 million total); and (3) sales of low-carbon electricity and marketable carbon offset/credits into carbon-constrained markets (approximately $11-$17 million/year) to West Coast consumers. These revenue estimates are broadly consistent with other studies that concluded arbitrage in JBP retail electricity sales between California and Wyoming could help to support the cost of deploying CCS at JBP. The models also indicate that maintaining a successful business case for long-term CCS -- wherein associated capture, transportation and operational costs are fully addressed at JBP -- is reliant on numerous external factors that complicate predictions. This study suggests that this site meets the CarbonSAFE program’s requirements of being able to feasibly store 50+ million metric tons of CO 2 over 25 years within the site’s stacked reservoirs, especially if coupled with pressure management strategies. This study suggests that the Nugget Sandstone and the Madison Limestone are two of the best reservoirs in the State with respect to overall storage capacity. In addition, the site’s proximity to existing CO 2 transportation networks, CO 2 -EOR opportunities, ability to sell to markets that value low-carbon electricity, Wyoming’s existing carbon regulatory framework, and a public that is well-educated with respect to energy markets and issues collectively increase the potential for implementing commercial-scale CCS adjacent to the JBP. With respect to the available technical and non-technical site data, and the reservoir quality at the study site, we suggest that this site is capable of moving directly into Phase III of the CarbonSAFE program.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2019

Kinetic particle simulations in a global toroidal geometry

The gyrokinetic toroidal code has been upgraded for global simulations by coupling the core and scrape-off layer regions across the separatrix with field-aligned particle-grid interpolations. A fully kinetic particle pusher for high frequency waves (ion cyclotron frequency and beyond) and a guiding center pusher for low frequency waves have been implemented using cylindrical coordinates in a global toroidal geometry. The two integrators correctly capture the particle orbits and agree well with each other, conserving energy and canonical angular momentum. As a verification and application of this new capability, ion guiding center simulations have been carried out to study ion orbit losses at the edge of the DIII-D tokamak for single null magnetic separatrix discharges. The ion loss conditions are examined as a function of the pitch angle for cases without and with a radial electric field. We report the simulations show good agreement with past theoretical results and with the experimentally observed feature in which high energy ions flow out along the ion drift orbits and then hit the divertor plates. A measure of the ion direct orbit loss fraction shows that the loss fraction increases with the ion energy for DIII-D in the initial velocity space. Finally, as a further verification of the capability of the new code, self-consistent simulations of zonal flows in the core region of the DIII-D tokamak were carried out. All DIII-D simulations were performed in the absence of turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2020

Neural Network Reflectance Prediction Model for Both Open Ocean and Coastal Waters

Remote sensing of global ocean color is a valuable tool for understanding the ecology and biogeochemistry of the worlds oceans, and provides critical input to our knowledge of the global carbon cycle and the impacts of climate change. Ocean polarized reflectance contains information about the constituents of the upper ocean euphotic zone, such as colored dissolved organic matter (CDOM), sediments, phytoplankton, and pollutants. In order to retrieve the information on these constituents, remote sensing algorithms typically rely on radiative transfer models to interpret water color or remote-sensing reflectance; however, this can be resource-prohibitive for operational use due to the extensive CPU time involved in radiative transfer solutions. In this work, we report a fast model based on machine learning techniques, called Neural Network Reflectance Prediction Model (NNRPM), which can be used to predict ocean bidirectional polarized reflectance given inherent optical properties of ocean waters. This supervised model is trained using a large volume of data derived from radiative transfer simulations for coupled atmosphere and ocean systems using the successive order of scattering technique (SOS-CAOS). The performance of the model is validated against another large independent test dataset generated from SOS-CAOS. The model is able to predict both polarized and unpolarized reflectances with an absolute error (AE) less than 0.004 for 99% of test cases. We have also shown that the degree of linear polarization (DoLP) for unpolarized incident light can be predicted with an AE less than 0.002 for 99% of test cases. In general, the simulation time of SOS-CAOS depends on optical depth, and required accuracy. When comparing the average speeds of the NNRPM against the SOS-CAOS model for the same parameters, we see that the NNRPM is able to predict the Ocean BRDF 6000 times faster than SOS-CAOS. Both ultraviolet and visible wavelengths are included in the model to help differentiate between dissolved organic material and chlorophyll in the study of the open ocean and the coastal zone. The incorporation of this model into the retrieval algorithm will make the retrieval process more efficient, and thus applicable for operational use with global satellite observations.

radiative transfer↗
DOE OSTI2021

A Review of Lattice-Boltzmann Models Coupled with Geochemical Modeling Applied for Simulation of Advanced Waterflooding and Enhanced Oil Recovery Processes

To maintain economic profit and improve the oil production efficiency after the primary and secondary production phase, advanced waterflooding techniques such as low salinity waterflooding in carbonate reservoirs have been investigated in numerical simulations, laboratory experiments, and field pilot tests. Multiple underlying mechanisms have been proposed based on these studies, and they are still under debate. Various numerical modeling approaches are introduced, but there exists a lack of a pore-scale comprehensive modeling scheme to fully understand the processes. Lattice-Boltzmann method (LBM) is a type of numerical fluid flow modeling technique that shows capabilities and flexibilities in modeling pore-scale fluid flow to integrate physical–chemical processes within complex structures. The intrinsic feature of LBM makes it a promising framework for simulating advanced waterflooding due to its flexibility, accuracy, and parallel efficiency. LBM works either by itself for solving reactive transport problems or by coupling with a third-party reaction solver. This review mainly introduces the LBM fluid flow and reactive transport capabilities and the concept and modeling approaches to simulate advanced waterflooding techniques. Meanwhile, an evaluation of the coupled LBM models for enhanced oil recovery (EOR) simulations is discussed with future research challenges and directions concluded.

02 PETROLEUM↗
NASA NTRS1996

Extinguishment of a Diffusion Flame Over a PMMA Cylinder by Depressurization in Reduced-Gravity

Extinction of a diffusion flame burning over horizontal PMMA (Polymethyl methacrylate) cylinders in low-gravity was examined experimentally and via numerical simulations. Low-gravity conditions were obtained using the NASA Lewis Research Center's reduced-gravity aircraft. The effects of velocity and pressure on the visible flame were examined. The flammability of the burning solid was examined as a function of pressure and the solid-phase centerline temperature. As the solid temperature increased, the extinction pressure decreased, and with a centerline temperature of 525 K, the flame was sustained to 0.1 atmospheres before extinguishing. The numerical simulation iteratively coupled a two-dimensional quasi-steady, gas-phase model with a transient solid-phase model which included conductive heat transfer and surface regression. This model employed an energy balance at the gas/solid interface that included the energy conducted by the gas-phase to the gas/solid interface, Arrhenius pyrolysis kinetics, surface radiation, and the energy conducted into the solid. The ratio of the solid and gas-phase conductive fluxes Phi was a boundary condition for the gas-phase model at the solid-surface. Initial simulations modeled conditions similar to the low-gravity experiments and predicted low-pressure extinction limits consistent with the experimental limits. Other simulations examined the effects of velocity, depressurization rate and Phi on extinction.

Goldmeer, Jeffrey Scott↗
DOE OSTI2016

ENSO-driven energy budget perturbations in observations and CMIP models

Various observation-based datasets are employed to robustly quantify changes in ocean heat content (OHC), anomalous ocean–atmosphere energy exchanges and atmospheric energy transports during El Niño-Southern Oscillation (ENSO). These results are used as a benchmark to evaluate the energy pathways during ENSO as simulated by coupled climate model runs from the CMIP3 and CMIP5 archives. The models are able to qualitatively reproduce observed patterns of ENSO-related energy budget variability to some degree, but key aspects are seriously biased. Area-averaged tropical Pacific OHC variability associated with ENSO is greatly underestimated by all models because of strongly biased responses of net radiation at top-of-the-atmosphere to ENSO. The latter are related to biases of mean convective activity in the models and project on surface energy fluxes in the eastern Pacific Intertropical Convergence Zone region. Moreover, models underestimate horizontal and vertical OHC redistribution in association with the generally too weak Bjerknes feedback, leading to a modeled ENSO affecting a too shallow layer of the Pacific. Vertical links between SST and OHC variability are too weak even in models driven with observed winds, indicating shortcomings of the ocean models. Furthermore, modeled teleconnections as measured by tropical Atlantic OHC variability are too weak and the tropical zonal mean ENSO signal is strongly underestimated or even completely missing in most of the considered models. In conclusion, results suggest that attempts to infer insight about climate sensitivity from ENSO-related variability are likely to be hampered by biases in ENSO in CMIP simulations that do not bear a clear link to future changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Future projection of mean and variability of the Asian Summer Monsoon and Indian Ocean Climate systems

The overall goal of this project is to assess the ability of the CMIP3/5 models to simulate the Indian-Ocean monsoon systems. The PI along with post-docs investigated research issues ranging from synoptic systems to long-term trends over the Asian monsoon region. The PI applied diagnostic tools such as moist static energy (MSE) to isolate: the moist and radiative processes responsible for extended monsoon breaks over South Asia, precursors in the ENSO-monsoon association, reasons for the drying tendency over South Asia and the possible effect on tropical Indian Ocean climate anomalies influencing certain aspects of ENSO characteristics. By diagnosing various observations and coupled model simulations, we developed working hypothesis and tested them by carrying out sensitivity experiments with both linear and nonlinear models. Possible physical and dynamical reasons for model sensitivities were deduced. On the teleconnection front, the ability of CMIP5 models in representing the monsoon-desert mechanism was examined recently. Further more, we have applied a suite of diagnostics and have performed an in depth analysis on CMIP5 integrations to isolate the possible reasons for the ENSO-monsoon linkage or lack thereof. The PI has collaborated with Dr. K.R. Sperber of PCMDI and other CLIVAR Asian-Australian monsoon panel members in understanding the ability of CMIP3/5 models in capturing monsoon and its spectrum of variability. The objective and process-based diagnostics aided in selecting models that best represent the present-day monsoon and its variability that are then employed for future projections. Two major highlights were an invitation to write a review on present understanding monsoons in a changing climate in Nature Climate Change, and identification of an east-west shift in observed monsoon rainfall (more rainfall over tropical western Pacific and drying tendency over South Asia) in the last six decades and attributing that shift to SST rise over the tropical western Pacific. On the training of post-doctoral scientists: the PI spent considerable amount of time and efforts in introducing the post-docs into climate modeling and designing the numerical experiments. With training provided and knowledge gained, post-docs worked in the project obtained long term positions elsewhere. The PI also enjoyed the experience in managing the works and educating work ethics to the younger generation. Based on the research achievements and publications, the PI gave invited talks in major international monsoon conferences/workshops, and gave lectures in various research organizations in the last six years. Finally, during the project period, the PI attended all the DOE organized PIs meeting and presented the major results. Some of the major implications of the project include: (i) Sustained observational efforts are necessary to monitor the three-dimensional moisture distribution over the Asian monsoon region that would aid in better understanding, modeling and predicting severe monsoons well in advance and (ii) process-based diagnostics lead pathways for model improvements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Physics-based multiscale coupling for full core nuclear reactor simulation

Numerical simulation of nuclear reactors is a key technology in the quest for improvements in efficiency, safety, and reliability of both existing and future reactor designs. Historically, simulation of an entire reactor was accomplished by linking together multiple existing codes that each simulated a subset of the relevant multiphysics phenomena. Recent advances in the MOOSE (Multiphysics Object Oriented Simulation Environment) framework have enabled a new approach: multiple domain-specific applications, all built on the same software framework, are efficiently linked to create a cohesive application. This is accomplished with a flexible coupling capability that allows for a variety of different data exchanges to occur simultaneously on high performance parallel computational hardware. Examples based on the KAIST-3A benchmark core, as well as a simplified Westinghouse AP-1000 configuration, demonstrate the power of this new framework for tackling—in a coupled, multiscale manner—crucial reactor phenomena such as CRUD-induced power shift and fuel shuffle. 2014 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-SA license

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

A Phase-Field Model for the U-Pu-Zr System

Modeling nuclear materials is a multi-scale, multi-physics endeavor. Researchers are currently investigating how bulk material properties are affected by the evolution of microscopic material structures, i.e. microstructures. Development of tools that can mechanistically predict the effects of processes like diffusion and grain growth on properties like thermal conductivity will transform nuclear material design, manufacturing, and use. Intermediate-scale, i.e. mesoscale, modeling is among the most promising tools available for this task. Expedited development of mesoscale models, however, will require structured collaboration between modelers and experimentalists. Coupling these fields promises to improve both model accuracy and experiment design. With these ideas in mind, a mesoscale model of the U-Pu-Zr system is being developed. Its development will accompany a series of experiments involving irradiation of metallic fuel samples.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Dynamics of a network of phase oscillators with plastic couplings

The processes of synchronization and phase cluster formation are investigated in a complex network of dynamically coupled phase oscillators. Coupling weights evolve dynamically depending on the phase relations between the oscillators. It is shown that the network exhibits several types of behavior: the globally synchronized state, two-cluster and multi-cluster states, different synchronous states with a fixed phase relationship between the oscillators and chaotic desynchronized state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Residual Stress Asymmetry in Thin Sheets of Double-Sided Shot Peened Aluminum

The local mechanical properties and residual stresses of a double-sided shot-peened aluminum alloy AA7050-T7451 sheet were measured by nanoindentation and simulated by the finite element method (FEM). The variation of the pileup area around the indentation due to residual stresses, as well as changes in hardness as a function of cold working, was assessed. Numerical simulations were conducted to investigate the residual stress profile. When the indentation contact area was corrected to account for pileup induced from residual stresses, the hardness increased about 10% over the range of depths at which compressive residual stresses exist, indicating some strain hardening did occur. The maximum compressive residual stress for the second shot-peened surface is on the order of 10% higher than the first surface; this is ascribed to the sequencing during the peening process, in which the first side peening influenced the stress throughout the part thickness for these thin plate structures. The ability to couple FEM simulations and rapid experimental assessments of residual stresses and hardness allows for future prediction of part thickness effects in double-sided peening processes.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Spray characterization for engine combustion network Spray G injector using high-fidelity simulation with detailed injector geometry

This article presents a computational fluid dynamics study of the engine combustion network Spray G, focusing on the transient characteristics of the spray during the start of injection and the impacts of nozzle geometry details derived from the manufacturing process. The large-eddy-simulation method, coupled with the volume-of-fluid method, was used to model the high-speed turbulent two-phase flow. A moving-needle boundary condition was applied to capture the internal flow boundary condition accurately. The injector geometry was measured with micron-level resolution using X-ray tomographic imaging at the Advanced Photon Source at Argonne National Laboratory, providing detailed machining tolerance and defects from manufacturing and a realistic rough surface. For comparison, a nominal geometry and a modified geometry incorporating measured large-scale geometric features but no surface details were also used in the simulations. Spray characteristics such as mass flow rate, injection velocity, and Sauter mean diameter were analyzed. Significantly distinct spray characteristics in terms of injection velocity, spray morphology, and primary breakup mechanism were predicted using different nozzle geometries, which is mainly attributable to the realistic surface finish and manufacturing defects. The measured high-resolution geometry predicts a lower injection velocity, a wider-spreading spray, and an overall slower breakup rate with evident injector tip wetting compared to the ideally smooth nozzle boundary. This result implies that the manufacturing details of the injector, which are usually ignored in fuel injection studies, have a significant impact on the spray development process and should be taken into account for design optimization.

Yue, Zongyu↗
DOE OSTI2017

Black holes in vector-tensor theories

We study static and spherically symmetric black hole (BH) solutions in second-order generalized Proca theories with nonminimal vector field derivative couplings to the Ricci scalar, the Einstein tensor, and the double dual Riemann tensor. We find concrete Lagrangians which give rise to exact BH solutions by imposing two conditions of the two identical metric components and the constant norm of the vector field. These exact solutions are described by either Reissner-Nordström (RN), stealth Schwarzschild, or extremal RN solutions with a non-trivial longitudinal mode of the vector field. We then numerically construct BH solutions without imposing these conditions. For cubic and quartic Lagrangians with power-law couplings which encompass vector Galileons as the specific cases, we show the existence of BH solutions with the difference between two non-trivial metric components. The quintic-order power-law couplings do not give rise to non-trivial BH solutions regular throughout the horizon exterior. The sixth-order and intrinsic vector-mode couplings can lead to BH solutions with a secondary hair. For all the solutions, the vector field is regular at least at the future or past horizon. The deviation from General Relativity induced by the Proca hair can be potentially tested by future measurements of gravitational waves in the nonlinear regime of gravity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Large-eddy simulation of a multi-injection flame in a diesel engine environment using an unsteady flamelet/progress variable approach

In this work, large-eddy simulations (LESs) are conducted for a multiple-injection flame in a diesel engine environment using an unsteady flamelet/progress variable (UFPV) approach in which differential diffusion is considered. The suitability of the UFPV tabulation approach is first evaluated through a priori analyses using the state-of-the-art direct numerical simulation (DNS) dataset. Both the instantaneous data and the conditional values for the major and minor species' mass fractions are compared between the UFPV and the DNS. The comparisons show that the proposed UFPV tabulation approach can give good predictions for the multiple-injection flame at different injection phases. While the gas temperature and major species mass fractions can be accurately predicted with or without differential diffusion being considered in the UFPV flamelet library, the prediction accuracy for the highly diffusive species (e.g., hydrogen) in the main injection phase can be noticeably improved when differential diffusion is taken into account. The fully coupled LES/UFPV simulations show that the overall structure of the multiple-injection flame can be predicted, and the conditional thermo-chemical values are close to the filtered DNS dataset. The reasons for the remaining discrepancies found in the a priori analyses and the a posteriori simulations using the UFPV approach are analyzed.

Wen, Xu (ORCID:0000000331334292)↗
DOE OSTI2025

Bayesian Calibration of Nuclear Graphite Property Models Accounting for Model Inadequacy and Impacts on Component Performance

Nuclear-grade structural graphite is extensively utilized in the core designs of various advanced nuclear reactors. In the reactor environment, graphite is subjected to prolonged exposure to extreme conditions, including high temperatures, radiation, and potentially molten salt and oxygen. Such exposure can induce several degradation mechanisms in graphite, including nonuniform volumetric strains caused by irradiation and thermal expansion, leading to stresses that may compromise the performance of graphite components. Assessing component integrity requires accurate models of graphite's thermomechanical response. This report documents the Bayesian calibration of thermomechanical properties for nuclear-grade graphite and their application to graphite component modeling and simulation using the Grizzly code. As part of this work, uncertainty-quantified models were developed for the elastic modulus, coefficient of thermal expansion, irradiation-induced dimensional change, and irradiation-induced creep for graphite grades IG-110, NBG-18, NBG-17, PCEA, and 2114. Using a hierarchical Bayesian approach, multiple experimental data sources were combined to develop Gaussian process models for the properties. Using the Kennedy O'Hagan framework, the uncertainties due to inadequacies in the model and the inherent spread in the experimental data were quantified for three different models. These uncertainty-quantified models, with a model-form correction, were subsequently applied to a coupled-physics simulation of representative graphite components, revealing that the uncertainties have a large impact on the components' deformation.

36 - MATERIALS SCIENCE↗
DOE OSTI2017

Thermo-mechanical performance assessment of selected plates from MP-1 high power experiments

Monolithic plate-type fuel is a fuel form that is being developed for the conversion of high performance research and test reactors to low-enrichment uranium fuels. These fuel-plates are comprised of a high density, low enrichment, U-Mo alloy based fuel foil encapsulated in an aluminum cladding. To benchmark this new design, number of plates has been irradiated with satisfactory performance. As a part of continuing evaluation efforts, a set of plates covering range of operational parameters is scheduled to be tested during MP-1 irradiation experiments. It is necessary to evaluate the thermo-mechanical performance of plates during irradiation. For this, selected plates with distinct operational histories; covering low power, high power and high fission density were simulated. Fully coupled three-dimensional models of plates with a capability to evolve mechanical and thermal properties of constituent materials with irradiation time and burn-up were developed. The models input used projected parameters, including plate geometry, irradiation history and coolant conditions as input. The model output included temperature, displacement and stresses in the fuel, cladding and diffusion barrier. The fuel behavioral model considered inelastic behavior including volumetric swelling due to solid and gaseous products, irradiation induced creep, thermal expansion, conductivity degradation and plasticity. A visco-plastic behavioral model was used for the cladding that included thermal creep, irradiation hardening, growth due to fast neutrons and Mises plasticity. The plates were then simulated by using projected irradiation parameters. The resulting temperature, displacement and stress-strains were comparatively evaluated on the selected paths. The results were then compared with those of plates from previous RERTR experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗