A passive scalar field convected by turbulence

Classical statistical mechanics is applied to the study of a passive scalar field convected by isotropic turbulence. A complete set of independent real parameters and dynamic equations are worked out to describe the dynamic state of the passive scalar field. The corresponding Liouville equation is solved by a perturbation method based upon a Langevin–Fokker–Planck model. The closure problem is treated by a variational approach reported in earlier papers. Two integral equations are obtained for two unknown functions: the scalar variance spectrum F(k) and the effective damping coefficient (k). The appearance of the energy spectrum of the velocity field in the two integral equations represents the coupling of the scalar field with the velocity field. As an application of the theory, the two integral equations are solved to derive the inertial-convective-range spectrum, obtaining F(k)=0.61 −1/3 k−5/3. Here is the dissipation rate of the scalar variance and is the dissipation rate of the energy of the velocity field. This theoretical value of the scalar Kolmogorov constant, 0.61, is in good agreement with experiments.

Nonequilibrium statistical mechanics of two-dimensional turbulence

The vorticity dynamics of two-dimensional turbulence are investigated analytically, applying the method of Qian (1983). The vorticity equation and its Fourier transform are presented; a set of modal parameters and a modal dynamic equation are derived; and the corresponding Liouville equation for the probability distribution in phase space is solved using a Langevin/Fokker-Planck approach to obtain integral equations for the enstrophy and for the dynamic damping coefficient eta. The equilibrium spectrum for inviscid flow is found to be a stationary solution of the enstrophy equation, and the inertial-range spectrum is determined by introducing a localization factor in the two integral equations and evaluating the localized versions numerically.

Variational approach to the closure problem of turbulence theory

A new method is proposed to solve the closure problem of turbulence theory and to drive the Kolmogorov law in an Eulerian framework. Instead of using complex Fourier components of velocity field as modal parameters, a complete set of independent real parameters and dynamic equations are worked out to describe the dynamic states of a turbulence. Classical statistical mechanics is used to study the statistical behavior of the turbulence. An approximate stationary solution of the Liouville equation is obtained by a perturbation method based on a Langevin-Fokker-Planck (LFP) model. The dynamic damping coefficient eta of the LFP model is treated as an optimum control parameter to minimize the error of the perturbation solution; this leads to a convergent integral equation for eta to replace the divergent response equation of Kraichnan's direct-interaction (DI) approximation, thereby solving the closure problem without appealing to a Lagrangian formulation. The Kolmogorov constant Ko is evaluated numerically, obtaining Ko = 1.2, which is compatible with the experimental data given by Gibson and Schwartz, (1963).

快速凝固/粉末冶金高阻尼铝合金挤压成形过程的数值模拟

采用有限元分析软件ANSYS,分析了快速凝固/粉末冶金(Rapidly Solidified/Powder Metallurgy,RS/PM)高阻尼铝合金FMS0714/10(Zn-30Al)的挤压成形过程,研究了模具与坯料间的摩擦条件和挤压比对阻尼铝合金成形过程的影响,探讨了挤压过程材料表面产生裂纹的机理,进行了成形工艺优化.研究结果表明:挤压过程中FMS0714/10(Zn-30Al)阻尼铝合金材料的流动与变形不均匀,导致了应力应变分布的不均匀:材料表面应力应变较大,芯部应力应变较小;表面过大的应变和轴向拉应力是表面裂纹的诱因;减小模具和坯料间的摩擦,增大挤压比,可以减小材料表面应力,使应变分布趋于均匀,从而减少材料表面损伤,优化材料表面质量,提高成品合格率.数值模拟的研究结果将为FMS0714/10(Zn-30Al)阻尼铝合金材料挤压工艺的制订以及新材料的设计和研制提供有益的参考.

弱粘性流体中垂直激励表面波的阻尼效应

在竖直振动的圆柱形容器中，将Navier-Stokes方程线性化，利用两时间尺度奇异摄动展开法研究了弱粘性流体的单一自由面驻波运动．整个流场被分为外部势流区和内部边界层区两部分，对两部分区域分别求解，得到包含阻尼项和外驱动影响的线性振幅方程．利用稳定性分析，得到形成稳定表面波的条件，给出了临界曲线．此外，还获得了阻尼系数的解析表达式．最后，将线性阻尼加到理想流体条件下所得到的色散关系中对其进行修正，理论结果证明修正后的驱动频率更加接近实验的结果．通过计算发现，当驱动的频率较低时，流体的粘性对表面波模式选择有重要影响，而表面张力的影响不明显；但当驱动频率较高时，流体的表面张力起主要作用，而流体的粘性影响甚小．

微梁谐振器中的空气阻尼

为了研究空气阻尼对谐振器工作性能的影响,基于Euler-Bernoulli假设、梁弯曲振动理论和气体阻尼理论,建立了考虑气体阻尼的梁振动方程,结合边界条件对方程进行了求解,得到了两种不同情况下气体阻尼对微梁谐振器的谐振频率和品质因子的影响.通过两种阻尼对谐振器影响的比较分析,指出在梁形状保持不变(长宽比和高宽比不变)情况下,微梁谐振器存在一个临界压膜厚度,此临界厚度可为器件的设计提供参考.

A New Model for Sediment Transport

The problem of predicting sediment transportation by water waves is treated analytically with the rate of wave energy dissipation or wave damping. With resorting to the theory of shallow water waves and the basis of Yamamoto’s Coulomb-damped poroelastic model, the Boussinesq-type equation has been derived over a variation depth bed. For convenience Cnoidal wave is just discussed, The Cnoidal wave with complex wave length and wave velocity, which are as a function of wave frequency, water depth, permeability, Poisson’s ratio and complex elastic moduli of bed soil, is applied to analyse the rate of sediment transportation. Considering the sediment transportation depended on the shear stress near-bed or the horizontal velocity, the conclusion of Yamamoto’s experiment in clay bed has been extended to general situation. It could be figured out that the model should provide a method to avoid the undistinguishable factors during sediment transport processes and relate mass transport with the sediment peculiarities.

Pool Boiling in Microgravity: Recent Results and Perspectives for the Project DEPA-SJ10

Two research projects on pool boiling in microgravity have been conducted aboard the Chinese recoverable satellites. Ground-based experiments have also been performed both in normal gravity and in short-term microgravity in the Drop Tower Beijing. Steady boiling of R113 on thin platinum wires was studied with a temperature-controlled heating method, while quasi-steady boiling of FC-72 on a plane plate was investigated with an exponentially increasing heating voltage. In the first case, slight enhancement of heat transfer is observed in microgravity, while diminution is evident for high heat flux in the second one. Lateral motions of bubbles on the heaters are observed before their departure in microgravity. The surface oscillation of the merged bubbles due to lateral coalescence between adjacent bubbles drives it to detach from the heaters. The Marangoni effect on the bubble behavior is also discussed. The perspectives for a new project DEPA-SJ10, which has been planned to be flown aboard the Chinese recoverable satellite SJ-10 in the future, are also presented.

Two bifurcation processes for onset of oscillatory thermocapillary convection in a floating half zone

An overview on the onset of thermocapillary oscillatory convection in a floating half zone is provided, and it is a typical subject in the microgravity sciences related to the space materials science, especially the floating zone processing, and also to the microgravity fluid physics. The main interests are focused around the process for onset of oscillatory thermocapillary convection, which is known also as the bifurcation transition from quasi-steady convection to oscillatory convection. The onset of oscillation depends on a set of critical parameters, such as the Marangoni number, Prandtl number, geometrical parameters, and heat transfer parameters. Recent studies show that, there exists the bifurcation transition from steady and axial symmetric convection to the steady and axial non-symmetric convection before the onset of oscillation in cases of small Prandtl number fluids and in cases of larger Prandtl number fluids of fat liquid bridge with small aspect ratio. The transition process is a strong non-linear process because the velocity deviation has the same order of magnitude as that of an average flow after the onset of oscillation, and unsteady 3-D numerical simulation is suitable to do in depth analysis on strong non-linear process, and leads generally to a better comparison with the experimental results.

Floating zone convection

The microgravity research, as a branch of the advanced sciences and a spe- cialized field of high technology, has been made in China since the late 1980's. The research group investigating microgravity fluid physics consisted of our col- leagues and the authors in the Institute of Mechanics of the Chinese Academy of Sciences (CAS), and we pay special attention to the floating zone convection as our first research priority. Now, the research group has expanded and is a part of the National Microgravity Laboratory of the CAS, and the research fields have been extended to include more subjects related to microgravity science. Howev- er, the floating zone convection is still an important topic that greatly holds our research interests.

均匀来流绕旋转及旋转振荡圆柱绕流研究

该文通过数值方法求解二维不可压Navier-Stokes方程，对均匀来流中静止、旋转和旋转振荡圆柱绕流进行了系统的数值模拟.该文采用有限体积法对控制方程进行离散，选用元结构化四边形网格剖分计算区域，关于速度-压力耦合的处理使用了SIMPLEC方法.经过了大量的数值模拟，分盺166L鸬玫搅苏饧钢秩屏鞯氖的Ｄ饨峁?该文重点是用快速傅里叶变换（FFT）方法对旋转振荡圆柱绕流中的频率耦合现象进行研究，并分析在不同频率耦合作用下涡形成、发展和脱落的规律.

Instability Analysis Of Torsional Mems/Nems Actuators Under Capillary Force

The static and dynamic instabilities of a torsional MEMS/NEMS actuator caused by capillary effects are studied, respectively. An instability number, eta, is defined, and the critical gap distance, g(cr), between the mainplate and the substrate is derived. According to the values of eta and g, the instability criteria of the actuator are presented. The dimensionless motion equation of the MEMS/NEMS torsional actuator is derived when it makes nonlinear oscillation under capillary force. The qualitative analysis of the nonlinear equation is made, and the phase portraits are presented on the phase plane. In addition, the bifurcation phenomena in the system are also analyzed. (C) 2008 Elsevier Inc. All rights reserved.

Effects of geometric shape on the hydrodynamics of a self-propelled flapping foil

The hydrodynamics of a free flapping foil is studied numerically. The foil undergoes a forced vertical oscillation and is free to move horizontally. The effect of chord-thickness ratio is investigated by varying this parameter while fixing other ones such as the Reynolds number, the density ratio, and the flapping amplitude. Three different flow regimes have been identified when we increase the chord-thickness ratio, i.e., left-right symmetry, back-and-forth chaotic motion, and unidirectional motion with staggered vortex street. It is observed that the chord-thickness ratio can affect the symmetry-breaking bifurcation, the arrangement of vortices in the wake, and the terminal velocity of the foil. The similarity in the symmetry-breaking bifurcation of the present problem to that of a flapping body under constraint is discussed. A comparison between the dynamic behaviors of an elliptic foil and a rectangular foil at various chord-thickness ratios is also presented.

AUTOROTATION OF TWO TANDEM TRIANGULAR CYLINDERS

The autorotation of two tandem triangular cylinders at different gap distances is investigated by numerical simulations. At the Reynolds number of 200， three distinct regimes are observed with the increase of gap distance: namely, angular oscillation, quasi-periodic autorotation and ‘chaotic’ autorotation. For various gap distances, the characteristic of vortex shedding and vortex interaction are discussed. The phase graphs (angular acceleration vs. angular velocity) and the power spectra of moment are analyzed to characterize the motion of the cylinder. The Lyapunov exponent is also calculated to identify the existence of chaos.

Two-dimensional Vortex-induced Vibrations of Submerged Floating Tunnel Tethers

A vortex-induced vibration (VIV) model is presented for predicting the nonlinear dynamic response of submerged floating tunnel (SFT) tethers which are subjected to wave, current and tunnel oscillatory displacements at their upper end in horizontal and vertical directions. A nonlinear fluid force formula is introduced in this model, and the effect of the nonlinearity of tether is investigated. First, the tunnel is stationary and the tether vibrates due to the vortices shedding. The calculated results show that the cross-flow amplitude of VIV decreases compared with the linear model. However the in-line amplitude of VIV increases. Next, the periodical oscillation of tunnel is considered. The oscillation caused by wave forces plays the roles of parametric exciter and forcing exciter to the VIV of tether. The time history of displacement of the tether mid-span is obtained by the proposed model. It is shown that the in-line amplitude increases obviously and the corresponding frequency is changed. The cross-flow amplitude exhibits a periodic behavior.