Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient Limiters”

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 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Experiments in Neural-Network Control of a Free-Flying Space Robot

Four important generic issues are identified and addressed in some depth in this thesis as part of the development of an adaptive neural network based control system for an experimental free flying space robot prototype. The first issue concerns the importance of true system level design of the control system. A new hybrid strategy is developed here, in depth, for the beneficial integration of neural networks into the total control system. A second important issue in neural network control concerns incorporating a priori knowledge into the neural network. In many applications, it is possible to get a reasonably accurate controller using conventional means. If this prior information is used purposefully to provide a starting point for the optimizing capabilities of the neural network, it can provide much faster initial learning. In a step towards addressing this issue, a new generic Fully Connected Architecture (FCA) is developed for use with backpropagation. A third issue is that neural networks are commonly trained using a gradient based optimization method such as backpropagation; but many real world systems have Discrete Valued Functions (DVFs) that do not permit gradient based optimization. One example is the on-off thrusters that are common on spacecraft. A new technique is developed here that now extends backpropagation learning for use with DVFs. The fourth issue is that the speed of adaptation is often a limiting factor in the implementation of a neural network control system. This issue has been strongly resolved in the research by drawing on the above new contributions.

NEURAL NETS↗
NASA NTRS2017

Vertical Profiles of Mars 1.27 Micrometer O2 Dayglow from MRO CRISM Limb Spectra: Seasonal/global Behaviors, Comparisons to LMDGCM Simulations, and a Global Definition for Mars Water Vapor Profiles

Since July of 2009, The Compact Reconnaissance Imaging Spectral Mapper (CRISM) onboard the Mars Re- connaissance Orbiter (MRO) has periodically obtained pole-to-pole observations (i.e., full MRO orbits) of limb scanned visible/near IR spectra (lambda= 0.4 - 4.0 micrometers, delta lambda approx. 10 nm- Murchie et al., 2007). These CRISM limb observations support the first seasonally and spatially extensive set of Mars 1.27 micometers O2 (1 delta(sub g)) day-glow profile retrievals (approx. 1100) over greater than or equal to 8-80 km altitudes. Their comparison to Laboratoire de Meteorologie Dynamique (LMD) global climate model (GCM) simulated O2 (1 delta(sub g)) volume emission rate (VER) profiles, as a function of altitude, latitude, and season (solar longitude, L(sub s), supports several key conclusions regarding Mars atmospheric water vapor (which is derived from O2 (1 delta(sub g)) emission rates), Mars O3, and the collisional de-excitation of O2 (1 delta(sub g)) in the Mars CO2 atmosphere. Current (Navarro et al., 2014) LMDGCM simulations of Mars atmospheric water vapor fall 2-3 times below CRISM derived water vapor abundances at 20-40 km altitudes over low-to-mid latitudes in northern spring (L(sub s) = 30-60 deg), and northern mid-to-high latitudes over northern summer (L(sub s) = 60-140 deg). In contrast, LMDGCM simulated water vapor is 2-5 times greater than CRISM derived values at all latitudes and seasons above 40 km, within the aphelion cloud belt (ACB), and over high-southern to mid-southern latitudes in southern summer (L(sub s) = 190-340 deg) at 15-35 km altitudes. Overall, the solstitial summer-to-winter hemisphere gradients in water vapor are reversed between the LMDGCM modeled versus the CRISM derived water vapor abundances above 10-30 km altitudes. LMDGCM-CRISM differences in water vapor profiles correlate with LMDGCM-CRISM differences in cloud mixing profiles; and likely reflect limitations in simulating cloud microphysics and radiative forcing, both of which restrict meridional transport of water from summer- to-winter hemispheres on Mars (Clancy et al., 1996; Montmessin et al., 2004; Steele et al., 2014; Navarro et al., 2014) and depend on uncertain cloud microphysical properties (Navarro et al., 2014). The derived low-to-mid latitude changes in Mars water vapor vertical distributions should reduce current model- data disagreements in column O3 and H2O2 abundances over low-to-mid latitudes (e.g., within the ACB; Lefevre et al., 2008; Encrenaz et al., 2015; Clancy et al., 2016). Lastly, the global/seasonal average com- parison of CRISM and LMDGCM O2 (1 delta(sub g)) VER below 20 km altitudes indicates a factor of approx. 3 times lower value (0.25 x 10(exp -20) cu cm sec(exp -1)) for the CO2 collisional de-excitation rate coefficient of O2 (1 delta(sub g)) than derived recently by Guslyakova et al. (2016).

Clancy, R. Todd↗
NASA NTRS2018

The Central Role of Air Quality Observations in NASA's GEOS Composition Forecasting Model

The NASA GEOS composition forecast model (GEOS-CF) provides global, high-resolution (25 km) air quality forecasts in near-real time. This system combines the operational GEOS-5 weather forecasting model with the state-of-the-science GEOS-Chem chemistry module to provide detailed chemical analysis of a wide range of air pollutants such as ozone, carbon monoxide, nitrogen oxides, and fine particulate matter (PM2.5). The resolution of the forecasts is the highest compared to current, publicly-available global composition forecasts.Air quality observations are an indispensable tool to evaluate the model's ability to capture the strong temporal and spatial gradients of air pollutants across the globe. We show how comparisons against near-real time observations available through OpenAQ (www.openaq.org) demonstrate the model's overall success in reproducing surface concentrations of ozone, nitrogen dioxide, and PM2.5. This analysis also helps identifying current limitations of the model, for example over South America. The model-observation mismatches are most likely caused by uncertainties in the emissions data. Using the example of Rio de Janeiro, we show how the model skill can be improved by using local, high-resolution emission inventories in combination with air quality data.

Keller, Christoph A.↗
NASA NTRS2004

Adjoint Formulation for an Embedded-Boundary Cartesian Method

Many problems in aerodynamic design can be characterized by smooth and convex objective functions. This motivates the use of gradient-based algorithms, particularly for problems with a large number of design variables, to efficiently determine optimal shapes and configurations that maximize aerodynamic performance. Accurate and efficient computation of the gradient, however, remains a challenging task. In optimization problems where the number of design variables dominates the number of objectives and flow- dependent constraints, the cost of gradient computations can be significantly reduced by the use of the adjoint method. The problem of aerodynamic optimization using the adjoint method has been analyzed and validated for both structured and unstructured grids. The method has been applied to design problems governed by the potential, Euler, and Navier-Stokes equations and can be subdivided into the continuous and discrete formulations. Giles and Pierce provide a detailed review of both approaches. Most implementations rely on grid-perturbation or mapping procedures during the gradient computation that explicitly couple changes in the surface shape to the volume grid. The solution of the adjoint equation is usually accomplished using the same scheme that solves the governing flow equations. Examples of such code reuse include multistage Runge-Kutta schemes coupled with multigrid, approximate-factorization, line-implicit Gauss-Seidel, and also preconditioned GMRES. The development of the adjoint method for aerodynamic optimization problems on Cartesian grids has been limited. In contrast to implementations on structured and unstructured grids, Cartesian grid methods decouple the surface discretization from the volume grid. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin e t al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation.

Nemec, Marian↗
DOE OSTI2016

Quasi-linear gyrokinetic predictions of the Coriolis momentum pinch in National Spherical Torus Experiment

This paper presents quasi-linear gyrokinetic predictions of the Coriolis momentum pinch for low aspect-ratio National Spherical Torus Experiment (NSTX) H-modes where previous experimental measurements were focused. Local, linear calculations predict that in the region of interest (just outside the mid-radius) of these relatively high-beta plasmas, profiles are most unstable to microtearing modes that are only effective in transporting electron energy. However, sub-dominant electromagnetic and electrostaticballooning modes are also unstable, which are effective at transporting energy, particles, and momentum. The quasi-linear prediction of transport from these weaker ballooning modes, assuming they contribute transport in addition to that from microtearing modes in a nonlinear turbulent state, leads to a very small or outward convection of momentum, inconsistent with the experimentally measured inward pinch, and opposite to predictions in conventional aspect ratio tokamaks. Additional predictions of a low beta L-mode plasma, unstable to more traditional electrostatic ion temperature gradient-trapped electron mode instability, show that the Coriolis pinch is inward but remains relatively weak and insensitive to many parameter variations. The weak or outward pinch predicted in NSTX plasmas appears to be at least partially correlated to changes in the parallel mode structure that occur at a finite beta and low aspect ratio, as discussed in previous theories. The only conditions identified where a stronger inward pinch is predicted occur either in the purely electrostatic limit or if the aspect ratio is increased. Lastly, as the Coriolis pinch cannot explain the measured momentum pinch, additional theoretical momentum transport mechanisms are discussed that may be potentially important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Numerical modeling of heat transfer in molten silicon during directional solidification process

Numerical investigation is performed for some of the thermal and fluid flow properties of silicon melt during directional solidification by numerical modeling. Dimensionless numbers are extremely useful to understand the heat and mass transfer of fluid flow on Si melt and control the flow patterns during crystal growth processes. The average grain size of whole crystal would increase when the melt flow is laminar. In the silicon growth process, the melt flow is mainly driven by the buoyancy force resulting from the horizontal temperature gradient. The thermal and flow pattern influences the quality of the crystal through the convective heat and mass transport. The computations are carried out in a 2D axisymmetric model using the finite-element technique. The buoyancy effect is observed in the melt domain for a constant Rayleigh number and for different Prandtl numbers. The convective heat flux and Reynolds numbers are studied in the five parallel horizontal cross section of melt silicon region. And also, velocity field is simulated for whole melt domain with limited thermal boundaries. The results indicate that buoyancy forces have a dramatic effect on the most of melt region except central part.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Counter-Current Flow Limitation Studies in Complex Geometries Utilizing Interface Capturing Simulations Coupled with PID Flow Rate Controller

In nuclear thermal-hydraulic studies, counter-current flow limitation (CCFL) typically refers to steam rising at a fast rate such that it prevents coolant from draining down within a confined channel. CCFL is a crucial issue in nuclear reactor safety analysis. This study investigates CCFL in debris bed channels using high-resolution interface-capturing simulations. A novel proportional-integral-derivative flow rate controller is developed to efficiently achieve the CCFL conditions. Verification studies confirm that CCFL occurs under the same conditions with or without the controller, demonstrating that PID control ensures accurate prediction. Three debris bed channel geometries were examined: a cylindrical channel, a channel with small obstacles, and a channel with large obstacles. Results show that obstacles significantly impact flow behavior, interfacial shear, wall shear, and pressure gradients required for CCFL. Furthermore, the comparison with experimental data confirmed that simulations incorporating geometric complexities align more closely with experimental CCFL conditions. A pressure gradient correlation was also developed for CCFL prediction.

Counter-current flow limitation↗
DOE OSTI2012

Special Analysis for Slit Trench Disposal of the Reactor Process Heat Exchangers

The Savannah River National Laboratory (SRNL), in response to a request from Solid Waste Management (SWM), conducted a Special Analysis (SA) to evaluate the performance of nineteen heat exchangers that are to be disposed in the E-Area low level waste facility Slit Trench 9 (ST 9). Although these nineteen heat exchangers were never decontaminated, the majority of the radionuclides in the heat exchanger inventory list were determined to be acceptable for burial because they are less than the 'generic' waste form inventory limits given in the 2008 Performance Assessment (PA) (WSRC, 2008). However, as generic waste, the H-3 and C-14 inventories resulted in unacceptable sum-of-fractions (SOFs). Initial scoping analyses performed by SRNL indicated that if alterations were made to certain external nozzles to mitigate various potential leak paths, acceptable SOFs could be achieved through the use of a 'Special' waste form. This SA provides the technical basis for this new 'Special' waste form and provides the inventory limits for H-3 and C-14 for these nineteen heat exchangers such that the nineteen heat exchangers can be disposed in ST 9. This 'Special' waste form is limited to these nineteen heat exchangers in ST 9 and applies for H-3 and C-14, which are designated as H-3X and C-14X, respectively. The SA follows the same methodology used in the 2008 PA and the 2008 SA except for the modeling enhancements noted below. Infiltration rates above the heat exchangers are identical to those used in the 2008 PA; however, flow through the heat exchangers is unique. Because it is unknown exactly how sealed heat exchanger openings will perform and how surface and embedded contaminants will be released, multiple base cases or scenarios were established to investigate a set of performances. Each scenario consists of flow options (based on the performance of sealed openings) and a near-field release of contaminants (based on corrosion and diffusion performance). Two disposal configurations were analyzed where heat exchangers were assumed to be disposed four across and five lengthwise (the 4x5 configuration, with one empty) and three across and seven lengthwise (the 3x7 configuration, with two empty). A large range of conditions was considered. For example, peak well concentrations at the 100-m boundary for H-3 are shown in Figure ES-1 for a wide range of configurations (i.e. release mechanism and degree of sealing options). The maximum contaminant level (MCL) and a 10% SOF goal for H-3 are also shown. The 10% goal was based on an estimated volume fraction that these nineteen heat exchangers would consume in ST 9 and was solely used for scoping purposes to assess disposal feasibility and sealing requirements. Because various line breaks and poor sealing greatly exceeded that 10% goal, the determination was made that mitigating activities were needed, such as protection from line breaks and better sealing. An initial set of scenarios was run to assess the requirements for sealing the heat exchanger openings and the need to ensure that the sealed heat exchangers stayed sealed during transit and disposal operations. After discovering that such mitigating activities were required, additional scenarios were run that included the mitigating activities. Scenarios deemed to have a very low probability of occurrence were excluded from consideration for calculating inventory limits (for example, those scenarios that assumed an instantaneous release of contaminants along with poor sealing). The SA used the most recent K d values for the C-14 analyses and the most recent Dose Conversion Factors for H-3 and C-14 which have been updated since the 2008 PA was issued. This SA took into account the location and the disposal timing of these heat exchangers. The disposal location is within a small area of the overall Slit Trench unit (about 6% of the total) and is behind a line that is 200 ft from the down-gradient edge of ST 9. The disposal timing is assumed to be after July 1, 2012 (because disposals cannot occur until this document is approved and mitigating activities are completed) which means that the disposal occurs after the first time period for the 2008 PA beta-gamma pathway (that time period is from December 1995 until December 2007), thus that pathway time period is not considered. Table ES-1 provides new 'Special' waste form groundwater pathway inventory limits for C-14X and H-3X in the heat exchangers. Inventory limits for generic C-14 and H-3 in the West Slit Trenches are included for comparison. The lowest limit for generic C-14 is 1.9E-1 Ci, while for C-14X it is 2.7E0, an increase of more than 14 times. Because time windows are employed, at later times C-14X exhibits lower limits than those for generic C-14 because with its smaller K d the C-14 moves much faster. The lowest limit for generic H-3 is 3.6E0 Ci, while for H-3X it is 1.7E3, an increase of almost 500 times.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRSDate not supplied

Multiple Probe Measurements at Uranus Motivated By Spatial Variability

Motivation: Spatial variations in the temperature field and composition of Uranus' atmosphere demonstrate a need for multiple entry probes to characterize vertical profiles in multiple locations. We will review variation of composition and temperature, which are produced by dynamical processes on a range of scales from global (polar anomalies, zonal bands) to regional (vortices, storms) [1–2]. In particular, the spatial variation of convective activity is not well understood based on existing remote sensing observations [3]. Understanding how these processes operate, and how they modulate variable composition, is key to constraining bulk atmospheric abundances. Abundances in turn provide cosmochemical constraints on planetary origins. Secondary probes at Uranus: Key measurements for secondary probes are temperature-pressure profiles, along with compositional profiles. Radio occultations are limited to shallow levels less than ~2 bar [4]. In situ temperature-pressure measurements (or atmospheric structure measurements) can unambiguously extend these results to deeper levels. Measurements with vertical resolution of at least 2 km are needed to characterize anomalies like the 1.2-bar feature from the Voyager 2 occultation at 2–6° S [4], which supports a range of temperature gradients depending on assumptions of composition. Simultaneous measurements of temperature, pressure, and composition are key to understanding this class of features. Species such as methane, hydrogen sulfide, and ammonia must be measured along a probe descent profile because they vary over several orders of magnitude due to the strong temperature dependence of their saturation vapor pressures [e.g., 5]. Results of these measurements can be interpreted to understand the potential for moist convective activity in the atmosphere. Lessons from the other giant planets: On Jupiter, the Galileo Probe's entry into a meteorologically distinct five-micron hot spot led many to interpret the local composition as column-stretched, so that well-mixed abundances were reached at deeper levels than in surrounding, unperturbed regions [6–8]. Ground-based microwave measurements and Juno data now indicate that the deep depletion of ammonia is a very widespread atmospheric characteristic not limited to 5-µm hot spots [9–11]. But in the absence of multiple probes, we do not know if the other volatiles H2S and H2O behave in the same way. The open questions for the Jupiter case strongly motivate sending multiple probes to Uranus. On Saturn, retrievals of NH3 and PH3 abundances at shallow levels vary with latitude [12], but it is unknown how deep these differences extend, which is why a Saturn probe (or probes) was a mission theme for NASA's New Frontiers 4 and 5 opportunities [13]. Challenges for secondary probes: Cost is an issue due to the perception that it involves sacrifices to other mission elements. Spacecraft trajectories may be constrained by needs for orbit insertion that limit probe deliveries to different latitudes, and additional limitations may be placed on communication windows for probe descent phases (particularly if multiple probes are released from the orbiter simultaneously) [14]. Finally, composition sensors for miniature secondary probes are not at the required technological maturity [15]. Mass spectrometers are typically too large, massive, and powerhungry, while smaller nanosensors are only beginning to be developed for planetary missions. Finally, the need for probe survival heating is most easily met by radioisotope heat sources, but these require regulatory approvals that are even more difficult to satisfy compared to standard environmental reviews [16], unless the secondary probe is designed as a core element of a mission. References: [1] Molter E.M., et al. (2021) PSJ, 2, 3. [2] Rowe-Gurney N., et al. (2021) Icar, 365, 114506. [3] Hueso R., Sánchez-Lavega A. (2019) SSRv, 215, 52. [4] Lindal G.F., et al. (1987) JGR, 92, 14987– 15001. [5] Simon A.A., et al. (2022) RemS, 14, 1518. [6] Atreya, S.K., et al. (1997) in The Three Galileos: The Man, the Spacecraft, the Telescope, pp. 249–260 (C. Barbieri et al., eds.). [7] Showman A.P., Ingersoll A.P. (1998) Icar, 132, 205–220. [8] Friedson A.J. (2005) Icar, 177, 1–17. [9] de Pater I., et al. (2001) Icar, 149, 66–78. [10] Li C., et al. (2017) GeoRL, 44, 5317–5325. [11] de Pater I., et al. (2019) Icar, 322, 168–191. [12] Fletcher L.N., et al. (2009) Icar, 202, 543–564. [13] National Research Council (2011) Vision and Voyages. [14] Sayanagi K.M., et al. (2020) SSRv, 216, 72. [15] Wong M.H. et al. (2021) BAAS, 53, 486. [16] Zide A., Mendoza-Hill A., Cheney D. (2022) COSPAR Abstracts H0.6-0012-22.

Michael H Wong↗
DOE OSTI2023

Colloidal State Machines as Smart Tracers for Chemical Reactor Analysis

A widely utilized tool in reactor analysis is passive tracers that report the residence time distribution, allowing estimation of the conversion and other properties of the system. Recently, advances in microrobotics have introduced powered and functional entities with sizes comparable to some traditional tracers. This has motivated the concept of Smart Tracers that could record the local chemical concentrations, temperature, or other conditions as they progress through reactors. Herein, the design constraints and advantages of Smart Tracers by simulating their operation in a laminar flow reactor model conducting chemical reactions of various orders are analyzed. It is noted that far fewer particles are necessary to completely map even the most complex concentration gradients compared with their conventional counterparts. Design criteria explored herein include sampling frequency, memory storage capacity, and ensemble number necessary to achieve the required accuracy to inform a reactor model. Cases of severe particle diffusion and sensor noise appear to bind the functional upper limit of such probes and require consideration for future design. The results of the study provide a starting framework for applying the new technology of microrobotics to the broad and impactful set of problems classified as chemical reactor analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

On the edge energy of lipid membranes and the thermodynamic stability of pores

To perform its barrier function, the lipid bilayer membrane requires a robust resistance against pore formation. Using a self-consistent field (SCF) theory and a molecularly detailed model for membranes composed of charged or zwitterionic lipids, it is possible to predict structural, mechanical, and thermodynamical parameters for relevant lipid bilayer membranes. We argue that the edge energy in membranes is a function of the spontaneous lipid monolayer curvature, the mean bending modulus, and the membrane thickness. An analytical Helfrich-like model suggests that most bilayers should have a positive edge energy. This means that there is a natural resistance against pore formation. Edge energies evaluated explicitly in a two-gradient SCF model are consistent with this. Remarkably, the edge energy can become negative for phosphatidylglycerol (e.g., dioleoylphosphoglycerol) bilayers at a sufficiently low ionic strength. Such bilayers become unstable against the formation of pores or the formation of lipid disks. In the weakly curved limit, we study the curvature dependence of the edge energy and evaluate the preferred edge curvature and the edge bending modulus. The latter is always positive, and the former increases with increasing ionic strength. These results point to a small window of ionic strengths for which stable pores can form as too low ionic strengths give rise to lipid disks. Higher order curvature terms are necessary to accurately predict relevant pore sizes in bilayers. The electric double layer overlap across a small pore widens the window of ionic strengths for which pores are stable.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

A multi-channel capacitive probe for electrostatic fluctuation measurement in the Madison Symmetric Torus reversed field pinch

A 40-channel capacitive probe has been developed to measure the electrostatic fluctuations associated with the tearing modes deep into Madison Symmetric Torus (MST) reversed field pinch plasma. The capacitive probe measures the ac component of the plasma potential via the voltage induced on stainless steel electrodes capacitively coupled with the plasma through a thin annular layer of boron nitride (BN) dielectric (also serves as the particle shield). When bombarded by the plasma electrons, BN provides a sufficiently large secondary electron emission for the induced voltage to be very close to the plasma potential. The probe consists of four stalks each with ten cylindrical capacitors that are radially separated by 1.5 cm. The four stalks are arranged on a 1.3 cm square grid so that at each radial position, there are four electrodes forming a square grid. Every two adjacent radial sets of four electrodes form a cube. The fluctuating electric field can be calculated by the gradient of the plasma potential fluctuations at the eight corners of the cube. The probe can be inserted up to 15 cm (r/a = 0.7) into the plasma. Furthermore, the capacitive probe has a frequency bandwidth from 13 Hz to 100 kHz, amplifier-circuit limit, sufficient for studying the tearing modes (5-30 kHz) in the MST reversed-field pinch.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Saccharomycotina yeasts defy long-standing macroecological patterns

The Saccharomycotina yeasts (“yeasts” hereafter) are a fungal clade of scientific, economic, and medical significance. Yeasts are highly ecologically diverse, found across a broad range of environments in every biome and continent on earth; however, little is known about what rules govern the macroecology of yeast species and their range limits in the wild. Here, we trained machine learning models on 12,816 terrestrial occurrence records and 96 environmental variables to infer global distribution maps at ~1 km2 resolution for 186 yeast species (~15% of described species from 75% of orders) and to test environmental drivers of yeast biogeography and macroecology. We found that predicted yeast diversity hotspots occur in mixed montane forests in temperate climates. Diversity in vegetation type and topography were some of the greatest predictors of yeast species richness, suggesting that microhabitats and environmental clines are key to yeast diversity. We further found that range limits in yeasts are significantly influenced by carbon niche breadth and range overlap with other yeast species, with carbon specialists and species in high-diversity environments exhibiting reduced geographic ranges. Finally, yeasts contravene many long-standing macroecological principles, including the latitudinal diversity gradient, temperature-dependent species richness, and a positive relationship between latitude and range size (Rapoport’s rule). These results unveil how the environment governs the global diversity and distribution of species in the yeast subphylum. These high-resolution models of yeast species distributions will facilitate the prediction of economically relevant and emerging pathogenic species under current and future climate scenarios.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Mechanisms for the clustering of inertial particles in the inertial range of isotropic turbulence

In this study, we consider the physical mechanism for the clustering of inertial particles in the inertial range of isotropic turbulence. We analyze the exact, but unclosed, equation governing the radial distribution function (RDF) and compare the mechanisms it describes for clustering in the dissipation and inertial ranges. We demonstrate that in the limit St r <<1, where St r is the Stokes number based on the eddy turnover time scale at separation r, the clustering in the inertial range can be understood to be due to the preferential sampling of the coarse-grained fluid velocity gradient tensor at that scale. When St r ≳O(1) this mechanism gives way to a nonlocal clustering mechanism. These findings reveal that the clustering mechanisms in the inertial range are analogous to the mechanisms that we identified for the dissipation regime. Further, we discuss the similarities and differences between the clustering mechanisms we identify in the inertial range and the “sweep-stick” mechanism developed by Coleman and Vassilicos. We show that the idea that initial particles are swept along with acceleration stagnation points is only approximately true because there always exists a finite difference between the velocity of the acceleration stagnation points and the local fluid velocity. This relative velocity is sufficient to allow particles to traverse the average distance between the stagnation points within the correlation time scale of the acceleration field. We also show that the stick part of the mechanism is only valid for St r <<1 in the inertial range. We emphasize that our clustering mechanism provides the more fundamental explanation since it, unlike the sweep-stick mechanism, is able to explain clustering in arbitrary spatially correlated velocity fields. We then consider the closed, model equation for the RDF given in Zaichik and Alipchenkov and use this, together with the results from our analysis, to predict the analytic form of the RDF in the inertial range for St r <<1, which, unlike that in the dissipation range, is not scale invariant. Finally, the results are in good agreement with direct numerical simulations, provided the separations are well within the inertial range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Methane production and emissions in trees and forests

Forest ecosystem methane (CH 4 ) research has focused on soils, but trees are also important sources and sinks in forest CH 4 budgets. Living and dead trees transport and emit CH 4 produced in soils; living trees and dead wood emit CH 4 produced inside trees by microorganisms; and trees produce CH 4 through an abiotic photochemical process. Here, we review the state of the science on the production, consumption, transport, and emission of CH 4 by living and dead trees, and the spatial and temporal dynamics of these processes across hydrologic gradients inclusive of wetland and upland ecosystems. Emerging research demonstrates that tree CH 4 emissions can significantly increase the source strength of wetland forests, and modestly decrease the sink strength of upland forests. Scaling from stem or leaf measurements to trees or forests is limited by knowledge of the mechanisms by which trees transport soil-produced CH 4 , microbial processes produce and oxidize CH 4 inside trees, a lack of mechanistic models, the diffuse nature of forest CH 4 fluxes, complex overlap between sources and sinks, and extreme variation across individuals. Understanding the complex processes that regulate CH 4 source–sink dynamics in trees and forests requires cross-disciplinary research and new conceptual models that transcend the traditional binary classification of wetland vs upland forest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

DP-TwoLevel: two-stage gradient subspace learning for differentially private federated learning

Federated learning (FL) enables collaborative model training across distributed data sources without sharing raw data, but faces fundamental challenges in communication efficiency and privacy. Differentially private (DP) training mitigates information leakage but introduces noise that degrades model performance, especially in high-dimensional settings. We propose DP-TwoLevel, a hierarchical gradient projection method that improves utility under fixed DP constraints by exploiting low-dimensional structure in model updates. Our approach learns a two-level PCA-based representation of gradients and applies DP noise in a reduced-dimensional subspace, thereby lowering the effective noise magnitude while preserving dominant signal components. We evaluate the method across three datasets (MNIST, Fashion-MNIST, CIFAR-10) and three privacy regimes (ϵ∈0.5, 1.0, 2.0). Across nine experimental settings, DP-TwoLevel consistently outperforms DP-FedAvg, achieving an average accuracy improvement of 9.44%, with larger gains observed in lower ϵ(higher-noise) regimes (up to +22.31%). We further analyze scalability across models ranging from 100K to 1.49M parameters and identify a variance-based success criterion: performance remains strong when the projection preserves more than 75% of gradient variance, degrades in a marginal regime (65–75%), and fails below this threshold. Our results demonstrate that structure-aware dimensionality reduction can significantly improve the privacy–utility tradeoff in FL without modifying formal privacy guarantees. We also provide empirical evidence of scaling limitations for global projections and motivate per-layer extensions for larger models.

Kotevska, Olivera [ORNL] (ORCID:0000000316772243)↗
DOE OSTI2021

Can Simple Machine Learning Tools Extend and Improve Temperature-Based Methods to Infer Streambed Flux?

Temperature-based methods have been developed to infer 1D vertical exchange flux between a stream and the subsurface. Current analyses rely on fitting physically based analytical and numerical models to temperature time series measured at multiple depths to infer daily average flux. These methods have seen wide use in hydrologic science despite strong simplifying assumptions including a lack of consideration of model structural error or the impacts of multidimensional flow or the impacts of transient streambed hydraulic properties. We performed a “perfect-model experiment” investigation to examine whether regression trees, with and without gradient boosting, can extract sufficient information from model-generated subsurface temperature time series, with and without added measurement error, to infer the corresponding exchange flux time series at the streambed surface. Using model-generated, synthetic data allowed us to assess the basic limitations to the use of machine learning; further examination of real data is only warranted if the method can be shown to perform well under these ideal conditions. We also examined whether the inherent feature importance analyses of tree-based machine learning methods can be used to optimize monitoring networks for exchange flux inference.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1983

Cumulus clouds - Numerical models, observations and entrainment

The first computer simulation of the organization phase of a buoyant atmospheric thermal is described. Although crude, it showed the spontaneous development of a rounded tight-gradient 'cap' and internal vortical circulation. The complexities involved in these 'field of motion' models in part motivated the development of entity models, based upon laboratory thermals. These one-dimensional models and their uses with observations are briefly described as well as their limitations. Finally, an application of Schlesinger's three-dimensional model to a GATE cumulus situation clarifies many apparently conflicting observations and postulates, thereby raising further challenging questions to be addressed jointly by the more sophisticated measuring and modeling tools available in the 1980's.

Simpson, J.↗