^{*}

^{*}

^{*}

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.

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. [

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.

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:

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.

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

where are the amplitude of the functions respectively.

Equations (1)-(3) become

where,.

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

where,.

The solution of Equation (8) has the form

where, and are the roots of the characteristic equation

The solution of Equation (8) is given by

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

where,

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

where,

On the surface at x = 0

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

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

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

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

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

Let us consider a 2-D simulation in which the pressure varies in the x and y directions.

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

the future to incorporate non-linear propagation effects.