A Partial Differential Inequality in Geological Models

Sedimentation and erosion processes in sedimentary basins can be modeled by a parabolic equation with a limiter on the fluxes and a constraint on the time variation.This limiter happens to satisfy a stationary scalar hyperbolic inequality,within a constraint,for which the authors prove the existence and the uniqueness of the solution.Actually,this solution is shown to be the maximal element of a convenient convex set of functions.The existence proof is obtained thanks to the use of a numerical scheme.     

Eulerian simulation of sedimentation flows in vertical and inclined vessels

Sedimentation of particles in inclined and vertical vessels is numerically simulated using a finite volume method where the Eulerian multiphase model is applied. The particulate phase as well as the fluid phase is regarded as a continuum while the viscosity and solid stress of the particulate phase are modelled by the kinetic theory of granular flows. The numerical results show an interesting phenomenon of the emergence of two circulation vortices of the sedimentation flow in a vertical vessel but only one in the inclined vessel. Several sensitivity tests are simulated to understand the factors that influence the dual-vortex flow structure in vertical sedimentation. Results show that a larger fluid viscosity makes the two vortex centres much closer to each other and the boundary layer effect at lateral walls is the key factor to induce this phenomenon. In the fluid boundary layer particles settle down more rapidly and drag the local carrier fluid to flow downward near the lateral walls and thus form the dual-vortex flow pattern.     

双颗粒沉降的动力学分析

采用格子Boltzmann-虚拟区域方法对雷诺数范围为50≤Re≤200的双颗粒自由沉降进行了直接数值模拟。首先,研究发现双颗粒最终沉降的位置分别在通道的1/4和3/4位置附近,且与颗粒的初始间距、雷诺数以及通道宽度无关。其次,重点研究了颗粒沉降过程中所受的侧向力(与沉降方向垂直),首次揭示了侧向力的振动频率与雷诺数呈二次关系,且单颗粒的结果始终小于双颗粒的结果;研究还发现侧向力的振动频率与通道宽度近似呈幂律函数关系,且幂律指数与雷诺数有关,雷诺数越大,幂律指数的绝对值越小。最后还研究了雷诺数及通道宽度对侧向力振动振幅的影响。     

MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTA-TION (Ⅰ)-THE CONTINUOUS-TIME CASE

IntroductionTheshallowwaterequationsareanimportantmathematicalmodelforavarietyofprobleminhydraulicengineering .Inrecentyears,therehasbeeninterestinthenumericalsolutionfortheshallowwaterequations.Thenumericalsimulationsfortheshallowwaterequationsystemcanbeappliedtomanypurposes .First,itcanserveasameansformodelingtidalfluctuationsforthosenterestedincapturingtidalenergyforcommercialpurposes.Secondly ,thesesimulationscanbeusedtocomputetidalrangesandsurgessuchashurricanesandtsunamiscausedbyextreme…     

A direct‐forcing pressure‐based lattice Boltzmann method for solving fluid–particle interaction problems

A direct‐forcing immersed boundary‐lattice Boltzmann method (IB–LBM) is developed to simulate fluid–particle interaction problems. This method uses the pressure‐based LBM to solve the incompressible flow field and the immersed boundary method to handle the fluid–particle interactions. The pressure‐based LBM uses the pressure distribution functions instead of the density distribution functions as the independent dynamic variables. The main idea is to explicitly eliminate the compressible effect due to the density fluctuation. In the IB method, a direct‐forcing method is introduced to capture the particle motion. It directly computes an IB force density at each lattice grid from the differences between the pressure distribution functions obtained by the LBM and the equilibrium pressure distribution functions computed from the particle velocity. By applying this direct‐forcing method, the IB–LBM becomes a purely LBM version. Also, by applying the Gauss theorem, the formulas for computing the force and the torque acting on the particle from the flows are derived from the volume integrals over the particle volume instead of from the surface integrals over the particle surface. The order of accuracy of the IB–LBM is demonstrated on the errors of velocity field, wall stress, and gradients of velocity and pressure. As a demonstration of the efficiency and capabilities of the new method, sedimentation of a large number of spherical particles in an enclosure is simulated. Copyright © 2010 John Wiley & Sons, Ltd.     

固体颗粒在热对流下沉降的直接数值模拟研究

在任意拉格朗日-欧拉算法模拟等温惰性颗粒两相流的基础上,增加对能量方程的联立求解,在热对流条件下对固体颗粒在不同雷诺数下的沉降规律进行了直接数值模拟.结果表明:热对流引起了流场流动的变化和不对称,颗粒在热流体中沉降,热对流产生的力加速了冷颗粒的运动,尾部形成了涡脱落;颗粒在冷流体中沉降,热对流产生的力阻碍了冷颗粒的运动,尾部形成了羽流. 随着雷诺数的增大,颗粒经历了稳定沉降、周期性摆动到不稳定、无规则运动3个阶段.      

NUMERICAL STUDY ON THE SEDIMENTED MOTION CHARACTERISTICS OF THREE ALIGNED CIRCULAR PARTICLES IN THE INCLINED CHANNELS

主要应用浸没边界的格子玻尔兹曼方法（immersed boundary-lattice Boltzmann method, IB–LBM） 对处于不同倾斜角度通道内的三个刚体圆形颗粒在重力作用下下落的动力学特性进行了计算研究. 首先分析通道倾斜角度的影响, 结果显示当通道倾斜角处于59°90°的范围时会发生后一个颗粒超越前一个颗粒的现象. 其次, 研究了Re对颗粒沉降特性的影响, 结果表明Re 越大, 颗粒间发生聚集的时间越早. 研究还发现当3 个颗粒的直径大小不均匀时, 颗粒由大到小纵向依次排列, 或者出现中间小球直径较相邻两个小球直径大的排列情况, 均能促使颗粒加快聚集. 本文的研究结果可为环境工程及地质学中的颗粒沉降问题提供有价值的参考.     

微尺度稀薄效应下颗粒沉降的格子Boltzmann方法模拟

基于多松弛格子Boltzmann模型,对竖直细长微通道内颗粒自由沉降过程进行模拟,分析气体稀薄效应、初始位置以及颗粒间相互作用对微颗粒沉降特性的影响.研究表明：随Knudsen数增大,微通道内气体稀薄效应增强,颗粒表面气体滑移速度增大,气相流体有效粘度减小,颗粒相同运动状态下受到气体阻力相应减小,颗粒沉降平衡速度明显增大;不同初始位置颗粒沉降过程存在明显差异,初始位置偏离中心线颗粒将发生水平方向位移且呈振荡趋势,最终稳定于中心线平衡位置;在微尺度双颗粒沉降DKT现象过程中,气体稀薄效应影响颗粒运动特性,后颗粒跟随过程明显增长.