Flow rate limitation in open wedge channel under microgravity

A study of flow rate limitation in an open wedge channel is reported in this paper. Under microgravity condition, the flow is controlled by the convection and the viscosity in the channel as well as the curvature of the liquid free surface. A maximum flow rate is achieved when the curvature cannot balance the pressure difference leading to a collapse of the free surface. A 1-dimensional theoretical model is used to predict the critical flow rate and calculate the shape of the free surface. Computa- tional Fluid Dynamics tool is also used to simulate the phenomenon. Results show that the 1-dimensional model overestimates the critical flow rate because extra pressure loss is not included in the governing equation. Good agreement is found in 3-dim- ensional simulation results. Parametric study with different wedge angles and channel lengths show that the critical flow rate increases with increasing the cross section area; and decreases with increasing the channel length. The work in this paper can help understand the surface collapsing without gravity and for the design in propellant management devices in satellite tanks.     

Constrained large-eddy simulation and detached eddy simulation of flow past a commercial aircraft at 14 degrees angle of attack

With the development of computational power and numerical algorithms,computational fluid dynamics(CFD) has become an important strategy for the design of aircraft,which significantly reduces the reliance on wind-tunnel and flight tests.In this paper,we conducted a numerical investigation on the flow past a full commercial aircraft at Mach number 0.2 and 14 degrees angle of attack by means of Reynolds-averaged Navier-Stokes(RANS),detached-eddy simulation(DES) and our newly developed constrained large-eddy simulation(CLES).The objective of this paper is to study the capability of these models in simulating turbulent flows.To our knowledge,this is the first large-eddy simulation method for full commercial aircraft simulation.The results show that the CLES can predict the mean statistical quantities well,qualitatively consistent with traditional methods,and can capture more small-scale structures near the surface of the aircraft with massive separations.Our study demonstrates that CLES is a promising alternative for simulating real engineering turbulent flows.     

Beam breakup simulation study for a high energy ERL

The maximum beam current that can be accelerated in an energy recovery linac (ERL) can be severely limited by the transverse multi-pass beam breakup instability (BBU), especially in future ERL light sources with multi-GeV high energy beam energy and more than 100 mA average current. In this paper, the multi-pass BBU of such a high energy ERL is studied based on the simulation of a 3-GeV ERL light source that is proposed by KEK. This work is expected to provide a reference for future high energy ERL projects.     

Simulations of Bingham plastic flows with the multiple-relaxation-time lattice Boltzmann model

Fresh cement mortar is a type of workable paste,which can be well approximated as a Bingham plastic and whose flow behavior is of major concern in engineering.In this paper,Papanastasiou’s model for Bingham fluids is solved by using the multiplerelaxation-time lattice Boltzmann model(MRT-LB).Analysis of the stress growth exponent m in Bingham fluid flow simulations shows that Papanastasiou’s model provides a good approximation of realistic Bingham plastics for values of m108.For lower values of m,Papanastasiou’s model is valid for fluids between Bingham and Newtonian fluids.The MRT-LB model is validated by two benchmark problems:2D steady Poiseuille flows and lid-driven cavity flows.Comparing the numerical results of the velocity distributions with corresponding analytical solutions shows that the MRT-LB model is appropriate for studying Bingham fluids while also providing better numerical stability.We further apply the MRT-LB model to simulate flow through a sudden expansion channel and the flow surrounding a round particle.Besides the rich flow structures obtained in this work,the dynamics fluid force on the round particle is calculated.Results show that both the Reynolds number Re and the Bingham number Bn afect the drag coefcients CD,and a drag coefcient with Re and Bn being taken into account is proposed.The relationship of Bn and the ratio of unyielded zone thickness to particle diameter is also analyzed.Finally,the Bingham fluid flowing around a set of randomly dispersed particles is simulated to obtain the apparent viscosity and velocity fields.These results help simulation of fresh concrete flowing in porous media.     

Simulation of natural convection under high magnetic field by means of the thermal lattice Boltzmann method

The thermal lattice Boltzmann method (TLBM), which was proposed by J. G. M. Eggels and J. A. Somers previously, has been improved in this paper. The improved method has introduced a new equilibrium solution for the temperature distribution function on the assumption that flow is incompressible, and it can correct the effect of compressibility on the macroscopic temperature computed. Compared to the previous method, where the half-way bounce back boundary condition was used for non-slip velocity and temperature, a non-equilibrium extrapolation scheme has been adopted for both velocity and temperature boundary conditions in this paper. Its second-order accuracy coincides with the ensemble accuracy of lattice Boltzmann method. In order to validate the improved thermal scheme, the natural convection of air in a square cavity is simulated by using this method. The results obtained in the simulation agree very well with the data of other numerical methods and benchmark data. It is indicated that the improved TLBM is also successful for the simulations of non-isothermal flows. Moreover, this thermal scheme can be applied to simulate the natural convection in a non-uniform high magnetic field. The simulation has been completed in a square cavity filled with the aqueous solutions of KCl (11wt%), which is considered as a diamagnetic fluid with electrically low-conducting, with Grashof number Gr=4.64× 10^4 and Prandtl number Pr=7.0. And three cases, with different cavity locations in the magnetic field, have been studied. In the presence of a high magnetic field, the natural convection is quenched by the body forces exerted on the electrically low-conducting fluids, such as the magnetization force and the Lorentz force. From the results obtained, it can be seen that the quenching efficiencies decrease with the variation of location from left, symmetrical line, to the right. These phenomena originate from the different distributions of the magnetic field strengths in the zones of the symmetrical central line of the magnetic fields. The results are also compared with those without a magnetic field. Finally, we can conclude that the improved TLBM will enable effective simulation of the natural convection under a high magnetic field.     

Study on the dynamic characteristics of flow over building cluster at high Reynolds number by large eddy simulation

In this paper,the dynamic characteristics of building clusters are simulated by large eddy simulation at high Reynolds number for both homogeneous and heterogeneous building clusters.To save the computational cost a channel-like flow model is applied to the urban canopy with free slip condition at the upper boundary.The results show that the domain height is an important parameter for correct evaluation of the dynamic characteristics.The domain height must be greater than 8h(h is the average building height)in order to obtain correct roughness height while displacement height and roughness sublayer are less sensitive to the domain height.The Reynolds number effects on the dynamic characteristics and flow patterns are investigated.The turbulence intensity is stronger inside building cluster at high Reynolds number while turbulence intensity is almost unchanged with Reynolds number above the building cluster.Roughness height increases monotonously with Reynolds number by 20%from Re*=103 to Re*=105 but displacement height is almost unchanged.Within the canopy layer of heterogeneous building clusters,flow structures vary between buildings and turbulence is more active at high Reynolds number.     

Direct numerical simulation methods of hypersonic flat-plate boundary layer in thermally perfect gas

High-temperature effects alter the physical and transport properties of air such as vibrational excitation in a thermally perfect gas,and this factor should be considered in order to compute the flow field correctly.Herein,for the thermally perfect gas,a simple method of direct numerical simulation on flat-plat boundary layer is put forward,using the equivalent specific heat ratio instead of constant specific heat ratio in the N-S equations and flux splitting form of a calorically perfect gas.The results calculated by the new method are consistent with that by solving the N-S equations of a thermally perfect gas directly.The mean flow has the similarity,and consistent to the corresponding Blasius solution,which confirms that satisfactory results can be obtained basing on the Blasius solution as the mean flow directly in stability analysis.The amplitude growth curve of small disturbance is introduced at the inlet by using direct numerical simulation,which is consistent with that obtained by linear stability theory.It verified that the equation established and the simulation method is correct.     

Coherent structures in wall turbulence and mechanism for drag reduction control

Herein is introduced the mechanism for active control influencing the generation of the near-wall streamwise vortices,which are closely related to the production of high skin friction in wall-bounded turbulent flows.A new opposition control scheme with adjusting control amplitude is proposed and evaluated in turbulent channel flow by direct numerical simulations.The maximum drag reduction rate can be greatly enhanced by the strengthened control.Finally the effectiveness of the control to the coherent structures at high Reynolds numbers is investigated by using a linear transient growth model.     

Modeling the heterogeneous traffic flow considering the effect of self-stabilizing and autonomous vehicles

With the increasing maturity of automatic driving technology,the homogeneous traffic flow will gradually evolve into the heterogeneous traffic flow,which consists of human-driving and autonomous vehicles.To better study the characteristics of the heterogeneous traffic system,this paper proposes a new car-following model for autonomous vehicles and heterogeneous traffic flow,which considers the self-stabilizing effect of vehicles.Through linear and nonlinear methods,this paper deduces and analyzes the stability of such a car-following model with the self-stabilizing effect.Finally,the model is verified by numerical simulation.Numerical results show that the self-stabilizing effect can make the heterogeneous traffic flow more stable,and that increasing the self-stabilizing coefficient or historical time length can strengthen the stability of heterogeneous traffic flow and alleviate traffic congestion effectively.In addition,the heterogeneous traffic flow can also be stabilized with a higher proportion of autonomous vehicles.     

Multiple solutions and hysteresis in the flows driven by surface with antisymmetric velocity profile

Multiple steady solutions and hysteresis phenomenon in the square cavity flows driven by the surface with antisymmetric velocity profile are investigated by numerical simulation and bifurcation analysis.A high order spectral element method with the matrix-free pseudo-arclength technique is used for the steady-state solution and numerical continuation.The complex flow patterns beyond the symmetry-breaking at Re■320 are presented by a bifurcation diagram for Re 2500.The results of stable symmetric and asymmetric solutions are consistent with those reported in literature,and a new unstable asymmetric branch is obtained besides the stable branches.A novel hysteresis phenomenon is observed in the range of2208 Re 2262,where two pairs of stable and two pairs of unstable asymmetric steady solutions beyond the stable symmetric state coexist.The vortices near the sidewall appear when the Reynolds number increases,which correspond to the bifurcation of topology structure,but not the bifurcation of Navier-Stokes equations.The hysteresis is proposed to be the result of the combined mechanisms of the competition and coalescence of secondary vortices.     

利用多松弛格子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的扁椭球区域,数值发现了气泡形状的周期性振荡行为.     

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

综合考虑声学边界层内的热损失和黏性损失,建立处于平面驻波声压波节位置二维球形颗粒外声流计算模型,利用分离时间尺度的数值方法对颗粒外声流流场特征进行模拟。将模拟结果与相应的解析解和实验结果对比,验证了数值模拟的可靠性。在此基础上,研究了雷诺数Re和斯特劳哈尔数Sr对球形颗粒声学边界层内二阶声流流场结构、涡流强度及范围的影响规律。结果表明,随Sr和Re增大,声学边界层内的涡流结构尺度呈指数形式减小,其涡流尺度与颗粒直径D和激励频率f成反比,与流体介质运动黏度v成正比;且满足低Sr和高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.     

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.     

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.     

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]是一致的。