Prediction of particle distribution in isotropic turbulence by large-eddy simulation

The recent application of large-eddy simulation (LES) to particle-laden turbulence requires that the LES with a subgrid scale (SGS) model could accurately predict particle distributions. Usually, a SGS particle model is used to recover the small-scale structures of velocity fields. In this study, we propose a rescaling technique to recover the effects of small-scale motions on the preferential concentration of inertial particles. The technique is used to simulate particle distribution in isotropic turbulence by LES and produce consistent results with direct numerical simulation (DNS). Key words: particle distribution, particle-laden turbulence, large-eddy simulation, subgrid scale model.

Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

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.

Space-time correlations of fluctuating velocities in turbulent shear flows

Space-time correlations or Eulerian two-point two-time correlations of fluctuating velocities are analytically and numerically investigated in turbulent shear flows. An elliptic model for the space-time correlations in the inertial range is developed from the similarity assumptions on the isocorrelation contours: they share a uniform preference direction and a constant aspect ratio. The similarity assumptions are justified using the Kolmogorov similarity hypotheses and verified using the direct numerical simulation DNS of turbulent channel flows. The model relates the space-time correlations to the space correlations via the convection and sweeping characteristic velocities. The analytical expressions for the convection and sweeping velocities are derived from the Navier-Stokes equations for homogeneous turbulent shear flows, where the convection velocity is represented by the mean velocity and the sweeping velocity is the sum of the random sweeping velocity and the shearinduced velocity. This suggests that unlike Taylor’s model where the convection velocity is dominating and Kraichnan and Tennekes’ model where the random sweeping velocity is dominating, the decorrelation time scales of the space-time correlations in turbulent shear flows are determined by the convection velocity, the random sweeping velocity, and the shear-induced velocity. This model predicts a universal form of the spacetime correlations with the two characteristic velocities. The DNS of turbulent channel flows supports the prediction: the correlation functions exhibit a fair good collapse, when plotted against the normalized space and time separations defined by the elliptic model.

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.

H-S管激励振荡流动的数值模拟及实验研究

Hartmann-Sprenger(H-S)管在一定条件下,可以使射流气体产生按一定频率激励的强烈振荡.对H-S管的振荡流动进行了数值模拟和实验研究.数值模拟结果显示出H-S管在吞吐模式下流场的周期变化,在对应的试验瞬态纹影图上,则观察到相应的吞吐射流气体,揭示了H-S管激励的工作机理.通过大量的实验研究,对小管径共振管的频率估算公式进行了修正,提高了振荡频率估算的准确度,为以后设计给定激励频率的射流喷嘴奠定了基础.

Bubble Behavior and Heat Transfer in Quasi-Steady Pool Boiling in Microgravity

Pool boiling of degassed FC-72 on a plane plate heater has been studied experimentally in microgravity. A quasi-steady heating method is adopted, in which the heating voltage is controlled to increase exponentially with time. Compared with terrestrial experiments, bubble behaviors are very different, and have direct effect on heat transfer. Small, primary bubbles attached on the surface seem to be able to suppress the activation of the cavities in the neighborhoods, resulting in a slow increase of the wall temperature with the heat flux. For the high subcooling, the coalesced bubble has a smooth surface and a small size. It is difficult to cover the whole heater surface, resulting in a special region of gradual transitional boiling in which nucleate boiling and local dry area can co-exist. No turning point corresponding to the transition from nucleate boiling to film boiling can be observed. On the contrary, the surface oscillation of the coalesced bubble at low subcooling may cause more activated nucleate sites, and then the surface temperature may keep constant or even fall down with the increasing heat flux. Furthermore, an abrupt transition to film boiling can also be observed. It is shown that heat transfer coefficient and CHF increase with the subcooling or pressure in microgravity, as observed in normal gravity.

Dimensionless Analysis of the Simulation of Bucket Foundations under Dynamic Load

Firstly, the main factors are obtained by use of dimensionless analysis. Secondly, the time scaling factors in centrifuge modeling of bucket foundations under dynamic load are analyzed based on dimensionless analysis and control- ling equation. A simplified method for dealing with the conflict of scaling factors of the inertial and the percolation in sand foundation is presented. The presented method is that the material for experiments is not changed while the effects are modified by perturbation method. Thirdly, the characteristic time of liquefaction state and the characteristic scale of affected zone are analyzed.

Two mechanisms for population inversion via two-colour coherent excitation in three-level systems

We analyse the physical origin of population inversion via continuous wave two-colour coherent excitation in three-level systems by dressing the inverted transition. Two different mechanisms are identified as being responsible for the population inversion. For V-configured systems and cascade (E) configured systems with inversion on the lower transition, the responsible mechanism is the selective trapping of dressed states, and the population inversion approaches the ideal value of 1. For Lambda-configured systems and Xi-configured systems with inversion on the upper transition, population inversion is based on the selective excitation of dressed states, with the population inversion tending towards 0.5. As the essential difference between these two mechanisms, the selective trapping of dressed states occurs in systems with strong decay into dressed states while the selective excitation appears in systems with strong decay out of dressed states.

Angular distribution and conversion of multi-keV L-shell X-ray sources produced from nanosecond laser irradiated thick-foil targets

An experimental study on the angular distribution and conversion of multi-keV X-ray sources produced from 2 ns-duration 527nm laser irradiated thick-foil targets on Shenguang II laser facility (SG-II) is reported. The angular distributions measured in front of the targets can be fitted with the function of f(theta) = alpha+ (1- alpha)cos(beta) theta (theta is the viewing angle relative to the target normal), where alpha = 0.41 +/- 0.014, beta = 0.77 +/- 0.04 for Ti K-shell X-ray Sources (similar to 4.75 keV for Ti K-shell), and alpha = 0.085 +/- 0.06, beta = 0.59 +/- 0.07 for Ag/Pd/Mo L-shell X-ray Sources (2-2.8 keV for Mo L-shell, 2.8-3.5 keV for Pd L-shell, and 3-3.8 keV for Ag L-shell). The isotropy of the angular-distribution of L-shell emission is worse than that of the K-shell emission at larger viewing angle (>70 degrees), due to its larger optical depth (stronger self-absorption) in the cold plasma side lobe Surrounding the central emission region, and in the central hot plasma region (emission region). There is no observable difference in the angular distributions of the L-shell X-ray emission among Ag, Pd, and Mo. The conversion efficiency of Ag/Pd/Mo L-shell X-ray sources is higher than that of the Ti K-shell X-ray sources, but the gain relative to the K-shell emission is not as high as that by using short pulse lasers. The conversion efficiency of the L-shell X-ray sources decrease, with increasing atomic numbers (or X-ray photon energy), similar to the behavior of the K-shell X-ray Source.