一维非定常流的格子Boltzmann模拟

本文建立了一维非定常流的一般格子Ｂｏｌｔｚｍａｎｎ模型，并利用２－速度模型研究了激波管和激波的形成，反射以及相互作用。给了了不同的初始条件和边界条件的模拟结果，它与理论分析放数据分析的结果十分吻合。     

一维Tyson反应扩散系统的格子Boltzmann方法模拟

建立基于格子Boltzmann模型的一维Tyson反应扩散系统的数值求解法.利用浓度分布的Chapmann-Enskogz展开及多尺度技术,获得激励介质在反应与扩散机制同时作用的一维反应扩散方程,用于求解Tyson反应扩散的反应物和催化剂随时间的浓度空间分布值.数值结果表明本文中所提供的方法是有效的.     

一类偏微分方程的格子Boltzmann模型

研究了对流扩散方程、Burgers方程和Modified-Burgers方程等具有相同形式的一类偏微分方程。并且构建了带修正函数项的D1Q3格子Boltzmann模型求解这类方程。为了能准确地恢复出此宏观方程,利用Chapman-Enskog展开和多尺度分析技术,推导出了各个方向的平衡态分布函数和修正函数的具体表达式。数值计算结果表明该模型是稳定、有效的。     

KdV方程带修正函数的格子Boltzmann模拟

采用一种带修正函数的新格子Boltzmann模型模拟了KdV方程,分析了由此得出的迭代格式的单调性和稳定性,得到了格式的单调性条件。在单调性条件下,迭代格式是[L1]稳定的。数值模拟结果表明该格式是可行的。     

求解二维对流扩散方程的格子Boltzmann方法

针对二维对流扩散方程,基于D2Q4格子速度,用Chapman-Enskog多尺度分析技术,将时间尺度取为二阶,空间尺度取为一阶,推导了各个速度方向上的平衡态分布函数所满足的条件,给出了简单且对称的平衡态分布函数表达式,所得到的平衡态分布函数能正确地恢复出二维对流扩散方程,从而构建了一种新的求解二维对流扩散方程的D2Q4格子Boltzmann（LB）模型。用所给LB模型对扩散方程和两个不同初边界条件的对流扩散方程进行了数值求解,数值实验结果表明数值解与精确解吻合较好,与相关文献结果比较边界误差要小得多,验证了模型的有效性。     

利用多尺度分析方法求解KdV方程和KdV—Burgers方程的格子Boltzmann模型

本文建立了用格点法解一般偏微分方程(PDE)的理论框架,构造出求解KdV方程及KdV—Burgers方程的三速格子BGK模型。引进三种时间尺度,利用多尺分析求出Boltzmann演化方程的平衡分布函数。     

后台阶流动的非均匀格子Boltzmann方法模拟

使用非均匀格子Boltzmann方法对后台阶流动进行了数值模拟.将流体流动区域划分为不同的子区域:对于每个子区域内部,分布函数使用均匀网格计算;对于区域边界,分布函数采用嵌套网格方法进行处理.数值计算结果与其它实验、数值结果相吻合.     

基于格子Boltzmann方法的换热器优化模拟

为优化换热器的结构设计,用格子Boltzmann方法(Lattice Boltzmann Method,LBM)结合多孔介质模型模拟换热器内的换热,研究雷诺数、普朗特数和热扩散率比的变化对温度场和换热性能的影响.模拟结果表明:在小雷诺数范围内,随着雷诺数的增加,努塞尔数先增加后减小,即存在一个使换热性能达到最好的雷诺数;随着普朗特数的增加,努塞尔数减小,换热性能降低;随着热扩散率比的增加,换热性能提高.分析不同管柱排列方式对换热性能的影响,结果表明:叉排的换热效果明显优于顺排,当横向节距等于2时,对于均匀顺排或叉排,努塞尔数均随纵向节距的增加而减小,这与实验结果相符;对于非均匀叉排,采用\"前密\"或\"中间密\"的排布方式有利于换热.     

二维带分数阶Laplacian算子的对流扩散方程的格子Boltzmann模型研究

带有分数阶Laplacian算子的对流扩散方程常被用来刻画自然界与社会系统中的反常扩散现象.本文提出了一种新的格子Boltzmann模型,用于求解二维带分数阶Laplacian算子的对流扩散方程.首先,基于分数阶Laplacian算子的Fourier变换和Gauss型求积公式,得到控制方程的近似方程.然后,将速度空间、时间和空间进行离散,并构造合适的平衡态分布函数和离散作用力,建立有效的格子Boltzmann-BGK模型.通过Chapman-Enskog分析,可由建立的格子Boltzmann-BGK模型恢复出宏观方程,从而证明了模型的有效性.最后,将模型应用于求解带有解析解的数值算例和Allen-Cahn方程,数值结果进一步验证了模型的正确性和有效性.     

多重网格格子Boltzmann方法的并行算法

针对复杂流动数值模拟中的格子Boltzmann方法存在计算网格量大、收敛速度慢的缺点,提出了基于三维几何边界的多重笛卡儿网格并行生成算法,并基于该网格生成方法提出了多重网格并行格子Boltzmann方法（LBM）。该方法结合不同尺度网格间的耦合计算,有效减少了计算网格量,提高了收敛速度;而且测试结果也表明该并行算法具有良好的可扩展性。     

Sailfish: A flexible multi-GPU implementation of the lattice Boltzmann method

We present Sailfish, an open source fluid simulation package implementing the lattice Boltzmann method (LBM) on modern Graphics Processing Units (GPUs) using CUDA/OpenCL. We take a novel approach to GPU code implementation and use run-time code generation techniques and a high level programming language (Python) to achieve state of the art performance, while allowing easy experimentation with different LBM models and tuning for various types of hardware. We discuss the general design principles of the code, scaling to multiple GPUs in a distributed environment, as well as the GPU implementation and optimization of many different LBM models, both single component (BGK, MRT, ELBM) and multicomponent (Shan–Chen, free energy). The paper also presents results of performance benchmarks spanning the last three NVIDIA GPU generations (Tesla, Fermi, Kepler), which we hope will be useful for researchers working with this type of hardware and similar codes.     

Large eddy simulation of turbulent open duct flow using a lattice Boltzmann approach

Large eddy simulations of turbulent open duct flow are performed using the lattice Boltzmann method (LBM) in conjunction with the Smagorinsky sub-grid scale (SGS) model. A smaller value of the Smagorinsky constant than the usually used one in plain channel flow simulations is used. Results for the mean flow and turbulent fluctuations are compared to experimental data obtained in an open duct of similar dimensions. It is found that the LBM simulation results are in good qualitative agreement with the experiments.     

Immersed boundary method and lattice Boltzmann method coupled FSI simulation of mitral leaflet flow

Coupling the immersed boundary (IB) method and the lattice Boltzmann (LB) method might be a promising approach to simulate fluid-structure interaction (FSI) problems with flexible structures and moving boundaries. To investigate the possibility for future IB-LB coupled simulations of the heart flow dynamics, an IB-LB coupling scheme suitable for rapid boundary motion and large pressure gradient FSI is proposed, and the mitral valve jet flow considering the interaction of leaflets and fluid is simulated. After analyzing the respective concepts, formulae and advantages of the IB and LB methods, we first explain the coupling strategy and detailed implementation procedures, and then verify the effectiveness and second-order accuracy of the scheme by simulating a benchmark case, the relaxation of a stretched membrane immersed in fluid. After that, the diastolic filling jet flow between mitral leaflets in a simplified 2D left heart model is simulated. The model consists of the simplified transmitral passage of the heart and two curvilinear leaflets. In the simulation, the atrial and ventricular pressure histories of normal human are specified as boundary conditions, and the leaflets are treated as fibers that interact with the fluid to define their deformations and movements. The resulting opening and closing movements of the leaflets and the flow patterns of the filling jet are qualitatively reasonable and compare well with existing numerical and measured data. It is shown that this IB-LB coupling method is feasible for treating flexible boundary FSI problems with rapid boundary motion and large pressure gradient, the results of the mitral leaflet flow are valuable for understanding the transmitral FSI dynamics, and it is possible to simulate the more realistic 3D heart flow by the scheme in the future.     

Application of lattice Boltzmann method in free surface flow simulation of micro injection molding

The filling flow in micro injection molding was simulated by using the lattice Boltzmann method (LBM). A tracking algorithm for free surface to handle the complex interaction between gas and liquid phases in LBM was used for the free surface advancement. The temperature field in the filling flow is also analyzed by combining the thermal lattice Boltzmann model and the free surface method. To simulate the fluid flow of polymer melt with a high Prandtl number and high viscosity, a modified lattice Boltzmann scheme was adopted by introducing a free parameter in the thermal diffusion equation to overcome the restriction of the thermal relaxation time. The filling flow simulation of micro injection molding was successfully performed in the study.     

基于下三角H矩阵的LDLC整形研究

针对低密度格码(LDLC)编码后不能适应在功率限制的加性高斯白噪声(AWGN)信道上传输的问题,研究了编码前的整形问题。构造了一种特殊结构的下三角H矩阵,并结合超立方整形和系统整形方案,分析并仿真了整形前后格点的变化及取得的整形增益。结果显示：整形后的码字沿格的Voronoi域均匀分布,并当码长为10000、误码率为10 -5时,可得到1.31dB的整形增益,相对于传统的整形技术提高了0.31dB。通过整形,有效地产生了功率受限的格点。     

Lattice Boltzmann based PDE solver on the GPU

In this paper, we propose a hardware-accelerated PDE (partial differential equation) solver based on the lattice Boltzmann model (LBM). The LBM is initially designed to solve fluid dynamics by constructing simplified microscopic kinetic models. As an explicit numerical scheme with only local operations, it has the advantage of being easy to implement and especially suitable for graphics hardware (GPU) acceleration. Beyond the Navier–Stokes equation of fluid mechanics, a typical LBM can be modified to solve the parabolic diffusion equation, which is further used to solve the elliptic Laplace and Poisson equations with a diffusion process. These PDEs are widely used in modeling and manipulating images, surfaces and volumetric data sets. Therefore, the LBM scheme can be used as an GPU-based numerical solver to provide a fast and convenient alternative to traditional implicit iterative solvers. We apply this method to several examples in volume smoothing, surface fairing and image editing, achieving outstanding performance on contemporary graphics hardware. It has the great potential to be used as a general GPU computing framework for efficiently solving PDEs in image processing, computer graphics and visualization.     

A front-tracking lattice Boltzmann method to study flow-induced deformation of three-dimensional capsules

In this paper, a hybrid method is proposed to study the flow-induced deformation of three-dimensional capsules. The capsules consist of Newtonian liquid drops enclosed by thin elastic membranes. In the proposed approach, the front-tracking method is coupled with the lattice Boltzmann method. The fluids inside and outside the capsule is treated as one fluid with varying physical properties, and is modeled by the lattice Boltzmann equation. The capsule membrane is explicitly tracked by the membrane nodes that are advected by the flow. The multi-block strategy of the lattice Boltzmann method is employed to refine the mesh near the capsule, which greatly increase the accuracy and efficiency of the three-dimensional computation. The capsule membrane is discretized into unstructured flat triangular elements, and a finite element model is incorporated to account for the membrane mechanics. With the present method, the transient deformation of initially spherical capsules with membrane following Neo-Hookean constitutive laws is simulated in shear flow, under various dimensionless shear rates and ratios of internal to surrounding liquid viscosities. The present results, including the Taylor shape parameter, the capsule inclination angle and the tank-treading frequency, agree well with previously published numerical results.     

Analysis of single droplet dynamics on striped surface domains using a lattice Boltzmann method

In this paper, the behavior of a micron-scale fluid droplet on a heterogeneous surface is investigated using a two-phase lattice Boltzmann method (LBM). The two-phase LBM permits the simulation of the time dependent three-dimensional motion of a liquid droplet on solid surface patterned with hydrophobic and hydrophilic strips. A nearest-neighbor molecular interaction force is used to model the adhesive forces between the fluid and solid walls. The solid heterogeneous wall is a uniform hydrophilic substrate painted with hydrophobic strips. The model is validated by demonstrating the consistency of the simulation results with an exact solution for capillary rise and through qualitative comparison of computed dynamic contact line behavior with experimentally measured surface properties and observed surface shapes of a droplet on a heterogeneous surface. The dependence of spreading behavior on wettability, the width of hydrophobic strip, initial location of the droplet relative to the strips, and gravity is investigated. A decrease in contact angle of the liquid on a hydrophilic surface may lead to breakup of the droplet for certain substrate patterns. The simulations suggest that the present lattice Boltzmann (LB) model can be used as a reliable way to study fluidic control on heterogeneous surfaces and other wetting related subjects.