Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable coefficients”

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 487 records · Page 27

TRW - Lunar Descent Engine

As we went through the program, what we determined, and what we all agreed on, was that the thrust coefficient (Cf) of the nozzle, after you get past a certain point, is really an engineering parameter. It s not a fundamental parameter that is going to be highly variable. Once we knew what the contour of the nozzle was, and once we knew what its characteristic was out to 2:1, we could calculate what the 48:1 thrust coefficient was going to be. In every case that we made a test, the calculation was precise. We weren't looking for a problem out at 48:1. Once we crushed the nozzle and said, "Yeah, we can land on the boulder," and once we had the thermal profile of that columbium nozzle, we did not require a lot of effort there. The real characterization was done in throttling over the 10:1 with the injector and controlling the mixture ratio on that - the whole head-end assembly - out to 2:1. I think everybody at NASA and Grumman agreed that flying like you test is great, particularly if you are using an aircraft engine. But, in this case, the thrust coefficient of the nozzle was not an issue. We had the tandem configuration of the service module, the command module, and the LEM sitting out there, and we were to fire the LEM. On Apollo 5, we were firing the LEM to show how it would work. There was a problem. I can t remember where the problem was, but something caused a problem before that engine had finished its burn. It was not in the engine, but there was some other problem, and NASA made a controlled shutdown. Then, they came to us and asked, "Hey, we re up there. We want to finish this test program. Is it okay if we restart that engine again in space with this tandem configuration?" We said, "As long as it has been more than forty minutes since you shut down, our analysis says that you will be okay in terms of the thermal characteristics of the inside of that chamber." They restarted it and pushed that system around in orbit on Apollo 5. It turned out, that when it came to Apollo 13, we went back into the record, and said, "Hey, we have pushed this system around up there on Apollo 5, and we have also restarted this tandem configuration." The requirements on Apollo 13 were to put it back into play. The spacecraft was out of free return to the earth at the time of the accident. It would not have come back. NASA said, "Okay, we ll use the descent engine to put the spacecraft in a free trajectory; it will go around the moon and be on free trajectory back to Earth." Then, as it came around the far side of the moon, the guys found out that they had an oxygen problem. As you remember, things were getting pretty bad in there. They said, "We ve got to get it back as fast as we can. Is it okay if we re-fire the engine? Now, we re in a free trajectory, so we want to put as much delta-v (or change in velocity) in as we can. Can we re-fire right now?" We said, "Yes, the data says it has been this period of time." We could re-fire the engine, run the rest of the duty cycle up as far as we needed while preserving enough fluids to make the final correction as the spacecraft got near Earth, and restart the engine. It was pretty fortuitous that we could give them those answers.

Elverum, Gerard W.↗

Transformations of the perturbed two-body problem to unperturbed harmonic oscillators

Singular, nonlinear, and Liapunov unstable equations are made regular and linear through transformations that change the perturbed planar problem of two bodies into unperturbed and undamped harmonic oscillators with constant coefficients, so that the stable solution may be immediately written in terms of the new variables. The use of arbitrary and special functions for the transformations allows the systematic discussion of previously introduced and novel anomalies. For the case of the unperturbed two-body problem, it is proved that if transformations are power functions of the radial variable, only the eccentric and the true anomalies (with the corresponding transformations of the radial variable) will result in harmonic oscillators. The present method significantly reduces computation requirements in autonomous space operations.

Szebehely, V.↗

Crystal-chemistry and partitioning of REE in whitlockite

Partitioning of Rare Earth Elements (REE) in whitlockite is complicated by the fact that two or more charge-balancing substitutions are involved and by the fact that concentrations of REE in natural whitlockites are sufficiently high such that simple partition coefficients are not expected to be constant even if mixing in the system is completely ideal. The present study combines preexisting REE partitioning data in whitlockites with new experiments in the same compositional system and at the same temperature (approximately 1030 C) to place additional constraints on the complex variations of REE partition coefficients and to test theoretical models for how REE partitioning should vary with REE concentration and other compositional variables. With this data set, and by combining crystallographic and thermochemical constraints with a SAS simultaneous-equation best-fitting routine, it is possible to infer answers to the following questions: what is the speciation on the individual sites Ca(B), Mg, and Ca(IIA) (where the ideal structural formula is Ca(B)18 Mg2Ca(IIA)2P14O56); how are REE's charge-balanced in the crystal; and is mixing of REE in whitlockite ideal or non-ideal. This understanding is necessary in order to extrapolate derived partition coefficients to other compositional systems and provides a broadened understanding of the crystal chemistry of whitlockite.

Colson, R. O.↗

Low Temperature Unbalanced Magnetron Deposition of Hard, Wear-Resistant Coatings for Liquid-Film Bearing Applications

The original program for evaluating the tribological properties several different hard coatings for liquid film bearing applications was curtailed when the time for the program was reduced from 3 years to 1. Of the several different coatings originally planned for evaluation, we decided to concentrate on one coating, carbon nitride. At BIRL, we have been instrumental in the development of reactively sputtered carbon nitride coatings, and we have found that it is a very interesting new material with very good tribological properties. In this program, we found that the reactively sputtered carbon nitride does not bond well directly to hardened 440C stainless steel; but if an interlayer of titanium nitride is added between the carbon nitride and the 440C, the adhesion of the dual coating combination is very good. Statistically designed experiments were run with the dual layer combination, and 3 variables were chosen for the Box-Benken design, which were the titanium nitride interlayer thickness, the nitrogen partial pressure during the reactive sputtering of the carbon nitride, and the carbon nitride substrate bias voltage. Two responses were studied from these three variables; the adhesion of the dual coating combination to the 440C substrate and the friction coefficient of the carbon nitride in dry sliding contact with 52100 steel in air. The best adhesion came with the thickness interlayer thickness studied, which was 4 micrometers, and the lowest coefficient of friction was 0.1, which was achieved when the bias voltage was in the range of -80 to - 120 V and the nitrogen partial pressure was 3 mTorr.

Sproul, William D.↗

Electric Propulsion Interactions Code (EPIC): Recent Enhancements and Goals for Future Capabilities

The Electric Propulsion Interactions Code (EPIC) is the leading interactive computer tool for assessing the effects of electric thruster plumes on spacecraft subsystems. EPIC, developed by SAIC under the sponsorship of the Space Environments and Effects (SEE) Program at the NASA Marshall Space Flight Center, has three primary modules. One is PlumeTool, which calculates plumes of electrostatic thrusters and Hall-effect thrusters by modeling the primary ion beam as well as elastic scattering and charge-exchange of beam ions with thruster-generated neutrals. ObjectToolkit is a 3-D object definition and spacecraft surface modeling tool developed for use with several SEE Program codes. The main EPIC interface integrates the thruster plume into the 3-D geometry of the spacecraft and calculates interactions and effects of the plume with the spacecraft. Effects modeled include erosion of surfaces due to sputtering, re-deposition of sputtered materials, surface heating, torque on the spacecraft, and changes in surface properties due to erosion and deposition. In support of Prometheus I (JIMO), a number of new capabilities and enhancements were made to existing EPIC models. Enhancements to EPIC include adding the ability to scale and view individual plume components, to import a neutral plume associated with a thruster (to model a grid erosion plume, for example), and to calculate the plume from new initial beam conditions. Unfortunately, changes in program direction have left a number of desired enhancements undone. Variable gridding over a surface and resputtering of deposited materials, including multiple bounces and sticking coefficients, would significantly enhance the erosion/deposition model. Other modifications such as improving the heating model and the PlumeTool neutral plume model, enabling time dependent surface interactions, and including EM1 and optical effects would enable EPIC to better serve the aerospace engineer and electric propulsion systems integrator. We review EPIC S overall capabilities and recent modifications, and discuss directions for future enhancements.

Gardner, Barbara M.↗

Unitary Plan Wind Tunnel Tests of Cook Technological Center Parachutes in the Wake of a Conical-Nosed Cylindrical Body Having a Base Diameter of 2.375-Inches (Part 5 of 6)

The film depicts two tests of a flat roof, conical inlet canopy parachute. The first test is a series of wind tunnel trials with a flat circular ribbon roof of 22 percent porosity. The second test is a single series of wind tunnel trials with a flat circular ribbon roof of 25 percent porosity. Variables for both trials include Mach number, dynamic pressure, longitudinal separation distances (x/d), and drag coefficient C(sub d).

Source record↗

Forecasting Ocean Chlorophyll in the Equatorial Pacific

Using a global ocean biogeochemical model combined with a forecast of physical oceanic and atmospheric variables from the NASA Global Modeling and Assimilation Office, we assess the skill of a chlorophyll concentrations forecast in the Equatorial Pacific for the period 2012-2015 with a focus on the forecast of the onset of the 2015 El Niño event. Using a series of retrospective 9-month hindcasts, we assess the uncertainties of the forecasted chlorophyll by comparing the monthly total chlorophyll concentration from the forecast with the corresponding monthly ocean chlorophyll data from the Suomi-National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite (S-NPP VIIRS) satellite. The forecast was able to reproduce the phasing of the variability in chlorophyll concentration in the Equatorial Pacific, including the beginning of the 2015-2016 El Niño. The anomaly correlation coefficient (ACC) was significant (p less than 0.05) for forecast at 1-month (R=0.33), 8-month (R=0.42) and 9-month (R=0.41) lead times. The root mean square error (RMSE) increased from 0.0399 microgram ch1 L(exp -1) for the 1-month lead forecast to a maximum of 0.0472 microgram ch1 L(exp -1) for the 9-month lead forecast indicating that the forecast of the amplitude of chlorophyll concentration variability was getting worse. Forecasts with a 3-month lead time were on average the closest to the S-NPP VIIRS data (23% or 0.033 microgram ch1 L(exp -1)) while the forecast with a 9-month lead time were the furthest (31% or 0.042 microgram ch1 L(exp -1)). These results indicate the potential for forecasting chlorophyll concentration in this region but also highlights various deficiencies and suggestions for improvements to the current biogeochemical forecasting system. This system provides an initial basis for future applications including the effects of El Niño events on fisheries and other ocean resources given improvements identified in the analysis of these results.

Anomaly Correlation Coefficient (ACC)↗

Performance evaluation of CMIP6 models on the Arctic-Siberian Plain teleconnection affecting the East Asian heat waves

The frequency and intensity of summer heat waves in East Asia have increased sharply in recent decades, significantly impacting public health and the economy. The Arctic-Siberian Plain (ASP) teleconnection pattern has been identified as a key driver, with ASP warming amplifying atmospheric circulation patterns conducive to extreme temperatures. This study evaluates the ability of Coupled Model Inter-comparison Project phase 6 models to simulate the ASP pattern across interannual variability (IAV) and intra-seasonal variability (ISV) timescales using the Common Basis Function method. The multi-model mean shows statistically significant pattern correlations with ERA5 reanalysis, with correlation coefficients of 0.90 and 0.99 for IAV and ISV, respectively. While the ASP pattern is generally well captured, models exhibit substantial inter-model diversity in the intensity and position of anticyclonic anomalies over the ASP and East Asia. Models with ASP pattern variability similar to reanalysis better reproduce extreme East Asian temperatures, whereas those over- or underestimating ASP variability exhibit lower skill. These performance differences are related to differences in simulating key variables associated with the development of the ASP pattern. Our findings highlight the role of the ASP pattern in modulating extreme heat events, as models with improved ASP simulations align more closely with observed temperature extremes. Refining ASP representations in models could enhance seasonal heat wave predictions, improving climate adaptation strategies.

Arctic-Siberian Plain (ASP)↗

NLSCIDNT user's guide maximum likehood parameter identification computer program with nonlinear rotorcraft model

A nonlinear, maximum likelihood, parameter identification computer program (NLSCIDNT) is described which evaluates rotorcraft stability and control coefficients from flight test data. The optimal estimates of the parameters (stability and control coefficients) are determined (identified) by minimizing the negative log likelihood cost function. The minimization technique is the Levenberg-Marquardt method, which behaves like the steepest descent method when it is far from the minimum and behaves like the modified Newton-Raphson method when it is nearer the minimum. Twenty-one states and 40 measurement variables are modeled, and any subset may be selected. States which are not integrated may be fixed at an input value, or time history data may be substituted for the state in the equations of motion. Any aerodynamic coefficient may be expressed as a nonlinear polynomial function of selected 'expansion variables'.

Source record↗

Sea surface height variability during 1986-1988 in the tropical Pacific Ocean

A description of the sea surface height variabilities during 1986-1988 in the tropical Pacific Ocean is presented on the basis of Geosat altimeter data and a tropical Pacific Ocean GCM (OGCM) forced with observed winds from Florida State University. Both data sets are validated against in situ sea-level measurements at 12 selected sea-level stations in the equatorial Pacific and were found to be in good agreement with these measurements with correlation coefficients of 0.69 for Geosat and 0.71 for OGCM. The Geosat altimeter data are capable of describing sea surface height variabilities during 1986-1988 in the tropical Pacific, in particular, the seasonal cycle, the 1987 El Nino and 1988 La Nina episodes, and the December 1986 equatorial Kelvin wave event. The OGCM-simulated surface dynamic height is found to be sensitive to different surface wind products.

Chao, YI↗

Interdecadal Variations in the Alaska Gyre

Climatic dynamic topography variations in the Alaska gyre during the period 1968-1990 are described with an objective analysis of more than 12000 STD and XBT stations, and COADS wind stress data Interannual the dynamic height and SST variations were correlated and were consistent with recently described large-scale climatic shifts in the North Pacific. The gyre was centered more to the east, circulation appeared stronger, and SST was lower during the early to mid-1970s than during the 1980s. The Aleutian low (NP and PNA indices) intensified during the interim, but the response did not appear as a gyre spinup. Instead, the associated wind stress anomalies forced a slowly varying dynamic height anomaly across the eastern and northern part of the gyre through Ekman convergence, which had the effect of displacing the gyre's low somewhat to the WSW in the 1980s. The wind curl spectrum was white, and the slow oceanic response was modeled as stochastic-forced climate variability with a simple first-order Markov autoregression process. Forcing was assumed to be Ekman pumping of the pycnocline, and the damping coefficient was estimated from the data to be approx. 1 yr. A hindcast with observed winds gave estimated dynamic height patterns similar to those observed, with a canonical correlation of 0.79 at 99% confidence. This response was weak in the western half of the gyre, where slow baroclinic variability may have been influenced by long Rossby wave propagation. A simple autoregression simulation using artificial white noise forcing shows the evolution of decadal variations similar in nature to those observed. This result, along with the low frequency correlation between dynamic height and SST, suggests that the upper-ocean climatic variability in this region is primarily wind forced.

Lagerloef, Gary S.E.↗

Stochastic Adaptive Droop Control in Frequency Regulation of Power Systems With Intermittent Generators

Modern power systems (MPSs), including microgrids (MGs), are increasingly incorporating multiple renewable energy sources (RESs) such as wind and solar power, as well as battery storage and controllable loads. While environmentally beneficial, these sources pose challenges for control and management due to their intermittent and stochastic nature, especially in maintaining frequency stability with multiple interconnected generators of varying capacities. Traditional droop control methods are effective in systems with generators that are dispatchable and have fixed generation capacities, but they fall short when applied to systems with RESs, where generation capacities are dynamic and affected by unpredictable environmental conditions. To address these challenges, this paper introduces a novel stochastic adaptive droop control (SADC) method for load frequency control (LFC). The proposed method adapts droop coefficients in real time, based on the measured stochastic data of power generation capacities, enabling more effective frequency regulation in systems with variable and intermittent power generation. Unlike traditional adaptive control methods, which assume constant or slowly-varying system parameters, this approach accounts for stochastic processes by modeling them as Markov chains, enabling robust performance under highly dynamic and unpredictable conditions. The key contributions of this work include the development of real-time droop coefficient adaptation algorithms, derivation of their stability and convergence properties, and the demonstration of the advantages of the method through simulations. Case studies highlight the improved performance of frequency regulation, particularly in addressing the impact of stochastic weather conditions and the benefits of reducing dependence on battery reserves in dealing with intermittency of RESs. Finally, this paper provides a comprehensive analysis of the theoretical foundations of the method, as well as practical implementation insights for future power systems with high penetration of RESs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Assessing the Skill of Chlorophyll Forecasts: Latest Development and Challenges Ahead Using the Case of the Equatorial Pacific

Using a global ocean biogeochemical model combined with a forecast of physical oceanic and atmospheric variables from the NASA Global Modeling and Assimilation Office, we assess the skill of a chlorophyll concentrations forecast in the Equatorial Pacific for the period 2012-2015 with a focus on the forecast of the onset of the 2015 El Nino event. Using a series of retrospective 9-month hindcasts, we assess the uncertainties of the forecasted chlorophyll by comparing the monthly total chlorophyll concentration from the forecast with the corresponding monthly ocean chlorophyll data from the Suomi-National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite (S-NPP VIIRS) satellite. The forecast was able to reproduce the phasing of the variability in chlorophyll concentration in the Equatorial Pacific, including the beginning of the 2015-2016 El Nino. The anomaly correlation coefficient (ACC) was significant (p less than 0.05) for forecast at 1-month (R=0.33), 8-month (R=0.42) and 9-month (R=0.41) lead times. The root mean square error (RMSE) increased from 0.0399 microgram chl L(exp -1) for the 1-month lead forecast to a maximum of 0.0472 microgram chl L(exp -1) for the 9-month lead forecast indicating that the forecast of the amplitude of chlorophyll concentration variability was getting worse. Forecasts with a 3-month lead time were on average the closest to the S-NPP VIIRS data (23% or 0.033 microgram chl L(exp -1)) while the forecast with a 9-month lead time were the furthest (31% or 0.042 microgram chl L(exp -1)). These results indicate the potential for forecasting chlorophyll concentration in this region but also highlights various deficiencies and suggestions for improvements to the current biogeochemical forecasting system. This system provides an initial basis for future applications including the effects of El Nino events on fisheries and other ocean resources given improvements identified in the analysis of these results.

El Nin↗

Evaluating Nighttime CALIOP 0.532 micron Aerosol Optical Depth and Extinction Coefficient Retrievals

NASA Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP) Version 3.01 5-km nighttime 0.532 micron aerosol optical depth (AOD) datasets from 2007 are screened, averaged and evaluated at 1 deg X 1 deg resolution versus corresponding/co-incident 0.550 micron AOD derived using the US Navy Aerosol Analysis and Prediction System (NAAPS), featuring two-dimensional variational assimilation of quality-assured NASA Moderate Resolution Imaging Spectroradiometer (MODIS) and Multi-angle Imaging Spectroradiometer (MISR) AOD. In the absence of sunlight, since passive radiometric AOD retrievals rely overwhelmingly on scattered radiances, the model represents one of the few practical global estimates available from which to attempt such a validation. Daytime comparisons, though, provide useful context. Regional-mean CALIOP vertical profiles of night/day 0.532 micron extinction coefficient are compared with 0.523/0.532 micron ground-based lidar measurements to investigate representativeness and diurnal variability. In this analysis, mean nighttime CALIOP AOD are mostly lower than daytime (0.121 vs. 0.126 for all aggregated data points, and 0.099 vs. 0.102 when averaged globally per normalised 1 deg. X 1 deg. bin), though the relationship is reversed over land and coastal regions when the data are averaged per normalised bin (0.134/0.108 vs. 0140/0.112, respectively). Offsets assessed within single bins alone approach +/- 20 %. CALIOP AOD, both day and night, are higher than NAAPS over land (0.137 vs. 0.124) and equal over water (0.082 vs. 0.083) when averaged globally per normalised bin. However, for all data points inclusive, NAAPS exceeds CALIOP over land, coast and ocean, both day and night. Again, differences assessed within single bins approach 50% in extreme cases. Correlation between CALIOP and NAAPS AOD is comparable during both day and night. Higher correlation is found nearest the equator, both as a function of sample size and relative signal magnitudes inherent at these latitudes. Root mean square deviation between CALIOP and NAAPS varies between 0.1 and 0.3 globally during both day/night. Averaging of CALIOP along-track AOD data points within a single NAAPS grid bin improves correlation and RMSD, though day/night and land/ocean biases persist and are believed systematic. Vertical profiles of extinction coefficient derived in the Caribbean compare well with ground-based lidar observations, though potentially anomalous selection of a priori lidar ratios for CALIOP retrievals is likely inducing some discrepancies. Mean effective aerosol layer top heights are stable between day and night, indicating consistent layer-identification diurnally, which is noteworthy considering the potential limiting effects of ambient solar noise during day.

aerosol optical depth↗

Exploring Environmental and Aerosol Impacts on Maritime Tropical Convection using Airborne Radiometer, Radar, Lidar, and Dropsondes

The field deployment phase of NASA’s Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP2Ex) took place around the Philippines during August–October 2019, with the primary goal of jointly investigating aerosols and tropical meteorology at the micro-β to meso-β scales. A suite of instruments was deployed on NASA’s P-3B Orion aircraft to accomplish this mission, including: the Advanced Microwave Precipitation Radiometer (AMPR), Airborne Precipitation and cloud Radar 3rd Generation (APR-3), High Spectral Resolution Lidar 2 (HSRL2), and Advanced Vertical Atmospheric Profiling System (AVAPS) dropsondes. To examine potential aerosol influences on maritime tropical convection, P-3 flight segments throughout CAMP2Ex were binned into similar environmental groups using “low,” “medium,” and “high” values of nine AVAPS-derived parameters with known physical connections to convective frequency and/or intensity. Aerosol concentrations in each flight segment were evaluated using three HSRL2 variables at 355 and 532 nm: aerosol backscatter, aerosol extinction, and aerosol optical thickness. A set of radiometer- and radar-derived variables directly related to convective frequency and/or intensity was used to characterize convection, which included: AMPR-derived integrated cloud liquid water; polarization-corrected temperatures at 10.7, 19.35, 37.1, and 85.5 GHz; peak equivalent radar reflectivity factor (ZH); peak height of 30-dBZ ZH; and the number of radar data columns with composite ZH > 30 dBZ. The ZH analyses were performed using both Ku- and Ka-band APR-3 data. For each flight segment, correlation coefficients were calculated between each convective parameter and aerosol concentrations within the “low,” “medium,” and “high” groups for each environmental variable. The environmental stratification thresholds were then varied in a series of sensitivity tests. Several noteworthy correlations were observed between the convective and aerosol parameters within the environmental groups, especially when stratifying the environments based on their 850–500-hPa temperature lapse rate, 700–500-hPa temperature lapse rate, and K-Index. The convective parameters were often correlated most strongly with 355-nm extinction, 532-nm extinction, and 532-nm backscatter, while the presence of precipitation-sized liquid and ice hydrometeors contributed to some unexpected negative correlations. In general, as environmental conditions became more favorable for convection, a stronger correlation was observed between the convective parameters and aerosol concentrations. However, a “Goldilocks” zone of medium aerosol concentration was correlated most strongly with the convective parameters in many cases. These results stress the importance of considering environmental and aerosol conditions together when evaluating their impacts on convection. This presentation will provide a detailed discussion of these results and their implications, a summary of the limitations associated with such an analysis, and suggestions for future work.

Corey G. Amiot↗

Kinetics of the Reaction C2H + O2 from 193 to 350 K Using Laser Flash Kinetic Infrared Absorption Spectroscopy

Rate coefficients for the reaction C2H + O2 are measured from 193 to 350 K by using transient absorption spectroscopy with an infrared color center laser. Ethynyl radicals are produced by pulsed laser photolysis of C2H2 in a temperature variable flow cell and a tunable color center laser probes the transient removal Of C2H ((bar-X(sup 2))Sigma(+))(0,0,0)) in absorption. The rate coefficient has a slight negative temperature dependence over the range 193-350 K which can be expressed as k(O2) = (1.5 +/- 0.3) x 10(exp -11) exp((230 +/- 36)/T) Cu cm/(Molecule.s). The lack of a pressure dependence under our experimental conditions strongly suggests that the HCCOO* complex is short lived with respect to encountering a collision and rapidly dissociates to products without competition from collisional stabilization to HCCOO or collisionally induced dissociation to reactants.

Opansky, Brian J.↗

Gas Transfer Across Air‐Water Interfaces in Inland Waters: From Micro‐Eddies to Super‐Statistics

In inland water covering lakes, reservoirs, and ponds, the gas exchange of slightly soluble gases such as carbon dioxide, dimethyl sulfide, methane, or oxygen across a clean and nearly flat air‐water interface is routinely described using a water‐side mean gas transfer velocity $\overline{k_{L}}$, where overline indicates time or ensemble averaging. The micro‐eddy surface renewal model predicts $\overline{k_{L}}$ = α o Sc -1/2 ($v\bar{ϵ}$) 1/4 , where Sc is the molecular Schmidt number, $v$ is the water kinematic viscosity, and $\bar{ϵ}$ is the waterside mean turbulent kinetic energy dissipation rate at or near the interface. While α o = 0.39 - 0.46 has been reported across a number of data sets, others report large scatter or variability around this value range. It is shown here that this scatter can be partly explained by high temporal variability in instantaneous ϵ around $\bar{ϵ}$, a mechanism that was not previously considered. As the coefficient of variation (CV e ) in ϵ increases, α o must be adjusted by a multiplier (1 = CV e 2 ) -3/32 that was derived from a log‐normal model for the probability density function of ϵ. Reported variations in α o with a macro‐scale Reynolds number can also be partly attributed to intermittency effects in ϵ. Such intermittency is characterized by the long‐range (i.e., power‐law decay) spatial auto‐correlation function of ϵ. That α o varies with a macro‐scale Reynolds number does not necessarily violate the micro‐eddy model. Instead, it points to a coordination between the macro‐ and micro‐scales arising from the transfer of energy across scales in the energy cascade.

Batchelor scale↗

Characterization of Arsenic and Selenium in Coal Fly Ash to Improve Evaluations for Disposal and Reuse Potential (Final Technical Report)

Coal fly ash is a high volume waste material that is discarded in landfills and surface water impoundments across the U.S. and is also widely recycled for a variety of applications. The leaching of potential of contaminants of concern, such as arsenic (As) and selenium (Se), is often the driver of risk assessments for coal ash disposal and reuse. The extent of leachable As and Se depends on several factors related to environmental conditions and fly ash characteristics. Previous studies employed various methods to delineate the concentration, chemical form, and distribution of As and Se in fly ash materials. However, few studies have attempted to directly correlate these properties to mobilization parameters relevant to disposal and reuse. Instead, the coal residuals industries often rely upon standardized leaching protocols that can be laborious or involve hazardous chemicals. The goals of the project were to: 1) Develop and evaluate a characterization protocol that can be used to screen fly ash samples for leachability of As and Se; 2) Characterize As, Se, and associated constituents of fly ash particles at multiple length scales (nanometer to micrometer) to determine if elemental associations differ as a function of the resolution of characterization; and 3) Establish a predictive model for the chemical composition of coal ash produced annually at major U.S. coal fired power facilities on 50-year national coal supply records. For the first objective, we performed leaching experiments with 52 fly ash samples collected from 15 different U.S. power plants and representing coal feedstocks from the three major domestic coal regions. For this work, we assessed the mobilization potential of As and Se in fly ash based on standardized leaching protocols and performed multivariate and lasso regression analyses to explore correlations of leachable As and Se contents with characteristics such as major element contents, loss on ignition (LOI) and pH. The results of regression models indicated that major elements (Fe, Ca, Al) for a wide range of fly ashes can serve as predictor variables for the leaching potential of As, but not for Se. LOI and pH were not important predictive variables in the models. Both regression approaches resulted in relatively strong fits for leachable As (correlation coefficient R 2 = 0.78 for both models) compared to models for leachable Se (R 2 = 0.49). Overall, these results suggest that correlation models combined with on-site elemental analysis with portable analyzers may enable a screening method for leachable As in coal ash. For the second objective, we utilized nanoscale 2-D imaging (30-50 nm spot size) with the Hard X-ray Nanoprobe (HXN) in combination with microprobe X-ray capabilities (~5 µm resolution) to determine As and Se elemental associations in fly ash particles. Speciation of As and Se was also measured at the nano- to microscale with X-ray absorption spectroscopy. The enhanced resolution of HXN showed As and Se that were diffusely located around or comingled with Ca- and Fe-rich particles. The results also showed nanoparticles of Se attached to the surface of fly ash grains. Overall, a comparison of As and Se species across scales highlights the heterogeneity and complexity of chemical associations for these trace elements of concern in coal fly ash. For the final objective, we developed a predictive model for major element composition of coal ash in reserve at disposal sites of major U.S. coal fired power plants. This model was constructed from coal purchase records of 705 power stations from 1973-2022 and was trained on coal ash composition data showing that coal ash elemental composition is strongly associated with the source of feedstock coal. The model showed regional shifts in the major element contents of ash produced by power plants in the last 50 years, particularly for calcium and iron (expressed as %CaO and %Fe 2 O 3 ), as coal-fired power stations changed their source of coal over this time frame. Our approach enables an estimation of coal ash chemical composition that is stored in waste impoundments at individual power stations. Such information can help delineate the regional market potential for material applications that would utilize coal ash harvested from disposal sites across the U.S.

01 COAL, LIGNITE, AND PEAT↗