Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “partial correlation”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evolution of Charge-Lattice Dynamics across the Kuramoto Synchronization Phase Diagram of Quantum Tunneling Polarons in Cuprate Superconductors

Because of its sensitivity to the instantaneous structure factor, S(Q,t = 0), Extended X-ray Absorption Fine Structure (EXAFS) is a powerful tool for probing the dynamic structure of condensed matter systems in which the charge and lattice dynamics are coupled. When applied to hole-doped cuprate superconductors, EXAFS has revealed the presence of internal quantum tunneling polarons (IQTPs). An IQTP arises in EXAFS as a two-site distribution for certain Cu–O pairs, which is also duplicated in inelastic scattering but not observed in standard diffraction measurements. The Cu–Sr pair distribution has been found to be highly anharmonic and strongly correlated to both the IQTPs and to superconductivity, as, for example, in YSr 2 Cu 2.75 Mo 0.25 O 7.54 (T c =84 K). In order to describe such nontrivial, anharmonic charge-lattice dynamics, we have proposed a model Hamiltonian for a prototype six-atom cluster, in which two Cu-apical-O IQTPs are charge-transfer bridged through Cu atoms by an O atom in the CuO 2 plane and are anharmonically coupled via a Sr atom. By applying an exact diagonalization procedure to this cluster, we have verified that our model indeed produces an intricate interplay between charge and lattice dynamics. Then, by using the Kuramoto model for the synchronization of coupled quantum oscillators, we have found a first-order phase transition for the IQTPs into a synchronized, phase-locked phase. Most importantly, we have shown that this transition results specifically from the anharmonicity. Finally, we have provided a phase diagram showing the onset of the phase-locking of IQTPs as a function of the charge-lattice and anharmonic couplings in our model. We have found that the charge, initially confined to the apical oxygens, is partially pumped into the CuO 2 plane in the synchronized phase, which suggests a possible connection between the synchronized dynamic structure and high-temperature superconductivity (HTSC) in doped cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

X-ray Photoelectron Spectroscopy Analysis of Nafion-Containing Samples: Pitfalls, Protocols, and Perceptions of Physicochemical Properties

X-ray photoelectron spectroscopy (XPS) is one of the most common techniques used to analyze the surface composition of catalysts and support materials used in polymer electrolyte membrane (PEM) fuel cells and electrolyzers, providing important insights for further improvement of their properties. Characterization of catalyst layers (CLs) is more challenging, which can be at least partially attributed to the instability of ionomer materials such as Nafion during measurements. This work explores the stability of Nafion during XPS measurements, illuminating and addressing Nafion degradation concerns. The extent of Nafion damage as a function of XPS instrumentation, measurement conditions, and sample properties was evaluated across multiple instruments. Results revealed that significant Nafion damage to the ion-conducting sulfonic acid species (>50% loss in sulfur signal) may occur in a relatively short time frame (tens of minutes) depending on the exact nature of the sample and XPS instrument. This motivated the development and validation of a multipoint XPS data acquisition protocol that minimizes Nafion damage, resulting in reliable data acquisition by avoiding significant artifacts from Nafion instability. The developed protocol was then used to analyze both thin film ionomer samples and Pt/C-based CLs. Comparison of PEM fuel cell CLs to Nafion thin films revealed several changes in Nafion spectral features attributed to charge transfer due to interaction with conductive catalyst and support species. This study provides a method to reliably characterize ionomer-containing samples, facilitating fundamental studies of the catalyst-ionomer interface and more applied investigations of structure-processing-performance correlations in PEM fuel cell and electrolyzer CLs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Participation in the TOMS Science Team

Because of the nominal funding provided by this grant, some of the relevant research is partially funded by other sources. Research performed for this funding period included the following items: We have investigated errors in TOMS ozone measurements caused by the uncertainty in wavelength calibration, coupled with the ozone cross sections in the Huggins bands and their temperature dependence. Preliminary results show that 0.1 nm uncertainty in TOMS wavelength calibration at the ozone active wavelengths corresponds to approx. 1% systematic error in O3, and thus potential 1% biases among ozone trends from the various TOMS instruments. This conclusion will be revised for absolute O3 Measurements as cross sections are further investigated for inclusion in the HITRAN database at the SAO, but the potential for relative errors remains. In order to aid further comparisons among TOMS and GOME ozone measurements, we have implemented our method of direct fitting of GOME radiances (BOAS) for O3, and now obtain the best fitting precision to date for GOME O3 Columns. This will aid in future comparisons of the actual quantities measured and fitted for the two instrument types. We have made comparisons between GOME ICFA cloud fraction and cloud fraction determined from GOME data using the Ring effect in the Ca II lines. There is a strong correlation, as expected, but there are substantial systematic biases between the determinations. This study will be refined in the near future using the recently-developed GOME Cloud Retrieval Algorithm (GOMECAT). We have improved the SAO Ring effect determination to include better convolution with instrument transfer functions and inclusion of interferences by atmospheric absorbers (e.g., O3). This has been made available to the general community.

Chance, Kelly↗
DOE OSTI2019

Turbulent flow and vortex characteristics in a blocked subchannel of a helically wrapped rod bundle

The study of flow and heat transfer through a bundle of rods is important, as this type of configuration is widely used in many engineering applications, for example, in heat exchangers, clusters of support structures, steam generators, and nuclear power reactors. This work experimentally investigated the flow-field characteristics of rod bundles with helically wrapped wires (wire-wrapped rod assemblies) proposed for advanced burner reactors. Subchannel blockage is one of the most interesting accident scenarios for liquid metal fast reactors as different mechanisms can lead to a partial or total subchannel blockage. It is of paramount important to understand the thermal-hydraulic behavior of the flow in the vicinity of such blockages. The flow mixing characteristics in two interior subchannels located near the core of the bundle, one without and one with the presence of a blockage, were experimentally investigated. The velocity fields within the spatial gaps of two subchannels were obtained by applying the matched-index-of-refraction and time-resolved particle image velocimetry (TR-PIV) technique for Reynolds numbers of 4000 and 17,000. The first- and second-order statistics, including the mean velocity, root-mean-square fluctuating velocity, and Reynolds stress, were computed from the obtained TR-PIV velocity vector fields. In the blocked sub-channel, the TR-PIV results revealed a large recirculation flow region in the vicinity of the blockage that was created by the interaction between inflows from neighboring subchannels. Spectral analysis was performed and the vortex shedding frequency indicated by the Strouhal number was found to be St = 0.16 for the two Reynolds numbers investigated. The flow characteristics of the interior subchannels were investigated via the two-point cross-correlation of fluctuating velocities. In addition, proper orthogonal decomposition (POD) analysis was applied to the instantaneous vorticity fields to extract the statistically dominant flow structures. Lastly, the statistical characteristics of the vortex, such as the vortex populations, spatial distributions, and vortex strengths, were acquired by combining POD analysis and the vortex identification algorithm.

42 ENGINEERING↗
DOE OSTI2023

Evaluation of liquid cloud albedo susceptibility in E3SM using coupled eastern North Atlantic surface and satellite retrievals

The impact of aerosol number concentration on cloud albedo is a persistent source of spread in global climate predictions due to multi-scale, interactive atmospheric processes that remain difficult to quantify. We use 5 years of geostationary satellite and surface retrievals at the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) eastern North Atlantic (ENA) site in the Azores to evaluate the representation of liquid cloud albedo susceptibility for overcast cloud scenes in the DOE Energy Exascale Earth System Model version 1 (E3SMv1) and provide possible reasons for model–observation discrepancies. The overall distribution of surface 0.2 % CCN concentration values is reasonably simulated, but simulated liquid water path (LWP) is lower than observed and layer mean droplet concentration (N d ) comparisons are highly variable depending on the N d retrieval technique. E3SMv1's cloud albedo is greater than observed for given LWP and N d values due to a lower cloud effective radius than observed. However, the simulated albedo response to N d is suppressed due to a correlation between the solar zenith angle (SZA) and N d created by the seasonal cycle that is not observed. Controlling for this effect by examining the cloud optical depth (COD) shows that E3SMv1's COD response to CCN concentration is greater than observed. For surface-based retrievals, this is only true after controlling for cloud adiabaticity because E3SMv1's adiabaticities are much lower than observed. Assuming a constant adiabaticity in surface retrievals as done in top-of-atmosphere (TOA) retrievals narrows the retrieved ln N d distribution, which increases the cloud albedo sensitivity to ln N d to match the TOA sensitivity. The greater sensitivity of COD to CCN is caused by a greater Twomey effect in which the sensitivity of N d to CCN is greater than observed for TOA-retrieved N d , and once model–observation cloud adiabaticity differences are removed, this is also true for surface-retrieved N d . The LWP response to N d in E3SMv1 is overall negative as observed. Despite reproducing the observed LWP–N d relationship, observed clouds become much more adiabatic as N d increases, while E3SMv1 clouds do not, associated with more heavily precipitating clouds that are partially but not completely caused by deeper clouds and weaker inversions in E3SMv1. These cloud property differences indicate that the negative LWP–N d relationship is likely not caused by the same mechanisms in E3SMv1 and observations. The negative simulated LWP response also fails to mute the excessively strong Twomey effect, highlighting potentially important confounding factor effects that likely render the LWP–N d relationship non-causal. N d retrieval scales and assumptions, particularly related to cloud adiabaticity, contribute to substantial spreads in the model–observation comparisons, though enough consistency exists to suggest that aerosol activation, drizzle, and entrainment processes are critical areas to focus E3SMv1 development for improving the fidelity of aerosol–cloud interactions in E3SM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The LSST AGN Data Challenge: Selection Methods

Abstract Development of the Rubin Observatory Legacy Survey of Space and Time (LSST) includes a series of Data Challenges (DCs) arranged by various LSST Scientific Collaborations that are taking place during the project's preoperational phase. The AGN Science Collaboration Data Challenge (AGNSC-DC) is a partial prototype of the expected LSST data on active galactic nuclei (AGNs), aimed at validating machine learning approaches for AGN selection and characterization in large surveys like LSST. The AGNSC-DC took place in 2021, focusing on accuracy, robustness, and scalability. The training and the blinded data sets were constructed to mimic the future LSST release catalogs using the data from the Sloan Digital Sky Survey Stripe 82 region and the XMM-Newton Large Scale Structure Survey region. Data features were divided into astrometry, photometry, color, morphology, redshift, and class label with the addition of variability features and images. We present the results of four submitted solutions to DCs using both classical and machine learning methods. We systematically test the performance of supervised models (support vector machine, random forest, extreme gradient boosting, artificial neural network, convolutional neural network) and unsupervised ones (deep embedding clustering) when applied to the problem of classifying/clustering sources as stars, galaxies, or AGNs. We obtained classification accuracy of 97.5% for supervised models and clustering accuracy of 96.0% for unsupervised ones and 95.0% with a classic approach for a blinded data set. We find that variability features significantly improve the accuracy of the trained models, and correlation analysis among different bands enables a fast and inexpensive first-order selection of quasar candidates.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Harvest Intensity Effects on Carbon Stocks and Biodiversity Are Dependent on Regional Climate in Douglas-Fir Forests of British Columbia

Temperate forests provide crucial ecosystems services as living sinks for atmospheric carbon (C) and repositories of biodiversity. Applying harvesting at intensities that minimize losses offers one means for mitigating global change. However, little is known of overstory retention levels that best conserve ecosystem services in different regional climates, and likewise as climate changes. To quantify the effect of harvest intensity on C stocks and biodiversity, we compared five harvesting intensities (clearcutting, seedtree retention, 30% patch retention, 60% patch retention, and uncut controls) across a climatic aridity gradient that ranged from humid to semi-arid in the Douglas-fir (Pseudotsuga menziesii) forests of British Columbia. We found that increased harvesting intensity reduced total ecosystem, above ground, and live tree C stocks 1 year post-harvest, and the magnitude of these losses were negatively correlated with climatic aridity. In humid forests, total ecosystem C ranged from 50% loss following clearcut harvest, to 30% loss following large patch retention harvest. In arid forests this range was 60 to 8% loss, respectively. Where lower retention harvests are sought, the small patch retention treatment protected both C stocks and biodiversity in the arid forests, whereas the seedtree method performed as well or better in the humid forests. Below ground C stocks declined by an average of 29% after harvesting, with almost all of the loss from the forest floor and none from the mineral soil. Of the secondary pools, standing and coarse deadwood declined in all harvesting treatments regardless of cutting intensity or aridity, while C stocks in fine fuels and stumps increased. The understory plant C pool declined across all harvesting intensities in the humid forests, but increased in arid forests. Shannon’s diversity and richness of tree and bryoid species declined with harvesting intensity, where tree species losses were greatest in the humid forests and bryoid losses greatest in arid forests. Shrub and herb species were unaffected. This study showed that the highest retention level was best at reducing losses in C stocks and biodiversity, and clearcutting the poorest, and while partial retention of canopy trees can reduce losses in these ecosystem services, outcomes will vary with climatic aridity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Predicting Antimicrobial Resistance Using Partial Genome Alignments

Antimicrobial resistance (AMR) is an important global health threat that impacts millions of people worldwide each year. Developing methods that can detect and predict AMR phenotypes can help to mitigate the spread of AMR by informing clinical decision making and appropriate mitigation strategies. Many bioinformatic methods have been developed for predicting AMR phenotypes from whole-genome sequences and AMR genes, but recent studies have indicated that predictions can be made from incomplete genome sequence data. In order to more systematically understand this, we built random forest-based machine learning classifiers for predicting susceptible and resistant phenotypes for Klebsiella pneumoniae (1,640 strains), Mycobacterium tuberculosis (2,497 strains), and Salmonella enterica (1,981 strains). We started by building models from alignments that were based on a reference chromosome for each species. We then subsampled each chromosomal alignment and built models for the resulting subalignments, finding that very small regions, representing approximately 0.1 to 0.2% of the chromosome, are predictive. In K. pneumoniae, M. tuberculosis, and S. enterica, the subalignments are able to predict multiple AMR phenotypes with at least 70% accuracy, even though most do not encode an AMR-related function. We used these models to identify regions of the chromosome with high and low predictive signals. Finally, subalignments that retain high accuracy across larger phylogenetic distances were examined in greater detail, revealing genes and intergenic regions with potential links to AMR, virulence, transport, and survival under stress conditions. IMPORTANCE Antimicrobial resistance causes thousands of deaths annually worldwide. Understanding the regions of the genome that are involved in antimicrobial resistance is important for developing mitigation strategies and preventing transmission. Machine learning models are capable of predicting antimicrobial resistance phenotypes from bacterial genome sequence data by identifying resistance genes, mutations, and other correlated features. They are also capable of implicating regions of the genome that have not been previously characterized as being involved in resistance. In this study, we generated global chromosomal alignments for Klebsiella pneumoniae, Mycobacterium tuberculosis, and Salmonella enterica and systematically searched them for small conserved regions of the genome that enable the prediction of antimicrobial resistance phenotypes. In addition to known antimicrobial resistance genes, this analysis identified genes involved in virulence and transport functions, as well as many genes with no previous implication in antimicrobial resistance.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1990

A study of the topology of dissipating motions in direct numerical simulations of time-developing compressible and incompressible mixing layers

A preliminary investigation of the geometry of flow patterns in numerically simulated compressible and incompressible mixing layers was carried out using 3-D critical point methodology. Motions characterized by high rates of kinetic energy dissipation and/or high enstrophy were of particular interest. In the approach the partial derivatives of the velocity field are determined at every point in the flow. These are used to construct the invariants of the velocity gradient tensor and the rate-of-strain tensor (P, Q, R, and P(sub s), Q(sub s), R(sub s) respectively). For incompressible flow the first invariant is zero. For the conditions of the compressible simulation, the first invariant is found to be everywhere small, relative to the second and third invariants, and so in both cases the local topology at a point is mainly determined by the second and third invariants. The data at every grid point is used to construct scatter plots of Q versus R and Q(sub s) versus R(sub s). Most points map to a cluster near the origin in Q-R space. However, fine scale motions, that is motions which are characterized by velocity derivatives which scale with the square root of R(sub delta), tend to map to regions which lie far from the origin. Definite trends are observed for motions characterized by high enstrophy and/or high dissipation. The observed trends suggest that, for these motions, the second and third invariants of the velocity gradient and rate-of-strain tensors are strongly correlated. The second and third invariants of the rate-of-strain tensor are related by K(-Q(sub s))(exp 3/2), which is consistent with the above scaling of velocity derivatives. The quantity K appears to depend on Reynolds number with an upper limit K = 2(the square root of 3)/9 corresponding to locally axisymmetric flow. For both the compressible and incompressible mixing layer, regions corresponding to high rates of dissipation are found to be characterized by comparable magnitudes of R(sub ij)R(sub ij) and S(sub ij)S(sub ij). For the incompressible mixing layer, regions characterized by the highest values of enstrophy are found to have relatively low strain rates.

Chen, J. H.↗
DOE OSTI2018

Alleviation of fatty liver in a rat model by enhancing N 1-methylnicotinamide bioavailability through aldehyde oxidase inhibition

Highlights: • A slight positive correlation between human liver MNA and NNMT protein level. • Increase of liver MNA concentrations in rat alleviate fatty liver. • Inhibition of liver MNA metabolism would lead to fatty liver amelioration. • Hydralazine might be suitable for the combined treatment of hypertension and NAFLD. Nonalcoholic fatty liver disease (NAFLD) has increased worldwide in recent years. NAFLD is classified into two types, nonalcoholic fatty liver (NAFL), with few complications, and nonalcoholic steatohepatitis (NASH), which leads to liver cirrhosis or cancer. This study was based on previous reports that N{sup 1}-methylnicotinamide (MNA) can stabilise sirtuin 1 protein, leading to decreased lipid levels in the liver. We hypothesised that fatty liver improvement by MNA would be further enhanced by suppressing its rapid metabolism by aldehyde oxidase in the liver. To test this, hydralazine (HYD), a potent aldehyde oxidase inhibitor, was administered orally to NAFL model rats. Liver triglyceride (TG) levels in the model were nearly unchanged by administration of MNA alone. In contrast, TG levels were marked decreased in NAFL rats treated with a combination of MNA and HYD. In addition, TG levels were decreased even in NAFL rats treated with only HYD. These findings supported our hypothesis that maintaining MNA concentrations in the liver, by suppressing MNA metabolism, would at least partially ameliorate fatty liver.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1994

Lick optical spectra of quasar HS 1946+7658 at 10 kilometers per second resolution Lyman-alpha forest and metal absorption systems

We present optical spectra of the most luminous known quasi stellar object (QSO) HS 1946+7658 (z(sub em) = 3.051). Our spectra have both full wavelength coverage, 3240-10570 A, and in selected regions, either high signal-to-noise ratio, SNR approximately equals 40-100, or unusually high approximately 10 km/sec resolution, and in parts of the Ly alpha forest and to the red of Ly alpha emission they are among the best published. We find 113 Ly alpha systems and six metal-line systems, three of which are new. The metal systems at z(sub abs) = 2.844 and 3.050 have complex velocity structure with four and three prominent components, respectively. We find that the system at z(sub abs) = 2.844 is a damped Ly alpha absorption (DLA) system, with a neutral hydrogen column density of log N(H I) = 20.2 +/- 0.4, and it is the cause of the Lyman limit break at lambda approximately equals 3520 A. We believe that most of the H I column density in this system is in z(sub abs) = 2.8443 component which shows the strongest low-ionization absorption lines. The metal abundance in the gas phase of the system is (M/H) approximately equals -2.6 +/- 0.3, with a best estimate of (M/H) = -2.8, with ionizaion parameter log gamma = -2.75, from a photoionization model. The ratios of the logarithmic abundances of C, O, Al, and Si are all within a factor of 2 of solar, which is important for two reasons. First, we believe that the gas abundances which we measure are close to the total abundances, because the ratio of aluminum to other elements is near cosmic, and Al is a refractory element which depletes very readily like chromium, in the interstellar medium. Second, we do not see the enhancement of O with respect to C of (O/C) approximately equals 0.5-0.9 reported in three partial Lyman limit systems by Reimers et al. (1992) and Vogel & Reimers (1993); we measure (O/C) = -0.06 for observed ions and (O/C) approximately equals 0.2 after ionization corrections, which is consistent with solar abundances. We see C II*(lambda 1335) offset by 15 km/sec with respect to C II(lambda 1334), presumably because the gas density varies from 2 to 8 cm(exp -3) with changing velocity in the DLA system. These densities imply that the damped component is 6-25 pc thick, which is reasonable for a single cloud in a cold spiral disk. They also imply that the cloud is relatively highly ionized with more C III than C II, more O III than O I, and log N(H I) = 20.72, which is 3 times the H I column. The system at z(abs) = 1.7382 is also believed to be damped with N(H I) approximately equals 10(exp 21) cm(exp -2), because we see Cr II, but its Ly alpha line will never be seen because it is below the Lyman limit of the other DLA system. We see a 2.6 sigma lack of Ly alpha forest lines well away from the QSO redshift, which may be a chance fluctuation. We also see a correlation between column density N(H I) and Doppler parameter b for 96 unsaturated Ly Alpha forest absorption lines, and although this correlation persists in the 36 Ly alpha lines which lie in regions where the SNR approximately equals 8-16, we agree with Rauch et al. (1993) that it is probably a bogus effect of low supernova remnant (SNR). The same applies to lines with very low b values: in regions where SNR less than or equal to 8 we see many Ly alpha lines which appear to have 10 less than or equal to b less than or equal to 20, but when 8 less than or equal to SNR less than or equal to 16 we see only one line with b less than or equal to 15 km/sec, and two others which we believe have b less than or equal to 20, with values of 20 and 16 km/sec. Traditional Ly alpha line samples which include all lines which have W/sigma(W) greater than or equal to 4 are not adequate to explore the distribution of the properties of individual clouds, because we need much higher (W/sigma(W)) and SNR to avoid the strong biases.

Fan, Xiao-Ming↗
DOE OSTI2021

Airborne hyperspectral imaging of nitrogen deficiency on crop traits and yield of maize by machine learning and radiative transfer modeling

Nitrogen is an essential nutrient that directly affects plant photosynthesis, crop yield, and biomass production for bioenergy crops, but excessive application of nitrogen fertilizers can cause environmental degradation. To achieve sustainable nitrogen fertilizer management for precision agriculture, there is an urgent need for nondestructive and high spatial resolution monitoring of crop nitrogen and its allocation to photosynthetic proteins as that changes over time. Here, we used visible to shortwave infrared (400–2400 nm) airborne hyperspectral imaging with high spatial (0.5 m) and spectral (3–5 nm) resolutions to accurately estimate critical crop traits, i.e., nitrogen, chlorophyll, and photosynthetic capacity (CO 2 -saturated photosynthesis rate, V max,27 ), at leaf and canopy scales, and to assess nitrogen deficiency on crop yield. We conducted three airborne campaigns over a maize (Zea mays L.) field during the growing season of 2019. Physically based soil-canopy Radiative Transfer Modeling (RTM) and data-driven approaches i.e. Partial-Least Squares Regression (PLSR) were used to retrieve crop traits from hyperspectral reflectance, with ground truth of leaf nitrogen, chlorophyll, V max,27 , Leaf Area Index (LAI), and harvested grain yield. To improve computational efficiency of RTMs, Random Forest (RF) was used to mimic RTM simulations to generate machine learning surrogate models RTM-RF. The results show that prior knowledge of soil background and leaf angle distribution can significantly reduce the ill-posed RTM retrieval. RTM-RF achieved a high accuracy to predict leaf chlorophyll content (R 2 = 0.73) and LAI (R 2 = 0.75). Meanwhile, PLSR exhibited better accuracy to predict leaf chlorophyll content (R 2 = 0.79), nitrogen concentration (R 2 = 0.83), nitrogen content (R 2 = 0.77), and V max,27 (R 2 = 0.69) but required measured traits for model training. We also found that canopy structure signals can enhance the use of spectral data to predict nitrogen related photosynthetic traits, as combining RTM-RF LAI and PLSR leaf traits well predicted canopy-level traits (leaf traits × LAI) including canopy chlorophyll (R 2 = 0.80), nitrogen (R 2 = 0.85) and V max,27 (R 2 = 0.82). Compared to leaf traits, we further found that canopy-level photosynthetic traits, particularly canopy V max,27 , have higher correlation with maize grain yield. This study highlights the potential for synergistic use of process-based and data-driven approaches of hyperspectral imaging to quantify crop traits that facilitate precision agricultural management to secure food and bioenergy production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Understanding latent timescales in neural ordinary differential equation models of advection-dominated dynamical systems

The neural ordinary differential equation (ODE) framework has shown considerable promise in recent years in developing highly accelerated surrogate models for complex physical systems characterized by partial differential equations (PDEs). For PDE-based systems, state-of-the-art neural ODE strategies leverage a two-step procedure to achieve this acceleration: a nonlinear dimensionality reduction step provided by an autoencoder, and a time integration step provided by a neural-network based model for the resultant latent space dynamics (the neural ODE). This work explores the applicability of such autoencoder-based neural ODE strategies for PDEs in which advection terms play a critical role. More specifically, alongside predictive demonstrations, physical insight into the sources of model acceleration (i.e., how the neural ODE achieves its acceleration) is the scope of the current study. Such investigations are performed by quantifying the effects of both autoencoder and neural ODE components on latent system time-scales using eigenvalue analysis of dynamical system Jacobians. To this end, the sensitivity of various critical training parameters – de-coupled versus end-to-end training, latent space dimensionality, and the role of training trajectory length, for example – to both model accuracy and the discovered latent system timescales is quantified. Furthermore, this work specifically uncovers the key role played by the training trajectory length (the number of rollout steps in the loss function during training) on the latent system timescales: larger trajectory lengths correlate with an increase in limiting neural ODE time-scales, and optimal neural ODEs are found to recover the largest time-scales of the full-order (ground-truth) system. Demonstrations are performed across fundamentally different unsteady fluid dynamics configurations influenced by advection: (1) the Kuramoto–Sivashinsky equations (2) Hydrogen-Air channel detonations (the compressible reacting Navier–Stokes equations with detailed chemistry), and (3) 2D Atmospheric flow.

Advection-dominated dynamical systems↗
DOE OSTI2023

Modeling Advanced Strippers for CO2 Capture from Gas-Fired Power Plants using Aqueous Piperazine

Energy performance of 5 m piperazine was evaluated with an advanced stripper at the pilot scale, accounting for heat loss. Modeling studies indicated that heat loss at locations other than the heat source impacted the heat duty. For strippers with excess packing, the column was the most important source of heat loss, and values as low as 0.03 GJ/hr can cause pinch and reabsorption of CO2. Solvents like 5 m PZ are more affected by heat loss in the column due to a top-side temperature pinch at high lean loading. Heat loss impacts NGCC CO2 capture more than coal-based CO2 capture. Surface temperature measurements pilot plants showed that heat loss was 50-75% controlled by natural convection. Measured heat loss was correlated with steam flowrate and wind speed and ranged from 8126 to 219668 Btu/hr, representing 35% of the measured heat rate on average. Net heat duty linearly varied with CO2 removal at low lean loading and was similar for coal and NGCC conditions at fixed removal. At 90% removal, net heat duty was about 2.5 GJ/t with 4 and 12% CO2 in the flue gas. 7 At NGCC conditions, strippers with finite packing can benefit from a flashing feed to the top of the column by using a hot bypass. A flashing feed was linked to a reduction in irreversibility of the stripper from temperature driving forces. It also reduces heat duty by up to 6% at fixed packing height and reduces packing requirement for a fixed steam heater size compared to a warm bypass. Lean loading and lean solvent rate were effective handles to maximize profitability of a fixed stripper design at low gas price. High ambient temperature operation can benefit from low pressure stripping to about 0.18 lean loading at 150 ℃ to maintain the cyclic capacity at reduced rich loading. At low gas price, the capture plant was able to maximize profitability even with a high heat duty of 2.7-3.1 GJ/t at a lean loading of 0.2-0.22 mol/mol. Heat recovery by partial water vapor condensation in a gas-liquid exchanger can be replaced by a direct contact condenser (DCC). The DCC improved gas cooling and reduced heat duty compared to the CO2 exchanger at 0.2 lean loading. The DCC when used with a 150 ℃/5.5 bar stripper can reduce the cost of capture compared to the base case by $3- 4/tonne at 15 ft of packing (optimum) but increases cost of capture by $6/tonne with a 120 ℃/2 bar stripper and 10 ft of packing (optimum). The DCC can work effectively with a solvent with a high heat of absorption and thermal stability.

Suresh Babu, Athreya↗
NASA NTRS2010

An Experimental Study of Small-Scale Variability of Raindrop Size Distribution

An experimental study of small-scale variability of raindrop size distributions (DSDs) has been carried out at Wallops Island, Virginia. Three Joss-Waldvogel disdrometers were operated at a distance of 0.65, 1.05, and 1.70 km in a nearly straight line. The main purpose of the study was to examine the variability of DSDs and its integral parameters of liquid water content, rainfall, and reflectivity within a 2-km array: a typical size of Cartesian radar pixel. The composite DSD of rain events showed very good agreement among the disdrometers except where there were noticeable differences in midsize and large drops in a few events. For consideration of partial beam filling where the radar pixel was not completely covered by rain, a single disdrometer reported just over 10% more rainy minutes than the rainy minutes when all three disdrometers reported rainfall. Similarly two out of three disdrometers reported5%more rainy minutes than when all three were reporting rainfall. These percentages were based on a 1-min average, and were less for longer averaging periods. Considering only the minutes when all three disdrometers were reporting rainfall, just over one quarter of the observations showed an increase in the difference in rainfall with distance. This finding was based on a 15-min average and was even less for shorter averaging periods. The probability and cumulative distributions of a gamma-fitted DSD and integral rain parameters between the three disdrometers had a very good agreement and no major variability. This was mainly due to the high percentage of light stratiform rain and to the number of storms that traveled along the track of the disdrometers. At a fixed time step, however, both DSDs and integral rain parameters showed substantial variability. The standard deviation (SD) of rain rate was near 3 mm/h, while the SD of reflectivity exceeded 3 dBZ at the longest separation distance. These standard deviations were at 6-min average and were higher at shorter averaging periods. The correlations decreased with increasing separation distance. For rain rate, the correlations were higher than previous gauge-based studies. This was attributed to the differences in data processing and the difference in rainfall characteristics in different climate regions. It was also considered that the gauge sampling errors could be a factor. In this regard, gauge measurements were simulated employing existing disdrometer dataset. While a difference was noticed in cumulative distribution of rain occurrence between the simulated gauge and disdrometer observations, the correlations in simulated gauge measurements did not differ from the disdrometer measurements.

Tokay, Ali↗
NASA NTRS1995

A method for combining passive microwave and infrared rainfall observations

Because passive microwave instruments are confined to polar-orbiting satellites, rainfall estimates must interpolate across long time periods, during which no measurements are available. In this paper the authors discuss a technique that allows one to partially overcome the sampling limitations by using frequent infrared observations from geosynchronous platforms. To accomplish this, the technique compares all coincident microwave and infrared observations. From each coincident pair, the infrared temperature threshold is selected that corresponds to an area equal to the raining area observed in the microwave image. The mean conditional rainfall rate as determined from the microwave image is then assigned to pixels in the infrared image that are colder than the selected threshold. The calibration is also applied to a fixed threshold of 235 K for comparison with established infrared techniques. Once a calibration is determined, it is applied to all infrared images. Monthly accumulations for both methods are then obtained by summing rainfall from all available infrared images. Two examples are used to evaluate the performance of the technique. The first consists of a one-month period (February 1988) over Darwin, Australia, where good validation data are available from radar and rain gauges. For this case it was found that the technique approximately doubled the rain inferred by the microwave method alone and produced exceptional agreement with the validation data. The second example involved comparisons with atoll rain gauges in the western Pacific for June 1989. Results here are overshadowed by the fact that the hourly infrared estimates from established techniques, by themselves, produced very good correlations with the rain gauges. The calibration technique was not able to improve upon these results.

Kummerow, Christian↗
NASA NTRS1992

Computational studies of metal-metal and metal-ligand interactions

The geometric structure of Cr(CO)6 is optimized at the modified coupled-pair functional (MCPF), single and double excitation coupled-cluster (CCSD) and CCSD(T) levels of theory (including a perturbational estimate for connected triple excitations), and the force constants for the totally symmetric representation are determined. The geometry of Cr(CO)5 is partially optimized at the MCPF, CCSD and CCSD(T) levels of theory. Comparison with experimental data shows that the CCSD(T) method gives the best results for the structures and force constants, and that remaining errors are probably due to deficiencies in the one-particle basis sets used for CO. A detailed comparison of the properties of free CO is therefore given, at both the MCPF and CCSD/CCSD(T) levels of treatment, using a variety of basis sets. With very large one-particle basis sets, the SSCD(T) method gives excellent results for the bond distance, dipole moment and harmonic frequency of free CO. The total binding energies of Cr(CO)6 and Cr(CO)5 are also determined at the MCPF, CCSD and CCSD(T) levels of theory. The CCSD(T) method gives a much larger total binding energy than either the MCPF or CCSD methods. An analysis of the basis set superposition error (BSSE) at the MCPF level of treatment points out limitations in the one-particle basis used here and in a previous study. Calculations using larger basis sets reduced the BSSE, but the total binding energy of Cr(CO)6 is still significantly smaller than the experimental value, although the first CO bond dissociation energy of Cr(CO)6 is well described. An investigation of 3s3p correlation reveals only a small effect. The remaining discrepancy between the experimental and theoretical total binding energy of Cr(CO)6 is probably due to limitations in the one-particle basis, rather than limitations in the correlation treatment. In particular an additional d function and an f function on each C and O are needed to obtain quantitative results. This is underscored by the fact that even using a very large primitive se (1042 primitive functions contracted to 300 basis functions), the superposition error for the total binding energy of Cr(CO)6 is 22 kcal/mol at the MCPF level of treatment.

Barnes, Leslie A.↗
NASA NTRS2012

Decorrelation Times of Photospheric Fields and Flows

We use autocorrelation to investigate evolution in flow fields inferred by applying Fourier Local Correlation Tracking (FLCT) to a sequence of high-resolution (0.3 "), high-cadence (approx = 2 min) line-of-sight magnetograms of NOAA active region (AR) 10930 recorded by the Narrowband Filter Imager (NFI) of the Solar Optical Telescope (SOT) aboard the Hinode satellite over 12 - 13 December 2006. To baseline the timescales of flow evolution, we also autocorrelated the magnetograms, at several spatial binnings, to characterize the lifetimes of active region magnetic structures versus spatial scale. Autocorrelation of flow maps can be used to optimize tracking parameters, to understand tracking algorithms f susceptibility to noise, and to estimate flow lifetimes. Tracking parameters varied include: time interval Delta t between magnetogram pairs tracked, spatial binning applied to the magnetograms, and windowing parameter sigma used in FLCT. Flow structures vary over a range of spatial and temporal scales (including unresolved scales), so tracked flows represent a local average of the flow over a particular range of space and time. We define flow lifetime to be the flow decorrelation time, tau . For Delta t > tau, tracking results represent the average velocity over one or more flow lifetimes. We analyze lifetimes of flow components, divergences, and curls as functions of magnetic field strength and spatial scale. We find a significant trend of increasing lifetimes of flow components, divergences, and curls with field strength, consistent with Lorentz forces partially governing flows in the active photosphere, as well as strong trends of increasing flow lifetime and decreasing magnitudes with increases in both spatial scale and Delta t.

Welsch, B. T.↗