凡凤仙a，b， 陈超凡a， 罗伟乐a
上海理工大学 a. 能源与动力工程学院， b. 中国机械工业多相流测试技术重点实验室， 上海 200093
引用格式：
凡凤仙， 陈超凡， 罗伟乐. 初始位置对圆管泊肃叶流中颗粒惯性迁移的影响［J］. 中国粉体技术， 2026， 32（5）： 1-12.
Fan Fengxian， Chen Chaofan， Luo Weile. Influence of initial position on inertial migration of particles in circular pipe Poiseuille flow［J］. China Powder Science and Technology， 2026， 32（5）： 1-12.
DOI：10.13732/j.issn.1008-5548.2026.05.002
收稿日期： 2026-03-21， 修回日期： 2026-05-10，上线日期： 2026-06-18。
基金项目：国家自然科学基金项目，编号：52476157；中国机械工业多相流测试技术重点实验室基金，编号：2024SA-05-22。
第一作者：凡凤仙（1982—），女，博士，教授，博士生导师，研究方向为多相流动与传热。E-mail：fanfengxian@usst.edu.cn。
摘要：【目的】探究初始位置对圆管泊肃叶流中单个中性浮力球形颗粒运动的影响，以准确预测颗粒全过程动力学行为。【方法】采用计算流体力学与离散元方法，建立数学模型，描述流体和颗粒的运动方程以及流体和颗粒的相互作用方程，在管道进口施加充分发展的泊肃叶流速度边界条件，设置出口为零压力边界，管道壁面无滑移；开展网格无关性验证与模型验证，验证网格的高精度和模型的可靠性，从4种升力作用下分析典型雷诺数下颗粒的惯性迁移行为及其机制。【结果】随着雷诺数的增大，颗粒迁移轨迹呈现单调、过冲、振荡3种模式；颗粒的径向平衡位置与初始位置无关，而随着雷诺数的增大先向壁面而后向中心轴线移动。初始位置偏离平衡位置时，颗粒在单调区的径向速度先增大后减小，在过冲区的径向速度先减小至0而后反向增大，继而再减小至0；在振荡区呈衰减振荡特性。平均径向迁移速度由初始位置和雷诺数共同决定；颗粒旋转速度的演变规律与径向迁移轨迹一致，最终旋转速度随着雷诺数的增大而增大，且与初始位置无关。【结论】雷诺数主导迁移模式与平衡位置，初始位置决定迁移动态特性。
关键词：液固两相流；泊肃叶流；惯性迁移；颗粒动力学；计算流体力学-离散元方法
Abstract
Objective The inertial migration of particles in circular pipe Poiseuille flow （Segré-Silberberg effect） holds important application value in microfluidic technologies such as particle enrichment， separation， and detection. Existing studies have mainly focused on the final radial equilibrium position of particles， while the understanding of the dynamic characteristics throughout the entire migration process （such as the evolution patterns of trajectory， migration velocity， and rotational velocity） and their dependence on the initial position remains insufficient. To investigate the influence mechanism of the initial position on the inertial migration dynamics of particles， this study conducts numerical simulations of the motion of a single neutrally buoyant spherical particle in a circular pipe Poiseuille flow.
Methods Using the fictitious domain method based on the coupling of computational fluid dynamics and the discrete element method， a three-dimensional circular pipe flow model was established. The pipe had a diameter of 8 mm and a length of 0.5 m. The particle diameter was 1.3 mm. Both the particle and fluid densities were 1 050 kg·m-3， and the kinematic viscosity of the fluid was 1.428 6×10-6 m²·s-1. A fully developed Poiseuille-flow velocity boundary was imposed at the pipe inlet， a zero-pressure boundary was set at the outlet， and a no-slip condition was applied at the pipe wall. The initial radial positions of the particle covered three cases： between the equilibrium position and the pipe center， near the equilibrium position， and between the equilibrium position and the pipe wall. The reliability of the numerical model was validated through mesh independence verification （using 655 583 static meshes with three-level dynamic refinement in the particle domain） and by comparison with results from the literature. Simulations were conducted under Reynolds numbers Re=200， 400， 600， and 800 to analyze the time evolution and final equilibrium characteristics of particle radial position， radial migration velocity， and rotational velocity.
Results and Discussion As the Reynolds number increased， the inertial migration trajectory of particles exhibited three typical modes： monotonic， overshooting， and oscillatory migration. At Re=200， particles monotonically approached the radial equilibrium position. At Re=400， overshooting occurred， and particles returned after crossing the equilibrium position. When Re exceeded 600， particles oscillated with decaying amplitude after reaching the equilibrium position and eventually became stabilized. The radial equilibrium position of particles was independent of the initial position. However， with increasing Reynolds number， the equilibrium position first moved toward the pipe wall （Re=200-400） and then toward the centerline （Re = 400-800）. The dimensionless equilibrium position r*/R varied from 0.687 to 0.716， with a change rate of less than 2.6%. In the monotonic migration regime （Re=200）， when the initial position deviated from the equilibrium position， the radial migration velocity first increased and then decreased to 0. In the overshooting regime （Re=400）， the radial velocity first decreased to 0， then increased in the opposite direction， and finally decreased to 0 again. In the oscillatory regime （Re=600 and 800）， the radial velocity showed a decaying oscillation. The average radial migration velocity for particles to reach the equilibrium position for the first time was jointly affected by the initial position and the Reynolds number. When the initial position was close to the pipe center， the average velocity increased monotonically with Re. However， when the initial position was close to the pipe wall， the average velocity at Re=500 was significantly lower than that at Re=400 because oscillatory migration at high Reynolds numbers prolonged the time required to first reach the equilibrium position. The temporal evolution patterns of the particle rotational velocity were highly consistent with the radial migration trajectory， and the final rotational velocity increased monotonically with the Reynolds number and was independent of the initial position. The fundamental reason was that particle rotation was dominated by the velocity gradient of the local flow field， and this gradient increased with Reynolds number.
Conclusion The Reynolds number governs the inertial migration mode and final equilibrium position of particles， whereas the initial position determines the migration dynamic characteristics （such as the evolution path of radial velocity and the average migration velocity）. The findings provide a theoretical basis for achieving precise regulation and efficient separation of particles through inertial migration.
Keywords： liquid-solid two-phase flow； Poiseuille flow； inertial migration； particle dynamics； computational fluid dynamics-discrete element method （CFD-DEM）
参考文献（References）
［1］叶挺， 凡凤仙. 泊肃叶流下中性浮力颗粒横向迁移的研究进展［J］. 上海理工大学学报， 2021， 43（2）： 118-126.
Ye Ting， Fan Fengxian. Research progress on the lateral migration of neutrally buoyant particles in a Poiseuille flow［J］. Journal of University of Shanghai for Science and Technology， 2021， 43（2）： 118-126.
［2］彭德其， 刘建， 付翁， 等. 缩放管-液固两相流换热器的流场特征影响因素［J］. 中国粉体技术，2025，31（3）： 67-77.
Peng Deqi， Liu Jian， Fu Weng， et al. Influencing factors of flow field characteristics of liquid-solid two-phase flow in converging-diverging tube heat exchangers［J］. China Powder Science and Technology， 2025， 31（3）： 67-77.
［3］张鑫镝， 张卫义， 李兆亭， 等. 液-固流化床最小流态化速度的数值模拟与验证［J］. 中国粉体技术， 2023， 29（1）： 80-89.
Zhang Xindi， Zhang Weiyi， Li Zhaoting， et al. Numerical simulation and verification of the minimum fluidization velocity in a liquid-solid fluidized bed［J］. China Powder Science and Technology， 2023， 29（1）： 80-89.
［4］Huang L R， Cox E C， Austin R H， et al. Continuous particle separation through deterministic lateral displacement［J］. Science， 2004， 304（5673）： 987-990.
［5］Di Carlo D， Irimia D， Tompkins R G， et al. Continuous inertial focusing， ordering， and separation of particles in microchannels［J］. Proceedings of the National Academy of Sciences， 2007， 104（48）： 18892-18897.
［6］Park W， Lee J， Kim M. Recent progress of inertial microfluidics for high-throughput and scalable separation of microscale bioparticles［J］. Sensors and Actuators A： Physical， 2025， 382： 116112.
［7］Segré G， Silberberg A. Radial particle displacements in poiseuille flow of suspensions［J］. Nature， 1961， 189（4760）： 209-210.
［8］Matas J P， Morris J F， Guazzelli É. Inertial migration of rigid spherical particles in Poiseuille flow［J］. Journal of Fluid Mechanics， 2004， 515： 171-195.
［9］Morita Y， Itano T， Sugihara-Seki M. Equilibrium radial positions of neutrally buoyant spherical particles over the circular cross-section in Poiseuille flow［J］. Journal of Fluid Mechanics， 2017， 813： 750-767.
［10］Aanad P， Subramanian G. Inertial migration in a pressure-driven channel flow： beyond the segre-silberberg pinch［J］. Physical Review Letters， 2024， 132（5）： 054002.
［11］Luo Weile， Parteli E J R， Pöschel T， et al. Inertial migration regimes of a neutrally buoyant sphere in pipe Poiseuille flow［J］. Journal of Fluid Mechanics， 2026， 1029： A53.
［12］Jeffrey R C， Pearson J R A. Particle motion in laminar vertical tube flow［J］. Journal of Fluid Mechanics， 1965， 22（4）： 721-735.
［13］Karnis A， Goldsmith H L， Mason S G. The flow of suspensions through tubes： V. Inertial effects［J］. The Canadian Journal of Chemical Engineering， 1966， 44（4）： 181-193.
［14］Nakayama S， Yamashita H， Yabu T， et al. Three regimes of inertial focusing for spherical particles suspended in circular tube flows［J］. Journal of Fluid Mechanics， 2019， 871： 952-969.
［15］Wang Jing， Joseph D D. Lift forces on a cylindrical particle in plane Poiseuille flow of shear thinning fluids［J］. Physics of Fluids， 2003， 15（8）： 2267-2278.
［16］Shao Xueming， Yu Zhaosheng， Sun Bo. Inertial migration of spherical particles in circular Poiseuille flow at moderately high Reynolds numbers［J］. Physics of Fluids， 2008， 20（10）： 103307.
［17］Pan T W， Li A， Glowinski R. Numerical study of equilibrium radial positions of neutrally buoyant balls in circular Poiseuille flows［J］. Physics of Fluids， 2021， 33（3）： 033301.
［18］Aouane O， Sega M， Bäuerlein B， et al. Inertial focusing of a dilute suspension in pipe flow［J］. Physics of Fluids， 2022， 34（11）： 113312.
［19］Liu Wenwei， Wu Chuanyu. Inertial migration of a neutrally buoyant circular particle in a planar Poiseuille flow with thermal fluids［J］. Physics of Fluids， 2021， 33（6）： 063315.
［20］Prohm C， Zöller N， Stark H. Controlling inertial focussing using rotational motion［J］. The European Physical Journal E， 2014， 37（5）： 36.
［21］Alexeev D， Litvinov S， Economides A， et al. Inertial focusing of spherical particles： the effects of rotational motion［J］. Physical Review Fluids， 2025， 10（5）： 054202.
［22］Ren Wanlong， Zhang Xuhui， Zhang Yan， et al. Investigation of motion characteristics of coarse particles in hydraulic collection［J］. Physics of Fluids， 2023， 35（4）： 043322.
［23］Li Chengzhi， Zhang Yu， Shen Jiahe， et al. Coupled simulation of fluid-particle interaction for large complex granules： a resolved CFD-DEM method for modelling the airflow in a vertical fixed bed of irregular sinter particles［J］. Particuology， 2024， 90： 292-306.
［24］Hager A， Kloss， C， Pirker， S， et al. Parallel resolved open source CFD-DEM： method， validation and application［J］. The Journal of Computational Multiphase Flows， 2014， 6（1）： 13-27.
［25］Kloss C， Goniva. C， Hager A，et al. Models， algorithms and validation for opensource DEM and CFD-DEM［J］. Progress in Computational Fluid Dynamics， an International Journal， 2012， 12（2/3）： 140-152.
