Mitigating instability by actuating the swirler in a combustor

Thermoacoustic instability in a lean premixed combustor is a major impediment towards reliable operation of gas turbine engines for both aerospace and land based applications. In this communication, we investigate the following concept: in a laboratory combustor, could the otherwise static swirler be actuated to a rotary motion, such that the higher intensity turbulence and higher swirl number generated in the flame stabilization region might alter the flame position, structure and thereby assist in mitigating thermoacoustic instabilities? Results obtained using microphone and high speed imaging, show prominent reductions in the amplitudes of the first mode of the thermoacoustically unstable flame, with increased rotation rate of the swirler. (C) 2015 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Structure-rheology relationship in a sheared lamellar fluid

The structure-rheology relationship in the shear alignment of a lamellar fluid is studied using a mesoscale model which provides access to the lamellar configurations and the rheology. Based on the equations and free energy functional, the complete set of dimensionless groups that characterize the system are the Reynolds number (rho gamma L-2/mu), the Schmidt number (mu/rho D), the Ericksen number (mu(gamma)/B), the interface sharpness parameter r, the ratio of the viscosities of the hydrophilic and hydrophobic parts mu(r), and the ratio of the system size and layer spacing (L/lambda). Here, rho and mu are the fluid density and average viscosity, (gamma) over dot is the applied strain rate, D is the coefficient of diffusion, B is the compression modulus, mu(r) is the maximum difference in the viscosity of the hydrophilic and hydrophobic parts divided by the average viscosity, and L is the system size in the cross-stream direction. The lattice Boltzmann method is used to solve the concentration and momentum equations for a two dimensional system of moderate size (L/lambda = 32) and for a low Reynolds number, and the other parameters are systematically varied to examine the qualitative features of the structure and viscosity evolution in different regimes. At low Schmidt numbers where mass diffusion is faster than momentum diffusion, there is fast local formation of randomly aligned domains with ``grain boundaries,'' which are rotated by the shear flow to align along the extensional axis as time increases. This configuration offers a high resistance to flow, and the layers do not align in the flow direction even after 1000 strain units, resulting in a viscosity higher than that for an aligned lamellar phase. At high Schmidt numbers where momentum diffusion is fast, the shear flow disrupts layers before they are fully formed by diffusion, and alignment takes place by the breakage and reformation of layers by shear, resulting in defects (edge dislocations) embedded in a background of nearly aligned layers. At high Ericksen number where the viscous forces are large compared to the restoring forces due to layer compression and bending, shear tends to homogenize the concentration field, and the viscosity decreases significantly. At very high Ericksen number, shear even disrupts the layering of the lamellar phase. At low Ericksen number, shear results in the formation of well aligned layers with edge dislocations. However, these edge dislocations take a long time to anneal; the relatively small misalignment due to the defects results in a large increase in viscosity due to high layer stiffness and due to shear localization, because the layers between defects get pinned and move as a plug with no shear. An increase in the viscosity contrast between the hydrophilic and hydrophobic parts does not alter the structural characteristics during alignment. However, there is a significant increase in the viscosity, due to pinning of the layers between defects, which results in a plug flow between defects and a localization of the shear to a part of the domain.

Numerical study of driven flows of shear thinning viscoelastic fluids in rectangular cavities

In this work, we present a numerical study of flow of shear thinning viscoelastic fluids in rectangular lid driven cavities for a wide range of aspect ratios (depth to width ratio) varying from 1/16 to 4. In particular, the effect of elasticity, inertia, model parameters and polymer concentration on flow features in rectangular driven cavity has been studied for two shear thinning viscoelastic fluids, namely, Giesekus and linear PTT. We perform numerical simulations using the symmetric square root representation of the conformation tensor to stabilize the numerical scheme against the high Weissenberg number problem. The variation in flow structures associated with merging and splitting of elongated vortices in shallow cavities and coalescence of corner eddies to yield a second primary vortex in deep cavities with respect to the variation in flow parameters is discussed. We discuss the effect of the dominant eigenvalues and the corresponding eigenvectors on the location of the primary eddy in the cavity. We also demonstrate, by performing numerical simulations for shallow and deep cavities, that where the Deborah number (based on convective time scale) characterizes the elastic behaviour of the fluid in deep cavities, Weissenberg number (based on shear rate) should be used for shallow cavities. (C) 2016 Elsevier B.V. All rights reserved.

A nano-scale instability in the beta phase of dilute Ti-Mo alloys

A nano-scale instability in the beta phase resulting in the formation of the disordered orthorhombic O' phase has been discovered in a fairly dilute binary Ti-Mo alloy, using selected area electron diffraction and high resolution scanning transmission electron microscopy. The O' phase informed in the alloy when the Mo content exceeds a critical value. The instability occurs in beta-solutionized samples that have been quenched to room temperature and is found to co-exist with athermal omega to phase. Interestingly, this nano-scale instability, involving the {110}<1<(1)over bar>0> soft-phonon shuffle, occurs in the beta phase without deliberate additions of either interstitial or substitutional solutes. (C) 2016 Elsevier Ltd. All rights reserved.

Effect of particle size of sand and surface asperities of reinforcement on their interface shear behaviour

This paper investigates the effect of particle size of sand and the surface asperities of reinforcing material on their interlocking mechanism and its influence on the interfacial shear strength under direct sliding condition. Three sands of different sizes with similar morphological characteristics and four different types of reinforcing materials with different surface features were used in this study. Interface direct shear tests on these materials were performed in a specially developed symmetric loading interface direct shear test setup. Morphological characteristics of sand particles were determined from digital image analysis and the surface roughness of the reinforcing materials was measured using an analytical expression developed for this purpose. Interface direct shear tests at three different normal stresses were carried out by shearing the sand on the reinforcing material fixed to a smooth surface. Test results revealed that the peak interfacial friction and dilation angles are hugely dependent upon the interlocking between the sand particles and the asperities of reinforcing material, which in turn depends on the relative size of sand particles and asperities. Asperity ratio (AS/D-50) of interlocking materials, which is defined as the ratio of asperity spacing of the reinforcing material and the mean particle size of sand was found to govern the interfacial shear strength with highest interfacial strength measured when the asperity ratio was equal to one, which represents the closest fitting of sand particles into the asperities. It was also understood that the surface roughness of the reinforcing material influences the shear strength to an extent, the influence being more pronounced in coarser particles. Shear bands in the interface shear tests were analysed through image segmentation technique and it was observed that the ratio of shear band thickness (t) to the median particle size (D-50) was maximum when the AS/D-50 was equal to one. (C) 2015 Elsevier Ltd. All rights reserved.

Morphological transitions and buckling characteristics in a nanoparticle-laden sessile droplet resting on a heated hydrophobic substrate

In this work, we have established the evaporation-liquid flow coupling mechanism by which sessile nanofluid droplets on a hydrophobic substrate evaporate and agglomerate to form unique morphological features under controlled external heating. It is well understood that evaporation coupled with internal liquid flow controls particle transport in a spatiotemporal sense. Flow characteristics inside the heated droplet are investigated and found to be driven by the buoyancy effects. Velocity magnitudes are observed to increase by an order at higher temperatures with similar looking flow profiles. The recirculating flow induced particle transport coupled with collision of particles and shear interaction between them leads to the formation of dome shaped viscoelastic shells of different dimensions depending on the surface temperature. These shells undergo sol-gel transition and subsequently undergo buckling instability leading to the formation of daughter cavities. With an increase in the surface temperature, droplets exhibit buckling from multiple sites over a larger sector in the top half of the droplet. Irrespective of the initial nanoparticle concentration and substrate temperature, growth of a daughter cavity (subsequent to buckling) inside the droplet is found to be controlled by the solvent evaporation rate from the droplet periphery and is shown to exhibit a universal trend.

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.

Maintaining local displacive disorders in Na0.5Bi0.5TiO3 piezoceramics by K0.5Bi0.5TiO3 substitution

A comprehensive global and local structural characterization using X-ray and neutron diffraction, and EXAFS have been performed to investigate the effect of electrical poling on the structure of (1-x)Na0.5Bi0.5TiO3-(x)K0.5Bi0.5TiO3 compositions. While the annealed samples showed definite indication of presence of local displacive disorders in all the compositions studied (0 local displacive disorder decreased with increasing K0.5Bi0.5TiO3 content in the samples. (C) 2016 Elsevier Ltd. All rights reserved.

Deformation, Phase Transformation and Recrystallization in the Shear Bands Induced by High-Strain Rate Loading in Titanium and its Alloys

alpha-titanium and its alloys with a dual-phase structure (alpha+beta) were deformed dynamically under strain rate of about 10(4) s(-1). The formation and microstructural evolution of the localized shear bands were characterized by scanning electron microscopy (SEM) and transmission electron microscopy (TEM). The results reveal that both the strain and strain rate should be considered simultaneously as the mechanical conditions for shear band formation, and twinning is an important mode of deformation. Both experimental and calculation show that the materials within the bands underwent a superhigh strain rate (9 x 10(5) s(-1)) deformation, which is two magnitudes of that of average strain rate required for shear band formation; the dislocations in the bands can be constricted and developed into cell structures; the phase transformation from alpha to alpha(2) within the bands was observed, and the transformation products (alpha(2)) had a certain crystallographic orientation relationship with their parent; the equiaxed grains with an average size of 10 mu m in diameter observed within the bands are proposed to be the results of recrystallization.

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

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

Stress-strain fields and the effectiveness shear properties for three-phase composites with imperfect interface

Based on the 'average stress in the matrix' concept of Mori and Tanaka (:Mori, T., Tanaka, K., 1973. Average stress in matrix and average elastic energy of materials with misfitting inclusion. Acta Metall. 21, 571-580) a micromechanical model is presented for the prediction of the elastic fields in coated inclusion composites with imperfect interfaces. The solutions of the effective elastic moduli for this kind of composite are also obtained. In two kinds of composites with coated particulates and fibers, respectively, the interface imperfections are takes to the assumption that the interface displacement discontinues are linearly related to interface tractions like a spring layer of vanishing thickness. The resulting effective shear modulus for each material and the stress fields in the composite are presented under a transverse shear loading situation.

Influence of Local Magnetic Fields on P-Doped Si Floating Zone Melting Crystal Growth in Microgravity

A two-dimensional axisymmetric numerical model is presented to study the influence of local magnetic fields on P-doped Si floating zone melting crystal growth in microgravity. The model is developed based on the finite difference method in a boundary-fitted curvilinear coordinate system. Extensive numerical simulations are carried out, and parameters studied include the curved growth interface shape and the magnetic field configurations. Computed results show that the local magnetic field is more effective in reducing the impurity concentration nonuniformity at the growth interface in comparison with the longitudinal magnetic field. Moreover, the curved growth interface causes more serious impurity concentration nonuniformity at the growth interface than the case with a planar growth interface.

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.

An exact solution for the three-phase piezoelectric cylinder model under antiplane shear and its applications to piezoelectric composites

A three-phase piezoelectric cylinder model is proposed and an exact solution is obtained for the model under a farfield antiplane mechanical load and a far-field inplane electrical load. The three-phase model can serve as a fiber/interphase layer/matrix model, in terms of which a lot of interesting mechanical and electrical coupling phenomena induced by the interphase layer are revealed. It is found that much more serious stress and electrical field concentrations occur in the model with the interphase layer than those without any interphase layer. The three-phase model can also serve as a fiber/matrix/composite model, in terms of which a generalized self-consistent approach is developed for predicting the effective electroelastic moduli of piezoelectric composites. Numerical examples are given and discussed in detail.

Material response and failure mechanism of unidirectional metal matrix composites under impulsive shear loading

The material response and failure mechanism of unidirectional metal matrix composite under impulsive shear loading are investigated in this paper. Both experimental and analytical studies were performed. The shear strength of unidirectional C-f/A356.0 composite and A356.0 aluminum alloy at high strain rate were measured with a modified split Hopkinson torsional bar technique. The results indicated that the carbon fibers did not improve the shear strength of aluminum matrix if the fiber orientation aligned with the shear loading axis. The microscopic inspection of the fractured surface showed a multi-scale zigzag feature which implied a complicated shear failure mechanism in the composite. In addition to testing, the micromechanical stress field in the composite was analyzed by the generalized Eshelby equivalent method (GEEM). The influence of cracking in matrix on the micromechanical stress field was investigated as well. The results showed that the stress distribution in the composite is quite nonhomogeneous and very high shear stress concentrations are found in some regions in the matrix. The high shear stress concentration in the matrix induces tensile cracking at 45 degrees to the shear direction. This in turn aggravates the stress concentration at the fiber/matrix interface and finally leads to a catastrophic failure in the composite. From the correlation between the analysis and experimental results, the shear failure mechanism of unidirectional C-f/A356.0 composite can be elucidated qualitatively.