利用多松弛格子Boltzmann方法预测多孔介质的渗透率

多孔介质内的流动问题在工程热物理领域有着重要的研究价值和应用背景。本文利用多松弛格子Boltzmann方法详细预测了两种二维多孔介质中的渗透率。研究结果表明:一方面,多松弛模型可以用来克服由单松弛模型带来的一些不足;另一方面,借助于达西定律,多松弛模型可以准确预测多孔介质的渗透率,并将计算结果与已有文献做了对比。     

非稳态横掠管束周期性充分发展对流换热的数值模拟

本文采用非稳态数学模型对横掠管束的周期性充分发展流动与换热进行了数值模拟。结果显示,在所考虑的参数条件下,Re≥100以后流场和温度场随时间发生振荡,并且振幅随Re增加而增大;进一步增加Re到一定值后,在管束间会有交替出现的涡。计算得到的平均Nusselt数与现有文献中的实验关联式基本是吻合的。平均Nusselt数的计算值随流动方向管间距的增加而增加。     

空调蒸发器用管翅式换热器开缝翅片的数值设计

通过三维层流数值模拟,与平片计算结果对比,设计了适用于低流速和Re数(迎面风速1-3 m／s,对应Re数 901-2702)下的空调蒸发器用管翅式换热器的开缝翅片形式。研究表明:在上述迎面风速和Re数范围内,开缝翅片性能曲线和平片性能曲线不可避免地有一个交叉点,对应Re数称为转折Re数,在交叉点前的速度和Re数范围内,平片换热器的综合性能优于开缝翅片,而在高于转折Re数后,随着迎面风速的增大,开缝翅片的综合换热性能将越来越优于平片;依据“前疏后密”原则适当减少开缝翅片换热器开缝的条数,可以有效降低转折Re数,并显著改善低流速和Re数下换热器的综合性能。另外,再次验证了场协同理论-温度场和速度场的协同性与换热器换热量间的必然联系。     

黏性液体中单个气泡上升的形状特性

采用基于Level Set方法的直接数值模拟技术对黏性液体中单个气泡的上升运动进行三维模拟.数值模拟采用拟单相流模型处理气泡内外的气液两相流动,应用Level Set方法捕捉运动气泡的变形.针对Eo数从O(0)～O(2),Mo数从O(-11)～O(2)的流动范围,重点研究了上升气泡的形状特性,并与经典的气泡形状图谱进行了比较.模拟结果表明,上升气泡的形状与无量纲参数(Eo、Mo和Re)密切相关.在高Re的扁椭球区域,数值发现了气泡形状的周期性振荡行为.     

Multi-relaxation-time lattice Boltzmann front tracking method for two-phase flow with surface tension

<正>In this paper,an improved incompressible multi-relaxation-time lattice Boltzmann-front tracking approach is proposed to simulate two-phase flow with a sharp interface,where the surface tension is implemented.The lattice Boltzmann method is used to simulate the incompressible flow with a stationary Eulerian grid,an additional moving Lagrangian grid is adopted to track explicitly the motion of the interface,and an indicator function is introduced to update the fluid properties accurately.The interface is represented by using a four-order Lagrange polynomial through fitting a set of discrete marker points,and then the surface tension is directly computed by using the normal vector and curvature of the interface.Two benchmark problems,including Laplace’s law for a stationary bubble and the dispersion relation of the capillary wave between two fluids are conducted for validation.Excellent agreement is obtained between the numerical simulations and the theoretical results in the two cases.     

声波作用下球形颗粒外声流分布的数值模拟

综合考虑声学边界层内的热损失和黏性损失,建立处于平面驻波声压波节位置二维球形颗粒外声流计算模型,利用分离时间尺度的数值方法对颗粒外声流流场特征进行模拟.将模拟结果与相应的解析解和实验结果对比,验证了数值模拟的可靠性.在此基础上,研究了雷诺数Re和斯特劳哈尔数Sr对球形颗粒声学边界层内二阶声流流场结构、涡流强度及范围的影...     

Thermodynamic of collapsing cavitation bubble investigated by pseudopotential and thermal MRT-LBM

The thermodynamic of cavitation bubble collapsing is a complex fundamental issue for cavitation application and prevention. The pseudopotential and thermal multi-relaxation-time lattice Boltzmann method (MRT-LBM) is adopted to investigate the thermodynamic of collapsing cavitation bubble in this paper. The simulation results satisfy the maximum temperature equation of the bubble collapse, which derived from the Rayleigh-Plesset (R-P) equation. The validity of thermal MRT-LBM in simulating the collapse process of cavitation bubble is verified. It shows that the temperature evolution of vapor-liquid phase is well captured. Furthermore, the two-dimensional (2D) temperature, velocity and pressure field of the bubble near a solid wall are analyzed. The maximum temperature inside the bubble and wall temperature under different position offset parameters are discussed in details.     

A novel method of evaluating the lift force on the bluff body based on Noca's flux equation

The influence of experimental error on lift force evaluated by Noca's flux equation is studied based on adding errors into the direct numerical simulation data for flow past cylinder at Re = 100. As Noca suggested using the low-pass filter to get rid of the high-frequency noise in the evaluated lift force, we verify that his method is inapplicable for dealing with the dataset of 1% experimental error, although the precision is acceptable in practice. To overcome this defect, a novel method is proposed in this paper. The average of the lift forces calculated by using multiple control volume is taken as the evaluation before applying the low-pass filter. The method is applied to an experimental data for flow past a cylinder at approximately Re = 900 to verify its validation. The results show that it improves much better on evaluating the lift forces.     

Pseudopotential multi-relaxation-time lattice Boltzmann model for cavitation bubble collapse with high density ratio

The dynamics of the cavitation bubble collapse is a fundamental issue for the bubble collapse application and prevention. In the present work, the modified forcing scheme for the pseudopotential multi-relaxation-time lattice Boltzmann model developed by Li Q et al. [ Li Q, Luo K H and Li X J 2013 Phys. Rev. E 87 053301] is adopted to develop a cavitation bubble collapse model. In the respects of coexistence curves and Laplace law verification, the improved pseudopotential multi-relaxation-time lattice Boltzmann model is investigated. It is found that the thermodynamic consistency and surface tension are independent of kinematic viscosity. By homogeneous and heterogeneous cavitation simulation, the ability of the present model to describe the cavitation bubble development as well as the cavitation inception is verified. The bubble collapse between two parallel walls is simulated. The dynamic process of a collapsing bubble is consistent with the results from experiments and simulations by other numerical methods. It is demonstrated that the present pseudopotential multirelaxation-time lattice Boltzmann model is applicable and efficient, and the lattice Boltzmann method is an alternative tool for collapsing bubble modeling.     

圆柱绕流问题的初期流动数值模拟

本文继Zheng-huan Teng[1]中介绍的解Navier-stokes方程的椭圆涡团法,研究了一种新的变形涡团法,用以模拟不可压粘性流体绕圆柱的不定常流动。圆柱在静止流体中突然起动并做匀速直线运动。对整个流动区域构造完全Navier-Stokes方程的解并不容易,近十年出现很多数值研究,本文对算法有所推进。把本文的方法称作变形涡团法是因为圆柱边界附近的流体中用椭圆涡团,远离边界时用圆形涡团。计算圆柱绕流比平板绕流在满足附着条件上更为困难,本文分析了怎样在圆柱边界上给出适当的附着条件的数值方法。在算例中雷诺数分别取200、550、3000,得到了不定常边界层分离,二次涡等复杂的物理现象,这些数值结果与近年实验结果[2]是一致的。     

Direct numerical simulation of transitional flow in a finite length curved pipe

Transition to turbulent flow in a curved pipe has been well studied through experiments and numerical simulations. Numerical simulations often use a helical pipe with an infinite length such that the inlet and outlet boundary conditions can be modelled as periodic which greatly reduces computational time. In this study, we examined a finite length curved pipe with Poiseuille flow imposed at the inlet and a stress-free boundary condition at the outlet. Direct numerical simulation of the Navier-Stokes equations for rigid walls and a Newtonian fluid was performed using nek5000. Straight extensions were added to the inlet and outlet such to diminish the impact of boundary conditions on the flow field in the region with curvature. The examined model has a pipe radius of curvature that is three times the pipe radius. The model has ～355 million nodes and required an order of magnitude greater computational time when compared with an infinite length curved pipe. Results show that the critical Reynolds number, the lowest value with instabilities present in the flow, is much greater than that of a straight pipe and occurs near Re=5000–5200. This is larger than the critical Reynolds number typically reported for an infinite length curved pipe (Re=4200–4300).     

纳米尺度圆柱绕流现象的分子动力学模拟研究

本文采用分子动力学模拟方法,利用Lennard-Jones(L-J)势能模型,模拟了Re为28时Ar流体流过Pt纳米圆柱的绕流现象.采用小步长时均值作为涡的瞬时值的方法,在纳米尺度、纳秒量级下得到了涡的周期性产生、组合、发展和脱落现象;在大步长时均条件下,得到了稳定的对称涡,表现出了绕流现象在不同时间范围内的不同特征.绕流过程还体现了流体的密度变化,圆柱上游密度大于下游密度、对称的两侧离轴线越远密度越大.结果表明,分子动力学模拟方法可以有效地细节刻画圆柱绕流现象.     

聚变堆液态包层准二维磁流体湍流模型研究

为研究聚变堆氚增殖包层中液态金属湍流磁流体动力学(MHD)效应,开发了一种新的准二维单方程 MHD 湍流模型,并进行了相关数值模拟程序的编制及验证。对于矩形管道中的准二维 MHD 湍流流动,三维流 动主要发生在哈德曼层中,中心的主流区呈现出二维流动。为了反映这种特殊的流动特征,新湍流模型在标准 k-ε 模型的基础之上去掉了传统的耗散项,代之以电磁耗散项来模拟湍流 MHD 效应。同时,采用 Bradshaw 假设来对 湍流涡粘系数进行模化。为验证该湍流模型是否合理,编制了相关数值模拟程序,并利用直接数值模拟(DNS)结 果对该程序进行了校正,数值模拟结果与 DNS 结果吻合较好。计算结果表明,该湍流模型可应用于聚变堆液态 包层 MHD 湍流流动的数值模拟。     

Hydrodynamic binary coalescence of droplets under air flow in a hydrophobic microchannel

Based on the volume of fluid(VOF) method, we conduct a numerical simulation to study the hydrodynamic binary coalescence of droplets under air flow in a hydrophobic rectangular microchannel. Two distinct regimes, coalescence followed by sliding motion and that followed by detaching motion, are identified and discussed. Additionally, the detailed hydrodynamic information behind the binary coalescence is provided, based on which a dynamic mechanical analysis is conducted to reveal the hydrodynamic mechanisms underlying these two regimes. The simulation results indicate that the sliding motion of droplets is driven by the drag force and restrained by the adhesion force induced by the interfacial tension along the main flow direction. The detachment(i.e., upward motion) of the droplet is driven by the lift force associated with an aerodynamic lifting pressure difference imposed on the coalescent droplet, and also restrained by the adhesion force perpendicular to the main flow direction. Especially, the lift force is mainly induced by an aerodynamic lifting pressure difference imposed on the coalescent droplet. Two typical regimes can be quantitatively recognized by a regime diagram depending on Re and We. The higher Re and We respectively lead to relatively larger lift forces and smaller adhesion forces acting on the droplet, both of which are helpful to detachment of the coalesced droplet.     

二维顶盖驱动半圆腔内流动MRT-LBM研究

本文采用多弛豫时间格子玻尔兹曼方法(multi relaxation time lattice Boltzmann method—MRT LBM)对二维顶盖驱动半圆腔内流动进行了数值模拟,得到了雷诺数为500～50000范围内半圆腔内流场分布情况。在二维顶盖驱动半圆腔流场中,随雷诺数的增大,流场内旋涡的数目逐渐增加,且流动依次呈现出稳定流、周期流、混沌流等状态。本文计算结果表明,MRT-LBM模型可显著提高计算的稳定性,适用于大范围的雷诺数流动情况。     

交叉缩放椭圆管换热与流阻实验研究及分析

对交叉缩放椭圆管进行了实验和数值研究,给出了换热和沿程阻力系数实验拟合关联式。交叉缩放椭圆管管内 截面交叉变化诱导产生强烈的二次流和纵向涡流,改善了速度场与温度场之间的协同关系。实验和数值模拟结果表明,交 叉缩放椭圆管管内的流动在Re≥500即表现为湍流,换热强化效果显著。     

离体射流对S809翼型气动性能的影响

在Re=5×106的条件下,分别在S809翼型前缘点附近不同位置处设置离体射流装置,改变射流动量的大小和射流口宽度,探究其对S809翼型气动性能的影响.并通过流场分析,研究这种流动控制手段有效的物理机理.结果表明:在射流装置位置和射流口宽度固定时,射流动量的大小对控制效果影响显著;在S809翼型表面附近设置微小离体射流...     

变截面折流杆换热器的流动与传热分析

本文在折流杆换热器的基础上,采用数值计算的方法,以水为流动介质,对不同折流杆间隔的折流杆式换热器进行流动与传热特性的分析。同时,本文设计了一种新型的变截面折流杆式换热器,对换热器进行优化,结果表明,变截面折流杆换热器比等截面折流杆换热器的综合性能最多高13%～14%。     

基于LBM方法的风力发电机尾流结构仿真

针对传统CFD数值计算方法难以实现风力机动态旋转及其旋转状态下的流固耦合计算,本文结合格子玻尔兹曼(LBM)方法易于处理动态复杂边界的特点及大涡模拟(LES)方法在非稳态涡流结构捕捉上的优势,采用LBM-LES联合方法进行三维风力发电机整机气动性能及尾流结构仿真研究,同时采用尺度自适应方法对尾涡结构进行跟踪和精细化计算。针对NREL PhaseⅥ型试验机进行模拟,得到了与实验结果吻合的流动形态及尾流结构演变规律,分析了尾流区速度演变规律并对比了不同亚格子湍流模型对计算结果的影响.     

过渡区微尺度流动的有效黏性多松弛系数格子Boltzmann模拟

为提高采用二维九速离散速度模型的格子Boltzmann方法 (LBM)模拟微尺度流动中非线性现象的精度和效率,引入Dongari等提出的有效平均分子自由程对黏性进行修正(Dongari N,Zhang Y H,Reese J M2011 J.Fluids Eng.133 071101);并针对以往研究微尺度流动时采用边界处理格式含有离散误差的问题,采用多松弛系数格子Boltzmann方法结合二阶滑移边界条件,对微尺度Couette流动和周期性Poiseuille流动进行模拟,并将速度分布以及质量流量等模拟结果与直接模拟蒙特卡罗方法模拟数据、线性Boltzmann方程的数值解以及现有的LBM模型模拟结果进行对比.结果表明,相对于现有的LBM模型,引入新的修正函数所建立的有效黏性多松弛系数LBM模型有效提高了LBM模拟过渡区的微尺度流动中的非线性现象的能力.