Search NASA⌕ Search

SEARCH · Search NASA

Results for “Homogenization”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Rapid quantification of whole seed fatty acid amount, composition, and shape phenotypes from diverse oilseed species with large differences in seed size

Seed oils are widely used in the food, biofuel, and industrial feedstock industries, with their utility and value determined by total oil content and fatty acid composition. Current high throughput seed oil analysis methods either lack accuracy in total fatty acid profiling or require extensive labor for lipid extraction prior to derivatization to fatty acid methyl esters (FAME) and quantification by gas chromatography (GC). Alternatively, direct whole seed FAME production methods have been developed for the very small seeds in the model species Arabidopsis thaliana but these have generally not been adapted to larger seeds of most oilseed crops. High-throughput direct whole seed FAME production methods were optimized for seeds up to 5 mg each utilizing acid-catalyzed esterification. For the oilseed species Camelina sativa, Thlaspi avernse (pennycress), Cuphea viscosissima, and Brassica napus (var. Canola), the total seed fatty acid content and composition from direct seed esterification to FAME matched that of lipid extract derivatization demonstrating the accuracy of the methods. In combination with seed phenotyping using GridFree, this approach enabled the development of a rapid pipeline for simultaneous seed weight, count, size/shape phenotyping, and oil analysis. For the larger and tougher seeds produced by Limnanthes alba (Meadowfoam) and Cannabis sativa L. (hemp) the whole seed acid-based method proved insufficient, and prior laborious homogenization of seeds was required. Therefore, a rapid one-tube bead homogenization and base catalyzed-esterification method was developed. Base-derived fatty acid esterification cannot derivatize free fatty acids leading to slightly lower total seed fatty acid than acid-catalyzed methods, however the seed oil content and fatty acid composition that is valuable for screening large numbers of samples in research populations was accurately measured. New rapid whole seed fatty acid esterification and phenotyping protocols were developed to accurately assess oilseed lipid content. These methods are particularly valuable in oilseed research, breeding, and engineering applications where efficient analysis of large numbers of samples and accurate oil fatty acid profiling is essential. While having been developed for current and emerging oilseed crops, these methods also provide a foundation from which protocols might be established for new and emerging crop species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Data for Impact of Vertical and Seasonal Variation in Leaf Traits on Simulating Soybean Canopy Photosynthesis via 1D and 3D Modeling

Accurate modeling of photosynthesis is crucial for predicting crop productivity and quantifying the carbon cycle in agroecosystems. Leaf traits are essential inputs for modeling canopy photosynthesis. Yet, many existing models still use fixed plant functional type (PTF)-based values to parameterize leaf traits under a big-leaf or two-big-leaf assumption, neglecting their vertical profiles and seasonal changes. This simplification may introduce significant uncertainties in estimating gross primary productivity (GPP). In this study, we simulated soybean GPP and tested the effects of vertical and seasonal variation in three key leaf photosynthetic traits: the maximum carboxylation rate at 25 °C (Vcmax25), leaf chlorophyll content (LCC), and leaf mass per area (LMA) in the 1D-SCOPE and 3D-Helios models. Weekly field measurements were conducted during the growing season of 2024 to support the simulation. We designed ten leaf trait parameterization schemes by incorporating different combinations of vertical profiles and seasonal changes, while assuming homogeneous canopy architecture in both models. Our results revealed that Vcmax25 vertical and seasonal variation had the strongest influence on simulated GPP in both 1D and 3D models, while LCC and LMA effects were minimal. Particularly, the scheme with an empirically parameterized Vcmax25 profile achieved comparable performance to the scheme with the measured Vcmax25 profile. Both 1D-SCOPE and 3D-Helios accurately modeled GPP (SCOPE: R2 = 0.87, Bias = 0.55 µmol m⁻² s⁻¹; Helios: R2 = 0.9, Bias = 0.22 µmol m⁻² s⁻¹) under the most complex scheme, and their responses to vertical and seasonal variation in leaf traits were consistent, demonstrating the robustness of our findings. Based on our findings, we propose a scalable framework for parameterizing leaf traits to improve GPP simulations. This study contributes to improving the representation of leaf trait dynamics in canopy-level photosynthesis models, potentially enhancing our ability to predict crop productivity and understand agroecosystem carbon dynamics.

Photosynthesis↗
DOE OSTI2025

Ink-sprayed transparent ceramic multilayer planar waveguides for ribbon lasers

Multilayer transparent ceramic planar waveguides were fabricated using a newly developed particle-loaded ink spraying method, characterized, and lased. Waveguides with five Yb:YAG gain regions and four Lu:YAG index-matched regions were produced for application as a “ribbon laser,” as well as a homogeneous Yb:YAG waveguide. Elemental mapping revealed that the waveguide thicknesses varied from 40 to 69 μm, comprised of layers of less than 10 µm. The Yb concentration profiles were fit to a diffusion profile. The homogeneous and multilayer waveguides were tested in a laser cavity and produced similar slope efficiencies, with a maximum of 31% achieved. In conclusion, their beam profiles provided evidence of the ribbon structure affecting the preferred mode.

Osborne, Ross [Lawrence Livermore National Laborat↗
DOE OSTI2024

Solidification of SRPPF Aqueous Recovery Liquid: Process Disruptions

Savannah River National Laboratory has identified a magnesium oxysulfate grout formulation to solidify the liquid effluent from the Savannah River Plutonium Processing Facility (SRPPF) Aqueous Recovery System (ARS). The formulation uses reactive, light burnt MgO, anhydrous MgSO 4 , and dead burnt MgO and provides good mixability, similar density to the original Portland Cement-based mix, and a leachate pH of 9.39, within the assumed Waste Isolation Pilot Plant (WIPP) brine pH range. The formulation will solidify the effluent using a solid first, lost paddle design, within a 55-gallon drum. A series of seven different process disruptions, or irregular processing during the expected ARS liquid effluent solidification process, were completed to assess the severity of the impact of the process disruption and how the formulation behaves under these conditions. Seven process disruptions were investigated including: No Mixing, Partial Mixing, Compacted/Stratified, Sealed, Overfilled, Underfilled, and Undermixed. The impact on solidification, leachate pH, Er (a Pu/Am surrogate) distribution, and density of the resulting grout were established. Overall, four of the process disruptions, No Mixing, Partial Mixing, Compacted/Stratified, and Overfilled, resulted in basic, unabsorbed liquid that would require additional solidification to be acceptable for WIPP. Most solids generated in this testing met the expected WIPP leachate pH requirements. Solids from the Compacted/Stratified process disruption test, however, produced a leachate pH slightly above the assumed range in the WIPP performance assessment. The distribution of Er for most of the tests was relatively consistent. The No Mixing process disruption test, however, resulted in non-homogeneous Er distribution, indicating that radioactive material would not be evenly distributed throughout this solid’s matrix. The Overfilled test is the only mix with a density slightly below the acceptable value. The Sealed and Overmixed tests produced a solid closest to the standard mix with no unabsorbed liquid, expected leachate pH levels, consistent Er concentrations, and expected density, indicating these process disruptions do not have a negative impact on the final grout form. Based on the other process disruptions, it is critical to have premixed dry materials and a well-defined mixing process to produce a solid, homogeneous magnesium oxysulfate cement, with no unabsorbed liquid. These process disruption tests show that the lack of adequate mixing is fundamental to creating an improperly solidified material that may require remediation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

VARI3D & PERSENT: Perturbation and Sensitivity Analysis (Revision 5)

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry-based transport code.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Managing Subsurface Pressure Buildup and Interference in Commercial-Scale CO 2 Storage Project with Proximal Injection Wells

Large-scale decarbonization using carbon capture and storage (CCS) is likely to involve many commercial-scale CO 2 storage projects located in close proximity to each other. This close proximity raises concerns over pressure interference among the storage projects. Pressure interference between injection and storage efforts can reduce the practicable CO 2 storage resource and force wells to inject CO 2 at a lower rate to avoid the fracture pressure thresholds per United States Environmental Protection Agency (EPA) Class VI well regulations to preserve injection and confining zone integrity and potentially mitigate against inducing seismic activity. These analyses employ numerical full-physics reservoir modeling to evaluate how pressure buildup fronts and CO 2 plumes evolve under commercial-scale injection volumes of CO 2 in which multiple storage sites located in close proximity occur in tandem. The simulation models mimic injection at pseudo basin-scale and assume homogeneous saline formation(s) as storage targets with a pair of upper and lower homogeneous seal layer/s. These analyses specifically investigate the efficacy of two basin-wide reservoir pressure management strategies in addressing the technical challenges associated with pressure buildup and CO 2 plume commingling. The strategies explored include: 1) enlarging the area of injection well spacing (WS) and 2) storing CO 2 in a stacked sequence (SSS) of saline formations compared to a single formation. The storage and confining zones properties assumed were held common across the scenarios, unless specified otherwise. Analyses results show that after injecting 4 million tons per year for 30 years using 4 separate wells (each injecting 1 million metric tons per year), the radius of CO 2 plume extends to a mere 3 km or less from injection wells. Meanwhile, the radius of pressure buildup ranges on the order of tens to a few hundreds of kilometers, depending on the magnitude of pressure buildup threshold that one would use to define the front. CO 2 plume commingling from different injection wells appears to occur 50 years post-injection, especially under scenarios with narrowly spaced (i.e., < 5 km apart) injection well locations. Findings from sensitivity cases on the well spacing suggest that storage formations modeled would require different well spacing to avoid fracture pressure thresholds. For instance, modeled storage formations with high fracture gradients (i.e., 0.8 psi/ft) would need less than 5–km well spacing, whereas those with lower fracture gradients (i.e., 0.7 psi/ft) would need approximately 20–km well spacing, based on assumed modeling parameters. Under stacked injection, the pressure challenges (described above) still exist but are more alleviated due to distributing the same injection volume across more available reservoir volume. These analyses demonstrate that stacked-sequence storage can effectively address the challenges, while still providing the same target CO 2 storage volumes and allowing a large number of storage projects to be deployed in the same basin by better utilizing the available storage resource across different reservoir depths. Among cases modeled, the resulting pressure buildup front is most suppressed when each storage project distributes injection volumes over several wells, each of which injects a portion of the total CO 2 across the stacked sequence. This strategy results in the smallest CO 2 aerial footprint amongst scenarios evaluated but also shows the largest reduction in the pressure buildup at the top of perforation at the injection wells (upwards of approximately 42 percent compared to the commercial-scale single-formation storage), the result of which is crucial to maintain caprock integrity. The findings presented by this research draw attention to the importance of greater coordination among storage operators and regulatory stakeholders to foster the upscaling and deployment of CCS. These analyses provide insights into required decision-making when considering multi-project deployment in a shared basin. Because these analyses evaluate a very specific geologic situation, they bear further investigations across other geologic situations.

42 ENGINEERING↗
DOE OSTI2025

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Discovery, Design, Synthesis and Testing of High Performance Structural Alloys (Final Technical Report)

The overarching goal of this project is to understand the phase stability and mechanical behavior of non-stoichiometric multi-principal element alloy (MPEA) materials. In order to identify suitable alloys, we plan to use a combinatorial thin film screening approach, in collaboration with scientists at Lawrence Berkeley National Laboratory who are performing computational work as well as complementary experimental work. Specific tasks within the scope of this project include the fabrication, using thin film deposition from six sputtering targets, of combinatorial samples with multi-dimensional gradients in composition and microstructure. These samples are studied to screen MPEA systems for promising candidate alloys with specific composition(s), based on characterization of composition, structure and mechanical behavior across the thin film. We want to produce single-phase MPEAs with chemical homogeneity in a given thin film region, simple grain structures, and no intermetallic phases present. Gradient films facilitate first-pass screening for desirable characteristics and inform the next stage of work that involves fabrication of bulk MPEA specimens for (tensile) mechanical testing and characterization. To make the bulk alloys, metal (elemental) pieces are melted to form MPEAs, followed by heat treatment to homogenize the composition and microstructure. A subset of alloys is also cast, using vacuum arc melting, to yield larger samples (diameter ~1 cm and length ~5-10 cm) and these allow us to assess viability of scale-up for the alloys in structural applications. Further processing plans include rolling and heat treatment to recrystallize selected bulk MPEAs and grow grains to different extents, in order to investigate size effects in the mechanical behavior of MPEAs. Microspecimen testing will be performed (primarily in tension) to assess the mechanical behavior over a range of temperatures. The deformation microstructure of mechanically tested alloys will be characterized using transmission electron microscopy (TEM) to provide a scientific basis for understanding the structure-property relationships in MPEA mechanical behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗
DOE OSTI2025

"Isotopics by nuclides" tool performance in InterSpec

This report provides a performance summary of the “Isotopics by Nuclides” tool in InterSpec. The primary objective of the NA-241 FY25 project was the detection and characterization of non-homogeneous uranium samples. However, the tool also demonstrates effectiveness in determining the enrichment of homogeneous uranium or plutonium from single gamma spectra. This paper presents a subset of evaluated data to avoid distribution limitations while offering potential users’ insight into the tool's performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Rapid Gamma Simulations of TRISO Fuel Elements

As energy demand rises, nuclear energy, particularly from reactors that use tristructural isotropic (TRISO) fuels, has gained attention due to the fuel’s enhanced resistance to radiation damage and high temperatures. This report investigates the modeling capabilities of the Gamma Detector Response and Analysis Software (GADRAS) for TRISO fuels, focusing on the gamma signatures of TRISO particles, which have not been extensively explored. Using the Monte Carlo N-Particle (MCNP) code as a benchmark, we developed both homogeneous and heterogeneous models of TRISO pebbles to analyze gamma spectra. Our findings reveal that the homogeneous and heterogeneous models produced different gamma signatures. Additionally, the GADRAS heterogeneous model significantly reduces computation times compared to MCNP, enabling effective modeling of gamma signatures for safeguards applications. This advancement is essential for the International Atomic Energy Agency (IAEA) in detecting anomalies and potential smuggling attempts in TRISO reactor fuel elements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Subplane Decusping for BWRs in MPACT

Control blade cusping can introduce significant error in boiling water reactor (BWR) calculations with MPACT when blade tips fall partway within an axial method of characteristics (MOC) plane, requiring homogenization of controlled and uncontrolled regions. This work implements subplane decusping for BWRs in MPACT by enabling BWR-compatible subplane coarse mesh finite difference (CMFD) and extending the existing decusping framework to represent between-assembly control blades that insert from the bottom of the core. The method resolves axial heterogeneity on a refined subplane mesh in the low-order solve and uses the resulting subplane flux shape to form flux-volume homogenized transport cross sections for the partially rodded MOC plane. The capability is evaluated using a single physics General Electric (GE)-14 assembly and a multiphysics Peach Bottom Unit 2 Type 1 assembly (PB2T1A) with thermal hydraulic feedback. In both cases, coarse axial meshes with and without subplane decusping are compared against fine mesh reference solutions over the full range of blade withdrawal positions. Subplane decusping reduces maximum/average eigenvalue errors from 3,369/220 pcm to 303/36 pcm for GE-14 and from 6,689/585 pcm to 1,582/109 pcm for PB2T1A. Additionally, it reduces pin power root mean square errors from 5.2%/0.8% to 1.8%/0.2% for GE-14 and from 6.9%/0.6% to 2.3%/0.1% for PB2T1A. These results demonstrate an effective, practical correction for BWR blade cusping in VERA-MPACT.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Mapping Spiking Neural Networks to Heterogeneous Crossbar Architectures using Integer Linear Programming

Advances in novel hardware devices and architectures allow Spiking Neural Network (SNN) evaluation using ultra-low power, mixed-signal, memristor crossbar arrays. As individual network sizes quickly scale beyond the dimensional capabilities of single crossbars, networks must be mapped onto multiple crossbars. Crossbar sizes within modern Memristor Crossbar Architectures (MCAs) are determined predominately not by device technology but by network topology; more, smaller crossbars consume less area thanks to the high structural sparsity found in larger, brain-inspired SNNs. Motivated by continuing increases in SNN sparsity due to improvements in training methods, we propose utilizing heterogeneous crossbar sizes to further reduce area consumption. This approach was previously unachievable as prior compiler studies only explored solutions targeting homogeneous MCAs. Our work improves on the state-of-the-art by providing Integer Linear Programming (ILP) formulations supporting arbitrarily heterogeneous architectures. By modeling axonal interactions between neurons, our methods produce better mappings while removing inhibitive a priori knowledge requirements. We first show a 16.7-27.6% reduction in area consumption for square-crossbar homogeneous architectures. Then, we demonstrate 66.9-72.7% further reduction when using a reasonable configuration of heterogeneous crossbar dimensions. Next, we present a new optimization formulation capable of minimizing the number of inter-crossbar routes. When applied to solutions already near-optimal in area, an 11.9-26.4% routing reduction is observed without impacting area consumption. Finally, we present a profile-guided optimization capable of minimizing the number of runtime spikes between crossbars. Compared to the best-area-then-route optimized solutions, we observe a further 0.5-14.8% inter-crossbar spike reduction while requiring 1–3 orders of magnitude less solver time.

Pohl, Devin [ORNL] (ORCID:0009000040149027)↗
DOE OSTI2026

Improving the Confidence in Retrievals of Vertical Distributions of Cloud Condensation Nuclei Number Concentration from ARM Supported by Aircraft In Situ Observations

Accurate quantification of the vertical distribution of cloud condensation nuclei (CCN) number concentrations is critical for improving our understanding of aerosol–cloud interactions. Ground-based Raman lidars operated by the Atmospheric Radiation Measurement (ARM) program, together with surface CCN measurements, are used to retrieve vertically resolved CCN number concentrations (Retrieved Number concentration of CCN, RNCCN). These retrievals rely on several assumptions, including that aerosol composition is vertically homogeneous. To assess this assumption, we developed and tested a framework to infer the dominant aerosol classes/types at different altitudes. This was done by applying a k-Nearest-Neighbors (kNN) algorithm to lidar ratio and linear depolarization ratio measurements from Raman lidar. We evaluated the framework using aircraft aerosol and CCN measurements from the ARM Holistic Interactions of Shallow Clouds, Aerosols, and Land Ecosystems (HI-SCALE) field campaign. The results show that RNCCN performance degrades as vertical aerosol complexity increases, i.e., RNCCN agrees with the aircraft CCN in vertically homogeneous conditions, but closure decreases in layered aerosol structures. To generalize beyond individual examples, we introduce a metric (heterogeneity index) that quantifies the vertical complexity by assessing the variation of inferred aerosol classes/types. Case-level statistics show a tendency for RNCCN and aircraft differences to increase with this metric. By detecting retrievals that are likely compromised by aerosol vertical heterogeneity, the proposed framework improves the interpretability and effective use of RNCCN used for long-term evaluation of models and aerosol–cloud interactions.

Tian, Jingjing↗
DOE OSTI2026

Investigation of onshore wind farm wake recovery with in situ aircraft measurements during AWAKEN

The share of wind power for electricity supply is increasing worldwide. This highly variable resource requires the improved prediction of power output for network stability. The interaction between wind farm wakes and the atmospheric boundary layer (ABL) introduces uncertainties in power production that warrant detailed investigation. The flow downwind of wind farms is characterized by a reduction in wind speed and an increase in turbulence, which both vary with atmospheric conditions. During the American WAKE experimeNt (AWAKEN), the Technische Universität Braunschweig conducted measurement flights with a research aircraft upwind and downwind of onshore wind farms in the southern Great Plains in Oklahoma in the USA. This study utilizes data from 20 flights conducted at approximately hub height in September 2023 to investigate the wind field variability downwind of the wind farms and vertical profiles to observe atmospheric stratification. The flights were aligned perpendicular to the main wind direction downwind of the King Plains and Armadillo Flats wind farms. Additionally, lidar data from both upwind and downwind ground-based measurement sites and sonic anemometer data were used for comprehensive analysis. Results indicate that under stable ABL conditions, the wake persists at greater downwind distances with a higher velocity deficit in the wake relative to the undisturbed flow compared to unstable stratification. In homogeneous terrain under stable conditions, wake recovery to 95 % occurs between a distance of 4.5 and 9 km downwind of the wind farm. In the semi-complex terrain characterized by shallow hills, slopes, and valleys, the wake exhibits a higher velocity deficit compared to homogeneous terrain, while in some cases the wake was amplified by the terrain resulting in higher velocity deficit 10 km downwind of the wind farm compared to the measurements closer to the wind farm. The turbulent kinetic energy (TKE) and “TKE difference” was found to be a valuable measure in understanding wakes in a semi-complex terrain, showing a clear wake recovery and formation depending on the stratification of the ABL.

17 WIND ENERGY↗
DOE OSTI2025

Evaluating Chemical Kinetics Predictions for Propane Using 3-D and 0-D Models in a Boosted Spark-Ignited Engine

Propane has been shown to be a promising alternative fuel to reduce emissions while simultaneously achieving high efficiencies in medium- and heavy-duty engines. These high-power density applications require boosted engines which, combined with high compression ratio, can lead to auto-ignition and knock. While three-dimensional (3-D) computational fluid dynamics (CFD) models are often used for resolving the complex fluid flow in engines, these models can become computationally expensive when simulating detailed chemical kinetics. Likewise, zero-dimensional (0-D) models are computationally concise enough for kinetics development, but lack any flow-field information which governs the flame propagation processes in spark ignition (SI) engines. This work presents a comprehensive comparison between 3-D and 0-D closed cycle simulations at knocking conditions in a high compression ratio high stroke-to-bore ratio propane engine. In order to initialize the flow-field for the 3-D closed cycle (intake valve closing, (IVC) to exhaust valve opening, (EVO)) simulation, a motored multi-cycle 3-D model was run using Converge to create a map at IVC, reducing the computational time. The map allowed a non-homogeneous 3-D closed cycle simulation to be satisfactorily validated against experiments, while a homogeneous case using only the turbulence field mapping was also simulated, mimicking 0-D modeling. The 3-D simulations were used to prescribe the initial conditions (e.g., IVC thermodynamics, speciation, burn-rate profile) for a 2-zone 0-D SI engine model in Chemkin Pro for both cases. It was found that 2-zone 0-D modeling underpredicted the knock onset timing, likely due to the lack of thermal stratification in the unburned gas region. Future work will carry multi-zone 0-D modeling to capture the fuel auto-ignition in the unburned region.

Douvry-Rabjeau, Julien [Oakland University, Roches↗
DOE OSTI2025

Micropolar Elastoplasticity Using a Fast Fourier Transform‐Based Solver

ABSTRACT This work presents a micromechanical spectral formulation for obtaining the full‐field and homogenized response of elastoplastic micropolar composites. A closed‐form radial‐return mapping is derived from thermodynamics‐based micropolar elastoplastic constitutive equations to determine the increment of plastic strain necessary to return the generalized stress state to the yield surface, and the algorithm implementation is verified using the method of numerically manufactured solutions. Then, size‐dependent material response and micro‐plasticity are shown as features that may be efficiently simulated in this micropolar elastoplastic framework. The computational efficiency of the formulation enables the generation of large datasets in reasonable computing times.

42 ENGINEERING↗
DOE OSTI2024

Enforcing global constraints for the dispersion closure problem: τ 2 -SIMPLE algorithm

Permeability and effective dispersion tensors are critical parameters to characterize flow and transport in porous media at the continuum scale. Homogenization theory defines a framework in which such effective properties are first computed from solving a closure problem in a repeating unit cell of the periodic microstructure and then used in a macroscopic formulation for efficient computation. The closure problem is formulated as a local boundary value problem subjected to global constraints, which guarantee the uniqueness of the solution and can be difficult to satisfy for complex geometries and at high flow conditions. These constraints also ensure that pore-scale pressure, velocity, and concentration fields can be accurately reconstructed from the closure variable. Building on a previous work, here we present a framework that allows to satisfy global constraints associated to both the permeability and the dispersion closure problems by introducing two artificial time scales. The algorithm, called τ 2 -SIMPLE, computes both permeability and effective dispersion given an arbitrarily complex geometry and flow condition. Furthermore, this algorithm is demonstrated to be accurate for both 2D and 3D geometries across varying flow conditions, and thus it can be used to quickly characterize effective properties from porous media images in many applications.

97 MATHEMATICS AND COMPUTING↗