Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel state estimation”

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 109 records · Page 6

Antiparallel interface coupling evidenced by negative rotatable anisotropy in IrMn/NiFe bilayers

Negative rotatable anisotropy is estimated via ferromagnetic resonance measurements in as-made, annealed, and ion-irradiated IrMn{sub 3}/Ni{sub 81}Fe{sub 19} bilayers. Opposite to previous observations, inverse correlation between rotatable anisotropy and coercivity is observed. The exchange-bias field, determined from hysteresis loop measurements, is higher than that obtained from ferromagnetic resonance for all samples. The results are discussed in terms of majority antiparallel coupling and magnetic-field-induced transitions from antiparallel to parallel states of uncompensated spins at ferromagnet/antiferromagnet interface. We affirm that an observation of negative rotatable anisotropy evidences antiparallel coupling even in systems presenting conventional exchange bias.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Editorial to the special issue on the 2015 TOUGH Symposium

This Special Issue of Transport in Porous Media features selected papers from the 2015 TOUGH Symposium, which took place on September 28–30, 2015, in Berkeley, California, USA. Such symposia are held every 3 years and are a forum where users of the TOUGH family of numerical simulators meet for an open exchange on the code applications and recent enhancements. The TOUGH codes (TOUGH is an acronym for Transport of Unsaturated Groundwater and Heat) were originally developed in the early 1980s for geothermal reservoir simulation by a team of researchers at the Lawrence Berkeley National Laboratory (LBNL) headed by Karsten Pruess (see http://esd.lbl.gov/TOUGH). The TOUGH codes have a modular structure in which different fluid components can be incorporated into the code by means of a replaceable equation of state (EOS) package. Since its first release, a variety of important EOS modules have been developed, enabling the modeling of volatile organic compounds, decay chains of multiple radionuclides, gas hydrates, tight and shale gas behavior, enhanced oil recovery, and supercritical water and CO 2 . Derivatives and descendants of the basic TOUGH code include coupled multi-phase flow and reactive chemical transport; coupled thermo-hydrologic-mechanical processes; inverse modeling for parameter estimation, sensitivity analysis, and uncertainty quantification; and parallelization for multi-core PCs, workstation clusters, and supercomputers. The TOUGH codes have been applied to a large variety of problems, including geothermal and petroleum reservoir engineering, nuclear waste and CO 2 geologic sequestration, subsurface environmental remediation, methane gas hydrate production, compressed air energy storage, and many other energy and environmental applications. Finally, the international TOUGH user community currently numbers almost 400 organizations (research laboratories, government agencies, private companies, universities, etc.) in over 60 countries.

99 GENERAL AND MISCELLANEOUS↗

A Trichomonas vaginalis C2-XYPPX-repeat protein with a structured C2 domain displaying dampened flexibility upon binding calcium

C2 domains are ubiquitous membrane-binding modules of ∼130 residues in eukaryotes that are often associated with proteins involved in membrane trafficking and lipid modification. The genome of Trichomonas vaginalis, the most common, non-viral, sexually transmitted human pathogen, encodes eight genes that contain a N-terminal C2 module linked to a XYPPX-repeat domain of more than four XYPPX repeats (C2-XYPPX). While the function of the XYPPX-repeat domain remains unknown, its multiple association with C2 domains in T. vaginalis suggests it is important. Here, the C2 domain from one of these C2-XYPPX-repeat proteins, Tv-C2-1, was structurally and physically characterized using X-ray crystallography and NMR spectroscopy. The crystal structure for Tv-C2-1 shows that this domain shares a fold common to all C2 domains, a compact Greek-key motif composed of eight anti-parallel β-strands in the type-2 topology. An NMR chemical shift perturbation study with Ca 2+ showed that Tv-C2-1 bound two Ca 2+ atoms primarily via two loops (loop-1 and loop-3) on the predicted calcium binding face of the protein with K d s of 58.0 ± 0.1 μM and 232 ± 6 μM. Estimations of the overall rotational correlation time, τ c , in the apo (11.1 ns) and Ca 2+ -bound (9.2 ns) state suggests the protein becomes more compact upon Ca 2+ binding, consistent with a decrease in dynamics in loop-3 and marginally in loop-1 suggested by amide 15 N heteronuclear steady-state { 1 H}- 15 N NOEs. Showing Tv-C2-1 binds calcium and adopts a compact Greek-key motif structure, two primary features of C2 domains, suggests understanding the function of the XYPPX-repeat domain may be warranted.

NMR spectroscopy↗

Stochastic Vector Techniques in Ground-State Electronic Structure

Herein we review a suite of stochastic vector computational approaches for studying the electronic structure of extended condensed matter systems. These techniques help reduce algorithmic complexity, facilitate efficient parallelization, simplify computational tasks, accelerate calculations, and diminish memory requirements. While their scope is vast, we limit our study to ground-state and finite temperature density functional theory (DFT) and second-order many-body perturbation theory. More advanced topics, such as quasiparticle (charge) and optical (neutral) excitations and higher-order processes, are covered elsewhere. We start by explaining how to use stochastic vectors in computations, characterizing the associated statistical errors. Next, we show how to estimate the electron density in DFT and discuss effective techniques to reduce statistical errors. Finally, we review the use of stochastic vectors for calculating correlation energies within the second-order Møller-Plesset perturbation theory and its finite temperature variational form. Example calculation results are presented and used to demonstrate the efficacy of the methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparing Parallel Plastic‐to‐X Pathways and Their Role in a Circular Economy for PET Bottles

Abstract The United States generates the most plastic waste of any country and is a top contributor to global plastic pollution. Multiple end‐of‐life strategies must be implemented to minimize environmental impacts and retain valuable plastic material, but it is challenging to compare options that generate products with different lifetimes and utilities. Herein, they present a material flow model equipped with consequential life cycle assessment, cost analysis, and a plastic circularity indicator that considers product quality and lifetime. The model is used to estimate the greenhouse gas (GHG) emissions, circularity, and cost of polyethylene terephthalate (PET) bottle mechanical downcycling to lower‐quality resin, closed‐loop glycolysis to food‐grade PET, upcycling to glass fiber‐reinforced plastic, and conversion to non‐plastic products (electricity, oil) on a United States economy‐wide basis for the year 2020. A brute force algorithm suggests that a combination of 68% glycolysis, 11% mechanical recycling, 6% upcycling, 9% landfilling, and 5% incineration can minimize the cost and GHG emissions and maximize the circularity of the current PET economy. However, uncertainty around transportation distances, materials recovery facility efficiencies, and recycling yields can result in different “optimal” pathway mixes. This flexible framework enables informed decision‐making to move toward a cost‐ and environment‐conscious circular economy for plastic.

09 BIOMASS FUELS↗

GRAVITATIONAL WAVE SIGNATURES OF HYPERACCRETING COLLAPSAR DISKS

By performing two-dimensional special relativistic (SR) magnetohydrodynamic simulations, we study possible signatures of gravitational waves (GWs) in the context of the collapsar model for long-duration gamma-ray bursts. In our SR simulations, the central black hole is treated as an absorbing boundary. By doing so, we focus on the GWs generated by asphericities in neutrino emission and matter motions in the vicinity of the hyperaccreting disks. We compute nine models by adding initial angular momenta and magnetic fields parametrically to a precollapse core of a 35 M{sub Sun} progenitor star. As for the microphysics, a realistic equation of state is employed and the neutrino cooling is taken into account via a multi-flavor neutrino leakage scheme. To accurately estimate GWs produced by anisotropic neutrino emission, we perform a ray-tracing analysis in general relativity by a post-processing procedure. By employing a stress formula that includes contributions from both magnetic fields and SR corrections, we also study the effects of magnetic fields on the gravitational waveforms. We find that the GW amplitudes from anisotropic neutrino emission show a monotonic increase with time, whose amplitudes are much larger than those from matter motions of the accreting material. We show that the increasing trend of the neutrino GWs stems from the excess of neutrino emission in the direction near parallel to the spin axis illuminated from the hyperaccreting disks. We point out that a recently proposed future space-based interferometer like Fabry-Perot-type DECIGO would permit the detection of these GW signals within Almost-Equal-To 100 Mpc.

79 ASTRONOMY AND ASTROPHYSICS↗

Two-Tower Quantum Matrix Chain Multiplication: Trading Qubits for Depth

Matrix chain multiplication -- computing $\mathcal{W} = M^{(0)}\cdots M^{(K-1)}$ where $M^{(k)} \in \mathbb{R}^{P_k \times P_{k+1}}$-- arises in scientific computing, machine learning, and graph analysis. Despite the importance of this problem, for chains of distinct matrices, the classical number of operations grows linearly with the chain length $K$ and polynomially in the matrix dimensions. We present \emph{Two-Tower Matrix Multiplication}, a quantum subroutine that encodes the product $\mathcal{W}$ of the $K$ matrices into a quantum state in circuit depth $\mathcal{O}(\max_{k} \mathrm{polylog} (P_k P_{k+1}))$, which is independent of~$K$ within the QRAM-based state-preparation model, whereas the qubit count is $\mathcal{O}\bigl(\sum_{k} \log P_k \bigr)$; the total gate count remains linear in $K$, so the gain is in the circuit depth. The construction interleaves state-preparation operators across two layers; within each layer, all operators act on disjoint registers and execute in parallel. This subroutine can be specialized for the chain-vector case, which computes the product of $K-1$ matrices applied to a vector. We prove the correctness of the subroutine for all $K$ and provide two implementations using the Qiskit and QCLAB frameworks. The subroutine is applicable to any downstream quantum algorithm that operates on a matrix encoded in the statevector, including norm estimation, graph-matrix powers, linear system solving, and quantum machine learning kernels.

Antonioli, Giacomo [Pisa U.] (ORCID:00090000668703↗

CO 2 Capture by Cold Membrane Operation with Actual Power Plant Flue Gas

The main objective of the project was to develop a post-combustion CO 2 capture process based on the hybrid cold temperature membrane operation. The CO 2 in the flue gas from coal fired power plant is pre-concentrated to >60% CO 2 in the first stage membrane operation followed by further liquefaction of permeate stream to achieve >99% CO 2 purity. The aim of the project was based on DOE program goal of 90% CO 2 capture with >95% CO 2 purity from Pulverized Coal (PC) fired power plants with $40/tonne of carbon capture cost by 2025. The project moves the technology from TRL 4 to TRL 5. The project involved optimization of Air Liquide commercial 12” PI-1 bundle to improve the bundle productivity by >30% compared to the previous baseline (DE-FE0004278) using computational fluid dynamics (CFD) modeling and bundle testing with synthetic flue gas at 0.1 MWe bench scale skid located at Delaware Research and Technology Center (DRTC). In parallel, the next generation polyimide based novel PI-2 membrane was developed with 10 times CO 2 permeance compared to the commercial PI-1 membrane. The novel PI-2 membrane was scaled from mini-permeator to 1” permeator and 1” bundle for testing. Bundle development was conducted with a Development Spin Unit (DSU) installed at MEDAL. Air Liquide’s cold membrane technology was demonstrated with real coal fired flue gas at the National Carbon Capture Center (NCCC) with a 0.3 MWe field-test unit (FTU). The FTU was designed to incorporate testing of two PI-1 commercial membrane bundles (12” or 6” diameter) in parallel or series. A slip stream was sent to the next generation PI-2 membrane for testing with real flue gas. The system exceeded performance targets with stable PI-1 membrane operation for over 500 hours of single bundle, steady state testing. The 12” PI-1 bundle exceeded the productivity target by achieving ~600 Nm3/hr, where the target was set at ~455 Nm3/hr at 90% capture rate. The cost of 90% CO 2 capture from a 550 MWe net coal power plant was estimated between 40 and $45/tonne. A 6” PI-1 bundle exhibited superior bundle performance compared to the 12” PI-1 bundle. However, the carbon capture cost was not lower with the 6” PI-1 bundle due to the higher bundle installed cost. A 1” PI-1 bundle was tested to compare bundles with different length / diameter ratios. This bundle exhibited the lowest performance due to the different fiber winding pattern and increased bundle non-ideality. Several long-term and parametric tests were conducted with 3,200 hours of total run-time at NCCC. Finally, the new PI-2 membrane fiber was tested at a small scale (1” modules) in real flue gas and exhibited up to 10 times the CO 2 permeance and slightly lower CO 2 /N 2 selectivity as the commercial PI-1 fiber. This corresponded to a projected 4 - 5 times increase in the productivity per bundle and a potential cost reduction of $3/tonne for CO2 capture, as compared with PI-1. An analytical campaign was conducted to trace different impurities such as NOx, mercury, Arsenic, Selenium in gas and liquid samples through the carbon capture system. An Environmental, Health and Safety (EH&S) analysis was completed to estimate emissions from a 550 MWe net power plant with carbon capture using cold membrane. A preliminary design and cost analysis was completed for 550 tpd (~25 MWe) plant to assess the capital investment and carbon capture cost for PI-1 and PI-2 membrane solutions from coal fired flue gas. A comparison was made with an amine based solution with significant cost advantage for the membrane at this scale. Additional preliminary design and cost analysis was completed between coal, natural gas and SMR flue gas for carbon capture at 550 tpd (~25 MWe) plant.

01 COAL, LIGNITE, AND PEAT↗

Properties of the fast forward shock driven by the 2012 July 23 extreme coronal mass ejection

Late on 2012 July 23, the STEREO-A spacecraft encountered a fast forward shock driven by a coronal mass ejection (CME) launched from the Sun earlier that same day. The estimated travel time of the disturbance (∼20 hr), together with the massive magnetic field strengths measured within the ejecta (>100 nT), made it one of the most extreme events observed during the space era. In this study, we examine the properties of the shock wave. Because of an instrument malfunction, plasma measurements during the interval surrounding the CME were limited, and our approach has been modified to capitalize on the available measurements and suitable proxies, where possible. We were able to infer the following properties. First, the shock normal was pointing predominantly in the radial direction (n=0.97e{sub r}−0.09e{sub t}−0.23e{sub n}). Second, the angle between n and the upstream magnetic field, θ{sub Bn}, was estimated to be ≈34°, making the shock “quasi-parallel,” and supporting the idea of an earlier “preconditioning” ICME. Third, the shock speed was estimated to be ≈3300 km s{sup −1}. Fourth, the sonic Mach number, M{sub s}, for this shock was ∼28. We support these results with an idealized numerical simulation of the ICME. Finally, we estimated the change in ram pressure upstream of the shock to be ∼5 times larger than the pressure from the energetic particles, suggesting that this cosmic-ray modified shock had not reached steady-state, but instead, had been caught in an early, transient phase in its evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

A novel adaptive time stepping variant of the Boris–Buneman integrator for the simulation of particle accelerators with space charge

We show that adaptive time stepping in particle accelerator simulation is an enhancement for certain problems. The new algorithm has been implemented in the OPAL (Object Oriented Parallel Accelerator Library) framework. The idea is to adjust the frequency of costly self-field calculations, which are needed to model Coulomb interaction (space charge) effects. In analogy to a Kepler orbit simulation that requires a higher time step resolution at the close encounter, we propose to choose the time step based on the magnitude of the space charge forces. Inspired by geometric integration techniques, our algorithm chooses the time step proportional to a function of the current phase space state instead of calculating a local error estimate like a conventional adaptive procedure. Building on recent work, a more profound argument is given on how exactly the time step should be chosen. An intermediate algorithm, initially built to allow a clearer analysis by introducing separate time steps for external field and self-field integration, turned out to be useful by its own, for a large class of problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mitigating Cascading Outages in Severe Weather Using Simulation-Based Optimization

Severe weather events can trigger cascading power outages and lead to significant losses. In this work, we investigate cascading outage mitigation under severe weather conditions. Given day-ahead weather forecasts and component failure models, we aim to identify a set of power lines that can be hardened to minimize the expected impact of potential cascading outages. Since the expected load shedding cannot be expressed as an explicit function of line hardening decisions and system states, we developed a cascading outage simulator to estimate the expected value of load shedding under various initial weather-related disruption scenarios generated using a weather forecast. To avoid massive enumeration of all possible combinations of line hardening decisions and reduce the simulation efforts, we employed an efficient simulation-based optimization approach that quickly identifies the (near) optimal line hardening decisions in the presence of both large simulation noises due to the highly variable initial disturbances and system states, and significant randomness in the subsequent cascades. Furthermore, the algorithm is also able to utilize parallel computing to dramatically reduce computation time to support decision making in preparation for severe weather conditions. We performed a case study on the Northeast Power Coordinating Council (NPCC) 140-bus system model to demonstrate that our approach can significantly improve power grid resilience to adverse weather events.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Direct Demonstration of Unified Brønsted-Evans-Polanyi Relationships for Proton-Coupled Electron Transfer Reactions on Transition Metal Surfaces

Brønsted-Evans-Polanyi (BEP) relationships, which relate elementary reaction barriers to reaction thermodynamics, have long been discussed in electrochemical science, but demonstration of their existence in heterogeneous electrocatalysis, across many different catalyst surfaces and voltages, is lacking. In this contribution, the BEP model is demonstrated to describe both multiple catalyst surfaces and variable voltages in electrochemical environments. For proton-coupled electron transfer (PCET) to surface nitrogen (N*) and nitric oxide (NO*), reaction energies and activation barriers are calculated using Density Functional Theory (DFT) on a parallel plate capacitor model at three different potentials and for nine different transition metal surfaces. Linear BEP relations that describe all potentials and catalyst surfaces are obtained for these elementary reactions, and the slopes of the correlations are shown to be directly related to the fractional coordinate of the transition states (FCTS) of the reactions. The results, which are explained in terms of Marcus Theory, prove a direct equivalence between unified BEP coefficients, describing both variable catalyst surfaces and voltages, and electrochemical symmetry factors and provide a straightforward means of estimating this quantify for PCET reactions on transition metal surfaces. The resulting relationships, in turn, could lead to predictions of electrocatalytic reactivity trends of enhanced accuracy and efficiency.

Electrochemistry↗

Transport Estimations of Final State Interaction Effects on Short-range Correlation Studies Using the $(e,e'p)$ and $(e,e'pp)$ Reactions

Short range correlated (SRC) nucleon-nucleon pairs in nuclei are typically studied using measurements of electron-induced hard nucleon-knockout reactions (e.g. $(e,e'p)$ and $(e,e'pN)$), where the kinematics of the knocked-out nucleons are used to infer their initial state prior to the interaction. The validity of this inference relies on our understanding of the scattering reaction, most importantly how rescattering of the detected nucleons (final state interactions or FSI) distort their kinematical distributions. Recent SRC measurements on medium to heavy nuclei have been performed at high-$x_B$ (i.e., anti-parallel kinematics) where calculations of light nuclei indicate that such distortion effects are small. Here we study the impact of FSI on recent $^{12}$C$(e,e'p)$ and $^{12}$C$(e,e'pp)$ measurements using a transport approach. We find that while FSI can significantly distort the measured kinematical distributions of SRC breakup events, selecting high-$x_B$ anti-parallel events strongly suppresses such distortions. In addition, including the effects of FSI improves the agreement between Generalized Contact Formalism-based calculations and data and can help identify those observables that have minimal sensitivity to FSI effects. This result helps confirm the interpretation of experimental data in terms of initial-state momentum distributions and provides a new tool for the study of SRCs using lepton-scattering reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Impact of Different Thermal Gradients on the Dynamics of Cylindrical Lithium-ion Cells Subject to Accelerated Aging and on Module Performance

This study investigates the impacts of applying different thermal gradient patterns to cylindrical lithium-ion cells in a module on cell dynamics (temperatures, current flows, state of charge), module performance (evolution of resistance, capacity, and energy versus cycle number), and module lifetime. The thermal gradients were generated using cooling plates (CPs) with three different flow-field designs, namely, straight, perpendicular, and U-turn. The study uses computational fluid dynamics (CFD), the pseudo-two-dimensional (P2D) battery model, capacity loss and increased impedance due to the growth of a solid-electrolyte-interphase, and the electric current distribution from module terminals to cells that depends on the series-parallel electrical connections among the cells. The impact of the thermal gradient (resulting from the CP designs) on the variability in resistance, current, state of charge, and voltage among the cells was analyzed and linked to differences in the module's performance. Applying a thermal gradient to parallel-connected strings of series-connected cells led to variation in the current through each parallel string and an imbalance in the voltage of series-connected cells. Module performance is poorer when the thermal gradient causes a voltage imbalance than when it causes a current imbalance. Module performance becomes the worst when both current variation and voltage imbalance happen together. For instance, the module's lifetime (estimated as reaching 80% of its initial capacity) varied by 5% to 17.5%, depending on the magnitude and pattern of the imposed thermal gradient. As the relative orientation between thermal gradients and cells' electrical connectivity influences the module's performance, appropriate consideration should be given to the choice of the CP, especially if large thermal gradients are allowed.

Battery thermal management↗

Properties of Lower Hybrid Drift Waves and Energy Transfer Near and Inside the Magnetic Reconnection Electron Diffusion Regions

Here, we investigate properties of lower hybrid drift waves (LHDWs) near and inside the electron diffusion regions in 17 magnetopause and 9 magnetotail reconnection events using Magnetospheric MultiScale (MMS) mission observations. Our analysis show that LHDW type depend on the electron beta, as electron beta increases LHDWs become more electromagnetic in nature. The energy transfer from electromagnetic fields to particles is higher in electrostatic LHDWs and it is largely in parallel direction with respect to the local magnetic field. Linear dispersion analysis shows that electrostatic LHDWs are perpendicular propagating while electromagnetic waves propagate in oblique directions and the normalized wavenumber of all LHDW types falls within 0.5–0.8 range. A simple estimate on the LHDW nonlinear saturation suggests a possibly important roles played by these waves in supporting the reconnection electric field.

79 ASTRONOMY AND ASTROPHYSICS↗

PyPWA: A Partial-Wave/Amplitude Analysis Software Framework

The PyPWA project aims to develop a software framework for Partial Wave and Amplitude Analysis of data; providing the user with software tools to identify resonances from multi-particle final states in photoproduction. Most of the code is written in Python. The software is divided into two main branches: one general-shell where amplitude's parameters (or any parametric model) are to be estimated from the data. This branch also includes software to produce simulated data-sets using the fitted amplitudes. A second branch contains a specific realization of the isobar model (with room to include Deck-type and other isobar model extensions) to perform PWA with an interface into the computer resources at Jefferson Lab. We are currently implementing parallelism and vectorization using the Intel's Xeon Phi family of coprocessors.

Salgado, Carlos↗

Effect of magneto-mechanical synergism in the process-structure correlation in Fe–C alloys: A phase-field modeling approach

Applied magnetic fields can alter phase equilibria and kinetics in steels; however, quantitatively resolving how magnetic, chemical, and elastic driving forces jointly influence the microstructure remains challenging. We develop a quantitative magneto-mechanically coupled phase-field model for the Fe–C system that couples a CALPHAD-based chemical free energy with demagnetization-field magnetostatics and microelasticity. Here, the model reproduces single- and multi-particle evolution during the α → γ inverse transformation at 1023 K under external fields up to 20 T, including ellipsoidal morphologies observed experimentally at 8 T. Chemically driven growth is isotropic; a magnetic interaction introduces an anisotropic driving force that elongates γ precipitates along the field into ellipsoids, while elastic coherency promotes faceting, yielding elongated cuboidal or “brick-like” particles under combined magneto-elastic coupling. Growth kinetics increase with C content, and decrease with field strength and misfit strain. Multi-particle simulations reveal dipolar interaction-mediated coalescence for field-parallel neighbors and ripening for field-perpendicular neighbors. Incorporating field-dependent diffusivity from experiment slows kinetics as expected; a first-principles-motivated anisotropic diffusivity correction is estimated to be small (<2%). These results establish a process-structure link for magnetically assisted heat treatments of Fe–C alloys and provide guidance for microstructure control via chemo-magneto-mechanical synergism.

Magnetic field↗

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 2—Practice and Demonstration

Summary This study continues the work of presenting a novel approach for making petrophysical assessments of tight core samples. This method, the full-immersion pressure-pulse decay (Hannon 2019), involves applying a rapid increase in pressure in a chamber surrounding the entire outer surface area of a cylindrical sample, shutting the system in, and monitoring the pressure decay in the chamber as it reaches a new equilibrium. A precursor article covered the numerical simulator designed to model flow through the sample, demonstrating its performance and accuracy in addition to providing a first-order comparison between the speed and shape of the pressure-decay responses of the full-immersion method with those of other similar transient methods. This study covers the parameter-estimation procedure and experimental verification through a proof-of-concept laboratory investigation. The investigations provided here demonstrate that under appropriate, achievable experimental conditions, the pressure data can be analyzed in such a way that returns an estimate of the porosity and apparent permeabilities both parallel and perpendicular to bedding from a single test performed on a single cylindrical sample. After determining these experimental conditions (the uniqueness window), this report outlines a data-inversion strategy to estimate the petrophysical properties (porosity, horizontal permeability, and vertical permeability) from each test. This strategy is put to the test through comparisons with measurements performed by a commercial core laboratory. A common set of samples recovered from an outcrop of a tight-gas sandstone formation were investigated using the full-immersion method, and their results are compared with those from conventional steady-state measurement procedures performed by the commercial laboratory. Comparisons between petrophysical characterizations of these samples, which had permeabilities between 25 nd and 2.3 μd, demonstrated close agreement in most cases. However, whereas steady-state measurements performed at the professional laboratory required 4 to 5 hours of testing time per measurement of a single permeability, similar assessments using the full-immersion technique, requiring approximately 5 to 10 minutes to complete, returned estimates of the horizontal and vertical permeability simultaneously. Additional analyses are provided to determine principal reasons of discrepancies in instances where agreement was not as strong. Based on lessons learned from these experiences, the report closes with suggestions on areas of improvement in the experimental approach. Once complete, these developments should propel this technology to fill a critical need to determine petrophysical properties (porosity and permeability) of tight rocks in a time-efficient manner and in a way that does not compromise their accuracy.

Energy & Fuels↗