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)     

Wavy structures in compressible mixing layers

Semi-periodic structures namely inclined wavy structures(IWS) are experimentally observed in compressible mixing layers at two convective Mach numbers(Mc=0.11 and 0.47). Flow structures are visualized by the laserinduced planar laser Mie scattering(PLMS) technique. Two methods are developed to investigate the spatial distribution and geometry of IWS:(1) the dominant mode extraction(DME) method, to extract the dominant modes of IWS from the streamwise gray-level fluctuation, and(2) the phase tracking(PT) method, to identify the shape of IWS. The results suggest that pressure perturbations account for the formation of IWS in the initial mixing region and the joint efect of dilatation and coherent vortices enhances IWS in the welldeveloped region. The large transverse(cross-flow) scale of the IWS and their relation to coherent vortices(CV) indicate that the disturbance originated from CV in the mixing center propagates far into the free streams. The DME and the PT method are shown to be the efective tools to study the geometrical features of wavy structures in compressible shear flows.     

An experimental study of the large scale coherent structures in a forced free shear layer

The dynamic characteristics of the large scale coherent structures in a forced free shear layer are experimentally studied by means of flow visualization. The quantitative measurements are acquired by the use of a LDV. It is shown that the development of the coherent structures can be greatly influenced by upstream artificial perturbations and as a result the mixing in the layer can be controlled. Like vortex merging, vortex splitting is also a common evolution pattern in the development of the coherent structures.     

Numerical research on the coherent structure in the viscoelastic second-order mixing layers

I.IntroductionViscoelastict'olwsareprevalentinstiring,mixingandchemicalreactionofdilutepolymersolutions.SinceTomsdiscoveredthedragreducingeffectofpolylllerudditivesontul.bulent11owsin1948,1llilnystudiesonitsmechanisms11acebeendoneinviewofitsenormouspr'tcticillalld[lleorcticill\'al[Ics.UPtodata,1ilorcasonilbicilltcrl,rclilti')lls,II()wcvcr,ill'cstillun'lvtlilabjedLlctotiledilllcultiesillfoeviscoelasticdyllalnicsalldthetLlrbulellcetllcory.Al'cwpolymericadditiveshavinglittleeffectonthesolvent…     

Numerical simulations of compressible mixing layers with a discontinuous Galerkin method

Discontinuous Galerkin(DG) method is known to have several advantages for flow simulations,in particular,in fiexible accuracy management and adaptability to mesh refinement. In the present work,the DG method is developed for numerical simulations of both temporally and spatially developing mixing layers. For the temporally developing mixing layer,both the instantaneous fiow field and time evolution of momentum thickness agree very well with the previous results. Shocklets are observed at higher convective Mach numbers and the vortex paring manner is changed for high compressibility. For the spatially developing mixing layer,large-scale coherent structures and self-similar behavior for mean profiles are investigated. The instantaneous fiow field for a three-dimensional compressible mixing layer is also reported,which shows the development of largescale coherent structures in the streamwise direction. All numerical results suggest that the DG method is effective in performing accurate numerical simulations for compressible shear fiows.     

Formation and evolution of the streamwise vortices in mixing layers

The evolution of the three-dimensional time-developing mixing layer is simulated numerically using the pseudo-spectral method. The initial perturbations used in this study consisted of the two-dimensional fundamental wave and the streamwise-invariant three-dimensional disturbance. A comparison of the formations of the streamwise vortices with different amplitude functions for three-dimensional disturbances is made. In one case the results are similar to that of Rogers & Moser^[1], whereas a different way in which the quadrupole forms and sudden expansion of the rib are observed in another case. The simulation also confirms that the stretching by the forming roller rather than Rayleigh centrifugal instability is responsible for the formation of the rib. Finally, numerical flow visualization results are presented. The project supported by the Zhejiang Province Natural Science Special Fund for Youth Scientistis' Cultivation.     

Spanwise domain effects on the evolution of the plane turbulent mixing layer

Large Eddy Simulation is used to simulate a series of plane mixing layers. The influence of the spanwise domain on the development of the mixing layer, and the evolution of the coherent structures, are considered. The mixing layers originate from laminar conditions, and an idealised inflow condition is found to produce accurate flow predictions when the spanwise computational domain extent is sufficient to avoid confinement effects. Spanwise domain confinement of the flow occurs when the ratio of spanwise domain extent to local momentum thickness reaches a value of ten. Flow confinement results in changes to both the growth mechanism of the turbulent coherent structures, and the nature of the interactions that occur between them. The results demonstrate that simulations of the two-dimensional mixing layer flow requires a three-dimensional computational domain in order that the flow will evolve in a manner that is free from restraints imposed by the spanwise domain.     

A comparative study of inflow conditions for two‐ and three‐dimensional spatially developing mixing layers using large eddy simulation

A series of spatially developing mixing layers are simulated using the large eddy simulation (LES) technique. A hyperbolic tangent function and data derived from boundary layer simulations are used to generate the inflow condition, and their effects on the flow are compared. The simulations are performed in both two and three dimensions. In two‐dimensional simulations, both types of inflow conditions produce a layer that grows through successive pairings of Kelvin–Helmholtz (K–H) vortices, but the composition ratio is lower for the hyperbolic tangent inflow simulations. The two‐dimensional simulations do not undergo a transition to turbulence. The three‐dimensional simulations produce a transition to turbulence, and coherent structures are found in the post‐transition region of the flow. The composition ratio of the three‐dimensional layers is reduced in comparison to the counterpart two‐dimensional runs. The mechanisms of growth are investigated in each type of simulation, and amalgamative pairing interactions are found in the pre‐transition region of the three‐dimensional simulations, and throughout the entire computational domain of those carried out in two‐dimensions. The structures beyond the post‐transition region of the three‐dimensional simulations appear to behave in a much different manner to their pre‐transition cousins, with no pairing‐type interactions observed in the turbulent flow. In order to accurately simulate spatially developing mixing layers, it is postulated that the inflow conditions must closely correspond to the conditions present in the reference experiment. Copyright © 2007 John Wiley & Sons, Ltd.     

Secondary instability of large scale structure in free turbulent shear layer

The secondary instability theory is used to study the behavior of spatially growingdisturbance in free turbulent shear layer.The numerical results indicate that secondaryinstability of subharmonic mode shows a strong choice of spanwise wavenumber andthe maximum growth two dimensional case.In contrast to that secondary instabilities of the fundamental mode occur in a wide scope of spanwise wavenumber.We have found so called translative at β=0 and bifurcation phenomenon for an amplitude of the KH wave larger than 0.06.Dey words instability,large scale structure,bifurcation     

Coherent fine scale eddies in turbulence transition of spatially-developing mixing layer

To investigate the relationship between characteristics of the coherent fine scale eddy and a laminar–turbulent transition, a direct numerical simulation (DNS) of a spatially-developing turbulent mixing layer with Re_ω,0 = 700 was conducted. On the onset of the transition, strong coherent fine scale eddies appears in the mixing layer. The most expected value of maximum azimuthal velocity of the eddy is 2.0 times Kolmogorov velocity (u_k), and decreases to 1.2u_k, which is an asymptotic value in the fully-developed state, through the transition. The energy dissipation rate around the eddy is twice as high compared with that in the fully-developed state. However, the most expected diameter and eigenvalues ratio of strain rate acting on the coherent fine scale eddy are maintained to be 8 times Kolmogorov length (η) and :β:γ = −5:1:4 in the transition process. In addition to Kelvin–Helmholtz rollers, rib structures do not disappear in the transition process and are composed of lots of coherent fine scale eddies in the fully-developed state instead of a single eddy observed in early stage of the transition or in laminar flow.     

On decentralized stabilization of linear large scale systems with symmetric circulant structure

The decentralized stabilization of continuous and discrete linear large scale systems with symmetric circulant structure was studied. A few sufficient conditions on decentralized stabilization of such systems were proposed. For the continuous systems, by introducing a concept called the magnitude of interconnected structure, a very important property that the decentralized stabilization of such systems is fully determined by the structure of each isolated subsystem that is obtained when the magnitude of interconnected structure of the overall system is given. So the decentralized stabilization of such systems can be got by only appropriately designing or modifying the structure of each isolated subsystem, no matter how complicated the interconnected structure of the overall system is. A algorithm for obtaining decentralized state feedback to stabilize the overall system is given. The discrete systems were also discussed. The results show that there is a great dfference on decentralized stabilization between continuous case and discrete case.     

Eulerian-Lagranigan simulation of aerosol evolution in turbulent mixing layer

The formation and evolution of aerosol in turbulent flows are ubiquitous in both industrial processes and nature. The intricate interaction of turbulent mixing and aerosol evolution in a canonical turbulent mixing layer was investigated by a direct numerical simulation (DNS) in a recent study (Zhou, K., Attili, A., Alshaarawi, A., and Bisetti, F. Simulation of aerosol nucleation and growth in a turbulent mixing layer. Physics of Fluids, 26, 065106 (2014)). In this work, Monte Carlo (MC) simulation of aerosol evolution is carried out along Lagrangian trajectories obtained in the previous simulation, in order to quantify the error of the moment method used in the previous simulation. Moreover, the particle size distribution (PSD), not available in the previous works, is also investigated. Along a fluid parcel moving through the turbulent flow, temperature and vapor concentration exhibit complex fluctuations, triggering complicate aerosol processes and rendering complex PSD. However, the mean PSD is found to be bi-modal in most of the mixing layer except that a tri-modal distribution is found in the turbulent transition region. The simulated PSDs agree with the experiment observations available in the literature. A different explanation on the formation of such PSDs is provided.     

超声速平面混合层小激波的形成与演变

为了揭示超声速混合层中小激波形成机理及其与涡相互作用的演变过程,本文基于大涡模拟（LES）方法,结合五阶精度混合TCD／WENO格式,对超声速平面混合层在对流马赫数为Mc=0．65条件下的流场结构进行了数值模拟,数值结果详细描述了超声速混合层中小激波的形成过程。研究了小激波形成后,随涡运动而产生的变形、脱落及发展过程。同时,对混合层双涡合并过程中,小激波与相邻涡相互作用所产生的变形与演变过程进行了讨论。     

粘性可压混合层时间稳定性对称紧致差分求解

基于可压扰动方程组的一阶改型 ,将高精度对称紧致格式引入边值法数值线性稳定性分析。对所获非线性离散特征值问题给出了一个通用形式二阶迭代局部算法 ,实现了时间模式和空间模式的统一求解 ,并将扰动特征值及其特征函数同时得到。据此分析了可压平面自由混合层时间稳定性 ,涉及二维 /三维扰动波、粘性 /无粘扰动波、第一 /第二模态、特征函数、伪特征值谱等。研究表明 ,压缩性效应和粘性效应对最不稳定扰动波数和增长率呈相似的减抑作用 ;在 Mc=1附近 ,从高波数段开始 ,粘性效应可强化二维不稳定扰动波由第一模态向第二模态的过渡     

On problems in the weakly nonlinear theory of hydrodynamic stability and its improvement

There are three main problems in the weakly nonlinear theory of hydrodynamic stability: (1)The radius of convergence with respect to the perturbation parameter is too small and there is no concrete estimation for it. (2)The solution has a special structure, thus in general, it can not satisfy the initial condition posed by many practical problems. (3)When the linear part of its solution does not correspond to a neutral case, there are more than one way in determining the Landau constants, and practically no one knows which is the best way. In this paper, problems(1) and (2)are solved theoretically, and ways for its improvement have been proposed. By comparing the theoretical results with those obtained by numerical simulations, problem(3)has also been clarified. Project supported by the National Natural Science Foundation of China     

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

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

混合层流动拟序结构的大涡模拟

采用大涡模拟方法对空间发展的二维平面混合层进行了数值模拟 ,动量方程采用分步投影法求解 ,亚格子项采用标准Smagorinsky亚格子模式模拟 ,压力泊松方程采用修正的循环消去法快速求解 ,同时求解了标志物输运方程以实现数值流场显示。模拟结果给出了混合层流动的瞬态发展过程以及流动中拟序结构的发展演变过程 ,成功地模拟了混合层发展中的各种瞬态细节过程 ,如涡的卷起、增长 ,涡与涡之间的配对、合并过程 ,以及大涡破碎为小涡的级联过程 ,为各种以混合层流动为原型流动的射流、尾流等工业流动的控制和优化提供了理论基础。     

A parallel monolithic algorithm for the numerical simulation of large‐scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large‐scale fluid structure interaction problems. The governing incompressible Navier–Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian–Eulerian formulation‐based side‐centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant–Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large‐scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned sub‐domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd.     

基于修正型罗德里格斯参数的三维刚体摆姿态控制

3D 刚体摆是研究地球静止轨道航天器的一个力学简化模型, 它绕一个固定、无摩擦的支点旋转, 具有3 个转动自由度. 文章给出基于修正型罗德里格斯(Rodrigues) 参数描述的3D 刚体摆的姿态动力学方程, 针对3D 刚体摆姿态和角速度稳定的非线性控制设计问题, 基于无源性控制理论利用能量法设计了3D 刚体摆的系统控制器, 并证明了系统满足无源性. 构造了系统的李雅普诺夫(Lyapunov) 函数, 利用能量法设计出3D 刚体摆的姿态控制律, 并由拉萨尔(LaSalle) 不变集原理证明了该控制律的渐近稳定性. 仿真实验给出了3D 刚体摆在倒立平衡位置的姿态和角速度的渐近稳定性, 仿真实验结果表明基于能量方法的3D 刚体摆姿态控制是有效的.     

新修正偶应力层合板理论中的几个基本假设对自由振动计算的影响

基于新修正偶应力理论,在对微细观尺度的复合材料层合梁/板进行力学响应计算时,往往采用一系列假设来简化模型。现有文献都全部或部分应用了这些假设,但对这些假设是否会对计算结果造成影响尚未进行充分讨论分析。本文建立了未经简化的新修正偶应力Reddy层合板模型,并对其自由振动进行了分析。通过数值算例的对比,讨论了常用的几个简化假设对微细观复合材料四边简支方板自振频率的影响以及适用范围。算例结果表明,常用的几个简化假设对于微尺度层合薄板自由振动的影响很小,对于厚板的低阶频率影响也很小,但对厚板的高阶频率影响显著。