Search NASA⌕ Search

SEARCH · Search NASA

Results for “model structure uncertainty”

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 631 records · Page 35

Life Prediction Issues in Thermal/Environmental Barrier Coatings in Ceramic Matrix Composites

Issues and design requirements for the environmental barrier coating (EBC)/thermal barrier coating (TBC) life that are general and those specific to the NASA Ultra-Efficient Engine Technology (UEET) development program have been described. The current state and trend of the research, methods in vogue related to the failure analysis, and long-term behavior and life prediction of EBCITBC systems are reported. Also, the perceived failure mechanisms, variables, and related uncertainties governing the EBCITBC system life are summarized. A combined heat transfer and structural analysis approach based on the oxidation kinetics using the Arrhenius theory is proposed to develop a life prediction model for the EBC/TBC systems. Stochastic process-based reliability approach that includes the physical variables such as gas pressure, temperature, velocity, moisture content, crack density, oxygen content, etc., is suggested. Benefits of the reliability-based approach are also discussed in the report.

Shah, Ashwin R.↗

Understanding Model-Observation Discrepancies in Satellite Retrievals of Atmospheric Temperature Using GISS ModelE

We examine multiple factors in the representation of satellite-retrieved atmospheric temperature diagnostics in historical simulations of climate change during the satellite era (specifically 1979–2021) using GISS ModelE contributions to the Coupled Model Intercomparison Project (Phase 6) (CMIP6). The tropospheric and stratospheric trends in these diagnostics are affected by greenhouse gases (notably carbon dioxide and ozone), coupling with the ocean, volcanic aerosols, solar activity and compositional and dynamic feedbacks. We explore the impacts of internal variability, changing forcing specifications, composition interactivity, the quality of the stratospheric circulation, vertical resolution, and possible impacts of the mis-specification of volcanic aerosol optical depths. Overall temperature trends throughout the satellite period are well captured, but discrepancies at all levels exist and have multiple distinct causes. We find that stratospheric comparisons (using Stratospheric Sounding Unit (SSU) retrievals and successor instruments) are most affected by variations in the representation of ozone depletion and feedbacks, followed by the volcanic signals. Tropospheric skill (using the Microwave Sounding Unit (MSU) retrievals) is affected by the trends in ocean heat uptake and tropospheric aerosols, but also by the representation of stratospheric processes through the impact of the Brewer-Dobson circulation on the height of the tropical tropopause. We demonstrate that no single factor is the dominant cause of the discrepancies and that almost all observations lie within the broad envelope of structural uncertainty.

satellite retrievals of atmospheric temperature↗

Multiscale-Informed Modeling of High Temperature Component Response with Uncertainty Quantification

This report summarizes a joint effort between Argonne National Laboratory, Idaho National Laboratory, and Los Alamos National Laboratory to develop and deploy constitutive models targeted at predicting the life of Grade 91 alloy components subjected to high temperature environments typical of those that structural components in advanced nuclear reactors would experience. Two distinct, but complementary constitutive modeling approaches have been taken here. The first employs a phenomenological viscoplastic model for which parameters have been calibrated based on experimental data for a wide range of Grade 91 alloy that has undergone a variety of processing. A Bayesian approach was used to derive distributions of uncertain parameters for this model based on this data set. The second approach is a reduced order model suitable for engineering-scale analysis that is based on the results of a large set of mesoscale simulations. Mesoscale models allow for the microstructure and composition of a particular alloy to be directly taken into account in the computation of the viscoplastic response, but are computationally expensive, which makes it impractical to directly call those models for the material constitutive response in an engineering-scale simulation. The reduced-order representation of the response of the underlying model used here allows for an engineering-scale model to take into account the characteristics of the underlying microstructure, while only incurring a reasonable computational expense. Both of these approaches have different strengths, and are applicable for different parts of the design/analysis process. The phenomenological models can be readily parameterized based on a set of experimental data for a given class of materials and used for scoping calculations. Once a specific material is chosen and adequately characterized, the reduced order models can accurately predict the response of that specific alloy, and because the models are based on predictive models of the underlying microstructure, they can be used to more confidently predict the response under conditions in regions where there is limited experimental data. Both of these models have been integrated in the Grizzly code, which is used here to perform proof-of-concept uncertainty quantification analyses of a simple component under prototypical conditions. The built- in stochastic analysis capabilities in the MOOSE framework that Grizzly is built on are used here to run large sets of simulations for this uncertainty quantification analysis. As would be expected, because the reduced order models are developed for a much more tightly defined alloy, they predict tighter distributions of the time to failure than the phenomenological models, which are calibrated to a broader set of data. Also important is that these simulations demonstrate that a reduced order modeling approach can be successfully deployed to propagate uncertainties from the material scale to practical engineering-scale component simulations.

42 ENGINEERING↗

Using High Frequency Passive Microwave, A-train, and TRMM Data to Evaluate Hydrometer Structure in the NASA GEOS-5 Data Assimilation System

Validating water vapor and prognostic condensate in global models remains a challenging research task. Model parameterizations are still subject to a large number of tunable parameters; furthermore, accurate and representative in situ observations are very sparse, and satellite observations historically have significant quantitative uncertainties. Progress on improving cloud / hydrometeor fields in models stands to benefit greatly from the growing inventory ofA-Train data sets. ill the present study we are using a variety of complementary satellite retrievals of hydrometeors to examine condensate produced by the emerging NASA Modem Era Retrospective Analysis for Research and Applications, MERRA, and its associated atmospheric general circulation model GEOS5. Cloud and precipitation are generated by both grid-scale prognostic equations and by the Relaxed Arakawa-Schubert (RAS) diagnostic convective parameterization. The high frequency channels (89 to 183.3 GHz) from AMSU-B and MRS on NOAA polar orbiting satellites are being used to evaluate the climatology and variability of precipitating ice from tropical convective anvils. Vertical hydrometeor structure from the Tropical Rainfall Measuring Mission (TRMM) and CloudSat radars are used to develop statistics on vertical hydrometeor structure in order to better interpret the extensive high frequency passive microwave climatology. Cloud liquid and ice water path data retrieved from the Moderate Resolution Imaging Spectroradiometer, MODIS, are used to investigate relationships between upper level cloudiness and tropical deep convective anvils. Together these data are used to evaluate cloud / ice water path, gross aspects of vertical hydrometeor structure, and the relationship between cloud extent and surface precipitation that the MERRA reanalysis must capture.

Robertson, Franklin↗

Goal-oriented real-time Bayesian inference for linear autonomous dynamical systems with application to digital twins for tsunami early warning

We present a goal-oriented framework for constructing digital twins with the following properties: (1) they employ discretizations of high-fidelity partial differential equation (PDE) models governed by autonomous dynamical systems, leading to large-scale forward problems; (2) they solve a linear inverse problem to assimilate observational data to infer uncertain model components followed by a forward prediction of the evolving dynamics; and (3) the entire end-to-end, data-to-inference-to-prediction computation is carried out without approximation and in real time through a Bayesian framework that rigorously accounts for uncertainties. Several challenges must be overcome to realize this framework, including the large scale of the forward problem, the high dimensionality of the parameter space, and for a class of problems including those we target, the slow decay of the singular values of the parameter-to-observable map. Here we introduce a methodology to overcome these challenges by exploiting the autonomous structure of the forward model to decompose the solution of the inverse problem into a one-time-only offline phase in which the PDE model is solved a limited number of times (equal to the number of sensors), and an online phase that maps well onto GPUs and computes the parameter inference and prediction of quantities of interest in real time, given observational data. Our ultimate goal is to apply this framework to construct digital twins for subduction zones, including Cascadia, to provide early warning for tsunamis generated by megathrust earthquakes. To this end, we demonstrate how our methodology can be used to employ seafloor pressure observations, along with the coupled acoustic–gravity wave equations, to infer the earthquake-induced spatiotemporal seafloor motion (discretized with $\mathscr{O}$ (10 9 ) parameters) and forward predict the tsunami propagation. We present results of an end-to-end inference, prediction, and uncertainty quantification for a representative test problem with $\mathscr{O}$ (10 8 ) inversion parameters for which goal-oriented Bayesian inference is accomplished exactly and in real time, that is, in a matter of seconds.

97 MATHEMATICS AND COMPUTING↗

Wall Modeled Large Eddy Simulations for NASA’s Jet Noise Consensus Database of Single-Flow, Round, Convergent Jets

A campaign of wall-modeled large-eddy simulations (WMLES) using structured curvilinear overlapping grids has been performed with the Launch Ascent and Vehicle Aerodynamics(LAVA) computational fluid dynamics (CFD) software to predict jet noise for single-stream axisymmetric round jets. The simulations address the new Prediction Uncertainty Reduction(PUR) technical challenge within the context of NASA’s Commercial Supersonic Technology(CST) project. The goal of PUR is to quantify and reduce uncertainties from scale-resolving simulations to assess noise characteristics of next generation quiet supersonic commercial jets during takeoff and landing conditions where the noise from the exhaust jet dominates. The focus of this effort is to generate a simulation database for single-stream axisymmetric round nozzles at several conditions both for static (no ambient co-flow), which is the focus of this article, and in-flight (flight stream co-flow) conditions, which are beyond the current scope. Nine different flow conditions ranging in jet exit Mach number from 0.38 to 1.0 with nozzle temperature ratios (NTR) from 0.84 to 2.7 have been conducted. Details of the structured overset grids, numerical discretization and wall-model are provided. Near-field comparisons to PIV data show great agreement for both velocity and normal stresses, however a systematic TKE overshoot at the nozzle exit is seen in the lip line shear-layer. A permeable Ffowcs Williams Hawkings (FWH) surface, enclosing the jet, is used to predict far-field noise from the simulated flow-field. Comparison of CFD predictions to microphone array measurements demonstrate excellent agreement within the resolved frequency range. A systematic under-prediction of far-aft observer angles larger than 150 degrees has been observed across all simulations. We achieved a cost reduction of an order of magnitude for these WMLES compared to an earlier study of this configuration due to algorithmic and software improvements. The accuracy of the results and short turnaround time demonstrate that WMLES within the LAVA framework is a cost-effective approach for jet noise predictions that could soon be incorporated into the design cycle of jet noise reduction technologies.

CST↗

Bayesian Structural Time Series for Behind-the-Meter Photovoltaic Disaggregation: Preprint

Distributed photovoltaic (PV) generation often occurs ``behind the meter": a grid operator can only observe the net load, which is the sum of the gross load and distributed PV generation. This lack of observability poses a challenge to system operation at both bulk level and distribution level. The lack of real-time or near-future disaggregated estimates of gross load and PV generation will lead to over scheduling of energy production and regulation reserves, reliability constraints violations, wear and tear of controller devices, and potentially cascading failures of a system. In this paper we propose the use of a Bayesian Structural Time Series (BSTS) model with local solar irradiance measurements to disaggregate the summed PV generation and gross load signals at a downstream measurement site. BSTSs are a highly expressive model class that blends classic time series models with the powerful Bayesian state space estimation framework. Disaggregation is done probabilistically, which automatically quantifies the uncertainties of the estimated PV generation and gross load consumption. Depending on the data availability in real-time, it can be used to disaggragate PV and gross load at customer site, or can be used at the feeder level. In this paper, we focus on solving the problem at feeder level. We compare the performance of a BSTS model as well as a handful of state-of-the-art methods on a Pecan Street AMI dataset, using the National Solar Radiation Database (NSRDB) to estimate local irradiance.

Bayesian structural time series↗

Quantification of spatial and seasonal trends in the atmosphere and construction of statistical models for infrasonic propagation

SUMMARY Infrasonic waves are influenced by variations in the density, pressure and temperature as well as the ambient winds. Modelling infrasonic propagation can be challenging due to the dynamic nature of the atmosphere as well as the sparseness of measurements which result in variability and notable uncertainty. A framework is presented to quantify spatial and seasonal trends in atmospheric structure via analysis of the effective sound speed profile and identification of temporal trends in the middle atmospheric waveguide produced by the circumpolar vortex winds. Seasonal definitions identifying typical atmospheric structures during the summer, winter and spring/fall transition periods are identified using atmospheric data from 2010 through 2020. Seasonal trend analysis is conducted for a number of locations across the contiguous United States to quantify spatial variations in atmospheric structure that impact infrasonic propagation. For each season and location, empirical orthogonal function analysis is used to reduce the historical archive of atmospheric data into a smaller representative set that can be analysed using numerical tools more efficiently. Infrasonic ray tracing and finite-frequency modal propagation analyses are applied to construct propagation path geometry and transmission loss statistics which are useful in localization and yield estimation for infrasonic sources, respectively. An example application is detailed in which transmission loss statistics are combined with an explosive source model and noise statistics to quantify the capability of a network to detect nearby sources.

58 GEOSCIENCES↗

Tidal disruption of Periodic Comet Shoemaker-Levy 9 and a constraint on its mean density

The apparent tidal disruption of Periodic Comet Shoemaker-Levy 9 (1993e) during a close encounter within approximately 1.62 planetary radii of Jupiter can be used along with theoretical models of tidal disruption to place an upper bound on the density of the predisruption body. Depending on the theoretical model used, these upper bounds range from rho(sub c) less than 0.702 +/- 0.080 g/cu cm for a simple analytical model calibrated by numerical smoothed particle hydrodynamics (SPH) simulations to rho(sub c) less than 1.50 +/- 0.17 g/cu cm for a detailed semianalytical model. The quoted uncertainties stem from an assumed uncertainty in the perijove radius. However, the uncertainty introduced by the different theoretical models is the major source of error; this uncertainty could be eliminated by future SPH simulations specialized to cometary disruptions, including the effects of initially prolate, spinning comets. If the SPH-based upper bound turns out to be most appropriate, it would be consistent with the predisruption body being a comet with a relatively low density and porous structure, as has been asserted previously based on observations of cometary outgassing. Regardless of which upper bound is preferable, the models all agree that the predisruption body could not have been a relatively high-density body, such as an asteroid with rho approximately = 2 g/cu cm.

Boss, Alan P.↗

Use of the Long Duration Exposure Facility's thermal measurement system for the verification of thermal models

Results are presented of the comparison between the Long Duration Exposure Facility (LDEF) Thermal Measurements Systems (THERM) recorded temperature data and the predicted values as calculated prior to the LDEF deployment. The postflight thermal model was verified and calculated temperature uncertainties were reduced to under + or - 18 F from the preflight uncertainties of + or - 40 F. The THERM consisted of 8 temperature sensors, a shared tape recorder, a standard LDEF flight battery, and an electronics control box. The temperatures were measured at selected locations on the LDEF structure interior during the first 490 days of flight and recorded for postflight analysis. After the recorder was recovered from space, the tape recorder was recovered and the data reduced for comparison to the LDEF predicted temperatures. By comparing the calculated values to the measured data, a verified thermal model that presents the best agreement with the THERM data was obtained. The THERM experiment provided an economical way of performing a postflight verification of the LDEF Thermal Model by recording a limited number of flight temperatures on typical locations of the LDEF structure.

Berrios, William M.↗

Application of Stereo PIV on a Supersonic Parachute Model

The Mars Science Laboratory (MSL) is the next step in NASA's Mars Exploration Program, currently scheduled for 2011. The spacecraft's descent into the Martian atmosphere will be slowed from Mach 2 to subsonic speeds via a large parachute system with final landing under propulsive control. A Disk-Band-Gap (DBG) parachute will be used on MSL similar to the designs that have been used on previous missions, however; the DBG parachute used by MSL will be larger (21.5 m) than in any of the previous missions due to the weight of the payload and landing site requirements. The MSL parachute will also deploy at higher Mach number (M 2) than previous parachutes, which can lead to instabilities in canopy performance. Both the increased size of the DBG above previous demonstrated configurations and deployment at higher Mach numbers add uncertainty to the deployment, structural integrity and performance of the parachute. In order to verify the performance of the DBG on MSL, experimental testing, including acquisition of Stereo Particle Imaging Velocimetry (PIV) measurements were required for validating CFD predictions of the parachute performance. A rigid model of the DBG parachute was tested in the 10x10 foot wind tunnel at GRC. Prior to the MSL tests, a PIV system had never been used in the 10x10 wind tunnel. In this paper we discuss some of the technical challenges overcome in implementing a Stereo PIV system with a 750x400 mm field-of-view in the 10x10 wind tunnel facility and results from the MSL hardshell canopy tests.

Wernet, Mark P.↗

Dynamic Transmission Line Switching Amid Wildfire-Prone Weather Under Decision-Dependent Uncertainty

During dry and windy seasons, environmental conditions significantly increase the risk of wildfires, exposing power grids to disruptions caused by transmission line failures. Wildfire propagation exacerbates grid vulnerability, potentially leading to prolonged power outages. To address this challenge, we propose a multistage optimization model that dynamically adjusts transmission grid topology in response to wildfire propagation, aiming to develop an optimal response policy. By accounting for decision-dependent uncertainty, where line survival probabilities depend on usage, we employ distributionally robust optimization to model uncertainty in line survival distributions. We adapt the stochastic nested decomposition algorithm and derive a deterministic upper bound for its finite convergence. To enhance computational efficiency, we exploit the Lagrangian dual problem structure for a faster generation of Lagrangian cuts. Using realistic data from the California transmission grid, we demonstrate the superior performance of dynamic response policies against two-stage alternatives through a comprehensive case study. In addition, after solving the multistage formulation, we construct easy-to-implement policies that significantly reduce computational burden while maintaining good performance in real-time deployment. History: Accepted by Russell Bent, Area Editor for Network Optimization: Algorithms and Applications. Funding: This work was supported by the U.S. Department of Energy, Office of Electricity [Grant DE-AC02-05CH11231]. The work of R. Jiang was supported in part by the U.S. National Science Foundation, Division of Electrical, Communications and Cyber Systems [Grant ECCS-1845980] and the U.S. Air Force Office of Scientific Research [Grant FA9550-23-1-0323]. Supplemental Material: The software that supports the findings of this study is available within the paper and its Supplemental Information ( https://pubsonline.informs.org/doi/suppl/10.1287/ijoc.2025.1210 ) as well as from the IJOC GitHub software repository ( https://github.com/INFORMSJoC/2025.1210 ). The complete IJOC Software and Data Repository is available at https://informsjoc.github.io/ .

Estrada-Garcia, Juan-Alberto↗

Towards a Probabilistic Criterion for Preliminary Shell Design

Thin-walled stiffened or unstiffened, metallic or composite shells are widely used structural elements in aeronautical and space applications. Buckling strength, which is a major concern in all these areas, is affected by the uncertainties in the definition of loads, material properties, geometric variables, engineering models and the accuracy of the analysis tools used in the design phase. The NASA design criteria monographs from the late sixties account for these design uncertainties by the use of a lump sum safety factor or the so-called "knockdown" factor gamma, which usually results in an overly conservative design. In the present paper a new, reliability based, probabilistic design procedure for buckling critical imperfect isotropic shells is proposed. It essentially consists of a stochastic approach based on a new, improved "knockdown" factor lambda(sub a) that is not as conservative as the traditional one. It is felt that by quantifying and understanding the "problem uncertainties" such as initial imperfections and their influence on the design variables, one can develop a better engineered, better designed and safer system.

Arbocz, Johann↗

Kriging as a Means of Improving WAAS Availability

The Wide Area Augmentation System (WAAS), an augmentation of the Global Positioning System (GPS), provides safe and reliable use of GPS signals for airline navigation over much of North America. Currently the largest source of positioning error in the system is signal delay caused by the ionosphere. To allow the user to take account of such error, WAAS computes and broadcasts ionospheric vertical delays at a set of regularly-spaced grid points. In addition, WAAS computes and broadcasts a safety-critical integrity bound at each ionospheric grid point (IGP) called the Grid Ionospheric Vertical Error (GIVE). GIVEs are constructed to be sufficiently large to protect the user against positioning error due to the presence of ionospheric irregularity. In the initial operating capability (IOC) of WAAS, the vertical delay estimate at each IGP is determined from a planar fit of neighboring slant delay measurements, projected to vertical using an obliquity factor specified by the standard thin-shell model of the ionosphere. In WAAS Follow-On (WFO) Release 3, however, the vertical delay will be estimated by an established, geo-statistical technique known as kriging. Compared to the planar fit model, the kriging model is found, in general, to match better the observed random structure of the vertical delay. This paper presents the kriging methodology that will be used to estimate the vertical delay and its uncertainty at each IGP, and it assesses the subsequent improvement in WAAS availability enabled by kriging.

GPS↗

Near Earth Asteroid Scout Thrust and Torque Model

The Near Earth Asteroid (NEA) Scout is a solar sail mission whose objective is to scout at least one Near Earth Asteroid in preparation for manned missions to asteroids. NEA Scout will use a solar sail as the primary means of propulsion. Thus it is important for mission planning to accurately characterize the thrust of the sail. Additionally, the solar sail creates a relatively large solar disturbance torque that must be mitigated. For early mission design studies a flat plate model of the solar sail with a fixed center of pressure was adequate, but as mission concepts and the sail design matured, greater fidelity was required. Here we discuss the progress to a three-dimensional sail model that includes the effects of tension and thermal deformation that has been derived from a large structural Finite Element Model (FEM) developed by the Langley Research Center. We have found that the deformed sail membrane affects torque relatively much more than thrust. We have also found that other than uncertainty over the precise shape, the effect of small (approximately millimeter scale) wrinkles on the diffusivity of the sail is the leading remaining source of uncertainty. We demonstrate that millimeter-scale wrinkles can be modeled analytically as a change in the fraction of specular reflection. Finally we discuss the implications of these results for the NEA Scout mission.

Heaton, Andrew↗

Evaluation of PBR Spent Fuel Criticality and Dose Rate Compliance for Storage and Transportation

Spent tri-structural isotropic (TRISO)–based fuels have a strong track record in storage and transportation without documented incidents. This work seeks to reduce uncertainty to aid in more informed spent fuel management of TRISO-based fuels by modeling both fresh and spent pebble bed reactor (PBR) fuel and comparing the results to the regulatory standards from 10 CFR 71. SCALE was used for all modeling due to it having fast and accurate methods for handling PBR fuel modeling, as well as having an efficient method for shielding calculations in monaco with automated variance reduction using importance calculations (MAVRIC), which utilizes the consistent adjoint-driven importance sampling (CADIS) and the forward-weighted consistent adjoint-driven importance sampling (FW-CADIS) methods. KENO-VI was used for all criticality calculations, TSUNAMI was used for uncertainty quantification on k-effective, TRITON and the Oak Ridge isotope generation code (ORIGEN) were both used for depletion of the fuel, and MAVRIC was used for shielding calculations. For criticality assessments, this study focused on the requirement that the value of the neutron multiplication factor, k-effective (k-eff), would not exceed a peak value of 0.95, including uncertainty, with 95% confidence. Criticality was initially examined by modeling fresh fuel from three different designs—HTR-10 fuel, PBMR-400 fuel, and demonstration fuel representative of a TRISO-fueled modern high-temperature gas reactor (HTGR) design, henceforth referred to as Demo HTGR—and placing them into various sized containers with conditions described in 10 CFR 71 to quantify the peak k-eff state. When the peak value of 0.95 k-eff was exceeded, mitigation methods were examined in those scenarios. Burnup credit, pebble displacement in areas of strong neutron multiplication, and random pebble replacement using pebbles of various compositions and replacement fractions were examined. In summary, the criticality of PBR fuels can be well accounted for by restricting container size, taking credit for burnup, or by displacing/replacing pebbles. Uncertainty of the k-eff due to nuclear data uncertainties was recorded at ~0.6644%Δk/k, or roughly 664% mil (pcm). The nuclear data–induced uncertainty was relatively small and should not require significant modification in the design to be accounted for. Revisions to the evaluated nuclear data file values have been shown to have a larger impact than nuclear data–induced uncertainty. For dose rate aspects, U.S. Nuclear Regulatory Commission regulations require a maximum dose rate of 10 millirem per hour (mrem/h) at 2 meters. In examining the dose rate behavior of spent PBR fuel, the representative Demo HTGR fuel was modeled exclusively due to it possessing the highest target burnup of the examined fuels. Equilibrium cycle modeling methods were used to produce a higher-fidelity discharge isotopic composition than simple assumptions, such as reflected pebbles. The discharge composition was used as a source term in the fixed-source transport shielding calculations, and dose rates were calculated at 2 m for the shortest possible cooling time. The low concentration of fuel material led to dose rates that were in line with regulatory limits, despite the high burnup when compared to traditional light water reactor fuels. In conclusion, the methods employed in this study would require more work to further verify and validate and are limited to the criticality and dose rate analyses performed.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Prediction of DIII-D Pedestal Structure from Externally Controllable Parameters

The sharp increase of pressure at the edge of a high confinement mode (H-mode) plasma, the pedestal, strongly impacts overall plasma performance. Predicting the pedestal is a necessity to control and optimize tokamak operations. An experimental data-driven machine learning (ML) approach is presented that predicts the pedestal heights and widths of electron density (ne) and electron temperature (Te) profiles as well as the separatrix ne from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The OMFIT framework was used with DIII-D data to efficiently, robustly, and automatically build a database of pedestal parameters to train machine learning models. Database creation was enabled by the search engine tool for DIII-D data, TokSearch, which parallelizes data fetching, enabling fast searches through basic signals of thousands of DIII-D shots and selection of relevant time intervals. Principal Component Analysis (PCA) separated the database into three clusters that represent classes of plasma shapes that are regularly used in DIII-D. The most important parameters for setting the pedestal structure were plasma current (Ip), toroidal magnetic field (Bφ), neutral beam heating power (PNBI) and shaping quantities. The Deep Jointly Informed Neural Networks (DJINN) algorithm was applied to identify suitable neural network (NN) architectures that appropriately capture the features of the pedestal database. Separate NNs were implemented for each pedestal parameter, and ensembling methods were used to improve the prediction accuracy and allowed estimation of the prediction uncertainty. The pedestal predictions of the test dataset lie within the measurement uncertainties of the pedestal parameters. The NN outperformed simple Linear Regression (LR) analysis, indicating non-linear dependencies in the pedestal structure. The presented achievements illustrate a promising path for future research, using feature extraction to infer experimental trends and thereby improve pedestal models as well as deploying NN for a fast pedestal prediction in DIII-D scenario development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Prediction of DIII-D Pedestal Structure From Externally Controllable Parameters

The sharp increase of pressure at the edge of a high confinement mode (H-mode) plasma, the pedestal, strongly impacts overall plasma performance. Predicting the pedestal is a necessity to control and optimize tokamak operations. Here, an experimental data-driven machine learning (ML) approach is presented that predicts the pedestal heights and widths of electron density (n e ) and electron temperature (T e ) profiles as well as the separatrix ne from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The OMFIT framework was used with DIII-D data to efficiently, robustly, and automatically build a database of pedestal parameters to train machine learning models. Database creation was enabled by the search engine tool for DIII-D data, TokSearch, which parallelizes data fetching, enabling fast searches through basic signals of thousands of DIII-D shots and selection of relevant time intervals. Principal Component Analysis (PCA) separated the database into three clusters that represent classes of plasma shapes that are regularly used in DIII-D. The most important parameters for setting the pedestal structure were plasma current (I p ), toroidal magnetic field (B Φ ), neutral beam heating power (P NBI ) and shaping quantities. The Deep Jointly Informed Neural Networks (DJINN) algorithm was applied to identify suitable neural network (NN) architectures that appropriately capture the features of the pedestal database. Separate NNs were implemented for each pedestal parameter, and ensembling methods were used to improve the prediction accuracy and allowed estimation of the prediction uncertainty. The pedestal predictions of the test dataset lie within the measurement uncertainties of the pedestal parameters. The NN outperformed simple Linear Regression (LR) analysis, indicating non-linear dependencies in the pedestal structure. The presented achievements illustrate a promising path for future research, using feature extraction to infer experimental trends and thereby improve pedestal models as well as deploying NN for a fast pedestal prediction in DIII-D scenario development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗