Acoustic radiation and scattering from cylindrical bodies and analysis of transducer arrays

New mathematical methods to analytically investigate linear acoustic radiation and scattering from cylindrical bodies and transducer arrays are presented. Three problems of interest involving cylinders in an infinite fluid are studied. In all the three problems, the Helmholtz equation is used to model propagation through the fluid and the beam patterns of arrays of transducers are studied. In the first problem, a method is presented to determine the omni-directional and directional far-field pressures radiated by a cylindrical transducer array in an infinite rigid cylindrical baffle. The solution to the Helmholtz equation and the displacement continuity condition at the interface between the array and the surrounding water are used to determine the pressure. The displacement of the surface of each transducer is in the direction of the normal to the array and is assumed to be uniform. Expressions are derived for the pressure radiated by a sector of the array vibrating in-phase, the entire array vibrating in-phase, and a sector of the array phase-shaded to simulate radiation from a rectangular piston. It is shown that the uniform displacement required for generating a source level of 220 dB ref. μPa @ 1m that is omni directional in the azimuthal plane is in the order of 1 micron for typical arrays. Numerical results are presented to show that there is only a small difference between the on-axis pressures radiated by phased cylindrical arrays and planar arrays. The problem is of interest because cylindrical arrays of projectors are often used to search for underwater objects. In the second problem, the errors, when using data-independent, classical, energy and split beam correlation methods, in finding the direction of arrival (DOA) of a plane acoustic wave, caused by the presence of a solid circular elastic cylindrical stiffener near a linear array of hydrophones, are investigated. Scattering from the effectively infinite cylinder is modeled using the exact axisymmetric equations of motion and the total pressures at the hydrophone locations are computed. The effect of the radius of the cylinder, a, the distance between the cylinder and the array, b, the number of hydrophones in the array, 2H, and the angle of incidence of the wave, α, on the error in finding the DOA are illustrated using numerical results. For an array that is about 30 times the wavelength and for small angles of incidence (α<10), the error in finding the DOA using the energy method is less than that using the split beam correlation method with beam steered to α; and in some cases, the error increases when b increases; and the errors in finding the DOA using the energy method and the split beam correlation method with beam steered to α vary approximately as a7 / 4 . The problem is of interest because elastic stiffeners – in nearly acoustically transparent sonar domes that are used to protect arrays of transducers – scatter waves that are incident on it and cause an error in the estimated direction of arrival of the wave. In the third problem, a high-frequency ray-acoustics method is presented and used to determine the interior pressure field when a plane wave is normally incident on a fluid cylinder embedded in another infinite fluid. The pressure field is determined by using geometrical and physical acoustics. The interior pressure is expressed as the sum of the pressures due to all rays that pass through a point. Numerical results are presented for ka = 20 to 100 where k is the acoustic wavenumber of the exterior fluid and a is the radius of the cylinder. The results are in good agreement with those obtained using field theory. The directional responses, to the plane wave, of sectors of a circular array of uniformly distributed hydrophones in the embedded cylinder are then computed. The sectors are used to simulate linear arrays with uniformly distributed normals by using delays. The directional responses are compared with the output from an array in an infinite homogenous fluid. These outputs are of interest as they are used to determine the direction of arrival of the plane wave. Numerical results are presented for a circular array with 32 hydrophones and 12 hydrophones in each sector. The problem is of interest because arrays of hydrophones are housed inside sonar domes and acoustic plane waves from distant sources are scattered by the dome filled with fresh water and cause deterioration in the performance of the array.

THE PARAMETERS OF THE LEWIS METRIC FOR THE WEYL CLASS

We provide physical interpretation for the four parameters of the stationary Lewis metric restricted to the Weyl class. Matching this spacetime to a completely anisotropic, rigidly rotating, fluid cylinder, we obtain from the junction conditions that one of these parameters is proportional to the vorticity of the source. From the Newtonian approximation a second parameter is found to be proportional to the energy per unit of length. The remaining two parameters may be associated to a gravitational analog of the Aharanov-Bohm effect. We prove, using the Cartan scalars, that the Weyl class metric and static Levi-Civita metric are locally equivalent, i.e., indistinguishable in terms of its curvature.

On the parameters of the Lewis metric for the Lewis class

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

3D CFD Simulation of Vortex-induced Vibration of Cylinder

The 3D flow around a circular cylinder free to oscillate transversely to the free stream was simulated using Computational Fluid Dynamics (CFD) and the Spalart-Allmaras Detached Eddy Simulation (DES) turbulence model for a Reynolds number Re = 10(4). Simulations were carried out for a small mass-damping parameter m*zeta = 0.00858, where m* = 3.3 and zeta = 0.0026. We found good agreement between the numerical results and experimental data. The simulations predicted the high observed amplitudes of the upper branch of vortex-induced vibrations for low mass-damping parameters.

Discrete approximation to the global spectrum of the tangent operator for flow past a circular cylinder

This paper deals with the calculation of the discrete approximation to the full spectrum for the tangent operator for the stability problem of the symmetric flow past a circular cylinder. It is also concerned with the localization of the Hopf bifurcation in laminar flow past a cylinder, when the stationary solution loses stability and often becomes periodic in time. The main problem is to determine the critical Reynolds number for which a pair of eigenvalues crosses the imaginary axis. We thus present a divergence-free method, based on a decoupling of the vector of velocities in the saddle-point system from the vector of pressures, allowing the computation of eigenvalues, from which we can deduce the fundamental frequency of the time-periodic solution. The calculation showed that stability is lost through a symmetry-breaking Hopf bifurcation and that the critical Reynolds number is in agreement with the value presented in reported computations. (c) 2007 IMACS. Published by Elsevier B.V. All rights reserved.

Turbulent flow around a rotating stepped cylinder

Direct numerical simulation (DNS) of turbulent flow around a rotating cylinder with two backward-facing steps axisymmetrically mounted in the circumferential direction was performed and compared with DNS of plane backward-facing step flow (PBSF) of Le [J. Fluid Mech. 330, 349 (1997)]. The original motivation of this work stemmed from the efforts to design a simple device which can generate flows of high turbulence intensity at low cost for corrosion researchers. It turned out that the current flow shows flow structures quite similar to those of PBSF downstream of the step, even though configurations of the two flows are totally different from one another. The stepped cylinder appears to be a cost-effective tool in the generation of flow structures similar to those of PBSF. (C) 2002 American Institute of Physics.

Horizontal rotating cylinder - A compact apparatus for studying the effect of water wetting on carbon dioxide corrosion of mild steel

Water wetting is a crucial issue in carbon dioxide (CO.) corrosion of multiphase flow pipelines made from mild steel. This study demonstrates the use of a novel benchtop apparatus, a horizontal rotating cylinder, to study the effect of water wetting on CO2 corrosion of mild steel in two-phase flow. The setup is similar to a standard rotating cylinder except for its horizontal orientation and the presence of two phases-typically water and oil. The apparatus has been tested by using mass-transfer measurements and CO2 corrosion measurements in single-phase water flow. CO2 corrosion measurements were subsequently performed using a water/hexane mixture with water cuts varying between 5% and 50%. While the metal surface was primarily hydrophilic under stagnant. conditions, a variety of dynamic water wetting situations was encountered as the water cut and fluid velocity were altered. Threshold velocities were identified at various water cuts when the surface became oil-wet and corrosion stopped.

Influence of fluid concentration on freezing point depression and thermal conductivity of frozen orange juice

The freezing point depression (FPD) of orange juice at different concentrations was measured by using a simple apparatus. Results showed that the initial freezing point decreased approximately 90% with the increase of juice concentration between 46degrees and 66degrees Brix (water content respectively between 52.8 and 32.8% w/w). The thermal conductivity of orange juice as a function of fluid concentration was also investigated by using a coaxial dual-cylinder apparatus. Below the freezing point, the thermal conductivity was strongly affected by both the orange juice concentration and temperature. Simple equations in terms of water content and temperature could be adjusted to experimental data of FPD and thermal conductivity.

Rheological properties and fluid dynamics of egg yolk

The rheological behavior of egg yolk was studied at a range of temperatures (277-333 K) using a concentric cylinder viscometer. Rheological behavior was pseudoplastic and flow curves fitted by the power law model. The consistency and behavior indexes, dependent on temperature, were expressed by an Arrhenius-type equation. The rheological parameters, together with experimental values of pressure loss in tube flow were used to calculate friction factors. The good agreement between predicted and observed values confirmed the reliability of the equations proposed for describing the flow behavior of the egg yolk. (c) 2005 Published by Elsevier Ltd.

Rheological properties and fluid dynamics of coffee extract

The rheological behavior of coffee extract with different water contents (49 to 90%) was studied at a wide range of temperatures (274 to 365 K) using a concentric cylinder rheometer. The flow curves followed different models depending on the concentration and temperature level. Newtonian behavior was observed at high values of water content and temperature, changing to power law as these values were decreased. The Newtonian viscosity as well as the consistency and behavior index could be well correlated by functions simultaneously dependent on temperature and water content. The rheological parameters, together with experimental values of pressure loss in tube flow, were used to calculate friction factors. These showed to be in good agreement with those resulting from classical theoretical and empirical equations, thus confirming the reliability of the rheological measurements.

On nonlinear dynamics in a free fluid surface of a tank excited by a non-ideal power source

We investigate the nonlinear oscillations in a free surface of a fluid in a cylinder tank excited by non-ideal power source, an electric motor with limited power supply. We study the possibility of parametric resonance in this system, showing that the excitation mechanism can generate chaotic response. Additionally, the dynamics of parametrically excited surface waves in the tank can reveal new characteristics of the system. The fluid-dynamic system is modeled in such way as to obtain a nonlinear differential equation system. Numerical experiments are carried out to find the regions of chaotic solutions. Simulation results are presented as phase-portrait diagrams characterizing the resonant vibrations of free fluid surface and the existence of several types of regular and chaotic attractors. We also describe the energy transfer in the interaction process between the hydrodynamic system and the electric motor. Copyright © 2011 by ASME.

The role of wake stiffness on the wake-induced vibration of the downstream cylinder of a tandem pair

When a pair of tandem cylinders is immersed in a flow the downstream cylinder can be excited into wake-induced vibrations (WIV) due to the interaction with vortices coming from the upstream cylinder. Assi, Bearman & Meneghini ( J. Fluid Mech. , vol. 661, 2010, pp. 365–401) concluded that the WIV excitation mechanism has its origin in the unsteady vortex–structure interaction encountered by the cylinder as it oscillates across the wake. In the present paper we investigate how the cylinder responds to that excitation, characterising the amplitude and frequency of response and its dependency on other parameters of the system. We introduce the concept of wake stiffness , a fluid dynamic effect that can be associated, to a first approximation, with a linear spring with stiffness proportional to Re and to the steady lift force occurring for staggered cylinders. By a series of experiments with a cylinder mounted on a base without springs we verify that such wake stiffness is not only strong enough to sustain oscillatory motion, but can also dominate over the structural stiffness of the system. We conclude that while unsteady vortex–structure interactions provide the energy input to sustain the vibrations, it is the wake stiffness phenomenon that defines the character of the WIV response

Multilevel Domain Decomposition Algorithms for Monolithic Fluid-Structure Interaction Problems with Application to Haemodynamics

Finite element techniques for solving the problem of fluid-structure interaction of an elastic solid material in a laminar incompressible viscous flow are described. The mathematical problem consists of the Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian formulation coupled with a non-linear structure model, considering the problem as one continuum. The coupling between the structure and the fluid is enforced inside a monolithic framework which computes simultaneously for the fluid and the structure unknowns within a unique solver. We used the well-known Crouzeix-Raviart finite element pair for discretization in space and the method of lines for discretization in time. A stability result using the Backward-Euler time-stepping scheme for both fluid and solid part and the finite element method for the space discretization has been proved. The resulting linear system has been solved by multilevel domain decomposition techniques. Our strategy is to solve several local subproblems over subdomain patches using the Schur-complement or GMRES smoother within a multigrid iterative solver. For validation and evaluation of the accuracy of the proposed methodology, we present corresponding results for a set of two FSI benchmark configurations which describe the self-induced elastic deformation of a beam attached to a cylinder in a laminar channel flow, allowing stationary as well as periodically oscillating deformations, and for a benchmark proposed by COMSOL multiphysics where a narrow vertical structure attached to the bottom wall of a channel bends under the force due to both viscous drag and pressure. Then, as an example of fluid-structure interaction in biomedical problems, we considered the academic numerical test which consists in simulating the pressure wave propagation through a straight compliant vessel. All the tests show the applicability and the numerical efficiency of our approach to both two-dimensional and three-dimensional problems.

Apatite as an indicator of fluid salinity: An experimental study of chlorine and fluorine partitioning in subducted sediments

In order to constrain the salinity of subduction zone fluids, piston-cylinder experiments have been conducted to investigate the partitioning behaviour of Cl and F in subducted sediments. These experiments were performed at H2O-undersaturated conditions with a synthetic pelite starting composition containing 800 ppm Cl, over a pressure and temperature range of 2.5–4.5 GPa and 630–900 °C. Repetitive experiments were conducted with 1900 ppm Cl + 1000 ppm F, and 2100 ppm Cl. Apatite represents the most Cl-abundant mineral phase, with Cl concentration varying in the range 0.1–2.82 wt%. Affinity for Cl decreases over the following sequence: aqueous fluid > apatite ⩾ melt > other hydrous minerals (phengite, biotite and amphibole). It was found that addition of F to the Cl-bearing starting composition significantly lowers the Cl partition coefficients between apatite and melt (DClAp–melt) and apatite and aqueous fluid (DClAp–aq). Cl–OH exchange coefficients between apatite and melt (KdCl–OHAp–melt) and apatite and aqueous fluid (KdCl–OHAp–aq) were subsequently calculated. KdCl–OHAp–melt was found to vary from 1 to 58, showing an increase with temperature and a decrease with pressure and displaying a regular decrease with increasing H2O content in melt. Mole fractions of Cl and OH in melt were calculated based on an ideal mixing model for H2O, OH, O, Cl and F. The Cl contents of other hydrous minerals (phengite, biotite and amphibole) fall between 200 and 800 ppm, with resultant Cl partition coefficients from 0.02 to 0.49, appearing independent of the bulk Cl and F content. Preliminary data from this study show that the partitioning behaviour of F is strongly in favour of apatite relative to melt and phengite, with DFAp–melt = 15–51. Apatites from representative eclogite facies metasediments were examined and found to have low Cl contents close to ∼100 ppm. Calculations using our experimentally determined KdCl–OHAp–aq of 0.004 at 2.5 GPa, 630 °C indicate a low salinity character (0.5–2 wt% NaCleq) for the fluid formed during dehydration of subducted oceanic sediment at ∼80 km depth.

The scenario of two-dimensional instabilities of the cylinder wake under EHD forcing: A linear stability analysis

We propose to study the stability properties of an air flow wake forced by a dielectric barrier discharge (DBD) actuator, which is a type of electrohydrodynamic (EHD) actuator. These actuators add momentum to the flow around a cylinder in regions close to the wall and, in our case, are symmetrically disposed near the boundary layer separation point. Since the forcing frequencies, typical of DBD, are much higher than the natural shedding frequency of the flow, we will be considering the forcing actuation as stationary. In the first part, the flow around a circular cylinder modified by EHD actuators will be experimentally studied by means of particle image velocimetry (PIV). In the second part, the EHD actuators have been numerically implemented as a boundary condition on the cylinder surface. Using this boundary condition, the computationally obtained base flow is then compared with the experimental one in order to relate the control parameters from both methodologies. After validating the obtained agreement, we study the Hopf bifurcation that appears once the flow starts the vortex shedding through experimental and computational approaches. For the base flow derived from experimentally obtained snapshots, we monitor the evolution of the velocity amplitude oscillations. As to the computationally obtained base flow, its stability is analyzed by solving a global eigenvalue problem obtained from the linearized Navier–Stokes equations. Finally, the critical parameters obtained from both approaches are compared.