LES of Chemical and Acoustic Forcing of a Premixed Dump Combustor

This paper describes the first steps in the development of a large eddy simulation (LES) code able to compute combustion instabilities in gas turbines. This code was used to compute the forcing of an experimentally investigated premixed dump combustor. It is shown that the main effect of acoustic waves entering the combustion chamber is to create large vortices and unsteady heat release when these vortices burn. Another effect of waves entering the combustor is to modulate the fuel and air flow rates produced by the feeding lines. In this case the equivalence ratio of the mixture entering the combustor may also vary. This was investigated in a “chemical effect” simulation where the inlet equivalence ratio fluctuates but the total flow rate remains constant. For perturbations from stoichiometric burning, this mechanism was shown to induce less destabilizing effects than the purely aerodynamical mechanism due to vortex formation and combustion. It is shown that the LES methodology developed is able to reproduce the experimentally observed phase shift between acoustic excitation and total reaction rate in the chamber. This revised version was published online in July 2006 with corrections to the Cover Date.     

Modified Multidimensional Dissipation Scheme on Unstructured Meshes for High-speed Compressible Flow Analysis

A Roe's flux-difference splitting scheme, combining with the entropy fix method according to Van Leer et al., and the H-correction entropy fix method by Pandolfi and D'Ambrosio, is proposed. The presented scheme eliminates unphysical flow behaviors such as a spurious bump of the carbuncle phenomenon that occurs on the bow shock from flow over a blunt body, and the expansion shock generated from flow over a forward facing step. The proposed scheme is further extended to obtain high-order spatial and temporal solution accuracy. The scheme is, in addition, combined with an adaptive meshing technique that generates unstructured triangular meshes to resemble the flow phenomena for reducing computational effort. The entire procedure is evaluated by solving several benchmarks as well as complex steady-state and transient high-speed compressible flow problems.     

Numerical study of bubble rise and interaction in a viscous liquid

In this article, the flow instabilities during the rise of a single bubble in a narrow vertical tube are studied using a transient two-dimensional/axisymmetric model. To predict the shape of the bubble deformation, the Navier-Stokes equations in addition to an advection equation for liquid volume fraction are solved. A modified volume-of-fluid technique based on Youngs' algorithm is used to track the bubble deformation. To validate the model, the results of simulations for terminal rise velocity and bubble shape are compared with those of the experiments. The effect of different parameters such as initial bubble radius, channel height, liquid viscosity and surface tension on the shape and rise velocity of the bubble is investigated.     

Instabilities of micellar systems under homogeneous and non-homogeneous flow conditions

The rheological behavior of a cetylpyridinium chloride 100 mmol l^–1/sodium salicylate 60 mmol l^–1 aqueous solution was studied in this work under homogeneous (cone and plate) and non-homogeneous flow conditions (vane-bob and capillary rheometers), respectively. Instabilities consistent with non-monotonic flow curves were observed in all cases and the solution exhibited similar behavior under the different flow conditions. Hysteresis and the sigmoidal flow curve suggested as characteristic of systems that show constitutive instabilities were observed when running cycles of increasing and decreasing stress or shear rate, respectively. This information, together with a detailed determination of steady states at shear stresses close to the onset of the instabilities, allowed one to show unequivocally that "top and bottom jumping" are the mechanisms to trigger the instabilities in this micellar system. It is shown in addition that there is not a true plateau region in between the "top and bottom jumping". Finally, the flow behavior beyond the upturn seemed to be unstable and was found accompanied by an apparent violation of the no-slip boundary condition.     

Bistable behaviour in squeezed vacua: II. Stability analysis and chaos

Linear stability analysis and (numerical) investigation of the periodic and chaotic self-pulsing behaviour are presented for the Maxwell-Bloch equations of a bistable model in contact with a squeezed vacuum field. Effect of the squeeze phase parameter on the period doubling bifurcation that preceeds chaos is examined for the adiabatic and non-adiabatic regimes.     

“Slow” active control of combustion instabilities by modification of liquid fuel spray properties

This paper describes an experimental investigation of the feasibility of using “slow” active control approaches, which “instantaneously” change liquid fuel spray properties, to suppress combustion instabilities. The objective of this control approach was to break up the feedback between the combustion process heat release and combustor pressure oscillations that drive the instability by changing the characteristics of the combustion process (e.g., the characteristic combustion time). To demonstrate the feasibility of such control, this study used a proprietary fuel injector (Nanomiser^TM), which can vary its fuel spray properties, to investigate the dependence of acoustics–combustion process coupling, i.e., the driving of combustion instabilities, upon the fuel spray properties. This study showed that by changing the spray characteristics it is possible to significantly damp combustion instabilities. Furthermore, using combustion zone chemiluminescence distributions, which were obtained by Abel’s deconvolution synchronized with measured acoustic data, it has been shown that the instabilities were mostly driven midway between the combustor centerline and wall, a short distance downstream from the flame holder, where the mean axial flow velocity is approximately zero in the vortex near the flame holder. The results of this study strongly suggest that a “slow” active control system that employs controllable fuel injectors could be effectively used to prevent the onset of detrimental combustion instabilities.     

合肥光源逐束团相位测量及纵向不稳定性诊断

为了研究储存环束流多束团纵向运动特性,在合肥光源（HLS-II）上研制了一套逐束团相位测量系统。该系统利用示波器直接采集BPM和信号,采用过零点检测法、时间差分（Temporal Difference, TD）法相结合的方法从BPM和信号中提取出逐束团相位。介绍了逐束团相位测量系统的系统构架、相位提取方法以及在HLS-II上的一些实验结果。通过对该系统记录的5 ms时间长度的逐束团多圈相位数据的离线分析,得到多束团纵向运动的同步振荡的频率、纵向工作点、束团振荡的模式信息以及振荡模式增长率等特征信息,诊断出HLS-II在top-off恒流运行期间存在2个较强的纵向耦合束团不稳定模式,并提取出了2个振荡模式的增长率。该系统的逐束团相位测量结果及相关纵向不稳定性分析可为机器研究、纵向反馈系统调试评估及高频RF系统的性能评估等提供参考。     

Massively parallel LES of azimuthal thermo-acoustic instabilities in annular gas turbines

Increasingly stringent regulations and the need to tackle rising fuel prices have placed great emphasis on the design of aeronautical gas turbines, which are unfortunately more and more prone to combustion instabilities. In the particular field of annular combustion chambers, these instabilities often take the form of azimuthal modes. To predict these modes, one must compute the full combustion chamber, which remained out of reach until very recently and the development of massively parallel computers. In this article, full annular Large Eddy Simulations (LES) of two helicopter combustors, which differ only on the swirlers' design, are performed. In both computations, LES captures self-established rotating azimuthal modes. However, the two cases exhibit different thermo-acoustic responses and the resulting limit-cycles are different. With the first design, a self-excited strong instability develops, leading to pulsating flames and local flashback. In the second case, the flames are much less affected by the azimuthal mode and remain stable, allowing an acceptable operation. Hence, this study highlights the potential of LES for discriminating injection system designs. To cite this article: P. Wolf et al., C. R. Mecanique 337 (2009).     

Combination of two basic types of synchronization in a coupled semiconductor laser system

Combination of two basic types of synchronization, anticipatory synchronization and lagged synchronization, is investigated numerically between two coupled semiconductor lasers. It is found that lagged synchronization produced by a backward coupling with a suitable delay can combine with the originally hidden anticipatory synchronization and produce a type of synchronization overcoming the original lagged synchronization produced by a forward coupling. We study the combination synchronization phenomenon when the delay of the backward coupling is different from that of the original anticipatory synchronization. Our results suggest that the synchronization combination phenomenon might allow an interpretation of an experimental observation by Sivaprakasam et al. [Phys. Rev. Lett. 87, 154101 (2001)] that the anticipating time is irrespective of the external-cavity round trip time, which to date remains to be understood.