Instabilities in a fluid overlying an inclined anisotropic and inhomogeneous porous layer

In this paper, linear stability analysis on a Newtonian fluid film flowing under the effect of gravity over an inclined porous medium saturated with the same fluid in isothermal condition is carried out. The focus is placed on the effect of the anisotropic and inhomogeneous variations in the permeability of the porous medium on the shear mode and surface mode instabilities. The fluid-porous system is modelled by a coupled two-dimensional Navier-Stokes/Darcy problem. The perturbation equations are solved numerically using the Chebyshev collocation method. Detailed stability characteristics as a function of the depth ratio (the ratio of the depth of the fluid layer to that of the porous layer), the anisotropic parameter (the ratio of the permeability in the direction of the basic flow to that in the direction transverse to the basic flow) and the inhomogeneity functions are presented.

Absolute/Convective Instability Transition in a Backward Facing Step Combustor: Fundamental Mechanism and Influence of Density Gradient

Hydrodynamic instabilities of the flow field in lean premixed gas turbine combustors can generate velocity perturbations that wrinkle and distort the flame sheet over length scales that are smaller than the flame length. The resultant heat release oscillations can then potentially result in combustion instability. Thus, it is essential to understand the hydrodynamic instability characteristics of the combustor flow field in order to understand its overall influence on combustion instability characteristics. To this end, this paper elucidates the role of fluctuating vorticity production from a linear hydrodynamic stability analysis as the key mechanism promoting absolute/convective instability transitions in shear layers occurring in the flow behind a backward facing step. These results are obtained within the framework of an inviscid, incompressible, local temporal and spatio-temporal stability analysis. Vorticity fluctuations in this limit result from interaction between two competing mechanisms-(1) production from interaction between velocity perturbations and the base flow vorticity gradient and (2) baroclinic torque in the presence of base flow density gradients. This interaction has a significant effect on hydrodynamic instability characteristics when the base flow density and velocity gradients are colocated. Regions in the space of parameters characterizing the base flow velocity profile, i.e., shear layer thickness and ratio of forward to reverse flow velocity, corresponding to convective and absolute instability are identified. The implications of the present results on understanding prior experimental studies of combustion instability in backward facing step combustors and hydrodynamic instability in other flows such as heated jets and bluff body stabilized flames is discussed.

Wall-mode instability in plane shear flow of viscoelastic fluid over a deformable solid

The linear stability analysis of a plane Couette flow of an Oldroyd-B viscoelastic fluid past a flexible solid medium is carried out to investigate the role of polymer addition in the stability behavior. The system consists of a viscoelastic fluid layer of thickness R, density rho, viscosity eta, relaxation time lambda, and retardation time beta lambda flowing past a linear elastic solid medium of thickness HR, density rho, and shear modulus G. The emphasis is on the high-Reynolds-number wall-mode instability, which has recently been shown in experiments to destabilize the laminar flow of Newtonian fluids in soft-walled tubes and channels at a significantly lower Reynolds number than that for flows in rigid conduits. For Newtonian fluids, the linear stability studies have shown that the wall modes become unstable when flow Reynolds number exceeds a certain critical value Re c which scales as Sigma(3/4), where Reynolds number Re = rho VR/eta, V is the top-plate velocity, and dimensionless parameter Sigma = rho GR(2)/eta(2) characterizes the fluid-solid system. For high-Reynolds-number flow, the addition of polymer tends to decrease the critical Reynolds number in comparison to that for the Newtonian fluid, indicating a destabilizing role for fluid viscoelasticity. Numerical calculations show that the critical Reynolds number could be decreased by up to a factor of 10 by the addition of small amount of polymer. The critical Reynolds number follows the same scaling Re-c similar to Sigma(3/4) as the wall modes for a Newtonian fluid for very high Reynolds number. However, for moderate Reynolds number, there exists a narrow region in beta-H parametric space, corresponding to very dilute polymer solution (0.9 less than or similar to beta < 1) and thin solids (H less than or similar to 1.1), in which the addition of polymer tends to increase the critical Reynolds number in comparison to the Newtonian fluid. Thus, Reynolds number and polymer properties can be tailored to either increase or decrease the critical Reynolds number for unstable modes, thus providing an additional degree of control over the laminar-turbulent transition.

REDUCED ORDER MODELLING OF COMBUSTION INSTABILITY IN A BACKWARD FACING STEP COMBUSTOR

This paper develops a fully coupled time domain Reduced Order Modelling (ROM) approach to model unsteady combustion dynamics in a backward facing step combustor The acoustic field equations are projected onto the canonical acoustic eigenmodes of the systems to obtain a coupled system of modal evolution equations. The heat release response of the flame is modelled using the G-equation approach. Vortical velocity fluctuations that arise due to shear layer rollup downstream of the step are modelled using a simplified 1D-advection equation whose phase speed is determined from a linear, local, temporal stability analysis of the shear layer just downstream of the step. The hydrodynamic stability analysis reveals a abrupt change in the value of disturbance phase speed from unity for Re < Re-crit to 0.5 for Re > Re-crit, where Remit for the present geometry was found to be approximate to 10425. The results for self-excited flame response show highly wrinkled flame shapes that are qualitatively similar to those seen in prior experiments of acoustically forced flames. The effect of constructive and destructive interference between the two contributions to flame surface wrinkling results in high amplitude wrinkles for the case when K-c -> 1.

ABSOLUTE/CONVECTIVE INSTABILITY TRANSITION IN A BACKWARD FACING STEP COMBUSTOR: FUNDAMENTAL MECHANISM AND INFLUENCE OF DENSITY GRADIENT

Hydrodynamic instabilities of the flow field in lean premixed gas turbine combustors can generate velocity perturbations that wrinkle and distort the flame sheet over length scales that are smaller than the flame length. The resultant heat release oscillations can then potentially result in combustion instability. Thus, it is essential to understand the hydrodynamic instability characteristics of the combustor flow field in order to understand its overall influence on combustion instability characteristics. To this end, this paper elucidates the role of fluctuating vorticity production from a linear hydrodynamic stability analysis as the key mechanism promoting absolute/convective instability transitions in shear layers occurring in the flow behind a backward facing step. These results are obtained within the framework of an inviscid, incompressible, local temporal and spatio-temporal stability analysis. Vorticity fluctuations in this limit result from interaction between two competing mechanisms - (1) production from interaction between velocity perturbations and the base flow vorticity gradient and (2) baroclinic torque in the presence of base flow density gradients. This interaction has a significant effect on hydrodynamic instability characteristics when the base flow density and velocity gradients are co-located. Regions in the space of parameters characterizing the base flow velocity profile, i.e. shear layer thickness and ratio of forward to reverse flow velocity, corresponding to convective and absolute instability are identified. The implications of the present results on prior observations of flow instability in other flows such as heated jets and bluff-body stabilized flames is discussed.

Bimodality and re-entrant behaviour in the hierarchical self-assembly of polymeric nanoparticles

We show that a film of a suspension of polymer grafted nanoparticles on a liquid substrate can be employed to create two-dimensional nanostructures with a remarkable variation in the pattern length scales. The presented experiments also reveal the emergence of concentration-dependent bimodal patterns as well as re-entrant behaviour that involves length scales due to dewetting and compositional instabilities. The experimental observations are explained through a gradient dynamics model consisting of coupled evolution equations for the height of the suspension film and the concentration of polymer. Using a Flory-Huggins free energy functional for the polymer solution, we show in a linear stability analysis that the thin film undergoes dewetting and/or compositional instabilities depending on the concentration of the polymer in the solution. We argue that the formation via `hierarchical self-assembly' of various functional nanostructures observed in different systems can be explained as resulting from such an interplay of instabilities.

四边固定载流矩形薄板的磁弹性稳定性分析

本文针对四边固定载流矩形薄板，利用Mathieu方程解的稳定性，研究其在电磁场与机械荷载共同作用下的磁弹性稳定性问题。首先在载流薄板的磁弹性非线性运动方程、物理方程、几何方程、洛仑兹力表达式及电动力学方程的基础上，导出了载流薄板在电磁场与机械荷载共同作用下的磁弹性动力稳定方程，然后应用Galerkin方法将稳定方程整理为Mathieu方程的标准形式，并将薄板的动力稳定性问题归结为对Mathieu方程的求解。利用Mathieu方程的稳定解区域与非稳定解区域的分界，即方程系数λ和η的本征值关系，得出了磁弹性问题失稳临界状态的判别方程。通过具体算例，给出了四边固定载流矩形薄板的磁弹性动力失稳临界状态与相关参量之间的关系曲线，并对计算结果及其变化规律进行了分析讨论。

边坡稳定性分析极限平衡法的简化条件

边坡稳定性分析极限平衡法在求解安全系数时,需要对条块间内力引入假定条件或某些平衡条件在求解时得不到满足.分析这些不同的假设条件对于计算结果的影响,建立一个满足极限平衡原理和合理性条件的数值优化模型.引入载荷系数,利用载荷系数与安全系数之间的单调关系,将模型的求解转化为求解一系列线性规划问题,算法中不存在数值收敛问题,且可以求出安全系数的全局最优解.通过对模型各个约束条件的变化来模拟引入的假定条件,以此探讨各种假定条件或静力平衡条件的缺失对计算结果影响及某些常用方法引入的假定条件可能带来的误差.

三维简化Janbu法分析边坡稳定性的扩展

对三维简化的Janbu法进行了扩展,使得离散后在同一行上的各条块在水平横向的安全系数相等以及在同一列上的条块在水平纵向的安全系数相等,得到不同行(列)的安全系数,最后通过各条块底滑面的几何特性与受力分析给出其独立的安全系数以及各条块潜在的滑动方向.该方法拓展了传统极限平衡方法只给出坡体总体安全系数的思路,可对坡体局部稳定性和潜在滑动方向进行判定.该方法可适用于任意形状的滑面,还可分析包括考虑孔隙水压力、地质体分层、上覆载荷以及地震力等其他形式的外载.

Novel Investigation on Principle of Similarity for Microgravity Two-Phase Flow

根据两流体同心环状流线性稳定性分析的结果 ,对微重力气 /液两相流地面模拟实验所应遵循的相似准则进行了探讨 ,得到了一个新的重力无关性准则 ,即Bond数和环形区流体相的毛细数之比的绝对值不大于 1 .此外 ,微重力气 /液两相流模拟实验还必须满足两个条件 ,即流量比和气相表观Weber数应与所模拟的流动中对应数值相等 . In the present paper, the principle of similarity for two phase flows at microgravity is studied based on the results of the linear stability analysis of the two fluid concentric annular flow configuration. A new criterion of gravity independence, namely the absolute value of the ratio between the Bond number and the capillary number of the phase flowing in the annulus is no more than one, is achieved. It is also pointed out that the flowrate ratio and the gas superficial Weber number must have the same ...

Sedimentation Process of Saturated Sand under Impact Loading

The initial small inhomogeneity of saturated sand could be amplified during the sedimentation process after liquefaction, and cracks could be observed in the sand column. Layers of fine sand could also be found at the exact place where cracks developed and disappeared. The phenomena and the whole process were experimentally shown by X-rays images. To account for the phenomena, a linearized stability analysis of the sedimentation of saturated sand was conducted; however, it did not produce a satisfactory result. A three-phase flow model describing the transportation of fine sand is presented in this paper. It is shown that such a kind of erosion/deposition model was qualitatively in good agreement with the experimental observation.

各向异性非线性强度条件下的边坡稳定性

解释了土体强度各向异性、非线性的物理本质,结合常规直剪试验、三轴试验结果,在前人工作基础上建立了边坡稳定性分析中强度各向异性、非线性的描述方法,其中特别提出了一个各向异性函数.基于Janbu普遍条分法(GPS),运用SPREADSHEET模板程序,提出了一个能将各向异性、非线性强度准则逐点等效到Mohr-Coulomb直线强度准则处理上的迭代方法,准确方便地获得了非线性强度下的边坡稳定性分析.最后的算例展示了方法的使用过程.

Oscillatory Instabilities of Two-Layer Rayleigh-Marangoni-Benard Convection

The Rayleigh–Marangoni–Bénard convective instability (R–M–B instability) in the two-layer systems such as Silicone oil (10cSt)/Fluorinert (FC70) and Silicone oil (2cSt)/water liquids are studied. Both linear instability analysis and nonlinear instability analysis (2D numerical simulation) were performed to study the influence of thermocapillary force on the convective instability of the two-layer system. The results show the strong effects of thermocapillary force at the interface on the time-dependent oscillations at the onset of instability convection. The secondary instability phenomenon found in the real two-layer system of Silicone oil over water could explain the difference in the comparison of the Degen’s experimental observation with the previous linear stability analysis results of Renardy et al.

边坡稳定分析的三维Spencer法

将边坡稳定分析中的二维Spencer法拓展到了三维。该方法对所有条块满足力的平衡以及整体力矩平衡，克服了其他方法中只适用对称问题的缺点，不需已知滑动方向，还可根据滑面的几何特征，进一步得到各条块局部的稳定性系数及其潜在的滑动方向。此外，给出了该方法的实现步骤，其算法的收敛性也较好。最后以不对称问题为算例，并与其他文献中的结果进行比较，验证了该方法的可行性和有效性。

细粒尾矿堆积坝的地震稳定性分析

在分析地震对尾矿坝造成破坏特点的基础上，介绍了尾矿坝地震稳定性分析方法。以龙都尾矿库细粒尾矿堆积坝为例，分析了在当地发生七级地震情况下，尾矿坝在不同状态下的结果，为矿山进行尾矿库的安全管理提供了科学依据。