袋型阻尼密封动力学特性双控制体Bulk Flow模型

为快速准确预测袋型阻尼密封泄漏特性和动力学特性,针对传统单控制体Bulk Flow模型预测精度低、无法预测交叉动力系数的问题,提出了袋型阻尼密封双控制体Bulk Flow模型和动力学特性数值预测方法,并开发了计算程序。首先,依据边界层理论,将袋型密封腔室划分为两个控制体,推导了控制体的连续性、周向动量和能量方程,引入Swamee-Jain和Takahashi方程,计算流体-壁面间和流体-流体间的周向黏性摩擦力;其次,采用牛顿-拉夫森算法和摄动分析法分别求解0阶和1阶控制方程,获得各刚度、阻尼动力特性系数;然后,通过与袋型阻尼密封泄漏量和动力特性系数的实验值、单控制体Bulk Flow模型和非定常计算流体动力学(CFD)数值结果进行比较,验证了模型和方法的准确性和可靠性;最后,研究了转子转速(10 000、15 000、20 000 r/min)和预旋比(0.067、0.724、0.997)对袋型阻尼密封动力学特性的影响。结果表明:所发展的模型和方法具有计算速度快、预测精度高(泄漏量预测误差小于6%,动力特性系数预测误差小于38%)的优点;转子转速和进口预旋的增大均会导致袋型阻尼密封有效阻尼显著减小,穿越频率显著增大,易诱发轴系失稳。

肝硬化患者4D flow MRI血流动力学参数与中医证型的相关性研究

目的探讨肝硬化患者4D flow MRI血流动力学参数与中医证型的相关性。方法将118例肝硬化患者依据中医辨证分为肝气郁结证、湿热蕴结证、肝肾阴虚证、脾肾阳虚证、瘀血阻络证5个证型,所有患者均行门静脉4D flow MRI检查,统计不同证型肝硬化患者分布情况,观察门静脉系统(门静脉主干、脾静脉、肠系膜上静脉)的血流动力学参数,包括血流量、流速、壁剪切力等,比较不同证型患者门静脉血流动力学参数差异。结果肝硬化代偿期以肝气郁结证、湿热蕴结证为主,肝硬化失代偿期以脾肾阳虚、瘀血络阻证为主;A级以肝气郁结证、湿热蕴结证为主,B、C级以瘀血络阻证为主。瘀血络阻证肝硬化患者门静脉主干及脾静脉血流量明显高于肝气郁结证、湿热蕴结证患者(P<0.05);脾肾阳虚证门静脉主干血流量明显高于肝气郁结证肝硬化患者(P<0.05);瘀血络阻证肝硬化患者门静脉主干流速及剪切力较肝气郁结证和湿热蕴结证低。结论肝硬化患者中医辨证分型与门静脉血流动力学参数具有一定相关性,4D flow MRI可为肝硬化的中医辨证提供血流动力学参考。

A semi-analytical model for coupled flow in stress-sensitive multi-scale shale reservoirs with fractal characteristics

A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.

Evaluating the stability and volumetric flowback rate of proppant packs in hydraulic fractures using the lattice Boltzmann-discrete element coupling method

The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.

An inverse analysis of fluid flow through granular media using differentiable lattice Boltzmann method

This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.

基于FLOW3D的集成式水下基盘泥沙冲淤三维数值模拟

泥面下集成式水下基盘是为开采渤海通航区等海域油气资源而提出的新型基盘,其基坑周围局部冲淤是工程实践关注的问题之一。基于不可压缩粘性流体运动的Navier-Stokes方程建立泥面下集成式水下基盘基坑周围三维水动力数学模型,对不同粒径和不同流速情况下的局部冲淤进行了模拟。结果表明:泥沙粒径为0.005 mm时,由于泥沙较难起动,基坑附近局部冲淤较小。粒径分别为0.05 mm和0.1 mm时,在典型流速作用下,基盘附近可分别形成1 m左右和4 m左右的淤积。

V Flow技术测量颈动脉壁面剪应力的一致性研究

目的评价血流向量成像(V Flow)技术在测量健康成年人颈动脉壁面剪应力(WSS)中的一致性。方法于2021年2月至2021年3月招募健康成年志愿者20人,由2名不同年资的超声医师使用配备V Flow功能的Mindray Resona 7超声仪和3~9 MHz线阵探头进行双侧颈动脉扫查,分别采集双侧颈总动脉远段、颈总动脉分叉至颈内动脉起始部两段动脉的动态V Flow图像,测量两侧颈总动脉远段的近心端、远心端、分叉处及颈内动脉起始部的前、后壁的WSS,使用组内相关系数(ICC)和Bland-Altman图评估组内一致性及组间一致性。结果双侧颈动脉前、后壁的4个不同节段之间WSS值差异均有统计学意义(P<0.05)。高年资超声医师2次测量结果的一致性结果显示,左侧颈总动脉远段的远心端一致性极好(ICC 0.779),右侧颈总动脉远段的近心端(ICC 0.605)、远心端(ICC 0.585)、颈内动脉起始部(ICC 0.457)、左侧颈总动脉分叉处(ICC 0.606)及颈内动脉起始部(ICC 0.702)一致性均较好;不同年资超声医师的测量结果显示,仅右侧颈总动脉分叉处(ICC 0.486)及左侧颈总动脉远段的远心端(ICC 0.576)一致性较好。结论V Flow技术可显示不同位点间颈动脉WSS之间的差异,其组内一致性较好,但在不同年资超声医师间存在一定的差异。

Performance characteristics of the airlift pump under vertical solid-water-gas flow conditions for conveying centimetric-sized coal particles

In this study,the installation of an airlift pump with inner diameter of 102 mm and length of 5.64 m was utilized to consider the conveying process of non-spherical coal particles with density of 1340 kg/m3 and graining 25-44.5 mm.The test results revealed that the magnitude of increase in the solid transport rate due to the changes in the three tested parameters between compressed air velocity,submergence ratio,and feeding coal possibility was not the same,which are stand in range of 20%,75%,and 40%,respectively.Hence,creating the optimal airlift pump performance is highly dependent on submergence ratio.More importantly,we measured the solid volume fraction using the method of one-way valves in order to minimize the disadvantages of conventional devices,such as fast speed camera and conductivity ring sensor.The results confirmed that the volume fraction of the solid phase in the transfer process was always less than 12%.To validate present experimental data,the existing empirical correlations together with the theoretical equations related to the multiphase flow was used.The overall agreement between the theory and experimental solid delivery results was particularly good instead of the first stage of conveying process.This drawback can be corrected by omitting the role of friction and shear stress at low air income velocity.It was also found that the model developed by Kalenik failed to predict the performance of our airlift operation in terms of the mass flow rate of the coal particles.

Analysis and experimental study on resistance-increasing behavior of composite high efficiency autonomous inflow control device

Bottom water coning is the main reason to reduce the recovery of horizontal bottom water reservoir. By water coning, we mean the oil-water interface changes from a horizontal state to a mound-shaped cone and breaks through to the wellbore. Autonomous inflow control device(AICD) is an important instrument maintain normal production after bottom water coning, however, the resistance increasing ability of the swirl type AICD is insufficient at present, which seriously affects the water control effect. Aiming this problem, this paper designs a multi-stage resistance-increasing and composite type AICD. The separation mechanism of oil-water two phases in this structure, the resistance form of oil-water single phase and the resistance-increasing principle of water phase are analyzed. Establishing the dual-phase multi-stage separation and resistance-increasing model, and verified by measuring the throttling pressure drop and oil-water volume fraction of the AICD, it is found that the composite type AICD has the effect of ICD and AICD at the same time, which can balance the production rate of each well section at the initial stage of production, delay the occurrence of bottom water coning. In the middle and later stages of production, water-blocking can be effectively increased to achieve water control and stable production.After structural sensitivity analysis, the influence law of various structural parameters on the water control performance of composite AICD was obtained. The simulation calculation results show that,compared with the existing swirl type AICD, composite AICD has higher sensitivity to moisture content,the water phase throttling pressure drop is increased by 4.5 times on average. The composite AICD is suitable for the entire stage of horizontal well production.

A transient production prediction method for tight condensate gas wells with multiphase flow

Considering the phase behaviors in condensate gas reservoirs and the oil-gas two-phase linear flow and boundary-dominated flow in the reservoir,a method for predicting the relationship between oil saturation and pressure in the full-path of tight condensate gas well is proposed,and a model for predicting the transient production from tight condensate gas wells with multiphase flow is established.The research indicates that the relationship curve between condensate oil saturation and pressure is crucial for calculating the pseudo-pressure.In the early stage of production or in areas far from the wellbore with high reservoir pressure,the condensate oil saturation can be calculated using early-stage production dynamic data through material balance models.In the late stage of production or in areas close to the wellbore with low reservoir pressure,the condensate oil saturation can be calculated using the data of constant composition expansion test.In the middle stages of production or when reservoir pressure is at an intermediate level,the data obtained from the previous two stages can be interpolated to form a complete full-path relationship curve between oil saturation and pressure.Through simulation and field application,the new method is verified to be reliable and practical.It can be applied for prediction of middle-stage and late-stage production of tight condensate gas wells and assessment of single-well recoverable reserves.

Numerical and experimental study on the falling film flow characteristics with the effect of co-current gas flow in hydrogen liquefaction process

Liquid hydrogen storage and transportation is an effective method for large-scale transportation and utilization of hydrogen energy. Revealing the flow mechanism of cryogenic working fluid is the key to optimize heat exchanger structure and hydrogen liquefaction process(LH2). The methods of cryogenic visualization experiment, theoretical analysis and numerical simulation are conducted to study the falling film flow characteristics with the effect of co-current gas flow in LH2spiral wound heat exchanger.The results show that the flow rate of mixed refrigerant has a great influence on liquid film spreading process, falling film flow pattern and heat transfer performance. The liquid film of LH2mixed refrigerant with column flow pattern can not uniformly and completely cover the tube wall surface. As liquid flow rate increases, the falling film flow pattern evolves into sheet-column flow and sheet flow, and liquid film completely covers the surface of tube wall. With the increase of shear effect of gas-phase mixed refrigerant in the same direction, the liquid film gradually becomes unstable, and the flow pattern eventually evolves into a mist flow.

Depositional process of hyperpycnal flow deposits:A case study on Lower Cretaceous Sangyuan outcrop in the Luanping Basin,Northeast China

Sedimentary process research is of great significance for understanding the distribution and characteristics of sediments.Through the detailed observation and measurement of the Sangyuan outcrop in Luanping Basin,this paper studies the depositional process of the hyperpycnal flow deposits,and divides their depositional process into three phases,namely,acceleration,erosion and deceleration.In the acceleration phase,hyperpycnal flow begins to enter the basin nearby,and then speeds up gradually.Deposits developed in the acceleration phase are reverse.In addition,the original deposits become unstable and are taken away by hyperpycnal flows under the eroding force.As a result,there are a lot of mixture of red mud pebbles outside the basin and gray mud pebbles within the basin.In the erosion phase,the reverse deposits are eroded and become thinner or even disappear.Therefore,no reverse grading characteristic is found in the proximal major channel that is closer to the source,but it is still preserved in the middle branch channel that is far from the source.After entering the deceleration phase,normally grading deposits appear and cover previous deposits.The final deposits in the basin are special.Some are reverse,and others are normal.They are superimposed with each other under the action of hyperpycnal flow.The analysis of the Sangyuan outcrop demonstrates the sedimentary process and distribution of hyperpycnites,and reasonably explain the sedimentary characteristics of hyperpycnites.It is helpful to the prediction of oil and gas exploration targets in gravity flow deposits.

Extending homogeneous fluidization flow regime of Geldart-A particles by exerting axial uniform and steady magnetic field

The homogeneous/particulate fluidization flow regime is particularly suitable for handling the various gas–solid contact processes encountered in the chemical and energy industry.This work aimed to extend such a regime of Geldart-A particles by exerting the axial uniform and steady magnetic field.Under the action of the magnetic field,the overall homogeneous fluidization regime of Geldart-A magnetizable particles became composed of two parts:inherent homogeneous fluidization and newly-created magnetic stabilization.Since the former remained almost unchanged whereas the latter became broader as the magnetic field intensity increased,the overall homogeneous fluidization regime could be extended remarkably.As for Geldart-A nonmagnetizable particles,certain amount of magnetizable particles had to be premixed to transmit the magnetic stabilization.Among others,the mere addition of magnetizable particles could broaden the homogeneous fluidization regime.The added content of magnetizable particles had an optimal value with smaller/lighter ones working better.The added magnetizable particles might raise the ratio between the interparticle force and the particle gravity.After the magnetic field was exerted,the homogeneous fluidization regime was further expanded due to the formation of magnetic stabilization flow regime.The more the added magnetizable particles,the better the magnetic performance and the broader the overall homogeneous fluidization regime.Smaller/lighter magnetizable particles were preferred to maximize the magnetic performance and extend the overall homogeneous fluidization regime.This phenomenon could be ascribed to that the added magnetizable particles themselves became more Geldart-A than-B type as their density or size decreased.

A multichannel thermal bubble-actuated impedance flow cytometer with on-chip TIA based on CMOS-MEMS

Electrochemical impedance spectroscopy(EIS)flow cytometry offers the advantages of speed,affordability,and portability in cell analysis and cytometry applications.However,the integration challenges of microfluidic and EIS read-out circuits hinder the downsizing of cytometry devices.To address this,we developed a thermal-bubble-driven impedance flow cytometric application-specific integrated circuit(ASIC).The thermal-bubble micropump avoids external piping and equipment,enabling high-throughput designs.With a total of 36 cell counting channels,each measuring 884×220μm^(2),the chip significantly enhances the throughput of flow cytometers.Each cell counting channel incorporates a differential trans-impedance amplifier(TIA)to amplify weak biosensing signals.By eliminating the parasitic parameters created at the complementary metal-oxidesemiconductor transistor(CMOS)-micro-electromechanical systems(MEMS)interface,the counting accuracy can be increased.The on-chip TIA can adjust feedback resistance from 5 to 60 kΩto accommodate solutions with different impedances.The chip effectively classifies particles of varying sizes,demonstrated by the average peak voltages of 0.0529 and 0.4510 mV for 7 and 14μm polystyrene beads,respectively.Moreover,the counting accuracies of the chip for polystyrene beads and MSTO-211H cells are both greater than 97.6%.The chip exhibits potential for impedance flow cytometer at low cost,high-throughput,and miniaturization for the application of point-of-care diagnostics.

Numerical simulation on the multiphase flow and reoxidation of the molten steel in a two-strand tundish during ladle change

A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.

Effect of aspect ratio of elliptical stirred vessel on mixing time and flow field characteristics in the absence of baffles

Elliptical tanks were used as an alternative to circular tanks in order to improve mixing efficiency and reduce mixing time in unbaffled stirred tanks(USTs). Five different aspect ratios of elliptical vessels were designed to compare their mixing time and flow field. Computational fluid dynamics(CFD) simulations were performed using the k–ε model to calculate the mixing time and simulate turbulent flow field features, such as streamline shape, velocity distribution, vortex core region distribution, and turbulent kinetic energy(TKE) transfer. Visualization was also carried out to track the tinctorial evolution of the liquid phase. Results reveal that elliptical stirred tanks can significantly improve mixing performance in USTs. Specifically, the mixing time at an aspect ratio of 2.00 is only 45.3% of the one of a circular stirred tank. Furthermore, the secondary flow is strengthened and the vortex core region increases with the increase of aspect ratio. The axial velocity is more sensitive to the aspect ratio than the circumferential and radial velocity. Additionally, the TKE transfer in elliptical vessels is altered. These findings suggest that elliptical vessels offer a promising alternative to circular vessels for enhancing mixing performance in USTs.

Debris flow runout behaviors considering the influences of densely populated buildings

Debris flows pose serious risks to communities in mountainous areas,often resulting in large losses of human life and property.The impeding presence of urban buildings often affects the runout behavior and deposition of debris flows.But the impact of different building densities and sizes on debris flow dynamics has yet to be quantified to guide urban planning in debris flow risk zones.This study focused on a debris flow that occurred in Zhouqu County,Gansu Province,China on August 7th,2010,which was catastrophic and destroyed many buildings.The FLO-2D software was used to simulate this debris flow in two scenarios,i.e.the presence and the absence of buildings,to obtain debris-flow intensity parameters.The developed model was then used to further analyze the influence of large buildings and narrow channels within the urban environment.The simulation results show that considering the presence of buildings in the simulation is essential for accurate assessment of debris flow intensity and deposition distribution.The layout of buildings in the upstream urban area,such as large buildings or parallel buildings which form narrow channels,can affect the flow velocity and depth of debris flow heading towards downstream buildings.To mitigate damage to downstream buildings,the relative spacing(d/a)between upstream and downstream buildings should not exceed a value of two and should ideally be even lower.These findings provide valuable insights for improving the resistance of mountainous cities to urban debris flows.

Enhancing wood efficiency through comprehensive wood flow analysis:Methodology and strategic insights

Wood,an essential natural resource in human civilization,remains widely used despite advances in technology and material substitution.The surge in greenhouse gas emissions and environmental concerns accentuates the need for optimizing wood utilization.Material flow analysis is a powerful tool for tracking material flows and stocks,aiding resource management and environmental decision-making.However,the full extent of its methodological dimensions,particularly within the context of the wood supply chain,remains relatively unexplored.In this study,we delve into the existing literature on wood flow analysis,discussing its primary objectives,materials involved,temporal and spatial scales,data sources,units,and conversion factors.Additionally,data uncertainty,data reconciliation and crucial assumptions in material flow analysis are highlighted in this paper.Key findings reveal the significance of wood cascading and substitution effects by replacing non-wood materials,where they can reduce greenhouse gas emissions more than the natural carbon sink of forests and wood products.The immediate impact of short-term wood cascading might not be as robust as the substitution effect,with energy substitution showcasing better results than material substitution.However,it's crucial to note that these conclusions could experience significant reversal from a long-term and global perspective.Strategies for improving wood efficiency involve maximizing material use,advancing construction technologies,extending product lifespans,promoting cascade use,and optimizing energy recovery processes.The study underscores the need for standardized approaches in wood flow analysis and emphasizes the potential of wood efficiency strategies in addressing environmental challenges.

Protective effects of baffles with different positions,row spacings,heights on debris flow impact

The baffle effectively slowed down debris flow velocity,reduced its kinetic energy,and significantly shortened the distance of debris flow movement.Consequently,they are widely used for protection against natural hazards such as landslides and mudslides.This study,based on the threedimensional DEM(Discrete Element Method),investigated the impact of different baffle positions on debris flow protection.Debris flow velocity and kinetic energy variations were studied through single-factor experiments.Suitable baffle positions were preliminarily selected by analyzing the influence of the first-row baffle position on the impact force and accumulation mass of debris flow.Subsequently,based on the selected baffle positions and four factors influencing the effectiveness of baffle protection(baffle position(P),baffle height(h),row spacing(S_(r)),and angle of transit area(α)),an orthogonal design was employed to further explore the optimal arrangement of baffles.The research results indicate that the use of a baffle structure could effectively slow down the motion velocity of debris flows and dissipate their energy.When the baffle is placed in the transit area,the impact force on the first-row baffle is greater than that when the baffle is placed in the deposition area.Similarly,when the baffle is placed in the transit area,the obstruction effect on debris flow mass is also greater than that when the baffle is placed in the deposition area.Through orthogonal experimental range analysis,when the impact on the first row of baffles is used as the evaluation criterion,the importance of each influencing factor is ranked asα>P>S_(r)>h.When the mass of debris flow behind the baffle is regarded as the evaluation criterion,the rank is changed to P>α>S_(r)>h.The experimental simulation results show that the optimal baffle arrangement is:P_(5),S_(r)=16,α=35°,h=9.

A new interacting capillary bundle model on the multiphase flow in micropores of tight rocks

Surfactants are widely used in the fracturing fluid to enhance the imbibition and thus the oil recovery rate. However, current numerical models cannot capture the physics behind capillary imbibition during the wettability alteration by surfactants. Although the interacting capillary bundle(ICB) model shows potential in characterizing imbibition rates in different pores during wettability alteration, the existing ICB models neglect the influence of wettability and viscosity ratio on the imbibition behavior, making it difficult to accurately describe the oil-water imbibition behavior within the porous media. In this work,a new ICB mathematical model is established by introducing pressure balance without assuming the position of the leading front to comprehensively describe the imbibition behavior in a porous medium under different conditions, including gas-liquid spontaneous imbibition and oil-water imbibition.When the pore size distribution of a tight rock is known, this new model can predict the changes of water saturation during the displacement process in the tight rock, and also determine the imbibition rate in pores of different sizes. The water saturation profiles obtained from the new model are validated against the waterflooding simulation results from the CMG, while the imbibition rates calculated by the model are validated against the experimental observations of gas-liquid spontaneous imbibition. The good match above indicates the newly proposed model can show the water saturation profile at a macroscopic scale while capture the underlying physics of the multiphase flow in a porous medium at a microscopic scale. Simulation results obtained from this model indicate that both wettability and viscosity ratio can affect the sequence of fluid imbibition into pores of different sizes during the multiphase flow, where less-viscous wetting fluid is preferentially imbibed into larger pores while more-viscous wetting fluid tends to be imbibed into smaller pores. Furthermore, this model provides an avenue to calculate the imbibition rate in pores of different sizes during wettability alteration and capture the non-Darcy effect in micro-and nano-scale pores.