Transport in Porous Media最新文献

{"title":"Adaptive Time-Stepping and Machine Learning Approach for Pore-Scale Reactive Transport Simulations","authors":"Micha P. Baur,&nbsp;Sergey V. Churakov,&nbsp;Nikolaos I. Prasianakis","doi":"10.1007/s11242-026-02287-4","DOIUrl":"10.1007/s11242-026-02287-4","url":null,"abstract":"<div><p>Pore-scale reactive transport simulations are typically computationally very expensive, which limits their application to complex heterogeneous systems. To solve the computational bottleneck, an adaptive time-stepping algorithm is developed and combined with machine learning-derived surrogate models for geochemical calculations. The time-step is adapted by monitoring the evolution of the diffusion field and the precipitation reactions and by exploiting intermediate stationary states of the system. The algorithm is benchmarked on a system relevant to cement–claystone interaction with the geochemical reaction being the precipitation of Calcium-Silicate-Hydrates (C-S-H) in the pore space of a claystone. The precipitation of C-S-H is modeled as a solid solution, and it is possible to calculate and trace the local Ca/Si ratio of C-S-H, as well as the local amount of gel porosity. In the reactive transport simulations presented here, the geochemical surrogate models alone lead to acceleration factors of up to two orders of magnitude. The adaptive time-stepping algorithm leads to an additional acceleration of three to five orders of magnitude maintaining the relative deviation below 2%. The overall combined acceleration is demonstrated to be six to seven orders of magnitude having a profound impact and opening new computational avenues.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"153 3","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147338591","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":"Penetrative Electrothermal Convection in Anisotropic Porous Media: Role of Velocity and Thermal Boundary Conditions","authors":"B. Arpitha Raju,&nbsp;C. E. Nanjundappa,&nbsp;I. S. Shivakumara","doi":"10.1007/s11242-026-02286-5","DOIUrl":"10.1007/s11242-026-02286-5","url":null,"abstract":"<div><p>This study examines penetrative electrothermal convection driven by uniform internal heat generation in a horizontal layer of dielectric fluid-saturated anisotropic Brinkman porous medium subjected to a uniform AC electric field. Anisotropy in permeability and thermal diffusivity is incorporated, with the vertical permeability fixed at twice the horizontal value, while the thermal anisotropy parameter is varied. The boundaries are taken to be either rigid or stress-free, with an adiabatic lower surface and an upper surface subject to a Robin thermal condition. A linear stability analysis yields a generalized eigenvalue problem, which is solved numerically using a Galerkin method. The principle of exchange of stabilities is shown to hold for rigid–rigid, rigid–free and free–free boundary combinations. Increasing the AC electric Rayleigh number and the Darcy number advances the onset of convection, whereas larger thermal anisotropy and Biot numbers delay instability. Velocity boundary conditions exert a significant quantitative influence, with rigid boundaries being the most stabilizing. As a limiting case, pure electroconvection in the absence of buoyancy is also considered. Although the qualitative trends remain similar to those in the coupled electrothermal regime, substantially higher critical AC electric Rayleigh numbers and wavenumbers are required to trigger instability, reflecting the need for electric forcing alone to overcome viscous and porous dissipation.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"153 3","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146090860","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":"On the Stability of Triple-Diffusive Convection in a Voigt Fluid-Saturated Darcy–Brinkman Porous Medium","authors":"K. R. Raghunatha,&nbsp; Sangamesh,&nbsp;C. Pragathi,&nbsp;I. S. Shivakumara","doi":"10.1007/s11242-026-02289-2","DOIUrl":"10.1007/s11242-026-02289-2","url":null,"abstract":"<div><p>The linear and weakly nonlinear instability of triple-diffusive convection in a Darcy–Brinkman porous layer saturated with a Navier–Stokes–Voigt fluid is investigated analytically. The instability is driven by competing thermal and solutal buoyancy forces in the presence of viscoelastic stresses arising from a Voigt-type strain-rate regularization. Linear instability analysis reveals that the Kelvin–Voigt parameter plays a decisive role in oscillatory convection, exerting either stabilizing or destabilizing effects depending on the solute Darcy–Rayleigh number. A distinctive feature of the system is the emergence, under certain parametric conditions, of a disconnected closed oscillatory instability branch separate from the stationary one. In this regime, the neutral stability boundary is defined by three distinct thermal Darcy–Rayleigh numbers, in sharp contrast to the single threshold that characterizes the corresponding double-diffusive system. A weakly nonlinear stability analysis yields a cubic complex Landau amplitude equation governing the evolution of the disturbance amplitude, suggesting the possibility of subcritical bifurcations under certain parametric conditions. Heat and mass transport are quantified using time- and area-averaged Nusselt and Sherwood numbers, which increase with the Darcy–Prandtl number but decrease with increasing Darcy number and the effective heat capacity group, while variations in the Kelvin–Voigt parameter induce crossover behavior in the transport characteristics. These results clarify the subtle interplay between buoyancy, diffusion, and viscoelastic effects in regulating convection in triply diffusive porous media.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"153 3","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146090929","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":"Time-Resolved X-ray and Neutron Imaging of Brine Percolation and Liquefaction in an Ultra-Soft Sandstone","authors":"Fazel Mirzaei,&nbsp;Mukul Jaiswal,&nbsp;Jessica Ann Zeman,&nbsp;Kim Robert Tekseth,&nbsp;Benoît Cordonnier,&nbsp;Nicolaine Agofack,&nbsp;Pierre Cerasi,&nbsp;François Renard,&nbsp;Basab Chattopadhyay,&nbsp;Dag Werner Breiby","doi":"10.1007/s11242-025-02284-z","DOIUrl":"10.1007/s11242-025-02284-z","url":null,"abstract":"<div><p>The inherent instability and risk of liquefaction associated with unstable glaciomarine clay are both scientifically intriguing and societally important. Stabilizing sensitive soils is increasingly necessary with the unfolding climate crisis, which leads to wilder and wetter weather. The relevant length scales extend from nanometer-sized pores to kilometer-sized geological features. Laboratory experiments are important for gaining a better understanding of these phenomena. Here, we present a study of the ultra-soft sandstone <i>Saltwash South</i>, which essentially consists of quartz particles glued together by clay, as a proxy for nanoscopically fine-grained sensitive soil to facilitate time-resolved imaging of the onset of activation and liquefaction. By employing in situ time-resolved combined X-ray and neutron computed tomography (CT), we visualize the structural deformation caused by the presence of water and salts. While neutron imaging was sensitive to the presence of normal and heavy water, simultaneous X-ray imaging was used to measure the porous structure, swelling, and initial liquefaction response of the consolidated rock. Uniform expansion was observed in regions exposed to water, reflecting clay swelling and disaggregation. In tightly confined samples, the swelling and disaggregation were suppressed. Finally, we discuss future perspectives of this promising approach to studying liquefaction phenomena in porous media.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"153 3","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11242-025-02284-z.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146090959","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":"SLM-Fabricated CuCrZr TPMS Heat Exchangers: Influence of Structural Design on Fluid Transport Characteristics","authors":"Xiaoqing Yang,&nbsp;Yaqin Yang,&nbsp;Chao Suo,&nbsp;Bin Liu,&nbsp;Xueli Chen,&nbsp;Zezhou Kuai,&nbsp;Yaoxian Zhang,&nbsp;Yuan Gao,&nbsp;Zhonghua Li,&nbsp;Rui Wang","doi":"10.1007/s11242-026-02288-3","DOIUrl":"10.1007/s11242-026-02288-3","url":null,"abstract":"<div><p>The hydrodynamic performance of porous structures is paramount for the efficiency of compact heat exchangers. This study investigates the fluid transport characteristics of three types of triply periodic minimal surface (TPMS) structures (Diamond, Primitive, and IWP) fabricated from CuCrZr alloy via selective laser melting (SLM). Combining computational fluid dynamics (CFD) simulations with experimental validation, we systematically elucidate the influence of structure and unit cell density on key performance metrics, including pressure drop, friction factor, and permeability. Results reveal that the Diamond structure, with the lowest inlet porosity (23.07%), exhibits the most complex flow patterns—characterized by significant flow separation, enhanced vorticity, and superior flow uniformity—leading to the highest pressure drop and friction factor, but the lowest permeability. In contrast, the primitive and IWP structures, with higher inlet porosities (46.34% and 50.07%, respectively), offer lower flow resistance. Furthermore, increasing the unit cell count for a fixed overall porosity amplifies these effects, resulting in a higher pressure drop per unit length and reduced permeability due to more frequent flow disturbances. The experimental measurements of pressure drop align well with numerical predictions, with discrepancies attributed to a reduction in effective functional porosity caused by adherent unmelted particles during the SLM process. This work provides fundamental insights into the structure-transport relationships in TPMS architectures and offers practical guidelines for designing high-performance heat exchangers for applications under extreme conditions, such as deep-sea thermal management.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"153 3","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146090930","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":"Stability of Double-diffusive Mixed Convection in a Vertical Pipe Filled by Porous Medium","authors":"Saurabh Kapoor,&nbsp;Durgaprasad Nayak","doi":"10.1007/s11242-025-02276-z","DOIUrl":"10.1007/s11242-025-02276-z","url":null,"abstract":"<div><p>Thermal gradients in many natural environments, such as undersea hot springs and hydrothermal vents, are often accompanied by variations in the concentration of chemical compounds carried into the surrounding seawater. This interaction gives rise to a phenomenon known as thermosolutal mixed convection. To investigate the underlying physical mechanisms in such systems, a vertical pipe filled with saline water through a porous medium is considered. The non-Darcy Brinkman–Forchheimer extended model is employed, assuming that the saturated porous medium is in a thermal equilibrium state. The instability boundary curve reveals three distinct regimes: (i) the Rayleigh–Taylor (R–T) phenomenon, (ii) a nonlinear variation of mixed convection on a log–log scale and (iii) a linear variation on a log–log scale in the pipe. These regimes occur at specific Reynolds numbers. The instability mechanism of the base flow is examined using the spectral collocation technique. The present study focuses on the instability of a dense porous medium, where the Darcy–Forchheimer drag is incorporated in the momentum equation. A linear stability analysis is performed for a wide range of Darcy numbers (<span>(10^{-8})</span> to <span>(10^{-5})</span>). Similar to cross-diffusive free convection in a vertical slot, a relationship is observed between the solutal Rayleigh number (<span>(Ra_S)</span>) and the thermal Rayleigh number (<span>(Ra_T)</span>). A hyperbolic relationship between the threshold values of <span>(Ra_S)</span> and <i>Da</i> is found at <span>(Re = 1000)</span> for the upper limit of <span>(Ra_S)</span> in the first regime, expressed as <span>(Ra_S)</span> <i>Da</i> =4.9 x <span>(10^{-4})</span>. Simulations of secondary flow profiles are also presented at the critical state for all three regimes.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"153 3","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11242-025-02276-z.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146090931","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":"Effects of Near-surface Winds and Soil Pore Air Pressure Fluctuations on Intra-soil Gas Transport Mechanisms","authors":"Junguo Hu,&nbsp;Taolve Chen,&nbsp;Lingxia Feng,&nbsp;Junjie Jiang,&nbsp;Chao Zhu,&nbsp;Wei Jiang","doi":"10.1007/s11242-026-02292-7","DOIUrl":"10.1007/s11242-026-02292-7","url":null,"abstract":"<div><p>Soil carbon flux is the main source of atmospheric CO₂. An accurate understanding of the soil CO₂ transport process is of substantial significance for the study of global carbon flux and the prediction of future climate change. In this study, data on CO₂ concentrations over time in five media under the influence of surface winds and pore air pressure fluctuations were analysed using (1) an analytical solution of the dispersive transport model and (2) a numerical solution of the advection transport model to accurately delineate the transport mechanisms that control soil gas transport. For loam and sandy soils, gas transport is mainly dominated by diffusion, and in highly permeable sand, gas transport changes from diffusion–dispersion transport to advection transport dominated by wind speed and pore air pressure fluctuations. The effects of surface wind and pore pressure fluctuations on advection velocity differed with depth. The pore pressure fluctuations had a stronger penetrating power. The ratio of the enhanced diffusion coefficient (<i>D</i>_<i>e</i>) to advection velocity(<i>v</i>) was effective in determining the transport mechanism driving gas transport within the medium, which is conducive to accurately modelling the release of CO₂ from soils to improve the understanding and assessment of the global carbon cycle.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"153 3","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146090960","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":"Backward Similarity Solution of the Boussinesq Groundwater Equation","authors":"Shuntaro Togo,&nbsp;Koichi Unami","doi":"10.1007/s11242-025-02285-y","DOIUrl":"10.1007/s11242-025-02285-y","url":null,"abstract":"<div><p>Groundwater flow in an unconfined aquifer resting on a horizontal impermeable layer with a boundary condition of a rapid increase in the source water level is considered in this work. The newly introduced condition, referred to as the backward power law head condition, represents a situation where the water level in the adjacent water body increases more rapidly than do conventional problems, which can only represent a situation akin to a traveling wave under rising water level conditions, given its consideration of infinite time. This problem admits the similarity transformation which allows the nonlinear partial differential equation to be converted into a nonlinear ordinary differential equation via the Boltzmann transformation. The reduced boundary value problem is interpreted as the initial value problem for a system of ordinary differential equations, which can be numerically solved via Shampine’s method. The numerical solutions are in good agreement with Kalashinikov’s special solution, which is also introduced into the Boussinesq equation. We find that the solution is consistent with the limit of the forward power law head condition. The new approximate analytical solution and the associated wetting front position are derived by assuming that the solution has the form of quadratic polynomials, which enables the description of the time progression of a real front position. The obtained approximation is compared to Shampine’s solution to check the accuracy. Furthermore, the finite element method is applied to the original partial differential equation, which validates Shampine’s solution for use as a benchmark.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"153 3","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11242-025-02285-y.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146090961","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":"Analysis and Comparison of Natural Shear and Induced Tensile Fractures for Caprock Leakage Assessment","authors":"Sahyuo Achuo Dze,&nbsp;Tomos Phillips,&nbsp;Reza Najafi-Silab,&nbsp;Sarah Perez,&nbsp;Tom Bultreys,&nbsp;Vladimir Novak,&nbsp;Christian M. Schlepütz,&nbsp;Veerle Cnudde,&nbsp;Florian Doster,&nbsp;Kamaljit Singh,&nbsp;Kevin Bisdom,&nbsp;Andreas Busch","doi":"10.1007/s11242-025-02283-0","DOIUrl":"10.1007/s11242-025-02283-0","url":null,"abstract":"<p>The comparison between laboratory-induced and subsurface fractures, and their corresponding flow is still unclear. Here, we examine three natural shear fractures and two induced tensile fractures from the same low-permeability lithology. Using high-resolution synchrotron imaging, we extracted three-dimensional fracture void geometries to analyse aperture distributions, surface roughness, and spatial correlation patterns. We then compared measured fracture transmissivities against theoretical predictions from parallel-plate models (cubic law) and direct numerical simulations (DNS) of flow to evaluate consistency and uncertainties. We find that despite differences in heterogeneity, induced tensile and natural shear fractures can generally yield similar flow in caprocks. Our comparisons further indicate that the choice of flow estimation method can introduce more uncertainty than the fracture opening mode.</p>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"153 2","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11242-025-02283-0.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145983244","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":"Comparative Study of Three Fine Grains Infiltration Techniques to Re–Fill Granular Materials Eroded by Suffusion","authors":"Abhijit Hegde,&nbsp;Nadia Benahmed,&nbsp;Antoine Wautier,&nbsp;Pierre Philippe","doi":"10.1007/s11242-025-02271-4","DOIUrl":"10.1007/s11242-025-02271-4","url":null,"abstract":"<div><p>With the aim of remediating soils eroded by suffusion, three different techniques were explored and compared at the laboratory scale in this study to effectively infiltrate fine grains into a porous granular filter. In the first technique, which serves as a reference owing to its simple implementation in the field conditions, a pure layer of fine grains is deposited on the upstream side of the filter. However, this technique is found inefficient due to a collective effects between the fine grains which induces formation of surface cake, preventing thus their penetration into the filter. A second technique is therefore proposed in which the fine grains are mixed homogeneously with coarse particles and deposited on the upstream of the filter. Besides, a third technique was designed, where the fine grains are put in suspension in the water seeping through the filter in order to limit their collective clogging. The influence of size ratio and fine grains shape for both the mixture and suspension techniques, and fine grains shape in the case of suspension technique were also explored. For the latter, the cumulative efficiency was further investigated through multiple rounds of infiltration. The results showed that, for a given size ratio between the coarse and fine grains, there is a significant increase in the amount of infiltrated fine grains for the mixture and suspension techniques, that is respectively twice and eight times larger than the value obtained with the reference one, respectively. In the mixture technique, fine grains infiltration is found to increase with the size ratio up to a limiting value whereas, in the suspension method, it decreases as the size ratio increases. Moreover, rounded fine grains are found to infiltrate more easily than angular ones. This work should contribute to improve engineering protocols for achieving deep infiltration through the selection of an appropriate fine grains infiltration technique, with the objective of developing an effective remediation method.</p></div>","PeriodicalId":804,"journal":{"name":"Transport in Porous Media","volume":"153 2","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145983249","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}