Local stability analysis and eigenvalue sensitivity of reacting bluff-body wakes

This paper presents an experimental and theoretical investigation of high-Reynolds-number low-density reacting wakes near a hydrodynamic Hopf bifurcation. This configuration is applicable to the wake flows that are commonly used to stabilize flames in high-velocity flows. First, an experimental study is conducted to measure the limit-cycle oscillation of this reacting bluff-body wake. The experiment is repeated while independently varying the bluff-body lip velocity and the density ratio across the flame. In all cases, the wake exhibits a sinuous oscillation. Linear stability analysis is performed on the measured time-averaged velocity and density fields. In the first stage of this analysis, a local spatiotemporal stability analysis is performed on the measured time-averaged velocity and density fields. The stability analysis results are compared to the experimental measurement and demonstrate that the local stability analysis correctly captures the influence of the lip-velocity and density-ratio parameters on the sinuous mode. In the second stage of the analysis, the linear direct and adjoint global modes are estimated by combining the local results. The sensitivity of the eigenvalue to changes in intrinsic feedback mechanisms is found by combining the direct and adjoint global modes. This is referred to as the eigenvalue sensitivity throughout the paper for reasons of brevity. The predicted global mode frequency is consistently within 10 % of the measured value, and the linear global mode shape closely resembles the measured nonlinear oscillations. The adjoint global mode reveals that the oscillation is strongly sensitive to open-loop forcing in the shear layers. The eigenvalue sensitivity identifies a wavemaker in the recirculation zone of the wake. A parametric study shows that these regions change little when the density ratio and lip velocity change. In the third stage of the analysis, the stability analysis is repeated for the varicose hydrodynamic mode. Although not physically observed in this unforced flow, the varicose mode can lock into longitudinal acoustic waves and cause thermoacoustic oscillations to occur. The paper shows that the local stability analysis successfully predicts the global hydrodynamic stability characteristics of this flow and shows that experimental data can be post-processed with this method in order to identify the wavemaker regions and the regions that are most sensitive to external forcing, for example from acoustic waves.Matthew Juniper acknowledges the financial support of the European Research Council under grant 2590620. Tim Lieuwen acknowledges the support of the University Turbine Systems Research (contract #DE-FC21-92MC29061) program under contract monitor Dr. Mark Freeman, and the Air Force Office of Scientific Research (contract #FA9550- 12-1-0107/RC657), under contract monitor Dr. Chiping Lee.This is the author accepted manuscript. The final version is available from Cambridge University Press via http://dx.doi.org/10.1017/jfm.2015.72

The structural sensitivity of open shear flows calculated with a local stability analysis

The structural sensitivity shows where an instability of a fluid flow is most sensitive to changes in internal feedback mechanisms. It is formed from the overlap of the flow's direct and adjoint global modes. These global modes are usually calculated with 2D or 3D global stability analyses, which can be very computationally expensive. For weakly non-parallel flows the direct global mode can also be calculated with a local stability analysis, which is orders of magnitude cheaper. In this theoretical paper we show that, if the direct global mode has been calculated with a local analysis, then the adjoint global mode follows at little extra cost. We also show that the maximum of the structural sensitivity is the location at which the local k+ and k- branches have the same imaginary value. Finally, we use the local analysis to derive the structural sensitivity of two flows: a confined co-flow wake at Re = 400, for which it works very well, and the flow behind a cylinder at Re = 50, for which it works reasonably well. As expected, we find that the local analysis becomes less accurate when the flow becomes less parallel.This is the preprint version of the manuscript as submitted to the journal. The final version will be published by Elsevier

Matrix factorizations for quantum complete intersections

We introduce twisted matrix factorizations for quantum complete intersections of codimension two. For such an algebra, we show that in a given dimension, almost all the indecomposable modules with bounded minimal projective resolutions correspond to such matrix factorizations.Comment: 13 page

State-space realization of a describing function

The describing function is a powerful tool for characterizing nonlinear dynamical systems in the frequency domain. In some cases, it is the only available description of a nonlinear operator characterizing a certain subcomponent of the system. This paper presents a methodology to provide a state-space realization of one given describing function, in order to allow the study of the system in the time domain as well. The realization is based on Hammerstein models and Fourier–Bessel series. It can be embedded in time domain simulations of complex configurations with many nonlinear elements interacting, accurately describing the nonlinear saturation of the system. The technique is applied to an example application in the field of combustion instability, featuring self-excited thermoacoustic oscillations. We benchmark the performance of the tool comparing the results with a frequency domain analysis of the same system, obtaining good agreement between the two formulations.This work was supported by the European Research Council through project ALORS N.259620 and by the German Research Association for Combustion Engines (Forschungsvereinigung Verbrennungskraftmaschinen e.V. FVV).This is the accepted manuscript of a paper published in Nonlinear Dynamics (Ghirardo G, Ćosić B, Juniper MP, Moeck JP, Nonlinear Dynamics, 2015, doi:10.1007/s11071-015-2134-x). The final version is available at http://dx.doi.org/10.1007/s11071-015-2134-

Cost-effectiveness of asthma control: an economic appraisal of the GOAL study

<i>Background</i>: The Gaining Optimal Asthma ControL (GOAL) study has shown the superiority of a combination of salmeterol/fluticasone propionate (SFC) compared with fluticasone propionate alone (FP) in terms of improving guideline defined asthma control. <i>Methods</i>: Clinical and economic data were taken from the GOAL study, supplemented with data on health related quality of life, in order to estimate the cost per quality adjusted life year (QALY) results for each of three strata (previously corticosteroid-free, low- and moderate-dose corticosteroid users). A series of statistical models of trial outcomes was used to construct cost effectiveness estimates across the strata of the multinational GOAL study including adjustment to the UK experience. Uncertainty was handled using the non-parametric bootstrap. Cost-effectiveness was compared with other treatments for chronic conditions. <i>Result</i>: Salmeterol/fluticasone propionate improved the proportion of patients achieving totally and well-controlled weeks resulting in a similar QALY gain across the three strata of GOAL. Additional costs of treatment were greatest in stratum 1 and least in stratum 3, with some of the costs offset by reduced health care resource use. Cost-effectiveness by stratum was £7600 (95% CI: £4800–10 700) per QALY gained for stratum 3; £11 000 (£8600–14 600) per QALY gained for stratum 2; and £13 700 (£11 000–18 300) per QALY gained for stratum 1. <i>Conclusion</i>: The GOAL study previously demonstrated the improvement in total control associated with the use of SFC compared with FP alone. This study suggests that this improvement in control is associated with cost-per-QALY figures that compare favourably with other uses of scarce health care resources

A viscoelastic deadly fluid in carnivorous pitcher plants

Background : The carnivorous plants of the genus Nepenthes, widely distributed in the Asian tropics, rely mostly on nutrients derived from arthropods trapped in their pitcher-shaped leaves and digested by their enzymatic fluid. The genus exhibits a great diversity of prey and pitcher forms and its mechanism of trapping has long intrigued scientists. The slippery inner surfaces of the pitchers, which can be waxy or highly wettable, have so far been considered as the key trapping devices. However, the occurrence of species lacking such epidermal specializations but still effective at trapping insects suggests the possible implication of other mechanisms. Methodology/Principal Findings : Using a combination of insect bioassays, high-speed video and rheological measurements, we show that the digestive fluid of Nepenthes rafflesiana is highly viscoelastic and that this physical property is crucial for the retention of insects in its traps. Trapping efficiency is shown to remain strong even when the fluid is highly diluted by water, as long as the elastic relaxation time of the fluid is higher than the typical time scale of insect movements. Conclusions/Significance : This finding challenges the common classification of Nepenthes pitchers as simple passive traps and is of great adaptive significance for these tropical plants, which are often submitted to high rainfalls and variations in fluid concentration. The viscoelastic trap constitutes a cryptic but potentially widespread adaptation of Nepenthes species and could be a homologous trait shared through common ancestry with the sundew (Drosera) flypaper plants. Such large production of a highly viscoelastic biopolymer fluid in permanent pools is nevertheless unique in the plant kingdom and suggests novel applications for pest control

On acceleration of Krylov-subspace-based Newton and Arnoldi iterations for incompressible CFD: replacing time steppers and generation of initial guess

We propose two techniques aimed at improving the convergence rate of steady state and eigenvalue solvers preconditioned by the inverse Stokes operator and realized via time-stepping. First, we suggest a generalization of the Stokes operator so that the resulting preconditioner operator depends on several parameters and whose action preserves zero divergence and boundary conditions. The parameters can be tuned for each problem to speed up the convergence of a Krylov-subspace-based linear algebra solver. This operator can be inverted by the Uzawa-like algorithm, and does not need a time-stepping. Second, we propose to generate an initial guess of steady flow, leading eigenvalue and eigenvector using orthogonal projection on a divergence-free basis satisfying all boundary conditions. The approach, including the two proposed techniques, is illustrated on the solution of the linear stability problem for laterally heated square and cubic cavities