基于格子玻尔兹曼法的三维溃坝流数值模拟与验证
|
摘要
溃坝问题一直是验证数值模拟方法有效性的经典问题。利用基于格子玻尔兹曼离散方法研发的XFlow软件来模拟三维溃坝问题,以此验证其有效性并提出进一步的改进方向。为了对比前人所做的实验和模拟结果,计算了无障碍三维溃坝流相应监测点的压力、水位等要素,得到了水流演进过程,并提出了用增强固壁函数法处理固壁边界以增强计算的稳定性,另外还考虑了有矩形障碍物的三维溃坝问题。研究结果表明:该方法能捕捉到某一时间点的自由液面,相关数据计算结果准确,且无需划分网格,能减少工作量,但结果稳定性较差,合理运用边界处理方法能有效解决这一问题。
Dam-break flood simulation is a classical issue to verify the effectiveness of numerical simulation. In this paper,3 D dam-break flood is simulated and studied by XFlow software based on Lattice Boltzmann Method to verify its validity and find out improvement approach. In order to compare the results of previous experiments and simulations,the pressure,water level etc. at the corresponding monitoring points for barrier-free 3 D dam-break flow are calculated and the evolution process of the flooding is obtained. In the calculation,solid boundary is treated by Enhanced Solid wall to increase calculation stability. In addition,3 D dam-break flood discharge with rectangular obstacle is also considered. The research results show that the proposed method can capture the free liquid surface at a certain time point,the relevant calculation results are accurate and the work load is reduced because of meshing needless. However,the stability of the results is poor and the problem can be effectively solved by reasonably using the boundary treatment method.
Dam-break flood simulation is a classical issue to verify the effectiveness of numerical simulation. In this paper,3 D dam-break flood is simulated and studied by XFlow software based on Lattice Boltzmann Method to verify its validity and find out improvement approach. In order to compare the results of previous experiments and simulations,the pressure,water level etc. at the corresponding monitoring points for barrier-free 3 D dam-break flow are calculated and the evolution process of the flooding is obtained. In the calculation,solid boundary is treated by Enhanced Solid wall to increase calculation stability. In addition,3 D dam-break flood discharge with rectangular obstacle is also considered. The research results show that the proposed method can capture the free liquid surface at a certain time point,the relevant calculation results are accurate and the work load is reduced because of meshing needless. However,the stability of the results is poor and the problem can be effectively solved by reasonably using the boundary treatment method.
引文
[1]Colagrossi A,Landrini M.Numerical simulation of interfacial flows by smoothed particle hydrodynamics[J].Journal of Computational Physics,2003,2(191):448-475.
[2]曹志先,郑力,钱忠东.溃坝水流三维湍流的试验与数值分析[J].武汉大学学报:工学版,2010,43(2):137-142.
[3]张雨新,万德成.MPS方法在三维溃坝问题中的应用[J].中国科学:物理学力学天文学,2011(2):140-154.
[4]Cherfils J M,Pinon G,Rivoalen E.JOSEPHINE:A parallel SPH code for free-surface flows[J].Computer Physics Communications,2012,183(7):1468-1480.
[5]曹洪建,万德成,杨驰.三维溃坝波绕方柱剧烈流动的数值模拟[J].水动力学研究与进展(A辑),2013,28(4):414-422.
[6]许晓阳,任胜章,邓方安.三维溃坝流的光滑粒子动力学方法模拟[J].计算力学学报,2016,55(5):676-682.
[7]何雅玲,王勇,李庆.格子Boltzmann方法的理论及应用[M].北京:科学出版社,2009.
[8]Zhou Z Q,Kat J,Buchner B.A Nonlinear 3-D approach to simulate green water dynamics on deck[A].Proceeding of the 7th International Conference on Numerical Ship Hydrodynamics[C]∥France:Nantes,1999.
[9]Mcnamara G R,Zanetti G.Use of the Boltzmann equation to simulate lattice automata[J].Physical Review Letters,1988,61(20):2332-2335.
[10]Bhatnagar P L,Gross E P,Krook M.A model for collision processes in gases.I.Small amplitude processes in charged and neutral onecomponent systems[J].Physical Review,1954(94):511-525.
[11]李慎龙,刘树成,邢庆坤,等.基于LBM-LES模拟的离合器摩擦副流致运动效应[J].吉林大学学报:工学版.2017,47(2):490-497.
[12]Qian H Y,D Humières D,Lallemand P.Lattice BGK models for Navier-Stokes equation[J].Europhysics Letters,1992,6(17):479.
[13]罗金辉.煤层气运移LBM模型与井间干扰模拟研究[D].徐州:中国矿业大学,2012.
[14]钱忠东,黄社华.四种湍流模型对空化流动模拟的比较[J].水科学进展,2006(2):203-208.
[15]赵雪莹,王昭升,盛金保.梯级水库溃坝洪水模拟[J].人民长江,2017,48(11):32-35.
[16]葛杭建,邓波,赵容娇.亭下水库溃坝洪水分析[J].人民长江,2017,48(S1):8-13.
[2]曹志先,郑力,钱忠东.溃坝水流三维湍流的试验与数值分析[J].武汉大学学报:工学版,2010,43(2):137-142.
[3]张雨新,万德成.MPS方法在三维溃坝问题中的应用[J].中国科学:物理学力学天文学,2011(2):140-154.
[4]Cherfils J M,Pinon G,Rivoalen E.JOSEPHINE:A parallel SPH code for free-surface flows[J].Computer Physics Communications,2012,183(7):1468-1480.
[5]曹洪建,万德成,杨驰.三维溃坝波绕方柱剧烈流动的数值模拟[J].水动力学研究与进展(A辑),2013,28(4):414-422.
[6]许晓阳,任胜章,邓方安.三维溃坝流的光滑粒子动力学方法模拟[J].计算力学学报,2016,55(5):676-682.
[7]何雅玲,王勇,李庆.格子Boltzmann方法的理论及应用[M].北京:科学出版社,2009.
[8]Zhou Z Q,Kat J,Buchner B.A Nonlinear 3-D approach to simulate green water dynamics on deck[A].Proceeding of the 7th International Conference on Numerical Ship Hydrodynamics[C]∥France:Nantes,1999.
[9]Mcnamara G R,Zanetti G.Use of the Boltzmann equation to simulate lattice automata[J].Physical Review Letters,1988,61(20):2332-2335.
[10]Bhatnagar P L,Gross E P,Krook M.A model for collision processes in gases.I.Small amplitude processes in charged and neutral onecomponent systems[J].Physical Review,1954(94):511-525.
[11]李慎龙,刘树成,邢庆坤,等.基于LBM-LES模拟的离合器摩擦副流致运动效应[J].吉林大学学报:工学版.2017,47(2):490-497.
[12]Qian H Y,D Humières D,Lallemand P.Lattice BGK models for Navier-Stokes equation[J].Europhysics Letters,1992,6(17):479.
[13]罗金辉.煤层气运移LBM模型与井间干扰模拟研究[D].徐州:中国矿业大学,2012.
[14]钱忠东,黄社华.四种湍流模型对空化流动模拟的比较[J].水科学进展,2006(2):203-208.
[15]赵雪莹,王昭升,盛金保.梯级水库溃坝洪水模拟[J].人民长江,2017,48(11):32-35.
[16]葛杭建,邓波,赵容娇.亭下水库溃坝洪水分析[J].人民长江,2017,48(S1):8-13.