Porosity Of Medium Research Articles

Algebraic growth of 2D optimal perturbation of a plane Poiseuille flow in a Brinkman porous medium

The impact of various factors on the stability of a plane porous Poiseuille flow is investigated in detail. Both modal and non-modal linear stability analyses are used to study the effects of the porosity of the porous media and the ratio of effective viscosity to the fluid viscosity. The present analysis includes solving the linearized Navier-Stokes equation in the form of the Orr-Sommerfeld (O-S) type equation by applying the 2D perturbation to the basic mean flow. The Chebyshev collocation method is used to solve the Orr-Sommerfeld equation numerically. Through modal analysis, the accurate values of the critical triplets (αc,Rc,cc), the eigen-spectrum, the growth rate curves, and the marginal stability curves are studied. Then, by using non-modal analysis, the transient energy growth G(t) of two-dimensional optimal perturbations, the ε-pseudospectrum of the non-normal O-S operator (L), and the regions of stability, instability, and potential instability of the fluid flow system are investigated in detail. The collective results of both modal and non-modal analysis show that the porous parameter and the ratio of effective viscosity to the fluid viscosity have stabilizing effects on the fluid system due to the postponement of the onset of stability.

Thermal investigation and parametric analysis of cascaded latent heat storage system enhanced by porous media

Cascaded latent heat storage (CLHS) has become a promising solution to the recovery and utilization of thermal energy. Many studies optimized the inlet temperature and velocity of heat transfer fluid to improve the heat transfer rate for the charging process of CLHS systems with pure phase change materials (PCMs), but the effect is not significant, mainly because the thermal resistance of the system comes from the PCMs side due to the low thermal conductivity of pure PCMs. Porous media composite phase change materials (CPCMs) have been shown to have promising prospects in a single PCM system, but no research precedent has been found for the application of porous media CPCM in the CLHS system. Besides, when CPCMs are encapsulated in the cavities of the CLHS system instead of pure PCMs, the difference in the degree of impact between the porosity of a porous media and boundary conditions on the thermal performance of a CLHS system is still unclear. The purpose of this article is to significantly improve the heat transfer rate by applying the CPCM to the CLHS system and to comprehensively compare the impact of porosity and HTF inlet boundary conditions on the charging/discharging time and energy/exergy charged/discharged rates of a CLHS system with CPCMs through two sensitivity analysis methods. When the porosity reduces from 0.97 to 0.92, the charging time and the discharging time decrease by 74.50 % and 73.39 % respectively; the average energy charge rate and average exergy charge rate increase by about 4.11 times; the average energy and exergy discharge rates increase by 5.76 times and 5.8 times respectively. Sensitivity analysis shows that the porosity of a porous media is the most important factor affecting the charging and discharging rates of a CLHS system, but its impact on energy and exergy charge/discharge rates is not monotonic. The effect of the inlet temperature and velocity of the HTF on a CLHS system is complex, and simply changing a single variable may not be effective in optimizing thermal performance parameters. This study provides a comprehensive performance enhancement strategy for the thermal energy storage and utilization processes of CLHS systems, which have great potential in fields such as solar energy collection, waste heat recovery, and building heating.

Some remarks on simplified double porosity model of immiscible incompressible two-phase flow

The paper is devoted to the derivation, by linearization, of simplified homogenized models of an immiscible incompressible two-phase flow in double porosity media in the case of thin fissures. In a simplified double porosity model derived previously by the authors the matrix-fracture source term is approximated by a convolution type source term. This approach enables to exclude the cell problem, in form of the imbibition equation, from the global double porosity model. In this paper we propose a new linear version of the imbibition equation which leads to a new simplified double porosity model. We also present numerical simulations which show that the matrix-fracture exchange term based on this new linearization procedure gives a better approximation of the exact one than the corresponding exchange term obtained earlier by the authors.

A POROUS MEDIA APPROACH FOR NUMERICAL OPTIMISATION OF THERMAL WHEEL

The experimental investigations of rotating heat exchangers are usually too costly and provide limited understanding for the phenomena of heat and fluid flow within them; hence, a less expensive and more comprehensive method is required to investigate what can affect their overall performance. In the current study, a porous media concept is presented as an alternative way to numerically analyse the fluid flow and heat transport through a rotary thermal regenerator. An aluminum core formed of multi-packed square passages is simulated as a porous medium of an orthotropic porosity in order to allow the counter-flowing streams to flow in a way similar to that inside the regenerator core. The geometric properties of the core were transformed into the conventional porous media parameters such as the permeability and inertial coefficient based on empirical equations; so, the core has been dealt with as a porous medium of known features. Fluid and solid phases are assumed to be in a local thermal non-equilibrium state with each other. A commercial CFD code "STAR CCM+" was used to solve the current problem numerically, where heat is allowed to be exchanged between the two phases and tracked by creating a heat exchanger interface in the core region. The results are presented by means of overall thermal effectiveness, pressure drop, and coefficient of performance COP. Using porous media approach has been found to be sufficient to simulate the current problem, where the currently computed data were found to deviate by up to 2.7% only from the corresponding analytical and experimental data. The data obtained reveal an obvious impact of the core geometrical parameters on both the heat restored and pressure loss; and hence, the overall efficiency of the regenerator system.

Impact of rotation on thermal instability of Darcy–Brinkman porous layer filled with a Jeffrey nanofluid

The novelty of this research lies in its investigation of the initiation of thermal instability within a Darcy–Brinkman porous layer containing a Jeffrey nanofluid under the influence of rotation. While non-Newtonian fluids complex flow patterns and heat transfer enhancements due to thermal gradients are known, this study addresses a specific and previously unexplored aspect, the impact of rotation on Darcy–Brinkman porous layer containing a Jeffrey nanofluid. When a fluid is subjected to rotation, it can induce additional forces and flow patterns within the fluid. The linear stability theory and normal mode analysis method are used to analyze the stability of the system analytically and are computed numerically using the software Mathematica version 11.3. The investigation focuses on understanding the influence of various factors, namely, the Taylor number, the Lewis number, the Jeffrey parameter, the modified diffusivity ratio, nanoparticles Rayleigh number, the Darcy–Brinkman number, and medium porosity on the onset of thermal instability in the physical system. This novel perspective can contribute to a deeper understanding of heat transfer processes in non-Newtonian fluids and has potential applications in improving the efficiency of industrial systems, such as cooling systems, heat exchangers, and microfluidic devices where efficient heat transfer is crucial for performance improvement.

Implementation and validation of effective thermal conductivity model of coal‐rock mass with rough capillary pore structure under CT three‐dimensional reconstruction

AbstractOn the strength of the capillary physical equivalent model and heat transfer law of series‐parallel combined structural unit, the calculation model of effective thermal conductivity of coal‐rock mass is obtained. The model considers the effect of roughness, porosity factor and other factors on effective thermal conductivity. CT digital core technology is applied to reconstruct the 3D model of coal‐rock mass, and one‐way heat transfer numerical simulation is carried out according to the model. It indicates that the internal temperature of coal‐rock mass is non‐uniformly distributed under the influence of pore structure, and the isothermal surface of pore region bulges towards the boundary of heat flux. Comparing the calculated and simulated values under the same conditions, the average error percentage of the theoretical calculation model is 6.95%, and the model has outstanding prediction performance in porous media with medium porosity and regular pore distribution, which provides a reference for predicting the effective thermal conductivity of coal‐rock mass.

Surface motion of P/SV wave scattering by a semicircular canyon in saturated half-space by complex function and conformal mapping

Based on the complex function theory of elastic dynamics, a complex function solution of P/SV-wave scattering by a semicircular canyon in an elastic porous saturated half-space is obtained. Firstly, through Biot theory and Helmholtz decomposition, the wave potential functions with unknown coefficients in the model are obtained based on the wave function expansion method. Secondly, the complex function expressions of stress and displacement in the model are obtained by using the complex function method and conformal mapping. According to the boundary conditions of the horizontal ground surface and semicircular canyon surface, this boundary value problem can be converted into a series of algebraic equations and solved numerically by truncation. Then, the correctness of the complex function solution in this paper is verified by comparing it with the published results through degenerating the scattering problem in the saturated porous elastic half-space in this paper into that in the elastic half-space. Finally, the influence of the boundary drainage conditions, the medium porosity, and the incident angle and frequency of the incident wave on the ground surface displacement is investigated. Although the P/SV-wave scattering problem is a traditional problem, there is currently no accepted analytical solution that strictly satisfies the zero-stress condition on the horizontal ground surface. This paper is an attempt at this problem, which may provide insight into the search for the strictly analytical solution for P/SV-wave scattering.

Homogenization of fluid saturated double porosity media with a new type of the contrast in the Biot mesoscopic model

This paper presents a mathematical model of a double poro-elastic medium derived by the homogenization of a two-component periodically heterogeneous Biot continuum characterized by strong heterogeneities in the permeability and poroelastic coefficients of mesoscopic structure. A new scaling of the mesoscopic material parameters with respect to the small parameter which is involved in the asymptotic analysis is introduced to capture this high contrast quality of the mesoscopic model. It is shown that such a scaling ansatz of the mesoscopic material parameters can be justified using different micromodels associated with the two mesoscopic components. While the “matrix” is a little permeable stiff phase, the “conductive channels” are made of a very soft fibrous structure equivalently represented by an network of helical springs. The unfolding method of the homogenization is employed to derive the macroscopic model involving the frequency-dependent effective parameters. Semipermeable interfaces between hard dual porosity and soft primary porosity are considered. The wave dispersion of two shear wave modes S1,S2, and two pressure wave modes P1,P2, is illustrated in an numerical example and validated, in a part, using the reference dispersion analysis based on the Bloch wave decomposition.

Significance role of dual porosity and interfacial nanolayer mechanisms on hybrid nanofluids flow: A symmetry flow model

The investigation of heat and mass transfer between two porous disks containing porous media is the focus of this paper. In this research, (Ag–Al2O3/water)-based hybrid nanofluid is used. Governing partial differential equations are converted to nonlinear ordinary differential equations by similarity transformations. Numerical solutions to nonlinear systems of equations are achieved using the shooting technique. It is examined how the particle’s diameter and nanolayer thickness affect the thermal conductivity of a hybrid nanofluid. Due to the increment of nanolayer thickness, the Nusselt number increases and vice versa for the radius of the particle. Tables and graphs show how the governing parameters influence velocity, temperature, and mass concentration profiles, as well as physical quantities such as skin friction coefficient and Nusselt number. As the values of the porous medium parameters ([Formula: see text]–13) and nanolayer thickness ([Formula: see text]–12) elevate, the flow of skin friction coefficient and Nusselt number escalate in both porous disks. Moreover, as the values of the particle radius ([Formula: see text]% to 12%) and Reynolds number ([Formula: see text]1 to 2) amplify, the flow of heat transfer rate is diminished in the lower disk. Porous materials have a high surface area, which increases the contact area between the nanofluids and the porous medium. Porous materials can enhance the permeability of the fluid, reducing the pressure drop in the system. This can be advantageous in applications where high flow rates and low-pressure drops are desired, such as in heat exchangers or cooling systems. Porous materials can provide trapping sites or adhesion points for the nanoparticles, enhancing their dispersion and stability within the hybrid nanofluids. The porosity of the flow medium can therefore influence the purification and filtration efficiency of the hybrid nanofluids, ensuring their cleanliness and functionality. Understanding and controlling the porosity can help optimize the flow characteristics of the nanofluids in various applications. Overall this field will focus on refining hybrid nanofluid compositions, improving modeling techniques, validating experimental findings, optimizing system designs, and exploring new applications.

Double stratification impacts on MHD flow of micropolar nanofluid over a stretching surface inserted in porous medium

A comprehensive analysis of double stratification impacts on the MHD flow of a nanofluid produced by a porous stretching sheet has been conducted by using a numerical approach. The convective boundary conditions (CBCs) and Newtonian heating (NH) cases have been particularly considered in this problem. The governing equations defined in partial differential equations (PDEs) form are converted into ordinary differential equations (ODEs) by using appropriate transformations to meet the requirement of the used numerical method which is bvp4c in the MATLAB software package. Moreover, the numerical findings of the ODEs are displayed via graphical and tabular forms. The impacts of different flow parameters used in particular problem are discussed for the profiles of the velocity, temperature and concentration of micropolar nanofluid along the values of coefficient of skin-friction [Formula: see text], Nusselt number [Formula: see text], the Sherwood number [Formula: see text] and in relation with the variation parameters in general, CBC and NH are calculated to get a clear insight into the flow situation. Some of the findings show that the skin friction increases by increase in magnetic, permeability and material parameters in all cases but the skin friction is reduced by increasing the thermophoresis parameter for CBC and NH cases, respectively. Microrotation profile decreases when material parameter increases, whereas microrotation profiles enhance with increase of the magnetic field and porosity of porous medium.

Analysis of slip condition in MHD nanofluid flow over stretching sheet in presence of viscous dissipation: Keller box simulations

In present paper the mathematical study of Slip conditions on electrically conducting nanofluid over a vertically stretching sheet is presented with effects of viscous dissipation, thermal radiation and Soret and Dufor. The state of flow has been explained mathematically with help of partial differential equations. Model equations are made convenient for calculation by suitable transformations. With the help of Kellerbox technique unknown functions are approximated in flow situation. The influence of medium porosity, magnetic field strength, thermal radiation, Soret and Dufour effects, Buoyancy forces, viscous dissipation wall slip parameter, and heat source are investigated in detail. Effects of emerging parameters on velocity, concentration, temperature and profiles are investigated. In relation to the parameters involved, we can observe the rates of momentum, heat, and mass transfer near the stretching surface. It can be concluded from the study that velocity slip parameter accelerates with fluid motion, temperature of nanoparticles is enhanced due to rise in Eckert number and Biot number, while the Soret effect enhances the nanoparticles concentration close to the stretching sheet. Also, under common assumptions, the analytical approximations for the solution of the current model that were produced by using the Keller box approach were shown to be in very excellent agreement with certain early publications.

STUDY ON HYDROGEOLOGICAL ASPECT OF GALANG SUB-DISTRICT, BATAM CITY FOR ACCEPTANCE OF PROSPECTIVE NPP SITES

This is a preliminary study in the selection of Nuclear Power Plant (NPP) site in Batam's Barelang area to support industrial growth in the area in the future. The initial site selection was conducted in 2015 and 2017 in the pre-survey phase, considering hydrogeological aspect. The results of previous research showed four potential areas, namely Pasir Panjang Beach, Tanjung Batu, Dapur 3 and Tanjung Rame. The hydrogeological aspect plays an important role in the consideration of site acceptance. This is related to the consideration of potential flow of radioactively contaminated groundwater seepage in the site area, in the event of a potential release. The acceptance of the NPP site from the hydrogeological aspect is based on the site’s permeable geological formation and porosity condition, where groundwater can be stored. The purpose of the study was to assess the potential site in Galang Sub-district, Batam City based on hydrogeological aspects which include surface geology, groundwater productivity and rock porosity. Research methods include primary and secondary data collection, literature review and ranking analysis. The results showed that Tanjung Batu, Dapur 3 and Tanjung Rame could be chosen as the potential sites for NPP based on the hydrogeological aspect. The three regions have surface geology in the form of sandstone, clay and claystone rock with medium rock porosity level and medium groundwater productivity. Pasir Panjang Beach is less preferable because it has a high porosity of rocks with the productivity of aquifers is being spread widely.

Numerical Study on Heat and Mass Transfer of Evaporated Binary Zeotropic Mixtures in Porous Structure

This study investigates the heat and mass transfer characteristics of a binary mixture (R134a/R245fa) evaporated in a porous medium. The Eulerian model coupled with the multiphase VOF model and species transport equations is employed to establish a multi-component evaporation model. The effects of heat flux ranging from 200 kW/m2 to 500 kW/m2, porosity ranging from 0.4 to 0.6, and mass fraction ratios (R134a/R245fa) of 3:7, 5:5, and 7:3 are explored. The results indicate that an increase in heat flux contributes to an increase in the evaporation rate. For the overall evaporation rate, the evaporation rates of R134a and R245fa improve by 11.3%, 6.9%, and 16.3%, respectively, while the maximum improvement in heat transfer coefficient is only 1.4%. The maximum evaporation rate is achieved at intermediate porosity in the porous medium, and the highest heat transfer coefficient is obtained at a porosity of 0.4. With the increase in mass fraction, the evaporation rate of the corresponding species also increases, while the overall evaporation rate and heat transfer coefficient remain almost unchanged.

Effect of variable gravity on thermal convection in Jeffrey nanofluid: Darcy-Brinkman model

The subject under consideration in this research has various geophysical and astrophysical applications. Specifically, we investigate the impact of variable gravity on the onset of thermal instability within a layer of Jeffrey nanofluid confined in a Darcy-Brinkman porous medium. The solution of the fluid layer, which is positioned between two free-free boundaries, is determined using a linear stability analysis employing the normal mode technique. The Rayleigh number on the onset of convection is derived by using the Galerkin method. For stationary convection, the effects of different variable gravity parameters on the Jeffrey parameter, Darcy-Brinkman number, Lewis number, moderated diffusivity ratio, porosity of porous media and nanoparticle Rayleigh number are analyzed and presented graphically. The choices of using Jeffrey nanofluid as the base fluid for the study add novelty. Non-Newtonian fluids encompass a diverse range of substances, such as engine oils, oil extraction, wallpaper paste, and various biological liquids like blood etc.

Effect of Hump Configurations of Porous Square Cavity on Free Convection Heat Transfer

Free convection is widely used in engineering applications, including solar energy, electronic devices, nuclear energy, and heat exchangers. A computational simulation utilizing Ansys Fluent-CFD was employed to examine the natural convection heat transfer inside a square cavity filled with pure water and saturated metal foam as a porous medium (porosity ɛ =0.9). The enclosure's lower wavy wall exhibits a high temperature (Th), while the side and upper walls have a low temperature (Tc). For different Rayleigh numbers, the study examines hump configuration and the bottom wall hump number (N). The predominant design of heat transmission was improved using the circular hump design parameters of ɛ=0.9, N=4 and Tc= 25C˚ for different Ra. This resulted in significant improvements in heat transfer enhancement and energy enhancement which were enhanced by 1.13 times, for both. The authenticity research included determining the optimal design for the square enclosure. This involved estimating the effects of hump configure and number of humps for bottom wall of enclosure. These parameters have not been studied yet. The optimum case showed the highest heat transfer coefficient (h) at circular hump, N=4 and Ra = 30´103. While the standard case had N=0 and Ra = 5´103. The CFD simulation results indicate that the primary objective of the study was achieved through the optimal design, which resulted in a significant enhancement of hydrothermal performance for both heat transfer enhancement and energy enhancement 1.13 times compared to standard case.

Simulation investigation of the effect of heating temperature and porosity of porous media on the water evaporation process

Abstract In the study, carbon fiber as porous media, water and vapor as liquid and gas phases, respectively, the evaporation process of multiphase fluid in carbon fiber was explored based on the pore network model. The Micro-CT technology was adopted to reconstruct the structure of carbon fiber. The applicability of the CFD model adopted in the simulation was validated by experiment. Effect of heating temperature and porosity on the evaporation process was investigated in detail. The study found that at various heating temperatures, fluid temperature increases gradually with time, and tends to be constant when it reaches a certain temperature. When heating temperature is high, the constant temperature is large, the time to arrive the constant temperature is short. The liquid evaporation rate V e and the heat absorbed by the fluid every 0.5 s Q 0.5 increase first, and decrease gradually when a maximum value is reached. When the heating temperature is high, the maximum that V e and Q 0.5 can be obtained is also large, the required time is also short. In porous media with different porosity, fluid temperature rises to an approximate constant. V e and Q 0.5 increase first, and decrease gradually after reaching a certain maximum value. The maximum values of V e and Q 0.5 in porous media with porosity of 0.598 are the largest.

Effects of fluid flow in triple porosity medium on fracture width and its propagation during lost circulation control

During oil and gas well construction, lost circulation is one of the major challenges, and to migrate the problem, lost circulation materials (LCMs) should be used to resist fracture reopening with the hypothesis of fracture tip isolation. A triple porosity medium Darcy flowing model was developed to capture pressure distribution along the fracture surface, coupled with a dislocation-based fracture mechanics model to calculate fracture width and assess its propagation. In the model, mud cake on wellbore and fracture surface was considered as the third type of pore medium. The rest dual pore medium models incorporate the fluid flow through LCMs within fracture and the fluid leakoff from fracture into matrix. Sensitivity analyses of the key parameters, such as triple pore medium permeability, viscosity, wellbore pressure, plugging location are implemented, and the results emphasize that fracturing resistance effect of LCMs was controlled not only by the fluid flow within fractures, but also by the flow through mud cake on wellbore into matrix. To inhibit fracture reopening, in high permeability formations, both reduction of mud cake permeability and LCMs permeability were required while in lower permeability formations, fluid flow within fractures dominates and should be minimized. Besides, tip isolation can be achieved by extending sealing fracture length. However, high wellbore pressure or long operation time will result in elevated distances and unallowable fluid loss. Then, wellbore pressure must be reduced to cure the lost circulation. More viscosity of the fluid and total integrity plugging of LCMs are another major factor extending the time for pressure transfer to fracture tip, resisting unexpected fracture propagation. The method proposed in this study clarifies that tip isolation can be realized by fracture plugging and reduction of fluid filtration from the wellbore to matrix. Moreover, operational insights on LCMs design can also be derived.

Establishing a transient model of double-layer porous media interfacial evaporation system to reveal the relationship between porous media parameters and evaporation performance

The solar interfacial evaporation system can generate steam efficiently by localizing solar thermal energy to heat air–water interface, which has great application potential in seawater desalination, sewage treatment, and other fields. At present, the solar absorption layer and substrate of the interfacial evaporation system are mainly made of carbon materials, sponges, foams, and other porous media materials. However, the transient interfacial evaporation heat and mass transfer model and the influence of the physical parameters of porous media on evaporation performance have rarely been investigated. Based on heat and mass transfer theory of porous media, a transient theoretical model of porous media interfacial evaporation system is established in this paper. Moreover, the effects of physical parameters such as the porosity and thermal conductivity of porous media on evaporation performance are analyzed, and the accuracy of the model is verified through outdoor experiments. Based on the model calculation, as the porosity of the substrate increases, the evaporation rate and efficiency of the evaporator first increase and then decrease. The lower the thermal conductivity of the substrate and the higher the porosity and thermal conductivity of the solar absorption layer are, the better the evaporation performance is, but the growth rate gradually decreases. Therefore, the relationship between the improvement of evaporation performance and the increased cost caused by the improvement of material properties needs to be considered. The model and related theoretical analysis results can provide a theoretical basis for system structure optimization, material selection and other aspects of future research on double-layer porous media interfacial evaporation system.

Effects of viscous dissipation and activation energy for the MHD Eyring-powell fluid flow with Darcy-Forchheimer and variable fluid properties

The study of flow, heat and mass transference of Eyring-Powell fluid flowing in the existence of micro cantilever sensor is useful to create enhancement in the field of non-Newtonian fluid. The Darcy-Forchheimer porosity medium is considered for the flow analysis. The impressions of thermal radiation, viscous dissipation, heat generation, activation energy, Soret and Dufour effects are taken here to investigate the transference analysis of heat and mass rate. The flow phenomena are inspected under the influence of variable fluid properties. The motivation arises to work in this type of study is due to their applications in various industries such as steam generators, chemical reactors, system of dying package beds, several cooling and heating process and treatment of various diseases. The flow equations are occurred in the form of PDEs and then we utilize suitable similarity transformations to convert them into ODEs. These ODEs are numerically solved by the numerical technique namely, ‘shooting method’. The result reveals that the fluid velocity upsurges due to Hartman number and fluid coefficient whereas opposite impact is illustrious due to porosity coefficient, inertia coefficient and viscosity parameter. The finding indicated that we can upturns the temperature within the flow system due to thermal radiation, thermal conduction coefficient, Eckert number, heat generation and Dufour number. The rate of mass transference enhances for higher inputs of heat conduction coefficient, activation energy and Soret number while reverse behavior is seen due to temperature difference coefficient.

A mathematical model for pre-Darcy flow in low permeability porous media with stress sensitivity and the boundary-layer effect

Pre-Darcy flow is widely found in many naturally porous media, such as low-permeability clay. Nonflowing liquid boundary layers strongly influence fluid flow in low-permeability porous media. There is evidence in the literature that the porosity and permeability of low-permeability porous media change with variations in the effective stress, which is called stress sensitivity. In this paper, a new mathematical model for pre-Darcy flow in low-permeability porous media is developed and validated with experimental data from the literature. The nonflowing boundary layer thickness can be obtained from the proposed model instead of the empirical equation. First, an equation for calculating the pore radius with stress sensitivity is derived. Two parameters εp and εr are defined to calculate the nonflowing boundary layer thickness. Second, based on the capillary bundle model, the apparent flow velocity and apparent liquid permeability of the porous media are determined. Finally, the stress sensitivity and boundary-layer effects on the apparent flow velocity and apparent liquid permeability are analyzed. The research results show that the effect of stress sensitivity on apparent flow velocity increases with increasing pore compressibility. When the pore compressibility Cp>0 MPa−1, the apparent liquid permeability increases continuously with increasing pressure gradient. When the pore compressibility Cp=0 MPa−1, the apparent liquid permeability continues to increase with increasing pressure gradient. The apparent liquid permeability eventually tends to a constant value. The apparent flow velocity decreases with εp, and εr decreases at a given pressure gradient. The larger εp and εr are, the faster the rate of increase in the apparent liquid permeability with increasing pressure gradient. The proposed model can estimate the apparent flow velocity at different pressure gradients with a high accuracy. This work is important for understanding flow phenomena in low permeability porous media, such as landslides caused by seepage of water in clay.