An analytical connection between temporal and spatio-temporal growth rates in linear stability analysis

We derive an exact formula for the complex frequency in spatio-temporal stability analysis that is valid for arbitrary complex wave numbers. The usefulness of the formula lies in the fact that it depends only on purely temporal quantities, which are easily calculated. We apply the formula to two model dispersion relations: the linearized complex Ginzburg--Landau equation, and a model of wake instability. In the first case, a quadratic truncation of the exact formula applies; in the second, the same quadratic truncation yields an estimate of the parameter values at which the transition to absolute instability occurs; the error in the estimate decreases upon increasing the order of the truncation. We outline ways in which the formula can be used to characterize stability results obtained from purely numerical calculations, and point to a further application in global stability analyses.Comment: 36 pages, 16 figures; Article has been tweaked and reduced in size but essential features remain the same; Supplementary material (16 pages) is also include

Absolute linear instability in laminar and turbulent gas/liquid two-layer channel flow

We study two-phase stratified flow where the bottom layer is a thin laminar liquid and the upper layer is a fully-developed gas flow. The gas flow can be laminar or turbulent. To determine the boundary between convective and absolute instability, we use Orr--Sommerfeld stability theory, and a combination of linear modal analysis and ray analysis. For turbulent gas flow, and for the density ratio r=1000, we find large regions of parameter space that produce absolute instability. These parameter regimes involve viscosity ratios of direct relevance to oil/gas flows. If, instead, the gas layer is laminar, absolute instability persists for the density ratio r=1000, although the convective/absolute stability boundary occurs at a viscosity ratio that is an order of magnitude smaller than in the turbulent case. Two further unstable temporal modes exist in both the laminar and the turbulent cases, one of which can exclude absolute instability. We compare our results with an experimentally-determined flow-regime map, and discuss the potential application of the present method to non-linear analyses.Comment: 33 pages, 20 figure

Inertial coalescence of droplets on a partially wetting substrate

We consider the growth rate of the height of the connecting bridge in rapid surface-tension-driven coalescence of two identical droplets attached on a partially wetting substrate. For a wide range of contact angle values, the height of the bridge grows with time following a power law with a universal exponent of 2/3, up to a threshold time, beyond which a 1/2 exponent results, that is known for coalescence of freely-suspended droplets. In a narrow range of contact angle values close to 90°, this threshold time rapidly vanishes and a 1/2 exponent results for a 90° contact angle. The argument is confirmed by three-dimensional numerical simulations based on a diffuse interface method with adaptive mesh refinement and a volume-of-fluid method

Finite-Weber-Number Motions of Bubbles Through a Nearly Inviscid Liquid

A method is described for computing the motion of bubbles through a liquid under conditions of large Reynolds and finite Weber numbers. Ellipsoidal harmonics are used to approximate the shapes of the bubbles and the flow induced by the bubbles, and a method of summing flows induced by groups of bubbles, using a fast multipole expansion technique is employed so that the computational cost increases only linearly with the number of bubbles. Several problems involving one, two and many bubbles are examined using the method. In particular, it is shown that two bubbles moving towards each other in an impurity-free, inviscid liquid touch each other in a finite time. Conditions for the bubbles to bounce in the presence of non-hydrodynamic forces and the time for bounce when these conditions are satisfied are determined. The added mass and viscous drag coefficients and aspect ratio of bubbles are determined as a function of bubble volume fraction and Weber number

Viscous simulations of shock-bubble interaction and Richtmyer-Meshkov instability

Viscous simulations of shock-bubble interaction and Richtmyer-Meshkov instability are performed using an explicit high-order computational method. The simulations are performed by solving the Navier-Stokes equations associated with two convection equations governing the interface between two fluids. The stiffened equation of state is used to relate the pressure to the total energy of a liquid or a gas. Two-dimensional two-phase flows are considered. The first flow concerns the Richtmyer-Meshkov instability developed on a post-shocked interface between air and sulphur hexafluoride (SF6). The influence of the grid refinement on the instability shape is studied. The second problem deals with a shock wave propagating in air and hitting a cylindrical bubble filled with helium or chlorodifluoromethane (R22). A spatial-time diagram represents the locations of the various pressure waves generated from the interaction between the shock wave and the interface. For both simulations, the numerical results are in agreement with experimental data and visualizations

Determination of Particle Size Distributions from Acoustic Wave Propagation Measurements

The wave equations for the interior and exterior of the particles are ensemble averaged and combined with an analysis by Allegra and Hawley @J. Acoust. Soc. Am. 51, 1545 ~1972!# for the interaction of a single particle with the incident wave to determine the phase speed and attenuation of sound waves propagating through dilute slurries. The theory is shown to compare very well with the measured attenuation. The inverse problem, i.e., the problem of determining the particle size distribution given the attenuation as a function of frequency, is examined using regularization techniques that have been successful for bubbly liquids. It is shown that, unlike the bubbly liquids, the success of solving the inverse problem is limited since it depends strongly on the nature of particles and the frequency range used in inverse calculations

Attenuation of Sound in Concentrated Suspensions: Theory and Experiments

Ensemble-averaged equations are derived for small-amplitude acoustic wave propagation through non-dilute suspensions. The equations are closed by introducing effective properties of the suspension such as the compressibility, density, viscoelasticity, heat capacity, and conductivity. These effective properties are estimated as a function of frequency, particle volume fraction, and physical properties of the individual phases using a self-consistent, effective-medium approximation. The theory is shown to be in excellent agreement with various rigorous analytical results accounting for multiparticle interactions. The theory is also shown to agree well with the experimental data on concentrated suspensions of small polystyrene particles in water obtained by Allegra & Hawley and for glass particles in water obtained in the present study

An artificial neural network stratifies the risks of reintervention and mortality after endovascular aneurysm repair:a retrospective observational study

Background Lifelong surveillance after endovascular repair (EVAR) of abdominal aortic aneurysms (AAA) is considered mandatory to detect potentially life-threatening endograft complications. A minority of patients require reintervention but cannot be predictively identified by existing methods. This study aimed to improve the prediction of endograft complications and mortality, through the application of machine-learning techniques. Methods Patients undergoing EVAR at 2 centres were studied from 2004-2010. Pre-operative aneurysm morphology was quantified and endograft complications were recorded up to 5 years following surgery. An artificial neural networks (ANN) approach was used to predict whether patients would be at low- or high-risk of endograft complications (aortic/limb) or mortality. Centre 1 data were used for training and centre 2 data for validation. ANN performance was assessed by Kaplan-Meier analysis to compare the incidence of aortic complications, limb complications, and mortality; in patients predicted to be low-risk, versus those predicted to be high-risk. Results 761 patients aged 75 +/- 7 years underwent EVAR. Mean follow-up was 36+/- 20 months. An ANN was created from morphological features including angulation/length/areas/diameters/ volume/tortuosity of the aneurysm neck/sac/iliac segments. ANN models predicted endograft complications and mortality with excellent discrimination between a low-risk and high-risk group. In external validation, the 5-year rates of freedom from aortic complications, limb complications and mortality were 95.9% vs 67.9%; 99.3% vs 92.0%; and 87.9% vs 79.3% respectively (p0.001) Conclusion This study presents ANN models that stratify the 5-year risk of endograft complications or mortality using routinely available pre-operative data

Non-isothermal droplet spreading/dewetting and its reversal

International audienceAxisymmetric non-isothermal spreading/dewetting of droplets on a substrate is studied , wherein the surface tension is a function of temperature, resulting in Marangoni stresses. A lubrication theory is first extended to determine the drop shape for spread-ing/dewetting limited by slip. It is demonstrated that an apparent angle inferred from a fitted spherical cap shape does not relate to the contact-line speed as it would under isothermal conditions. Also, a power law for the thermocapillary spreading rate versus time is derived. Results obtained with direct numerical simulations (DNSs), using a slip length down to O(10 −4) times the drop diameter, confirm predictions from lubrication theory. The DNS results further show that the behaviour predicted by the lubrication theory that a cold wall promotes spreading, and a hot wall promotes dewetting, is reversed at sufficiently large contact angles and/or viscosity of the surrounding fluid. This behaviour is summarized in a phase diagram, and a simple model that supports this finding, is presented. Although the key results are found to be robust when accounting for heat conduction in the substrate, a critical thickness of the substrate is identified above which wall conduction significantly modifies wetting behaviour