多孔介质内溶解与沉淀过程的格子Boltzmann方法模拟

本文采用格子Boltzmann方法模拟了多孔介质内的溶解和沉淀现象, 并分析了雷诺数、施密特数、达姆科勒数对多孔介质孔隙结构及浓度分布的影响. 结果表明: 对于多孔介质内的溶解(沉淀)过程, 当雷诺数越大时, 孔隙率越大(小), 平均浓度值越小(大); 当达姆科勒数或施密特数较小时, 溶解和沉淀过程均受反应控制, 此时反应在多孔介质的固体表面较为均匀的发生; 当达姆科勒数或施密特数较大时, 溶解和沉淀过程均受扩散控制, 此时反应主要发生在上游及大孔隙区域.     

用格子玻尔兹曼方法研究流动-反应耦合的非线性渗流问题

根据格子玻尔兹曼计算技术以及相应渗流理论，对多孔介质内流动-反应(矿物介质的溶解等)耦合这一非线性渗流问题进行了数值研究，计算结果与解析解基本符合.数字图像重构技术反映的结果表明流体流动和反应之间可以发生强烈的耦合和反耦合作用，同时可以形成条带结构这一自组织现象，与实验和其他理论分析结果符合也很好. 关键词： 非线性渗流 耦合反应 数值模型     

热对流条件下颗粒沉降的格子Boltzmann方法模拟

采用格子Boltzmann方法模拟了在热对流条件下的颗粒沉降问题, 在研究单颗粒在等温流体、热流体和冷流体中运动的基础上, 进一步模拟了两个不同温度的颗粒在流体中的沉降.结果表明:两等温颗粒的沉降方式与雷诺数Re以及格拉晓夫数Gr密切相关, 而两不同温度的颗粒与两等温颗粒的沉降规律有显著不同.无论初始位置如何, 冷颗粒最终总位于热颗粒下方运动, Re较大时, 发生连续的拖曳、接触现象, 而Re较小时, 冷颗粒会以较大的沉降速度远离热颗粒. 关键词： 格子Boltzmann方法 颗粒沉降 热对流     

Lattice Boltzmann simulation of 3-dimensional natural convection heat transfer of CuO/water nanofluids

The present study investigated fluid flow and natural convection heat transfer in an enclosure embedded with isothermal cylinder. The purpose was to simulate the three-dimensional natural convection by thermal lattice Boltzmann method based on the D3Q19 model. The effects of suspended nanoparticles on the fluid flow and heat transfer analysis have been investigated for different parameters such as particle volume fraction, particle diameters, and geometry aspect ratio. It is seen that flow behaviors and the average rate of heat transfer in terms of the Nusselt number (Nu) are effectively changed with different controlling parameters such as particle volume fraction (5 % ≤ φ ≤ 10 %), particle diameter (d _p = 10 nm to 30 nm) and aspect ratio (0.5 ≤ AR ≤ 2) with fixed Rayleigh number, Ra = 10⁵. The present results give a good approximation for choosing an effective parameter to design a thermal system.     

Lattice Boltzmann simulation of solid particles suspended in fluid

The lattice Boltzmann method, an alternative approach to solving a fluid flow system, is used to analyze the dynamics of particles suspended in fluid. The interaction rule between the fluid and the suspended particles is developed for real suspensions where the particle boundaries are treated as no-slip impermeable surfaces. This method correctly and accurately determines the dynamics of single particles and multi-particles suspended in the fluid. With this method, computational time scales linearly with the number of suspensions,N, a significant advantage over other computational techniques which solve the continuum mechanics equations, where the computational time scales asN ³. Also, this method solves the full momentum equations, including the inertia terms, and therefore is not limited to low particle Reynolds number.     

Lattice Boltzmann simulation for temperature-sensitive magnetic fluids in a porous square cavity

A lattice Boltzmann method is developed to simulate temperature-sensitive magnetic fluids in a porous cavity. In the simulation, the magnetic force, efficient gravity, viscous loss term and geometric loss term in porous medium are imported to the momentum equation. To test the reliability of the method, a validation with water in porous cavity is carried out. Good agreements with the previous results verify that the present lattice Boltzmann method is promising for simulation of magnetic fluids in porous medium. In this study, we investigate the change of magnetization with external magnetic field, and we present numerical results for the streamlines, isotherms, and magnetization at vertical or horizontal mid-profiles for different values of Ram. In addition, Nusselt numbers changing with magnetic Rayleigh numbers are also investigated.     

Lattice Boltzmann simulation of fluid flows in two-dimensional channel with complex geometries

Boundary conditions (BCs) play an essential role in lattice Boltzmann (LB) simulations. This paper investigates several most commonly applied BCs by evaluating the relative L₂-norm errors of the LB simulations for two-dimensional (2-D) Poiseuille flow. It is found that the relative L₂-norm error resulting from FHML's BC is smaller than that from other BCs as a whole. Then, based on the FHML's BC, it formulates an LB model for simulating fluid flows in 2-D channel with complex geometries. Afterwards, the flows between two inclined plates, in a pulmonary blood vessel and in a blood vessel with local expansion region, are simulated. The numerical results are in good agreement with the analytical predictions and clearly show that the model is effective. It is expected that the model can be extended to simulate some real biologic flows, such as blood flows in arteries, vessels with stenosises, aneurysms and bifurcations, etc.     

多孔介质中盐指现象的数值模拟

运用基于杂交网格的高精度数值方法研究了多孔介质中的盐指现象.该算法将基于边界拟合坐标下的高精度有限差分法和高精度的泊松方程快速求解器有效地结合在一起,从而达到提高整体的计算精度、计算效率和稳定性的目的.通过比较不同孔隙率的多孔介质对盐指对流的传热传质效应的影响,发现在标准孔隙率较低的多孔介质中,盐度扩散的速度明显比热扩散的速度快,盐指很快触及上下壁面,使得上下层的盐度梯度迅速减小,这是与非多孔介质具有明显差异之处. 关键词： 多孔介质 双扩散对流 盐指     

耦合界面张力的三维流体界面不稳定性的格子Boltzmann模拟

采用介观格子Boltzmann方法模拟界面张力作用下三维流体界面的Rayleigh-Taylor (RT)不稳定性的增长过程,主要分析表面张力对流体界面动力学行为及尖钉和气泡后期增长的影响机制.首先发现三维RT不稳定性的发生存在临界表面张力(σ_c),其值随着流体Atwood数的增大而增大,且数值预测值与理论分析结果σ_c=(ρ_h-ρ₁)g/k~2一致.另外,随着表面张力的增大,不稳定性演化过程中界面卷吸程度和结构复杂性逐渐减弱,系统中界面破裂形成离散液滴的数目也显著减少.相界面的后期动力学行为也从非对称发展转向始终保持关于中轴线对称.尖钉与气泡振幅在表面张力较小时对其变化不显著,当表面张力增大到一定值后,可以有效地抑制尖钉与气泡振幅的增长.进一步发现,高雷诺数三维RT不稳定性在不同表面张力下均经历4个不同的发展阶段:线性阶段、饱和速度阶段、重加速和混沌混合阶段.尖钉与气泡在饱和速度阶段以近似恒定的速度增长,其渐进速度的值与修正的势流理论模型结果一致.受非线性Kelvin-Helmholtz旋涡的剪切作...     

Lattice Boltzmann simulation of transverse wave travelling in Maxwell viscoelastic fluid

A nine-velocity lattice Boltzmann method for Maxwell viscoelastic fluid is proposed. Travelling of transverse wavein Maxwell viscoelastic fluid is simulated. The instantaneous oscillating velocity, transverse shear speed and decay rateagree with theoretical results very well.     

Lattice Boltzmann method on quadtree grids for simulating fluid flow through porous media: A new automatic algorithm

During the past two decades, the lattice Boltzmann (LB) method has been introduced as a class of computational fluid dynamic methods for fluid flow simulations. In this method, instead of solving the Navier Stocks equation, the Boltzmann equation is solved to simulate the flow of a fluid. This method was originally developed based on uniform grids. However, in order to model complex geometries such as porous media, it can be very slow in comparison with other techniques such as finite differences and finite elements. To eliminate this limitation, a number of studies have aimed to formulate the lattice Boltzmann on the unstructured grids. This paper deals with simulating fluid flow through a synthetic porous medium using the LB method and on the quadtree grid structure. To this end, the LB method was used on nonuniform grids coupled with a technique for image reconstruction which resulted in the quadtree grids for simulation of fluid flow through porous media. Accuracy and efficiency of this algorithm is compared against the conventional LB method based on uniform grids. While the decrease in computational time in the proposed LB method on nonuniform grids is found to be significant regarding the size of the initial and reconstructed images, the same level of accuracy is obtained when compared with the conventional LB method on uniform grids.     

Lattice Boltzmann simulation for hydrothermal analysis of free convection within dumbbell-shaped heat exchanger

The lattice Boltzmann simulation of nanofluid flow and heat transfer during natural convection within a dumbbell-shaped heat exchanger is carried out. The heat exchanger is filled with CuO–water. The KKL model is employed to predict the thermo-physical properties of nanofluid. In order to perform a comprehensive hydrothermal investigation, different post-processing approaches such as heatline visualization, total entropy generation, local entropy generation based on local fluid friction irreversibility and heat transfer irreversibility, average and local Nusselt variation are employed. In the present investigation, it is tried to present the impact of different influential parameters like Rayleigh number, solid volume fraction of nanofluid and thermal arrangement of internal fins-bodies on the fluid flow, heat transfer rate and entropy generation.     

Magnetic field stabilization of a two-dimensional fluid jet: a multiple relaxation Lattice Boltzmann simulation

The stabilization of a two-dimensional fluid jet from the Kelvin–Helmholtz (KH) instability by an external parallel magnetic field is examined by lattice Boltzmann techniques. For sufficiently strong magnetic fields, the jet does not break up into large-scale vortices but retains the major features of the jet, albeit somewhat expanded. There are time-dependent striations within the expanded jet.     

用晶格玻尔兹曼方法研究微结构表面的疏水性能

将固体表面分别近似为具有简单的周期性矩形、三角形和半圆形微粗糙结构表面,建立了两相流的晶格玻尔兹曼模型.通过测量不同微粗糙结构表面上液滴的接触角,探讨微结构形状和尺寸的改变对固体材料表面疏水性能的影响.最后,由流体在各种糙壁管中的速度滑移,验证了结论的正确性.     

Radiative effects on magnetohydrodynamic natural convection flows saturated in porous media

A boundary layer solution is presented to study the effects of joule heating on magnetohydrodynamic natural convection flow. The Rosseland approximation is used to describe the radiative heat flux in the energy equation. Four different cases of flows have been studied namely an isothermal surface, a uniform heat flux surface, a plane plume and flow generated from—a horizontal line energy source a vertical adiabatic surface. Numerical results presented for the perturbation analysis four boundary conditions with various parameters are tabulated.     

用晶格玻尔兹曼方法研究颗粒在涡流中的运动

用晶格玻尔兹曼方法研究小颗粒在涡流中的运动.涡流由流经空腔的流体产生.用动量交换法和压力张量积分法计算颗粒在涡流中的运动轨迹、速度和角速度.最后用张量积分法计算两个不同半径的颗粒在涡流中的运动. 关键词： 晶格玻尔兹曼方法 涡流 颗粒     

Lattice Boltzmann simulation of a binary fluid with different phase viscosities and its application to fingering in two dimensions

A thermodynamically consistent lattice Boltzmann scheme for simulating the flow of a binary fluid is extended to allow the fluid components to have different viscosities. The approach is tested for the shear and Poiseuille flow of layered immiscible fluids and for the dispersion relation and the damping of a capillary wave. We then consider the fingering that results when a fluid is displaced by a less viscous fluid in a two-dimensional channel. The finger widths obtained match the results of Reinelt and Saffman [#!Reinelt85!#], but differ somewhat from those of Halpern and Gaver [#!Halpern94!#] for capillary numbers above 2. A limiting finger width close to 1/2 is obtained for high capillary numbers and high viscosity ratios. Received 25 May 1999 and Received in final form 19 November 1999     

Theoretical and numerical study of a thermal convection problem with temperature-dependent viscosity in an infinite layer

A convection problem with temperature-dependent viscosity in an infinite layer is presented. This problem has important applications in mantle convection. The existence of a stationary bifurcation is proven together with a condition to obtain the critical parameters at which the bifurcation takes place. A numerical strategy has been developed to calculate the critical bifurcation curves and the most unstable modes for a general dependence of viscosity on temperature. An exponential dependence of viscosity on temperature has been considered in the numerical calculations. Comparisons with the classic Rayleigh-Bénard problem with constant viscosity indicate that the critical temperature difference threshold decreases as the exponential rate parameter increases. The vertical velocity of the marginal mode exhibits motion concentrated in the region where viscosity is smaller.