﻿Analytical Solution for Acoustic Waves Propagation in Fluids

World Journal of Mechanics
Vol.1 No.5(2011), Article ID:8162,4 pages DOI:10.4236/wjm.2011.15030

Analytical Solution for Acoustic Waves Propagation in Fluids

Mohamed I. A. Othman1, Mohamed Galal Sayed Ali2, Roushdi Mohamed Farouk3

1,3Department of Mathematics, Faculty of Science, Zagazig University, Zagazig, Egypt

2Department of Physics, Faculty of Science, Minia University, Minia, Egypt

E-mail:{m_i_othman, mgalal09, rmfarouk1}@yahoo.com

Received July 12, 2011; revised August 20, 2011; accepted September 2, 2011

Keywords: Acoustic waves, normal mode analysis, fluid medium

ABSTRACT

This paper presents a mathematical model of linear acoustic wave propagation in fluids. The benefits of a mathematical model over a normal mode analysis are first discussed, then the mathematical model for acoustic propagation in the test medium is developed using computer simulations. The approach is based on a analytical solution to the homogeneous wave equation for fluid medium. A good agreement between the computational presented results with published data.

1. Introduction

In recent years, physical acoustic wave modeling has become a successful tool in diagnostic and therapeutic ultrasound application. There are several wave equations available for describing acoustic wave propagation [1-4]. Numerical methods can be used as a tool for sound field simulation. Discrete-time simulation algorithms for wave propagation can be derived by numerically solving a acoustic wave equation in terms of the variables for sound pressure and particle velocity. Initial conditions for time derivatives and boundary conditions for space derivatives are necessary to provide a complete set of solutions of the wave equation. These equations are most commonly solved by propagation in time. However, when propagating over large distances, such methods are expensive in terms of memory and computational costs .

The normal mode method analysis gives exact solutions without any assumed restrictions on pressure and velocity components distributions. It is applied to wide range of problems in different branches (Othman [6-8], Sharma et al. , Othman and Kumar , Othman and Singh  and Othman et al. ). It can be applied to boundary-layer problems, which are described by the linearized Navier-stokes equations in electrohydrodynamic (Othman ).

In this paper, the normal mode analysis can be employed to solve linear acoustic wave equation analytically. The technique focuses on description of a linear model and discuses the conditions under which using this technique. The propagation of acoustic pressure wave by the normal mood analysis in a medium with two-dimensional spatially-variable acoustic properties has been explained.

2. Acoustic Wave Equation

Consider sound waves propagating in the water. Instead of the wave equation, we base our work on the basic Euler’s equation and the equation of continuity. For simplicity, the discussion is confined to a two-dimensional space. In a 2-D Cartesian coordinate system, the sound pressure and the particle velocity v satisfy the following linear equations: , (1) , (2) . (3)

where is the particle velocity, p(x, y, t) is the pressure and is the density of the fluid with wave number where , is the angular frequency, c and are the speed of sound and attenuation in inhomogeneous medium, respectively.

3. Normal Mode Analysis

The solution of considered physical variable can be decomposed in terms of normal modes as the following form (4)

where are the amplitude of the functions respectively.

Equations (1)-(3) become (5) , (6) . (7)

where, .

Equations (5)-(7) form a coupled system Eliminating and between Equations (5)-(7) we obtain (8)

where, .

The solution of Equation (8) has the form . (9)

where , and are the roots of the characteristic equation (10)

The solution of Equation (8) is given by (11)

From Equations (6) and (11) we can obtain (12)

where, (13)

From Equations (7) and (11) we can obtain (14)

where, (15)

4. Boundary Conditions

On the surface at x = 0 (16)

Substituting from (4) into (16) then Equations (11) and (14) (17) (18)

By adding Equations (17) and (18) we obtain (19)

By subtracting Equations (17) and (18) we get (20)

By substituting from Equations (19) and (20) into Equations (11), (12) and (14) (21) (22) (23)

5. Computational Results

To study the wave propagation phenomenon in viscous medium and under different frequencies, we can apply the theoretical acoustic viscous wave Equation (21). Using water as the medium, the parameters are given as following: and . Let the wave peak amplitude be Pa and m/sec at the source (x = 0), we simulate the pressure wave peak amplitude, in Equation (21), vs. the distance from the source at various frequencies. The results are shown in Figure 1 (in dB and linear scale). As expected, the peak wave amplitude becomes smaller as we move further from the source. We can also notice that as the wave frequency goes higher, more attenuation can be observed at any given location. The peak pressure as function of frequency shown in Figure 2 for fixed distance at x = 2 cm. It is clear from this figure that the magnitude of peak of pressure little changes with the frequency. The predicted results are very agreement to the recorded by Wang .

Let us consider a 2-D simulation in which the pressure varies in the x and y directions. Figure 3 shows the 2-D pressure computational as function in the plane x-y. From this figure, it can be seen that the pressure amplitude becomes smaller when moving in the x-y plane further from the source.

6. Conclusions

A normal mode analysis which accurately the pressure acoustic wave equation, has been developed. This technique has a number of attractive features, foremost of which is the speed and simplicity with which it can be designed and implemented. The model could be used in Figure 1. Pressure amplitude (dB) as function of the distance from the origin. Figure 2. The pressure amplitude as function of frequency at distance 2 cm from the source. Figure 3. The pressure amplitude as function of plane x-y direction.

the future to incorporate non-linear propagation effects.

REFERENCES

1. L. E. Kinsler, A. R. Frey, A. B. Coppens and J. V. Sanders, “Fundamentals of Acoustics,” 4th Edition, Wiley, New York, 2000.
2. S. W. Rienstra and A. Hirschberg, “An Introduction to Acoustics,” Eindhoven University of Technology, Eindhoven, 2009.
3. S. W. Rienstra, “Sound Propagation in Slowly Varying Lined Flow Ducts of Arbitrary Cross Section,” Journal of Fluid Mechanics, Vol. 495, 2003, pp. 157-173. doi:10.1017/S0022112003006050
4. D. Daniel and Joseph, “Review Potential Flow of Viscous Fluids: Historical Notes,” International Journal of Multiphase Flow, Vol. 32, No. 3, 2006, pp. 285-310. doi:10.1016/j.ijmultiphaseflow.2005.09.004
5. J. G. Maloney and K. E. Cummings, “Adaptation of FDTD Techniques to Acoustic Modeling,” 11th Annual Review of Progress in Applied Computational Electromagnetics, Vol. 2, 1995, pp. 724-731.
6. M. I. A. Othman, “Effect of Rotation in Case of 2-D Problems of Generalized Thermoelasticity with Thermal Relaxation,” Mechanics & Mechanical Engineering, Vol. 8, 2005, pp. 111-126.
7. M. I. A. Othman, “Effect of Rotation on Plane Waves in Generalized Thermo-Elasticity with Two Relaxation Times,” International Journal of Solids and Structures, Vol. 41, No. 11-12, 2004, pp. 2939-2956. doi:10.1016/j.ijsolstr.2004.01.009
8. M. I. A. Othman, “Lord-Shulman Theory under the Dependence of the Modulus of Elasticity on the Reference Temperature in Two-Dimensional Generalized Thermoelasticity,” Journal of Thermal Stresses, Vol. 25, No. 11, 2002, pp. 1027-1045. doi:10.1080/01495730290074621
9. J. N. Sharma, R. Chand and M. I. A. Othman, “On the Propagation of Lamb Waves in Visco-Thermoelastic Plates under Fluid Loadings,” International Journal of Engineering Science, Vol. 47, No. 3, 2009, pp. 391-404. doi:10.1016/j.ijengsci.2008.10.008
10. M. I. A. Othman and R. Kumar, “Reflection of Magneto-Thermoelastic Waves under the Effect of Temperature Dependent Properties in Generalized Thermo-Elasticity with Four Theories,” International Communications in Heat and Mass Transfer, Vol. 36, No. 5, 2009, pp. 513-520. doi:10.1016/j.icheatmasstransfer.2009.02.002
11. M. I. A. Othman and B. Singh, “The Effect of Rotation on Generalized Micropolar Thermoelasticity for a Halfspace under Five Theories,” International Journal of Solids and Structures, Vol. 44, No. 9, 2007, pp. 2748- 2762. doi:10.1016/j.ijsolstr.2006.08.016
12. M. I. A. Othman, Kh. Lotfy and R. M. Farouk, “Generalized Thermo-Micro-Stretch Elastic Medium with Temperature Dependent Properties for Different Theories,” Engineering Analysis of Boundary Element, Vol. 34, No. 3, 2010, pp. 229-237. doi:10.1016/j.enganabound.2009.10.003
13. M. I. A. Othman, “Electrohydrodynamic Stability in a Horizontal Viscoelastic Fluid Layer in the Presence of a Vertical Temperature Gradient,” International Journal of Engineering Science, Vol. 39, No. 11, 2001, pp. 1217- 1232. doi:10.1016/S0020-7225(00)00092-6
14. S. Wang, “Finite-Difference Time-Domain Approach to Underwater Acoustic Scattering problems,” Journal of the Acoustical Society of America, Vol. 99, No. 4, 1996, pp. 1924-1931. doi:10.1121/1.415375