Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic Modeling”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multiphysics modeling of a critical dual-purpose canister in a saturated geological repository

Dual purpose dry cask storage canisters for spent nuclear fuel are designed for storage and transportation, but are not licensed for permanent disposal in a geological repository. If dual purpose canisters were to be used to dispose of spent nuclear fuel in a geological repository, they would be expected to eventually breach and be flooded with groundwater, and it is shown that some fraction of these canisters will achieve criticality. In order to evaluate the consequences of canisters going critical in a repository, an initial capability has been developed for estimating the quasi-static power level of a critical canister using loosely coupled multiphysics simulations. The low power level in a critical canister enables coupling through precomputed physics proxies. This calculated power level is then used to compute the change in the critical canister’s isotopic inventory as a function of time. Three as-loaded canisters are evaluated and two were found to have power levels below 4 kW, with a modest effect on the radiological inventory over time. This effort also shows that although some DPCs will have extremely peaked power shapes, the relatively low power and long-time scales result in relatively homogenized thermohydraulic properties in the water within the DPC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1985

Geologic mapping using Landsat MSS and TM images - Removing vegetation by modeling spectral mixtures

Semi-arid regions of the island of Hawaii and of the Tucson Mountains, AZ, having less than 50 percent green vegetation cover were mapped into rock/soil types by first identifying spectral mixtures of vegetation and rock/soil, and then removing the vegetation information. Landsat MSS images were used for Hawaii, and a Landsat TM image was used for the Tucson Mountains. Following analysis of the images field transects were done in test areas to determine percent vegetation cover and vegetation community characteristics, and to observe rock and soil types. For both Hawaii and the Tucson Mountains 20 percent or more of green vegetation caused potential misidentification of the rock/soil units unless spectral mixing was taken into account.

Adams, J. B.↗
NASA NTRS1999

Strategies and Technologies for In Situ Mineralogical Investigations on Mars

Surface landers on Mars (Viking and Pathfinder) have not revealed satisfying answers to the mineralogy and lithology of the planet's surface. In part, this results from their prime directives: Viking focused on exobiology, Pathfinder focused on technology demonstration. The analytical instruments on board the landers made admirable attempts to extract the mineralogy and geology of Mars, as did countless modeling efforts after the missions. Here we suggest a framework for elucidating martian, or any other planetary geology, through an approach that defines (a) type of information required, (b) explorational strategy harmonious with acquisition of these data, (c) interpretation approach to the data, (d) compatible mission architecture, (e) instrumentation for interrogating rocks and soil. (a) Data required: The composition of a planet is ordered at scales ranging from molecules to minerals to rocks, and from geological units to provinces to planetary-scale systems. The largest ordering that in situ compositional instruments can attempt to interrogate is rock type "aggregate" information. This is what the geologist attempts to identify first. From this, mineralogy can be either directly seen or inferred. From mineralogy can be determined elemental abundances and perhaps the state of the compounds as being crystalline or amorphous. Knowledge of rock type and mineralogy is critical for elucidating geologic process. Mars landers acquired extremely valuable elemental data, but attempted to move from elements to aggregates, but this can only be done by making many assumptions and sometimes giant leaps of faith. Data we believe essential are elements, minerals, degree of ordering of compounds, and the aggregate or rock type that these materials compose. (b) Explorational strategy: A lander should function as a surrogate geologist. Of the total landscape, a geologist sees much, but gives detailed attention to an infinitesimally small amount of what is seen. To acquire samples worth detailed scrutiny, as many samples as possible need examining at a cursory or reconnaissance level. A representative, statistically-meaningful sample number cannot be overemphasized. This maxim still applies to geological exploration of our own planet of which we have abundant knowledge. Analysis of many samples mandates low-power consumption per sample. (c) Data interpretation: No single instrument can analyze the full spectrum of the x-axis. An instrument is optimized for detecting certain material characteristics and must therefore affix itself to some point on the x-axis. Any conclusions drawn about data to the left or right of the instrument's position on this axis must necessarily be derived by inference. Hence, it seems logical to include on a mission, instruments that are not closely spaced in their x-axis-position, and if only two analytical methods are used, as shown, they should start at opposite ends of the axis and work towards the center. As examples, we depict a high-resolution camera to evaluate rock type ("aggregate" state) and mineralogy, and an x-ray diffractometer-fluorescence spectrometer (XRD-XRF) to determine elements, minerals, and the degree of order of materials. (d) Mission architecture: No instrument or suite of instruments can be relied upon to always give truly unequivocal analyses. The suite of instruments should therefore permit conclusions of one instrument to be checked against those of another through closed analytical loops. These "loops" can be structured by a combination of orbital imagery, descent imagery, broad-band site viewing/analysis, and data that cover both x and y axes. For example, the detection of a basaltic-looking rock with a microscope should be checked against the elements detected, the appearance of the rock as a lava flow from descent imagery, and so forth. (e) Instrumentation: To satisfy the above criteria, it is necessary to: (i) See the rock or soil with high resolution + magnification, (ii) Examine many samples, (iii) Consume little power per analysis, (iv) Determine elemental species, (v) Determine mineralogy directly (not inferentially) and the degree of ordering of compounds, (vi) Start analyzing from both ends of the x-axis. Every geologist wants to see the hand sample first, and apply a hand lens to its surface. This has not been the starting point for missions to Mars. Thus, our technology satisfies all these criteria . This XRD-XRF-Optical instrument currently being developed, analyses rock or soil surfaces without the need for sample acquisition or preparation; this satisfies the power criterion, and enables many analyses. The device acquires direct mineralogy and determines elemental species. The embedded endoscopic camera satisfies the critical criterion of close inspection of samples; the fiber optic cable can also be used for IR, LTV, or laser sample analysis. Additional information is contained in the original (Figures).

Marshall, J. R.↗
NASA NTRS1992

Modeling orbital changes on tectonic time scales

Geologic time series indicate significant 100 ka and 400 ka pre-Pleistocene climate fluctuations, prior to the time of such fluctuations in Pleistocene ice sheets. The origin of these fluctuations must therefore depend on phenomena other than the ice sheets. In a previous set of experiments, we tested the sensitivity of an energy balance model to orbital insolation forcing, specifically focusing on the filtering effect of the Earth's geography. We found that in equatorial areas, the twice-yearly passage of the sun across the equator interacts with the precession index to generate 100 ka and 400 ka power in our modeled time series. The effect is proportional to the magnitude of land in equatorial regions. We suggest that such changes may reflect monsoonal variations in the real climate system, and the subsequent wind and weathering changes may transfer some of this signal to the marine record. A comparison with observed fluctuations of Triassic lake levels is quite favorable. A number of problems remain to be studied or clarified: (1) the EBM experiments need to be followed up by a limited number of GCM experiments; (2) the sensitivity to secular changes in orbital forcing needs to be examined; (3) the possible modifying role of sedimentary processes on geologic time series warrants considerably more study; (4) the effect of tectonic changes on Earth's rotation rate needs to be studied; and (5) astronomers need to make explicit which of their predictions are robust and geologists and astronomers have to agree on which of the predictions are most testable in the geologic record.

Crowley, Thomas J.↗
DOE OSTI2020

Sensor Recommendations for Long Term Monitoring of the F-Area Seepage Basins

In mid-2018, a new paradigm for long-term monitoring was developed after of decade of applied research projects funded by the Department of Energy’s office of Environmental Management Technology Development program. The program at SRNL was focused on transitioning complex environmental waste sites from active to passive remediations strategies. A key result of these studies was that the use of enhanced attenuation approaches at radiologically contaminated sites will result in the creation of secondary source areas in the subsurface that will require monitoring for decades. Alternative monitoring approaches are being developed and tested at the Savannah River Site’s F-Area Hazardous Waste Management Facility, the new paradigm provides innovative solutions that will significantly lower costs of monitoring through the coupling of data collection, machine learning and deterministic groundwater modeling. The foundation of this approach is a well-optimized network of sensors for measuring hydrogeochemical master variables that control, and therefore act as indicators of groundwater contaminant transport. By monitoring changes in the controlling master variables over time arising from geological and environmental shifts, predictive modelling can assist with identifying new strategies for ensuring regulatory requirements are met if trends toward conditions for potential remobilization of attenuated contaminants are detected. In this report, we evaluated commercially available single parameter sensor platforms (e.g., temperature/depth) and configurable multi-parameter sensor platforms (e.g., pH, oxidation-reduction potential, temperature, depth, dissolved oxygen, and conductivity). Each was scored using an optimization function based on how well the system supports the proposed long-term monitoring paradigm, in general, and the site-specific conditions at F-Area, in particular. Several viable sensor systems were identified. Of these, a combined platform including the In-Situ Aqua TROLL 500 multi-parameter sensor platform and the In-Situ temperature/depth sensor had the highest rating and was identified as the most suitable candidate for installation and monitoring of the master variables and potentiometric surface that control groundwater contaminant plumes emanating from the F-Area Seepage Basins. The discussion of recommended potential deployment locations builds upon recommendations made by Denham et al (2019).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1990

Remote sensing science for the Nineties; Proceedings of IGARSS '90 - 10th Annual International Geoscience and Remote Sensing Symposium, University of Maryland, College Park, May 20-24, 1990. Vols. 1, 2, & 3

Various papers on remote sensing (RS) for the nineties are presented. The general topics addressed include: subsurface methods, radar scattering, oceanography, microwave models, atmospheric correction, passive microwave systems, RS in tropical forests, moderate resolution land analysis, SAR geometry and SNR improvement, image analysis, inversion and signal processing for geoscience, surface scattering, rain measurements, sensor calibration, wind measurements, terrestrial ecology, agriculture, geometric registration, subsurface sediment geology, radar modulation mechanisms, radar ocean scattering, SAR calibration, airborne radar systems, water vapor retrieval, forest ecosystem dynamics, land analysis, multisensor data fusion. Also considered are: geologic RS, RS sensor optical measurements, RS of snow, temperature retrieval, vegetation structure, global change, artificial intelligence, SAR processing techniques, geologic RS field experiment, stochastic modeling, topography and Digital Elevation model, SAR ocean waves, spaceborne lidar and optical, sea ice field measurements, millimeter waves, advanced spectroscopy, spatial analysis and data compression, SAR polarimetry techniques. Also discussed are: plant canopy modeling, optical RS techniques, optical and IR oceanography, soil moisture, sea ice back scattering, lightning cloud measurements, spatial textural analysis, SAR systems and techniques, active microwave sensing, lidar and optical, radar scatterometry, RS of estuaries, vegetation modeling, RS systems, EOS/SAR Alaska, applications for developing countries, SAR speckle and texture.

Source record↗
NASA NTRS2003

Numerical Modeling of Shatter Cones Development in Impact Craters

Shatter cones are the characteristic forms of rock fractures in impact structures. They have been used for decades as unequivocal fingerprints of meteoritic impacts on Earth. The abundant data about shapes, apical angles, sizes and distributions of shatter cones for many terrestrial impact structures should provide insights for the determination of impact conditions and characteristics of shock waves produced by high-velocity projectiles in geologic media. However, previously proposed models for the formation of shatter cones do not agree with observations. For example, the widely accepted Johnson-Talbot mechanism requires that the longitudinal stress drops to zero between the arrival of the elastic precursor and the main plastic wave. Unfortunately, observations do not support such a drop. A model has been also proposed to explain the striated features on the surface of shatter cones but can not invoked for their conical shape. The mechanism by which shatter cones form thus remains enigmatic to date. In this paper we present a new model for the formation of shatter cones. Our model has been tested by means of numerical simulations using the hydrocodes SALE 2D enhanced with the Grady-Kipp-Melosh fragmentation model.

Baratoux, D.↗
DOE OSTI2022

Forecasting induced seismicity in Oklahoma using machine learning methods

Oklahoma earthquakes in the past decade have been mostly associated with wastewater injection. Here we use a machine learning technique—the Random Forest to forecast induced seismicity rate in Oklahoma based on injection-related parameters. We split the data into training (2011.01–2015.05) and test (2015.06–2020.12) periods. The model forecasts seismicity rate during the test period based on input features, including operational parameters (injection rate and pressure), geological information (depth to basement), and modeled pore pressure and poroelastic stress. The results show overall good match with observed seismicity rate (adjusted R 2 of 0.75). The model shows that pore pressure rate and poroelastic stressing rates are the two most important features in forecasting. The absolute values of pore pressure and poroelastic stress, and the injection rate itself, are less important than the stressing rates. These findings further emphasize that temporal changes of stressing rates would lead to significant changes in seismicity rates.

58 GEOSCIENCES↗
NASA NTRS2021

Post-Landing Major Element Quantification Using SuperCam Laser Induced Breakdown Spectroscopy

The SuperCam instrument on the PerseveranceMars 2020 rover uses a pulsed 1064 nm laser to ablate targets at a distance and conduct laser induced breakdown spectroscopy (LIBS) by analyzing the light from the resulting plasma. SuperCam LIBS spectra are preprocessed to remove ambient light, noise, and the continuum signal present in LIBS observations. Prior to quantification, spectra are masked to remove noisier spectrometer regions andspectra are normalized to minimize signal fluctuations and effectsof target distance.In some cases, the spectra are also standardized or binned prior to quantification. To determine quantitative elemental compositionsof diverse geologic materials at Jezero crater, Mars, we use a suite of 1198 laboratory spectra of 334 well-characterized reference samples. The samples were selected to span a wide range of compositions and include typical silicate rocks, pure minerals (e.g.,silicates, sulfates, carbonates, oxides),more unusual compositions (e.g.,Mn oreand sodalite), andreplicates of the sintered SuperCam calibration targets (SCCTs) onboardthe rover. For each major element (SiO2, TiO2, Al2O3, FeOT, MgO, CaO, Na2O, K2O), the database was subdivided into five“folds” with similar distributions of the element of interest. One fold was held out as an independent test set, and the remaining fourfolds were used to optimize multivariate regression models relating the spectrum to the composition. We considered a variety of models, and selected several for further investigation for each element, based primarily on the root mean squared error of prediction (RMSEP) on the test set, when analyzed at 3m. In cases with several models of comparable performance at 3 m, we incorporated the SCCT performance at different distances to choose the preferred model. Shortly after landing on Mars and collecting initial spectra of geologic targets, we selected one model per element. Subsequently, with additional data from geologic targets, some models were revised to ensure results that are more consistent with geochemical constraints. The calibration discussed here is a snapshot of an ongoing effort to deliver the most accurate chemical compositions with SuperCam LIBS.

Mars 2020↗
NASA NTRS1983

An investigation of MAGSAT and complementary data emphasizing precambrian shields and adjacent areas of West Africa and South America

There is lack of acceptance of quantitative crustal models by specialists on regional geology because of the problem of east-west banding on MAGSAT anomaly maps. The possibility of some of the Earth's magnetic field being created in the mantle below the Curie isotherm is discussed as well as the need to remove spherical harmonics of degrees less than 15 from the anomaly map. Two dimensional models were produced for northern South America that might be considered as working models. These models suggest that the presumed Amazon River depression is better developed near Belem and in the upper Amazon basin than in the central Amazon basin. An attempt is proposed for separating the core field from the crustal field using a flat-Earth model.

Hastings, D. A.↗
DOE OSTI2023

Time-lapse VSP integration and calibration of subsurface stress field utilizing machine learning approaches: A case study of the morrow B formation, FWU

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse vertical seismic profiling (VSP) integration. The selected study site is a region around the injector well located within Farnsworth field unit (FWU), where there is an ongoing CO 2 -enhanced oil recovery (EOR) operation. In our study, a site-specific rock physics model was created from extensive geological, geophysical, and geomechanical characterization through 3D seismic data, well logs, and core assessed as part of the 1D MEM conducted on the characterization well within the study area. Here, the Biot-Gassmann workflow was utilized to combine the rock physics and reservoir simulation outputs to determine the seismic velocity change due to fluid substitution. Modeled seismic velocities attributed to mean effective stress were determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from ultrasonic seismic velocity measurements. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer (PSO) was utilized to minimize a penalty function created between the modeled seismic velocities and the observed time-lapse VSP dataset. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.

58 GEOSCIENCES↗
DOE OSTI2021

Post-landing major element quantification using SuperCam laser induced breakdown spectroscopy

The SuperCam instrument on the Perseverance Mars 2020 rover uses a pulsed 1064 nm laser to ablate targets at a distance and conduct laser induced breakdown spectroscopy (LIBS) by analyzing the light from the resulting plasma. SuperCam LIBS spectra are preprocessed to remove ambient light, noise, and the continuum signal present in LIBS observations. Prior to quantification, spectra are masked to remove noisier spectrometer regions and spectra are normalized to minimize signal fluctuations and effects of target distance. In some cases, the spectra are also standardized or binned prior to quantification. To determine quantitative elemental compositions of diverse geologic materials at Jezero crater, Mars, we use a suite of 1198 laboratory spectra of 334 well-characterized reference samples. The samples were selected to span a wide range of compositions and include typical silicate rocks, pure minerals (e.g., silicates, sulfates, carbonates, oxides), more unusual compositions (e.g., Mn ore and sodalite), and replicates of the sintered SuperCam calibration targets (SCCTs) onboard the rover. For each major element (SiO 2 , TiO 2 , Al 2 O 3 , FeO T , MgO, CaO, Na 2 O, K 2 O), the database was subdivided into five “folds” with similar distributions of the element of interest. One fold was held out as an independent test set, and the remaining four folds were used to optimize multivariate regression models relating the spectrum to the composition. We considered a variety of models, and selected several for further investigation for each element, based primarily on the root mean squared error of prediction (RMSEP) on the test set, when analyzed at 3 m. In cases with several models of comparable performance at 3 m, we incorporated the SCCT performance at different distances to choose the preferred model. Shortly after landing on Mars and collecting initial spectra of geologic targets, we selected one model per element. Subsequently, with additional data from geologic targets, some models were revised to ensure results that are more consistent with geochemical constraints. The calibration discussed here is a snapshot of an ongoing effort to deliver the most accurate chemical compositions with SuperCam LIBS.

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

UNESE Data Analysis-Disko Elm Gas Transport Characteristics (NA-22 Quarterly Report)

Recently published papers have demonstrated the impact of sorption mechanisms on gas transport in the subsurface following an underground nuclear explosion. To include sorption processes in the multi-physics models of noble gas production and transport, researchers have conducted experiments of gas sorption in geologic materials and calibrated a linear (Henry’s law) sorption model for all gases in various geologic materials. As a part of Task #4, we requested the experimental data from Sandia Laboratories and derived a dual-mode (Henry’s absorption and Langmuir adsorption) sorption model. The derived model better describes the nonlinear sorption mechanism with much higher fitness with the experimental data. When gas concentrations are relatively low (in the conditions of field experiments or underground nuclear explosions) compared to those used sorption experiments, the nonlinear sorption model can be approximated as a linear sorption model. The equivalent Henry’s constant, which equals the product of Langmuir capacity and affinity, is systematically higher than the literature values. Further analyses and comparison of sorption models are planned in FY2020 Q4 and the sorption model will be coupled in models of nonisothermal multiphase transport for studying noble gas detectability.

58 GEOSCIENCES↗
NASA NTRS2013

Challenges in Quantifying Pliocene Terrestrial Warming Revealed by Data-Model Discord

Comparing simulations of key warm periods in Earth history with contemporaneous geological proxy data is a useful approach for evaluating the ability of climate models to simulate warm, high-CO2 climates that are unprecedented in the more recent past. Here we use a global data set of confidence-assessed, proxy-based temperature estimates and biome reconstructions to assess the ability of eight models to simulate warm terrestrial climates of the Pliocene epoch. The Late Pliocene, 3.6-2.6 million years ago, is an accessible geological interval to understand climate processes of a warmer world4. We show that model-predicted surface air temperatures reveal a substantial cold bias in the Northern Hemisphere. Particularly strong data-model mismatches in mean annual temperatures (up to 18 C) exist in northern Russia. Our model sensitivity tests identify insufficient temporal constraints hampering the accurate configuration of model boundary conditions as an important factor impacting on data- model discrepancies. We conclude that to allow a more robust evaluation of the ability of present climate models to predict warm climates, future Pliocene data-model comparison studies should focus on orbitally defined time slices.

Salzmann, Ulrich↗
NASA NTRS1987

Wind ripples in low density atmospheres

The effect of varying fluid density (rho) on particle transport was examined by conducting tests at atmospheric pressures between 1 and 0.004 bar in the Martian Surface Wind Tunnel (MARSWIT). This study specifically concerns the effect of varying rho on the character of wind ripples, and elicits information concerning generalized ripple models as well as specific geological circumstances for ripple formation such as those prevailing on Mars. Tests were conducted primarily with 95 micron quartz sand, and for each atmospheric pressure chosen, tests were conducted at two freestream wind speeds: 1.1 U*(t) and 1.5 U*(t), where U*(t) is saltation threshold. Preliminary analysis of the data suggests: (1) ballistic ripple wavelength is not at variance with model predictions; (2) an atmospheric pressure of approximately 0.2 bar could represent a discontinuity in ripple behavior; and (4) ripple formation on Mars may not be readily predicted by extrapolation of terrestrial observations.

Miller, J. S.↗
NASA NTRS2004

Numerical Modeling of Glaciers in Martian Paleoclimates

Numerous geologic features suggest the presence of ice flow on the surface of mars. These features include lobate debris aprons, concentric crater fill, and lineated valley fill. The lateral extent of these features can range from 100 meters to over 20 km. Previous work has demonstrated that these features could not have formed in current Martian conditions. It has long been speculated that changes in Mars orbital properties, namely its obliquity, eccentricity, and argument of perihelion, can result in dramatic changes to climate. Recent climate model studies have shown that at periods of increased obliquity north polar water ice is mobilized southward and deposited at low ad mid latitudes. Mid latitude accumulation of ice would provide the necessary conditions for rock glaciers to form. A time-marching, finite element glacier model is used to demonstrate the ability of ice and ice-rock mixtures to flow under Martian paleoclimate conditions. Input to this model is constrained by the NASA Ames Mars General Circulation Model (MGCM).

Colaprete, A.↗
DOE OSTI2020

Multi-Model and Multi-Scale Global Sensitivity Analysis for Identifying Controlling Processes of Complex Systems

An environmental model consists of multiple process level sub-models, and each sub-model represents a process that is key to the operation of the simulated system. Global sensitivity analysis methods have been widely used to identify important processes for system model development and improvement. The existing methods of global sensitivity analysis only consider parametric uncertainty, and are not capable of handling model uncertainty caused by multiple process models that arise from competing hypotheses about one or more processes. To address this problem, this project develops a new method to probe model output sensitivity to competing process models by integrating model averaging methods with variance-based global sensitivity analysis to address uncertainty in process models and parameters. The new method yields three process sensitivity indices. The first one is called first-order process sensitivity index, and it is derived as a single summary measure of relative process importance. Evaluating the index is computationally expensive, because it relies in a Monte Carlo scheme that requires thousands and even millions of model executions. To reduce computational cost, this project develops a computationally efficient, quasi Monte Carlo method, and this method is presented in Chapter 2 of this report with and a numerical example for demonstration. The numerical example shows that the results of the quasi Monte Carlo method are substantially close to those of the full Monte Carlo method, but the computational cost of the quasi Monte Carlo method is only 0.7% of that of the full Monte Carlo method. The second index is called total-effect process sensitivity index, and it measures interactions between different processes. Therefore, this sensitivity index includes the first-order process sensitivity index, and can be used to identify influential processes. On the other hand, the total-effect process sensitivity index can also be used to screen non-influential processes. This is demonstrated by two numerical examples using the Sobol-G* functions and groundwater flow models that consider recharge process, geological process, and snowmelt process. The numerical examples shows that the total-effect process sensitivity index is more informative than the first-order process sensitivity. The derivation of the process sensitivity index and the numerical examples are discussed in Chapter 3. Chapter 4 presents two computationally efficient methods for screening non-influential processes to exclude them from further investigation. The two methods are the multi-model difference-based sensitivity (MMDS) analysis method, which can be implemented using the Latin Hypercube Sampling. The second one is the implementation of MMDS method using a binning method. The numerical example for the Sobol-G* function indicates the two methods are capable of identifying non-influential models, and the numerical examples for the groundwater flow and reactive transport show that the two methods are effective for groundwater problems. However, it should be noted that the two methods are numerical approximations, and they can only be used for screening non-influential processes, not for ranking importance of system processes. All the sensitivity analysis methods are implemented by developing python codes, and the codes are in a software called SAMMPY: a python package for process sensitivity analysis under multiple models. The SAMMPY design and structure are discussed in Chapter 5, and the package is released to the public for free download.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Physics-informed semantic inpainting: Application to geostatistical modeling

A fundamental problem in geostatistical modeling is to infer the heterogeneous geological field based on limited measurements and some prior spatial statistics. Semantic inpainting, a technique for image processing using deep generative models, has been recently applied for this purpose, demonstrating its effectiveness in dealing with complex spatial patterns. However, the original semantic inpainting framework incorporates only information from direct measurements, while in geostatistics indirect measurements are often plentiful. In this work, to overcome this limitation, we propose a physics-informed semantic inpainting framework, employing the Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) and jointly incorporating the direct and indirect measurements by exploiting the underlying physical laws. Our simulation results for a high-dimensional problem with 512 dimensions show that in the new method, the physical conservation laws are satisfied and contribute in enhancing the inpainting performance compared to using only the direct measurements.

54 ENVIRONMENTAL SCIENCES↗