Fluid-solid Coupling Model Research Articles

Effects of Injection Pressure and Duration on Alternate High-Pressure Water-Gas Sequestration of Coalbed Methane

To solve the problems of high ground stress, low permeability and low coalbed methane recovery, this paper proposes an optimization technique for improving coal seam permeability, i.e., water-gas alternating displacement enhances coalbed methane recovery (WGA-ECBM). The seepage field of gas-water two-phase flow is derived and a multiphase fluid-solid coupling model between physical parameters is established. This paper studies the effects of three displacement methods (water injection, air injection and water-air alternating injection) on methane recovery. The influence of injection pressure and injection duration on the WGA-ECBM technique is also studied. The V C H 4 (methane volume fraction) of the three methods decrease by 17.84%, 12.72% and 44.62% in the borehole, respectively. The influence ranges (diameter) are 24.68 m, 149.73 m and 135.21 m, respectively. Pore pressures and displacement area also increase as injection pressure increases in a 5 MPa gradient. And V C H 4 in the borehole decrease by 35.57%, 43.82%, 46.78% and 51.72% from 10 MPa to 25 MPa, respectively. With the increase of injection duration, the influence range expands significantly but the pore pressure distribution changes little. And V C H 4 decreases slightly in the displaced area and the displacement velocity gradually slows down. This technique promotes multiphase flow migration in coal seam, and increases the displacement range and methane recovery efficiency.

Numerical Simulation of Boulder Fluid–Solid Coupling in Debris Flow: A Case Study in Zhouqu County, Gansu Province, China

Boulders mixed with debris flows roll downstream under interactions with debris flow slurry, which poses a great threat to the people, houses, bridges, and other infrastructure encountered during their movement. The catastrophic debris flow in Zhouqu County, which occurred on 7 August 2010, was used as an example to study the motion and accumulation characteristics of boulders in debris flows. In this study, a fluid–solid coupling model utilizing the general moving objects collision model and the renormalization group turbulent model was used in the FLOW-3D software, treating boulders with different shapes in the Zhouqu debris flow as rigid bodies and the debris flow as a viscous flow. Numerical simulation results show that this method can be used to determine the motion parameters of boulders submerged in debris flows at different times, such as the centroid velocity, angular velocity, kinetic energy, and motion coordinates. The research method employed herein can provide a reference for studying debris flow movement mechanisms, impact force calculations, and aid in designing engineering control structures.

Study on craniocerebral dynamics response and helmet protective performance under the blast waves

Blast traumatic brain injury (BTBI) is the major type of injury, existing helmets are mainly designed to resist bullets and fragments, so it is necessary to study the explosion- proof performance of helmets. In this study, the validated finite element models including head, Advanced Combat Helmet (ACH), fully-covered helmet based on the prototype of Iron Man helmet were used to establish the blast wave-helmet-head fluid–solid coupling model, dynamic responses such as blast wave flow field pressure, skull stress, brain tissue pressure, displacement and energy absorption were obtained through simulations.Quantitatively analyzed the protective performance of helmets show that wearing ACH can’t reduce the severity of brain injury, but rather increased about 5% compared to without wearing helmet. Wearing fully-covered helmet with same thickness as ACH can reduces the degree of brain injury by about 5%. Fully-covered helmet with same weight as ACH provides better protection, which can reduce the degree of brain injury to 65% than without wearing helmet. Research show that optimizing the shape and thickness of helmet can greatly reduce the brain injury, lighter fully-covered helmet can provide better protection for the brain under blast wave environment. The results can provide reference for new design of helmets.

Comparative Analysis of the Drainage Effect Based on Improved Hydraulic Flushing Technology.

The composite structure coal seam composed of coal seams with strong and weak deformation is widely developed in the coalfields in Southwest and North China. The implementation of hydraulic flushing in soft coal stratification can effectively relieve the pressure and increase the permeability of coal seams, so it can improve the effect of underground gas drainage. Taking Xin'an Coal Mine, West Henan province, China, as an example, this paper discussed the influence of hydraulic flushing on the gas drainage effect in soft coal seam by means of field investigation, numerical simulation analysis, similar simulation, and field verification. As the same time, a fluid-solid coupling model of gas drainage based on layered hydraulic flushing technology is established, and the differences and reasons of the gas drainage effect after layered hydraulic flushing in soft and hard composite coal seams are discussed. The results show that (1) the gas concentration change of layered flushing boreholes intermittently rises rapidly, and there is no obvious law, while that of conventional flushing boreholes shows a regular attenuation trend. (2) According to the numerical simulation results, after 30 days of drainage, the effective drainage range of layered hydraulic flushing is 8.33 m3 higher than that of conventional hydraulic flushing. (3) There are not a few points where the single borehole concentration data of each row of layered flushing boreholes are in the high level. Meanwhile, the maximum concentration of all conventional flushing boreholes is less than 40%, except that one borehole can reach 70% at the beginning. (4) All of the average gas drainage concentration of single-row boreholes using layered flushing technology is greater than 20% in the vast majority of days, while that of the conventional flushing boreholes is only higher than 20% in row 5 and has only lasted for 2 days. The difference of the drainage effect is obvious.

Numerical Simulation of Hydraulic Fracturing and Penetration Law in Continental Shale Reservoirs

The vertical heterogeneity of continental shale reservoirs is strong, the difference between lithology and stress between layers is large, the weak interface between layers develops, and the hydraulic fracture penetration and expansion are difficult, resulting in poor fracturing transformation effect. In view of this, based on the finite element and cohesive element method, this paper established a fluid-solid coupling model for the hydraulic fracture propagation through the continental shale and studied the control mechanism and influence law of various geological and engineering parameters on the hydraulic fracture propagation through the continental shale reservoir using single factor and orthogonal test analysis methods. Interfacial cementation strength between high layers, high vertical stress difference, low interlaminar stress difference, low tensile strength difference, low elastic modulus difference, high pressure fracturing fluid viscosity, and high injection displacement are conducive to the penetration and expansion of hydraulic fractures. The primary and secondary order of influence degree of each factor is: interlaminar interface cementation strength > interlaminar stress difference/tensile strength difference > fracturing fluid viscosity/injection displacement > vertical stress difference > elastic modulus. In addition, engineering application research has also been carried out, and it is recommended that the injection displacement during early construction should not be less than 3 m3/min, and the fracturing viscosity should not be less than 45 mPa·s. The field application effect is good, which verifies the engineering application value of the model.

Effect of elastic throttling structure on the load performance of aerostatic radial bearing

PurposeThis research aims to combine the throttling structure with the elastic element to enhance the load performance of aerostatic radial bearing.Design/methodology/approachIn this research, a fluid–solid coupling model of the elastic throttling structure is established while considering the interaction between the elastic element and the flow field. The effects of elastic element structural parameters on the stiffness and load capacity of aerostatic radial bearing are then researched. Finally, the effect of elastic element modulus on air film load performance and elastic element deformation is analyzed.FindingsThe results indicate that the aerostatic radial bearing with elastic element can significantly improve the load capacity and stiffness when compared to the common aerostatic bearing. By choosing the proper combination of parameters, the load performance can be improved by at least 16%.Originality/valueThe throttling structure of aerostatic bearing is optimized in this work, which significantly enhances the load performance of the aerostatic bearing.

The jetting process and spreading characteristics of the power-law fluids for material jetting process

Materials jetting, known as one of the 3D printing technologies, is widely applied in microelectronics packaging, biology and ceramic 3D printing due to its ability to print multi-materials by drop-on-demand. However, most of the materials are power-law fluids in 3D printing applications, the generation of satellites during the jetting process and droplet spreading characteristics are unclear and they have a great effect on the quality of the printout. In this paper, a common electromechanical and fluid-solid coupling model of the jet dispenser and observation platform of the jetting process are established. This modeling method is also suitable for other needle-driven jet dispensers. A commercial UV resin is adopted to study the jetting process of power-law fluid. To reveal the mechanism of satellite generation, the effects of input signals (rising time and falling time) on the dynamic characteristics of the needle and the jetting process are analyzed. On the basis thereof, the effectiveness of the optimal control parameters is demonstrated to eliminate satellites. In addition, the simulation and experimental results show that the falling time and fluid pressure can be controlled to adjust the spreading diameter and height of the droplet. Subsequently, the minimum line width of 0.276 mm is successfully printed with a nozzle of 0.07 mm.

A 3D peridynamic fluid–solid coupling damage model of hydraulic fracture propagation

AbstractPeridynamics provides a new model and method for describing material failure and fracture development. In this study, a three‐dimensional peridynamic fluid–solid coupling model of hydraulic fracture propagation was established. The model considered the effective stress and transverse isotropic characteristics of shale and was used to simulate the propagation of the hydraulic fracturing of bedding shale. The damage evolution influence on the fracture propagation behavior, as well as the interactions between the fractures and bedding, were investigated through a simulation. The damage evolution had discontinuous characteristics. In addition to the incomplete damage zone adjacent to the complete damage zone, there were nonadjacent and incomplete damage zones. The fracture and propagation of hydraulic fractures followed the same trends of damage evolution. Moreover, the simulation and related experimental results confirmed the existence of a large amount of mesodamage in the region far from the fracture, and the concept of “remote damage” in hydraulic fracturing was proposed. Finally, the damage evolution of the propagating fractures and bedding was simulated. The simulation showed that bedding hinders the propagation of main fractures in the height direction but promotes the complexity of fracture propagation. Overall, this study provides a new theoretical understanding and research ideas for an in‐depth understanding of hydraulic fracture geometry and the accurate calculation of Stimulated Reservoir Volume.

Numerical simulation of rock erosion performance of a high-speed water jet using an immersed-body method

Water jet drilling (WJD) is an effective technique for drilling micro-holes in the subsurface for reservoir stimulation . This study aims to investigate numerically the failure mechanism of rock during WJD and to assess the WJD performance before drilling. A 3D fluid-solid coupling model is developed for simulating WJD by coupling a mechanical solver based on the combined finite-discrete element method (FDEM) with a fluid solver using an immersed-body method. The new numerical model is capable of simulating crack initiation and propagation and fragment removal under the impact load of a high-speed water jet. The poroelastic effect is implemented via Biot's theory of poroelasticity . The numerical results show that: (1) rock failure is only observed in the Gildehaus sandstone with the lowest strengths among the three types of rock tested, (2) most of the cracks are tensile failures and pure shear cracks are rare, mixed mode cracks account for 15%∼40% of the total crack number depending on the mechanical boundary conditions , (3) increased water back pressure significantly suppresses jet erosion. The poroelastic effect on the rock failure is insignificant in the mesoscale simulations and will be further investigated using a microscale model in future.

Numerical research of lubricant flow behaviors in rolling-slider bearing under the microscopic scale

In order to reveal the fluid flow behaviors of rolling-slider bearing under the microscopic scale, a full-basin fluid-solid coupling model with different slider shapes is established by CFD and multiple mesh technology. The flow trajectory of fluid is monitored by a slip mesh model, and the dynamics performance is also compared. A reasonable model of rolling-slider bearing is derived. Moreover, this paper analyzes the influence of bearing speed, fluid viscosity, and fluid velocity on fluid behavior. The results show that there is a recirculation flow of fluid and vortex shedding in rolling-slider bearing during the process of lubrication. The distribution of maximum pressure on rolling body is in the orientation of bearing rotation. The fluid flow rate near inner ring is greater than that at outer ring. The established model is suitable for simulating the lubricant flow behaviors of rolling-slider bearings under different working conditions. It also provides a new method for designing the parameters of rolling-slider bearings and the function of working conditions.

Stability Analysis of a High-Steep Dump Slope under Different Rainfall Conditions

The existing slope stability research, which is based on the fluid–solid coupling theory, is mainly focused on the slopes of central and eastern China. The impact of rainfall on the stability of the dump slope has often been ignored. It is worthwhile to reveal the mechanism of the fluid–solid coupling mechanics of dump slopes in the arid desertification area of northwest China under the maximum precipitation. The method of combining the seepage mechanics theory with the geomechanics theory was adopted. Darcy’s law and the mass conservation law were introduced to derive and establish the fluid–solid coupling analysis method. Taking the Xinxing Coal Mine in Wuhai City, China, as an example, the finite element software ABAQUS was used to construct the fluid–solid coupling model for slope stability analysis with unsaturated soil. The equivalent rainfall intensity of 68 mm/h for 1 h and 18 mm/h for 24 h was designed in the simulation, respectively. Four different types of initial water content (i.e., 1.72%, 7.34%, 14.69%, and 22.03%) of the dump slopes were defined as the initial conditions. The high-steep slope was compared to the standard slope. Therefore, a set of sixteen rainfall schemes was proposed. The variation regularity of slope stability was thoroughly discussed in regards to four areas: vertical deformation, pore water pressure distribution, equivalent plastic strain, and safety factor. As was expected, the research showed that the slope height and angle have a significant effect on the slope stability. When high-intensity rainfall occurs for a short duration, the slope tends to be more stable as the initial water content increases on the slope. When low-intensity rainfall occurs over a long period, the slope stability reduces if the initial water content is too high or too low in the slope.

Numerical simulation on fluid flow and temperature prediction of motorcycles based on CFD

The temperature distributions of solid parts in motorcycle is of significance to give an quantitative cooling performance evaluation for assuring the feasible running of motorcycles. Therefore, the present work mainly tend to systematically investigate the flow velocity distributions and temperature distributions of motorcycle by constructing a fluid-solid coupling analysis model via CFD using commercial software STAR-CCM+. Based on computational fluid dynamics method (CFD), fluid-solid coupling analysis model of motorcycle is realized by integration of porous media model, heat exchanger model and fluid-solid co-simulation module and numerical investigations on fluid flow and temperature distribution characteristics at maximum running speed 120 km/h is conducted. The thermal management experiment test considering three pressure measuring points, fifty-eight temperature measuring points and two flowrate measuring points is investigated to know its cooling performance of motorcycle experimentally and provide experimental data for validation of simulation model. Results show that CFD results are in good agreement with experimental results and the maximum relative error is less than 10%. Thermal management test results show that left and right cylinder temperature, oil temperature are 183 °C, 190 °C, 120 °C, respectively. Air temperature of radiator at inlet side are from 29 °C to 35 °C, Air temperature of radiator at outlet side are from 60 °C to 90 °C, water temperature of radiator inlet and radiator outlet are 96 °C and 89 °C, respectively. The air temperature around left thigh and right thigh are 31 °C and 34 °C, respectively. CFD results show that there is no flow dead-zone and high temperature area in engine main components and muffler front cover. The temperature around motorcycle covering and both two thigh is not high due to the radiation of high-temperature components and the outlet temperature of exhaust muffler is also not high, which will not cause enough heat feelings to the riders for being uncomfortable. From the temperature field analysis results, there is no obvious poor temperature distributions for thermal comfort of this motorcycle. The results of this work can provide simulation data support and theoretical reference for the design and development of motorcycle and its crucial components. Through engineering applications of fluid-solid coupling model in a two-cylinder engine and electric generator, it verifies the universality of fluid-solid coupling model and it can be an effective way for the temperature predictions of industrial products.

The Method of Determining Layer in Bottom Drainage Roadway Taking Account of the Influence of Drilling Angle on Gas Extraction Effect

The pre-drainage of coalbed methane through boreholes in the bottom drainage roadway (BDR) is the key measure to prevent and control coal and gas outburst. Different arrangement layers in the BDR will make a difference in the range of drilling angle and affect the gas extraction effect. In this paper, the mathematical model of the rock loose circle area around elliptical drilling was constructed. Meanwhile, the fluid–solid coupling model is constructed by using COMSOL software, the dynamic response of coal permeability and volumetric strain with gas pressure and the Klinkenberg effect of gas are considered, and the effect of the change of the elliptical drilling angle on the pressure relief effect of the coal seam is studied. The results showed that the distance between the layer in the BDR and the pre-drainage coal seam would decrease, and the effective extraction length at the same point of gas extraction in the coal seam increases. The area of the rock loose circle and permeability around the drilling decayed negatively and exponentially with the increase in drilling angle. As the drilling angle decreased, the stress in the major axis of the ellipse at the drilling cross-section increased, so the drilling was prone to collapse, and the gas extraction was hindered. Finally, an optimal method of determining the layer in the BDR under the coupling effect of multiple factors was established by combining the measured ground stress. Through field measurement, the drilling extraction rate of the optimized scheme is 60% higher than that of the original scheme.

Numerical simulation of hydraulic fracturing process with consideration of fluid–solid interaction in shale rock

As a critical technology for the efficient exploitation of shale gas, hydraulic fracturing is a two-way fluid-solid coupling interaction process of water-rock. The accurate water pressure calculation in single shale rock fracture during hydraulic fracturing is critical to the fluid-solid coupling analysis. This paper selected the fractured shale rock from the shale reservoir for analysis based on the Saint-Venant principle. Moreover, the water pressure in a single shale rock fracture was calculated by solving the Navier-Stokes equations. Meanwhile, a fluid-solid coupling model of fractured shale rock during hydraulic fracturing was established and evaluated the changing process of hydraulic fracturing. Our results show that in the hydraulic fracturing process, the changing process of water pressure in a single shale fracture with or without fluid-solid coupling is similar, and both have a water hammer effect. Single shale rock fracture opened instantly in an elliptical shape. The maximum instantaneous water pressure at the fracture tip considering fluid-solid coupling was 3.01 times the initial water pressure. The Mises stress at the fracture tip caused by the maximum instantaneous water pressure was 2.97 times the Mises stress at the fracture tip caused by the water pressure in the fracture after stable water pressure change. Mode Ⅰ fracture occurred in the fractured shale rock when high-pressure water rapidly invaded 2/3 of the length of the fracture. The findings in this study may help to understand the hydraulic fracturing mechanism in shale rock.

Simulation of Fracturing and Well Pattern Optimization of Fractured Tight Sandstone Reservoirs

The fracturing and well pattern optimization of fractured reservoirs is a hot spot in petroleum engineering. In this study, taking the tight oil reservoirs of the Yanchang Formation in the Ordos Basin as an example, an optimized plan of hydraulic fracturing and well pattern deployment of the fractured tight reservoir were systematically studied using a fluid-solid coupling model. The results show that the Chang 6 Member of Yanchang Formation mainly develop feldspar sandstone. A large number of vertical fractures are developed in Chang 6 Member, and they have significant shearing properties. The target sandstone reservoir is in the mid-diagenetic A stage. Natural fractures have a significant effect on the extension of hydraulic fractures. When the fracturing direction is parallel to the main strike of natural fractures, the extension distance of the hydraulic fractures is long and the fracturing scale is large; however, when the fracturing direction is orthogonal or oblique to the main strike of natural fractures, the extension distance of the hydraulic fractures is limited. Statistics show that the average length of hydraulic fractures when fracturing parallel to natural fractures is 1.3 times of that when fracturing perpendicular to natural fractures. In addition, after optimization of fracturing parameters, we found that the displacement should be controlled at 2–3 cubic meters/min, and the sand ratio should be controlled at 25%. The smaller the horizontal minimum principal stress difference between the fracturing layer and the shielding layer, the lower the elastic modulus, and the easier it is to control the fracture height. Microseismic monitoring results show that the half-fracture length and fracture height are mainly distributed in 80–140 m and 15–25 m, respectively. The optimized fracture half-length was 120 m and the half-fracture height was 20 m. Finally, the well pattern of the target layer is optimized considering natural fractures. Through this study, we found that the reverse five-spot well pattern is suitable for the efficient development of tight oil reservoirs of the Chang 6 Member in this area.

Similarity Research on the Action of a Centrifuge Scaled Near-Field Underwater Explosion on a Concrete Gravity Dam

Concrete gravity dams are a common form of dams. The dynamic response analysis of concrete gravity dams under the action of underwater explosions is indispensable for survivability assessments. Due to the constraints of the environment and other conditions, an underwater explosion test of a prototype concrete gravity dam is difficult to achieve. The centrifuge scaled-down test that satisfies similarity theory provides a new way to study the dynamic response of underwater explosions on concrete gravity dams. In a near-field explosion, the shock wave of an underwater explosion will cause direct damage to the gravity dam, and the detonation products after the shock wave will aggravate the damage. Under the premise of verifying the concrete dynamic response material parameters and the fluid–solid coupling model of the underwater explosion, numerical calculations of the scaled underwater explosion model of the centrifuge and the prototype underwater explosion were carried out, and the similarity between the scaled dam response model of the underwater explosion centrifuge and the prototype was discussed. The analysis shows that the scaled-down centrifuge model can better reflect the local damage of the prototype dam subjected to shock waves. The difference between the scaled centrifuge model and the prototype model of underwater explosion damage to the dam is concentrated on the overall damage. The damage effect of detonation products on the prototype dam body is mainly concentrated in the middle and lower parts of the dam body, while the damage effect of detonation products on the dam body in the scaled-down model is mainly concentrated in the middle and upper parts of the dam body.

Mechanism of the traumatic brain injury induced by blast wave using the energy assessment method

Although blast-induced traumatic brain injury (b-TBI) is well recognized for its significance in the military population, the unique mechanisms of primary b-TBI remain undefined. The aim of the present study is to unveil the mechanism of the traumatic brain injury induced by the shock wave using the energy assessment method. First, the magnetic resonance images of the human head were processed to establish the finite element (FE) model of the human head including the skin, the skull, the brain and other structures. The simulation of the shock wave was implemented using the Coupled Eulerian-Lagrange (CEL) method, a fluid-solid coupling model of explosive shock wave-head. The coupled model was used to simulate the situations of head subjected to the frontal direction, the side direction and the rear direction by explosive shock wave. The predictions of the FE results were validated by comparing them with previous studies. Then, the degree of b-TBI in three different directions i.e., the frontal direction, the side direction and the rear direction, was assessed to investigate the influence of the direction of the shock wave on the b-TBI. Finally, the various energies including strain energy, kinetic energy, viscous energy and creep energy were extracted from the FE calculations for investigating the mechanism of b-TBI. The results showed that under the scenario of the shock wave, the brain injury was caused by the combined actions of acceleration of the head and the deformation of the skin& facial muscle tissues, the brain and the cerebrospinal fluid (CSF). In addition, it was also found that the dissipation capacity of the head was not sensitive to the impact direction. The present study provides important data and theoretical basis for the mechanism of b-TBI.

Design of an Eccentric Vacuum Preloaded Aerostatic Bearing to Improve the Lateral Stiffness of Guideway

The conventional aerostatic guideway with two relatively large opposed aerostatic bearing pads (OABPs) has large thrust forces at the lateral carriages. The structure is deformed due to the thrust force, and the thickness of the gas film becomes large and inconsistent, resulting in insufficient gas film stiffness. In this study, a fluid–solid coupling model is established for the aerostatic guideway, and the vacuum preloaded (VPL) aerostatic bearing with an eccentric vacuum chamber was obtained through optimization design and computational fluid dynamics (CFD) simulation. The CFD results show that the eccentric VPL aerostatic bearing has a larger stiffness. Then an experiment on the static performance of a lateral carriage with the eccentric VPL aerostatic bearing was carried out. The experimental results show that the lateral stiffness of the aerostatic guideway using the optimized VPL design is about 65 times as large as the previous OABP and the validity of the design optimization is demonstrated.

Analysis on Water Inrush Prevention Mechanism of Paste-Filled Floor above Confined Water

To study the influence of filling step and advancing distance on the deformation and failure of a working face floor, a mechanical model based on elastic foundation beam theory is established. The deflection and bending moment curves of the floor under different filling steps and advancing distance are obtained by Maple. Then, a fluid–solid coupling model of paste-filling mining on confined water is established by FLAC3D. The effects of different filling steps and advancing distance on the floor displacement, stress, and plastic zone of the floor are analyzed. The results show that there is a “concave” quadratic relationship between the filling step and the maximum displacement of the floor, and there is a “convex” quadratic relationship between the advancing distance and the maximum displacement of the floor. The maximum stress of the floor increases linearly with the increase in filling distance and tends to be stable with the increase in advancing distance. Moreover, the increase in filling steps will lead to the continuous increase in longitudinal failure. This study could guide paste-filling mining above confined water.

Evaluation on Distribution Characteristics of Pore Water Pressure within Saturated Pavement Structure Based on the Proposed Tire-Fluid-Pavement Coupling Model

To investigate the flow characteristics of pore water in asphalt pavement and the variation law of the pore water pressure under vehicle loading, a novel method based on BA network and quartet structure generation set method was proposed to reconstruct the three-dimensional (3D) pavement model with pores. The permeability coefficient and the gradation curve were adopted to evaluate the reliability and stability of the random growth pavement model. Then, the tire-fluid-pavement coupling model was established with FLUENT 3D based on the fluid Mie–Gruneisen state equation. According to the built fluid-solid coupling model, the pressure-velocity coupled finite volume algorithm was applied to study the distribution of the pore water pressure in asphalt pavement. Results show that the pore water pressure in asphalt pavement decays periodically with time under vehicle loading. For different types of asphalt pavement, the pore water pressure in open-graded friction course (OGFC) pavement is the smallest during the whole process. Moreover, the peak values of the pore water pressure decrease in the order of asphalt concrete (AC) pavement, stone mastic asphalt (SMA) pavement, and OGFC pavement. The maximum negative value of the pore water pressure is generally less than 0.3 times the maximum positive values. As for saturated pavement pores, the pore water pressure is hardly affected by the water film thickness. The positive peak value of the pore water pressure increases on an approximate parabolic curve as the vehicle speed improves gradually, while the negative one remains largely unchanged. The results are expected to help reduce tire hydroplaning risk and provide guidance for the selection of asphalt mixtures of drainage asphalt pavement.