Search NASA⌕ Search

SEARCH · Search NASA

Results for “data flow statistics”

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 163 records · Page 9

Particle Loading Dynamics of SRNL Electrostatic Particle Collectors - 19398

The Department of Energy is interested in decontamination methods and tools for the removal of airborne radioactive particulate from facilities prior to decommissioning. Electrostatic particle collectors (EPCs) have historically been used industrially to scrub particles from air streams for large scale operations, and Savannah River National Laboratory (SRNL) researchers have expanded on the fundamental collection capabilities of EPCs to develop unique samplers focused on particulate effluent collection. One distinct advantage of an electrostatic particle collection over filtration or inertial impaction-based methods is the low pressure drop associated with prolonged use of EPCs. Electrostatic particle collectors previously developed by SRNL have proven effective in the collection of atmospheric aerosols in multiple field trials, and their applicability as a decontamination method for airborne particulates. Initial experiments using aerosolized NaCl were performed to determine EPC collection efficiencies, the times required to achieve a desired particle loading, and the size distribution of collected particulates. The amount of NaCl collected on Si substrates was found to increase linearly with time (∼ 2 mg per hour) and showed little sign of change after 24 hours of collection. Experiments showed that NaCl, or a similar aerosol, would likely continue to collect unperturbed until physical limitations within the EPC are reached. Collection efficiencies (i.e., the number of particles collected vs. the number of particles entering the EPC) were calculated independently using an assessment scheme based on Particle Depletion Dynamics (PDD) and by Computer Controlled Scanning Electron Microscope (CCSEM) analysis of collection substrates for a range of particle sizes. For the ACE 2.0 with a corona voltage of 15 kV, collection efficiencies were determined using PDD to be 54, 29, and 13% for the flow rates 100, 300, and 490 L/min respectively. While CCSEM data suggests a very high collection efficiency for the MiniACE, the static sampling rate of 28.1 L/min (1 ft{sup 3}/min) was found to be too low volume for PDD to provide statistically meaningful data. Number density maps of NaCl particles detected via CCSEM revealed that particles are preferentially collected near the EPC inlet (i.e., opposite the direction of airflow) when using slower flow rates. Particles collected using higher flow rates (e.g., 300 and 490 L/min) were collected more evenly across the substrate. When segregated by size, smaller (e.g., 200-500 nm) particles are collected evenly across the collection substrates. As particle sizes increase they are collected preferentially near the EPC inlet, but not in sufficient quantity to obstruct airflow and hinder long-term operations. Subsequent deployment of the SRNL developed system into plutonium-238 contaminated hot cells will be performed at Savannah River Site Building 235-F. Presented in this work are the testing and considerations utilized by SRNL to show this technique as a viable means of rapidly decontaminating airborne contaminated areas. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Generating Emissions Inventory for Carbon Capture and Storage Analysis for Carbon-Intensive Industrial Sectors

Decarbonizing the industrial sector is critical to achieve carbon dioxide (CO2) emissions reductions goals of the Biden Administration. Currently available decarbonization options include electrification, fuel switching to zero carbon fuels like green hydrogen (H2) and carbon capture and storage (CCS). Application of post-combustion carbon capture (PCCC) technology in the power sector, as well as research at the U.S. Department of Energy's Fossil Energy and Carbon Management (FECM) Office has shown that its application in the industrial sector could have co-benefits in the form of emissions reductions of non-CO2 regulated pollutants. For example, solvent based PCCC systems require pre-conditioning of flue gas to remove sulfur and particulate matter (PM) upstream of the CO2 absorber. However, there is a lack of understanding about the type of non-CO2 pollutants which can be captured and the amount of reduction possible. PCCC application differs across industrial sectors as it depends on the availability of decarbonization options, characteristics of industrial processes and the amount and composition of pollutant flows. Certain facilities can also have multiple effluent flows with or without a CO2 stream. As such, understanding industrial processes and their effluent flows in detail is required to quantify the co-benefits opportunities presented by PCCC. Considering this requirement, the goal of this analysis is to develop a high-resolution inventory of effluent flows from facilities of 8 industrial sectors in the U.S. These industrial sectors - ethanol, ammonia, cement, steel, natural gas processing, hydrogen, petroleum refining and wood and pulp products - have carbon-intensive effluent flows, and thus are prime candidates for PCCC applications. In this study, we map the composition of pollutant flow from flue stacks across the identified facilities. Using data available in three Environmental Protection Agency (EPA) databases - the Green House Gas Reporting Program (GHGRP), the National Emissions Inventory (NEI) and the Toxic Release Inventory (TRI), we create a combined inventory which lists the type, amount, and concentration of pollutant flows. Using total weight of the pollutant flow back calculated from observed data for CO2 concentrations in flue gas for individual sectors, we calculate the concentration of each pollutant in the flue gas stream. Thus, the resultant emissions inventory includes the following details for each facility in the sector: facility-level and if possible, process-level pollutant flows, concentrations of pollutants in the flue gas, and geographical coordinates of the facilities. A detailed statistical analysis and summary allows us to search for erroneous data and remove them from the final inventory. The generation of the inventory is achieved using a python-based framework which can recreate this inventory for other industrial sectors as well as using newer releases of emission inventories from EPA. The statistical analysis performed on the inventory is also calibrated and automated to identify outliers efficiently.

air pollutants↗

Data-driven Feature Analysis in Control Design for Series-Compensated Transmission Systems

One challenge in power-system control designs is the gap between numerical model-based analysis and complex real-world power systems. With increased data and measurements being collected from power systems, data-driven analysis (e.g., machine learning) may provide an alternative approach to reveal hidden information through learning from the real system data, and provide insights for better control scheme design during the utility planning process. Data-driven feature analysis was performed to evaluate the relationships between series compensation, power generation, and path flows in a real transmission system, as well as temporal patterns. The main data-driven analysis methods, including statistical cross-correlation, multinomial logistical regression, and classification and regression trees, were integrated for feature selection and developing predictive models of series compensation. Analysis results demonstrated the effectiveness of the proposed methodology in feature analysis and the potential to help improve power-system control scheme design.

42 ENGINEERING↗

A perspective of Middle-Atmosphere Dynamics (MAD) studies at the New International Equatorial Observatory (NIEO)

The equatorial region has attracted many MAD studies mainly based on data of limited locations and resolutions. Established at NIEO are: (1) Climatology of the equatorial middle atmosphere (all of the mean zonal flow, the meridional and/or east-west circulations and the planetary/gravity waves are described based on massive, reliable data statistics); (2) Troposphere-stratosphere coupling at the equator (the candidate location of NIEO is just at the stratospheric fountain area where the tracers and waves are pumped up into the middle atmosphere); and (3) Mesosphere-thermosphere coupling at the equator; thermospheric superrotation, which may be caused either by ion drag or by tidal breaking, is examined in detail by observations covering a wide altitude range from the mesosphere through the thermosphere.

Yamanaka, M. D.↗

Large-Eddy Simulation of Flow Over Boeing Gaussian Bump Using Multiagent Reinforcement Learning Wall Model: Preprint

We develop a wall model for large-eddy simulation (LES) that takes into account various pressure-gradient effects using multi-agent reinforcement learning. The model is trained using low-Reynolds-number flow over periodic hills with agents distributed on the wall at various computational grid points. It utilizes a wall eddy-viscosity formulation as the boundary condition to apply the modeled wall shear stress. Each agent receives states based on local instantaneous flow quantities at an off-wall location, computes a reward based on the estimated wall-shear stress, and provides an action to update the wall eddy viscosity at each time step. The trained wall model is validated in wall-modeled LES of flow over periodic hills at higher Reynolds numbers, and the results show the effectiveness of the model on flow with pressure gradients. The analysis of the trained model indicates that the model is capable of distinguishing between the various pressure gradient regimes present in the flow. To further assess the robustness of the developed wall model, simulations of flow over the Boeing Gaussian bump are conducted at a Reynolds number of 2 x 10^6, based on the free-stream velocity and the bump width. The results of mean skin friction and pressure on the bump surface, as well as the velocity statistics of the flow field, are compared to those obtained from equilibrium wall model (EQWM) simulations and published experimental data sets. The developed wall model is found to successfully capture the acceleration and deceleration of the turbulent boundary layer on the bump surface, providing better predictions of skin friction near the bump peak and exhibiting comparable performance to the EQWM with respect to the wall pressure and velocity field. We also conclude that the subgrid-scale model is crucial to the accurate prediction of the flow field, in particular the prediction of separation.

boundary layer↗

Constraining the local variance of H {sub 0} from directional analyses

We evaluate the local variance of the Hubble Constant H {sub 0} with low-z Type Ia Supernovae (SNe). Our analyses are performed using a hemispherical comparison method in order to test whether taking the bulk flow motion into account can reconcile the measurement of the Hubble Constant H {sub 0} from standard candles ( H {sub 0} = 73.8±2.4 km s{sup -1} Mpc {sup -1}) with that of the Planck's Cosmic Microwave Background data ( H {sub 0} = 67.8 ± 0.9km s{sup -1} Mpc{sup -1}). We obtain that H {sub 0} ranges from 68.9±0.5 km s{sup -1} Mpc{sup -1} to 71.2±0.7 km s{sup -1} Mpc{sup -1} through the celestial sphere (1 σ uncertainty), implying a Hubble Constant maximal variance of δ H {sub 0} = (2.30±0.86) km s{sup -1} Mpc{sup -1} towards the ( l,b ) = (315°,27°) direction. Interestingly, this result agrees with the bulk flow direction estimates found in the literature, as well as previous evaluations of the H {sub 0} variance due to the presence of nearby inhomogeneities. We assess the statistical significance of this result with different prescriptions of Monte Carlo simulations, obtaining moderate statistical significance, i.e., 68.7% confidence level (CL) for such variance. Furthermore, we test the hypothesis of a higher H {sub 0} value in the presence of a bulk flow velocity dipole, finding some evidence for this result which, however, cannot be claimed to be significant due to the current large uncertainty in the SNe distance modulus. Then, we conclude that the tension between different H {sub 0} determinations can plausibly be caused to the bulk flow motion of the local Universe, even though the current incompleteness of the SNe data set, both in terms of celestial coverage and distance uncertainties, does not allow a high statistical significance for these results or a definitive conclusion about this issue.

79 ASTRONOMY AND ASTROPHYSICS↗

Plug20 Test Report: Particle Image Velocimetry (PIV) Measurements

This report provides turbulent flow statistics made using particle image velocimetry (PIV) for several internally mixed exhaust systems with external plugs (Plug20 designs). The tests were conducted at the NASA Glenn Research Center’s Aero-Acoustic Propulsion Laboratory (AAPL) in May and June 2022. The corresponding acoustic data for the configurations and flows presented in this report can be found in an earlier Technical Memorandum, NASA/TM-20210010291, “Plug20 Test Report.” The Plug20 nozzle family was designed to represent commercial supersonic vehicles in the 2025 to 2035 timeframe with a bypass ratio around 4, having an internal mixer, with and without an external plug to improve nacelle boattail drag and shock. The flow conditions measured represent engine cycles for commercial supersonic aircraft at takeoff, having pressure ratios of 2.0 and heated fan and core streams, measured both with and without a flight stream surrounding the plume. A few unheated flow conditions were also tested to tie results to other baseline test results. The PIV measurements in both streamwise and cross-stream planes are described, along with brief analysis to visually catalog the contents and confirm their validity. In this analysis, the streamwise PIV is found to match benchmark datasets within its experimental error. However, the cross-stream turbulent quantities measured by two systems do not match each other or the streamwise measurements in amplitude—caution is required when using the cross-stream data. Data acquired and included in this report are (1) velocity statistics of the plumes in both a streamwise and several cross-stream planes, and (2) the hot nozzle geometry. Tabular forms of the data and computer-aided design (CAD) files of the hot geometry are presented as supplemental files described in Appendix A (TM-20230006779SUPPL1.zip) and Appendix B (TM-20230006779SUPPL2.zip).

Turbulence↗

Calibration of Direct Velocimetry Using Hot Wire Anemometry

We studied the functionality of the IFA300 Hot Wire Anemometry system and its feasibility for the calibration of Direct Estimator Velocimetry (DEV) using Laser Induced Florescence (LIF). The experimental setup consisted of a single-wire hot wire probe measuring one component of velocity, the anemometry computational hardware and software, an HP oscilloscope, and a table fan to produce a simple flow with fluctuations. Measurements were taken at several points in the stream wise and transverse directions in the flow and various parameters recorded such as mean velocities temperature, turbulence intensities, skewness coefficients and flatness coefficients. The IFA300 software also allowed us to perform statistical manipulations such as spectrum analysis on velocities samples and correlation. Utilization of data files, also produced by the anemometry software, and post analysis were performed to produce graphical representations of turbulent intensity versus probe position and a flow field velocity profile. We concluded that the IFA300 Hot Wire Anemometry system is a reliable and functional method for calibration of DEV using LIF. Our future intentions are to set up a test chamber such that both velocity measurement techniques can be applied simultaneously, thus the calibration.

Holland, Michael L.↗

A hybrid point-particle force model that combines physical and data-driven approaches

Here this study improves upon the physics-based pairwise interaction extended point-particle (PIEP) model. The PIEP model leverages our physical understanding to predict fluid mediated interactions between solid particles. By considering the relative location of neighboring particles, the PIEP model is able to provide better predictions than existing drag models. While the current physical PIEP model is a powerful tool, its assumptions lead to increased error in flows with higher particle volume fractions. To reduce this error, a regression algorithm makes direct use of the results of direct numerical simulations (DNS) of an array of monodisperse solid particles subjected to uniform ambient flow at varying Reynolds numbers. The resulting statistical model and the physical PIEP model are superimposed to construct a hybrid, physics-based data-driven PIEP model. It must be noted that the performance of a pure data-driven approach without the model-form provided by the physical PIEP model is substantially inferior. The hybrid model's predictive capabilities are analyzed using additional DNS data that was not part of training the data-driven model. In every case tested, the hybrid models resulting from the regression were capable of (1) improving upon the physical PIEP model's prediction and (2) recovering underlying relevant physics from the DNS data. As the particle volume fraction increases, the physical PIEP model loses the ability to approximate the forces experienced by the particles, but the statistical model continues to produce accurate approximations.

97 MATHEMATICS AND COMPUTING↗

Data-driven learning of Mori–Zwanzig operators for isotropic turbulence

Developing reduced-order models for turbulent flows, which contain dynamics over a wide range of scales, is an extremely challenging problem. In statistical mechanics, the Mori–Zwanzig (MZ) formalism provides a mathematically exact procedure for constructing reduced-order representations of high-dimensional dynamical systems, where the effects due to the unresolved dynamics are captured in the memory kernel and orthogonal dynamics. Turbulence models based on MZ formalism have been scarce due to the limited knowledge of the MZ operators, which originates from the difficulty in deriving MZ kernels for complex nonlinear dynamical systems. In this work, we apply a recently developed data-driven learning algorithm, which is based on Koopman's description of dynamical systems and Mori's linear projection operator, on a set of fully resolved isotropic turbulence datasets to extract the Mori–Zwanzig operators. With data augmentation using known turbulence symmetries, the extracted Markov term, memory kernel, and orthogonal dynamics are statistically converged and the generalized fluctuation–dissipation relation can be verified. The properties of the memory kernel and orthogonal dynamics, and their dependence on the choices of observables are investigated to address the modeling assumptions that are commonly used in MZ-based models. A series of numerical experiments are then constructed using the extracted kernels to evaluate the memory effects on prediction. The results show that the prediction errors are strongly affected by the choice of observables and can be further reduced by including the past history of the observables in the memory kernel.

97 MATHEMATICS AND COMPUTING↗

Statistical distributions for transient transport

Here, this paper introduces the use of statistical distributions based on transport differential equations for clear distinction of transport modes within transient kinetic experiments. More specifically, novel techniques are developed for the transient data obtained through the Temporal Analysis of Products (TAP) reactor and are applicable to experiments where pulse response into a gas flow is used. The methodology allows distinguishing between two domains of diffusion transport in heterogeneous catalytic systems, i.e., Knudsen and non-Knudsen diffusion, using statistical fingerprints, and finding the transition domain. Two distribution parameters were obtained that directly result in coefficients that correspond to the concentration and the rate of transport. Using a linear relationship between the rate and concentration coefficients, Knudsen diffusion is revealed when the rate of transport is constant and non-Knudsen diffusion is confirmed when the rate of transport coefficient is a function of the concentration coefficient. As a result, accurate transport information can be extracted from experimental data even in the presence of comprising instrument drift or noise particularly when analyzing higher pressure pulse responses with complex transport. This enables more direct investigation of experiments influenced by gas-phase reactions.

TAP reactor↗

Mountain Basin Controls on the Snow-to-Streamflow Signal: An AIC-Weighted Multiple Linear Regression Framework

A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

Reduction of flow-measurement uncertainties in laser velocimeters with nonorthogonal channels

An analysis of certain geometrical limitations inherent in the application of laser velocimeters with nonorthogonal channels has led to the development of advanced-LDA-calibration and data-acquisition techniques that minimize systematic and statistical errors, respectively. The data-acquisition technique optimizes the number of velocity samples collected from three velocimeter channels as a function of local turbulence intensity, vector direction, and prescribed confidence interval. Linear velocity surveys and streamline traces measured in a turbulent flow field with a three-dimensional laser velocimeter are presented and the validity and accuracy of the theoretical analysis are discussed.

Snyder, P. K.↗

Benefits and Limits of Phasing Alleles for Network Inference of Allopolyploid Complexes

Abstract Accurately reconstructing the reticulate histories of polyploids remains a central challenge for understanding plant evolution. Although phylogenetic networks can provide insights into relationships among polyploid lineages, inferring networks may be hindered by the complexities of homology determination in polyploid taxa. We use simulations to show that phasing alleles from allopolyploid individuals can improve phylogenetic network inference under the multispecies coalescent by obtaining the true network with fewer loci compared with haplotype consensus sequences or sequences with heterozygous bases represented as ambiguity codes. Phased allelic data can also improve divergence time estimates for networks, which is helpful for evaluating allopolyploid speciation hypotheses and proposing mechanisms of speciation. To achieve these outcomes in empirical data, we present a novel pipeline that leverages a recently developed phasing algorithm to reliably phase alleles from polyploids. This pipeline is especially appropriate for target enrichment data, where the depth of coverage is typically high enough to phase entire loci. We provide an empirical example in the North American Dryopteris fern complex that demonstrates insights from phased data as well as the challenges of network inference. We establish that our pipeline (PATÉ: Phased Alleles from Target Enrichment data) is capable of recovering a high proportion of phased loci from both diploids and polyploids. These data may improve network estimates compared with using haplotype consensus assemblies by accurately inferring the direction of gene flow, but statistical nonidentifiability of phylogenetic networks poses a barrier to inferring the evolutionary history of reticulate complexes.

Evolutionary Biology↗

Dimensionless Regime Mapping and Chaos Analysis of Solids Volume Fraction Data in a Circulating Fluidized Bed Riser

Presentation of experimental results of solids volume fraction data measurements in a 12-inch diameter cold flow circulating fluidized bed in relation to a dimensionless fluidization regime map. Statistical and chaotic system analysis results are shown to illustrate differences in hydrodynamic behavior at different operating conditions. Presentation was made at the ASTFE 10th Thermal and Fluids Engineering Conference, March 09-12, 2025.

chaos↗

Surface topography of the Greenland Ice Sheet from satellite radar altimetry

Surface elevation maps of the southern half of the Greenland subcontinent are produced from radar altimeter data acquired by the Seasat satellite. A summary of the processing procedure and examples of return waveform data are given. The elevation data are used to generate a regular grid which is then computer contoured to provide an elevation contour map. Ancillary maps show the statistical quality of the elevation data and various characteristics of the surface. The elevation map is used to define ice flow directions and delineate the major drainage basins. Regular maps of the Jakobshavns Glacier drainage basin and the ice divide in the vicinity of Crete Station are presented. Altimeter derived elevations are compared with elevations measured both by satellite geoceivers and optical surveying.

Bindschadler, Robert A.↗

NASA Goddard Space Flight Center

The NASA SLR Operational Center is responsible for: 1) NASA SLR network control, sustaining engineering, and logistics; 2) ILRS mission operations; and 3) ILRS and NASA SLR data operations. NASA SLR network control and sustaining engineering tasks include technical support, daily system performance monitoring, system scheduling, operator training, station status reporting, system relocation, logistics and support of the ILRS Networks and Engineering Working Group. These activities ensure the NASA SLR systems are meeting ILRS and NASA mission support requirements. ILRS mission operations tasks include mission planning, mission analysis, mission coordination, development of mission support plans, and support of the ILRS Missions Working Group. These activities ensure than new mission and campaign requirements are coordinated with the ILRS. Global Normal Points (NP) data, NASA SLR FullRate (FR) data, and satellite predictions are managed as part of data operations. Part of this operation includes supporting the ILRS Data Formats and Procedures Working Group. Global NP data operations consist of receipt, format and data integrity verification, archiving and merging. This activity culminates in the daily electronic transmission of NP files to the CDDIS. Currently of all these functions are automated. However, to ensure the timely and accurate flow of data, regular monitoring and maintenance of the operational software systems, computer systems and computer networking are performed. Tracking statistics between the stations and the data centers are compared periodically to eliminate lost data. Future activities in this area include sub-daily (i.e., hourly) NP data management, more stringent data integrity tests, and automatic station notification of format and data integrity issues.

Carter, David↗

Partnership for Edge Physics Simulation

A major goal of our participation in the Edge Physics Simulation project has been to contribute to the understanding of the self-organization of tokamak turbulence fluctuations resulting in the formation of a staircase structure in the ion temperature. A second important goal is to demonstrate how small scale turbulence in plasmas self-organizes with dynamically driven quasi-stationary flow shear. These goals have been accomplished through the analyses of the statistical properties of XGC1 flux driven Gyrokinetic electrostatic ion temperature gradient (ITG) turbulence simulation data in which neutrals are included. The ITG turbulence data, and in particular fluctuation data, were obtained from a massively parallel flux-driven gyrokinetic full-f particle-in-cell simulation of a DIII-D like equilibrium. Below some the findings are summarized. It was observed that the emergence of staircase structure is related to the variations in the normalized temperature gradient length (R/LT) and the poloidal flow shear. Average turbulence intensity is found to be large in the vicinity of minima in R/LTi, where ITG growth is expected to be lower. The distributions of the occurrences of potential fluctuation are found to be Gaussian away from the staircase-step locations, but they are found to be non-Gaussian in the vicinity of staircase-step locations. The results of analytically derived expressions for the distribution of the occurrences of turbulence intensity and intensity flux were compared with the corresponding quantities computed using XGC1 simulation data and good agreement is found. The derived expressions predicts inward and outward propagation of turbulence intensity flux in an intermittent fashion. The outward propagation of turbulence intensity flux occurs at staircase-step locations and is related to the change in poloidal flow velocity shear and to the change in the ion temperature gradient. The standard deviation, skewness and kurtosis for turbulence quantities were computed and found to be large in the vicinity of the staircase-step structures. Large values of skewness and kurtosis can be explained by a temporary opening and closing of the structure which allows turbulence intensity events to propagate. The staircase patterns may reduce the ion heat transport and a manipulation of these patterns may be used to optimize heat transport in tokamaks. An additional objective of the research in support of the Edge Physics Simulation initiative has been to improve the understanding of scrape-off layer thermal transport. In planning experiments and designing future tokamaks, it is important to understand the physical effects that contribute to divertor heat-load fluxes. The research accomplished will contribute to developing new models for the scrape-off layer region. The XGC0 code was used to compute the heat fluxes and the heat-load width in the outer divertor plates of C-Mod and DIII-D tokamaks. It was observed that the width of the XGC0 neoclassical heat-load was approximately inversely proportional to the total plasma current. Anomalous transport in the H-mode pedestal region of five Alcator C-Mod discharges, representing a collisionality scan, was analyzed. The understanding of anomalous transport in the pedestal region is important for the development of a comprehensive model for the H-mode pedestal slope. It was found that the electron thermal anomalous diffusivities at the pedestal top increase with the electron collisionality. This dependence can point to the DRIBM as the modes that drive the anomalous transport in the plasma edge of highly collisional discharges. The effects of plasma shaping on the H-mode pedestal structure was also investigated. The differences in the predicted H-mode pedestal width and height for the DIII-D discharges with different elongation and triangularities were discussed. For the discharges with higher elongation, it was found that the gradients of the plasma profiles in the H-mode pedestal reach semi-steady states. In these simulations, the pedestal slowly continued to evolve to higher pedestal pressures and bootstrap currents until the peeling ballooning stability conditions were satisfied. The discharges with lower elongation do not reach the semi-steady state, and ELM crashes were triggered at earlier times. The plasma elongation was found to have a stronger stabilizing effect than the plasma triangularity. For the discharges with lower elongation and lower triangularity, the ELM frequency was large, and the H-mode pedestal evolves rapidly. It was found that the temperature of neutrals in the scrape-off-layer region can affect the dynamics of the H-mode pedestal buildup. However, the final pedestal profiles were nearly independent of the neutral temperature. The elongation and triangularity affected the pedestal widths of plasma density and electron temperature profiles differently. This study illustrated a new mechanism for controlling the pedestal bootstrap current and the pedestal stability.

97 MATHEMATICS AND COMPUTING↗