Fluid-structure interaction with mean flow: Over-scattering and unstable resonance growth

It is known theoretically [1-3] that infinitely long fluid loaded plates in mean flow exhibit a range of unusual phenomena in the 'long time' limit. These include convective instability, absolute instability and negative energy waves which are destabilized by dissipation. However, structures are necessarily of finite length and may have discontinuities. Moreover, linear instability waves can only grow over a limited number of cycles before non-linear effects become dominant. We have undertaken an analytical and computational study to investigate the response of finite, discontinuous plates to ascertain if these unusual effects might be realized in practice. Analytically, we take a "wave scattering" [2,4] - as opposed to a "modal superposition" [5] - view of the fluttering plate problem. First, we solve for the scattering coefficients of localized plate discontinuities and identify a range of parameter space, well outside the convective instability regime, where over-scattering or amplified reflection/transmission occurs. These are scattering processes that draw energy from the mean flow into the plate. Next, we use the Wiener-Hopf technique to solve for the scattering coefficients from the leading and trailing edges of a baffled plate. Finally, we construct the response of a finite, baffled plate by a superposition of infinite plate propagating waves continuously scattering off the plate ends and solve for the unstable resonance frequencies and temporal growth rates for long plates. We present a comparison between our computational results and the infinite plate theory. In particular, the resonance response of a moderately sized plate is shown to be in excellent agreement with our long plate analytical predictions. Copyright © 2010 by ASME.

Non-modal instabilities of two-dimensional disturbances in plane Couette flow of a power-law fluid

Instabilities of fluid flows have traditionally been investigated by normal mode analysis, i.e. by linearizing the equations of flow and testing for unstable eigenvalues of the linearized problem. However, the results of eigenvalue analysis agree poorly in many cases with experiments, especially for shear flows. In this paper we study the instabilities of two-dimensional Couette flow of a polymeric fluid in the framework of non-modal stability theory rather than normal mode analysis. A power-law model is used to describe the polymeric liquid. We focus on the response to external excitations and initial conditions by examining the pseudospectra structures and the transient energy growths. For both Newtonian and non-Newtonian flows, the results show that there can be a rather large transient growth even though the linear operator of Couette flow has no unstable eigenvalue. The effects of non-Newtonian viscosity on the transient behaviors are examined in this study. The results show that the "shear-thinning/shear-thickening" effect increases/decreases the amplitude of responses to external excitations and initial conditions. (C) 2010 Elsevier B.V. All rights reserved.

Viscous Flow and Fluid Dynamics

Exam questions and solutions for a third year maths course. Diagrams for the questions are all together in the support.zip file, as .eps files

The unsteady flow of a weakly compressible fluid in a thin porous layer. I: Two-dimensional theory

We consider the problem of determining the pressure and velocity fields for a weakly compressible fluid flowing in a two-dimensional reservoir in an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting or extracting fluid. Numerical solution of this problem may be expensive, particularly in the case that the depth scale of the layer h is small compared to the horizontal length scale l. This is a situation which occurs frequently in the application to oil reservoir recovery. Under the assumption that epsilon=h/l<<1, we show that the pressure field varies only in the horizontal direction away from the wells (the outer region). We construct two-term asymptotic expansions in epsilon in both the inner (near the wells) and outer regions and use the asymptotic matching principle to derive analytical expressions for all significant process quantities. This approach, via the method of matched asymptotic expansions, takes advantage of the small aspect ratio of the reservoir, epsilon, at precisely the stage where full numerical computations become stiff, and also reveals the detailed structure of the dynamics of the flow, both in the neighborhood of wells and away from wells.

The unsteady flow of a weakly compressible fluid in a thin porous layer II: three-dimensional theory

We consider the problem of determining the pressure and velocity fields for a weakly compressible fluid flowing in a three-dimensional layer, composed of an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting and/or extracting fluid. Numerical solution of this three-dimensional evolution problem may be expensive, particularly in the case that the depth scale of the layer h is small compared to the horizontal length scale l, a situation which occurs frequently in the application to oil and gas reservoir recovery and which leads to significant stiffness in the numerical problem. Under the assumption that $\epsilon\propto h/l\ll 1$, we show that, to leading order in $\epsilon$, the pressure field varies only in the horizontal directions away from the wells (the outer region). We construct asymptotic expansions in $\epsilon$ in both the inner (near the wells) and outer regions and use the asymptotic matching principle to derive expressions for all significant process quantities. The only computations required are for the solution of non-stiff linear, elliptic, two-dimensional boundary-value, and eigenvalue problems. This approach, via the method of matched asymptotic expansions, takes advantage of the small aspect ratio of the layer, $\epsilon$, at precisely the stage where full numerical computations become stiff, and also reveals the detailed structure of the dynamics of the flow, both in the neighbourhood of wells and away from wells.

The unsteady flow of a weakly compressible fluid in a thin porous layer III: three-dimensional computations

We describe a novel method for determining the pressure and velocity fields for a weakly compressible fluid flowing in a thin three-dimensional layer composed of an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting and/or extracting fluid. Our approach uses the method of matched asymptotic expansions to derive expressions for all significant process quantities, the computation of which requires only the solution of linear, elliptic, two-dimensional boundary value and eigenvalue problems. In this article, we provide full implementation details and present numerical results demonstrating the efficiency and accuracy of our scheme.