Transport in Porous Media最新文献

{"title":"A Semi-Analytical Model to Simulate Fluid Flow in Fractured Reservoirs","authors":"V. L. Teja-Juárez,&nbsp;J. A. Torres-López,&nbsp;R. M. Mariscal-Romero,&nbsp;R. G. Camacho-Velázquez,&nbsp;A. P. Gómora-Figueroa,&nbsp;D. E. Solano-Rojas","doi":"10.1007/s11242-025-02203-2","DOIUrl":"10.1007/s11242-025-02203-2","url":null,"abstract":"<div><p>A semi-analytical simulator has been developed to describe single-phase gas flow in hydraulically fractured wells hosted by naturally fractured rock. The governing formulation embeds Streltsova’s unsteady-state matrix-transfer function to fracture directly in a dual-porosity finite difference framework and advances the solution with an implicit–explicit scheme. This strategy preserves the analytical description of transient exchange, while avoiding the convolution integrals and tight time-step controls that burden fully numerical unsteady-state models. Both pseudo steady-state and unsteady-state transfer functions can be invoked with different shape factor and we use the logarithmic spaced, local grid refinement methodology to consider the hydraulic fracturing effects that occur near well at minimal computational cost. The simulator was tested on homogeneous, discretely segmented, and fully heterogeneous synthetic reservoirs under constant rate and constant pressure production. In every case, the unsteady-state option reproduced the characteristic early-time pressure transients and matrix-fracture flow that the classical Warren–Root’s model with the pseudo-transitory transfer function may not capture. Compared with a published hybrid numerical–analytical workflow for same shape factor, the new implementation is faster than twelve times, while presenting good pressure behavior and maintaining numerical precision. Because the transfer function remains analytical, computational effort scales only with fracture count and not with matrix grid density, making the approach well suited to large, heterogeneous unconventional plays as well as geothermal and subsurface storage applications. The tests confirm that the method achieves the desired balance between physical fidelity and computational efficiency, providing a robust tool for transient flow analysis and production-strategy optimization in fractured reservoirs.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"152 9","pages":""},"PeriodicalIF":2.6,"publicationDate":"2025-07-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11242-025-02203-2.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145165504","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"The Research and Impact of the Works by Don Nield","authors":"Sid Becker,&nbsp;David Dempsey,&nbsp;Dan Zhao","doi":"10.1007/s11242-025-02199-9","DOIUrl":"10.1007/s11242-025-02199-9","url":null,"abstract":"<div><p>This is a commemoration of the work of Don Nield. His theoretical and applied mathematical models have had a tremendous impact in many fields including convection heat and mass transfer, stability analysis, nanofluid enhanced convection, and flows in porous media.  Here the major contributions of Don Nield’s Scopus-indexed outputs are reviewed. His research outputs and related impact are presented chronologically. The citation counts are presented with respect to each work published in that year. The work and impact of Don Nield as single author and in his collaborations are presented.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"152 9","pages":""},"PeriodicalIF":2.6,"publicationDate":"2025-07-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11242-025-02199-9.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145165503","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Assessing the Influence of Fluid Properties on Darcian Flow Dynamics in Porous Media: A Detailed Study of Liquids and Gases","authors":"Asseel M. Rasheed Al-Gaheeshi,&nbsp;Farhan Lafta Rashid,&nbsp;Mudhar A. Al-Obaidi,&nbsp;Atef Chibani,&nbsp;Moustafa Boukraa,&nbsp;Tawfiq Chekifi","doi":"10.1007/s11242-025-02201-4","DOIUrl":"10.1007/s11242-025-02201-4","url":null,"abstract":"<div><p>The study intends to investigate the effects of fluid characteristics on the Darcian flow through porous media, which would help to perceive the merits of its employment in different designs of porous media systems. To systematically conduct this aim, both theoretical modeling and numerical simulations are used while using different incompressible fluids of four liquids (water, oil, n-heptane, and n-hexane) and three gases (nitrogen, carbon dioxide, and methane). In this regard, an acceptance between the simulation results and the experimental findings is confirmed. The results demonstrate the actual effect of friction factor and pressure gradient on the fluid properties of viscosity and density. An increase in Reynolds numbers caused a decrease in the friction factor, which signify a shift in dominance from viscous to inertial forces. The pressure gradient rises with higher inlet velocities, proposing that a higher flow rate necessitates a larger pressure differential to overwhelmed resistance of the porous structure.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"152 9","pages":""},"PeriodicalIF":2.6,"publicationDate":"2025-07-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145165502","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Instabilities and Bifurcations in Turbulent Porous Media Flow","authors":"Vishal Srikanth,&nbsp;Andrey V. Kuznetsov","doi":"10.1007/s11242-025-02195-z","DOIUrl":"10.1007/s11242-025-02195-z","url":null,"abstract":"<p>Microscale turbulent flow in porous media is conducive to the development of flow instabilities due to strong vortical and shearing flow occurring within the pore space. When the flow instabilities around individual solid obstacles interact with numerous others within the porous medium, unique symmetry-breaking phenomena emerge as a result. This paper focuses on investigations of the vortex dynamics and flow instabilities behind solid obstacles in porous media, emphasizing how solid obstacle geometry and porosity influence both microscale and macroscale flow behavior. Two distinct symmetry-breaking mechanisms were identified in different porosity ranges. In low porosity media (&lt; 0.8), a “deviatory flow” phenomenon occurs, where the macroscale flow deviates from the direction of applied pressure gradient at Reynolds numbers above 500. Deviatory flow is a source of macroscale Reynolds stress anisotropy, which is counterbalanced by a diminished vortex core size. In the intermediate porosity regime (0.8–0.95), a “jetting flow” mechanism creates asymmetric microscale velocity channels in the pore space through temporally biased vortex shedding, occurring during the transition to turbulence. Both symmetry-breaking phenomena are critically influenced by solid obstacle shape, porosity, and Reynolds number. Circularity of solid obstacle geometry and an adequately high-Reynolds number provide critical conditions for symmetry-breaking, whereas porosity can be used to parametrize the degree of symmetry-breaking. This paper provides fundamental insights into the intricate flow dynamics in porous media, offering a comprehensive understanding of how microscale vortex interactions generate macroscale flow asymmetries across different geometric configurations.</p>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"152 8","pages":""},"PeriodicalIF":2.6,"publicationDate":"2025-07-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11242-025-02195-z.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145161382","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Reactive Transport and Péclet Number Analysis of Hydrogen Flux Pathways in Uniform Clay Matrix: Implications for Underground Storage","authors":"Samuel Bowman,&nbsp;Arkajyoti Pathak,&nbsp;Shikha Sharma","doi":"10.1007/s11242-025-02200-5","DOIUrl":"10.1007/s11242-025-02200-5","url":null,"abstract":"<div><p>Reactive transport modeling in unique chamosite, illite, and kaolinite monomineralic clay reservoirs was performed to predict spatial and mineralogical control on H_2(aq) distribution and flow dynamics at simulated underground hydrogen storage pressure and temperature conditions in high permeability depleted oil reservoirs and lower permeability sandstones. Results are normalized against simulations performed in pure quartz for comparative purposes and to serve as a benchmark for H₂ storage in sandstone or quartz-arenite reservoirs. Isobaric and isothermal condition results indicate that Péclet numbers increase with volume rate of injected fluid. Advection-mediated transport is ubiquitous for injection rates between 0.01 and 1.0 L per second, although diffusion-mediated transport is prevalent in the center of the RMT block at 0.001 L per second injection. Source/sink term and Péclet number analysis indicate that the effect of mineralogy on H_2(aq) transport is small. Flow velocities in kaolinite are typically the fastest, but chamosite Péclet numbers are greatest. This suggests that kaolinite favors diffusion, chamosite favors advection, and illite is intermediate. A more accurate reflection of underground hydrogen storage conditions incorporating temperature and pressure gradients, permeability anisotropy, and mineralogical heterogeneity shows a decrease in Péclet numbers proportional to distance from injection well. Thus, along the reservoir-caprock boundary and in the absence of cushion gas, H_2(aq) loss through diffusion is probable. Although quantitative flow regime analysis cannot be determined at many grid point locations due to uniform H_2(aq) concentration, these locations are very likely diffusion dominant.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"152 8","pages":""},"PeriodicalIF":2.6,"publicationDate":"2025-07-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145161381","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Macroscopic Orientation of Inertial Flows in Porous Media","authors":"Yanis Bendali,&nbsp;Morgan Chabanon,&nbsp;Quentin Holka,&nbsp;Benoît Goyeau","doi":"10.1007/s11242-025-02186-0","DOIUrl":"10.1007/s11242-025-02186-0","url":null,"abstract":"<div><p>Macroscopic models of inertial flows in porous media have many practical applications where direct numerical simulations are not feasible. The Forchheimer equation describes macroscopic momentum transport accounting for inertial effects at the pore scale through a nonlinear correction tensor <span>({textbf{F}}_beta)</span> to the permeability. The goal of this work is to study the effects of inertial flow orientation on the Forchheimer correction. Using up-scaling approaches such as the volume averaging method, <span>({textbf{F}}_beta)</span> can be determined. However, the procedure requires to deal with a nonlinear problem for the deviations of the local velocity field. This is commonly tackled by assuming that the inertial convective velocity is decoupled from the velocity deviations. Here, we propose an alternative approach based on regular perturbation expansion leading to a series of linear closure problems. The values of <span>({textbf{F}}_beta)</span> predicted by both approaches are compared for various values of the Reynolds number and flow orientation. Compared to the local inertial–convection approach, the proposed linearized closure problem has the advantage of being self-consistent, independent of the pore Reynolds number and of flow orientation. It is, however, limited in validity by Reynolds number below one and requires the solution of closure problems of higher dimensions. Then, macroscopic simulations are performed to evaluate the importance of varying pressure gradient orientation on the macroscopic inertial flow. Numerical results of the general macroscopic model obtained by the volume averaging method highlight the necessity to account for extra-diagonal terms as well as macroscopic gradient orientation in the determination of the Forchheimer tensor.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"152 8","pages":""},"PeriodicalIF":2.6,"publicationDate":"2025-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11242-025-02186-0.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145161225","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Convective Heat Transfer in Brinkman–Darcy–Kelvin–Voigt Fluid with Variable Gravity and Generalized Maxwell–Cattaneo Law","authors":"Amit Mahajan,&nbsp;Saravanan P","doi":"10.1007/s11242-025-02194-0","DOIUrl":"10.1007/s11242-025-02194-0","url":null,"abstract":"<div><p>This article investigates thermal convection in Kelvin–Voigt fluids saturating a Brinkman–Darcy-type porous medium under variable gravity effects. The stability analysis encompasses linear, nonlinear conditional, and nonlinear unconditional regimes, leveraging the generalized Maxwell–Cattaneo law for heat flux. Normal mode technique is applied for linear analysis, while nonlinear governing equations for conditional and unconditional stability are derived using energy methods. The compound matrix method is utilized to compute critical Rayleigh numbers and corresponding wave numbers. Numerical computations are performed in MATLAB to determine all critical values, with graphical results illustrating stability trends. Six gravity variation profiles are examined, revealing that the gravity modulation parameter <span>(epsilon)</span> influences stability depending on the gravity profile, either promoting or suppressing convection. Additionally, the parameter <span>(xi)</span> consistently enhances stability by increasing the critical Rayleigh numbers across all cases. These findings highlight the role of gravity variations and material parameters in shaping convection onset, contributing to a deeper understanding of thermal stability in viscoelastic porous systems.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"152 8","pages":""},"PeriodicalIF":2.6,"publicationDate":"2025-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145161156","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Self-Pumping Transpiration Cooling: A Joint Experimental and Numerical Study","authors":"Julian Härter,&nbsp;Maziar Veyskarami,&nbsp;Martin Schneider,&nbsp;Johannes C. Müller,&nbsp;Hanchuan Wu,&nbsp;Rainer Helmig,&nbsp;Bernhard Weigand,&nbsp;Grazia Lamanna,&nbsp;Rico Poser","doi":"10.1007/s11242-025-02198-w","DOIUrl":"10.1007/s11242-025-02198-w","url":null,"abstract":"<div><p>A joint experimental and numerical study is presented to close the current gap in fully coupled data and modeling capabilities for self-pumping transpiration cooling (SPTC). An experimental setup was developed to investigate the effects of the porous medium properties, the flow conditions, and the interactions between solid and coolant on SPTC. Additionally, a two-reference-point, locally emissivity-corrected evaluation methodology for analyzing infrared (IR) measurements was developed, which is valid for quasi-steady evaporation regimes and achieves a better repeatability. For the numerical simulations, we developed an upscaling workflow with pore-network models derived from micro computed tomography (CT) data to accurately describe effective representative elementary volume (REV)-scale parameters and relations. Using upscaled properties, we created a non-isothermal, two-phase Darcy-scale model for the porous medium and modeled free-flow with Reynolds-averaged Navier–Stokes equations, employing an shear stress transport (SST) <span>(ktext {-}omega)</span> turbulence closure to capture near-wall shear stress effects. Coupling conditions ensured mass, momentum, and energy transfer at the interface. The experimental results show a high reproducibility and new insights for the surface temperature at SPTC with the new IR method. The comparison between experimental and numerical results show good agreements. The developed simulation workflow is a major step toward creating a digital twin of an experimental SPTC system. This work lays the foundation for investigating the influence of parameters on SPTC systems and optimizing their efficiency.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"152 8","pages":""},"PeriodicalIF":2.6,"publicationDate":"2025-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11242-025-02198-w.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145161272","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Buoyancy-Induced Convection in a Partially Porous Chamber Filled with a Fluid of Variable Viscosity Under Periodically Heated Source","authors":"M. S. Astanina,&nbsp;M. A. Sheremet","doi":"10.1007/s11242-025-02191-3","DOIUrl":"10.1007/s11242-025-02191-3","url":null,"abstract":"<div><p>This study deals with the thermal gravitational convection of a variable viscosity fluid in a closed two-dimensional cavity in the presence of a porous layer. In addition, at the lower edge of the region there is a heat-conducting energy source with periodic heat release. The problem has been solved in dimensionless non-primitive parameters including stream function and vorticity using the finite difference schemes on a uniform structured mesh. The influence of the defining characteristics of heat transfer such as the Darcy and Rayleigh numbers, the height of the porous layer, the heater location, as well as the frequency of the heat dissipation from the heater, has been studied. The results have been demonstrated using isolines of the stream function and temperature combined with time dependences of the integral characteristics of heat transport for various values of dimensionless parameters. The conclusions obtained have demonstrated that the periodic nature of heat dissipation completely determines the periodicity of convective processes in the cavity. In addition, the presence of a porous layer and the ability to change the position of the heater make it possible to control heat removal.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"152 8","pages":""},"PeriodicalIF":2.6,"publicationDate":"2025-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145161155","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Improving Numerical Fluid Flow Simulation by Ring Artifact Removal in Micro-CT Images of Porous Media Using Attention Autoencoder–Decoders","authors":"Mehdi Mahdaviara,&nbsp;Maryam Mousavi,&nbsp;Yousef Rafiei,&nbsp;Amir Raoof,&nbsp;Mohammad Sharifi","doi":"10.1007/s11242-025-02190-4","DOIUrl":"10.1007/s11242-025-02190-4","url":null,"abstract":"<div><p>The emergence of micro-computed tomography has significantly enhanced our ability to examine the morphology of porous materials and the dynamics of fluid flow within pore spaces. However, image-based analyses can be compromised by various artifacts, particularly ring artifacts, which appear as concentric rings in the images. These artifacts can be misinterpreted as part of the pore space, artificially connecting pores and thus influencing numerical simulations. This study examines the influence of ring artifacts on pore network modeling (PNM), direct numerical simulation (DNS), and prominent numerical techniques, and presents a computing approach for their effective mitigation. For this purpose, a dataset was compiled from the literature that includes the images of Fontainebleau, Boise, and Belgian sandstones. Data augmentation was implemented by extracting real ring patterns from Fontainebleau samples and superimposing them onto clean images of the sandstones. Two U-Net autoencoder architectures (base and Attention U-Net) were trained for a regression task aimed at removing ring artifacts while reconstructing the underlying pore morphologies. The Attention U-Net outperformed the base model, achieving a mean squared error of 0.07 (calculated based on the grayscale values between 0 and 255). Visual evaluations confirmed the model’s effectiveness in artifact removal and pore morphology reconstruction. The model was further tested on unseen pore-scale data containing real ring artifacts, which indicated a high performance in removing the artifacts. DNS and PNM were performed on both original (with real rings) and improved 3D samples (200³ voxels) to assess the impact of artifact removal on transport properties. The results revealed that ring artifacts, identified as flow pathways, significantly influence the velocity profiles. While the presence of the artifact had a minimal effect on porosity (a 1.68% error) and the number of pores (1.45% error), it significantly increased the permeability by 34%.</p><h3>Graphical Abstract</h3>\u0000<div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"152 8","pages":""},"PeriodicalIF":2.6,"publicationDate":"2025-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11242-025-02190-4.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145161154","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}