Cold gas in cluster cores: global stability analysis and non-linear simulations of thermal instability

We perform global linear stability analysis and idealized numerical simulations in global thermal balance to understand the condensation of cold gas from hot/virial atmospheres (coronae), in particular the intracluster medium (ICM). We pay particular attention to geometry (e.g. spherical versus plane-parallel) and the nature of the gravitational potential. Global linear analysis gives a similar value for the fastest growing thermal instability modes in spherical and Cartesian geometries. Simulations and observations suggest that cooling in haloes critically depends on the ratio of the cooling time to the free-fall time (t(cool)/t(ff)). Extended cold gas condenses out of the ICM only if this ratio is smaller than a threshold value close to 10. Previous works highlighted the difference between the nature of cold gas condensation in spherical and plane-parallel atmospheres; namely, cold gas condensation appeared easier in spherical atmospheres. This apparent difference due to geometry arises because the previous plane-parallel simulations focused on in situ condensation of multiphase gas but spherical simulations studied condensation anywhere in the box. Unlike previous claims, our non-linear simulations show that there are only minor differences in cold gas condensation, either in situ or anywhere, for different geometries. The amount of cold gas depends on the shape of tcool/tff; gas has more time to condense if gravitational acceleration decreases towards the centre. In our idealized plane-parallel simulations with heating balancing cooling in each layer, there can be significant mass/energy/momentum transfer across layers that can trigger condensation and drive tcool/tff far beyond the critical value close to 10.

圆柱受限空间内喷雾燃烧振荡特性研究=Investigation of the Characteristics of Thermal Instability in A Confined Combustor

本文研究了圆柱受限空间内的喷雾火焰嫩烧压力振荡的特性。为了更清楚地了解火焰的构造, 首先测量了火焰的温度场，在较大的一次风和二次风变化范围内, 测量了压力的振荡特性。结果表明，火焰的稳定是由回流区完成的, 在较小的一次风燃料当量比和中等的二次风量时, 振荡最强, 达到100Pa左右的量级,其频率为200-230Hz左右,分析表明燃烧室中的振荡是轴向驻波振荡。

Numerical Investigation of Richtmyer-Meshkov Instability Driven by Cylindrical Shocks

In this paper, a numerical method with high order accuracy and high resolution was developed to simulate the Richtmyer-Meshkov(RM) instability driven by cylindrical shock waves. Compressible Euler equations in cylindrical coordinate were adopted for the cylindrical geometry and a third order accurate group control scheme was adopted to discretize the equations. Moreover, an adaptive grid technique was developed to refine the grid near the moving interface to improve the resolution of numerical solutions. The results of simulation exhibited the evolution process of RM instability, and the effect of Atwood number was studied. The larger the absolute value of Atwood number, the larger the perturbation amplitude. The nonlinear effect manifests more evidently in cylindrical geometry. The shock reflected from the pole center accelerates the interface for the second time, considerably complicating the interface evolution process, and such phenomena of reshock and secondary shock were studied.

The thermo-visco-plastic instability analysis of saturated soil

A theoretical analysis of instability of saturated soil is presented considering the simple shearing of a heat conducting thermo-visco-plastic material. It is shown that the instability is mainly the consequence of thermal softening which overcomes the strain hardening and the other type of instability is controlled by strain softening. The effects of other factors such as permeability to the instability are discussed in this paper.

On the shear instability of saturated soil

In this paper. the dynamic instability of simple shear of saturated soil is discussed. The governing equations are obtained based on mixture theory in which the inertia effect and the compressibility of grains are considered. Perturbation method is used to analyze and it is shown that two types of instability may exist. One of them is dominated by pore-pressure-softening, while the other by strain-softening.

Pull-in Instability Study of Carbon Nanotube Tweezers Under the Influence of Van der Waals Forces

The pull-in instability of two nanotubes under van der Waals force is studied. The cantilever beam with large deformation model is used. The influence of nanotube parameters such as the interior radius, the gap distance between the two nanotubes, etc, on the pull-in instability is studied. The critical nanotube length is determined for each specific set of nanotube parameters. The Galerkin method is applied to discretize the governing equations, and it shows good convergence.

Large-Eddy-Simulation of 3-Dimensional Rayleigh-Taylor Instability in Incompressible Fluids

The 3-dimensiqnal incompressible Rayleigh-Taylor instability is numerically studied through the large-eddy-simulation (LES) approach based on the passive scalar transport model. Both the instantaneous velocity and the passive scalar fields excited by sinu

Instability Analysis of Nonlinear Surface Waves in a Circular Cylindrical Container Subjected to a Vertical Excitation

Singular perturbation theory of two-time scale expansions was developed both in inviscid and weak viscous fluids to investigate the motion of single surface standing wave in a liquid-filled circular cylindrical vessel, which is subject to a vertical periodical oscillation. Firstly, it is assumed that the fluid in the circular cylindrical vessel is inviscid, incompressible and the motion is irrotational, a nonlinear evolution equation of slowly varying complex amplitude, which incorporates cubic nonlinear term, external excitation and the influence of surface tension, was derived from solvability condition of high-order approximation. It shows that when forced frequency is low, the effect of surface tension on mode selection of surface wave is not important. However, when forced frequency is high, the influence of surface tension is significant, and can not be neglected. This proved that the surface tension has the function, which causes free surface returning to equilibrium location. Theoretical results much close to experimental results when the surface tension is considered. In fact, the damping will appear in actual physical system due to dissipation of viscosity of fluid. Based upon weakly viscous fluids assumption, the fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates damping term and external excitation, was derived from linearized Navier-Stokes equation. The analytical expression of damping coefficient was determined and the relation between damping and other related parameters (such as viscosity, forced amplitude and depth of fluid) was presented. The nonlinear amplitude equation and a dispersion, which had been derived from the inviscid fluid approximation, were modified by adding linear damping. It was found that the modified results much reasonably close to experimental results. Moreover, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent. Finally, instability of the surface wave is analyzed and properties of the solutions of the modified amplitude equation are determined together with phase-plane trajectories. A necessary condition of forming stable surface wave is obtained and unstable regions are illustrated. (c) 2005 Elsevier SAS. All rights reserved.

Instability of Two-Layer Rayleigh-Benard Convection with Interfacial Thermocapillary Effect

The linear instability analysis of the Rayleigh-Allarangoni-Benard convection in a two-layer system of silicon oil 10cS and fluorinert FC70 liquids are performed in a larger range of two-layer depth ratios H, from 0.2 to 5.0 for different total depth H less than or equal to 12 mm. Our results are different from the previous study on the Rayleigh-Benard instability and show strong effects of thermocapillary force at the interface on the time-dependent oscillations arising from the onset of instability convection.

Effect of liquid bridge volume on the instability in small-Prandtl-number half zones

A perturbation method is used to examine the linear instability of thermocapillary convection in a liquid bridge of floating half-zone filled with a small Prandtl number fluid. The influence of liquid bridge volume on critical Marangoni number and flow features is analyzed. The neutral modes show that the instability is mainly caused by the bulk flow that is driven by the nonuniform thermocapillary forces acting on the free surface. The hydrodynamic instability is dominant in the case of small Prandtl number fluid and the first instability mode is a stationary bifurcation. The azimuthal wave number for the most dangerous mode depends on the liquid bridge volume, and is not always two as in the case of a cylindrical liquid bridge with aspect ratio near 0.6. Its value may be equal to unity when the liquid bridge is relatively slender.

Numerical and Analytical Study on the Pull-in Instability of Micro-Structure under Electrostatic Loading

The one-mode analysis method on the pull-in instability of micro-structure under electrostatic loading is presented. Taylor series are used to expand the electrostatic loading term in the one-mode analysis method, which makes analytical solution available. The one-mode analysis is the combination of Galerkin method and Cardan solution of cubic equation. The one-mode analysis offers a direct computation method on the pull-in voltage and displacement. In low axial loading range, it shows little difference with the established multi-mode analysis on predicting the pull-in voltages for three different structures (cantilever, clamped-clamped beams and the plate with four edges simply-supported) studied here. For numerical multi-mode analysis, we also show that using the structural symmetry to select the symmetric mode can greatly reduce both the computation effort and the numerical fluctuation.

Marangoni-Benard instability with the exchange of evaporation at liquid-vapour interface

A new two-sided model rather than the one-sided model in previous works is put forward. The linear instability analysis is performed on the Marangoni-Benard convection in the two-layer system with an evaporation interface. We define a new evaporation Biot number which is different from that in the one-sided model, and obtain the curves of critical Marangoni number versus wavenumber. The influence of evaporation velocity and Biot number on the system is discussed and a new phenomenon uninterpreted before is now explained from our numerical results.