Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model interpretability”

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Aerosol Thermodynamics (Aerosol Liquid Water Concentration and Aerosol pH) at S2 during CoURAGE

These data report aerosol liquid water content (ALWC) and aerosol pH for the S2 (Mt. Airy) site during CoURAGE for April - June 2025. The ISORROPIA-II aerosol thermodynamic equilibrium model was used to compute ALWC and pH. The model was run in forward mode, with solids formation disabled (metastable mode), according to Pye et al. (Atmospheric Chemistry and Physics, 2020). Information about the model can be found at https://www.epfl.ch/labs/lapi/models-and-software/isorropia/, and in Fountoukis and Nenes (Atmospheric Chemistry and Physics, 2007). Inputs to the model were: (1) measured Temperature and relative humidity, both from the DoE ARM “metwxt” data products at S2, (2) aerosol chemical composition from the ARM ACSM, and (3) gas-phase ammonia concentrations. For ISORROPIA-II, ambient RH values above 0.995 were input to the model as 0.995; therefore, caution should be exercised interpreting ALWC when RH was above 0.995. Aerosol composition or NH3(g) concentrations below the LOD were input into the model as 0.5*LOD. Model runs were only conducted for data points at 30-min resolution with all three inputs (meteorology, ACSM, and NH3). If any of the three model inputs were missing, aerosol thermodynamic outputs are flagged as -999. Similarly, aerosol thermodynamic outputs are flagged as -999 if measured NH3(g) and aerosol ammonium (NH4+) were simultaneously below the respective LODs. The ALWC includes liquid water from inorganic aerosol components, calculated directly by ISORROPIA-II, as well as liquid water from organics, estimated from the ACSM organics assuming a kappa value of 0.12 according to Guo et al. (Atmospheric Chemistry and Physics, 2015).

Aerosol Liquid Water Content↗
NASA NTRS1993

Inference of cirrus cloud properties using satellite-observed visible and infrared radiances. II - Verification of theoretical cirrus radiative properties

The study presents a methodology developed to apply a parameterization of radiative transfer calculations to satellite analyses of cirrus clouds. Cloud heights and optical depths are derived from visible and IR window measurements taken during FIRE when cirrus clouds were present. Geostationary satellite retrievals are compared to lidar-derived cloud heights and retrievals from a polar-orbiting satellite taken at different angles to determine which theoretical models of scattering phase function and single-scattering albedo best represent actual cirrus clouds. Models using small scattering phase function and single-scattering albedo best represent actual cirrus clouds. It is concluded that interpretation of cirrus reflectance with water-droplet models leads to biased results. The cloud-height and optical depth biases can be minimized with the aid of the C20 or cirrostratus models.

Minnis, Patrick↗
NASA NTRS2021

Mid-Pliocene Atlantic Meridional Overturning Circulation simulated in PlioMIP2

In the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2), coupled climate models have been used to simulate an interglacial climate during the mid-Piacenzian warm period (mPWP; 3.264 to 3.025 Ma). Here, we compare the Atlantic Meridional Overturning Circulation (AMOC), poleward ocean heat transport and sea surface warming in the Atlantic simulated with these models. In PlioMIP2, all models simulate an intensified mid-Pliocene AMOC. However, there is no consistent response in the simulated Atlantic ocean heat transport nor in the depth of the Atlantic overturning cell. The models show a large spread in the simulated AMOC maximum, the Atlantic ocean heat transport and the surface warming in the North Atlantic. Although a few models simulate a surface warming of ~8–12°C in the North Atlantic, similar to the reconstruction from Pliocene Research, Interpretation and Synoptic Mapping (PRISM) version 4, most models appear to underestimate this warming. The large model spread and model–data discrepancies in the PlioMIP2 ensemble do not support the hypothesis that an intensification of the AMOC, together with an increase in northward ocean heat transport, is the dominant mechanism for the mid-Pliocene warm climate over the North Atlantic.

Pliocene Model Intercomparison Project Phase 2 (Pl↗
DOE OSTI2024

A Science Gateway for the Repeatable Analysis of Machine Learning Predicted Gravity Anomalies

In recent years, deep learning has become an increasingly popular alternative for modeling in geoscience applications due to its scalability and efficiency. However, the interpretability, compute, data volume, and hyperparameter tuning requirements of deep learning models make development and monitoring difficult. Furthermore, model explainability and communicating results obtained by these models to users or domain experts is a challenge, as domain experts in geoscience also need to have a deep understanding of how those models function in order to support their scientific works. Here, we describe a science gateway and machine learning pipeline for predicting gravity anomalies from geophysical data. The gateway, built on open-source technologies, provides a holistic view of the pipeline through interactive visualizations aimed at enabling efficient exploratory data analysis. The repeatability, reproducibility, and monitoring capabilities of this overall system allow us to iterate and analyze at scale. Using this pipeline and gateway, we can repeatedly produce accurate high-resolution gravity anomaly datasets. By describing the underlying technologies, implementation, and results, here we provide a foundation for the broader adoption of science gateways into cross-cutting geoscience and machine learning research projects as a means to improve the scientific discovery and collaboration in the geophysics and computational sciences community.

58 GEOSCIENCES↗
DOE OSTI2022

Estimation of the NiCu Cycle Strength and Its Impact on Type I X-Ray Bursts

Abstract Type I X-ray bursts (XRBs) are powered by thermonuclear burning on proton-rich unstable nuclides. The construction of burst models with accurate knowledge of nuclear physics is required to properly interpret burst observations. Numerous studies that have investigated the sensitivities of burst models to nuclear inputs have commonly extracted the strength of the NiCu cycle in the rp process, determined by the 59 Cu( p , α ) 56 Ni and 59 Cu( p , γ ) 60 Zn thermonuclear reaction rates, as critical in the determination of reaction flow in the burst. In this study, the strength of the cycle at the XRB temperature range was estimated based on published experimental data. The nuclear properties of the compound nucleus 60 Zn were evaluated for the 59 Cu( p , α ) 56 Ni and 59 Cu( p , γ ) 60 Zn reaction rate calculations. Monte Carlo rate calculations were conducted to include the large uncertainties of nuclear properties in the calculations. In the current work, a weak NiCu cycle is expected, whereas the rates adopted by the previous studies suggest a strong NiCu cycle. Model simulations were performed with the new rates to assess the impact on Type I XRBs. The results show that the estimated cycle strength does not strongly influence the model predictions of the burst light curve or synthesized abundances.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Observational evidence of preferred flow regimes in the Northern Hemisphere winter stratosphere

Ten years of stratospheric geopotential height data are analyzed in an attempt to determine whether there are preferred flow regimes in the Northern Hemisphere winter stratosphere. The data are taken from Stratospheric Sounding Units on board NOAA satellites. The probability density estimate of the amplitude of the wavenumber 1 10-mb height is found to be bimodal. The density distribution is composed of a dominant large-amplitude mode and a less frequent low-amplitude mode. When the wavenumber 1 10-mb height data are projected onto the phase plane defined by the 10-mb zonal-mean winds and wavenumber 1 100-mb heights, three preferred regimes are evident. The small-amplitude mode separates into a strong zonal wind-weak wave regime and a weak zonal wind-weak wave regime. The large-amplitude mode is an intermediate zonal wind-strong wave regime. Transitions between the large-amplitude regime and the weak zonal wind-weak wave regime are found to be associated with major stratospheric warmings. The clustering of the stratospheric data into the preferred flow regimes is interpreted in light of the bifurcation properties of the Holton and Mass model. The interannual variability of the Northern Hemisphere winter stratosphere is interpreted in terms of the relative frequency of the observed preferred regimes.

Pierce, R. B.↗
NASA NTRS1994

Interpretation of satellite airglow observations during the March 22, 1979, magnetic storm, using the coupled ionosphere-thermosphere model developed at University College, London

The University of California, Berkeley, extreme ultraviolet spectrometer aboard the U.S. Air Force STP 78-1 satellite measured emission features in the Earth's dayglow due to neutral and ionized species in the atmosphere, in the 35 to 140-nm range. The spectrometer was operating between March 1979 and March 1980, including the period of the magnetic storm on March 22, 1979. Some of these measurements are interpreted using the predictions of the three-dimensional time-dependent coupled ionosphere-thermosphere model developed at University College, London. The observations show a reduction in the atomic oxygen 130.4-nm airglow emission at high northern latitudes following the storm. Model simulations show that this reduction in 130.4-nm emission is associated with an increase in the O2/O ratio. Analysis of model results using electron transport and radiative transport codes show that the brightness of 130.4-nm emission at high latitudes due to resonantly scattered sunlight is approximately twice that due to photoelectron impact excitation. However, the observed decrease in the brightness at high northern latitudes is mainly due to a change in the photoelectron impact source, which contributes approximately 75% of the total, as well as its multiple scattering component; for the photoelectron impact source at 70 deg latitude and 200 km altitude, the reduction in multiple scattering is 1.5 times greater than the reduction in the initial excitation. The reduction in the airglow emission is visible only in the norther n hemisphere because the south pole was not sunlit over the storm period. The comparison of model results with observations suggests that 130.4-nm emission may be useful as a tracer for global changes in the concentration of atomic energy.

Parish, H. F.↗
DOE OSTI2023

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Interpretable AI forecasting for numerical relativity waveforms of quasicircular, spinning, nonprecessing binary black hole mergers

We present a deep-learning artificial intelligence model (AI) that is capable of learning and forecasting the late-inspiral, merger and ringdown of numerical relativity waveforms that describe quasicircular, spinning, nonprecessing binary black hole mergers. We used the NRHybSur3dq8 surrogate model to produce train, validation and test sets of ℓ = |m| = 2 waveforms that cover the parameter space of binary black hole mergers with mass ratios q ≤ 8 and individual spins |s$^{z}_ {(1,2)}$| ≤ 0.8. These waveforms cover the time range t ∊ [-5000 M, 130 M], where t = 0M marks the merger event, defined as the maximum value of the waveform amplitude. We harnessed the ThetaGPU supercomputer at the Argonne Leadership Computing Facility to train our AI model using a training set of 1.5 million waveforms. We used 16 NVIDIA DGX A100 nodes, each consisting of 8 NVIDIA A100 Tensor Core GPUs and 2 AMD Rome CPUs, to fully train our model within 3.5 h. Our findings show that artificial intelligence can accurately forecast the dynamical evolution of numerical relativity waveforms in the time range t ∊ [-100 M, 130 M]. Sampling a test set of 190,000 waveforms, we find that the average overlap between target and predicted waveforms is ≳99% over the entire parameter space under consideration. We also combined scientific visualization and accelerated computing to identify what components of our model take in knowledge from the early and late-time waveform evolution to accurately forecast the latter part of numerical relativity waveforms. This work aims to accelerate the creation of scalable, computationally efficient and interpretable artificial intelligence models for gravitational wave astrophysics

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Downstream Separation of Formic Acid with Anion-Exchange Resin from Electrocatalytic Carbon Dioxide (CO 2 ) Conversion: Adsorption, Kinetics, and Equilibrium Modeling

The electrocatalytic carbon dioxide reduction reaction (CO 2 RR) that produces liquid formic acid (FA) is one potential route to minimize the CO 2 concentration in the atmosphere by utilizing CO 2 as a feedstock. However, one of the main disadvantages of this route is the high energy demand required for purification through conventional distillation. This method of purification has been shown to result in a large carbon footprint in life-cycle analysis (LCA). Therefore, finding an alternative purification technique that can integrate with an electrochemical CO 2 RR process with the potential to lower the carbon footprint would be of great interest. Under the experimental conditions of the CO 2 RR, FA remains in its basic form as the formate anion in the catholyte. Anion exchange resins with their basic properties can separate the formate anion from the catholyte mixture as an alternative to energy intensive distillation. In this work, three anion exchange resins with increasing basicity order: Amberlite IRA-96 < Amberlite IRA-910 < Ambersep 900 are tested to separate FA in the presence of aqueous solution of potassium bicarbonate catholyte with varying levels of resin and FA. Kinetics and equilibrium studies data for the FA adsorption are interpreted using several kinetics and isotherm models. The kinetics data fit well with a pseudo-first-order model at high initial FA concentrations and a pseudo-second-order model at low initial FA concentrations. The experimental data can be best explained with the Freundlich isotherm model. The fitted results show that strongly basic Ambersep 900 is more effective at separating FA, with a maximum adsorption capacity of 336.7 mg/g compared to Amberlite IRA-96 (275.2 mg/g) and Amberlite IRA-910 (209.2 mg/g) in neutral water. However, weakly basic Amberlite IRA-96 (110.8 mg/g) can separate FA more efficiently from potassium bicarbonate catholyte than Amberlite IRA-910 (99.9 mg/g) and Ambersep 900 (40.0 mg/g). As a result, competitive adsorption occurs in the presence of potassium bicarbonate catholyte, but the formic acid adsorption capacity is still sufficient for adequate separation.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRS1992

Geometry and mass model of ionizing radiation experiments on the LDEF satellite

Extensive measurements related to ionizing radiation environments and effects were made on the LDEF satellite during its mission lifetime of almost 6 years. These data, together with the opportunity they provide for evaluating predictive models and analysis methods, should allow more accurate assessments of the space radiation environment and related effects for future missions in low Earth orbit. The LDEF radiation dosimetry data is influenced to varying degrees by material shielding effects due to the dosimeter itself, nearby components and experiments, and the spacecraft structure. A geometry and mass model is generated of LDEF, incorporating sufficient detail that it can be applied in determining the influence of material shielding on ionizing radiation measurements and predictions. This model can be used as an aid in data interpretation by unfolding shielding effects from the LDEF radiation dosimeter responses. Use of the LDEF geometry/mass model, in conjunction with predictions and comparisons with LDEF dosimetry data currently underway, will also allow more definitive evaluations of current radiation models for future mission applications.

Colborn, B. L.↗
NASA NTRS2004

Discharge Chamber Primary Electron Modeling Activities in Three-Dimensions

Designing discharge chambers for ion thrusters involves many geometric configuration decisions. Various decisions will impact discharge chamber performance with respect to propellant utilization efficiency, ion production costs, and grid lifetime. These hardware design decisions can benefit from the assistance of computational modeling. Computational modeling for discharge chambers has been limited to two-dimensional codes that leveraged symmetry for interpretation into three-dimensional analysis. This paper presents model development activities towards a three-dimensional discharge chamber simulation to aid discharge chamber design decisions. Specifically, of the many geometric configuration decisions toward attainment of a worthy discharge chamber, this paper focuses on addressing magnetic circuit considerations with a three-dimensional discharge chamber simulation as a tool. With this tool, candidate discharge chamber magnetic circuit designs can be analyzed computationally to gain insight into factors that may influence discharge chamber performance such as: primary electron loss width in magnetic cusps, cathode tip position with respect to the low magnetic field volume, definition of a low magnetic field region, and maintenance of a low magnetic field region across the grid span. Corroborating experimental data will be obtained from mockup hardware tests. Initially, simulated candidate magnetic circuit designs will resemble previous successful thruster designs. To provide opportunity to improve beyond previous performance benchmarks, off-design modifications will be simulated and experimentally tested.

Steuber, Thomas J.↗
DOE OSTI2023

An Interpretable Index of Social Vulnerability to Environmental Hazards

Index-based measures of social vulnerability to environmental hazards are commonly modeled from composites of population-level risk factors. These models overlook individual context in communities' experiences of environmental hazards, producing metrics that may hinder spatial decision support for mitigating and responding to hazards. This paper introduces an interpretable, high-resolution model for generating an individual-oriented social vulnerability index (IOSVI) for the United States built on synthetic populations that couples individual and social determinants of vulnerability. The IOSVI combines an individual vulnerability index (IVI) that ranks individuals in an area’s synthetic population based on intersecting risk factors, with a social vulnerability index (SVI) based on the population’s cumulative distribution of IVI scores. Interpretability of the IOSVI procedure is demonstrated through examples of national, metropolitan, and neighborhood (census tract) level spatial variation in index scores and IVI themes, as well as an exploratory analysis examining risk factors affecting a specific sub-population (military veterans) in areas of high social and environmental vulnerability.

Tuccillo, Joe↗
NASA NTRS2010

Abundances of Jupiter's Trace Hydrocarbons from Voyager and Cassini

The flybys of Jupiter by the Voyager spacecraft in 1979, and over two decades later by Cassini in 2000, have provided us with unique datasets from two different epochs, allowing the investigation of seasonal change in the atmosphere. In this paper we model zonal averages of thermal infrared spectra from the two instruments, Voyager 1 IRIS and Cassini CIRS, to retrieve the vertical and meridional profiles of temperature, and the abundances of the two minor hydrocarbons, acetylene (C2H2) and ethane (C2H6). The spatial variation of these gases is controlled by both chemistry and dynamics, and therefore their observed distribution gives us an insight into both processes, We find that the two gases paint quite different pictures of seasonal change. Whilst the 2-D cross-section of C2H6 abundance is slightly increased and more symmetric in 2000 (northern summer solstice) compared to 1979 (northern fall equinox), the major trend of equator to pole increase remains. For C2H2 on tile other hand, the Voyager epoch exhibits almost no latitudinal variation, whilst the Cassini era shows a marked decrease polewards in both hemispheres. At the present time, these experimental findings are in advance of interpretation, as there are no published models of 2-D Jovian seasonal chemical variation available for comparison.

Nixon, C. A.↗
DOE OSTI2023

Development of fracture diagnostic methods for fluid distribution based on quantitative interpretation of distributed acoustic sensing and distributed temperature sensing

Multistage hydraulic fracturing design on horizontal wells has significantly evolved with larger fluid volume, more fracturing stages, and tighter perforation cluster spacing to efficiently stimulate unconventional reservoirs. From the published field observations, the recent fracturing design results in complex fracture networks or swarm of fractures. Fracture treatment evaluation is extremely challenging in such a case because of the large amount of variables in well completion and stimulation design. Combined measurements from different technologies can help in fracture diagnosis. Fluid distribution, either during fracture injection or during production, directly relates to the stimulation efficiency at the cluster level and at the stage level. Because it is unlikely in the real world to distribute the injected fluid uniformly among all the clusters, we need diagnostic techniques to generate the flow profile along a lateral. Fiber-optic measurements, such as distributed acoustic sensing (DAS) and distributed temperature sensing (DTS), are currently used to diagnose downhole flow conditions. This technology allows us to qualitatively confirm the fluid flow profile and other issues occurring downhole during fracturing such as leakage through plugs. For optimizing a fracturing design, we also need to understand how the design parameters are correlated with the stimulation efficiency. In this study, we combine two sets of models of DAS and DTS data interpretation for injected fluid volume distribution. The DAS is interpreted based on an empirical correlation between fluid flow rates and frequency band energy from the acoustic signals. The DTS is interpreted by performing temperature history match-based thermal energy conservation. Because of the completely different physics behind the interpretations, the confirmation of two interpretations provides confidence in fluid distribution.

Geochemistry & Geophysics↗
NASA NTRS1984

Aerosol backscatter lidar calibration and data interpretation

A treatment of the various factors involved in lidar data acquisition and analysis is presented. This treatment highlights sources of fundamental, systematic, modeling, and calibration errors that may affect the accurate interpretation and calibration of lidar aerosol backscatter data. The discussion primarily pertains to ground based, pulsed CO2 lidars that probe the troposphere and are calibrated using large, hard calibration targets. However, a large part of the analysis is relevant to other types of lidar systems such as lidars operating at other wavelengths; continuous wave (CW) lidars; lidars operating in other regions of the atmosphere; lidars measuring nonaerosol elastic or inelastic backscatter; airborne or Earth-orbiting lidar platforms; and lidars employing combinations of the above characteristics.

Kavaya, M. J.↗
DOE OSTI2023

A data-driven framework for permeability prediction of natural porous rocks via microstructural characterization and pore-scale simulation

Understanding the microstructure–property relationships of porous media is of great practical significance, based on which macroscopic physical properties can be directly derived from measurable microstructural informatics. However, establishing reliable microstructure–property mappings in an explicit manner is difficult, due to the intricacy, stochasticity, and heterogeneity of porous microstructures. In this paper, a data-driven computational framework is presented to investigate the inherent microstructure–permeability linkage for natural porous rocks, where multiple techniques are integrated together, including microscopy imaging, stochastic reconstruction, microstructural characterization, pore-scale simulation, feature selection, and data-driven modeling. A large number of 3D digital rocks with a wide porosity range are acquired from microscopy imaging and stochastic reconstruction techniques. A broad variety of morphological descriptors are used to quantitatively characterize pore microstructures from different perspectives, and they compose the raw feature pool for feature selection. Here high-fidelity lattice Boltzmann simulations are conducted to resolve fluid flow passing through porous media, from which reliable permeability references are obtained. The optimal feature set that best represents permeability is identified through a performance-oriented feature selection process, upon which a cost-effective surrogate model is rapidly fitted to approximate the microstructure-permeability mapping via data-driven modeling. This surrogate model exhibits great advantages over empirical/analytical formulas in terms of prediction accuracy and generalization capacity, which can predict reliable permeability values spanning four orders of magnitude. Besides, feature selection also greatly enhances the interpretability of the data-driven prediction model, from which new insights into the mechanism of how microstructural characteristics determine intrinsic permeability are obtained.

58 GEOSCIENCES↗
NASA NTRS1983

Estimation of clear-sky insolation using satellite and ground meteorological data

Ground based pyranometer measurements were combined with meteorological data from the Tiros N satellite in order to estimate clear-sky insolations at five U.S. sites for five weeks during the spring of 1979. The estimates were used to develop a semi-empirical model of clear-sky insolation for the interpretation of input data from the Tiros Operational Vertical Sounder (TOVS). Using only satellite data, the estimated standard errors in the model were about 2 percent. The introduction of ground based data reduced errors to around 1 percent. It is shown that although the errors in the model were reduced by only 1 percent, TOVS data products are still adequate for estimating clear-sky insolation.

Staylor, W. F.↗