Acoustic scattering from baffled membranes that are backed by elastic cavities

An elastic membrane backed by a fluid-filled cavity in an elastic body is set into an infinite plane baffle. A time harmonic wave propagating in the acoustic fluid in the upper half-space is incident on the plane. It is assumed that the densities of this fluid and the fluid inside the cavity are small compared with the densities of the membrane and of the elastic walls of the cavity, thus defining a small parameter . Asymptotic expansions of the solution of this scattering problem as →0, that are uniform in the wave number k of the incident wave, are obtained using the method of matched asymptotic expansions. When the frequency of the incident wave is bounded away from the resonant frequencies of the membrane, the cavity fluid, and the elastic body, the resultant wave is a small perturbation (the “outer expansion”) of the specularly reflected wave from a completely rigid plane. However, when the incident wave frequency is near a resonant frequency (the “inner expansion”) then the scattered wave results from the interaction of the acoustic fluid with the membrane, the membrane with the cavity fluid, and finally the cavity fluid with the elastic body, and the resulting scattered field may be “large”. The cavity backed membrane (CBM) was previously analyzed for a rigid cavity wall. In this paper, we study the effects of the elastic cavity walls on modifying the response of the CBM. For incident frequencies near the membrane resonant frequencies, the elasticity of the cavity gives only a higher order (in ) correction to the scattered field. However, near a cavity fluid resonant frequency, and, of course, near an elastic body resonant frequency the elasticity contributes to the scattered field. The method is applied to the two dimensional problem of an infinite strip membrane backed by an infinitely long rectangular cavity. The cavity is formed by two infinitely long rectangular elastic solids. We speculate on the possible significance of the results with respect to viscoelastic membranes and viscoelastic instead of elastic cavity walls for surface sound absorbers.     

The method of fundamental solutions for Stokes flow in a rectangular cavity with cylinders

Numerical solutions based on the method of fundamental solutions are discussed for Stokes flow inside a rectangular cavity in the presence of circular cylinders. The Stokeslets are used as the fundamental solutions to obtain the solution for the flow field by a linear combination of fundamental solutions. Flow results on the cellular structure of flow field resulting from the dynamics of cylinders and the horizontal walls of the cavity are reported for (i) one rotating cylinder in a rectangular cavity with two parallel horizontal sides moving in the same directions as well as in the opposite directions, (ii) two rotating cylinders kept apart in a rectangular cavity with two parallel horizontal sides moving in the same directions as well as in the opposite directions. The effect of aspect ratio of the rectangular cavity, direction of movement of the two parallel horizontal sides of the cavity and the diameter of the rotating cylinder on the flow structure are studied. The flow results obtained for the single cylinder case are in accordance with the results available in the literature. From the computational point of view, the present numerical procedure based on the method of fundamental solutions is efficient and simple to implement as compared to the mesh-dependent schemes, which needs complex mesh generation procedure for the multiply connected geometrical domains considered in this article.     

Finite Difference Method for the Acoustic Radiation of an Elastic Plate Excited by a Turbulent Boundary Layer: A Spectral Domain Solution

A finite difference method is developed to study, on a two-dimensional model, the acoustic pressure radiated when a thin elastic plate, clamped at its boundaries, is excited by a turbulent boundary layer. Consider a homogeneous thin elastic plate clamped at its boundaries and extended to infinity by a plane, perfectly rigid, baffle. This plate closes a rectangular cavity. Both the cavity and the outside domain contain a perfect fluid. The fluid in the cavity is at rest. The fluid in the outside domain moves in the direction parallel to the system plate/baffle with a constant speed. A turbulent boundary layer develops at the interface baffle/plate. The wall pressure fluctuations in this boundary layer generates a vibration of the plate and an acoustic radiation in the two fluid domains. Modeling the wall pressure fluctuations spectrum in a turbulent boundary layer developed over a vibrating surface is a very complex and unresolved task. Ducan and Sirkis [1] proposed a model for the two-way interactions between a membrane and a turbulent flow of fluid. The excitation of the membrane is modeled by a potential flow randomly perturbed. This potential flow is modified by the displacement of the membrane. Howe [2] proposed a model for the turbulent wall pressure fluctuations power spectrum over an elastomeric material. The model presented in this article is based on a hypothesis of one-way interaction between the flow and the structure: the flow generates wall pressure fluctuations which are at the origin of the vibration of the plate, but the vibration of the plate does not modify the characteristics of the flow. A finite difference scheme that incorporates the vibration of the plate and the acoustic pressure inside the fluid cavity has been developed and coupled with a boundary element method that ensures the outside domain coupling. In this paper, we focus on the resolution of the coupled vibration/interior acoustic problem. We compare the results obtained with three numerical methods: (a) a finite difference representation for both the plate displacement and the acoustic pressure inside the cavity; (b) a coupled method involving a finite difference representation for the displacement of the plate and a boundary element method for the interior acoustic pressure; (c) a boundary element method for both the vibration of the plate and the interior acoustic pressure. A comparison of the numerical results obtained with two models of turbulent wall pressure fluctuations spectrums - the Corcos model [3] and the Chase model [4] - is proposed. A difference of 20 dB is found in the vibro-acoustic response of the structure. In [3], this difference is explained by calculating a wavenumber transfer function of the plate. In [6], coupled beam-cavity modes for similar geometry are calculated by the finite difference method. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical investigation of roll convection stability

This study is devoted to the mathematical modeling of Rayleigh-Bénard convection in a rectangular cavity with rigid boundaries. The stability of the roll motion induced by an initial disturbance of special form is studied on the basis of two-dimensional and three-dimensional calculations. Different patterns of flow restructuring with respect to the wavenumber are analyzed on the Rayleigh number range Ra = 1708–7000 for the Prandtl numbers Pr = 1 and 0.71.     

Unsteady phenomena of an oscillating turbulent jet flow inside a cavity: Effect of aspect ratio

Self-sustained oscillatory phenomena in confined flow may occur when a turbulent plane jet is discharging into a rectangular cavity. An experimental set-up was developed and the flow analysis has been made using mainly hot-wire measurements, which were complemented by visualisation data. Previous studies confirmed that periodic oscillations may occur, depending on the location of the jet exit nozzle inside the cavity, and also the distance between the side-walls. The present study deals with the symmetrical interaction between a turbulent plane jet and a rectangular cavity and the influence of the geometrical characteristics of the cavity on the oscillatory motion. The size and aspect ratio of the cavity were varied together with the jet width compared to that of the cavity. The study is carried out both numerically and experimentally. The numerical method solves the unsteady Reynolds averaged Navier–Stokes equations (URANS) together with the continuity equation for an incompressible fluid. The closure of the flow equations system is achieved using a two-scale energy-flux model at high Reynolds number in the core flow coupled with a wall function treatment in the vicinity of the wall boundaries. The fundamental frequency of the oscillatory flow was found to be practically independent of the cavity length. Moreover, the oscillations are attenuated as the cavity width increases, until they disappear for a critical value of the cavity width. Contour maps of the instantaneous flow field are drawn to show the flow pattern evolution at the main phases of oscillation. They are given for several aspect ratios of the cavity, keeping constant values for the cavity width and the jet thickness. The proposed approach may help to investigate further the oscillation mechanisms and the entrainment process occurring in pressure driven jet–cavity interactions.     

参数多项式方法的扩展——参数分布形状特征式法

在深入分析参数多项式方法数学原理的基础上，提出了参数分布形状特征式法。根据流动问题的物理本质及大量试验数据，总结出参数分布的形状特征式，将其代入流动基本方程，就可以快速而准确地求解流场。它是参数多项式方法的一种扩展，具有重要的应用前景。     

A Computational Method for Free-Surface Viscous Flows Through Porous Media

In this paper we have investigated free-surface unsteady viscous flow through porous media using a finite difference technique via primitive formulation. A transformation has been used to convert the irregular flow domain into a rectangular domain. Two illustrations have been presented in case of a standing wave problem and a cavity flow problem with a free surface. The method is found satisfactory for numerical solution for Reynolds number upto 200.     

Flow on oscillating rectangular duct for Maxwell fluid

This paper presents an analysis for the unsteady flow of an incompressible Maxwell fluid in an oscillating rectangular cross section.By using the Fourier and Laplace transforms as mathematical tools,the solutions are presented as a sum of the steady-state and transient solutions.For large time,when the transients disappear,the solution is represented by the steady-state solution.The solutions for the Newtonian fluids appear as limiting cases of the solutions obtained here.In the absence of the frequency of oscillations,we obtain the problem for the flow of the Maxwell fluid in a duct of a rectangular cross-section moving parallel to its length.Finally,the required time to reach the steady-state for sine oscillations of the rectangular duct is obtained by graphical illustrations for different parameters.Moreover,the graphs are sketched for the velocity.     

Bending of an Elastic Rectangular Clamped Plate: Exact Versus ‘Engineering’ Solutions

This paper addresses the fascinating long history of the classical problem of bending of a thin rectangular elastic plate with clamped edges by uniform pressure. Among various mathematical and engineering approaches, a method of superposition proposed by Lamé (1852, 1859) and Mathieu (1881, 1890) and developed by the mathematician Koialovich (1902) and engineers Boobnoff (1902, 1914), Hencky (1913) and Inglis (1925) appears to be very useful for the analysis of distribution of stresses and deflection inside a plate. The object of this paper is both to clarify some purely mathematical questions connected with the solution of the infinite systems of linear algebraic equations and to provide a considerable simplification of the numerical algorithm. This revised version was published online in August 2006 with corrections to the Cover Date.     

封闭水平矩形腔内流体自然对流第一次逆叉形分岔的数值研究

本文采用二阶全展开ETG(Euler-Taylor-Galerkin)分裂步有限元方法,对长宽比为3.5(L/B=3.5,如图 1所示)的封闭矩形腔体内,三种Pr数条件下,定常层流范围内,流体自然对流叉形分岔随Rayleigh数的演化过程进行了数值模拟。研究结果表明,该矩形腔内对应三种Pr数条件下,流体的叉形分岔的演化过程中,在第二次模态Ⅱ型叉形分岔之后,均会出现两个较小尺度涡旋合并,突变为一个较大尺度涡旋的全新叉形分岔模态。即在某临界Ra数两侧,存在定常四涡结构和定常三涡结构两个定常解支,当系统控制参数Ra越过临界值,前者被后者突发性取代,这是完全不同于传统叉形分岔的逆叉形分岔。其数值预报,则采用分半法结合流动拓扑结构及典型截面处速度扩线上鞍点的变化来确定。计算结果表明,在计算的Pr数条件下,随Pr数的增加逆叉形分岔对应临界Ra数的取值也会提高。     

A novel hybrid finite element analysis of inplane singular elastic field around inclusion corners in elastic media

This paper deals with the inplane singular elastic field problems of inclusion corners in elastic media by an ad hoc hybrid-stress finite element method. A one-dimensional finite element method-based eigenanalysis is first applied to determine the order of singularity and the angular dependence of the stress and displacement field, which reflects elastic behavior around an inclusion corner. These numerical eigensolutions are subsequently used to develop a super element that simulates the elastic behavior around the inclusion corner. The super element is finally incorporated with standard four-node hybrid-stress elements to constitute an ad hoc hybrid-stress finite element method for the analysis of local singular stress fields arising from inclusion corners. The singular stress field is expressed by generalized stress intensity factors defined at the inclusion corner. The ad hoc finite element method is used to investigate the problem of a single rectangular or diamond inclusion in isotropic materials under longitudinal tension. Comparison with available numerical results shows the present method is an efficient mesh reducer and yields accurate stress distribution in the near-field region. As applications, the present ad hoc finite element method is extended to discuss the inplane singular elastic field problems of a single rectangular or diamond inclusion in anisotropic materials and of two interacting rectangular inclusions in isotropic materials. In the numerical analysis, the generalized stress intensity factors at the inclusion corner are systematically calculated for various material type, stiffness ratio, shape and spacing position of one or two inclusions in a plate subjected to tension and shear loadings.     

Multiplicity of Steady Two-Dimensional Flows in Two-Sided Lid-Driven Cavities

The two-dimensional steady incompressible flow in rectangular cavities is calculated numerically by a finite volume method. The flow is driven by two opposing cavity side walls which move with constant velocities tangentially to themselves. Depending on the cavity aspect ratio and the two side-wall Reynolds numbers different flow states exist. Their range of existence and the bifurcations between different states are investigated by a continuation method accurately locating the bifurcation points. When both side walls move in opposite directions up to seven solutions are found to exist for the same set of parameters. Three of these are point-symmetric and four are asymmetric with respect to the center of the cavity, if the side-wall Reynolds numbers have the same magnitude. When the walls move in the same direction, up to five different flow states are found. In this case only a single mirror symmetric solution exists for equal Reynolds numbers. Received 9 February 2000 and accepted 9 October 2000     

Pseudostress–velocity formulation for incompressible Navier–Stokes equations

This paper presents a numerical algorithm using the pseudostress–velocity formulation to solve incompressible Newtonian flows. The pseudostress–velocity formulation is a variation of the stress–velocity formulation, which does not require symmetric tensor spaces in the finite element discretization. Hence its discretization is greatly simplified. The discrete system is further decoupled into an H ( div ) problem for the pseudostress and a post‐process resolving the velocity. This can be done conveniently by using the penalty method for steady‐state flows or by using the time discretization for nonsteady‐state flows. We apply this formulation to the 2D lid‐driven cavity problem and study its grid convergence rate. Also, computational results of the time‐dependent‐driven cavity problem and the flow past rectangular problem are reported. Copyright © 2009 John Wiley & Sons, Ltd.     

Improved-Delayed-Detached-Eddy Simulation of cavity-induced transition in hypersonic boundary layer

Hypersonic flow transition from laminar to turbulent due to the surface irregularities, like local cavities, can greatly affect the surface heating and skin friction. In this work, the hypersonic flows over a three-dimensional rectangular cavity with length-to-width-to-depth ratio, L:W:D, of 19.9:3.57:1 at two angles of attack (AoA) were numerically studied with Improved-Delayed-Detached-Eddy Simulation (IDDES) method to highlight the mechanism of transition triggered by the cavity. The present approach was firstly applied to the transonic flow over M219 rectangular cavity. The results, including the fluctuating pressure and frequency, agreed with experiment well. In the hypersonic case at Mach number about 9.6 the cavity is seen as “open” at AoA of −10° but “closed” at AoA of −15° unconventional to the two-dimensional cavity case where the flow always exhibits closed cavity feature when the length-to-depth ratio L/D is larger than 14. For the open cavity flow, the shear layer is basically steady and the flow maintains laminar. For the closed cavity case, the external flow goes into the cavity and impinges on the bottom floor. High intensity streamwise vortices, impingement shock and exit shock are observed causing breakdown of these vortices triggering rapid flow transition.     

Simulation of partial closure of a crack-like cavity with cohesion between the faces in an isotropic medium

A mathematical model for the closure of a crack-like cavity with cohesive end zones in an isotropic medium is constructed using methods of elastic theory. It is assumed that the interaction between the surfaces of the crack-like cavity under the action of body and surface forces can lead to the formation of contact zones on their surfaces. Determination of the unknown parameters characterizing the closure of the crack-like cavity reduces to a system of singular integrodifferential equations. The integral equations are converted to a system of nonlinear algebraic equations which is solved by the method of successive approximations. The contact stresses, the interaction forces between the faces of the crack-like cavity, and the size of the contact zone in which the faces of the crack-like cavities are closed are determined.     

Lifting of Disperse Particles from a Cavity Behind the Front of an Unsteady Shock Wave with a Triangular Velocity Profile

The gas flow in plane shock waves slipping along an impermeable surface with a rectangular cavity where solid disperse particles are suspended is considered numerically. The motion of the gas and particles (gas suspension) is modeled by equations of mechanics of multiphase media. Some laws of the behavior of the dusty cloud in the cavity are established for the case of wave interaction with the cavity.     

3D dynamic stress intensity factors at three rectangular cracks in an infinite elastic medium subjected to a time-harmonic stress wave

Summary Dynamic stresses around three coplanar cracks in an infinite elastic medium are determined in the paper. Two of the cracks are equal, rectangular and symmetrically situated on either side of the centrally located rectangular crack. Time-harmonic normal traction acts on each surface of the three cracks. To solve the problem, two kind of solutions are superposed: one is a solution for a rectangular crack in an infinite elastic medium, and the other one is that for two rectangular cracks in an infinite elastic medium. The unknown coefficients in the combined solution are determined by applying the boundary conditions at the surfaces of the cracks. Finally, stress intensity factors are calculated numerically for several crack configurations. Received 14 July 1998; accepted for publication 2 December 1998     

NUMERICAL SIMULATION OF TWO–DIMENSIONAL BUOYANCY–DRIVEN TURBULENCE IN A TALL RECTANGULAR CAVITY

A computational study of natural convection of air in a tall rectangular cavity with 4:1 aspect ratio is conducted. In an effort to investigate the applicability of the Boussinesq approximation to turbulent flow simulation, the cavity is differentially heated from the sides and is insulated at the ends at a Rayleigh number of 10⁹. Starting from quiescent and isothermal flow conditions, the flow is driven to turbulence without any artificial perturbations. The computer programme developed integrates the two-dimensional, time-dependent Navier–Stokes equations with the Boussinesq approximation and the energy equation by a time-accurate method on a stretched, staggered grid. The simulation proceeds to a statistically steady solution in which large-scale structures are found in the mean. Both mean and fluctuating quantities provide good agreement with experimental results.     

Experimental measurement of flow past cavities of different shapes

Experiments are carried out in order to investigate the flow structure past a rectangular, triangular and semi-circular cavity of length-to-depth ratio of 2:1 using the Particle Image Velocimetry (PIV) technique. The experiments are performed in a large scale water channel with three different upstream velocities resulting in Reynolds numbers of 1230, 1460 and 1700, based on inflow momentum thickness, for each cavity type. Contours of constant averaged streamwise and transverse components of velocity, contours of constant averaged vorticity, Reynolds stress and streamline plots for each cavity type for the aforementioned three Reynolds numbers are presented. In addition, streamwise velocity, Reynolds stress and turbulence intensity values are compared for all cavity types. Effect of cavity shape on flow structure within the cavity is discussed in detail. Moreover, spectrum of instantaneous streamwise velocity fluctuations in shear layer near the downstream of the leading corner and the upstream of the trailing corner of the cavities are obtained and it was found that no organized oscillations are present in the flow; rectangular and triangular cavities have the largest amplitudes while semi-circular cavity has the smallest. Calculated turbulence intensities also reveal that the maximum turbulence intensities occur at cavity lid in the centerline section and rectangular and triangular cavities have larger turbulence intensity compared to semi-circular cavity.     

Calculation of a conductive-laminar thermo-convection regime of a Newtonian fluid in a rectangular cavity with isothermal vertical boundaries

A method of calculating a free-convection flow of a viscous incompressible medium in a rectangular cavity is proposed. The lateral boundaries of the cavity have the same constant temperature, the top and bottomboundaries are thermally insulated, and the half-width to height ratio is much smaller than unity. A decomposition of the flow domain into zones with ascending and descending flows made it possible to formulate and solve analytically a coupled parabolic system of linear equations based on the classic Oberbeck-Boussinesq assumptions.