Thermoacoustic Instability Model With Porous Media: Linear Stability Analysis and the Impact of Porous Media

Author(s):  
Cody S. Dowd ◽  
Joseph W. Meadows

Lean premixed (LPM) combustion systems are susceptible to thermoacoustic instability, which occurs when acoustic pressure oscillations are in phase with the unsteady heat release rates. Porous media has inherent acoustic damping properties, and has been shown to mitigate thermoacoustic instability; however, theoretical models for predicting thermoacoustic instability with porous media do not exist. In the present study, a 1-D model has been developed for the linear stability analysis of the longitudinal modes for a series of constant cross-sectional area ducts with porous media using a n-Tau flame transfer function. By studying the linear regime, the prediction of acoustic growth rates and subsequently the stability of the system is possible. A transfer matrix approach is used to solve for acoustic perturbations of pressure and velocity, stability growth rate, and frequency shift without and with porous media. The Galerkin approximation is used to approximate the stability growth rate and frequency shift, and it is compared to the numerical solution of the governing equations. Porous media is modeled using the following properties: porosity, flow resistivity, effective bulk modulus, and structure factor. The properties of porous media are systematically varied to determine the impact on the eigenfrequencies and stability growth rates. Porous media is shown to increase the stability domain for a range of time delays (Tau) compared to similar cases without porous media.

2018 ◽  
Vol 141 (4) ◽  
Author(s):  
Cody S. Dowd ◽  
Joseph W. Meadows

Lean premixed (LPM) combustion systems are susceptible to thermoacoustic instability, which occurs when acoustic pressure oscillations are in phase with the unsteady heat release rates. Porous media has inherent acoustic damping properties and has been shown to mitigate thermoacoustic instability; however, theoretical models for predicting thermoacoustic instability with porous media do not exist. In the present study, a one-dimensional (1D) model has been developed for the linear stability analysis of the longitudinal modes for a series of constant cross-sectional area ducts with porous media using a n-Tau flame transfer function (FTF). By studying the linear regime, the prediction of acoustic growth rates and subsequently the stability of the system is possible. A transfer matrix approach is used to solve for acoustic perturbations of pressure and velocity, stability growth rate, and frequency shift without and with porous media. The Galerkin approximation is used to approximate the stability growth rate and frequency shift, and it is compared to the numerical solution of the governing equations. Porous media is modeled using the following properties: porosity, flow resistivity, effective bulk modulus, and structure factor. The properties of porous media are systematically varied to determine the impact on the eigenfrequencies and stability growth rates. Porous media is shown to increase the stability domain for a range of time delays (Tau) compared to similar cases without porous media.


1982 ◽  
Vol 22 (05) ◽  
pp. 625-634 ◽  
Author(s):  
David A. Krueger

Abstract Downhole steam generation leads to consideration of reservoir fluid displacement by a mixture of steam and nitrogen. The linear stability analysis of the steam condensation front has been generalized to include a noncondensing gas. Roughly speaking, the addition of nitrogen increases the likelihood of having fingers, but, compared with the no-nitrogen case, the fingers will grow more slowly. Introduction The theory of the stability of flows through porous media has been a subject of interest for more than 25 years, dating back to the pioneering work of Dietz, Chuoke et al., and Saffman and Taylor. They considered injecting one fluid (e.g., water) to force a second fluid (e.g., oil) out of a porous medium. The primary result was that instabilities (fingering) occurred when the driving fluid was more mobile than the driven fluid. Hagoort included multiple fluid phases. Miller generalized the original work to include steam driving water (liquid). He showed that the thermodynamic phase transition (steam to water) introduces two stabilizing effects. The first effect introduces a water/steam velocity ratio as a multiplier of the mobility ratio. This factor is less than one because of the volume change upon condensation. The second effect is the cooling of incipient steam fingers by the surrounding water, which retards their growth. Baker anticipated these effects in a qualitative way to explain his experiments, which showed a more stable displacement by steam than was expected on the basis of mobility ratios alone. Armento and Miller also have considered the stability of the in-situ combustion front in porous media. Their work deals with a region where steam is generated. This paper reformulates Miller's results for a condensation front in a more useful form including general numerical results and extends the theory to include injection of a noncondensing gas (e.g., nitrogen) together with the steam. Depending on the particular situation, the presence of nitrogen can be either stabilizing or destabilizing. The motivation for the generalization comes from enhanced oil recovery projects where the exhaust gases from the steam generator are injected into the reservoir along with the steam. This paper considers perturbations on a flat condensation front that is perpendicular to its velocity. The gravitational force along this velocity is included, but the component of the gravitational force perpendicular to the velocity is not. Thus we include the effect of gravity on fingering, but we do not discuss the gravity override problem. In Stability Analysis we present two steps:determination of the motion of a flat condensation front (details are in the Appendix) andevaluation of the characteristic time for growth or decay of a perturbation of that front. In Results wegive the results for a specific reservoir;discuss the sensitivity of these results to the important reservoir parameters (flow velocities and absolute permeabilities),show that, if surface tension and gravitation are unimportant, the stability condition is independent of the absolute permeability and absolute flow rates, anddiscuss the longest wavelength for a stable perturbation. In the final section we discuss the main conclusions. Stability Analysis We consider a homogeneous porous medium with fluids in two regions as illustrated in Fig. 1. A steam/nitrogen mixture is injected at the left, and water (liquid) and nitrogen are produced at the fight. The linear stability analysis proceeds in two main stages and follows the general methods as discussed by Chandrasekhar and the specific application of Miller. First, we assume that the condensation front is flat, moves with constant velocity, v, and has properties that vary with z alone. SPEJ P. 625^


2001 ◽  
Author(s):  
Davide Valtorta ◽  
Khaled E. Zaazaa ◽  
Ahmed A. Shabana ◽  
Jalil R. Sany

Abstract The lateral stability of railroad vehicles travelling on tangent tracks is one of the important problems that has been the subject of extensive research since the nineteenth century. Early detailed studies of this problem in the twentieth century are the work of Carter and Rocard on the stability of locomotives. The linear theory for the lateral stability analysis has been extensively used in the past and can give good results under certain operating conditions. In this paper, the results obtained using a linear stability analysis are compared with the results obtained using a general nonlinear multibody methodology. In the linear stability analysis, the sources of the instability are investigated using Liapunov’s linear theory and the eigenvalue analysis for a simple wheelset model on a tangent track. The effects of the stiffness of the primary and secondary suspensions on the stability results are investigated. The results obtained for the simple model using the linear approach are compared with the results obtained using a new nonlinear multibody based constrained wheel/rail contact formulation. This comparative numerical study can be used to validate the use of the constrained wheel/rail contact formulation in the study of lateral stability. Similar studies can be used in the future to define the limitations of the linear theory under general operating conditions.


2013 ◽  
Vol 721 ◽  
pp. 268-294 ◽  
Author(s):  
L. Talon ◽  
N. Goyal ◽  
E. Meiburg

AbstractA computational investigation of variable density and viscosity, miscible displacements in horizontal Hele-Shaw cells is presented. As a first step, two-dimensional base states are obtained by means of simulations of the Stokes equations, which are nonlinear due to the dependence of the viscosity on the local concentration. Here, the vertical position of the displacement front is seen to reach a quasisteady equilibrium value, reflecting a balance between viscous and gravitational forces. These base states allow for two instability modes: first, there is the familiar tip instability driven by the unfavourable viscosity contrast of the displacement, which is modulated by the presence of density variations in the gravitational field; second, a gravitational instability occurs at the unstably stratified horizontal interface along the side of the finger. Both of these instability modes are investigated by means of a linear stability analysis. The gravitational mode along the side of the finger is characterized by a wavelength of about one half to one full gap width. It becomes more unstable as the gravity parameter increases, even though the interface is shifted closer to the wall. The growth rate is largest far behind the finger tip, where the interface is both thicker, and located closer to the wall, than near the finger tip. The competing influences of interface thickness and wall proximity are clarified by means of a parametric stability analysis. The tip instability mode represents a gravity-modulated version of the neutrally buoyant mode. The analysis shows that in the presence of density stratification its growth rate increases, while the dominant wavelength decreases. This overall destabilizing effect of gravity is due to the additional terms appearing in the stability equations, which outweigh the stabilizing effects of gravity onto the base state.


2010 ◽  
Vol 652 ◽  
pp. 5-44 ◽  
Author(s):  
S. DUPONT ◽  
F. GOSSELIN ◽  
C. PY ◽  
E. DE LANGRE ◽  
P. HEMON ◽  
...  

In order to investigate the possibility of modelling plant motion at the landscape scale, an equation for crop plant motion, forced by an instantaneous velocity field, is introduced in a large-eddy simulation (LES) airflow model, previously validated over homogeneous and heterogeneous canopies. The canopy is simply represented as a poroelastic continuous medium, which is similar in its discrete form to an infinite row of identical oscillating stems. Only one linear mode of plant vibration is considered. Two-way coupling between plant motion and the wind flow is insured through the drag force term. The coupled model is validated on the basis of a comparison with measured movements of an alfalfa crop canopy. It is also compared with the outputs of a linear stability analysis. The model is shown to reproduce the well-known phenomenon of ‘honami’ which is typical of wave-like crop motions on windy days. The wavelength of the main coherent waving patches, extracted using a bi-orthogonal decomposition (BOD) of the crop velocity fields, is in agreement with that deduced from video recordings. The main spatial and temporal characteristics of these waving patches exhibit the same variation with mean wind velocity as that observed with the measurements. However they differ from the coherent eddy structures of the wind flow at canopy top, so that coherent waving patches cannot be seen as direct signatures of coherent eddy structures. Finally, it is shown that the impact of crop motion on the wind dynamics is negligible for current wind speed values. No lock-in mechanism of coherent eddy structures on plant motion is observed, in contradiction with the linear stability analysis. This discrepancy may be attributed to the presence of a nonlinear saturation mechanism in LES.


2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard

This study seeks to characterise the breakdown of the steady two-dimensional solution in the flow around a 180-degree sharp bend to infinitesimal three-dimensional disturbances using a linear stability analysis. The stability analysis predicts that three-dimensional transition is via a synchronous instability of the steady flows. A highly accurate global linear stability analysis of the flow was conducted with Reynolds number $\mathit{Re}<1150$ and bend opening ratio (ratio of bend width to inlet height) $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 5$. This range of $\mathit{Re}$ and $\unicode[STIX]{x1D6FD}$ captures both steady-state two-dimensional flow solutions and the inception of unsteady two-dimensional flow. For $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 1$, the two-dimensional base flow transitions from steady to unsteady at higher Reynolds number as $\unicode[STIX]{x1D6FD}$ increases. The stability analysis shows that at the onset of instability, the base flow becomes three-dimensionally unstable in two different modes, namely a spanwise oscillating mode for $\unicode[STIX]{x1D6FD}=0.2$ and a spanwise synchronous mode for $\unicode[STIX]{x1D6FD}\geqslant 0.3$. The critical Reynolds number and the spanwise wavelength of perturbations increase as $\unicode[STIX]{x1D6FD}$ increases. For $1<\unicode[STIX]{x1D6FD}\leqslant 2$ both the critical Reynolds number for onset of unsteadiness and the spanwise wavelength decrease as $\unicode[STIX]{x1D6FD}$ increases. Finally, for $2<\unicode[STIX]{x1D6FD}\leqslant 5$, the critical Reynolds number and spanwise wavelength remain almost constant. The linear stability analysis also shows that the base flow becomes unstable to different three-dimensional modes depending on the opening ratio. The modes are found to be localised near the reattachment point of the first recirculation bubble.


1997 ◽  
Vol 352 ◽  
pp. 265-281 ◽  
Author(s):  
A. M. H. BROOKER ◽  
J. C. PATTERSON ◽  
S. W. ARMFIELD

A non-parallel linear stability analysis which utilizes the assumptions made in the parabolized stability equations is applied to the buoyancy-driven flow in a differentially heated cavity. Numerical integration of the complete Navier–Stokes and energy equations is used to validate the non-parallel theory by introducing an oscillatory heat input at the upstream end of the boundary layer. In this way the stability properties are obtained by analysing the evolution of the resulting disturbances. The solutions show that the spatial growth rate and wavenumber are highly dependent on the transverse location and the disturbance flow quantity under consideration. The local solution to the parabolized stability equations accurately predicts the wave properties observed in the direct simulation whereas conventional parallel stability analysis overpredicts the spatial amplification and the wavenumber.


2014 ◽  
Vol 574 ◽  
pp. 15-20
Author(s):  
Zhi Wei Guo ◽  
Si Min Shen ◽  
Wei Min Feng ◽  
Bo Fu Wang

Temporal linear stability of a compressible axisymmetric swirling jet is investigated. The present work extends a previous analysis to include the effects of swirl number on the stability of flow dynamics. Results obtained show that the optimal growth rate of disturbance for azimuthal wavenumber n = -1 is larger than that for n = -2 while the corresponding frequencies for both n increases as axial wavenumber increases. As swirl number q increases, the optimal growth rate of disturbance also increases. What is more, there is an optimal swirl number for small axial wavenumbers, which is different from the situation for medium and large axial wavenumbers.


2005 ◽  
Vol 74 (1) ◽  
pp. 8-12 ◽  
Author(s):  
Jérôme Colin

The linear stability analysis of the shape of a spherical cavity embedded in an infinite-size matrix under stress has been performed when infinitesimal perturbation from sphericity of the rod is assumed to appear by surface diffusion. Developing the perturbation on a basis of complete spherical harmonics, the growth rate of each harmonic Ylm(θ,φ) has been determined and the conditions for the development of the different fluctuations have been discussed as a function of the applied stress and the order l of the perturbation.


Sign in / Sign up

Export Citation Format

Share Document