黏弹性流体扩展及收缩流动的格子Boltzmann模拟分析

本文应用格子Boltzmann 方法(LBM)并结合Oldroyd-B 模型,讨论了不可压缩的 Navier-Stokes 方程和平流扩散本构方程的解耦及各自求解方法,以及两类问题的边界处理格式,实现了黏弹性流体在二维1:3 扩展流道以及3:1 收缩流道中的流动的数值模拟．获得了不同雷诺数Re 和维森伯格数Wi 以及黏度vs 下流动的流线分布,计算给出了漩涡的涡心位置和大小,并分析了参数Re、Wi 和vs 对流动特点的影响．模拟结果表明本文所采用模型和边界处理方法具有良好的精度和稳定性．     

基于格子Boltzmann两相模型的粘弹流体挤出胀大分析

挤出胀大的数值模拟是非牛顿流体研究中具有挑战性的问题．本文运用格子Boltzmann方法(LBM)分析Oldroyd-B和多阶松弛谱PTT粘弹流体的挤出胀大现象，采用颜色模型模拟出口处粘弹流体和空气的两相流动，通过重新标色获得两种流体的界面，并最终获得胀大的形状．Navier-Stokes方程和本构方程的求解采用双分布函数模型．将胀大的结果与解析解、实验解和单相自由面LBM结果进行了比较，发现格子Boltzmann两相模型结果与解析解和实验结果相吻合，相比于单相模型，收敛速度更快，解的稳定性更高．研究了流道尺寸对胀大率的影响，并对挤出胀大的内在机理进行了分析．     

可压缩流体中封闭薄球壳的自由振动

通过引入位移函数，成功地研究了薄球壳在可压缩流体中的自由振动。发现存在两类自由振动：第一类与外界流体无关；第二类则受到流体性质的影响．证明了第二类振动的频率方程具有多项式形式并只存在复频率（除ｎ＝１时有Ω＝０）．求解ｎ＝０，１和２时的频率方程，并讨论参数的影响及给出根轨迹图。最后就小阻尼系数法作了对比分析。     

近岸非平整海底上耗散动力系统的广义波作用量守恒方程

为刻划近岸波-流-海底相互作用耗散动力系统的多种复杂作用机制，着眼于波浪对近岸大尺度变化环境流作用和考虑多变海底地形(可典型地刻划为由慢变水深和快变水深构成)的影响，由基于黏性流体Navie-Stokes方程的平均流方程，建立了近岸耗散动力系统的广义波作用量守恒方程，从中提出垂向速度波作用量和耗散波作用量这两种新概念，使得它们和经典的波作用量相互间达成了一种互补、协调而又主次分明的更为广泛的守恒形式．从而把波作用量这一经典概念从理想的平均流守恒系统引申到实际的平均流耗散系统(即广义守恒系统)中去，为解释沿岸过程和应用于近海、海岸工程提供了一个理论基础．     

直角和贴体坐标系下两种格子Boltzmann方法的比较研究

采用格子Boltzmann方法(LBM)和改进的插值格子Boltzmann方法(GILBM)研究了45°斜方腔的顶盖驱动流和Roach通道内的流动特性,并与基准解进行了对比。结果表明,对于45°斜方腔的顶盖驱动流,当雷诺数较小时,两种方法的计算结果与基准解吻合较好;但当雷诺数较大时,采用LBM的计算结果准确性降低,而基于GILBM方法得到的结果准确度升高,且计算稳定性好。对于Roach通道内的流体流动而言,两种方法的计算精度和复杂边界的复杂程度与雷诺数大小有关。根据流场边界形状的复杂程度,网格划分与计算精确度的不同要求,两种方法各有利弊。     

非Четаев型非完整系统的Lie对称性与守恒量

研究非Четаев型非完整系统的Ｌｉｅ对称性．首先利用微分方程在无限小变换下的不变性建立Ｌｉｅ对称所满足的确定方程和限制方程，给出结构方程并求出守恒量；其次研究上述问题的逆问题：根据已知积分求相应的Ｌｉｅ对称性；最后举例说明结果的应用．     

基于格子Boltzmann方法的挤出胀大的数值模拟

将单相格子Boltzmann方法(lattice Boltzmann method, LBM)引入到粘弹流体的瞬态挤出胀大的数值模拟中,建立了基于双分布函数的自由面粘弹性流动格子Boltzmann模型．分析得到的流道中流动速度分布和构型张量结果与理论解十分吻合．对粘弹流体瞬态挤出胀大过程进行了模拟,并分析了运动粘度比和剪切速率对挤出胀大率的影响,得到的胀大率结果与理论分析和其它模拟结果基本一致．表明给出的LBM可以捕捉挤出胀大的瞬态效应．     

颗粒在黏弹性流体扩张和收缩流中沉降的格子Boltzmann模拟分析

对于Oldroyd-B型黏弹性流体，本文应用格子Boltzmann方法(LBM)，实现了流体在二维1:3扩张流道及3:1收缩流道中流动的数值模拟，获得了黏弹性流体在扩张和收缩流道中的流场分布．结合颗粒的受力和运动规则，基于点源颗粒模型，数值分析了颗粒在扩张流和收缩流中的沉降过程和特征，讨论了颗粒相对质量和起始位置以及雷诺数Re和威森伯格数Wi对颗粒沉降特征的影响．结果表明，颗粒相对质量和起始位置以及Re对颗粒沉降轨迹和落点影响较大，而Wi的影响则较小．     

基于分数阶模型的非完整系统的Mei对称性及其守恒量

为了进一步揭示非完整系统的对称性和守恒量之间的内在关系，提出并研究基于分数阶模型的非完整系统的Mei对称性及其守恒量．首先，根据分数阶d’Alembert-Lagrange原理建立基于分数阶模型的非完整系统的动力学方程．其次，根据动力学方程中的动力学函数经无限小变换后仍满足原方程的不变性，建立分数阶模型下非完整系统的Mei对称性定理，给出Mei守恒量．再次，讨论了几个特例：分数阶Hamilton系统、经典非完整系统和受非完整约束的分数阶Lagrange系统的Mei对称性定理．文末举例说明结果的应用．     

柔性微粒介电泳分离过程的多尺度模拟

介电泳分离是一种高效的微细颗粒分离技术,利用非均匀电场极化并操纵分离微流道中的颗粒. 柔性微粒在介电泳分离过程中同时受多种物理场、多相流和微粒变形等复杂因素的影响,仅用单一的计算方法对其进行模拟存在一定的难度,本文采用有限单元——格子玻尔兹曼耦合计算的方法处理这一难题.介观尺度的格子玻尔兹曼方法将流体看成由大量微小粒子组成,在离散格子上求解玻尔兹曼输运方程,易于处理多相流及大变形问题,特别适合模拟柔性颗粒在介电泳分离过程中的变形情况.另一方面,介电泳分离过程的模拟需求解流体、电场和微粒运动方程,计算量相当庞大,通过有限单元法求解介电泳力,提高计算效率.利用这种多尺度耦合计算方法,对一款现有的介电泳芯片分离过程进行了模拟.分析了微粒在电场作用下产生的介电泳力,揭示了介电泳力与电场变化率等因素之间的关系.对微粒运动轨迹及其变形的情况进行了研究,发现微粒的变形主要与流体剪切作用有关.这种多尺度耦合计算方法,为复杂微流体的计算提供了一种有效的解决方案.      

Lattice Boltzmann simulations of turbulent shear flow between parallel porous walls

The effects of two parallel porous walls are investigated, consisting of the Darcy number and the porosity of a porous medium, on the behavior of turbulent shear flows as well as skin-friction drag. The turbulent channel flow with a porous surface is directly simulated by the lattice Boltzmann method （LBM）. The Darcy-Brinkman- Forcheimer （DBF） acting force term is added in the lattice Boltzmann equation to simu- late the turbulent flow bounded by porous walls. It is found that there are two opposite trends （enhancement or reduction） for the porous medium to modify the intensities of the velocity fluctuations and the Reynolds stresses in the near wall region. The parametric study shows that flow modification depends on the Darcy number and the porosity of the porous medium. The results show that, with respect to the conventional impermeable wall, the degree of turbulence modification does not depend on any simple set of param- eters obviously. Moreover, the drag in porous wall-bounded turbulent flow decreases if the Darcy number is smaller than the order of O（10-4） and the porosity of porous walls is up to 0.4.     

Lattice Boltzmann method for simulating particle–fluid interactions

The lattice Boltzmann method （LBM） has gained increasing popularity in the last two decades as an alternative numerical approach for solving fluid flow problems. One of the most active research areas in the LBM is its application in particle-fluid systems, where the advantage of the LBM in efficiency and parallel scalability has made it superior to many other direct numerical simulation （DNS） techniques. This article intends to provide a brief review of the application of the LBM in particle-fluid systems. The numerical techniques in the LBM pertaining to simulations of particles are discussed, with emphasis on the advanced treatment for boundary conditions on the particle-fluid interface. Other numerical issues, such as the effect of the internal fluid, are also briefly described. Additionally, recent efforts in using the LBM to obtain closures for particle-fluid drag force are also reviewed.     

Lattice Boltzmann simulations of high-order statistics in isotropic turbulent flows

The lattice Boltzmann method (LBM) is coupled with the multiple-relaxation-time (MRT) collision model and the three-dimensional 19-discrete-velocity (D3Q19) model to resolve intermittent behaviors on small scales in isotropic turbulent flows. The high-order scaling exponents of the velocity structure functions, the probability distribution functions of Lagrangian accelerations, and the local energy dissipation rates are investigated. The self-similarity of the space-time velocity structure functions is explored using the extended self-similarity (ESS) method, which was originally developed for velocity spatial structure functions. The scaling exponents of spatial structure functions at up to ten orders are consistent with the experimental measurements and theoretical results, implying that the LBM can accurately resolve the intermittent behaviors. This validation provides a solid basis for using the LBM to study more complex processes that are sensitive to small scales in turbulent flows, such as the relative dispersion of pollutants and mesoscale structures of preferential concentration of heavy particles suspended in turbulent flows.     

A hybrid phase field multiple relaxation time lattice Boltzmann method for the incompressible multiphase flow with large density contrast

A hybrid phase field multiple relaxation time lattice Boltzmann method (LBM) is presented in this paper for simulation of multiphase flows with large density contrast. In the present method, the flow field is solved by a lattice Boltzmann equation. Concurrently, the interface of two fluids is captured by solving the macroscopic Cahn‐Hilliard equation using the upwind scheme. To be specific, for simulation of the flow field, an lattice Boltzmann equation (LBE) model developed in Shao et al. (Physical Review E, 89 (2014), 033309) for consideration of density contrast in the momentum equation is used. Moreover, in the present work, the multiple relaxation time collision operator is applied to this LBE to enable simulation of problems with large viscosity contrast or high Reynolds number. For the interface capturing, instead of solving another set of LBE as in many phase field LBMs, the macroscopic Cahn‐Hilliard equation is directly solved by using a weighted essentially non‐oscillatory scheme. In this way, the present hybrid phase field LBM shares full advantages of the phase field LBM while enhancing numerical stability. The ability of the present method to simulate multiphase flow problems with large density contrast is demonstrated by several numerical examples. Copyright © 2015 John Wiley & Sons, Ltd.     

Unstructured lattice Boltzmann method in three dimensions

Over the last decade, the lattice Boltzmann method (LBM) has evolved into a valuable alternative to continuum computational fluid dynamics (CFD) methods for the numerical simulation of several complex fluid‐dynamic problems. Recent advances in lattice Boltzmann research have considerably extended the capability of LBM to handle complex geometries. Among these, a particularly remarkable option is represented by cell‐vertex finite‐volume formulations which permit LBM to operate on fully unstructured grids. The two‐dimensional implementation of unstructured LBM, based on the use of triangular elements, has shown capability of tolerating significant grid distortions without suffering any appreciable numerical viscosity effects, to second‐order in the mesh size. In this work, we present the first three‐dimensional generalization of the unstructured lattice Boltzmann technique (ULBE as unstructured lattice Boltzmann equation), in which geometrical flexibility is achieved by coarse‐graining the lattice Boltzmann equation in differential form, using tetrahedrical grids. This 3D extension is demonstrated for the case of 3D pipe flow and moderate Reynolds numbers flow past a sphere. The results provide evidence that the ULBE has significant potential for the accurate calculation of flows in complex 3D geometries. Copyright © 2005 John Wiley & Sons, Ltd.     

Lattice Boltzmann modeling of shallow water flows over discontinuous beds

Numerical modeling of shallow water flows over discontinuous beds is presented. The flows are described with the shallow water equations and the equations are solved using the lattice Boltzmann method (LBM) with single relaxation time (Bhatnagar–Gross–Krook‐LBM (BGK‐LBM)) and the multiple relaxation time (MRT‐LBM). The weighted centered scheme for force term together with the bed height for a bed slope is described to improve simulation of flows over discontinuous bed. Furthermore, the resistance stress is added to include the local head loss caused by flow over a step. Four test cases, one‐dimensional tidal over regular bed and steps, dam‐break flows, and two‐dimensional shallow water flow over a square block, are considered to verify the present method. Agreements between predictions and analytical solutions are satisfactory. Furthermore, the performance and CPU cost time of BGK‐LBM and MRT‐LBM are compared and studied. The results have shown that the lattice Boltzmann method is simple and accurate for simulating shallow water flows over discontinuous beds. This demonstrates the capability and applicability of the lattice Boltzmann method in modeling shallow water flows on bed topography with a discontinuity in practical hydraulic engineering. Copyright © 2014 John Wiley & Sons, Ltd.     

MEMS稀薄气体内部流动模拟中的信息保存法

首先综述了处理低速稀薄气体流动的一些方法: 线化Boltzmann方程方法、Lattice Boltzmann方法(LBM)、加滑移边界的Navier-Stokes方程、以及DSMC方法, 并讨论它们在模拟MEMS中过渡领域低速流动特别是内部流动所遇到的困难, 其中表明了LBM现有方案不适合模拟过渡领域中的MEMS流动问题. 信息保存(IP)法通过保存一个模拟分子所代表的大量分子的平均信息,克服了流速低使得信息噪声比小而引起统计模拟的困难. 本文给出了方法的一些理论证实. MEMS中内部流动的特点, 即流速低和大的长宽比的特点, 引起椭圆性问题, 即出入口边界条件相互影响需要协调的问题. 通过对(长约几千微米的)微槽道流动应用IP方法的算例,演示了采用守恒形式的质量守恒方程和超松弛法可成功地解决这一问题. 借助同样的方法,用IP方法求解了真实长度(1\,000\,$\mu$m)硬盘驱动器读写头在过渡领域的薄膜支撑问题, 压力分布与具有严格气体动理论基础的概括化Reynolds方程完全相符, 而在此之前, DSMC方法只对短的读写头(5\,$\mu$m)与Reynolds方程做了校验. 作者建议将原来用于求解读写头润滑问题的Reynolds方程退化来求解过渡领域中的微槽道流动问题, 从而提供了一个有严格气体动理论品性的检验方法来验证求解MEMS内部流动的各种方法.      

Airfoil design optimization based on lattice Boltzmann method and adjoint approach

We present a new aerodynamic design method based on the lattice Boltzmann method (LBM) and the adjoint approach. The flow field and the adjoint equation are numerically simulated by the GILBM (generalized form of interpolation supplemented LBM) on non-uniform meshes. The first-order approximation for the equilibrium distribution function on the boundary is proposed to diminish the singularity of boundary conditions. Further, a new treatment of the solid boundary in the LBM is described particularly for the airfoil optimization design problem. For a given objective function, the adjoint equation and its boundary conditions are derived analytically. The feasibility and accuracy of the new approach have been perfectly validated by the design optimization of NACA0012 airfoil.