Specific heat ratio effects of compressible Rayleigh–Taylor instability studied by discrete Boltzmann method

Rayleigh–Taylor (RT) instability widely exists in nature and engineering fields. How to better understand the physical mechanism of RT instability is of great theoretical significance and practical value. At present, abundant results of RT instability have been obtained by traditional macroscopic methods. However, research on the thermodynamic non-equilibrium (TNE) effects in the process of system evolution is relatively scarce. In this paper, the discrete Boltzmann method based on non-equilibrium statistical physics is utilized to study the effects of the specific heat ratio on compressible RT instability. The evolution process of the compressible RT system with different specific heat ratios can be analyzed by the temperature gradient and the proportion of the non-equilibrium region. Firstly, as a result of the competition between the macroscopic magnitude gradient and the non-equilibrium region, the average TNE intensity first increases and then reduces, and it increases with the specific heat ratio decreasing; the specific heat ratio has the same effect on the global strength of the viscous stress tensor. Secondly, the moment when the total temperature gradient in y direction deviates from the fixed value can be regarded as a physical criterion for judging the formation of the vortex structure. Thirdly, under the competition between the temperature gradients and the contact area of the two fluids, the average intensity of the non-equilibrium quantity related to the heat flux shows diversity, and the influence of the specific heat ratio is also quite remarkable.     

Modeling of cascaded high isolation bidirectional amplification in long-distance fiber-optic time and frequency synchronization system

We propose a physical model of estimating noise and asymmetry brought by high isolation Bi-directional erbiumdoped fiber amplifiers(Bi-EDFAs), no spontaneous lasing even with high gain, in longdistance fiber-optic time and frequency(T/F) synchronization system. It is found that the Rayleigh scattering noise can be suppressed due to the high isolation design, but the amplified spontaneous emission(ASE) noise generated by the high isolation Bi-EDFA and the bidirectional asymmetry of the transmission link caused by the high isolation Bi-EDFA will deteriorate the stability of the system.The calculated results show that under the influence of ASE noise, the frequency instability of a 1200 km system composed of 15 high isolation Bi-EDFAs is 1.773 × 10-13/1 s. And the instability caused by asymmetry is 2.6064 × 10-16/30000–35000 s if the total asymmetric length of the bidirectional link length is 30 m. The intensity noises originating from the laser and detector, the transfer delay fluctuations caused by the variation in ambient temperature and the jitter in laser output wavelength are also studied. The experiment composed of three high isolation Bi-EDFAs is done to confirm the theoretical analysis. In summary, the paper shows that the short-term instability of the T/F synchronization system composed of high isolation Bi-EDFAs is limited by the accumulation of ASE noise of amplifiers and the laser frequency drift, while the long-term instability is limited by the periodic variation in ambient temperature and the asymmetry of the amplifiers.The research results are useful for pointing out the direction to improve the stability of the fiber-optic T/F synchronization system.     

Study of Coulomb collisional effects on the propagation of electromagnetic waves in a turbulent plasma

The presence of Weibel instability in laser-irradiated fuel could be detrimental to the process of ablative implosion, which is necessary for achieving thermonuclear fusion reactions. In this paper, the effect of the Coulomb collisional within the turbulent plasma on the Weibel instability growth rate has been investigated for linear and circular polarization. The results indicate that the Weibel instability growth rate at circular polarization near the ignition centre of the fuel fusion (collisional plasma) is about 10⁵ times higher than the collisional Weibel instability growth rate at linear polarization. The Weibel instability growth rate is observed near the critical density of the fuel fusion (collisionless plasma) at linear polarization and enhancement near the foot of the heat in front of the fuel fusion. By increasing the steps of the density gradient plasma in the low-density corona, electromagnetic instability occurs at a higher stress flow. Therefore, the deposition condition of electron beam energy in circular polarization of turbulent plasma can be shifted to the fuel core for suitable ignition.     

周期调制结构平面薄膜电爆炸实验研究

在金属层表面引入微结构以实现对Z箍缩等离子体形成和发展过程中不稳定性的调控具有重要研究价值.在“强光一号”装置上(峰值电流~1.4 MA,上升时间~100 ns),开展了针对具有一维周期性凹槽调制结构的金属薄膜的电爆炸实验研究.实验负载采用外推型平面结构,基底为30μm厚铝膜,刻蚀周期为2 mm,刻蚀深度约为10μm.通过激光阴影成像、激光干涉成像和可见光自辐射成像等系统进行联合诊断.实验结果表明刻蚀结构对等离子体发展过程的不稳定性特征产生了明显调制作用,原本征波长也受到抑制,微结构周期对不稳定结构波长产生趋同效应;未刻蚀一侧边界层同样受刻蚀层结构的影响,在不稳定结构上表现出相似形貌,且内外侧不稳定性特征的耦合关联性增强;刻蚀凹槽处在爆炸过程中膨胀更为迅速,形成的表面等离子体结构与初始结构反相;在刻蚀结构的几何突变处会形成细长的等离子体喷流,在二分之一刻蚀波长处出现波谱特征峰.理论分析表明电流密度调制造成电热不稳定性分布改变是调控作用产生的重要原因.     

Mechanism of wave formation at the interface in explosive welding

Some of the main progress on the investigation of the mechanism of the wave formation in explosive welding at the Institute of Mechanics is summarized and others' previous works are reviewed. Our systematic experiments and analysis do not substantiate the theory of wave formation based on Karman vortex-street analogy or Helmholtz instability. On the contrary, they show that material strength insensitive to strain rate plays an important role. A simple hydro-plastic model is presented to explain the main features regarding the interfacial wave formation and to estimate the magnitude of wave length. The result is in broad agreement with experiment.     

Characteristic Regimes of Transitional Separation Bubbles in Unsteady Flow

This experimental investigation deals with transition phenomena of a separated boundary layer under unsteady inlet flow conditions. The main purpose of this investigation is to understand the influence of the rotor-stator interaction in turbomachinery on the subsequent, highly loaded boundary layer. The research project is divided into two phases. In the first phase, which has been completed recently, only the variation of mean velocity caused by upstream blades was simulated in the experiments while the free-stream turbulence intensity was retained at a constant low level. The experiments are carried out in an Eifel-type wind tunnel to investigate the laminar separated boundary layer of a flat plate under oscillating inlet conditions. The adverse pressure gradient, similar to that of turbomachines, is generated by the contoured upper wall. The unsteadiness is produced by a rotating flap located downstream of the test section. The reduced frequency, the amplitude and the mean Reynolds number are varied to simulate the conditions prevailing in turbomachines. In addition to the Kelvin–Helmholtz instability of the separated shear layer, a lower frequency instability was observed. This is frequently referred to as `free shear layer flapping' and results in two distinctly different ways of re-attachment, depending primarily on the Reynolds number. For low momentum thickness Reynolds numbers at the separation point, large-scale vortices locked to the frequency of the unsteady main flow are identified. They originate nearly at the top of the separation bubble and are ejected downstream. A fully turbulent boundary layer develops after these vortices mix out. For higher Reynolds numbers, transition is completed within a short length of the free shear layer and there-attachment region. The characteristic momentum thickness Reynolds number separating these two regimes in unsteady flow is about 125. The Strouhal number (reduced frequency) does not appear to have any significant effect. Based on the experimental results, this behaviour is discussed in some detail. This revised version was published online in July 2006 with corrections to the Cover Date.     

Cellular structures in solid fuel combustion

The objective of this contribution is to investigate whether the mechanism of the thermal diffusion instability in gaseous flames causing cellular flame structures also occurs during the combustion of porous solid fuel. Based on conservation for mass and energy, the relevant set of differential equations was derived. Assuming thermal equilibrium between fuel and oxidiser, a global energy equation was valid for both solid and gaseous phase. The resulting set of differential equations was discretised by the Collocation method to arrive at a system of algebraic equations. In order to investigate into cellular flame structures, an infinitesimal disturbance was superimposed onto the plane conversion front. Carrying out a linear instability analysis, yielded eigenvalues dependent on the wave number of the disturbance. A critical wave number exists below which the real part of the eigenvalues is positive, thus, indicating a regime of instability. Within this region, eigenvalues with a not-vanishing imaginary part of the eigen value existed causing cellular flame structures. However, the growth rate of disturbances was found to be small, which may explain the difficulty to investigate this phenomena experimentally.     

Purely elastic interfacial instabilities in superposed flow of polymeric fluids

Purely elastic interfacial stability of superposed plane Poiseuille flow of polymeric liquids has been investigated utilizing both asymptotic and numerical techniques. It is shown that these instabilities are caused by an unfavorable jump in the first normal stress difference across the fluid interface. To determine the significance of these instabilities in finite experimental geometries, a comparison between the maximum growth rates of purely elastic instabilities with instabilities driven primarily by a viscosity or a combined viscosity and elasticity difference is made. Based on this comparison, it is shown that purely elastic interfacial instabilities can play a major role in superposed flow of polymeric liquids in finite experimental geometries.     

Einfluß nichtnewtonscher Stoffeigenschaften auf die Taylor-Wirbelströmung

Zusammenfassung Nach einer kurzen Beschreibung über die Anwendung der Taylor-Wirbelströmung bei der Filtration von nichtnewtonschen Flüssigkeiten werden polymere Modellflüssigkeiten charakterisiert, die als Partikel- und Netzwerklösungen in einer Wirbelströmungsapparatur eingesetzt wurden. Die Fließkurven der Polymerlösungen zeigen aufgrund der extrem hochmolekularen Polymerproben neben strukturviskosen Erscheinungen auch dilatantes Fließverhalten mit besonderer Wirkung auf die Taylor-Wirbelströmung.Die Versuchsflüssigkeiten offenbaren vier verschiedene Typen von Strömungsinstabilitäten: spiralförmige, schwingende und stationäre Instabilitäten sowie gedämpfte Turbulenz. Während die stark strukturviskosen Netzwerklösungen alle genannten Formen aufweisen, fehlt bei den Partikellösungen die spiralförmige Instabilität.Unter Zuhilfenahme des Hantelmolekülmodells zur Beschreibung viskoelastischer Strömungsphänomene gelingt es, durch Einführung einer kritischen Deborahzahl den Einsatzpunkt nichtnewtonscher Taylorströmungseffekte vorauszusagen. Die gefundene Beziehung steht in engem Zusammenhang mit dehnviskositätserhöhenden Polymerwirkungen in Porenströmungen und mit reibungsmindernden Polymereffekten in turbulenten Rohrströmungen.

---

Experimental investigations dealing with dilute polymer solutions are described after a short review of the application of Taylor-vortex flow in the filtration processes of non-Newtonian fluids. The test fluids represent both viscoelastic solutions with isolated macromolecules and network solutions with power law fluid behaviour.These solutions show four different types of flow instabilities: spiral-shaped, oscillatory, steady and turbulent phenomena. The Taylor-numbers which depend upon the polymer concentration are determined for the onset of these instability types. For isolated macromolecule solutions, the Deborah-number concept for dilute dumbbell solutions can be applied to describe the first appearance of irregular nonstationary Taylor vortices.The present data are compared to literature values. This fluid behaviour is related to extensional viscosity increases which are also observed in porous media flow and turbulent pipe flow of dilute macromolecular solutions.

---

Nomenklatur A, B Konstanten aus Randbedingungen der Ringspaltströmung - C Polymerkonzentration - D Schergeschwindigkeit - De Deborahzahl - l Länge einer Wirbelzelle - L Zylinderlänge - m Gesamtanzahl der Wirbelpaare zur Bestimmung der Wellenzahl - M Molmasse der Polymere - M _w Gewichtsmittel der Molmasse - M _v Viskositätsmittel der Molmasse - n Drehzahl des Rotors - universelle Gaskonstante - r Radius - R _a Radius des Außenzylinders - R _i Radius des Innenzylinders - s Spaltweite - s ^* dimensionslose Spaltweite - T Temperatur - Ta Taylorzahl - v Umfangsgeschwindigkeit - z Anzahl der Wirbelpaare zur Bestimmung der Wellenzahl Griechische Symbole Deformationskoeffizient von Makromolekülen - Wellenzahl - Dehnrate - dynamische Viskosität - [] Grenzviskositätszahl - Relaxationszeit - Dichte - Schubspannung - Winkelgeschwindigkeit - _a Winkelgeschwindigkeit des Außenzylinders - _i Winkelgeschwindigkeit des Innenzylinders Indices c kritisch (erstmaliges Auftreten von Taylorwirbel) - N newtonsch - o onset, Schwellwert - P polymer - r radial - Sch schwingend - Spir spiralförmig - Stat stationär - Turb turbulent - T Taylorströmung - Umfangsrichtung Herrn Prof. Dr. Heinz Harnisch zum 60. Geburtstag gewidmet     

A nonlinear PSE method for two-fluid shear flows with complex interfacial topology

A nonlinear stability method is developed for laminar two-fluid shear flows which undergo changes in the interface topology. The method is based on the nonlinear parabolized stability equations (PSE) and incorporates a scalar-based interface capturing (IC) scheme in order to track complex deformations of the fluid interface. In doing so, the formulation retains the flexibility and physical insight of instability-wave based methods, while providing hydrodynamic modeling capabilities similar to direct numerical calculations: the new formulation, referred to as the IC-PSE, can capture the nonlinear physical mechanisms responsible for generating large-scale, two-fluid structures, without incurring heavy computational costs. This approach is valid for spatially developing, laminar two-fluid shear flows which are convectively unstable, and can naturally account for the growth of finite amplitude interfacial waves, along with changes to the interfacial topology. We demonstrate the accuracy of the IC-PSE against direct Navier–Stokes calculations for two-fluid mixing layers with density and viscosity stratification. The comparisons show that the IC-PSE can predict the dynamics of the instability waves and capture the formation of Kelvin–Helmholtz vortex rolls and large scale liquid structures, at an order of magnitude less computational cost than direct calculations. The role of surface tension in the IC-PSE formulation is shown to be valid for flows in which Re/We ? 1, and the method accurately predicts the formation and non-linear evolution of flow structures in this limit. This is demonstrated for spatially developing mixing layers which lead to vortex roll-up and ligaments, prior to droplet formation. The pinch-off process itself is a high surface tension phenomenon and in not considered herein. The method also accurately captures the effect of interfacial waves on the mean flow, and the topology changes during the non-linear evolution of the two-fluid structures.