SEARCH · Engineering Papers
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.
To investigate fully balanced steady-state free precession (bSSFP) with optimized acquisition protocols for magnetic resonance imaging (MRI)-based postimplant quality assessment of low-dose-rate (LDR) prostate brachytherapy without an endorectal coil (ERC).
Here we consider the problem of radiation transport through purely absorbing particle clouds. The gold-standard solution of particle-resolved Monte Carlo ray-tracing method to this problem is computationally expensive and therefore solving the radiation transport equation (specifically, the Beer-Bouguer law (BB-law)) on a Eulerian mesh is often preferred. While the absorption coefficient in the real problem is infinite, the BB-law approximates it to be a finite number in the form of number density through a set of assumptions. For particle clouds that do not obey these assumptions, the BB-law predicts an incorrect exponential decay. Also, when the number density is computed using the nearest-neighbor approach, the BB-law solution diverges when the Eulerian mesh size becomes closer to or smaller than the particle size. This numerical divergence is due to the homogenization error. Although the filtering strategy for number density minimizes the homogenization error, it still converges to an incorrect exponential decay for particle clouds that break the BB-law constraints and the cost of the filtering is equivalent to that of the gold-standard solution. In this study, we develop a novel, highly accurate, verifiable and cost-effective solution to the radiation transport equation on a Eulerian domain using an energy balance approach where we derive an expression for the absorption coefficient as a function of particle and Eulerian mesh sizes. We apply our new method to Poisson and turbulent particle clouds that violate all the BB-law constraints and show that the solution is converging upon the mesh refinement, and eventually, we recover the same gold-standard solution for a much cheaper computational cost.
An extended population balance model (PBM) and a deep learning-based enhanced deep neural operator (DNO+) model are introduced for predicting particle size distribution (PSD) of comminuted biomass through a large knife mill. Experimental tests using corn stalks with varied moisture contents, mill blade speeds, and discharge screen sizes are conducted to support model development. A novel mechanism in the extended PBM allows for including additional input parameters such as moisture content, which is not possible in the original PBM. The DNO+ model can include influencing factors of different data types such as moisture content and discharge screen size, which significantly extends the engineering applicability of the standard DNO model that only admits feed PSD and outcome PSD. Test results show that both models are remarkably accurate in the calibration or training parameter space and can be used as surrogate models to provide effective guidance for biomass preprocessing design.
CD4{sup +}CD25{sup +} regulatory T (Treg) cells and Th17 cells play important roles in the progression of metabolic-associated fatty liver disease (MAFLD). However, the contribution of monokine induced by interferon-gamma (MIG)/CXCL9 to the Treg/Th17 imbalance in MAFLD is only partially understood. In the present study, we detected increased levels of MIG/CXCL9 and a Treg/Th17 imbalance in the setting of metabolic-associated steatohepatitis (MASH). Recombinant adeno-associated virus-mediated gene transfer and silencing of MIG/CXCL9 expression in mice alleviated MASH and increased the Treg/Th17 ratio. Furthermore, the percentage of Th17 cells, but not Treg cells, differentiated from splenic CD4{sup +} T cells was significantly increased by administration of MIG/CXCL9. MIG/CXCL9 also promoted Th17 cell proliferation, and its effects were dose dependent. Levels of phosphorylated c-Jun N-terminal kinase (JNK) decreased dramatically when MIG/CXCL9 was inhibited in a murine MASH model. In cultured Treg cells, phosphorylated JNK levels decreased dose-dependently in response to MIG/CXCL9 inhibition, but increased in cultured Th17 cells. This effect was blocked in the presence of a JNK inhibitor. These findings underline the fundamental importance of MIG/CXCL9 in maintaining the Treg/Th17 balance in MAFLD and provide the foundations for a novel approach to preventing and treating MAFLD.
The attached eddy hypothesis of Townsend (The Structure of Turbulent Shear Flow, 1956, Cambridge University Press) states that the logarithmic mean velocity admits self-similar energy-containing eddies which scale with the distance from the wall. Over the past decade, there has been a significant amount of evidence supporting the hypothesis, placing it to be the central platform for the statistical description of the general organisation of coherent structures in wall-bounded turbulent shear flows. Nevertheless, the most fundamental question, namely why the hypothesis has to be true, has remained unanswered over many decades. Under the assumption that the integral length scale is proportional to the distance from the wall y , in the present study we analytically demonstrate that the mean velocity is a logarithmic function of y if and only if the energy balance at the integral length scale is self-similar with respect to y , providing a theoretical basis for the attached eddy hypothesis. The analysis is subsequently verified with the data from a direct numerical simulation of incompressible channel flow at the friction Reynolds number Reτ≃5200 (Lee & Moser, J. Fluid Mech., vol. 774, 2015, pp. 395–415).
Here, the ongoing accumulation of plastic waste in landfills and in the environment is driving research on chemical processes and catalysts to recycle polymers. Traditional modeling strategies are not applicable to these processes because they involve too many reactants and intermediates, one for each molecular weight and each functionalization. To model the kinetics, we have developed population balance models (PBMs) that account for macromolecular reactants in the bulk and macromolecular catalytic intermediates. These PBMs couple to each other through polymer adsorption and desorption models and to traditional rate equations for small molecule products and co-reactants (like hydrogen or ethylene). The models, in combination with experimental data, are being used in many ways: (i) to test mechanistic hypotheses, (ii) to extract rate parameters, (iii) to quantitatively compare catalyst activities, (iv) to account for mass transfer and vapor–liquid partitioning in two-phase reactors, and (v) to design novel support architectures and catalysts that mimic the processive action of natural depolymerization enzymes. Some key theoretical advances allow PBMs to be constructed from elementary rates and mechanisms, as opposed to traditional formulations with pseudoelementary rate parameters invoked as fitting parameters. We discuss ways to build these models “bottom-up” from first-principles calculations and ways to extract model parameters from “top down” analyses of rate data. The combination provides a quantitative bridge between first-principles calculations and the kinetics of complex macromolecular transformations for polymer upcycling and beyond.
Magnetic circular dichroism of transition metal complexes and open-shell systems are challenging to simulate and analyze, mainly due to the interplay of spin–orbit couplings and finite-magnetic-field induced Zeeman effects with the complex selection rules dictated by the circularly polarized light. In this study, we introduce an ab initio relativistic two-component formalism based on the restricted magnetic-balanced Hamiltonian for simulating MCD spectra. Both homogeneous finite magnetic field and relativistic effects are included variationally in the ground state reference. Finite-field London orbitals are used to enforce the constrained gauge-origin independence in the calculation using localized atomic orbitals. Through benchmark studies of AuCl 4 – , Pt(CN) 4 2– , and Mo(CN) 8 3– , we discuss how relativistic effects are manifested in MCD for both closed-shell and open-shell molecular complexes and how the interplay between spin–orbit coupling and magnetic field modulates the MCD selection rules. Finally, an investigation on temperature-dependent MCD is carried out and compared to experiment.
A panchromatic triad and a charge-separation triad are joined in a crossbar architecture to capture solar energy. The panchromatic triad (T) is comprised of a central free-base porphyrin that is strongly coupled via direct ethyne linkages to two perylene-monoimides. Here, the charge-separation unit incorporates a free-base or zinc chlorin (C or ZnC) as a hole acceptor (or electron donor) and a perylene-diimide (PDI) as an electron acceptor, both attached to the porphyrin via diphenylethyne linkers. The free-base porphyrin is common to both light-harvesting and charge-separation motifs. The chlorin and PDI also function as ancillary light absorbers, complementing direct excitation of the panchromatic triad to produce the discrete lowest excited state of the array (T*). Attainment of full charge separation across the pentad entails two steps: (1) an initial excited-state hole/electron transfer process to oxidize the chlorin (and reduce the triad) or reduce the PDI (and oxidize the triad); and (2) subsequent ground-state electron/hole migration to produce oxidized chlorin and reduced PDI. Full charge separation for pentad ZnC-T-PDI to generate ZnC + -T-PDI¯ occurs with a quantum yield of ~30% and mean lifetime ~1 μs in dimethyl sulfoxide. For C-T-PDI, initial charge separation is followed by rapid charge recombination. The molecular designs and studies reported here reveal the challenges of balancing the demands for charge separation (linker length and composition, excited-state energies, redox potentials, and medium polarity) with the constraints for panchromatic absorption (strong electronic coupling of the porphyrin and two perylene-monoimide units) for integrated function in solar-energy conversion.
Abstract The surface energy imbalance, viz. that the turbulent dissipation does not fully account for available energy, has for long been an outstanding challenge in geophysical studies. In this study, we developed a novel approach based on nonequilibrium thermodynamics by representing the atmospheric boundary layer as a heat engine. In addition, an analytically tractable approach was used to estimate the ground heat flux based on Green's function approach, which in turn determines the available energy that drives the atmospheric heat engine. The proposed model was evaluated using heat fluxes measured by eight AmeriFlux eddy covariance towers with atmospheric temperature profiles recorded at adjacent radiosonde sites. The surface energy balance closure can be improved by ~11% over various landscapes, by including the estimated power production from the atmospheric heat engine.
Abstract The Hamburg Aerosol Module version 2.3 (HAM2.3) from the ECHAM6.3‐HAM2.3 global atmosphere‐aerosol model is coupled to the recently developed icosahedral nonhydrostatic ICON‐A (icon‐aes‐1.3.00) global atmosphere model to yield the new ICON‐A‐HAM2.3 atmosphere‐aerosol model. The ICON‐A and ECHAM6.3 host models use different dynamical cores, parameterizations of vertical mixing due to sub‐grid scale turbulence, and parameter settings for radiation balance tuning. Here, we study the role of the different host models for simulated aerosol optical thickness (AOT) and evaluate impacts of using HAM2.3 and the ECHAM6‐HAM2.3 two‐moment cloud microphysics scheme on several meteorological variables. Sensitivity runs show that a positive AOT bias over the subtropical oceans is remedied in ICON‐A‐HAM2.3 because of a different default setting of a parameter in the moist convection parameterization of the host models. The global mean AOT is biased low compared to MODIS satellite instrument retrievals in ICON‐A‐HAM2.3 and ECHAM6.3‐HAM2.3, but the bias is larger in ICON‐A‐HAM2.3 because negative AOT biases over the Amazon, the African rain forest, and the northern Indian Ocean are no longer compensated by high biases over the sub‐tropical oceans. ICON‐A‐HAM2.3 shows a moderate improvement with respect to AOT observations at AERONET sites. A multivariable bias score combining biases of several meteorological variables into a single number is larger in ICON‐A‐HAM2.3 compared to standard ICON‐A and standard ECHAM6.3. In the tropics, this multivariable bias is of similar magnitude in ICON‐A‐HAM2.3 and in ECHAM6.3‐HAM2.3. In the extra‐tropics, a smaller multivariable bias is found for ICON‐A‐HAM2.3 than for ECHAM6.3‐HAM2.3.
Abstract Quantifying flow and transport from hillslopes is vital for understanding water quantity and quality in rivers, but remains obscure because of limited subsurface measurements. Using measured hydraulic conductivity K profiles and water balance over a single year to calibrate a transmissivity feedback model for a hillslope in the East River watershed (Colorado) proved unsatisfactory for predicting flow over the subsequent years. Well‐constrained field‐scale K were obtained by optimizing subsurface flux predictions over years having large differences in recharge, and by including estimates of interannual transfer of excess snowmelt recharge. Water and solute exports during high snowmelt recharge occur predominantly via shallow groundwater flow through weathered rock and soil because of their enlarged transmissivities under saturated conditions. Conversely, these shallow pathways are less active in snow drought years when the water table remains deeper within the weathering zone. Hillslope soil water monitoring showed that rainfall does not infiltrate deeply during summer and fall months, and revealed water losses consistent with model ET predictions. By combining water table‐dependent fluxes with pore water chemistry in different zones, time‐dependent rates of solute exports become predictable. As an example, calibrated K were combined with dissolved nitrogen concentrations in pore waters to show the snowmelt‐dependence of reactive nitrogen exported from the hillslope, further supporting the recent finding that the weathering zone is the dominant source of reactive nitrogen at this site. Subsurface export predictions can now be obtained for wide ranges of recharge based on measurements of water table elevation and profiles of pore water chemistry.
Plant responses to water stress is a major uncertainty to predicting terrestrial ecosystem sensitivity to drought. Different approaches have been developed to represent plant water stress. Empirical approaches (the empirical soil water stress (or Beta) function and the supply-demand balance scheme) have been widely used for many decades; more mechanistic based approaches, that is, plant hydraulic models (PHMs), were increasingly adopted in the past decade. However, the relationships between them—and their underlying connections to physical processes—are not sufficiently understood. This limited understanding hinders informed decisions on the necessary complexities needed for different applications, with empirical approaches being mechanistically insufficient, and PHMs often being too complex to constrain. Here we introduce a unified framework for modeling transpiration responses to water stress, within which we demonstrate that empirical approaches are special cases of the full PHM, when the plant hydraulic parameters satisfy certain conditions. We further evaluate their response differences and identify the associated physical processes. Finally, we propose a methodology for assessing the necessity of added complexities of the PHM under various climatic conditions and ecosystem types, with case studies in three typical ecosystems: a humid Midwestern cropland, a semi-arid evergreen needleleaf forest, and an arid grassland. Notably, Beta function overestimates transpiration when VPD is high due to its lack of constraints from hydraulic transport and is therefore insufficient in high VPD environments. With the unified framework, we envision researchers can better understand the mechanistic bases of and the relationships between different approaches and make more informed choices.
Large spatio-temporal gradients in the Congo basin vegetation and rainfall are observed. However, its water-balance (evapotranspiration minus precipitation, or ET - P) is typically measured at basin-scales, limited primarily by river-discharge data, spatial resolution of terrestrial water storage measurements, and poorly constrained ET. We use observations of the isotopic composition of water vapor to quantify the spatio-temporal variability of net surface water fluxes across the Congo Basin between 2003 and 2018. These data are calibrated at basin scale using satellite gravity and total Congo river discharge measurements and then used to estimate time-varying ET - P over four quadrants representing the Congo Basin, providing first estimates of this kind for the region. We find that the multi-year record, seasonality, and interannual variability of ET - P from both the isotopes and the gravity/river discharge based estimates are consistent. Additionally, we use precipitation and gravity-based estimates with our water vapor isotope-based ET - P to calculate time and space averaged ET and net river discharge within the Congo Basin. These quadrant-scale moisture flux estimates indicate (a) substantial recycling of moisture in the Congo Basin (temporally and spatially averaged ET/P > 70%), consistent with models and visible light-based ET estimates, and (b) net river outflow is largest in the Western Congo where there are more rivers and higher flow rates. Our results confirm the importance of ET in modulating the Congo water cycle relative to other water sources.
During the 2015–2016 El Niño, the Amazon basin released almost one gigaton of carbon (GtC) into the atmosphere due to extreme temperatures and drought. The link between the drought impact and recovery of the total carbon pools and its biogeochemical drivers is still unknown. With satellite-constrained net carbon exchange and its component fluxes including gross primary production and fire emissions, we show that the total carbon loss caused by the 2015–2016 El Niño had not recovered by the end of 2018. Forest ecosystems over the Northeastern (NE) Amazon suffered a cumulative total carbon loss of ~0.6 GtC through December 2018, driven primarily by a suppression of photosynthesis whereas southeastern savannah carbon loss was driven in part by fire. We attribute the slow recovery to the unexpected large carbon loss caused by the severe atmospheric aridity coupled with a water storage deficit during drought. We show the attenuation of carbon uptake is three times higher than expected from the pre-drought sensitivity to atmospheric aridity and ground water supply. Our study fills an important knowledge gap in our understanding of the unexpectedly enhanced response of carbon fluxes to atmospheric aridity and water storage deficit and its impact on regional post-drought recovery as a function of the vegetation types and climate perturbations. Our results suggest that the disproportionate impact of water supply and demand could compromise resiliency of the Amazonian carbon balance to future increases in extreme events.
Abstract Forecasts of root growth and carbon sequestration under global change are compromised by uncertainty in how plants will allocate biomass between above and belowground pools. Here, we develop a simple model to assess whether functional balance theory can explain a complex biomass allocation response observed in a brackish marsh under experimental warming and elevated CO 2 . Our model shows how treatment‐driven changes in nitrogen supply and demand can explain divergent observations of root growth (i.e., maximum responses under intermediate warming and elevated CO 2 ). The model also reveals a surprising interaction between warming and eutrophication, where enhanced N loading to coastal marshes may reduce adverse impacts of warming on root growth. Our findings provide a mechanistic basis for incorporating biomass allocation into forecast models of marsh evolution. They also provide a general example of using ecological theory to decompose complex net responses observed in multi‐factor global change experiments into constituent processes.
Abstract Land surface models (LSMs) play a crucial role in elucidating water and carbon cycles by simulating processes such as plant transpiration and evaporation from bare soil, yet calibration often relies on comparing LSM outputs of landscape total evapotranspiration ( ET ) and discharge with measured bulk fluxes. Discrepancies in partitioning into component fluxes predicted by various LSMs have been noted, prompting the need for improved evaluation methods. Stable water isotopes serve as effective tracers of component hydrologic fluxes, but data and model integration challenges have hindered their widespread application. Leveraging National Ecological Observation Network measurements of water isotope ratios at 16 US sites over 3 years combined with LSM‐modeled fluxes, we employed an isotope‐enabled mass balance framework to simulate ET isotope values ( δET ) within three operational LSMs (Mosaic, Noah, and VIC) to evaluate their partitioning. Models simulating δET values consistent with observations were deemed more reflective of water cycling in these ecosystems. Mosaic exhibited the best overall performance (Kling‐Gupta Efficiency of 0.28). For both Mosaic and Noah there were robust correlations between bare soil evaporation fraction and error (negative) as well as transpiration fraction and error (positive). We found the point at which errors are smallest ( x ‐intercept of the multi‐site regression) is at a higher transpiration fraction than is currently specified in the models. Which means that transpiration fraction is underestimated on average. Stable isotope tracers offer an additional tool for model evaluation and identifying areas for improvement, potentially enhancing LSM simulations and our understanding of land‐surface hydrologic processes.
Integrating genome-scale metabolic networks with reactive transport models (RTMs) provides a detailed description of the dynamic changes in microbial growth and metabolism. Despite promising demonstrations in the past, computational inefficiency has been pointed out as a critical issue to overcome because it requires repeated application of linear programming (LP) to obtain flux balance analysis (FBA) solutions in every time step and spatial grid. To address this challenge, we propose a new simulation method where we train and validate artificial neural networks (ANNs) using randomly sampled FBA solutions and incorporate the resulting surrogate FBA model (represented as algebraic equations) into RTMs as source/sink terms. We demonstrate the efficiency of our method via a case study of Shewanella oneidensis MR-1. During aerobic growth on lactate, S. oneidensis produces metabolic byproducts (such as pyruvate and acetate), which are subsequently consumed as alternative carbon sources when the preferred nutrients are depleted. To effectively simulate these complex dynamics, we used a cybernetic approach that models metabolic switches as the outcome of dynamic competition among multiple growth options. In both zero-dimensional batch and one-dimensional column configurations, the ANN-based surrogate models achieved substantial reduction of computational time by several orders of magnitude compared to the original LP-based FBA models. Moreover, the ANN models produced robust solutions without any special measures to prevent numerical instability. These developments significantly promote our ability to utilize genome-scale networks in complex, multi-physics, and multi-dimensional ecosystem modeling.
How coastal forest productivity varies with local nutrient availability and water supply remains a knowledge gap under climate change. In a two-decade field experiment manipulating fertilization and density in coastal pine forests, we show that a decade of growth enhancement by simulated sedimentary nutrient inputs has resulted in a striking reversal in growth and increased mortality risk as drought and sea level rise progressed. Recent groundwater salinization has further triggered a shift from nutrient to water limitation, causing severe stomatal closure and decoupling of tree carbon–water balance, which induces a negative intrinsic water use efficiency (iWUE)–growth relationship. Higher tree iWUE predicted sharper tree growth declines, and both nutrient enrichment and high stand density amplified this feedback, increasing the risk of hydraulic failure and mortality. These results suggest that a transient nutrient-stimulatory effect could drive further coastal forest degradation due to heightened belowground saltwater stress under sea-level rise.