The problem of clear air turbulence: Changing perspectives in the understanding of the phenomenon

Due to rapid improvements in on-board instrumentation and atmospheric observation systems, in most cases, aircraft are able to steer clear of regions of adverse weather. However, they still encounter unexpected bumpy flight conditions in regions away from storms and clouds. This is the phenomenon of clear air turbulence (CAT), which has been a challenge to our understanding as well as efforts at prediction. While most of such cases result in mild discomfort, a few cases can be violent leading to serious injuries to passengers and damage to the aircraft. The underlying physical mechanisms have been sought to be explained in terms of fluid dynamic instabilities and waves in the atmosphere. The main mechanisms which have been proposed are: (i) Kelvin-Helmholtz instability of shear layers, (ii) waves generated from flow over mountains, (iii) inertia-gravity waves from clouds and other sources, (iv) spontaneous imbalance theory and (v) horizontal vortex tubes. This has also undergone a change over the years. We present an overview of the mechanisms proposed and their implications for prediction.

Mixed finite elements for electromagnetic analysis

The occurrence of spurious solutions is a well-known limitation of the standard nodal finite element method when applied to electromagnetic problems. The two commonly used remedies that are used to address this problem are (i) The addition of a penalty term with the penalty factor based on the local dielectric constant, and which reduces to a Helmholtz form on homogeneous domains (regularized formulation); (ii) A formulation based on a vector and a scalar potential. Both these strategies have some shortcomings. The penalty method does not completely get rid of the spurious modes, and both methods are incapable of predicting singular eigenvalues in non-convex domains. Some non-zero spurious eigenvalues are also predicted by these methods on non-convex domains. In this work, we develop mixed finite element formulations which predict the eigenfrequencies (including their multiplicities) accurately, even for nonconvex domains. The main feature of the proposed mixed finite element formulation is that no ad-hoc terms are added to the formulation as in the penalty formulation, and the improvement is achieved purely by an appropriate choice of finite element spaces for the different variables. We show that the formulation works even for inhomogeneous domains where `double noding' is used to enforce the appropriate continuity requirements at an interface. For two-dimensional problems, the shape of the domain can be arbitrary, while for the three-dimensional ones, with our current formulation, only regular domains (which can be nonconvex) can be modeled. Since eigenfrequencies are modeled accurately, these elements also yield accurate results for driven problems. (C) 2014 Elsevier Ltd. All rights reserved.

A nodal domain theorem for integrable billiards in two dimensions

Eigenfunctions of integrable planar billiards are studied - in particular, the number of nodal domains, nu of the eigenfunctions with Dirichlet boundary conditions are considered. The billiards for which the time-independent Schrodinger equation (Helmholtz equation) is separable admit trivial expressions for the number of domains. Here, we discover that for all separable and nonseparable integrable billiards, nu satisfies certain difference equations. This has been possible because the eigenfunctions can be classified in families labelled by the same value of m mod kn, given a particular k, for a set of quantum numbers, m, n. Further, we observe that the patterns in a family are similar and the algebraic representation of the geometrical nodal patterns is found. Instances of this representation are explained in detail to understand the beauty of the patterns. This paper therefore presents a mathematical connection between integrable systems and difference equations. (C) 2014 Elsevier Inc. All rights reserved.

Acoustic response of vortex breakdown modes in a coaxial isothermal unconfined swirling jet

The present experimental work is concerned with the study of amplitude dependent acoustic response of an isothermal coaxial swirling jet. The excitation amplitude is increased in five distinct steps at the burner's Helmholtz resonator mode (i.e., 100 Hz). Two flow states are compared, namely, sub-critical and super-critical vortex breakdown (VB) that occur before and after the critical conical sheet breakdown, respectively. The geometric swirl number is varied in the range 2.14-4.03. Under the influence of external pulsing, global response characteristics are studied based on the topological changes observed in time-averaged 2D flow field. These are obtained from high resolution 2D PIV (particle image velocimetry) in the longitudinal-mid plane. PIV results also illustrate the changes in the normalized vortex core coordinates (r(vcc)/(r(vcc))(0) (Hz), y(vcc)/(y(vcc))(0) (Hz)) of internal recirculation zone (IRZ). A strong forced response is observed at 100 Hz (excitation frequency) in the convectively unstable region which get amplified based on the magnitude of external forcing. The radial extent of this forced response region at a given excitation amplitude is represented by the acoustic response region (b). The topological placement of the responsive convectively unstable region is a function of both the intensity of imparted swirl (characterized by swirl number) and forcing amplitude. It is observed that for sub-critical VB mode, an increase in the excitation amplitude till a critical value shifts the vortex core centre (particularly, the vortex core moves downstream and radially outwards) leading to drastic fanning-out/widening of the IRZ. This is accompanied by similar to 30% reduction in the recirculation velocity of the IRZ. It is also observed that b < R (R: radial distance from central axis to outer shear layer-OSL). At super-critical amplitudes, the sub-critical IRZ topology transits back (the vortex core retracts upstream and radially inwards) and finally undergoes a transverse shrinkage ((r(vcc))/(r(vcc))(0 Hz) decreases by similar to 20%) when b >= R. In contrast, the vortex core of super-critical breakdown mode consistently spreads radially outwards and is displaced further downstream. Finally, the IRZ fans-out at the threshold excitation amplitude. However, the acoustic response region b is still less than R. This is explained based on the characteristic geometric swirl number (S-G) of the flow regimes. The super-critical flow mode with higher S-G (hence, higher radial pressure drop due to rotational effect which scales as Delta P similar to rho u theta(2) and acts inwards towards the center line) compared to sub-critical state imposes a greater resistance to the radial outward spread of b. As a result, the acoustic energy supplied to the super-critical flow mode increases the degree of acoustic response at the pulsing frequency and energizes its harmonics (evident from power spectra). As a disturbance amplifier, the stronger convective instability mode within the flow structure of super-critical VB causes the topology to widen/fan-out severely at threshold excitation amplitude. (C) 2015 AIP Publishing LLC.

A model supersonic buried-nozzle jet: instability and acoustic wave scattering and the far-field sound

We consider sound source mechanisms involving the acoustic and instability modes of dual-stream isothermal supersonic jets with the inner nozzle buried within an outer shroud-like nozzle. A particular focus is scattering into radiating sound waves at the shroud lip. For such jets, several families of acoustically coupled instability waves exist, beyond the regular vortical Kelvin-Helmholtz mode, with different shapes and propagation characteristics, which can therefore affect the character of the radiated sound. In our model, the coaxial shear layers are vortex sheets while the incident acoustic disturbances are the propagating shroud modes. The Wiener-Hopf method is used to compute their scattering at the sharp shroud edge to obtain the far-field radiation. The resulting far-field directivity quantifies the acoustic efficiency of different mechanisms, which is particularly important in the upstream direction, where the results show that the scattered sound is more intense than that radiated directly by the shear-layer modes.

Structure and Dynamics of Octamethylcyclotetrasiloxane Confined between Mica Surfaces

Using a molecular model for octamethylcydotetrasiloxane (OMCTS), molecular dynamics simulations are carried out to probe the phase state of OMCTS confined between two mica surfaces in equilibrium With a reservoir. Molecular dynamics simulations are carried out for elevations ranging from 5 to 35 K above the melting point for the OMCTS model used in this study. The Helmholtz free energy is, computed for a specific confinement using the :two-phase thermodynamic (2PT) method. Analysis of the in-plane pair correlation functions did not reveal signatures of freezing even under an extreme confinement of two layers. OMCTS is found to orient with a wide distribution of orientations with respect to the mica surface, with a distinct preference for the surface parallel configuration in the contact layers. The self-intermediate scattering function is found to decay with increasing relaxation times as the surface separation is decreased, and the two-step relaxation in the scattering function, a signature of glassy dynamics, distinctly evolves as the temperature is lowered. However, even at 5 K above the melting point, we did not observe a freezing transition and the self-intermediate scattering functions relax within 200 ps for the seven-layered confined system. The self diffusivity and relaxation times obtained from the Kohlrausch-Williams-Watts stretched exponential fits to the late alpha-relaxation exhibit power law scalings with the packing fraction as predicted by mode coupling theory. A distinct discontinuity in the Helmholtz free energy, potential energy, and a sharp change in the local bond order parameter, Q(4), was observed at 230 K for a five-layered system upon cooling, indicative of a first-order transition. A freezing point depression of about 30 K was observed for this five-layered confined system, and at the lower temperatures, contact layers were found to be disordered with long-range order present only in the inner layers. These dynamical signatures indicate that confined OMCTS undergoes a slowdown akin to a fluid approaching a glass transition upon increasing confinement, and freezing under confinement would require substantial subcooling below the bulk melting point of OMCTS.

Evaluation of the Cumulative Formation of High-Speed Liquid Jets

This paper describes the generation of pulsed, high-speed liquid jets using the cumulation method. This work mainly includes (1) the design of the nozzle assembly, (2) the measurement of the jet velocity and (3) flow visualization of the injection sequences. The cumulation method can be briefly described as the liquid being accelerated first by the impact of a moving projectile and then further after it enters a converging section. The experimental results show that the cumulation method is useful in obtaining a liquid jet with high velocity. The flow visualization shows the roles of the Rayleigh-Taylor and Kelvin-Helmholtz instabilities in the breakup of the liquid depend on the jet diameter and the downstream distance. When the liquid jet front is far downstream from the nozzle exit, the jet is decelerated by air drag. Meanwhile, large coherent vortex structures are formed surrounding the jet. The liquid will break up totally by the action of these vortices. Experimental results showing the effect of the liquid volume on the jet velocity are also included in this paper. Finally, a method for measuring the jet velocity by cutting two carbon rods is examined.

纳米铁热力学性质的分子动力学计算

应用分子动力学方法计算了包括纳米铁晶体的熵、Helmholtz自由能等一系列热力学性质，发现纳米晶体的熵比单晶体的熵值高得多，且这些差别随温度增高而迅速加大，这些结果可以用纳米晶体的特殊结构来解释。纳米晶体中大量不规则原子的存在是熵值增大的根本原因。因而，只要改变纳米晶体的尺寸(事实上也是改变了其中不规则原子的数目)将会改变它的熵值及其它的热力学量，从微观结构与热力学量的关系来说，它将有助于改善纳米晶体的性能。

Orientation-based Continuum Damage Models for Rocks

A general formulation of the Helmholtz free energy used in thermodynamics of damage process of rocks is derived within a multi-scale framework. Such a physically-based thermodynamic state potential has a hybrid, discrete/continuum, nature in the sense tha

Evolution of three-dimensional coherent structures in compressible axisymmetric jet

A high order difference scheme is used to simulate the spatially developing compressible axisymmetric jet. The results show that the Kelvin-Helmholtz instability appears first when the jet loses its stability, and then with development of jet the increase in nonlinear effects leads to the secondary instability and the formation of the streamwise vortices. The evolution of the three-dimensional coherent structure is presented. The computed results verify that in axisymmetric jet the secondary instability and formation of the streamwise vortices are the important physical mechanism of enhancing the flow mixing and transition occurring.

VISUALIZATION RESULTS OF AXISYMMETRICAL WAKE FLOW WITH/WITHOUT ACOUSTIC EXCITATION

Visualization results demonstrate the evolution of Kelvin-Helmholtz unstable waves into vortex pairing in a separated shear layer of a blunf circular. The results with acoustic excitation are quite different from that without acoustic excitation, and the phenomenon with excitation in a separated shear layer follows the rule of Devil s staircase, which always occurs in a non-linear dynamical system of two coupling vibrators.

地震条件下损伤–愈合模型的初步研究

脆性介质的损伤与破坏是力学中极为复杂、极具挑战性的难题之一，而地震是一类典型脆性介质的突发破 裂或失稳。在地震等灾变发生前都会出现一定的前兆现象，地震之后地壳介质又会缓慢愈合，地壳介质的损伤愈 合过程对地震预测具有重要的研究意义。首先介绍V. Lyakhovsky 等提出的损伤流变模型，且结合岩石破坏声发射 实验数据分析该模型的特点；进一步基于损伤力学，从Helmholtz 自由能出发，并考虑地质材料的可愈合效应， 建立一维损伤–愈合模型，采用数值计算的方法得到地质材料的损伤–愈合模式。用常数加载率的线性载荷叠加 一个正弦扰动来模拟固体潮引力，模拟计算该加载条件下材料的损伤–愈合过程；同时也考察了整个损伤愈合演 化过程中加卸载响应比值的变化情况，发现在灾变发生之前加卸载响应比值会出现明显异常、回落的现象。再次 验证了加卸载响应比作为一种重要的前兆现象，能够定量地刻画介质的损伤演化过程，这也为加卸载响应比方法 提供了更为坚实的物理基础。

平面混合流拟序结构的直接数值模拟

采用高阶精度差分格式,求解二维可压缩N-S方程,直接数值模拟了可压缩平面混合流的二维拟序结构.给出了流动失稳,Kelvin-Helmholtz不稳定波的发展.展向大涡的卷起和相邻两涡卷对并,包括3次对并的发展过程.研究了平面混合流时-空的发展和可压缩效应对其发展的影响.

径向收缩的等离子体的运动不稳定性

<正> 近几年来,人们详细地研究过等离子体某些运动状态的不稳定性,如Rayleigh-Taylor不稳定性和Kelvin-Helmholtz不稳定性等等。本文先根据雪铲模型讨论收缩效应的运动不稳定性,其次指出,在径向运动情形下,稳定性不只和加速度有关,而且也和速度等其它因素有关,因此,它并不相当于平面情形的Rayleig-Taylor不稳定性。 1.收缩效应的运动不稳定性在雪铲模型里,如果假定在某时刻t_1出现径向、周向和轴向的微扰动γ,θ,z,并以初始时刻to的位置θo、Zo以及时刻t为自变量,则在假定(γ,θ,z)的形式为(γ(t),-iθ(t),-iz(t))eimθo+ikzσ之后,就可得到线性化的扰动方程(采用高斯单位)为:

Numerical simulation of coherent structure in two-dimensional compressible mixing layers

The coherent structure in two-dimensional mixing layers is simulated numerically with the compressible Navier-Stokes equations. The Navier-Stokes equations are discretized with high-order accurate upwind compact schemes. The process of development of flow structure is presented: loss of stability, development of Kelvin-Helmholtz instability, rolling up and pairing. The time and space development of the plane mixing layer and influence of the compressibility are investigated.