SEARCH · Search NASA

Results for “LIFT DISTRIBUTION - WINGS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Characterization of Surface Heterogeneity-Induced Convection Using Cluster Analysis

Land-atmosphere interactions and boundary layer processes often control the formation of shallow clouds and subsequently deep convective precipitation over the Southern Great Plains. In this study, we examine the impacts of large-scale advection on the cloud populations and land-atmospheric coupling observed during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign in 2016. In this work, we performed two Large Eddy Simulations (LES) using the Weather Research and Forecasting Model for a day in which the transition from clear-sky to shallow precipitating clouds and isolated deep convection was observed. The control simulation reproduced the overall distribution of cloud populations by using realistic soil conditions with an interactive land model. In the sensitivity simulation where large-scale advection is removed, a strong relationship between the land and boundary layer is found. To study the timing, location, and intensity of convective initiation and the relationship of clouds with land surface properties, a cluster analysis of equivalent potential temperature is performed for the simulation without large-scale advection. That analysis shows that convective clouds first form over regions with higher surface sensible heat flux. Precipitation from those convective clouds likely triggers new updrafts nearby about two hours later through the lifting associated with cold pools. Furthermore, the cluster analysis also shows that in addition to the spatial pattern of soil moisture, land use and soil texture in western Oklahoma also influence the location of convective initiation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Correlating dislocation structures to basal and prismatic slip bands near grain boundaries in tensile-tested Mg–4Al using a multiscale electron microscopy approach

Dislocation structures where basal and prismatic slip bands meet grain boundaries in tensile-tested Mg with 4 wt% Al (Mg–4Al) were studied using a multiscale electron microscopy approach to explore the assumptions made for the dislocation pileup theory: a) slip bands consist of dislocation arrays piling up near grain boundaries, and b) stress concentration in the adjacent grain is solely caused by dislocation pileup. After post-testing electron backscatter diffraction (EBSD) scans, focused-ion-beam (FIB) lift-out specimens were prepared from regions of interest so that the specimen plane is parallel to the bulk sample surface to allow for direct correlation between mesoscale and microscale characterization. High dislocation density and the dislocation arrays in the slip bands corroborated with the first assumption. Evidence of plastic deformation in all grains, however, showed that the second assumption was false. Additionally, dislocation structures including crystal rotation boundaries and dislocations dissociating to the basal plane provide additional stress relief mechanisms near grain boundaries that are not accounted for in many crystal plasticity models, which requires auxiliary models to more accurately predict the stress and strain distribution in a microstructure.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Role of Surface Hydroxyls and Lattice Oxygen in Governing Selectivity and Restructuring During Acetic Acid Conversion on Fe 3 O 4 (001)

Understanding the reactivity of carboxylic acids on metal oxide surfaces is critical for elucidating ketonization mechanisms relevant to biomass upgrading. Here, we investigate the adsorption and thermal decomposition of acetic acid (CH 3 COOH) on Fe 3 O 4 (001) using scanning tunneling microscopy (STM), temperature-programmed reaction spectroscopy (TPRS), and X-ray photoelectron spectroscopy (XPS). At room temperature, acetic acid adsorbs dissociatively to form ordered bidentate acetate (CH 3 COO) overlayer that lifts the (√2 × √2)R 45° surface reconstruction. TPRS reveals ketene (CH 2 CO) as the dominant product, along with CO, CO 2 , and H 2 O, the latter evolving via a Mars–van Krevelen (MvK) mechanism. Isotopic labeling shows preferential CO 2 formation from the carboxyl carbon and a more balanced CO/CO 2 ratio from the methyl carbon, suggesting distinct oxidation pathways. STM imaging reveals embedded acetate intermediates filling surface oxygen vacancies created in MvK steps. Upon product formation completion (~700 K), extensive surface etching is observed, with pits elongated along the octahedral Fe rows. Approximately 20% of the surface oxygen is removed, consistent with vacancy formation stoichiometry inferred from product distributions. These findings demonstrate that carboxylate-induced restructuring of Fe 3 O 4 (001) involves both surface healing and reduction processes, offering mechanistic insights relevant to ketonization and broader carboxylic acid chemistry on metal oxides.

acetic acid↗

DOE OSTI2025

Vertical column dual-comb spectroscopy to a TBS

Open-path dual-frequency-comb spectroscopy (DCS) is a broadband, high spectral resolution, and high precision method for measuring gas concentrations over kilometer-scale paths. It has been used for detection and quantification of emissions of pollutants, hazardous gasses, and greenhouse gasses (GHGs). DCS has been shown to measure trace-gas mixing ratios with 0.14%-0.4% agreement between instruments. To achieve high signal-to-noise ratios (SNRs) with DCS, comb light is targeted onto a retroreflector at the end of the measurement path, which returns the signal to a detector co-located with the launch signal. Installing the retroreflector on a mobile platform such as a balloon or unmanned aerial vehicle (UAV) extends the capabilities of DCS by enabling variable path lengths, greater mobility, and access to higher altitudes. Mobile-target DCS has many uses in plume and leak detection, emissions modeling, and planetary boundary layer (PBL) studies. It is a promising method for observing vertical distributions of GHGs and mixing processes in the PBL, which are difficult to measure but important for pollution and climate monitoring as well as for understanding transport of gasses through the atmosphere. Tethered balloons are an intriguing platform because they enable longer flight durations and higher altitudes than easily obtainable with a UAV and thus allow for column measurements up to and above the PBL. New measurements completed in October/November 2024 reach the highest altitudes above ground level yet achieved by mobile-target DCS. For these measurements, the retroreflector was mounted on a 7 m diameter tethered helium balloon. An actively tracking gimbal on the ground holds the DCS launch telescope and keeps the 5 cm beam pointed onto the retroreflector while the balloon is lifted, lowered, and moved by wind and turbulence.

atmosphere↗

DOE OSTI2024

Quantitative Identification of Dopant Occupation in Li‐Rich Cathodes

Elemental doping is widely used to improve the performance of cathode materials in lithium‐ion batteries. However, macroscopic/statistical investigation on how doping sites are distributed in the material lattice, despite being a key prerequisite for understanding and manipulating the doping effect, has not been effectively established. Herein, to solve this predicament, a universal strategy is proposed to quantitatively identify the locations of Al and Mg dopants in lithium‐rich layered oxides (LLOs). Solid evidence confirms that Al prefers to occupy the transition metal (TM) layer, while Mg evenly occupies both TM and Li layers. As a result, Mg significantly reduces the thickness of LiO 2 slabs at room temperature, which will increase the energy barrier of oxygen activation and enhance the structure stability of LLOs. The suppressed oxygen activity in Mg‐doped LLO can be kinetically unlocked at 55 °C. The different characteristics of Al and Mg enlighten an Al/Mg co‐doping strategy to optimize LLOs, which significantly improves the cycle performance while lifting the capacity. In conclusion, these insights from the quantitative identification of doping sites shed light on the manipulation of doping effects toward better cathodes.

25 ENERGY STORAGE↗

DOE OSTI2022

Upper troposphere dust belt formation processes vary seasonally and spatially in the Northern Hemisphere

Dust aerosols impact global energy balance substantially by acting as efficient ice nuclei to alter cold cloud properties. However, the estimate of dust indirect effect remains uncertain due to simulating dust distributions poorly and lacking reliable dust observations, especially in the upper-troposphere. Here, we characterize and understand upper-troposphere dust sources and transport with an improved dust dataset derived from A-train satellite lidar and radar measurements and an air parcel trajectory model. The distinct upper-troposphere dust belt over the northern hemisphere has seasonally varying base and top heights of 3.65 ± 2.84 and 8.35 ± 1.50 km above mean sea level and its column loading is strongest during spring (March-April-May). The out-of-phase annual cycles of mid-level dust concentration and westerly wind over source regions control the seasonal upper-tropospheric dust loading variations. African deserts contribute the most (46.3%) to the upper-troposphere dust belt in spring and the synoptic trough is the leading (49%) dust lifting mechanism.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

GIS-Based Modeling of Contaminated Soil Volumes at Multiple Sites in the Formerly Utilized Sites Remedial Action Program - 20149

The remediation of hazardous, toxic, and radioactive waste (HTRW) sites produces cost-related risks associated with the estimation of contaminated soil or debris volumes. Historical risk-management techniques include cost contingencies to cover volume uncertainties that affect project budgeting and decision-making. The Buffalo District teamed with project partners to lessen volume uncertainty and reduce project risks at multiple HTRW sites managed under the Formerly Utilized Sites Remedial Action Program (FUSRAP). Historical remedial investigations under FUSRAP commonly identified the presence of radiological material in site media, the associated human health risk, and then areas of remediation. To manage remedial execution and reduce risk, pre-design or remediation-phase sampling essentially 'chased' contamination, which was not conducive to efficient predictive budgeting derived from Feasibility Study (FS) cost analyses. The Buffalo District first optimized their approach to better understand volume uncertainty by utilizing the Argonne National Laboratory's Bayesian Approaches for Adaptive Spatial Sampling (BAASS) software [1]. BAASS processed soft data (e.g., gamma walk-over data) and spatial sampling data to estimate the lateral extent of contaminated soil irrespective of depth (i.e., gross contamination extent) and define areas of contaminant uncertainty. The software performed a binary transformation of contaminant concentrations at all sampling points based upon remedial action goals or a sum of ratios approach (i.e., clean, impacted, or range of impacts in soil). The model produced two-dimensional (horizontal) contaminant probability contours and statistical uncertainty in the sampling coverage and resulting contaminant extents. This method was translated vertically by partitioning the sampling data into depth brackets that produced a stacked representation of contaminant extents and uncertainty in the subsurface (i.e., similar to construction lifts). The results commonly led to a better understanding of project uncertainty and the need for sampling strategies that produce high-confidence soil volumes, which control costs. The BAASS-based delineations were eventually replaced by Empirical Bayesian Kriging (EBK) methods available in ArcGIS Spatial or 3D Analysts [2]. The EBK method calculates contaminant probability zones derived from user-controlled semivariograms of the spatial datasets. The resulting probability zones (e.g., 50% or 80% of contaminant probability) represent the two-dimensional surface delineation of the overall horizontal remedial area, similarly to BAASS. However, unlike BAASS, the vertical sampling data within these probability zones became vertical control points to contour a subterranean surface that connects subsurface points to the land-surface delineations of contamination. The resulting representation of horizontal and vertical impacts within an enclosed envelop (volume) of soil included uncertainty distributions that are used to plan uncertainty-reduction sampling. These data-driven and math-based models of three-dimensional sampling results produced well-bounded remedial volumes for project planning and better uncertainty predictions during project budgeting. The EBK method was applied to several FUSRAP sites managed by the Buffalo District and compared to less rigorously modeled sites previously remediated by the District. The comparison of modeled to actual remediated volumes provide a basis for validating the volume-estimation method. This comparison is important to ensure modeled volumes match physical boundaries of site remediation. FUSRAP sites with denser investigative sampling and lesser volume uncertainty proved useful in remedial planning and contracting. The Buffalo District noted that historical sites with sparser sampling arrays had greater disparity between estimated volumes and final remedial volumes. The benefit achieved over the cost of detailed soil sampling appears positive for FUSRAP projects, especially where impacts vary widely and appear unbounded by investigation-phase sampling. The subsequent Empirical Bayesian Kriging of contamination coupled with vertical contouring for soil estimations reduces uncertainty in soil volumes or indicates where sampling is required to reduce uncertainty, which together optimize remedial planning and budgeting. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Drop size distribution broadening mechanisms in a bin microphysics Eulerian model

In this study, processes that broaden drop size distributions (DSDs) in Eulerian models with two-moment bin microphysics are analyzed. Numerous tests are performed to isolate the effects of different physical mechanisms that broaden DSDs in two- and three-dimensional Weather Research and Forecasting model simulations of an idealized ice-free cumulus cloud. Additionally, sensitivity of these effects to modifying horizontal and vertical model grid spacings is also examined. As expected, collision-coalescence is a key process broadening the modeled DSDs. In-cloud droplet activation also contributes substantially to DSD broadening, whereas evaporation has only a minor effect and sedimentation has little effect. Cloud dilution (mixing of cloud-free and cloudy air) also broadens the DSDs considerably, whether or not it is accompanied by evaporation. This mechanism involves the reduction of droplet concentration from dilution along the cloud’s lateral edges, leading to locally high supersaturation and enhanced drop growth when this air is subsequently lifted in the updraft. DSD broadening ensues when the DSDs are mixed with those from the cloud core. Decreasing the horizontal and vertical model grid spacings from 100 m to 30 m has limited impact on the DSDs. However, when these physical broadening mechanisms (in-cloud activation, collision-coalescence, dilution, etc.) are turned off, there is a reduction of DSD width by up to ~20-50% when the vertical grid spacing is decreased from 100 m to 30 m, consistent with effects of artificial broadening from vertical numerical diffusion. Nonetheless, this artificial numerical broadening appears to be relatively unimportant overall for DSD broadening when physically based broadening mechanisms in the model are included for this cumulus case.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A phenomenological rate of injection model for predicting fuel injection with application to mixture formation in light-duty diesel engines

Fuel injection rate laws are one of the most important pieces of information needed when modeling engine combustion with computational fluid dynamics. In this study, a simple phenomenological model of a common-rail injector was developed and calibrated for the Bosch CRI2.2 platform. The model requires three tunable parameter fits, making it relatively easy to calibrate and suitable for injector modeling when high-fidelity information about the internal injector’s geometry and electrical circuit details are not available. Each injection pulse is modeled as a sequence of up to four stages: an injection needle mechanical opening transient, a full-lift viscous flow inertial transient, a Bernoulli steady-state stage, and a needle descent transient. Parameters for each stage are obtained as polynomial fits from measured injection rate properties. The model enforces total injected mass, and the intermediate stages are only introduced if the injection pulse duration is long enough. Experimental rates of injection from two separate campaigns on the same injector were used to calibrate the model. The model was first validated against measured injection rate laws featuring pilot injections, short partially premixed combustion pulses, and conventional diesel combustion injection strategies. Then, it was employed as an input to engine computational fluid dynamics simulations, which were run to simulate experiments of mixture formation in an optically accessible light-duty diesel engine. It was found that, though simple, this model is capable of predicting both pilot and main injection pulse mass flow rates well: the simulations yielded accurate predictions of in-cylinder equivalence ratio distributions from injection strategies for both partially premixed combustion and pilot injections. Also, once calibrated, the model produced appropriate results for a wide range of injected mass and rail pressure values. Finally, it was observed that usage of such a relatively simple model can be a good choice when high-fidelity injection rate input and highly detailed information of the injector’s geometry and operation are not available, particularly as noticeable discrepancies can be present also among different experimental campaigns on similar hardware.

42 ENGINEERING↗

DOE OSTI2023

Measurement of needle and armature dynamics in a gasoline direct injector by high-speed neutron imaging

In modern spark ignition engines, precise delivery of fuel with gasoline direct injection has become increasingly important in the effort to meet ever stricter efficiency and emissions regulations. Use of multiple small close-coupled injections has become more common in attempt to precisely control fuel distribution in the cylinder, but these strategies are hindered by nonlinear injection effects due to operation in the ballistic region of the solenoid-operated valve and due to armature bounce at the end of injection. Understanding the internal dynamics of the injector is crucial to minimizing and controlling non-linearity and shot-to-shot variation, but the experimental techniques available to date are capable only of tracking the position of either the top of the needle (via laser sensors or by monitoring current and voltage in the solenoid coil) or the bottom of the needle (via transparent nozzles or high-speed x-ray imaging). Furthermore, a complete picture of the axial and radial motion of valve needle has until now remained elusive. In this work, we present high-speed ensemble neutron transmission imaging of an entire 8-hole gasoline direct injector operating at 200 bar, allowing for both visualization and quantification of the actuation dynamics including lift and wobble of the valve ball, oscillation and bending of the valve needle, lift, rocking, and bounce of the armature, compression of the springs, and radial swelling of the solenoid during energization. Because neutrons offer high penetration through the metal injector while also being sensitive to the 1 H nuclei in fuel molecules, it is also possible to simultaneously see the fluid dynamics of the injection process, including filling of the sac volume, emanation of the spray through the nozzle holes and into the downstream gas, and the formation and evolution of fuel films on the tip of the injector and the walls of the spray container.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Equation-Free Coarse Control of Distributed Parameter Systems via Local Neural Operators

The control of high-dimensional distributed parameter systems (DPS) remains a challenge when explicit coarse-grained equations are unavailable. Classical equation-free (EF) approaches rely on fine-scale simulators treated as black-box timesteppers. However, repeated simulations for steady-state computation, linearization, and control design are often computationally prohibitive, or the microscopic timestepper may not even be available, leaving us with data as the only resource. We propose a data-driven alternative that uses local neural operators, trained on spatiotemporal microscopic/mesoscopic data, to obtain efficient short-time solution operators. These surrogates are employed within Krylov subspace methods to compute coarse steady and unsteady-states, while also providing Jacobian information in a matrix-free manner. Krylov-Arnoldi iterations then approximate the dominant eigenspectrum, yielding reduced models that capture the open-loop slow dynamics without explicit Jacobian assembly. Both discrete-time Linear Quadratic Regulator (dLQR) and pole-placement (PP) controllers are based on this reduced system and lifted back to the full nonlinear dynamics, thereby closing the feedback loop.

93B52, 93C20, 47N70, 65J15, 65M32, 68T07, 68T20, 6↗

DOE OSTI2024

Atom probe datasets from neutron irradiated Fe-Cr alloys

A series of model Fe–Cr alloys containing 3–18 at.% Cr was neutron irradiated at a nominal temperature of 563 K to 1.82 dpa. Solute distributions were analyzed by atom probe tomography, which revealed α′ precipitation for alloys containing more than 9 at.% Cr. Both the Cr concentration dependence of α′ precipitation and the measured matrix compositions are in agreement with the recently published Fe–Cr phase diagrams. An irradiation-accelerated precipitation process is strongly suggested. Irradiation was carried out in the Advance Test Reactor (ATR) at Idaho National Laboratory. A series of six Fe–Cr alloys of nominal compositions 3, 6, 9, 12, 15 and 18 at.% Cr was irradiated at a neutron fluence (E > 1 MeV) of 1.1 × 1021 n cm−2 at 563 ± 15 K and to a damage level of 1.82 displacements per atom (dpa). Nominal neutron flux and dpa rate are 2.3 × 1014 n cm−2 s−1 and 3.4 × 10−7 dpa s−1, respectively. The microstructures of the Fe–Cr alloys were studied by atom probe tomography (APT) using a Cameca 4000X HR instrument. APT specimens were prepared by a standard lift-out process using a Quanta 3D 200i dual beam scanning electron microscope ensuring that the analyses were performed away from grain boundaries.

Bachhav, Mukesh↗

DOE OSTI2025

6cr

Bachhav, Mukesh↗

DOE OSTI2020

Impact of the adhesion-force lever-arm “a” on the Rock ‘n’ Roll resuspension model and how to compute it from contact mechanics

The paper proposes a methodology to fill a gap in the context of mechanistic particle resuspension models related to the adhesion and lift force lever arms ( a ): to date, only one experimental observation exists from a centrifuge experiment with a specific particle/surface pair which found r P / a = 100 and assumed that the contact points were symmetric with respect to the particle center. However, the Reeks and Hall’s (2001) model is highly sensitive to this parameter, which therefore cannot be used as-is for different particle diameters or materials. This paper proposes a methodology to compute the adhesion and lift force lever arms ( a A and a L , respectively) based on Atomic Force Microcopy (AFM) surface roughness measurements and contact mechanics arguments. In practice, a spherical, perfectly smooth particle is brought in contact (computationally) with the measured surface roughness and is then rotated by the point of contact to determine the two or more asperities over which the particle would sit on the surface. A method to compute the asperity radius of curvature is devised to compute the asperity deformation and the correct partitioning of the total deformation between the particle and the asperity is presented. For the glass surface under consideration, the asperity deformation was negligible, i.e., it was below the van der Waals radius, which is the limit of applicability of contact mechanics. The Reeks and Hall’s (2001) model is then extended to account for the joint probability of a A and a L . Significant differences in particle resuspended after a day in common atmospheric conditions are found between the approach that uses the lever-arm distributions or their mean values. Also, important differences are found when comparing the lever-arm distribution and r P / a = 100. Hence, one should be cautious of using r P / a = 100 for surface/particle pairs differing from those in the original study. Finally, limits of applicability of the theories used are presented together with considerations regarding how to extend the method to non-spherical particles and different materials.

42 ENGINEERING↗

DOE OSTI2024

Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets

This repository contains simulated outputs from superwell – a hydro-economic tool for long-term assessment of groundwater cost and supply – providing globally gridded groundwater extractable volumes and associated unit costs ($/km³) for accessible groundwater production, based on a variety of user-defined depletion and ponded depth scenarios. Key model documentation: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 Find the source code of the superwell model on GitHub: https://github.com/JGCRI/superwell Repository Overview Main output: superwell_outputs.7z contains 6 files (4.5 GB) named as superwell_py_deep_all_0.*PD_0.*DL.csv. These files present superwell outputs of global groundwater extraction volumes and cost estimates on a 0.5° scale for six scenarios with different Ponded Depth (PD; 0.3 and 0.6 m) and Depletion Limit (DL; 5%, 25%, and 40% of available volume) targets over the entire pumping lifetime of a grid cell superwell_py_deep_all_0.3PD_0.25DL_sample_100.csv contains superwell outputs for 100 data points sampled to match the global inputs' distribution superwell_py_deep_all_0.3PD_0.25DL_Grid_72548.csv contains superwell output for a single grid cell concept_v5.png provides an overview of the superwell workflow Outputs Description year_number: year of pumping depletion_limit: set depletion limit (DL) as a volume fraction of total available groundwater Mappings: continent, country, gcam_basin_id, Basin_long_name, grid_id: geographic identifiers and basin information Inputs: grid_area (km²): area of the grid cell whyclass: hydrogeological classification of the aquifer permeability (m/day), porosity (%), total_thickness (m), depth_to_water (m): aquifer properties. The geo-processed input data has been published separately: https://doi.org/10.57931/2307831 Model outputs: orig_aqfr_sat_thickness (m), aqfr_sat_thickness (m): original and remaining/instantaneous saturated thickness of the aquifer hydraulic_conductivity (m/day), transmissivity (m²/day): hydraulic properties of the aquifer radius_of_influence (m), areal_extent (km²): well radius and area of influence from the center of the well number_of_wells (-): number of wells in a grid cell determined by a ratio of well area and grid area max_drawdown (m), drawdown (m), drawdown_interference (m): well and aquifer drawdown during extraction total_head (m): total lift for the groundwater (depth to water plus drawdown) total_well_length (m): total depth of wells drilled well_yield (m³/day): pumping rate or well yield power (kW), energy (kWh): power and energy required for pumping groundwater Volume Outputs: volume_produced_perwell (m³), cumulative_vol_produced_perwell (m³): production volume metrics per well volume_produced_allwells (m³), cumulative_vol_produced_allwells (m³): aggregate extraction volumes for all wells in a grid cell available_volume (m³): available groundwater in storage for the grid cell as determined by aquifer properties depleted_vol_fraction: fraction of total volume pumped over available volumes in a grid cell (same as depletion limit) Cost Outputs: well_installation_cost ($): well installation cost based on the hydrogeological complexity of the aquifer annual_capital_cost, maintenance_cost, nonenergy_cost ($): nonenergy costs energy_cost_rate ($/kWh): electricity rate energy_cost ($): energy cost of pumping groundwater total_cost_perwell ($), total_cost_allwells ($): total annual energy and non-energy cost for each and all wells in a grid cell a unit_cost ($/m³), unit_cost_per_km3 ($/km³), unit_cost_per_acreft ($/acre-ft): total cost of pumping a unit of groundwater, indicated for different spatial units Key Resources Model documentation: Niazi, H., Ferencz, S., Graham, N., Yoon, J., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024; In-prep). Long-term Hydro-economic Assessment Tool for Evaluating Global Groundwater Cost and Supply: Superwell v1. Geoscientific Model Development. Input data: Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MSD-LIVE Data repository. https://doi.org/10.57931/2307831 superwell source code: https://github.com/JGCRI/superwell Cite as Niazi, H., Ferencz, S., Yoon, J., Graham, N., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024). Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets. MSD-LIVE Data repository. https://doi.org/10.57931/2307832 Contact Reach out to Hassan Niazi or Stephen Ferencz or open an issue in the superwell repository for questions or suggestions.

Earth Systems↗

DOE OSTI2024

LES and RANS Spray Combustion Analysis of OME 3-5 and n-Dodecane

Clean-burning oxygenated and synthetic fuels derived from renewable power, so-called e-fuels, are a promising pathway to decarbonize compression–ignition engines. Polyoxymethylene dimethyl ethers (PODEs or OMEs) are one candidate of such fuels with good prospects. Their lack of carbon-to-carbon bonds and high concentration of chemically bound oxygen effectively negate the emergence of polycyclic aromatic hydrocarbons (PAHs) and even their precursors like acetylene (C 2 H 2 ), enabling soot-free combustion without the soot-NO x trade-off common for diesel engines. The differences in the spray combustion process for OMEs and diesel-like reference fuels like n-dodecane and their potential implications on engine applications include discrepancies in the observed ignition delay, the stabilized flame lift-off location, and significant deviations in high-temperature flame morphology. For CFD simulations, the accurate modeling and prediction of these differences between OMEs and n-dodecane proved challenging. This study investigates the spray combustion process of an OME 3 − 5 mixture and n-dodecane with advanced optical diagnostics, Reynolds-Averaged Navier–Stokes (RANS), and Large-Eddy Simulations (LESs) within a constant-volume vessel. Cool-flame and high-temperature combustion were measured simultaneously via high-speed (50 kHz) imaging with formaldehyde (CH 2 O) planar laser-induced fluorescence (PLIF) representing the former and line-of-sight OH* chemiluminescence the latter. Both RANS and LES simulations accurately describe the cool-flame development process with the formation of CH 2 O. However, CH 2 O consumption and the onset of high-temperature reactions, signaled by the rise of OH* levels, show significant deviations between RANS, LES, and experiments as well as between n-dodecane and OME. A focus is set on the quality of the simulated results compared to the experimentally observed spatial distribution of OH*, especially in OME fuel-rich regions. The influence of the turbulence modeling is investigated for the two distinct ambient temperatures of 900 K and 1200 K within the Engine Combustion Network Spray A setup. The capabilities and limitations of the RANS simulations are demonstrated with the initial cool-flame propagation and periodic oscillations of CH 2 O formation/consumption during the quasi-steady combustion period captured by the LES.

CFD↗

DOE OSTI2023

TEAMER: Water Tunnel Data from Testing the Pterofin Skimmer Concept

Pterofin's Skimmer concept relies on a flapping and pitching hydrofoil to extract hydrokinetic energy from water flows. The concept aims to utilize unsteady fluid dynamics phenomena (added mass, shed vorticity, and unsteady boundary layer development) to achieve higher lift coefficients, enabling increased power density of the hydrokinetic device and a fundamental shift in the rpm/torque scaling of the power take off compared with turbines. The Applied Research Laboratory at Penn State, in collaboration with Pterofin, designed and built a proof-of-concept flapping/pitching mechanism which was subsequently tested in ARL's 12-inch water tunnel facility. The mechanical power supplied to or extracted from the mechanism was measured for a range of hydrofoils provided by Pterofin over operating conditions including reduced frequency, Reynolds number, and the ratio between pitching and flapping amplitudes. The power lost to friction in the mechanism was removed from the net power measurement by means of a bare hub tare, with the resultant hydrodynamic power being used to calculate a mechanism-independent and non-dimensional power coefficient. The product of this effort is a dataset describing the power coefficient of a hydrofoil having simultaneous pitching and flapping motions, both of which are approximately sinusoidal. Power coefficients were collected for a range of primary design variables including: - Reduced frequency: 0.01 to 0.95 - Pitching/flapping peak angle ratio: 1.5 to 3.0 - Chord-based Reynolds number: 60,000 to 560,000 Secondary design variables relating to the hydrofoil geometry were explored including: - Aspect ratio - Planform shape - Section thickness distribution - Hydrofoil position relative to the pitching axis - Hydrofoil sweep angle relative to the pitching axis Measured data are provided in mean and time series formats. MATLAB scripts are provided which can be used to generate figures of time-averaged and phase-averaged hydrodynamic power coefficients calculated from the measured data. A complete description of the experiment and data reduction can be found in the Post Access Report for the Pterofin Skimmer test effort which will be available on the TEAMER website. This work was supported by the Pacific Energy Ocean Trust via a TEAMER award.

16 TIDAL AND WAVE POWER↗

DOE OSTI2025

Compositing and Characterization of SE Quadrant Waste Exemplars

This report outlines the experimental investigation and characterization of transport properties in Hanford SE quadrant High-Level Waste (HLW). The goal of the study was to establish baseline behaviors of bulk composite rheology and settling characteristics to facilitate waste treatment process design for the Waste Treatment and Immobilization Plant (WTP) and avoid waste conditions and properties favorable to bubble cascade gas release events. The study focused on two major objectives: 1) identifying, obtaining, and preparing relevant Hanford tank waste samples for evaluation and 2) quantifying the “as-received” rheology and transport properties of the samples. Twenty-three centrifuged core segments originating from tanks AN-101, AN-106, and AW-105 were selected based on compositional relevance to SE quadrant PUREX cladding waste. These materials were composited into five waste composites enriched with target analytes: aluminum (Al), iron (Fe), phosphate (PO 4 ), uranium (U), and zirconium (Zr). Physical property and transport testing examined particle size distributions, bulk densities, settling behaviors, rheological properties, shear strengths, and just-suspended mixing speeds (NJS). Testing revealed two distinct composite classifications based on rheological characteristics: non-Newtonian composites (Fe and PO 4 ) and Newtonian composites (Al, U, and Zr). The Fe and PO 4 composites exhibited slow settling rates and reduced mobilization proclivity, attributable to strong particle-particle interactions and the formation of yield structures within non-Newtonian slurries. In contrast, the Al, U, and Zr composites displayed rapid settling and dense compaction behaviors, indicative of minimal structuring and interactions. Shear strengths for all composites were generally low relative to prior studies of SE quadrant waste, with only the U composite showing elevated strength approaching values reported in previous literature. Repeat shear strength measurements revealed contributions from dense granular material in the U composite and stronger cohesive properties in the Al composite. Settling data highlighted hindered settling behavior, with rates falling more than one order of magnitude below estimates based on Stokes’ law and rate decreasing as composite UDS content increased. NJS testing demonstrated different mobilization behaviors between cohesive and granular composites. The Fe composite required the highest mixing rate for resuspension, while the Al composite was the easiest to resuspend. Comparison of measured NJS against predictions made using the Zwietering correlation suggests non-Newtonian behavior alters resuspension mechanics, rendering non-Newtonian systems more stable against resuspension lift forces relative to their granular counterparts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗