SEARCH · Engineering Papers

Results for “Bed”

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 991 records · Page 55

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Chapter 14: Heat and Mass Transfer

Fluidized beds are used in a variety of processes due to their excellent heat and mass transfer properties. Heat transfer in fluidized beds can refer to interphase heat transfer, or more commonly to bed-to-surface heat transfer. Heat transfer between bubble phase and dense phase can be characterized by an interphase heat transfer coefficient. This chapter elaborates the derivation of the interphase heat transfer coefficient by an analogy to mass transfer. Mass transfer in fluidized beds usually refers to gas-to-particle mass transfer, or interphase mass transfer, where a component migrates from particles to fluids, or vice versa. The chapter describes the two major components of the interphase mass transfer: particle and fluid mass transfer in the dense phase; and bubble- to dense-phase mass transfer.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗

DOE OSTI2020

Thermal characterization of the build chamber in electron beam melting

Electron beam powder bed fusion, commonly termed electron beam melting (EBM), offers great versatility in multiple-part processes and can produce high quality as-built components due to low residual stresses. The EBM process is complex and requires careful thermal management, including uniform and consistent preheating of the powder bed, in order to ensure quality and consistency of the product. However, most of the simulations in the literature focus on the selective melting stage of the process. As of today, optimal conditions for initial pre-heating temperatures are only available for specific materials. Thus, in order to extend the EBM technology to other desired build materials, a much better understanding of the pre-heating stages is required. In this work, numerical and experimental approaches are combined in order to investigate the effects and sensitivities of heat removal from the build plate during EBM pre-heating stages. For this purpose, a carefully reconstructed numerical model of the build chamber of an ARCAM Q20 + machine is developed. It includes all main parts of the chamber and all relevant heat transfer mechanisms, whereas special attention is paid to radiation heat exchange between various bodies. In order to validate the model, dedicated experiments are performed, in which a system of thermocouples is installed in the build chamber, allowing direct measurement of the local temperatures of the start-plate and heat shields. A good agreement between the simulation and experimental findings is achieved, leading to a better understanding of the thermal processes characteristic to the pre-heating stages. This basic analysis is followed by a representative pre-heating case, where a powder bed is also considered. The energy required to obtain the desired pre-heating temperatures is evaluated, and the role of the powder bed in heat transfer within the chamber is assessed. The pre-heating stage, simulated in the present work, is supposed to create proper conditions for sintering and consequent melting of the powder, which are highly dependent on the local temperatures and heat transfer features. Thus, the findings of the reported work present a step towards a better understanding of the thermal processes that characterize EBM. The reported model can be further used to provide realistic boundary condition inputs for other meso- or macro-scale models as a function of time and geometry. The model can serve also for verification of machine settings, i.e., jump-safe and melt-safe ones which actually provide desired preheat, and for development of settings for new powders.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Effects of slope steepness and cross-slope wind speed on fire spreading behavior for various vegetation

Wildland fire behavior is significantly influenced by environmental factors such as slope steepness, wind speed, and fuel type. Understanding these interactions is critical for improving predictive models and fire management. This study explores how slope steepness and cross-slope wind speed influence fire spread dynamics in various fuel bed types. Simulations are conducted using a physics-based wildland fire model, HIGRAD/FIRETEC, across six slope angles (0–50 %), four cross-slope wind speeds (4–10 m s –1 ), and three fuel bed types (grass, shrubland, and forest). Representative cases are additionally compared with FARSITE fireline evolution. Fire behavior is categorized into distinct propagation types based on spread characteristics and analyzed. The fire propagation angle, which indicates deviation from the wind direction, generally increases with steeper slopes and decreases with stronger cross-slope winds. Secondary upslope propagation is observed in shrubland under moderate slopes, while secondary downwind propagation occurs in all fuel beds at higher wind speeds. These findings highlight fire spread characteristics that differ from predictions by traditional models like Rothermel’s. By capturing complex propagation patterns and dynamics, this study demonstrates the value of a physics-based, atmosphere-fire coupled model for accurate wildland fire prediction. Incorporating secondary propagations and the influence of fuel bed complexities into predictive models can improve the accuracy of fire spread forecasts, enabling more effective fire management and risk mitigation efforts.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Pore accessibility and trapping of methane in Marcellus shale

Accessibility of pore space in unconventional reservoirs is an important factor influencing both methane storage capacity and the kinetics of methane desorption. The determination of open (accessible) versus closed (inaccessible) porosity is therefore vital for the prediction of gas production potential. This study investigates accessibility of pores to methane in overmature middle Devonian Marcellus shale samples (cut parallel and perpendicular to bedding) using small and ultra-small neutron scattering (SANS and USANS) with contrast matching (CM), supplemented by other complementary techniques, such as mercury injection capillary pressure (MICP) and low pressure gas (N2 and CO2) adsorption. Our results demonstrate that for the samples studied, only about 6 % of pores with diameter 25-500 nm are accessible to methane. The accessibility fraction for pores larger than 500 nm is 35%. For nanopores smaller than 25 nm, pore accessibility could not be quantitatively determined due to increased methane density and condensation effects in confinement. Our observations indicate that methane penetrates the accessible small mesopores and micropores down to at least 1 nm in diameter, and the density of confined deuterated methane (CD4) is 0.68 g/cm3 for pores of diameter 25 nm and gradually increases with the decreased pore size. Moreover, elevated gas pressure causes formation of additional high-density methane nano-clusters. These clusters have a form of slightly anisotropic polydisperse discs oriented along bedding plane, about 1-12 nm in diameter and with average thickness of 3.6 nm. Utilizing samples cut parallel and perpendicular to the bedding, this study also briefly addresses anisotropy of pores. Based on the iso-size intensity ratio R, our SANS and USANS results demonstrate anisotropy in the out-of-bedding direction and suggest that the degree of anisotropy depends on the pore size. Specifically, for pore diameters ~2.5 to 250 nm, the extent of anisotropy is smaller than for pores ~500 nm to 6 µm in diameter. Finally, comparison of pore size distribution results calculated from SANS/USANS to those obtained using MICP shows good agreement at low pressures, but large difference at pressures above 1000 bar. Here, this discrepancy requires further testing; it is possible that the high mercury pressure used in MICP alters the mesopore and micropore structure of shales.

03 NATURAL GAS↗

DOE OSTI2021

A micromechanical hydro-mechanical-damage coupled model for layered rocks considering multi-scale structures

An anisotropic Hydro-Mechanical-Damage coupled model was established in the framework of thermodynamics for saturated layered rocks, which contain arbitrarily-distributed cracks of smaller scales and orientated bedding planes of much larger scales. Here, the proposed model considers the anisotropic damage growth and frictional sliding of cracks and shear sliding and dilatancy of bedding planes, as well as the interaction between cracks and bedding planes. Numerical implementation of the nonlinear coupled model was developed based on TOUGHREACT, which is a well-established simulator for multiphase fluid flow and reactive transport analysis. The proposed model was first validated with laboratory water injection test results on rock blocks, and then used to study the effects of injection rates, in situ stresses and bedding planes on fluid injection processes at the field scale, demonstrating the importance to consider the variation of multi-scale structures for better understanding the coupled Hydro-Mechanical responses in layered rocks.

58 GEOSCIENCES↗

DOE OSTI2020

Analysis of hematite attrition in a grid jet apparatus

Particulate attrition is of interest for novel carbon-capture processes such as chemical looping combustion because the makeup cost of oxygen carrier is a significant portion of operating cost. As such, models to study and predict attrition of various oxygen carriers in fluidized bed systems are being developed. One of the regions of concern in fluidized bed systems is the high-velocity jet region near gas distributors in a fluid bed. This work studies the attrition of hematite particles using a modified ASTM apparatus to measure the particle size distribution throughout the experiment. Bed weight and gas velocity were varied. Hematite particles above the corresponding threshold value had a severe variation of particle size distributions which decreased with time. Weight fractions of the sieves over time were fit to a linear, time-variant population balance model to offer insight into particle attrition. The first-order rate constant was modified as a decaying exponential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Microbe-Encapsulated Silica Gel Biosorbents for Selective Extraction of Scandium from Coal Byproducts

Establishing an inexpensive and environmentally friendly scandium (Sc) supply is critical for the development of clean energy technologies. However, refining Sc from Sc-bearing sources using current technologies poses economic and environmental challenges due to its low concentration relative to base metals and chemically similar lanthanides. In this work, we developed a biosorption-based flow-through process for extraction of Sc from low-grade feedstocks. A novel biosorbent material was synthesized by encapsulating Arthrobacter nicotianae, a bacterium that selectively adsorbs Sc, within a porous silica matrix. Through the use of several complementary microscopy techniques, we demonstrated a highly porous biosorbent structure with a high cell-loading density and a homogenous cell distribution. Batch adsorption assays revealed selective Sc adsorption over lanthanides and common based metals, with the exception of Fe(III), as evidenced by separation factors greater than 50. The MESG particles were packed into fixed-bed columns to enable adsorption under flow-through conditions, which demonstrated effective Sc extraction at flow rates up to 0.08 cm/s and high stability for reuse; greater than 95% of the adsorption capacity was maintained after 10 consecutive adsorption/desorption cycles. When applied to a lignite coal leachate, whose Fe content was depleted through pH-mediated precipitation, the MESG particles yielded Sc breakthrough at 30 bed volumes, whereas all other metals broke through after only a few bed volumes. Following desorption, an eluate with a 124-fold increase in Sc purity was achieved relative to the lignite leachate with Sc constituting 96.4% of the total REEs. This study established a rapid, facile, and scalable cell immobilization approach that enables the use of microbial biomass for Sc recovery from low-grade sources.

Dong, Ziye↗

DOE OSTI2020

Optimizing Ammonia Separation via Reactive Absorption for Sustainable Ammonia Synthesis

Metal halide salts such as magnesium chloride have been shown to be promising candidates for ammonia storage materials for energy storage and agriculture applications due to their ability to incorporate several moles of ammonia per mole of salt. Ammonia exiting a synthesis reactor can be separated from nitrogen and hydrogen by absorption into magnesium chloride. Such an absorption can be more complete and hotter than separation via ammonia condensation, the current standard in the Haber–Bosch process. Here, we discuss the optimal conditions for the cyclic uptake and release of ammonia from the supported magnesium chloride absorbents. An automated system was designed for measuring the nonequilibrium working capacity of the absorbent, as well as the impact of important operating conditions such as absorption and desorption temperature, pressure, and desorption time. Measurements of absorption and desorption kinetics provide insight into the mechanisms involved. The temperatures and pressures during absorption and desorption were designed to use minimal energy input to maximize the uptake and release of ammonia within a reasonable amount of time. In a laboratory-scale bed, absorption has a small unused bed length, so it is largely independent of mass transfer; it is dominated by how fast ammonia is fed to the bed. Yet, desorption is restricted both by the speed of heating the bed and by diffusion out of the absorbent. These measurements provide guidelines for ammonia separations and cycling sorbent materials on a larger scale.

10 SYNTHETIC FUELS↗

DOE OSTI2021

Selective butene formation in direct ethanol-to-C3+-olefin valorization over Zn-Y/Beta and single-atom alloy composite catalysts using in situ generated hydrogen

The selective production of C3+ olefins from renewable feedstocks, especially via C1 and C2 platform chemicals, is a critical challenge for obtaining economically viable low-carbon middle distillate transportation fuels (i.e., jet and diesel). Here, we report a multifunctional catalyst system composed of Zn-Y/Beta and “single-atom” alloy (SAA) Pt-Cu/Al2O3 which selectively catalyzes ethanol-to-olefins (C3+, ETO) valorization in the absence of cofed hydrogen, forming butenes as the primary olefin products. Beta zeolites containing predominately isolated Zn and Y metal sites catalyze ethanol upgrading steps (588 K, 3.1 kPa ethanol, ambient pressure) regardless of cofed hydrogen partial pressure (0-98.3 kPa H2), forming butadiene as the primary product (60% selectivity at 87% conversion). The Zn-Y/Beta catalyst possesses site-isolated Zn and Y Lewis acid sites (at ~7 wt% Y) and Brønsted acidic Y sites, the latter of which has been previously uncharacterized. A secondary bed of SAA Pt-Cu/Al2O3 selectively hydrogenates butadiene to butene isomers at a consistent reaction temperature using hydrogen generated in situ from ethanol-to-butadiene (ETB) conversion. This unique hydrogenation reactivity at near-stoichiometric hydrogen and butadiene partial pressures is not observed over monometallic Pt or Cu catalysts, highlighting these operating conditions as a critical SAA catalyst application area for conjugated diene selective hydrogenation at high reaction temperatures (>573 K) and low H2/diene ratios (e.g., 1:1). Single-bed steady state selective hydrogenation rates, associated apparent hydrogen and butadiene reaction orders, and DFT calculations of the Horiuti-Polanyi reaction mechanisms indicate that the unique butadiene selective hydrogenation reactivity over SAA Pt-Cu/Al2O3 reflects lower hydrogen scission barriers relative to monometallic Cu surfaces and limited butene binding energies relative to monometallic Pt surfaces. DFT calculations further indicate the preferential desorption of butene isomers over SAA Pt-Cu(111) and Cu(111) surfaces while Pt(111) surface favors subsequent butene hydrogenation reactions to form butane over butene desorption events. Under operating conditions without of hydrogen cofeeding, this combination of Zn-Y/Beta and SAA Pt-Cu catalysts can selectively form butenes (65% butenes, 78% C3+ selectivity at 94% conversion) and avoid butane formation using only in situ generated hydrogen, avoiding costly hydrogen cofeeding requirements that hinder many renewable energy processes.

Cordon, Michael↗

DOE OSTI2020

Performance Restoration of a Tritium-Aged LaNi 4.25 Al 0.75 Sample

Hydride beds containing LaNi 4.25 Al 0.75 (LANA.75) are used to store significant quantities of tritium. These hydride beds have a limited service life due to radiolytic decay of tritium to 3He within the metal matrix. The crystal structure of the hydride is altered by trapped 3 He, which has a very low solubility in the metal. The altered structure induces the formation of a heel of trapped hydrogen isotopes and diminishes the reversible capacity of the hydride. With sufficient tritium exposure, the bed loses the ability to deliver 3He-free tritium, and replacement is needed. Demonstration of a means to regenerate tritium-aged LANA.75 in situ would delay or even eliminate the need to replace lanthanum nickel aluminum (LANA) hydride beds. This paper presents test results obtained during regeneration testing. The efficacy of regeneration testing was evaluated by comparing tritium desorption isotherms collected on the hydride before and after exposure to regeneration conditions. Testing was performed on a bench-scale tritium-aged LANA.75 sample that was previously isotopically exchanged (from tritium to deuterium), passivated, and recovered. Once transferred to a high-temperature test cell, the deuterium heel of the sample was isotopically exchanged with tritium, and a baseline desorption isotherm was collected for comparison purposes. The sample was then heated under vacuum, and comparative isotherms were gathered between regeneration evolutions. Shifts in isotherms show progressive improvements with higher-temperature exposure over the tritium-aged baseline. The heel was significantly reduced, and the reversible capacity of the hydride was essentially restored to near virgin values. For all tested conditions, the plateau pressure remained higher than virgin LANA.75.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Shock interacting with a random array of stationary particles underwater

Accurately predicting the force experienced by particles within a random bed during shock traversal is a challenging problem. In this paper, particle-resolved inviscid simulations of shock propagation over randomly distributed beds are performed, with the goal of quantifying the force on individual particles. Water was considered the continuum material, and then results were compared with previous research completed in air and the differences resulting from these different mediums were examined. Simulations were conducted for four different combinations of incident Mach number and bed volume fraction with the particles remaining stationary. Additionally, time-resolved streamwise and transverse force coefficients were calculated for each particle in the bed. It was observed that while the average force coefficient was similar to that of an isolated particle, the force on individual particles was substantially different. An important observation was that the inviscid drag force on some particles remained consistently positive, but others consistently negative even long after the passage of the shock. This persistent force on the particle was attributed to the potential flow fields that result from neighbors' specific locations upstream and downstream. Similar behavior was observed in the transverse force, with particles experiencing sustained force in the same lateral direction. These persistent forces are inviscid quasisteady contributions which with the viscous counterpart will play an essential role in the particle's long-term dispersion if allowed to move in response to the force.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Resilience Metrics for Building-Level Electrical Distribution Systems with Energy Storage

The energy system infrastructure that delivers power to a building's loads needs to be resilient to withstand and recover from extreme outages (e.g., grid faults that leave millions of people without power during severe weather events). Building-level electrical distribution systems (BEDSs) distribute power from a building's energy sources - including the grid, solar photovoltaic (PV) panels, and batteries - to its loads, including lighting, HVAC, and plug loads. BEDS with storage can provide resilience by distributing local electricity supply to critical loads during an outage. Quantitative metrics are needed to assess the resilience improvements associated with new BEDS and storage system technologies. In this paper, we apply an existing metric, the probability of outage survival curve (POSC), to BEDS with storage and propose novel metrics that improve upon POSC. Through a simulation-based case study, we demonstrate how these metrics are impacted by the BEDS design and how they can be used to design a resilient system.

buildings↗

DOE OSTI2024

Development of a Hardware-in-The-Loop Testbed for a Decentralized, Data-Driven Electric Vehicle Charging Control Algorithm

This study presents the design of an electric vehicle (EV)-grid integration (EVGI) hardware test-bed to implement smart EV charging algorithms. Here, the proposed test-bed also allows to create different grid events via flexible integration of other power hardware (e.g., controllable loads and battery energy storage systems) and test their impacts on EV charging. The design uses a real-time digital simulator to realize a complex distribution grid model with primary and secondary networks. A grid simulator physically realizes the selected nodes of the simulated grid to power an actual EV, forming a hardware-in-the-loop (HIL) test setup. The EV-grid integration is demonstrated based on the custom hardware and software implementation of the J1772 charging protocol using dSPACE MicroLabBox, operating as a custom EV Supply Equipment (EVSE). The HIL test-bed features a novel testing platform for accurate implementation and analysis of scalable charging algorithms. To this end, a data-driven, decentralized, model-free charging controller based on the Additive Increase and Multiplicative Decrease (AIMD) algorithm is presented and validated on an EV using the HIL test-bed. We tested the proposed algorithm under various case studies, and presented a comparison study with an existing droop-based, decentralized charging solution. The results showed that the EV successfully performed charging commands generated by the EVSE and regulated its charging power to effectively reduce the system loading caused by high EV penetration.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2021

Experimental Study of Thermal Management and Enhancement of Adsorption-Based Onboard Hydrogen Storage System

Abstract Although hydrogen has one of the highest specific energies, its energy density in terms of volume is very poor compared to liquid fuels. Thus, to achieve attractive energy density for hydrogen, either high-pressure compression or a storage method is needed. For onboard (vehicles) hydrogen storage, up to 700 bars are needed for commercial fuel cell vehicles. This creates extreme requirements for material strength and thus safety concerns. A new metal-organic framework 5 (MOF-5) was selected as the adsorbent for H2 storage, as it provides promising storage capacity and is commercially available. Under the same H2 storage capacity and tank volume, the adsorption system is expected several folds reduction in pressure. Under the current study, a unique thermal management design using Modular Adsorbing Tank Insert (MATI) is paired with conduction-enhanced compressed MOF-5 beds. Compared to bare beds without conduction enhancement, all beds with conduction enhancement using either aluminum pins or expanded natural graphite (ENG) have shown various levels of improvement on bed thermal response, which can potentially help expedite system charge and discharge cycle times for real applications.

Electrochemistry↗

DOE OSTI2023

Granular Flow in Novel Octet Shape–Based Lattice Frame Material

Moving packed-bed heat exchangers in concentrated solar power (CSP) plants involves heat transfer between heated falling particles and supercritical carbon dioxide. The overall effective thermal conductivity of the moving packed bed and particle-side channel contact resistances are still the bottlenecks in achieving the desirable thermal transport levels. To this end, a novel moving packed bed heat exchanger consisting of an Octet lattice packed between the walls of the particle-side channel is proposed in this study. Granular flow analysis in Octet lattice moving packed bed heat exchanger (OLHX) was conducted through experiments and discrete element method (DEM)-based numerical simulations. The experimental images clearly demonstrated stagnation regions upstream of lattice fibers, void regions downstream of the fiber junctions, and wavy-type unobstructed flow on the lateral sides of the fibers. DEM simulations were successful in capturing all these critical flow phenomena. Larger flow velocities were observed on the lateral sides of the fibers in the simulations. Also, when the particles in the silo were emptied, the final images showed an accumulation of particles on the inter-fiber as well as fiber–channel wall junctions. Moreover, the fiber connections resulted in some regions devoid of particle contact on the channel endwall, which means that these regions would suffer from poor thermal exchange. Finally, the overall mass flowrate increased with increasing porosity for a fixed particle diameter.

14 SOLAR ENERGY↗

DOE OSTI2021

An Investigation of the Effects of Volume Fraction on Drag Coefficient of Non-Spherical Particles Using PR-DNS

Prediction of the drag coefficient is required in gas-particle multiphase flow modeling and simulation. Experimental data and correlations on the fixed-bed system of spherical particles with high volume fractions for various possible arrangements are available in the literature. However, the effect of volume fraction on the drag coefficient of non-spherical particles is not well studied. In solving the momentum equation, the volume fraction plays a vital role in determining the flow resistances. In this paper, we study the impact of volume fraction in the range of 0.069 to 0.65 on the drag coefficient using the computational fluid dynamics (CFD) simulation of air for Reynold number in the range of 10 to 10000 using particle resolved direct numerical solution (PR-DNS). Regular non-spherical particles such as a cube, tetrahedron, and spheroids are used in this study since their single particle’s drag coefficient data are available in the literature for comparison. For this work, the simulations are carried out in the Ansys Fluent using polyhedral mesh, which consumes significantly less computational time and power. The study showed the sphericity and volume fraction have significant impact on the bed pressure drop and average drag coefficient of the particles in the bed especially in high Reynolds number regime. The bed of the spheroid experiences the lowest drag being the most streamlined particle, and the particles with the edges result in a large drag coefficient due to flow separation at the discontinuity. The vector plots verify this behavior where large wake regions are observed behind the tetrahedron particle.

Mahyawansi, Pratik↗

DOE OSTI2022

Characterizing particle-based thermal storage performance using optical methods for use in next generation concentrating solar power plants

Concentrating Solar Power (CSP) generation is an attractive option for low-emission power generation; however, the high costs of thermal storage associated with concentrating solar create a large barrier for their use and adaptation into modern life. Lowering their operation costs, while maintaining high thermal storage and transfer performance is essential. Solid particle-based heat exchange systems can reduce CSP cost but are often less efficient. Efforts to increase their performance have led to use of binary size particle mixes. Presented is an optical-based thermal analysis technique used to measure near-wall thermal conductivity of particle beds essential in determining their heat exchanger efficiency. Modulated Photothermal Radiometry is used to make dynamic temperature measurements, allowing for the extraction of the most relevant thermal properties like thermal conductivity, specific heat, and effusivity. The system uses a modulated laser source causing a damped periodic heat flux, resulting in a frequency and thermal property dependent surface temperature, of which is measured using radiometry. Lock-In techniques are used to extrapolate the amplitude of the signal. Plotting the amplitude against the root angular frequency allows for effusivity measurement by ratio to a known sample. Using specific heat measurements from literature and density measurements, the thermal conductivity of the particle mixes can be calculated. The simplicity of MPTR to probe through the depth of the bed is ideal for use in CSP for dynamic thermal performance monitoring.

Corona, Javier↗

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↗