Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “root dynamics”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Plant myo -inositol transport influences bacterial colonization phenotypes

Plant microbiomes are assembled and modified through a complex milieu of biotic and abiotic factors. Despite dynamic and fluctuating contributing variables, specific host metabolites are consistently identified as important mediators of microbial interactions. We combine information from a large-scale metatranscriptomic dataset from natural poplar trees and experimental genetic manipulation assays in seedlings of the model plant Arabidopsis thaliana to converge on a conserved role for transport of the plant metabolite myo-inositol in mediating host-microbe interactions. While microbial catabolism of this compound has been linked to increased host colonization, we identify bacterial phenotypes that occur in both catabolism-dependent and -independent manners, suggesting that myo-inositol may additionally serve as a eukaryotic-derived signaling molecule to modulate microbial activities. Our data suggest host control of this compound and resulting microbial behavior are important mechanisms at play surrounding the host metabolite myo-inositol.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Enforcing detailed balance in the Borgnakke–Larsen redistribution method with temperature dependent relaxation models

For decades, it has been observed that the commonly used Borgnakke–Larsen method for energy redistribution in Direct Simulation Monte Carlo codes fails to satisfy the principle of detailed balance when coupled to a wide variety of temperature dependent relaxation models, while seemingly satisfying detailed balance when coupled to others. Many attempts have been made to remedy the issue, yet much ambiguity remains, and no consensus appears in the literature regarding the root cause of the intermittent compatibility of the Borgnakke–Larsen method with temperature dependent relaxation models. This paper alleviates that ambiguity by presenting a rigorous theoretical derivation of the Borgnakke–Larsen method's requirement for satisfying detailed balance. Specifically, it is shown that the Borgnakke–Larsen method maintains detailed balance if and only if the probability of internal-energy exchange during a collision depends only on collision invariants (e.g., total energy). Here, the consequences of this result are explored in the context of several published definitions of relaxation temperature, including translational, total, and cell-averaged temperatures. Of particular note, it is shown that cell-averaged temperatures, which have been widely discussed in the literature as a way to ensure equilibrium is reached, also fail in a similar, although less dramatic, fashion when the aforementioned relationship is not enforced. The developed theory can be used when implementing existing or new relaxation models and will ensure that detailed balance is satisfied.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRSDate not supplied

Systematic Errors in Simulated L-Band Brightness Temperature in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission has been providing L-band (1.4 GHz) brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter (EnKF), the SMAP Level-4 Soil Moisture (L4_SM) product provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications. The EnKF-based L4_SM analysis assumes unbiased forecast errors. Consequently, the seasonally varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations. The L4_SM system is thus designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (possibly wrong) climatology. In this paper, we examine the Tb observation-minus-forecast (O-F) residuals from the L4_SM Version 7 product (computed after rescaling the Tb observations to the mean seasonal cycle of the simulated Tb). The long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. The model forecast Tb, however, still exhibits undesirable systematic errors relative to the (rescaled) Tb observations. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to improve the L4_SM algorithm.

Rolf Reichle↗
DOE OSTI2020

DSCOVR/EPIC-derived global hourly and daily downward shortwave and photosynthetically active radiation data at 0.1° × 0.1° resolution

Downward shortwave radiation (SW) and photosynthetically active radiation (PAR) play crucial roles in Earth system dynamics. Spaceborne remote sensing techniques provide a unique means for mapping accurate spatiotemporally continuous SW–PAR, globally. However, any individual polar-orbiting or geostationary satellite cannot satisfy the desired high temporal resolution (sub-daily) and global coverage simultaneously, while integrating and fusing multisource data from complementary satellites/sensors is challenging because of co-registration, intercalibration, near real-time data delivery and the effects of discrepancies in orbital geometry. The Earth Polychromatic Imaging Camera (EPIC) on board the Deep Space Climate Observatory (DSCOVR), launched in February 2015, offers an unprecedented possibility to bridge the gap between high temporal resolution and global coverage and characterize the diurnal cycles of SW–PAR globally. In this study, we adopted a suite of well-validated data-driven machine-learning models to generate the first global land products of SW–PAR, from June 2015 to June 2019, based on DSCOVR/EPIC data. The derived products have high temporal resolution (hourly) and medium spatial resolution (0.1°×0.1°), and they include estimates of the direct and diffuse components of SW–PAR. We used independently widely distributed ground station data from the Baseline Surface Radiation Network (BSRN), the Surface Radiation Budget Network (SURFRAD), NOAA's Global Monitoring Division and the U.S. Department of Energy's Atmospheric System Research (ASR) program to evaluate the performance of our products, and we further analyzed and compared the spatiotemporal characteristics of the derived products with the benchmarking Clouds and the Earth's Radiant Energy System Synoptic (CERES) data. We found both the hourly and daily products to be consistent with ground-based observations (e.g., hourly and daily total SWs have low biases of -3.96 and -0.71 W m -2 and root-mean-square errors (RMSEs) of 103.50 and 35.40 W m -2 , respectively). The developed products capture the complex spatiotemporal patterns well and accurately track substantial diurnal, monthly, and seasonal variations in SW–PAR when compared to CERES data. They provide a reliable and valuable alternative for solar photovoltaic applications worldwide and can be used to improve our understanding of the diurnal and seasonal variabilities of the terrestrial water, carbon and energy fluxes at various spatial scales.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2005

MESSENGER: Exploring Mercury's Magnetosphere

The MESSENGER mission to Mercury offers our first opportunity to explore this planet s miniature magnetosphere since the brief flybys of Mariner 10. Mercury s magnetosphere is unique in many respects. The magnetosphere of Mercury is among the smallest in the solar system; its magnetic field typically stands off the solar wind only - 1000 to 2000 km above the surface. For this reason there are no closed drift paths for energetic particles and, hence, no radiation belts. The characteristic time scales for wave propagation and convective transport are short and kinetic and fluid modes may be coupled. Magnetic reconnection at the dayside magnetopause may erode the subsolar magnetosphere allowing solar wind ions to impact directly the regolith. Inductive currents in Mercury s interior may act to modify the solar wind interaction by resisting changes due to solar wind pressure variations. Indeed, observations of these induction effects may be an important source of information on the state of Mercury s interior. In addition, Mercury s magnetosphere is the only one with its defining magnetic flux tubes rooted in a planetary regolith as opposed to an atmosphere with a conductive ionospheric layer. This lack of an ionosphere is probably the underlying reason for the brevity of the very intense, but short-lived, - 1-2 min, substorm-like energetic particle events observed by Mariner 10 during its first traversal of Mercury s magnetic tail. Because of Mercury s proximity to the sun, 0.3 - 0.5 AU, this magnetosphere experiences the most extreme driving forces in the solar system. All of these factors are expected to produce complicated interactions involving the exchange and re-cycling of neutrals and ions between the solar wind, magnetosphere, and regolith. The electrodynamics of Mercury s magnetosphere are expected to be equally complex, with strong forcing by the solar wind, magnetic reconnection at the magnetopause and in the tail, and the pick-up of planetary ions all driving field-aligned electric currents. However, these field-aligned currents do not close in an ionosphere, but in some other manner. In addition to the insights- into magnetospheric physics offered by study of the solar wind - Mercury system, quantitative specification of the "external" magnetic field generated by magnetospheric currents is necessary for accurate determination of the strength and multi-polar decomposition of Mercury s intrinsic magnetic field. MESSENGER S highly capable instrumentation and broad orbital coverage will greatly advance our understanding of both the origin of Mercury s magnetic field and the acceleration of charged particles in small magnetospheres. In. this article, we review what is known about Mercury s magnetosphere and describe the MESSENGER science team s strategy for obtaining answers to the outstanding science questions surrounding the interaction of the solar wind with Mercury and its small, but dynamic, magnetosphere.

Slavin, James A.↗
DOE OSTI2016

Staggered scheduling of sensor estimation and fusion for tracking over long-haul links

Networked sensing can be found in a multitude of real-world applications. Here, we focus on the communication-and computation-constrained long-haul sensor networks, where sensors are remotely deployed over a vast geographical area to perform certain tasks. Of special interest is a class of such networks where sensors take measurements of one or more dynamic targets and send their state estimates to a remote fusion center via long-haul satellite links. The severe loss and delay over such links can easily reduce the amount of sensor data received by the fusion center, thereby limiting the potential information fusion gain and resulting in suboptimal tracking performance. In this paper, starting with the temporal-domain staggered estimation for an individual sensor, we explore the impact of the so-called intra-state prediction and retrodiction on estimation errors. We then investigate the effect of such estimation scheduling across different sensors on the spatial-domain fusion performance, where the sensing time epochs across sensors are scheduled in an asynchronous and staggered manner. In particular, the impact of communication delay and loss as well as sensor bias on such scheduling is explored by means of numerical and simulation studies that demonstrate the validity of our analysis.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Numerical Investigation of Wall-Cooling Effect on Aero-Optical Distortions for Hypersonic Boundary Layer

Compressible wall-modeled large-eddy simulations of Mach 8 turbulent boundary-layer flows over a flat plate were carried out for the conditions of the hypersonic wind tunnel at Sandia National Laboratories. The simulations provide new insight into the effect of wall cooling on the aero-optical path distortions for hypersonic turbulent boundary-layer flows. Four different wall-to-recovery temperature ratios, 0.3, 0.48, 0.71, and 0.89, are considered. Despite the much lower grid resolution, the mean velocity, temperature, and resolved Reynolds stress profiles from the simulation for a temperature ratio of 0.48 are in good agreement with those from a reference direct numerical simulation. The normalized root-mean-square optical path difference obtained from the present simulations is compared with that from reference direct numerical simulations, Sandia experiments, as well as predictions obtained with a semi-analytical model by Notre Dame University. Here the present analysis focuses on the effect of wall cooling on the wall-normal density correlations, on key underlying assumptions of the aforementioned model such as the strong Reynolds analogy, and on the elevation angle effect on the optical path difference. Wall cooling is found to increase the velocity fluctuations and decrease the density fluctuations, resulting in an overall reduction of the normalized optical path distortion. Compared to the simulations, the basic strong Reynolds analogy overpredicts the temperature fluctuations for cooled walls. Also different from the strong Reynolds analogy, the velocity and temperature fluctuations are not perfectly anticorrelated. Finally, as the wall temperature is raised, the density correlation length, away from the wall but inside the boundary layer, increases significantly for beam paths tilted in the downstream direction.

42 ENGINEERING↗
DOE OSTI2017

Effects of strain rate on mechanical properties and deformation behavior of an austenitic Fe-25Mn-3Al-3Si TWIP-TRIP steel

The effects of quasi-static and low-dynamic strain rate (ε = 10 –4 /s to ε = 10 2 /s) on tensile properties and deformation mechanisms were studied in a Fe-25Mn-3Al-3Si (wt%) twinning and transformation-induced plasticity [TWIP-TRIP] steel. The fully austenitic microstructure deforms primarily by dislocation glide but due to the room temperature stacking fault energy [SFE] of 21 ± 3 mJ/m 2 for this alloy, secondary deformation mechanisms such as mechanical twinning (TWIP) and epsilon martensite formation (TRIP) also play an important role in the deformation behavior. The mechanical twins and epsilon-martensite platelets act as planar obstacles to subsequent dislocation motion on non-coplanar glide planes and reduce the dislocation mean free path. A high-speed thermal camera was used to measure the increase in specimen temperature as a function of strain, which enabled the use of a thermodynamic model to predict the increase in SFE. The influence of strain rate and strain on microstructural parameters such as the thickness and spacing of mechanical twins and epsilon-martensite laths was quantified using dark field transmission electron microscopy, electron channeling contrast imaging, and electron backscattered diffraction. The effect of sheet thickness on mechanical properties was also investigated. Increasing the tensile specimen thickness increased the product of ultimate tensile strength and total elongation, but had no significant effect on uniform elongation or yield strength. The yield strength exhibited a significant increase with increasing strain rate, indicating that dislocation glide becomes more difficult with increasing strain rate due to thermally-activated short-range barriers. Finally, a modest increase in ultimate tensile strength and minimal decrease in uniform elongation were noted at higher strain rates, suggesting adiabatic heating, slight changes in strain-hardening rate and observed strain localizations as root causes, rather than a significant change in the underlying TWIP-TRIP mechanisms at low values of strain.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Aboveground biomass density models for NASA’s Global Ecosystem Dynamics Investigation (GEDI) lidar mission

NASA's Global Ecosystem Dynamics Investigation (GEDI) is collecting spaceborne full waveform lidar data with a primary science goal of producing accurate estimates of forest aboveground biomass density (AGBD). This paper presents the development of the models used to create GEDI's footprint-level (~25 m) AGBD (GEDI04_A) product, including a description of the datasets used and the procedure for final model selection. The data used to fit our models are from a compilation of globally distributed spatially and temporally coincident field and airborne lidar datasets, whereby we simulated GEDI-like waveforms from airborne lidar to build a calibration database. We used this database to expand the geographic extent of past waveform lidar studies, and divided the globe into four broad strata by Plant Functional Type (PFT) and six geographic regions. GEDI's waveform-to-biomass models take the form of parametric Ordinary Least Squares (OLS) models with simulated Relative Height (RH) metrics as predictor variables. From an exhaustive set of candidate models, we selected the best input predictor variables, and data transformations for each geographic stratum in the GEDI domain to produce a set of comprehensive predictive footprint-level models. We found that model selection frequently favored combinations of RH metrics at the 98th, 90th, 50th, and 10th height above ground-level percentiles (RH98, RH90, RH50, and RH10, respectively), but that inclusion of lower RH metrics (e.g. RH10) did not markedly improve model performance. Second, forced inclusion of RH98 in all models was important and did not degrade model performance, and the best performing models were parsimonious, typically having only 1-3 predictors. Third, stratification by geographic domain (PFT, geographic region) improved model performance in comparison to global models without stratification. Fourth, for the vast majority of strata, the best performing models were fit using square root transformation of field AGBD and/or height metrics. There was considerable variability in model performance across geographic strata, and areas with sparse training data and/or high AGBD values had the poorest performance. These models are used to produce global predictions of AGBD, but will be improved in the future as more and better training data become available.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2022

Aboveground biomass density models for NASA’s Global Ecosystem Dynamics Investigation (GEDI) lidar mission

NASA’s Global Ecosystem Dynamics Investigation (GEDI) is collecting spaceborne full waveform lidar data with a primary science goal of producing accurate estimates of forest aboveground biomass density (AGBD). This paper presents the development of the models used to create GEDI’s footprint-level (~25 m) AGBD (GEDI04_A) product, including a description of the datasets used and the procedure for final model selection. The data used to fit our models are from a compilation of globally distributed spatially and temporally coincident field and airborne lidar datasets, whereby we simulated GEDI-like waveforms from airborne lidar to build a calibration database. We used this database to expand the geographic extent of past waveform lidar studies, and divided the globe into four broad strata by Plant Functional Type (PFT) and six geographic regions. GEDI’s waveform-to-biomass models take the form of parametric Ordinary Least Squares (OLS) models with simulated Relative Height (RH) metrics as predictor variables. From an exhaustive set of candidate models, we selected the best input predictor variables, and data transformations for each geographic stratum in the GEDI domain to produce a set of comprehensive predictive footprint-level models. We found that model selection frequently favored combinations of RH metrics at the 98th, 90th, 50th, and 10th height above ground-level percentiles (RH98, RH90, RH50, and RH10, respectively), but that inclusion of lower RH metrics (e.g. RH10) did not markedly improve model performance. Second, forced inclusion of RH98 in all models was important and did not degrade model performance, and the best performing models were parsimonious, typically having only 1-3 predictors. Third, stratification by geographic domain (PFT, geographic region) improved model performance in comparison to global models without stratification. Fourth, for the vast majority of strata, the best performing models were fit using square root transformation of field AGBD and/or height metrics. There was considerable variability in model performance across geographic strata, and areas with sparse training data and/or high AGBD values had the poorest performance. These models are used to produce global predictions of AGBD, but will be improved in the future as more and better training data become available.

Laura Duncanson↗
NASA NTRS2002

Observer-Based Magnetic Bearing Controller Developed for Aerospace Flywheels

A prototype of a versatile, observer-based magnetic bearing controller for aerospace flywheels was successfully developed and demonstrated on a magnetic bearing test rig (see the photograph) and an actual flywheel module. The objective of this development included a fast, yet low risk, control development process, and a robust, high-performance controller for a large variety of flywheels. This required a good system model, an efficient development procedure, and a model-based controller that addressed the key problems associated with flywheel and bearing imbalance, sensor error, and vibration. The model used in this control development and tuning procedure included the flexible rotor dynamics and motor-induced vibrations. Such a model was essential for low-risk scheduling of speed-dependent control parameters and for reliable evaluation of novel control strategies. The successfully tested control prototype utilized an extended Kalman filter to estimate the true rotor principal-axis motion from the raw sensor position feedback. For control refinement, the extended Kalman filter also estimated and eliminated the combined effects of mass-imbalance and sensor runouts from the input data. A key advantage of the design based on the extended Kalman filter is its ability to accurately estimate both the rotor's principal-axis position and gyroscopic rates with the least amount of phase lag. This is important for control parameter scheduling to dampen the gyroscopic motions. Because of large uncertainties in the magnetic bearing and imbalance characteristics, this state-estimation scheme alone is insufficient for containing the rotor motion within the desired 1-mil excursion radius. A nonlinear gain adjustment based on an estimation of the principal-axis orbit size was needed to provide a coarse (nonoptimal), but robust, control of the orbit growth. Control current minimization was achieved with a (steepest gradient) search of synchronous errors in the principal-axis position input data. Actual flywheel tests of this observer-based controller (developed entirely in-house) at the NASA Glenn Research Center showed that the model correctly predicted the rotor orbit growth as a function of rotational speed, and it demonstrated the capability of gain adjustments to arrest this growth. Data from these tests on an actual flywheel module spun to 26,000 rpm proved that the controller was able to contain the shaft motion to within much less than 0.5 mils of radial excursion with axis currents less than 300 mA in root-mean-square estimate. The test speed range was limited because of thermal expansion concerns for this particular flywheel unit, not because of any deficiency in the controller. Simulations for this unit indicated that the controller should be robust up to its top operating speed of 60,000 rpm. Aside from these important achievements, and most significantly, it took less than 1 week to adapt this controller from the simple test rig to the actual flywheel and to demonstrate full five-axis levitation and control. This demonstration showed that both the controller and the model-based development and tuning framework are easily adaptable to a wide range of rotors and bearing configurations and, hence, are capable of reducing design risks and costs for many future flywheel technology developments.

Le, Dzu K.↗
DOE OSTI2016

Combined measurement and modeling of the hydrological impact of hydraulic redistribution using CLM4.5 at eight AmeriFlux sites

Effects of hydraulic redistribution (HR) on hydrological, biogeochemical, and ecological processes have been demonstrated in the field, but the current generation of standard earth system models does not include a representation of HR. Though recent studies have examined the effect of incorporating HR into land surface models, few (if any) have done cross-site comparisons for contrasting climate regimes and multiple vegetation types via the integration of measurement and modeling. Here, we incorporated the HR scheme of Ryel et al. (2002) into the NCAR Community Land Model Version 4.5 (CLM4.5), and examined the ability of the resulting hybrid model to capture the magnitude of HR flux and/or soil moisture dynamics from which HR can be directly inferred, to assess the impact of HR on land surface water and energy budgets, and to explore how the impact may depend on climate regimes and vegetation conditions. Eight AmeriFlux sites with contrasting climate regimes and multiple vegetation types were studied, including the Wind River Crane site in Washington State, the Santa Rita Mesquite savanna site in southern Arizona, and six sites along the Southern California Climate Gradient. HR flux, evapotranspiration (ET), and soil moisture were properly simulated in the present study, even in the face of various uncertainties. Our cross-ecosystem comparison showed that the timing, magnitude, and direction (upward or downward) of HR vary across ecosystems, and incorporation of HR into CLM4.5 improved the model-measurement matches of evapotranspiration, Bowen ratio, and soil moisture particularly during dry seasons. Lastly, our results also reveal that HR has important hydrological impact in ecosystems that have a pronounced dry season but are not overall so dry that sparse vegetation and very low soil moisture limit HR.

58 GEOSCIENCES↗
DOE OSTI2021

Vegetation Warming Experiment: 15N Uptake Experiment Arctagrostis latifolia Biomass and Chemistry, Utqiagvik (Barrow), Alaska, 2018

This dataset consists of measured tissue specific biomass, %C, % N and APE (atom percentage in excess) of 15N of harvested A. latifolia in vegetation warming experiment enclosures and paired control plots located on the BEO, Utqiagvik, Alaska. Vegetation warming chambers (Zero Power) were deployed on the Barrow Environmental Observatory (BEO), Utqiagvik, Alaska. These chambers (Figure 1) consistently elevated air temperatures by approximately 4°C using a self-venting system described by Lewin et al (2017). Five chambers were deployed from June 17, 2018 to September 24, 2018 on the BEO within a 1 km2 area centered on 71.275N, -156.641W. Each chamber was co-located with an ambient plot where temperatures were not manipulated on patches of tundra containing the target species Arctagrostis latifolia. An intensive field campaign in late July investigated the impact of warming had on A. latifolia biomass, chemistry, and uptake of 15N labeled ammonia that was injected into the surface soils for one week. Initial measurements were taken on July 21, 2018. Harvest occurred on July 27, 2018. Water-extractable nutrients in soils were measured in July following harvests of A. latifolia plants and underlying soils. Availability of ammonia, nitrate, and phosphate throughout the growing season was measured by extracting nutrients bound to anion and cation binding resins deployed from July through September. Environmental variables (thaw depth, surface soil temperatures, surface soil moisture) were measured. Leaf traits and root traits of A. latifolia were also measured. This dataset includes one csv file and one pdf file. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

Stormtime ring current and radiation belt ion transport: Simulations and interpretations

We use a dynamical guiding-center model to investigate the stormtime transport of ring current and radiation-belt ions. We trace the motion of representative ions' guiding centers in response to model substorm-associated impulses in the convection electric field for a range of ion energies. Our simple magnetospheric model allows us to compare our numerical results quantitatively with analytical descriptions of particle transport, (e.g., with the quasilinear theory of radial diffusion). We find that 10-145-keV ions gain access to L approximately 3, where they can form the stormtime ring current, mainly from outside the (trapping) region in which particles execute closed drift paths. Conversely, the transport of higher-energy ions (approximately greater than 145 keV at L approximately 3) turns out to resemble radial diffusion. The quasilinear diffusion coefficient calculated for our model storm does not vary smoothly with particle energy, since our impulses occur at specific (although randomly determined) times. Despite the spectral irregularity, quasilinear theory provides a surprisingly accurate description of the transport process for approximately greater than 145-keV ions, even for the case of an individual storm. For 4 different realizations of our model storm, the geometric mean discrepancies between diffusion coefficients D(sup sim, sub LL) obtained from the simulations and the quasilinear diffusion coefficient D(sup ql, sub LL) amount to factors of 2.3, 2.3, 1.5, and 3.0, respectively. We have found that these discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) can be reduced slightly by invoking drift-resonance broadening to smooth out the sharp minima and maxima in D(sup ql, sub LL). The mean of the remaining discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) for the 4 different storms then amount to factors of 1.9, 2.1, 1.5, and 2.7, respectively. We find even better agreement when we reduce the impulse amplitudes systematically in a given model storm (e.g., reduction of all the impulse amplitudes by half reduces the discrepancy factor by at least its square root) and also when we average our results over an ensemble of 20 model storms (agreement is within a factor of 1.2 without impulse-amplitude reduction). We use our simulation results also to map phase-space densities f in accordance with Liouville's theorem. We find that the stormtime transport of approximately greater than 145-keV ions produces little change in f-bar the drift-averaged phase-space density on any drift shell of interest. However, the stormtime transport produces a major enhancement from the pre-storm phase-space density at energies approximately 30-145 keV, which are representative of the stormtime ring current.

Lyons, Larry R.↗
NASA NTRS1999

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗
DOE OSTI2025

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Duke Forest FACE (FACTS-I): Plant and Soil Response Data

This dataset describing the responses of plant and soil pools and fluxes to elevated atmospheric CO2 concentration and increased nitrogen supply was collected from Duke Forest Free Air CO2 Enrichment (FACE) – Forest-Atmosphere Carbon Transfer and Storage (FACTS-I) experiment from 1996 to 2012. The dataset includes data files for allometry (diameter at breast height, tree height, and height to live crown base), leaf area index, biomass (stem, branch, foliage, and root biomass, tree density, and basal area), net primary productivity (stem, branch, foliage, reproductive, and coarse root NPP), sap flux density, soil CO2 efflux, and stem temperature. Data files were formatted as .csv (Microsoft Excel or other spreadsheet programs can be used to read the format) and file descriptions, including variable name, unit, and data range, can be found in ‘FileDescription_[data_name].txt’ files. The Duke FACE experiment was in a loblolly pine (Pinus taeda L.) plantation established in 1983. Naturally regenerated broadleaved species including sweetgum (Liquidambar styraciflua L.) and tulip poplar (Liriodendron tulipifera L.), mostly in the overstory, and winged elm (Ulmus alata Michx.) and red maple (Acer rubrum L.) were common in the understory. The FACE experiment commenced with two plots (plots 7-8) in 1994 (Oren et al. 2001), with six additional plots (plots 1-6) coming online on 27 August 1996. CO2 enrichment was terminated on 31 October 2010 and post-enrichment data collection continued through 2012. Complete fertilization was applied annually to half of plots 7-8 from 1998 to 2004. The nutrient addition experiment expanded to half of plots 1-6 with a common protocol of N-fertilization in 2005 and continued until 2012. The levels of treatment in this dataset were expressed as ambient CO2 (AMB) or elevated CO2 (ELE) for CO2 treatment and control soil (CONT) or fertilized soil (FERT) for N treatment, respectively.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1998

Diurnal Reflectance Changes in Vegetation Observed with AVIRIS

Among the most important short-term dynamic biological processes are diurnal changes in canopy water relations. Plant regulation of water transport through stomatal openings affects other gaseous transport processes, often dramatically decreasing photosynthetic fixation of carbon dioxide during periods of water stress. Water stress reduces stomatal conductance of water vapor through the leaf surface and alters the diurnal timing of stomatal opening. Under non-water stressed conditions, stomates typically open soon after dawn and transpire water vapor throughout the daylight period. During stress periods, stomates may close for part of the day, generally near mid-day. Under prolonged stress conditions, stomatal closure shifts to earlier times during the day; stomates may close by mid-morning and remain closed until the following morning - or remain closed entirely. Under these conditions the relationship between canopy greenness (e.g., measured with a vegetation index or by spectral mixture analysis) and photosynthetic fixation of carbon is lost and the remotely sensed vegetation metric is a poor predictor of gas exchange. Prediction of stomatal regulation and exchange of water and trace gases is critical for ecosystem and climate models to correctly estimate budgets of these gases and understand or predict other processes like gross and net ecosystem primary production. Plant gas exchange has been extensively studied by physiologists at the leaf and whole plant level and by biometeorologists at somewhat larger scales. While these energy driven processes follow a predictable if somewhat asymmetric diurnal cycle dependent on soil water availability and the constraints imposed by the solar energy budget, they are nonetheless difficult to measure at the tree and stand levels using conventional methods. Ecologists have long been interested in the potential of remote sensing for monitoring physiological changes using multi-temporal images. Much of this research has focused on day-to-day changes in water use, especially for agricultural applications. Ustin et al. showed seasonal changes in canopy water content in chaparral shrub could be estimated using optical methods. Vanderbilt et al. followed asymmetric diurnal changes in the reflectance of a walnut orchard, but could not attribute specific reflectance changes to specific changes in canopy architecture or physiology. Forests and shrub lands in California experience prolonged periods of drought, sometimes extending six months without precipitation. The conifer and evergreen chaparral communities common to the foothill region around the central valley of California retain their foliage throughout the summer and have low transpiration rates despite high net radiation and temperature conditions. In contrast, grasslands and drought resistant deciduous species in the same habitat are seasonally dormant in summer. Because of differences in the mechanisms of drought tolerance, rooting depth and physiology between different plant communities in the region, it is likely that they display differences in diurnal water relations. The presence of diverse plant communities provides an opportunity to investigate possible diurnal landscape patterns in water relations that could be observed by an airborne hyperspectral scanner. This investigation of AVIRIS data collected over forest and shrub land represents the continuation of a prior investigation involving spectral mixture analysis of diurnal effects in the same AVIRIS data set.

Vanderbilt, V. C.↗