SEARCH · Search NASA

Results for “VERTICAL MIGRATION”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Plastic shrinkage of mortars cured with a paraffin-based compound – Bimodal neutron/X-ray tomography study

Early-age drying (immediately after casting) of mortars and the corresponding plastic shrinkage were studied using bimodal neutron/X-ray computed tomography. This novel, correlative 3D imaging mode enabled studying simultaneously and without any source of spurious perturbation the water migration and loss processes together with the corresponding deformations due to plastic shrinkage. Bimodal imaging opens up new possibilities for studying dynamic processes of coupled water transport and deformations in porous solids. The measurements were carried out on model systems (cylindrical mortar specimens with height of 19 mm). The study focused on the effect of a paraffin-based curing compound. Our results confirm that when the curing compound was applied directly onto the drying surface in a sufficient amount, both the evaporation rate and the rate of vertical displacement (settlement) were substantially reduced. The results shed a new light on the mechanisms of plastic shrinkage and the action of curing compounds.

36 MATERIALS SCIENCE↗

NASA NTRS1989

Thermocapillary migration of liquid droplets in a temperature gradient in a density matched system

An experimental investigation of thermocapillary flow in droplets of a vegetable oil (partially hydrogenated soybean oil) immersed in silicone oil was conducted in a test cell with a heated top wall and a cooled bottom wall. The liquids are nearly immiscible and have equal densities at a temperature below the room temperature, thus providing a simulation of low-gravity conditions by reducing the buoyancy forces. The interfacial tension between the two oils was measured in the temperature range 20 to 50 C using a capillary tube and (d sigma)/(d T) was determined to be negative. Droplets ranging in sizes from 3 mm to 1 cm diameter were injected into the silicone oil. The vertical temperature profile in the bulk liquid (silicone oil) produces temperature variations along the interface which induce variations in the interfacial tension. The flow inside the droplet driven by the resulting interfacial shear stresses was observed using a laser light-sheet flow visualization technique. The flow direction is consistent with the sign of (d sigma)/(d T). The observed maximum surface velocities are compared to the theoretical predictions of Young et al. (1959).

Rashidnia, N.↗

NASA NTRS1991

Thermocapillary migration of liquid droplets in a temperature gradient in a density matched system

Rashidnia, N.↗

NASA NTRS1991

Factors controlling the structures of magma chambers in basaltic volcanoes

The depths, vertical extents, and lateral extents of magma chambers and their formation are discussed. The depth to the center of a magma chamber is most probably determined by the density structure of the lithosphere; this process is explained. It is commonly assumed that magma chambers grow until the stress on the roof, floor, and side-wall boundaries exceed the strength of the wall rocks. Attempts to grow further lead to dike propagation events which reduce the stresses below the critical values of rock failure. The tensile or compressive failure of the walls is discussed with respect to magma migration. The later growth of magma chambers is accomplished by lateral dike injection into the country rocks. The factors controlling the patterns of growth and cooling of such dikes are briefly mentioned.

Wilson, L.↗

NASA NTRS2009

Five-Year (2004-2009)Observations of Upper Tropospheric Water Vapor and Cloud Ice from MLS and Comparisons with GEOS-5 Analyses

This paper gives an overview of August 2004 through July 2009 upper tropospheric (UT) water vapor (H2O) and ice water content (IWC) from the Aura Microwave Limb Sounder (MLS) and comparisons with outputs from the NASA Goddard Earth Observing System Version 5 (GEOS-5) data assimilation system. Both MLS and GEOS-5 show that high values of H2O and IWC at 215 to 147 hPa are associated with areas of deep convection. They exhibit good (within approximately 15%) agreement in IWC at these altitudes, but GEOS-5 H2O is approximately 50% (215 hPa) to approximately 30% (147 hPa) larger than MLS, possibility due to its higher temperatures at these altitudes. GOES-5 produces a weaker intertropical convergence zone than MLS, while a seasonally-migrating band of tropical deep convection is clearly evident in both the MLS and GEOS-5 UT H2O and IWC. MLS and GEOS-5 both show spatial anti-correlation between IWC and H2O at 100 hPa, where less H2O is associated with low temperatures in regions of tropical convection. At 100 hPa, GEOS-5 produces 50% less IWC and 15% less H2O in the tropics, and approximately 20% more H2O in the extra-tropics, than does MLS. Behavior of the 100 hPa H2O, which exhibits a quasi-biennial oscillation, appears consistent with it being controlled by temperature. The seasonal cycle in the vertical transport of tropical mean H2O from approximately 147 hPa to approximately 10 hPa appears much stronger in MLS than in GEOS-5. The UT IWC and H2O interannual variations, from both MLS and GEOS-5, show clear imprints of the El Nino-Southern Oscillation.

Jiang, Jonathan H.↗

NASA NTRS1971

Plate convergence, transcurrent faults and internal deformation adjacent to Southeast Asia and the western Pacific

A model for oblique convergence between plates of lithosphere is proposed in which at least a fraction of slip parallel to the plate margin results in transcurrent movements on a nearly vertical fault which is located on the continental side of a zone of plate consumption. In an extreme case of complete decoupling only the component of slip normal to the plate margin can be inferred from underthrusting. Recent movements in the western Sunda region provide the most convincing evidence for decoupling of slip, which in this region is thought to be oblique to the plate margin. A speculative model for convergence along the margins of the Philippine Sea is constructed from an inferred direction of oblique slip in the Philippine region. This model requires that the triple point formed by the junction of the Japanese and Izu-Bonin trenches and the Nankai trough migrate along the Sagami trough.

Fitch, T. J.↗

DOE OSTI2021

Performance and reliability of β-Ga 2 O 3 Schottky barrier diodes at high temperature

Beta-gallium oxide (ß-Ga 2 O 3 ) is an ultrawide bandgap semiconductor that has potential for power electronic applications and devices operating at high temperatures. Particularly important for these applications are its 4.9 eV bandgap, facile electron doping, and the ability to grow ß-Ga 2 O 3 crystals from the melt. In this work, vertical ß-Ga 2 O 3 Schottky barrier diodes were fabricated using Pt Schottky and Ti-based Ohmic contacts and Au contact pads on unintentionally doped n-type, (2¯01)-oriented single crystal substrates. The diode’s temperature-dependent electrical properties up to 400°C were investigated, and the Pt/Ga 2 O 3 Schottky barrier height was determined to be close to 1.2 eV. The degradation of the contacts over multiple cycles up to 400°C was observed, resulting in a significant increase in series resistance of the diodes by 1000× at ambient temperature after they were cycled. According to electron microscopy measurements, this degradation is likely due in part to the migration and oxidation of Ti at the top surface of the Au contact pads. This degradation highlights the need for further research and development to ensure stable Ohmic and Schottky contacts to Ga 2 O 3 at temperatures above 400°C.

36 MATERIALS SCIENCE↗

DOE OSTI2023

New Insights from Legacy Seismic Data regarding Basalt Elevations and Variability on the Hanford Site

Migration of groundwater contaminants in the Gable Gap area of the Hanford Site in southeastern Washington State is strongly influenced by the distribution and permeability of basalts that lie beneath an unconfined aquifer. Locally, folding and faulting of the Columbia River Basalt associated with the Yakima fold and thrust belt followed by erosion due to the Lake Missoula floods resulted in a complex basalt surface that represents either an impermeable lower boundary to the unconfined aquifer system or localized regions of increased permeability that potentially promote communication between the unconfined aquifer system and deeper, confined aquifer systems. Paleo-channels carved into the basalt by floodwaters are thought to provide preferential flow paths for groundwater contaminants. In 2011, a seismic landstreamer campaign was carried out to image the basalt surface and produced pre-stack depth migrated p-wave reflection images. The reflection images identified two large troughs that may represent paleo-channels and several areas of possible faulting. Here, the streamer data are re-analyzed using refraction travel-time and Rayleigh wave dispersion analyses to obtain images of compressional and shear wave velocities within the suprabasalt sediment sections and the upper basalt surface. The combined interpretation of reflection and seismic velocity images shows complexity in the basalt velocity and elevation, which varies by 50 m or more within the study area. These results, along with other ongoing geophysical investigations, will be used to inform the site geologic model and potentially guide placement of future boreholes needed to quantify vertical flow between the confined and unconfined aquifers.

58 GEOSCIENCES↗

NASA NTRS2024

Alaska Ecological Conservation: Using NASA Earth Observations to Identify Recent Changes in River Ice Phenology and Its Impacts on Caribou Migration

Each fall, caribou (Rangifer tarandus) in the Western Arctic Herd migrate hundreds of kilometers across northwestern Alaska to acquire seasonal resources and reach wintering grounds. Various rivers intersect migration paths, and caribou can only cross over open water or high-ice extents but are unable to cross rivers in stages of partial freezing. Recent temperature increases in Alaska can alter the timing and duration of ice formation periods, impeding migratory patterns. The Massachusetts NASA DEVELOP team partnered with the National Park Service in Alaska to detect river ice onset dates and formation periods using Landsat 8 Operational Land Imager and Landsat 9 Operational Land Imager-2, Sentinel-1 C-band Synthetic Aperture Radar and Sentinel-2A/B MultiSpectral Instrument imagery in Google Earth Engine. This feasibility analysis measured ice coverage using the Normalized Difference Infrared Index, Relative Difference River Ice, and Vertical-Vertical/Vertical-Horizontal backscatter values based on the spectral and surface characteristics of rivers. The team produced annual freezing timelines, time series plots, and maps at three river stretches to analyze river ice phenology changes. Radar imagery assessed ice coverage more accurately, while optical imagery better identified ice onset dates. Meanwhile indices were unable to robustly establish ice formation thresholds across the study period. These study results can help the National Park Service to better evaluate spatiotemporal migratory shifts and contextualize recent regional caribou declines.

river phenology↗

DOE OSTI2020

Mechanistic and predictive understanding of needle litter decay in semi-arid mountain ecosystems experiencing unprecedented vegetation mortality. Final report

This research compared needle litter decay rates and underlying mechanisms between litter released from insect-killed trees to that of naturally senesced litter from healthy trees across climate variables with elevation and snowmelt rates as variables. Our research was designed to address questions of whether there is a significant impact to C and N cycling at the terrestrial-atmospheric interface in association with pine needle decomposition and the role of changes in soil moisture and temperature as a function of these processes. Field activities were primarily conducted within the Rocky Mountain Biological Laboratory’s (RMBL) East River watershed across a spectrum of ecoregions that spanned approximately 700 vertical meters. Forested catchments at these elevations in southwestern Colorado are threatened by heat-induced drought, altered snowmelt patterns, and invasive insect infestation that have all been associated with climate change. Needle litter from two members of the Pinaceae family experiencing large-scale insect infestations (spruce and lodgepole pine) were collected, redeployed and then sampled over the course of three years. We experimentally tested three tractable parameters that are known to impact litter decay and biogeochemical cycling and are relevant to bark beetle-induced tree mortality: (1) moisture content, (2) temperature by using elevation as a surrogate, and (3) changes to litter compositional chemistry that result from beetle infestation. We examined decay processes and monitored the migration of soluble and gaseous nitrogen and carbon species as well as microbiological signatures of decay such as biodegradable fractions, metal availability and the temporal succession of associated microorganisms at the terrestrial interface.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

STILGAR End-of-Project Report

The Subsurface Tunnel Imaging LeveraGed by Analysis of Rayleigh wave ellipticity (STILGAR) project demonstrated an integrated geophysical approach for detecting, locating, and characterizing underground structural changes using dense seismic arrays and advanced inversion techniques. Field campaigns were conducted at two operational mines—the Redmond salt mine (Utah) and Graymont Pleasant Gap limestone mine (Pennsylvania)—providing real-world testbeds for monitoring anthropogenic subsurface activity. At the Redmond salt mine, seismic interferometry combined with back-projection inversion successfully identified continuous, low-amplitude signals from mining operations. The approach differentiated stationary from migrating anthropogenic sources, captured daily operational cycles, and validated the potential of passive seismic monitoring for remote detection of underground activity. At the Graymont Pleasant Gap mine, two dense seismic deployments in the spring and fall of 2023 generated over 4 TB of high-resolution data. Key outcomes included the relocation of 199 underground and 8 surface explosions with accuracies within tens of meters and the development of a 3D P-wave velocity model using the triple-difference tomography algorithm (tomoTD) that resolved major structural features such as the mine entrance, low-velocity tunnels, and roof-collapse areas. Ambient noise cross-correlation and back-projection analyses revealed persistent sources linked to ongoing mining activity, whereas horizontal-to-vertical spectral ratio (HVSR) and ellipticity studies confirmed stable site responses across seasons and identified soil thickness trends consistent with regional erosional and depositional processes. Checkerboard and sensitivity tests further validated the robustness of the tomographic results. Overall, the findings emphasize that although significant progress has been made in subsurface imaging, further work is needed to enhance the detection and localization of underground structures. Accurate imaging requires higher frequencies, yet anthropogenic sources tend to dominate the seismic record at those frequencies, and high-frequency surface waves are affected by higher modes that complicate interpretation. The improved detection and localization of human-induced signals enabled detailed temporal and spatial mapping of daily mine operations, demonstrating the feasibility of continuous anthropogenic source monitoring. Sensitivity to signals from nontraditional sources, such as fan operations, highlights the broader applicability of this approach to other industrial environments in which continuous and impulsive signals are present. The field campaigns produced a substantial volume of high-quality seismic data, supporting the development and testing of new methods for seismic source characterization and subsurface imaging. Future deployments should include sensors capable of recording lower frequencies to probe deeper structures, increase bandwidth to enhance resolution and sensitivity to both shallow and deep targets, and collect additional large-scale datasets to refine imaging and source characterization techniques. Moreover, conducting 3D modeling studies of seismic wavefields at higher frequencies will provide a better understanding of wave scattering and cavity–wavefield interactions in complex underground environments. In conclusion, the STILGAR project demonstrated that integrated seismic monitoring can effectively characterize underground operations, capturing both natural and anthropogenic signals. The approaches developed provide a foundation for improved detection, localization, and imaging of subsurface structures and are directly transferable to broader industrial monitoring applications.

58 GEOSCIENCES↗

NASA NTRS1988

Confined turbulent fluid-particle flow modeling using multiple-realization particle trajectory schemes

A multiple-realization particle trajectory scheme has been developed and applied to the numerical prediction of confined turbulent fluid-particle flows. The example flows investigated include the vertical pipe upflow experimental data of Tsuji et al. and the experimental data of Leavitt for a coaxial jet flow, comprising a particle-laden central jet and a clean annular jet, into a large recirculation chamber. The results obtained from the numerical scheme agree well with the experimental data, lending confidence to the modeling approach. The multiple-realization particle trajectory turbulent flow modeling scheme is believed to be a more elegant and accurate approach to the extension of single-particle hydrodynamics to dilute multi-particle systems than the more commonly employed two-fluid modeling approach. It is also better able to incorporate additional force items such as lift, virtual mass and Bassett history terms directly into the particle equation of motion as appropriate. This makes it a suitable candidate for particle migration studies and an extension to situations involving liquid particulate phases with possible propulsion applications, such as in spray combustion, follows naturally.

Adeniji-Fashola, A. A.↗

NASA NTRS2018

Quantifying Methane Leak Emissions by Fused Airborne Imaging Spectroscopy with in Situ Surface Mobile and Airborne Observations of a California Producing Oil Field

Methane (CH4) emission budgets remain uncertain and are projected to grow as oil and gas production from short-lived wells increases and their subsequent transport through aging gas distribution networks. Orders-of-magnitude variations in temporal, spatial, and emission scales present a key challenge to leak detection and quantification. Also, the probability distributions for large and stochastic, leaky systems such as geological reservoirs (by natural migration-seeps) and petroleum production from those reservoirs remain largely unknown, needed to address current approach limitations. The scale of many petroleum systems favors remote sensing, but the sensitivity of such systems often precludes detection of weak emissions. Consequently, an accurate evaluation requires that the relative contribution from the emission "tails" of small leaks also be quantified, which is best carried out using high-sensitivity in situ methods. Fusion of remote sensing and in situ approaches leverages complementary capabilities to address these limitations. We show results from mobile surface (AMOG) and airborne in situ (AJAX) and thermal-infrared (TIR) hyperspectral imaging spectroscopy (Mako) data applied to a producing oil field in the California Central Valley near Bakersfield. AMOG is an automobile-based mobile lab that measures 13 trace gases, aerosol size distributions and vertical profiles, 3D winds and other meteorology, and atmospheric column measurements by solar spectroscopy at highway speeds. AJAX measures 5 trace gases and 3D winds at ~140 m/s. Mako is a broad-area TIR imaging spectrometer that can discriminate multiple gases present in each pixel acquired. In situ-derived, total field emissions were 31±16 Gg/yr CH4. This was compared with Mako-derived emissions from all plumes identified across the study site. We found that super-emitters were not the dominant emissions mode and the spatial pattern of plume locations from production infrastructure was correlated to geological structures.

Leifer, Ira↗

NASA NTRS2002

Diurnal Cycle of Convection in the East Pacific ITCZ during EPIC-2001

During the last three weeks of September 2001, the EPIC-2001 intensive field campaign focused on studies of deep convection in the ITCZ over the Mexican warm pool region (10N, 95W) of the East Pacific. This study focuses on the pronounced observed diurnal cycle of environmental and convective parameters within the experiment domain. Data from three primary sources are examined: the R/V Ronald H. Brown C-band weather radar, 4-hourly soundings from the Brown and the Global Atmospherics, Inc. National Lightning Detection Network (long range product). Satellite data from TRMM, GOES and OV-1 are also used. The domain boundary layer shows a robust daily evolution of moist enthalpy (as reflect by equivalent potential temperature, theta-e, or wet bulb potential temperature, theta-w), with contributions from changes in both dry and moist entropy. Peak theta-w is found after local nightfall; the average diurnal range of theta-w is approximately 1 deg C. A composite diurnal cycle of convective properties was derived from the C-band volume scans, sampled continuously through the experiment at 10 minute updates. Products derived from the volumetric data include a surface PPI, 15 and 30 dBZ echo top height, vertically integrated liquid, and 6 km (mixed phase region) reflectivity CAPPIs. For almost all products, the parameter means showed virtually no diurnal cycle. However, for the upper-level products, the parameter spectra showed a clear peak in the occurrence of deep/vigorous convection (the "tail end of the distribution") between 7-9 UTC (1-3 AM local), while overall frequency of occurrence peaked later, from 12-15 UTC (6-9 AM local). This represents a daily "outbreak" of isolated deep cells a couple of hours after sunset and subsequent growth, organization and decay through the nighttime hours. The coherence of the diurnal cycle of the convective spectrum is impressive given the wide variety of convective organization observed during the experiment, and given the modulation by passage of 3-5 day easterly waves. While earlier satellite OLR composites suggested an offshore coastal migration of storms into the domain at night, examination of the 150 km and 300 km range radar products showed little evidence of such organization; almost all convection developed "in-place" within the analysis domain. Consistent with the diurnal thermodynamic and microphysical evolution, a clear cycle in cloud-to-ground (CG) lightning occurrence was observed. The local CG diurnal cycle is significantly stronger than the satellite-derived tropical ocean diurnal cycle of total (IC+CG) lightning. Flash rates of 3-4 fl/min were often visually observed after nightfall; these are fairly 'healthy' flash rates for tropical ocean storms, and the domain was qualitatively noted to be unusually lightning-productive by the R/V Brown crew (also consistent with satellite-based climatologies).

Boccippio, Dennis J.↗

NASA NTRS2007

Landsat-based Analysis of Mountain Forest-tundra Ecotone Response to Climate Trends in Sayan Mountains

observations of temperatures Siberia has shown a several degree warming over the past 30 years. It is expected that forest will respond to warming at high latitudes through increased tree growth and northward or upward slope migration. migration. Tree response to climate trends is most likely observable in the forest-tundra ecotone, where temperature mainly limits tree growth. Making repeated satellite observations over several decades provides an opportunity to track vegetation response to climate change. Based on Landsat data of the Sayan Mountains, Siberia, there was an increase in forest stand crown closure and an upward tree-line shift in the of the forest-tundra ecotone during the last quarter of the 2oth century,. On-ground observations, supporting these results, also showed regeneration of Siberian pine in the alpine tundra, and the transformation of prostrate Siberian pine and fir into arboreal (upright) forms. During this time period sparse stands transformed into closed stands, with existing closed stands increasing in area at a rate of approx. 1 %/yr, and advancing their upper border at a vertical rate of approx. 1.0 m/yr. In addition, the vertical rate of regeneration propagation is approx. 5 m/yr. It was also found that these changes correlated positively with temperature trends

Kharuk, Viatcheslav I.↗

DOE OSTI2022

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗

DOE OSTI2020

Time-Lapse Vertical Seismic Profiling (VSP) for CO 2 Storage in a Depleted Oil Field in Northern Michigan

The goal of the Dover 33 VSP study was to test the effectiveness of time-lapse VSP for detecting and delineating a plume of more than 271,000 tonnes of CO 2 injected into the Brown Niagaran and A-1 Carbonate formations within the Dover 33 reef between March 2013 and September 2016. Five 2D walkaway VSP (WVSP) source lines were acquired by SIGMA3 in September 2016 to investigate the possible time-lapse response in both P-wave and PS-wave seismic data. The data was compared to the same survey geometry acquired in March 2013 by SR2020.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 1—Theory

Summary Core-scale measurements are considered the ground truth that oil and gas or electric utility companies use to predict the migration of fluids such as oil, natural gas, carbon dioxide, or brine deep underground during their extraction or injection operations. To provide a greater understanding of petrophysical properties of low-permeability geologic formations such as shales and tight gas sandstones, this study introduces a novel core-analysis procedure. The technique follows conventional pressure-pulse-decay permeametry, where the pressure in an inlet chamber adjacent to a cylindrical core plug undergoes a rapid pressurization, the system is shut in, and the pressure reaches a new equilibrium. However, unlike a standard unidirectional pressure-pulse decay, the full-immersion pressure-pulse decay (Hannon 2019) applies a pressure disturbance to the entire outer surface area of the sample. This article covers the numerical simulator designed to model flow through an anisotropic porous sample in this scenario. The model assumes distinct but uniform permeabilities along the radial and axial directions of the cylindrical plug sample. When extracting a plug vertically (or perpendicular to bedding), the permeability along the radial direction associates with the horizontal permeability (i.e., parallel to bedding), whereas in the axial direction, flow occurs perpendicular to bedding (a vertical permeability). Investigations of these model outputs demonstrate an approximately 20-fold decrease in time to complete a full-immersion experiment compared with conventional pressure-pulse decay. Furthermore, the pressure-decay curves resulting from the full-immersion method have slightly different shapes than those resulting from other unidimensional transient methods. These differences begin to demonstrate that, under achievable experimental conditions, the analysis of pressure data from one full-immersion test could enable the simultaneous estimation of the apparent permeabilities parallel and perpendicular to bedding of a cylindrical sample in addition to its porosity. A follow-up article finalizes the proof of this capability with a parameter-estimation procedure and presents experimental verification through a proof-of-concept study. Described in greater detail in that article, the parameter estimator requires multiple forward simulations to analyze the data, which behooves minimizing the compute time for each simulation. By using an alternating direction implicit time-marching scheme and a structured but variably spaced grid, the numerical simulator built for this purpose provides a forward model output with suitable accuracy in approximately 0.5 seconds.

Energy & Fuels↗