摘要:
针对滑移区复杂气-固边界存在速度滑移现象,提出了一种基于格子Boltzmann方法的非平衡态外推与有限差分相结合的曲边界处理新格式.该格式具有可考虑实际物理边界与网格线偏移量的优势,较传统half-way DBB (diffusive bounce-back)格式更能准确反映实际边界情况,同时还可获取壁面处气体宏观量及其法向梯度等信息.采用本文所提曲边界处理格式模拟分析了滑移区气体平直/倾斜微通道Poiseuille流、微圆柱绕流和同心微圆柱面旋转Couette流问题.研究结果表明,采用曲边界处理新格式所得结果与理论值以及文献结果符合良好,适用于滑移区气体流动的复杂边界处理,且比half-way DBB格式具有更高的精度,较修正DBB格式具有更好的适应性.
Abstract:
A new curved boundary treatment in lattice Boltzmann method is developed for micro gas flow in the slip regime.The proposed treatment is a combination of the nonequilibrium extrapolation scheme for curved boundary with noslip velocity condition and the counter-extrapolation method for the velocity and its normal gradient on the curved boundary.Taking into consideration the effect of the offset between the physical boundary and the closest grid line,the new treatment is proved to be more accurate than the traditional half-way diffusive bounce-back (DBB) scheme.The present treatment is also more applicable than the modified DBB scheme because the specific gas-wall interaction parameters need to be determined to ensure the validation of the modified DBB scheme.The proposed boundary treatment is implemented to simulate the benchmark problems,which include a Poiseuille flow in the aligned/inclined micro-channel,a flow past a inicrocylinder and a microcylindrical Couette flow.The results and conclusions are summarized as follows.l) The force-driven Poiseuille flow in an aligned microchannel is simulated separately with different values of wallgrid offset qδx (q =0.25,0.5,0.75,1.0).With the consideration of the wall-grid offset,the numerical results with the new boundary treatment show good agreement with the analytical solut7ons.However,the results obtained by using the half-way DBB scheme only accord well with the analytical solutions under the condition of a fixed wall-grid offset (q =0.5).2) To demonstrate the capability of the present treatment in dealing with gas flow in a more complex geometry,the force-driven Poiseuille flow in a micro-channel is investigated separately with different inclined angles.The present numerical results fit well with the analytical solutions.However,the discrepancy between the results obtained with the half-way DBB scheme and the analytical solutions can be clearly observed near the inclined boundaries.3) The gas flow past a microcylinder is simulated to prove that the present treatment can deal with the curved boundary.The slip velocity profile along the micro cylinder periphery obtained with the present treatment accords well with the available data in the published literature.However,the results obtained with the half-way DBB scheme show lower values than the results from the published work.4) In the simulations of the microcylindrical Couette flow between two coaxial rotating cylinders for different Knudsen numbers the results obtained by using the present treatment show excellent agreement with the analytical solutions.However,the results obtained with the half-way DBB scheme and the modified DBB scheme deviate obviously from the analytical solutions near the inner and outer cylindrical walls,respectively.In summary,the new boundary treatment proposed in this work is capable of dealing with the complex gas-solid boundary in the slip regime.The new treatment has a higher accuracy than the half-way DBB scheme and shows a better applicability than the modified DBB scheme.