Filtration Model Research Articles

Estimation of Infiltration Parameters: The Role of Pedotransfer Functions and Initial Moisture Conditions

Infiltration is one of the key components of rainfall-runoff processes at small catchment scale and is therefore a vital component of hydrological models. However, small catchments typically lack direct infiltration measurements and associated soil hydraulic characteristics, which are being substituted with soil databases and pedotransfer functions built upon these databases. We utilized two sets of pedotransfer functions Rosetta 3 and EUPTFv2 derived from two well known soil databases: UNSODA, based on US and global soil data, and EU-HYDI, based on EU soil data. Our study evaluates the impact of choosing one of these pedotransfer functions on modeled infiltration in relation to variable initial moisture content- It investigates the behavior of the simplified infiltration models of Green-Ampt and Philip under different initial conditions, in comparison with Richards’ complex solution, which is used as a benchmark. Our evaluation was conducted using soil property data available in the Czech Republic, covering ten USDA soil textures that are prevalent within the European Union. The applied initial moisture conditions were derived for specific suction pressure values that correspond to hydrological initial conditions ranging from near saturation to near the permanent wilting point. Variation in the derived parameters for the Green-Ampt and Philip models across the pedotransfer function sources results in considerable variation in the modeled infiltration. Our results imply different sensitivities of the simplified infiltration models to the initial moisture for different pedotransfer functions. The Green-Ampt and Philip infiltration models demonstrate better alignment with Richards’ complex solution with a low and consistent RMSE when using the Rosetta 3 pedotransfer function. By contrast, the EUPTFv2 model scenarios, show a clear trend of increasing RMSE towards dry initial conditions. As the European soil database is not an open data source, we only hypothesize that the distributions of soil characteristics in the European soil database differ systematically from those in UNSODA, and that this would cause different results of pedotransfer functions for different initial moisture conditions. The findings of this study are important for selecting appropriate infiltration models and, appropriate pedotransfer functions for estimating model parameters, and for assessing uncertainties connected with varying initial conditions.

Comparison of the mathematical modelling results of the relationship between cerebral ventricular size and capillary pressure based on experimental and clinical data

Aim of the study was to compare the results of mathematical modelling of the dependence between brain ventricle size and capillary pressure for humans and animals based on the equations of multicomponent poroelastic filtration for brain parenchyma. Material and methods. The study included two groups of animals - 4 male mice of each inbred line C57Bl/6 and BALB/C at the age of 12 weeks – and 4 healthy volunteers. The brain and cerebrospinal fluid system images of mice were obtained using an 11.7 T horizontal MR scanner, group of humans were examined using the Ingenia 3.0 T MRI scanner. An axial section at the level of –0.5 mm from bregma in the mouse groups and a frontal slice at the level of the middle of the bodies of the lateral and third ventricles, posterior to the foramen of Monroe in the human group were chosen as the geometry for mathematical modelling. Mathematical modelling is based on the stationary mathematical model of multicomponent poroelastic filtration. Multiple linear regression of mean ventricular wall displacement on fluid media interaction parameters was constructed to compare results obtained. Regression coefficients were compared using nonparametric analysis of variance based on the Kraskell–Wallis criterion and post-hoc Dunn’s criterion with Hill’s correction Results. A qualitative coincidence in the behavior of capillary pressure and mean ventricular wall displacement was demonstrated for the human and mouse groups. No significant differences were found between the two animal lines. For the animals characterized by small ventricular size (BALB/c), greater similarity to humans is observed than for the genetic line with hypertrophied ventricles (C57Bl/6). A significant difference between humans and mice is observed only for capillary-venous interaction. Conclusions. The low variance within groups and insignificant discrepancy between groups indicate the possibility of further accumulation of empirical data to establish correction coefficients of the animal model, which will bring it more in line with the model for humans. Thus, the analyzed models are sufficiently comparable with each other.

An Evaluation of the Relationship between Membrane Properties and the Fouling Mechanism Based on a Blocking Filtration Model

Microfiltration plays an increasingly important role in various fields. Consequently, elucidating the mechanism of membrane fouling has emerged as a pivotal issue that needs to be resolved. In this study, a blocking filtration model was employed to evaluate the effects of membrane properties on the fouling mechanism during the microfiltration of representative polysaccharides, namely sodium alginate, pectin, and xanthan gum. Microfiltration membranes composed of hydrophilic and hydrophobic PVDF, mixed cellulose ester, as well as hydrophilic and hydrophobic PTFE were used as filter media. The flux decline behavior was significantly affected by the membrane properties, with hydrophilic membranes exhibiting a slower decrease in filtration rate. The model analysis revealed a correlation between the blocking characteristic values and the membrane properties. Although the blocking index n showed membrane material dependence, the values of this parameter remained consistent across various filtration conditions, including the wettability of the membrane surface, solute concentration, and pressure (pectin: n = 1.86, 1.85, 1.50, and 1.50 for hydrophilic PVDF, hydrophobic PVDF, hydrophilic PTFE, and hydrophobic PTFE, respectively). The resistance coefficient k was influenced by the characteristics of the membrane surface; the k values of the hydrophobic membranes were higher than those of the hydrophilic ones (pectin: k = 0.00084, 0.00725, 0.00714, and 0.0384 s1−n/cm2−n for hydrophilic PVDF, hydrophobic PVDF, hydrophilic PTFE, and hydrophobic PTFE, respectively). The model calculations, based on the values of n and k, demonstrated a relatively good agreement with the experimental data.

An Optimal Model for Determination Shut-In Time Post-Hydraulic Fracturing of Shale Gas Wells: Model, Validation, and Application

The global shale gas resources are huge and have good development prospects, but shale is mainly composed of nanoscale pores, which have the characteristics of low porosity and low permeability. Horizontal drilling and volume fracturing techniques have become the effective means for developing the shale reservoirs. However, a large amount of mining data indicate that the fracturing fluid trapped in the reservoir will inevitably cause hydration interaction between water and rock. On the one hand, the intrusion of fracturing fluid into the formation causes cracks to expand, which is conducive to the formation of complex fracture networks; on the other hand, the intrusion of fracturing fluid into the formation causes the volume expansion of clay minerals, resulting in liquid-phase trap damage. At present, the determination of well closure time is mainly based on experience without theoretical guidance. Therefore, how to effectively play the positive role of shale hydration while minimizing its negative effects is the key to optimizing the well closure time after fracturing. This paper first analyzes the shale pore characteristics of organic pores, clay pores, and brittle mineral pores, and the multi-pore self-absorption model of shale is established. Then, combined with the distribution characteristics of shale hydraulic fracturing fluid in the reservoir, the calculation model of backflow rate and shut-in time is established. Finally, the model is validated and applied with an experiment and example well. The research results show that the self-imbibition rate increases with the increase in self-imbibition time, and the flowback rate decreases with the increase in self-imbibition time. The self-imbibition of slick water is the maximum, the self-imbibition of breaking fluid is the minimum, and the self-imbibition of mixed fluid is the middle, and the backflow rates of these three liquids are in reverse order. It is recommended the shut-in time of Longmaxi Formation shale is 17 days according to the hydration and infiltration model.

Modeling of filtration and scattered fracture in fractured zones with abnormally high formation pressure

A model of a double-porosity medium with a damageable matrix under the action of anomalously high formation pressure is discussed. It is assumed that damage in the matrix enhances mass transfer between the subsystems of the medium. The problem of the coupled process of fracture and fluid flow in a fractured porous medium in one-dimensional (1D) statement is solved numerically. The limited spatial extent of the zone of conductive fractures is taken into account.

Rainfall Infiltration through Stratified Colluvial Deposits: Analytical Approach vs. Numerical Modelling

This work investigates the rainfall infiltration process within homogeneous and stratified colluvial deposits caused by short (1–3 h) and intense (40–90 mm/h) rainfall, using both analytical and numerical infiltration modelling. The findings of the investigation demonstrate that the classic Green–Ampt model can be employed effectively to study homogeneous colluvial covers with permeability equal to or lower than kw = 10−5 m/s and that are subject to a 1 h rainfall with intensity I ≥ 45–50 mm/h. In these circumstances, a top-down saturation front forms within the colluvial deposit, leading to the saturation of a 70–100 cm-thick layer. This critical condition occurs every 5–10 years in the mountain area of the Friuli Venezia Giulia Region (NE Italy), which corresponds to a lower return period of critical hydrologic events when compared with other mountain basins in the Alps due to the higher initial degree of saturation characterising colluvial covers in this area (70–95%). When analysing stratified colluvial covers, the Dagan–Bresler approximate model, as well as the numerical modelling, emphasised the strong influence that abrupt variations in the permeability of the various soil layers have on the infiltration process at depth. In particular, the presence of a top organic soil horizon that is rich in macro-pores and is characterised by a higher permeability (k = 10−4 m/s) actually reduces the possibility of surficial ponding, which is the basic condition of the “piston” models. The highly permeable top soil allows for a rapid downward infiltration up to contact with the underlying colluvial material, which is less permeable (k = 10−5 m/s). Therefore, a perched water table forms starting from the organic soil–colluvium interface, originating pore–water overpressures within the colluvial deposit, with maximum values in the order of 5–10 kPa.

An improved analytical model for rainfall infiltration in embankment slope constructed with red-bed solid waste

Rainfall infiltration is one of the key triggering factors of landslide in embankment slope constructed with red-bed solid waste, so rainfall infiltration analytical model plays a critical role in the landslide prediction. However, the complexity of the red-bed solid waste slope’s components has not be considered by the conventional Green-Ampt infiltration model, which has also ignored the influence of the dynamic variation of unsaturated section in the wetted zone of slope under rainfall infiltration. This study proposed an improved rainfall infiltration analytical model for red-bed solid waste slope considering the properties of rock fragments and the dynamic development of saturated–unsaturated proportion in the wetted zone, and developed a solution program for the improved model. The proposed model was validated by comparing the rainfall infiltration experiment, finite element numerical simulation, and traditional Green-Ampt method. The proposed model was then applied to parameter sensitivity analyses to reveal the rainfall infiltration characteristics of red-bed solid waste slope. The results indicate that the rainfall infiltration model proposed in this work perform better to predict the rainfall infiltration process of red-bed solid waste slope, and it will lay the groundwork for more reliably utilizing the red-bed solid waste in transportation geotechnics.

Evaluation of infiltration models in clay loam and laterite soils under field conditions

The purpose of the investigation is to calculate soil infiltration rates with the help of infiltration models. The infiltration model helps to design and evaluate surface irrigation systems. The study calculated constant infiltration for two types of soils (clay loam soil and laterite soil) under field conditions (Unploughed and Ploughed). The double-ring infiltrometer has been implemented to experiment. The value of various constants of the models was calculated using the approach of averages counselled through a graphical technique. Fitting infiltration test data to prominent infiltration models such as Philip’s, Horton's and Kostiakov’s and The Nash- Sutcliffe efficiency (NSE), coefficient of determination (R2) and root mean square error (RMSE) statistics are used to evaluate the effectiveness of the model. The results indicate that Philip's model is the most reliable, with R2, NSE, and RMSE values ranging from 0.9044-0.9677, 0.294-0.957 and 1.2647-5.7129, respectively. Therefore, under identical circumstances and without any kind of infiltration information, the above model can be employed to artificially produce infiltration information.

Infiltration Models in EnergyPlus: Empirical Assessment for a Case Study in a Seven-Story Building

The current decarbonization transition to be achieved by 2050 according to the European Council has given great prominence to the use of Digital Twins as tools for energy management. For their correct operation, it is essential to control the uncertainties of the energy models, which lead to differences between the measured and predicted data. One of the key parameters that is most difficult to assess numerically is air leakage. The existent infiltration models available in EnergyPlus were developed to be applied in low-rise residential buildings with fewer than three stories. Therefore, it is common to rely on air leakage equations employing predefined coefficients. This research presents an empirical assessment of the performance of two EnergyPlus air leakage models, the “Effective Leakage Area” and the “Flow Coefficient”, in predicting dynamic infiltration within the attic of a seven-story building. Blower door tests, along with the application of CO2 tracer gas, were conducted to establish coefficients for the models. Then, they were evaluated in three independent periods according to the criteria established in the American Society for Testing Material D5157 Standard. Those models that only used in situ coefficients consistently met the standard across all three periods, demonstrating for both equations their accurate performance and reliability. For the best model derived from tracer gas data, the R2 and NMSE values are 0.94 and 0.019, respectively. In contrast, the model developed using blower door test data and EnergyPlus default values presented a 64% reduction in accuracy compared to the best one. This discrepancy could potentially lead to misleading energy estimates. Although other software options exist for estimating infiltration, this study specifically targets EnergyPlus users. Therefore, these findings offer valuable insights to make more informed decisions when implementing the infiltration models into energy simulations for high-rise buildings using EnergyPlus.

The liver and muscle secreted HFE2-protein maintains central nervous system blood vessel integrity

Liver failure causes breakdown of the Blood CNS Barrier (BCB) leading to damages of the Central-Nervous-System (CNS), however the mechanisms whereby the liver influences BCB-integrity remain elusive. One possibility is that the liver secretes an as-yet to be identified molecule(s) that circulate in the serum to directly promote BCB-integrity. To study BCB-integrity, we developed light-sheet imaging for three-dimensional analysis. We show that liver- or muscle-specific knockout of Hfe2/Rgmc induces BCB-breakdown, leading to accumulation of toxic-blood-derived fibrinogen in the brain, lower cortical neuron numbers, and behavioral deficits in mice. Soluble HFE2 competes with its homologue RGMa for binding to Neogenin, thereby blocking RGMa-induced downregulation of PDGF-B and Claudin-5 in endothelial cells, triggering BCB-disruption. HFE2 administration in female mice with experimental autoimmune encephalomyelitis, a model for multiple sclerosis, prevented paralysis and immune cell infiltration by inhibiting RGMa-mediated BCB alteration. This study has implications for the pathogenesis and potential treatment of diseases associated with BCB-dysfunction.

Modeling of Filtration and Scattered Fracture in Fractured Zones with Abnormally High Formation Pressure

Modeling of Filtration and Scattered Fracture in Fractured Zones with Abnormally High Formation Pressure

Fabrication and characterization of low-cost kaolin based tubular ceramic membrane for microalgal harvesting

Tubular membranes with natural kaolin as a key precursor were fabricated in this study. Membranes with five different compositions were prepared, and the optimized membrane (77% kaolin, 2% boric acid, 2% sodium metasilicate, 4% sodium carbonate, and 15% calcium carbonate) was used to harvest microalgal culture. The pore size, porosity, water permeability, chemical and mechanical strength of the membranes were investigated. An assessment of the cost related to the membrane fabrication was also conducted. The fabricated membranes had porosity of ∼26–47%, a pore diameter of 0.123–0.182 µm, water permeability of 4.2 × 10−8 – 17.1 × 10−8 m3 m−2 s−1 kPa−1, along with good mechanical and chemical strength. The optimized membrane was tested for microfiltration of microalgae Monoraphidium sp. KMC4 with 1.5 g L−1 of initial concentration at a persistent cross-flow rate (1.11 × 10−5 m3 s−1) and various transmembrane pressures (69 kPa - 345 kPa). The separation results yielded an average permeate flux of 1.85 × 10−5 m3 m−2 s−1 at 276 kPa transmembrane pressure. Complete algal cell recovery and substantial nutrient passage (>88%) were observed within the pressure range of 69 kPa to 345 kPa. Fouling mechanism was explained by fitting four distinct pore-blocking models, of which the cake filtration model provided the most accurate fit as compared to the complete, intermediate and standard pore-blocking models. Additionally, the total organic carbon varied in the range of 31.6–63.2 mg L−1. This essentially explained the source of pore blocking.

Deep bed filtration model for cake filtration and erosion

Deep bed filtration model for cake filtration and erosion

Effects of varying tillage practices and weed control methods on the efficacy of infiltration models.

Agricultural land preparation and weed control techniques are essential farm management tools that affect the dynamics of soil water infiltration and the estimation accuracy of infiltration models. To analyse the interaction effect of tillage and weed control methods on the changes in soil physical properties and the efficacy of infiltration models, an experiment was conducted on a sandy clay loam forest ochrosol at Hodzo near Ho in Ghana. Four tillage systems (No Tillage [NT], Reduced Tillage [RT], Plough + Harrow + Ridging [PHR], and Deep Tillage + Plough + Harrow + Ridging [DPHR]) and three weed control methods (Hoeing [H], Machete [MAT] and No Weeding [NW]) were employed. The study also tested the reliability of the models (Kostiakov, Philip, and Horton) using the goodness of fit statistical criteria: Root mean squared error (RMSE), Mean absolute error (MAE), Coefficient of determination (R2), and Nash-Sutcliffe efficiency (NSE). The results show that conservation tillage systems (CsT) and conventional tillage systems (CT) with MAT weeding treatments recorded the highest moisture content across the studied soil profile, especially for NT x MAT (11.189%) which was significant (p < 0.05) in the 15-30 cm layer; the lowest were observed in the CsT and CT with H weeding interactions, especially for the DPHR x H (8.086%). Comparing the interaction effect on the soil infiltration, the highest mean infiltration rate was significant (p < 0.05) under the NT X H treatment combination whilst the lowest infiltration rate was recorded in the DPHR X H and PHR X NW treatment combinations. The efficiency of the fitting models (Kostiakov > Horton > Philip) highly prioritised the soil tillage operations and weed management under the treatments DPHR x MAT > DPHR x NW > DPHR x H > RT x MAT > PHR x NW > PHR x MAT > NT x NW > RT x MAT > PHR x H > RT x H > NT x MAT > RT x NW > NT x H in that order. The trend shows that the increase in tillage intensity and the decrease in weed management intensity induce the quality of the estimation process and vice versa. The study, therefore, identified the use of machete (MAT) with DPHR under the Kostiakov model as the efficient land management for modelling the cumulative infiltration characteristics of the sandy clay loam ochrosols of the study area.

The effect of soil texture, layering and water head on the infiltration rate and infiltration model accuracy

AbstractInfiltration is one of the most important physical characteristics of soil and depends on various factors. This study investigated the influence of soil texture, layering and water head on the soil water infiltration rate. It also selected the most accurate infiltration models to determine the water infiltration rate in homogeneous and heterogeneous soil profiles. Experiments were carried out in four soil containers with a length, width and height of 20 × 20 × 70 cm. Treatments consisted of two soil textures (sandy loam, SL; clay loam, CL), four soil profiles (homogeneous texture, SL and CL; and heterogeneous texture, lighter texture on the top, SL/CL, and heavier texture on the top CL/SL) and three water head sizes (4, 7 and 10 cm). Several models were used to determine the water infiltration rate under homogeneous (Kostiakov, modified Kostiakov, Philip, Horton, traditional Green–Ampt, modified Green–Ampt and HYDRUS‐1D) and heterogeneous soils (traditional Green–Ampt, modified Green–Ampt and HYDRUS‐1D). According to the results, the infiltration rate decreased over time and along the soil profile. Nevertheless, it jumped at the interface of two‐layered soils when the heavier soil was in the bottom layer (SC treatments) due to the high potential of the second layer, and then it decreased. In the reverse layering, the infiltration rate in the interface was lowest (CS treatments) because of the higher hydraulic conductivity of the second layer. Additionally, the infiltration rate increased with increasing water head, but the rate of this increase was higher by changing the water head from 7 to 10 cm. The results of infiltration models showed that the accuracy of these models was higher in clay loam texture than in sandy loam texture. The modified Green–Ampt was the most accurate model in homogeneous and layered soils, with average RMSE of 0.0204 and 0.019, respectively. The Horton model had the weakest simulation in homogeneous soils, with an average RMSE of 0.1299. Additionally, the accuracy of HYDRUS‐1D in layered soils was less than that in homogeneous soils (NS of 0.95 and 0.85, respectively), and its accuracy decreased with increasing water head in most treatments.

Microscopy of Woven and Nonwoven Face Covering Materials: Implications for Particle Filtration.

A suite of natural, synthetic, and mixed synthetic-natural woven fabrics, along with nonwoven filtration layers from a surgical mask and an N95 respirator, was examined using visible light microscopy, scanning electron microscopy, and micro-X-ray computed tomography (µXCT) to determine the fiber diameter distribution, fabric thickness, and the volume of solid space of the fabrics. Nonwoven materials exhibit a positively skewed distribution of fiber diameters with a mean value of ≈3 μm, whereas woven fabrics exhibit a normal distribution of diameters with mean values roughly five times larger (>15 μm). The mean thickness of the N95 filtration material is 1093 μm and is greater than that of the woven fabrics that span from 420 to 650 μm. A new procedure for measuring the thickness of flannel fabrics is proposed that accounts for raised fibers. µXCT allowed for a quantitative nondestructive approach to measure fabric porosity as well as the surface area/volume. Cotton flannel showed the largest mean isotropy of any fabric, though fiber order within the weave is poorly represented in the surface electron images. Surface fabric isotropy and surface area/volume ratios are proposed as useful microstructural quantities to consider for future particle filtration modeling efforts of woven materials.

Pressurised sand bed filtration model: Set up and energy requirements for a filtration cycle

Rapid deep bed filtration is a common process of drip irrigation systems to prevent emitter clogging. The particle retention in the porous media increases the pressure difference between the filter's inlet and outlet. Commercial operational instructions preset a threshold value of this pressure difference to define the end of a filtration cycle. Accurate particle retention models may contribute in the determination of this setpoint to improve the energy efficiency of a filtration cycle. A two-step method was developed to calibrate a phenomenological particle retention model with transient data of the filter pressure drop, and retained mass in different media slices at the end of the filtration cycle. The first step used the vertical profile of accumulated mass in the media to fit the input parameters of the specific deposit rate equation. The second step applied the pressure data to find the input parameters of the pressure drop equation. The use of constant values of the input parameters provided reasonable results for all tests. When the end of the filtration cycle was set by applying a threshold pressure drop, both particle mass accumulation and removal efficiency decreased as the flow rate increased. By using the threshold pressure drop criterion, the energy consumption in a filtration cycle also decreased as the flow rate increased, but the energy consumption per unit of filtered water volume increased. For high setpoint pressure drop values, a minimum of the energy consumption per unit of filtered water volume as a function of the flow rate was found.

Using an ecosystem service model to inform restoration planning: A spatially explicit oyster filtration model for Pensacola Bay, Florida.

The development of science-based restoration goals that reflect the primary motivation of stakeholders is a key factor leading to large-scale, long-term restoration successes. The ability to predict the potential ecosystem service delivery from restoration can inform the setting of appropriate goals and facilitate the strategic planning of restoration activities. While recovery of the ecosystem services provided by oyster reefs is a regularly cited reason for undertaking restoration, few examples exist where large-scale oyster habitat restoration plans have been informed using ecosystem service functions. Such an approach is currently being implemented in the Pensacola Bay System, Florida, where a broad coalition of partners and community stakeholders are utilizing a watershed approach to restoring oysters with the aim of restoring oysters for multiple objectives including habitat, ecosystem services, and wild harvest and aquaculture. Through the process of developing a habitat management plan, water filtration was identified as a key ecosystem service by the stakeholders. To support restoration planning we derived a spatially explicit estimate of water filtration services provided by the eastern oyster in the Pensacola Bay system by linking an oyster habitat suitability map to a hydrodynamic-oyster filtration model. This spatially explicit model allowed us to identify the areas where restored oyster reefs have the potential to provide the greatest increase in filtration service as well as provide spatially explicit estimates of the potential filtration provided by oyster habitat restored. Such information is useful in restoration planning and management and for stakeholder engagement, outreach, and education programs.

Three-dimensional stability analysis of two-layer slope based on GIS under rainfall conditions

ABSTRACT Geographical Information System (GIS), as a computer-based tool, allows easy data processing, and complex geospatial data management and analysis, making it a suitable environment for 3D slope stability analysis. The study proposed an integral scheme for 3D stability analysis of a two-layer slope under rainfall conditions based on GIS. To construct the 3D model of the slope and the sliding surface in GIS, a method was proposed in which a grid is divided into four sections to improve the accuracy of calculation, and the extension of the Hovland model and the Janbu model in unsaturated soil was discussed. The transformation equations between the global coordinates and the ellipsoid local coordinates were derived to locate the critical sliding surface, and how to apply constraints of parameters to improve search efficiency was put forward. An infiltration model for a two-layer slope was built so that the changes in water content in slope at the both stages of rainwater infiltration and redistribution can be determined, and the model was integrated into the GIS-based scheme to assess the impact of rainfall on slope stability. Finally, two case studies were implemented to demonstrate the validity and correctness of the proposed scheme.

A multivariate model of deep slurry infiltration in slurry shield tunneling

There are multiple factors that influence slurry infiltration processes, making it necessary to consider these factors when comprehensively analyze the process of deep slurry infiltration. Firstly, this paper proposes a theoretical model for multivariate deep slurry infiltration and introduces a control equation governing the variation of the slurry infiltration domain. This equation accurately reflects changes in the infiltration boundary during the slurry infiltration process and meticulously responds to variations in pore water pressure, slurry retention rate, and pressure conversion efficiency-critical factors in the slurry infiltration process. Then, the Lax method was used to solve the corresponding control equations of the model, realizing the effective solution of the implicit partial differential equation system. The computational results of the theoretical model were compared with the indoor test results. Finally, the validated model was used to analyze the deep slurry infiltration characteristics. The results indicate that the theoretical model can reasonably simulate the deep slurry infiltration process, and the use of dynamic computational domain to describe deep slurry infiltration process is more in line with engineering practice. Furthermore, this model has guiding significance for the division of deep slurry infiltration stage and slurry surface infiltration stage.