Wavelet analysis of coherent structures in a three-dimensional mixing layer

Wavelet analysis is applied to the results obtained by the direct numerical simulation of a three-dimensional (3D) mixing layer in order to investigate coherent structures in dimension of scale. First, 3D orthonormal wavelet bases are constructed, and the corresponding decomposition algorithm is developed. Then the Navier-Stokes equations are transformed into the wavelet space and the architecture for multi-scale analysis is established. From this architecture, the coarse field images in different scales are obtained and some local statistical quantities are calculated. The results show that, with the development of a mixing layer, the energy spectrum densities for different wavenumbers increase and the energy is transferred from the average flow to vortex structures in different scales. Due to the non-linear interactions between different scales, cascade processes of energy are very complex. Because vortices always roll and pair at special areas, for a definite scale, the energy is obtained from other scales at some areas while it is transferred to other scales at other areas. In addition, energy dissipation and transfer always occur where an intense interaction between vortices exists. The project supported by the Research Fund for the Doctoral Program of Higher Education and the National Natural Science Foundation for Outstanding Youth of China (19925210)     

Nonlinear evolution of Görtler vortices in the near-wall part of a boundary layer

The evolution of Görtler vortices with wavelength smaller than the thickness of the boundary layer on a concave surface is modelled asymptotically at high Reynolds and Görtler numbers. It is known that in the initial linear stage of their evolution such vortices have the largest increment of amplitude growth. Numerical results demonstrate that taking the nonlinear interaction of the flow parameters into account considerably reduces the growth rate and leads to the forming of a perturbed vortex region core; profiles of the flow characteristics in the different stages of vortex evolution are presented.     

Three-dimensional numerical simulation of a vortex ring impacting a bump

A vortex ring impinging on a three-dimensional bump is studied using large eddy simulation for a Reynolds number Re = 4 × 10⁴ based on the initial translation speed and diameter of the vortex ring. The effects of bump height on the vortical flow phenomena and the underlying physical mechanisms are investigated. Based on the analysis of the evolution of vortical structures, two typical kinds of vortical structures, i.e., the wrapping vortices and the hair-pin vortices, are identified and play an important role in the flow state evolution. The circulation of the primary vortex ring reasonably elucidates some typical phases of flow evolution. Furthermore, the mechanism of flow transition from laminar to turbulent state has been revealed based on analysis of turbulent kinetic energy.     

Heat transfer mechanisms during short-pulse laser heating of two-layer composite thin films

 Using a simple perturbation technique, an analytical investigation is presented for the heat transfer mechanisms during ultrafast laser heating of two-layer composite thin slabs from a microscopic point of view. The composite slab consists of two thin metal films which are in perfect thermal contact. The microscopic parabolic two-step model is adopted to describe the behavior of the composite slab. In the microscopic two-step model, the heating process is modeled by the deposition of radiation energy on electrons, the transport of energy by electrons, and the heating of the material lattice through electron–phonon interactions. The proposed perturbation technique is used when the normalized temperature difference between the solid lattice and the electron gas is relatively a small perturbed quantity. Received on 20 September 2000 / Published online: 29 November 2001     

Metastability and Dynamics of Discrete Topological Singularities in Two Dimensions: A Γ-Convergence Approach

This paper aims at building a variational approach to the dynamics of discrete topological singularities in two dimensions, based on Γ-convergence. We consider discrete systems, described by scalar functions defined on a square lattice and governed by periodic interaction potentials. Our main motivation comes from XY spin systems, described by the phase parameter, and screw dislocations, described by the displacement function. For these systems, we introduce a discrete notion of vorticity. As the lattice spacing tends to zero we derive the first order Γ-limit of the free energy which is referred to as renormalized energy and describes the interaction of vortices. As a byproduct of this analysis, we show that such systems exhibit increasingly many metastable configurations of singularities. Therefore, we propose a variational approach to the depinning and dynamics of discrete vortices, based on minimizing movements. We show that, letting first the lattice spacing and then the time step of the minimizing movements tend to zero, the vortices move according with the gradient flow of the renormalized energy, as in the continuous Ginzburg–Landau framework.     

柱形汇聚几何中内爆驱动金属界面不稳定性

采用自研的高保真度爆轰与冲击动力学程序，对柱形汇聚几何中内爆驱动金属材料界面不稳定性的动力学行为，进行了数值模拟研究。结果表明，首次冲击后至约12 μs，界面发展以RM（Richtmyer-Meshkov）不稳定性为主；12 μs后至冲击波聚心反弹加载前，界面聚心运动处于加速减速状态，界面发展由RT (Rayleigh-Taylor)不稳定性主导；冲击波聚心反弹加载后，界面发展又由RM不稳定性主导。另外，还研究了初始条件（初始振幅、初始波长、钢壳初始厚度和几何构型）对柱形内爆驱动金属材料界面不稳定性的影响。结果显示:初始振幅较大时振幅增长也较大；初始波长较小（模数较大）时振幅增长较小，而且存在一个截止波长；钢壳厚度会抑制扰动增长，也存在一个截止厚度；几何汇聚效应会使扰动增长速度更快。     

哈特曼边界层的初级稳定性分析

导电流体在法向外置磁场的作用下,在贴近壁面处会形成哈特曼边界层.哈特曼边界层的稳定性研究对电磁冶金过程和热核聚变反应冷却系统等相关设备的设计和运行都有着十分重要的意义.本文采用非正则模态稳定性分析方法,对两无限大绝缘平行平板内导电流体流动的稳定性进行了研究.通过在时间上迭代求解扰动变量的控制方程组和伴随控制方程组,获得了在磁场作用下初级扰动的增长情况及其空间分布形式,分析了磁场强度对最优扰动增长倍数G_max、最优展向波数β_opt和最优时刻t_opt的影响,并考察了上下两个哈特曼边界层之间的相互作用.结果表明,最优初始扰动的空间分布形式为沿着流场方向的漩涡,关于法向方向对称或者反对称.当哈特曼数Ha较大时(Ha>10),对称漩涡和反对称漩涡形式的初始扰动增长倍数基本相等;上下两个哈特曼边界层可以认为是彼此独立的,不会相互影响,此时最优扰动增长倍数G_max与局部雷诺数R的平方成正比,相应的最优展向波数β_opt和最优时刻t_opt均正比于哈特曼数Ha.当哈特曼数Ha较小时(Ha<10),反对称漩涡形式的初始扰动更为不稳定,其增长倍数大于对称漩涡的增长倍数,且上下两个边界层之间存在着一定的相互作用,并对整个流场的稳定性产生一定的影响.      

Asymmetric vortices flow over slender body and its active control at high angle of attack

The studies of asymmetric vortices flow over slender body and its active control at high angles of attack have significant importance for both academic field and engineering area. This paper attempts to provide an update state of art to the investigations on the fields of forebody asymmetric vortices. This review emphasizes the correlation between micro-perturbation on the model nose and its response and evolution behaviors of the asymmetric vortices. The critical issues are discussed, which include the formation and evolution mechanism of asymmetric multi-vortices; main behaviors of asymmetric vortices flow including its deterministic feature and vortices flow structure; the evolution and development of asymmetric vortices under the perturbation on the model nose; forebody vortex active control especially discussed micro-perturbation active control concept and technique in more detail. However present understanding in this area is still very limited and this paper tries to identify the key unknown problems in the concluding remarks. The project supported by the National Natural Science Foundation of China (10172017), Aeronautical Science Foundation of China (02A51048) and Foundation of National Key Laboratory of Aerodynamic Design and Research (51462020504HK0101)     

圆射流中拟序涡配对的数值模拟

用三维涡丝法结合涡丝的增加与合并技术,对不可压圆射流场,数值模拟了在扰动中加入次谐波后流场拟序涡的配对,配对情况与实验结果一致,通过对不同时刻涡量图的比较与分析,给出了次谐波初始强度对涡配对的影响,初始强度越大,涡配对出现越早,涡配对位置越往上游移动;同时还得到了基波与次谐波初始相位差和涡配对的关系,没有相位差时,次谐波增长最快,相位差增加时,涡配对的速度降低。说明控制次谐波的初始强度和基波与次谐波的相位差,能起到控制圆射流剪切层乃至控制整个射流场的目的。     

Three-dimensional instability analysis of boundary layers perturbed by streamwise vortices

A parametric study is presented for the incompressible, zero-pressure-gradient flat-plate boundary layer perturbed by streamwise vortices. The vortices are placed near the leading edge and model the vortices induced by miniature vortex generators (MVGs), which consist in a spanwise-periodic array of small winglet pairs. The introduction of MVGs has been experimentally proved to be a successful passive flow control strategy for delaying laminar-turbulent transition caused by Tollmien–Schlichting (TS) waves. The counter-rotating vortex pairs induce non-modal, transient growth that leads to a streaky boundary layer flow. The initial intensity of the vortices and their wall-normal distances to the plate wall are varied with the aim of finding the most effective location for streak generation and the effect on the instability characteristics of the perturbed flow. The study includes the solution of the three-dimensional, stationary, streaky boundary layer flows by using the boundary region equations, and the three-dimensional instability analysis of the resulting basic flows by using the plane-marching parabolized stability equations. Depending on the initial circulation and positioning of the vortices, planar TS waves are stabilized by the presence of the streaks, resulting in a reduction in the region of instability and shrink of the neutral stability curve. For a fixed maximum streak amplitude below the threshold for secondary instability (SI), the most effective wall-normal distance for the formation of the streaks is found to also offer the most stabilization of TS waves. By setting a maximum streak amplitude above the threshold for SI, sinuous shear layer modes become unstable, as well as another instability mode that is amplified in a narrow region near the vortex inlet position.     

On strain-rate sensitivity and size effect of brittle solids: transition from cooperative phenomena to microcrack nucleation

An idealized brittle microscale system is subjected to dynamic uniaxial tension in the medium-to-high strain-rate range $(\dot \varepsilon \in \;[100\;s^{-1},\;1 \times 10^{7} \;s^{-1}])$ to investigate its mechanical response under constrained spatial and temporal scales. The setup of dynamic simulations is designed to ensure practically identical in-plane stress conditions on a system of continuum particles forming a two-dimensional, geometrically and structurally disordered, lattice. The rate sensitivity of size effects is observed as well as the ordering effect of kinetic energy. A simple phenomenological expression is developed to account for the tensile strength sensitivity of the small-sized brittle systems to the strain-rate and extrinsic size effects, which may serve as a guideline for formulation of constitutive relations in the MEMS design. The representative sample is defined as a square lattice size for which the tensile strength becomes rate-insensitive and an expression is proposed to model its evolution between two asymptotes corresponding to the limiting loading rates. The dynamics of damage accumulation is analyzed as a function of sample size and loading rate.     

用格子Boltzmann方法数值模拟混合层中旋涡的合并过程

用格子Boltzmann方法计算混合层中的流动问题。在流场的入口处加不同频率、振幅和相位的小扰动,观察混合层中旋涡的演进机理,模拟二维混合层中旋涡合并现象。在基本扰动波的基础上,又加入频率为基本波频率一半的亚谐波,得到了两个涡合并的计算结果,当加入的亚谐波频率为基本波频率的三分之一时,得到了三个涡合并的计算结果。这些计算结果与已有文献的结果基本一致,显示用格子Boltzmann方法模拟混合层问题是可行的。     

Numerical simulation on evolution of subharmonic low-speed streaks in minimal channel turbulent flow

The evolution of low-speed streaks in the turbulent boundary layer of the minimum channel flow unit at a low Reynolds number is simulated by the direct numer- ical simulation （DNS） based on the standard Fourier-Chebyshev spectral method. The subharmonic sinuous （SS） mode for two spanwise-aligned low-speed streaks is excited by imposing the initial perturbations. The possibilities and the physical realities of the turbulent sustaining in the minimal channel unit are examined. Based on such a flow field environment, the evolution of the low-speed streaks during a cycle of turbulent sus- taining, including lift-up, oscillation, and breakdown, is investigated. The development of streamwise vortices and the dynamics of vortex structures are examined. The results show that the vortices generated from the same streak are staggered along the streamwise direction, while the vortices induced by different streaks tilt toward the normal direction due to the mutual induction effect. It is the spatial variations of the streamwise vortices that cause the lift-up of the streaks. By resolving the transport dynamics of enstrophy, the strength of the vortices is found to continuously grow in the logarithmic layer through the vortex stretching mechanism during the evolution of streaks. The enhancement of the vortices contributes to the spanwise oscillation and the following breakdown of the low-speed streaks.     

NUMERICAL STUDY ON THE EFFECTS OF SMALL-AMPLITUDE INITIAL PERTURBATIONS ON RM INSTABILITY

采用Navier-Stokes 方程对入射激波及其反射激波连续诱导小振幅扰动界面的Richtmyer-Meshkov 不稳定性增长过程进行了二维数值模拟,分析了单模和随机多模两种不同初始形态的界面上钉结构和泡结构在反射激波作用前后的发展特性. 研究结果发现：单模扰动的初始界面形态对反射激波前、后界面的扰动增长都有影响,反射激波前的界面形态信息可以通过钉和泡结构之间的反转传递到反射激波过后. 扰动界面上钉结构的发展速度控制了界面混合区总体的发展速度,反射激波前界面上发展成具有完整冠部形态的钉,在反射激波后会反转成复杂的泡结构,此泡结构对反射激波后钉的发展不利. 随机多模界面显示了与单模界面类似的发展规律,但随机多模界面上的复杂泡结构分布的不对称性使得其对钉结构增长的拖曳效应相对要弱,这导致了相似扰动波长下多模随机界面的扰动发展相对单模界面扰动发展要快.     

Direct numerical simulation of free and forced square jets

Direct numerical simulation (DNS) of incompressible, spatially developing square jets in the Reynolds number range of 500–2000 is reported. The three-dimensional unsteady Navier–Stokes equations are solved using high order spatial and temporal discretization. The objective of the present work is to understand the evolution of free and forced square jets by examining the formation of large-scale structures. Coherent structures and related interactions of free jets suggest control strategies that can be used to achieve enhanced spreading and mixing of the jet with the surrounding fluid. The critical Reynolds number for the onset on unsteadiness in an unperturbed free square jet is found to be 875–900 while it reduces to the range 500–525 in the presence of small-scale perturbations. Disturbances applied at the flow inlet cause saturation of KH-instability and early transition to turbulence. Forced jet calculations have been carried out using varicose perturbation with amplitude of 15%, while frequency is independently varied. Simulations show that the initial development of the square jet is influenced by the four corners leading to the appearance hairpin structures along with the formation of vortex rings. Farther downstream, adjacent vortices strongly interact leading to their rapid breakup. Excitation frequencies in the range 0.4–0.6 cause axis-switching of the jet cross-section. Results show that square jets achieve greater spreading but are less controllable in comparison to the circular ones.     

小扰动界面形态对RM不稳定性影响的数值分析

采用Navier-Stokes 方程对入射激波及其反射激波连续诱导小振幅扰动界面的Richtmyer-Meshkov 不稳定性增长过程进行了二维数值模拟,分析了单模和随机多模两种不同初始形态的界面上钉结构和泡结构在反射激波作用前后的发展特性. 研究结果发现:单模扰动的初始界面形态对反射激波前、后界面的扰动增长都有影响,反射激波前的界面形态信息可以通过钉和泡结构之间的反转传递到反射激波过后. 扰动界面上钉结构的发展速度控制了界面混合区总体的发展速度,反射激波前界面上发展成具有完整冠部形态的钉,在反射激波后会反转成复杂的泡结构,此泡结构对反射激波后钉的发展不利. 随机多模界面显示了与单模界面类似的发展规律,但随机多模界面上的复杂泡结构分布的不对称性使得其对钉结构增长的拖曳效应相对要弱,这导致了相似扰动波长下多模随机界面的扰动发展相对单模界面扰动发展要快.      

Analytical solution for temperature field in thin film initially heated by a short-pulse laser source

Analytical solution for electron and lattice temperature distribution in the solid initially heated by a laser short-pulse is presented. Strained parameters method is introduced when formulating electron and lattice temperature distributions. Laser short pulse heating of gold film is simulated numerically and temperature data at the end of the heating pulse are adopted as initial condition to the governing equations of energy transport for analytical solutions. This enables to solve the governing equations of energy analytically in the cooling period. It is found that electron temperature decays sharply while lattice site temperature increases gradually in the surface regions during the cooling cycle. As the depth from the surface increases change in both temperatures become gradual.     

Non-linear and linear evolution of perturbation in stochastic basic flows

In this paper, we examine the non-linear and linear evolutions of perturbation in stochastic basic flows with two-dimensional quasi-geostrophic equations on a sphere. As the analytic solutions for the considered quasi-geostrophic equations are not available, the Fourier finite volume element method is used to perform numerical simulation. It is found that, the non-linear and linear evolutions of perturbation in stochastic basic flow will be consistent for a short period of time and small stochastic fluctuations when they are consistent in the deterministic basic flow. However, the tangent linear model will fail to approximate the original non-linear model when the time period is considerably long and stochastic fluctuation becomes large. Moreover, the global energy decays faster for stochastic basic flow with stronger fluctuations.     

Dynamics of Oscillating Drops with Thermocapillary Effects

The dynamics of oscillating drops warming up in a hot gas environment is investigated via numerical simulation. The formation of surface and internal flows, due to variation of surface tension with temperature, and their impact on the oscillations are discussed. Both surface and ambient temperature disturbances are considered in terms of spherical harmonics. The effects of various parameters including modes of surface and temperature disturbances on period and amplitude of oscillations, kinetic and surface energies, and temperature field are studied. The most obvious feature of thermocapillary flows is demonstrated by vortices whose number and strength varies with the mode of temperature disturbance. These vortices tend to modify the amplitude of oscillations and enhance the kinetic energy. It is also shown that the decrease of the surface tension with increasing temperatures results in the increase of the period of oscillations while decreasing the surface energy. Due to the presence of thermocapillary flows, at long times, the equilibrium shape of the drop is not spherical and the kinetic energy approaches nonzero asymptotic values. The average temperature shows a nearly linear increase in time while the root mean square temperature, used to indicate the spatial variation, levels off after a fast initial growth. Received 3 January 2000 and accepted 7 July 2000     

Receptivity of the Blunt-Leading-Edge Plate Boundary Layer to Unsteady Vortex Perturbations

The response of the boundary layer on a plate with a blunt leading edge to frozen-in vortex perturbations whose vorticity is normal to the plate surface is found. It is shown that these vortices generate an inhomogeneity of the streamwise velocity component in the boundary layer. This inhomogeneity is analogous to the streaky structure developing as the degree of free-stream turbulence increases. The dependence of the amplitude and shape of the boundary layer inhomogeneity on the distance from the leading edge, the streamwise and spanwise scales, and other parameters is found for periodic and local initial perturbations. It is shown that the receptivity of the boundary layer decreases with increase in the frequency and with decrease in the streamwise perturbation scale.