Open access peer-reviewed chapter

Green’s Function Method for Electromagnetic and Acoustic Fields in Arbitrarily Inhomogeneous Media

Written By

Vladimir P. Dzyuba and Roman Romashko

Submitted: July 22nd, 2020 Reviewed: October 30th, 2020 Published: December 4th, 2020

DOI: 10.5772/intechopen.94852

Chapter metrics overview

340 Chapter Downloads

View Full Metrics

Abstract

An analytical method based on the Green’s function for describing the electromagnetic field, scalar-vector and phase characteristics of the acoustic field in a stationary isotropic and arbitrarily inhomogeneous medium is proposed. The method uses, in the case of an electromagnetic field, the wave equation proposed by the author for the electric vector of the electromagnetic field, which is valid for dielectric and magnetic inhomogeneous media with conductivity. In the case of an acoustic field, the author uses the wave equation proposed by the author for the particle velocity vector and the well-known equation for acoustic pressure in an inhomogeneous stationary medium. The approach used allows one to reduce the problem of solving differential wave equations in an arbitrarily inhomogeneous medium to the problem of taking an integral.

Keywords

  • inhomogeneous media
  • Green’s function
  • electromagnetic field
  • acoustic field
  • analytical method

1. Introduction

The chapter discusses the procedure for using the Green’s function for the analytical description of electromagnetic and acoustic fields in a stationary isotropic and arbitrarily inhomogeneous medium. In the case of the electromagnetic field, the wave equation for the electric vector of the electromagnetic field in the inhomogeneous medium with conductivity, dielectric and magnetic permeability is used. In the case of the acoustic field, the wave equation proposed by the author for the vector of particle velocity and the well-known equation for acoustic pressure in an inhomogeneous stationary medium are used. Using the Green’s function and the method of successive approximations makes it possible to achieve the required accuracy of calculating the electric and magnetic vectors of the electromagnetic field, as well as to calculate the vectors of complex intensity and intensity, density of energy, acoustic pressure and the particle velocity vector of the acoustic field in media with arbitrary spatial variability of the parameters. The approach used allows one to reduce the problem of solving differential wave equations in an arbitrarily inhomogeneous medium to integration. The chapter is divided into two parts. At the beginning of each part, the corresponding wave equations are derived and next, a method of using the Green’s function is described and analytical expressions describing the fields are formulated. At the beginning we will describe the method as applied to the electromagnetic field, and then as applied to the acoustic field.

Research and modeling of the electromagnetic field in spatially inhomogeneous natural and composite media are actively developing in various fields of science and technology, ranging from systems of underground and underwater electromagnetic communication to photonics, metamaterials and metasurfaces [1]. Such a wide field of scientific research requires methods of mathematical modeling of the properties of the electromagnetic field in media with different spatial scales of conductivity, magnetic and dielectric permittivity. At present, analytical methods are applicable in a very limited range of environments. Among the methods of mathematical modeling of the electromagnetic field in the frequency range from fractions of the hertz to optical, various numerical methods and technologies are used [2, 3]. Numerical modeling uses a variety of methods and technologies, for example, parallel computing which are used in electrodynamic modeling programs. Among them, there are also direct and universal methods for solving boundary problems. The drawback of these methods is a large expenditure of computer resources, which leads to a significant simplification of physical models of the environment and mathematical approximations. There is a third class of methods, in which, at the initial stage, analytical methods are used, for example, the Green’s function method, which brings the problem to a form that can be solved by fairly simple numerical methods. Below we will use exactly this approach using the Green’s function. Green’s function is actively used in a wide range of problems [4, 5, 6] of describing electromagnetic and other physical fields in various multilayer, chiral and anisotropic media, including inhomogeneous ones. The proposed procedure is also applicable to media with boundaries and arbitrary dependence on the coordinates of conductivity, magnetic and dielectric permittivity. The source of the field in the environment can be the electric current or an external field. The electric current can be located also inside the medium and outside it. The problem of descriptions the electric vector in an inhomogeneous medium by using the Green’s function is formulated as the integral equation with its subsequent solution by the method of successive approximations. This procedure uses the equation for the vector of electric field strength in an inhomogeneous medium, with a certain conductivity, magnetic permeability, and dielectric constant.

The acoustic energy flux density vector (intensity vector), basically, until the beginning of the second half of the 20th century was only of theoretical interest. The second half of the 20th century brought about reliable means of synchronous measurement, practically at a single point, of the acoustic pressure and the components of particle velocity vector necessary to determine the intensity vector of acoustic field [7, 8, 9, 10, 11, 12, 13, 14]. However, this did not lead to a significant increase in the number and quality of theoretical research methods and modeling of the intensity vector in an inhomogeneous medium. For the complete theoretical description of the acoustic field, knowledge of its acoustic pressure and the particle velocity vector is required. These two quantities make it possible to find the field of the acoustic intensity vector, to describe the energy and phase structure of the acoustic field. Knowledge of these quantities is useful for solving fundamental and applied problems of acoustic tomography and sounding of the geosphere, applied and fundamental hydroacoustics, creation of acoustic metamaterials, technical and architectural acoustics, noise control, etc. [15, 16, 17]. The acoustic pressure Part and the particle velocity vector Vrt are interrelated. This connection is obvious for a plane wave and, in the approximation of a continuous medium, has the following form: Vrt=1ρ0rPartdt , where ρ0r is the density of the medium unperturbed by the acoustic field at the point r and t is time, and is the Nabla operator. This relationship largely determined the development of the theory of sound as a scalar field of acoustic pressure. Currently, there are several directions for the development of methods of calculation and theoretical analysis of the characteristics of the intensity vector. In the first direction, the relationship between the acoustic pressure and the particle velocity vector is used. This approach is applicable when there are mathematical expressions for the acoustic pressure field. As a rule, this is only possible in a homogeneous medium or for simple waveguides [18]. The second direction requires the use of the continuity equation and the equation of state of the inhomogeneous medium, as well as dynamic equations of motion of elementary volumes or particles of the inhomogeneous medium, for example: the Euler or Navier-Stokes equations. These equations are viewed as a system of equations for determining the pressure and the particle velocity vector. This approach is used to model the propagation of waves in various environments, including plasma and stellar atmospheres [19, 20, 21, 22]. These equations are widely known, but to find analytical wave solutions of such systems given an arbitrary dependence of the density and speed of sound on the coordinates is a very difficult task. The use of the acoustic energy transfer equation is the third approach [9]. This approach allows one to describe the energy structure of the acoustic field which makes it possible to study the statistical characteristics of the complex intensity vector in a Gaussian delta-correlated inhomogeneous medium with refraction [9]. It is a very difficult task to find solutions to the transport equation in an inhomogeneous media. In turn, numerical methods for modeling metamaterials and propagation of acoustic waves in a medium are usually limited to specific problems [23, 24, 25]. None of the listed approaches, including numerical ones, provides the possibility of a complete theoretical description of the characteristics of the acoustic field and their evolution during field propagation in an arbitrary inhomogeneous medium. One of the promising directions is to use two wave equations in an inhomogeneous medium: equations for the acoustic pressure and equations for the particle velocity vector. We use this very approach. It is based on the proposed by authors wave equation for the particle velocity vector and the well-known equation for acoustic pressure in an inhomogeneous stationary medium. The proposed wave equation for the vector of the particle velocity of the acoustic field in a stationary inhomogeneous and isotropic medium is much more complicated than for the acoustic pressure. This makes it difficult to find the analytical solution for inhomogeneous media with an arbitrary spatial dependence of the density of the medium and the speed of sound in it. However, in an inhomogeneous medium, in which the field of the acoustic intensity vector is weakly vortex, the use of the Green’s tensor together with the method of successive approximations makes it possible to find analytical solutions for an arbitrary spatial dependence of the speed of sound and density of the inhomogeneous medium.

Advertisement

2. Electromagnetic field

By an inhomogeneous medium, we mean a medium in which the conductivity σr, dielectric εr and magnetic μr permittivity, and the current density Jr have an arbitrary, but differentiable, in the ordinary and in the generalized sense, dependence on coordinates points of the medium. Below we will not point out the explicit dependence of these and other quantities on time and coordinates, where this will not lead to misunderstanding. In an isotropic inhomogeneous medium, ε, μ, σ are scalar functions of coordinates. To derive the wave equation of the electromagnetic field in such a medium, we use the following well known fundamental and material Maxwell equations in a continuous isotropic and stationary medium:

1. ×H=J+D∂t;2. ×E=+B∂t;3. J=σE;4. J=∂tρf+ρext;5. J=Jf+Jext;6. D=εE;7. B=uH;8. ε=ε0εr;9. μ=μ0μr,E1

where ρf is the density of free charges of the medium, and ρext is the density of external charges introduced into the medium, εr and μr are the relative dielectric and magnetic permeability of the medium, and Jext is the current density created by free and external charges. The wave equation for the electric vector in an inhomogeneous medium can be obtained, as for a homogeneous medium, excluding the vector of magnetic field strength from the system of Maxwell’s equations. For this, we use the well-known vector analysis formulas [26] and take the rotor from the 2nd equation in system (Eq. (1)):

××E=EΔE=t×BE2

In this case, the source of electromagnetic field is the electric current density J, so the divergence of the vector E , we need to associate with a current density in the medium. For this, we use (Eq. (5)) of the Maxwell system of equations (Eq. (1)) and obtain E=1σJEσ. Using the vector analysis formulas, we find the following expression:

E= Eσ2σ2lnσEElnσlnσ××E++J1σ+1σJ]E3

The expression for the rotor of the magnetic field induction vector has the form

×B=×μH=μ×H+μ×HE4

Using equations 1 and 2 of the system of Maxwell equations and expressions (Eq. (2)) and (Eq. (3)), we find that t×B=μσtEμtJεμ2t2E+lnμ××E , and the desired equation for an electric vector with a field source, in which the charge flux from the volume occupied by the current is not equal to zero, has the following form:

εμ2t2E+μσtEΔElnμσ××E+Eσσ2lnσEElnσ=μtJJ1σ1σJE5

If there is no injection of external charges into the medium, the value J is equal to zero and J=Jext=tρext otherwise. When deriving (Eq. (5)), no conditions on the field frequency were used. Therefore, the equation is valid up to frequencies that correspond to wavelengths λ larger than the sizes of atoms or molecules. The smallness of the ratio of the first and second terms of equation (Eq. (5)) corresponds to the condition of quasi-stationarity of the electromagnetic field. For a monochromatic field with the angular frequency ω, the modulus of their ratio is equal εσω and small under conditions of high conductivity, low dielectric constant, or low the angular frequency. In this case, the propagation of the field in the medium will have a predominantly diffusion character and will be described by the following equation:

μσtEΔElnμσ××E+Eσσ2lnσEElnσ==μtJJ1σ1σJE6

When εσω>>1 the field propagation in the media is of the wave-type mainly. At the present time, there are no methods for finding exact solutions of equations of the type (Eq. (5)) and (Eq. (6)). Solutions satisfying a given accuracy can be obtained in two ways. The first is to use numerical methods. The second, which we will follow, consists in passing from the differential equation (5) to the integral equation using the tensor Green’s function of the Helmholtz equation for the Fourier - the spectrum of the vector of the electric field strength. The solution to an integral equation can be written in the form of a sequence of approximate solutions, in which each subsequent term is more accurate. It is important that such a procedure for finding a solution is applicable for arbitrary differentiable, both in the usual and in the generalized sense, dependences of σ, ε, and μ on coordinates. For this, we express the vector of the electric field Ert and the current density Jrt through their Fourier spectra Erω and Jrω

Ert=12π+Erωeiωt,Jrt=12π+Jrωeiωt.E7

For high frequencies, when the dependence of ε and μ from the field frequency cannot be neglected, but spatial dispersion and nonlinear effects can be neglected, Jrω=σωrEωr, Dωr=εωrEωr и Bωr=μωrHωr. Spatial dispersion plays a minor role in comparison with temporal dispersion and is significant in media with the mean free path of the charge or its diffusion much longer than the field wavelength. Below we will not indicate the dependence of the conductivity and permeability on frequency. Let’s consider conductivity and permittivity as a sum of a constant and a space-dependent variable:

σr=σc+σ1r,μr=μc+μ1r,εr=εc+ε1rE8

Let’s introduce the following notations:

ω2εrμriωμrσrσrσr2=k2rω,
fextωr=iωμrJrω+Jrωσ1rσ2r+1σrJrω,E9

fωr=μ1rμr+σ1rσr××Erωσ1rσrErωErωσ1rσr. Substituting expressions (Eq. (8)) and (Eq. (9)) into equation (Eq. (5)), we arrive at the following equation:

ΔErω+k2rωErω=frω+fextrωE10

Eq. (10) must be supplemented with boundary conditions. In an inhomogeneous medium, the interface between the media can be considered as an inhomogeneity with its dependence on coordinates, described by the corresponding functions, for example: Heaviside step function, etc. Therefore, the boundary conditions will be the conditions at infinity, where the field and its divergence must be equal to zero. In Eq. (10) the field source is not only the external currents (term fextωr but also the heterogeneity of the environment). These sources are described by fωr. At present, there are no methods for the analytical solution of equations similar to (10) with an arbitrary dependence of the term fωron coordinates. Nevertheless, using the Green’s functions of the vector Helmholtz equation in a homogeneous isotropic medium we can reformulate (10) into the integral with respect to the vector Erω, the solution of which can be found in an iterative way, for example, by the method of successive approximations.

Using (Eq. (8)) one can formulate k2rωas the sum of independent on coordinates kc2 function on coordinates and k12rω

k2r=ω2εcμciωμcσc+ω2εcμ1r+ε1rμc+ε1rμrσcμ1r+σ1rμc+σ1rμrσ1rσr2=kc2+k12rE11

Thus equation (Eq. (10)) can be written as:

ΔErω+kc2Erω=flωr+fextωr,E12

where f1ωr=fωrk12rωErω.

Formally, we can consider Eq. (12) as an inhomogeneous Helmholtz equation and using the Green’s function for it, we can rewrite Eq. (12) in the form of an integral equation. For vector fields, the Green’s function [7, 8, 9, 10] is a tensor of the second rank. In an orthogonal coordinate system, Eq. (5) decomposes into a system of three scalar equations for the projections of the field E (r, ω) on the coordinate axis. This simplifies the form of the Green’s tensor and it has only diagonal elements that are not equal to zero. It can be represented as the vector

Grr1=i=13niGirr1, where Girr1 is the components of which are the Green’s functions of the one-dimensional Helmholtz equation and ni are the unit vectors of the coordinate axes. Let the area Ω in which we describe the field be large enough so that on its borders the field and its derivatives can be equated to zero. Using the Green tensor, we can rewrite Eq. (5) for the electric vector at the point rΩ of the in the form of the following integral equation

Erω=i=13niΩGirr1fiextωr1+fi1ωr1dr1E13

where rr1Ω, fiextωr1 and fi1ωr1 are the projections of vectors fextωr and f1ωr on to the coordinate axes. Integration is performed over the volume occupied by inhomogeneities, which are secondary sources of the field. In practice, the volume should be chosen such that secondary and higher order sources make a noticeable contribution to the field. Due to the rapid decrease in the amplitude of the Green’s function and, especially with a strong absorption of the electromagnetic field by the medium, the region of integration can be about 1/α where α is the absorption coefficient of the field.

The steps for finding Erω by the method of successive approximations can be as follows. We find the zeroth approximation E0rω for the field, which is valid in a homogeneous medium with parameters σc, μc, εc

E0rω=i13niΩ1Grr1nifextωr1dr1.E14

Integration is performed over the volume Ω1 occupied by the external current (primary source of the field). This solution describes the primary field created by an external current. Using obtained by Eq. (14) expression E0rω and expression (9), we find fi1ωr1. Using (Eq. (13)) and integrating, we obtain a more accurate first E1rω approximation for Erω, which takes into account the influence of medium inhomogeneities on the field. To find the second approximation, it is necessary to substitute E1rω into fi1ωr1 and using (Eq. (13)) to obtain the second more accurate approximation. Similarly, more accurate solutions are obtained that take into account multiple field scattering by medium inhomogeneities. At these stages, the integration is performed over the volume occupied by inhomogeneities, which are secondary sources of the field. If the source of the field in an inhomogeneous medium is an external field with an electric vector Eextrω it should be used as the vector E0rω.

For determining the magnetic field component one uses Maxwell’s equations and writes them in terms of magnetic and electric fields Fourier spectrums: ×Erω=ωBrω. Substituting (Eq. (13)) in this equation one obtains

Hrω=i1ωμrΩi=13ni×Grr1fiextωr1+fi1ωr1dr1.E14a

Using the Green’s function, as the experience of its use shows (e.g. [7]) in such tasks, significantly reduces the requirements for computing resources and reduces the computation time. Note that the proposed procedure can be effective in simulating the optical properties of metamaterials, nanocomposites, and nanostructures.

Advertisement

3. Acoustic field

The wave equation for acoustic pressure Part in a continuous inhomogeneous motionless and stationary medium is well known [16, 17]

1c2r2t2Part+ΔPart+Part+frtlnρ0r=frt.E15

where frt is the density of volumetric external forces that are the source of the acoustic field.

Eq. (1) is obtained by excluding the particle velocity vector from the linearized Euler equations, continuity and state of the medium. If we exclude acoustic pressure from these equations, then we get the equation for the vector of the particle velocity of the acoustic field. For this, we differentiate the equation of state in taking into account the smallness of the acoustic pressure, perturbation of the density of the medium by the ρart in comparison with the background values ρ0r and P0r the medium. In the inhomogeneous medium, the equation of state Pcρrt describes the relationship of the instantaneous local value of pressure and density of the medium. Therefore, it is necessary to use the total time derivative when differentiating the equation of state. Using it, we find in the linear approximation

ddtPcρrt=c2rtρart+ρ0rVrt,E16

where Cr=Pcρrtρrt is the local phase speed of sound for the acoustic pressure wave. We used the expansion ρrt=ρ0rt+ρart and the condition ρartVrt<<ρ0rVrt. With the help of expression (2) and representation, Pcρrt=P0r+Part the equation of continuitytρrt+ρrtVrt=0 is reduced to a linearized form

1ρ0rc2rtPrt+Vrt=0E17

In the inhomogeneity of the medium Vrt0 , even in the approximation of an incompressible medium. Let us Take the time derivative on the linearized Euler equation ρ0rtVrt+Prt+frt=0 and take the gradient of the equation of continuity (Eq. (17)). We exclude the acoustic pressure from the obtained expressions and find the equation for the particle velocity vector

1c2r2t2VrtΔVrtlnρ0rc2rVrt××Vrt==1ρ0rc2r∂tfrtE18

When frt=0 then ××Vrt=×lnρ0r×Vrt and can transformed equation (Eq. (18)) to the following form, which is valid in the absence of external forces

1c2r2t2VrtΔVrtlnρ0rc2rVrt++×lnρ0r×Vrt=0E19

Eqs. (18) and (19) are much more complicated than equation (Eq. (15)) due to the third and fourth vortex term. The estimate of their ratio is

lnρ0rc2rVrt×lnρ0r×Vrtlnρ0rc2rlnρ0rE20

In areas of the medium where this ratio is greater than unity, the fourth term in equations (Eqs. (18) and (19)) can be neglected. As a rule, these are media with a large relative gradient of the speed of sound. One of the examples of such media can be the marine environment, in which the local gradient of the speed of sound is determined less by the change in water density than salinity and temperature [12, 24]. Directly near the surface and the ocean floor or the interface between the media, the relative density gradient of the medium can be large. In these regions, the field of the particle velocity vector and acoustic intensity can have a significant rotational (vortex) component.

At present, the solution of these equations is possible only by numerical methods. If the fourth term in the equations is small, the equations for the vector of particle velocity and acoustic pressure allows one to find analytical expressions connecting the phases and moduli of vector of complex intensity and particle velocity vector, pressure, density of acoustic energy with the density of the medium and the speed of sound in it. Let us do it for equation 19. Using both scalar function Ψrt and vector Urt

Prt=ZprZp0ΨrtVrt=ZvrZv0Urt,E21

we rewrite equations (Eqs. (15) and (19)) in the following form:

1c2r2t2ΨrtΔΨrt+3ZprZpr4Zp2rΔZpr2ZprΨrt=0E22
1c2r2t2UrtΔUrt+3ZvrZvr4Zv2rΔZvr2ZvrUrt=0,E23

where Zpr=ρ0r and Zvr=1ρ0rc2r, and Zp0=ρ0r0, Zv0=1ρ0r0c2r0 are values Zpr and Zvr at some point in space r0. For the spectral components Ψrω and Urω using the Fourier transform of equations (Eqs. (22) and (23)) with respect to the time variable, we obtain the following equations

ΔΨrω+kψ2rΨrω=0,E24
ΔUrω+kU2rUrω=0,E25

where kψ2r=ω2c2r34ρ0rρ0r2+Δρ0r2ρ0r

kU2r=ω2c2r+54ρ0rρ0r2+ρ0rcrρ0rcr+3crcr2Δρ0rρ0r2ΔcrcrE26

From expressions (Eq. (26)) it follows that the gradient of the speed of sound affects only the vibrational speed in a medium with a small swirl of the particle velocity field. The acoustic pressure depends only on the gradient of the density of the medium and does not depend on the gradient of the speed of sound. This situation is valid for media in which the density gradient of the medium is less than the gradient of the speed of sound. These differences form the phase difference between the acoustic pressure and the vibrational velocity vector during the propagation of an acoustic wave in an inhomogeneous medium. Different field reactions Vrt and Prt to the density gradient of the medium and the gradient of the sound speed form the phase difference between the acoustic pressure Φprt and the particle velocity Φvrt vector during the propagation of an acoustic wave in an inhomogeneous medium. Solutions to equations (Eqs. (24) and (25)) can be found using the method of successive approximations. For this, we represent these equations in the following form:

ΔΨrω+k02Ψrω=k1Ψ2rΨrωE27
ΔUrω+k02Urω=k1v2rUrωE28

where k02=ω2c2r0, k1Ψ2r=k02kψ2r and k1U2r=k02kU2r

Similarly to the case of an electromagnetic field in the inhomogeneous medium, using the scalar Grr1 and vector Grr1 Green’s functions of the Helmholtz equation for a homogeneous unbounded medium. Eqs. (2) and (13) and Eq. (2.14) can be rewritten in the form of the following integral equations:

Ψrω=Ψ0rω+ΩGrr1k1Ψ2r1Ψr1ωdr1E29
Uirω=U0irω+ΩGirr1k1U2r1Uir1ωdr1E30

Here Uir1ω is the projections of the vector onto the coordinate axes and Ψ0rω and U0irω are the solutions of equations (Eq. (27)) and (Eq. (28)) with the right-hand side equal to zero, and Girr1 are the components of the vector Green’s function. The steps for finding a solution to equations (2.15) and (2.16) by the method of successive refinements can be as follows:

  1. We find Ψ0rω and U0irω which are the zeroth approximation for the field, valid in a homogeneous medium

  2. We find an explicit form of dependence k1Ψ2r=k02kψ2r and k1U2r=k02kU2r on coordinates

  3. Using Ψ0rω and U0irω and expressions for k1Ψ2r and k1U2rwith the help of (Eqs. (29) and (30)), we obtain more accurate first approximations that take into account single scattering of the primary field.

  4. To obtain the second more accurate approximation, it is necessary to substitute the first approximations in expressions (Eqs. (29) and (30)) and obtain the second more accurate approximation of the solution.

Similarly, you can get solutions that are more accurate. Integration is performed over the volume of the inhomogeneous medium, the inhomogeneities of which will be secondary, etc. field sources. In a real situation, the volume should be chosen such that its secondary sources make a noticeable contribution to the field.

Let us consider an example that shows how the parameters of an inhomogeneous medium affect the characteristics of the acoustic field. We represent the acoustic pressure, and the vector of the particle velocity of the monochromatic acoustic field by the frequency ω in the following form

Prt=P0rexpiωtΦpr and Vrt=V0rexpiωtΦvr. Complex intensity vector will be written as Ir=PrtVrt=I0rexpiΦvrΦpr. In a medium without absorption of the acoustic energy, the phases Φpr and Φvr, respectively are equal to the phases Ψrt and Urt. The wave vector of a wave is normal to its phase surface and is determined by the wave phase gradient and the wave number by the modulus of this gradient. For the wavenumbers of the pressure kpr and the particle velocity vector kvr, we can take, respectively, the quantities kψr and kUr if the inequalities Δψ0rψ0r<<kψ2rΦpr2 and ΔU0rU0r<<kU2rΦvr2. The refractive indices of the medium for the acoustic pressure and the particle velocity relative to the point r0 are different and accordingly, equal npr=kprk0 and nvr=kvrk0, where k0=ωcr0 and Cr0 is the phase velocity of sound for the acoustic pressure wave at the point r0. The phase velocities of the acoustic pressure wave and the particle velocity vector become different, which leads to the inequality of the phases of the acoustic pressure and the particle velocity vector when the acoustic wave propagates in an inhomogeneous medium. The absorption of acoustic energy by the medium is taken into account by assuming the speed of sound and the density of the medium to be complex quantities, in which the imaginary part is responsible for the absorption of the energy of the acoustic field. In this case, the phases and acquire an additive equal to the phases of the values ZprZp0 and ZvrZv0. To avoid cumbersome expressions, we restrict ourselves to the first approximation. The proposed example is often implemented in real measurements of the characteristics of the acoustic field. In practice, as a rule, the projections of the vibrational velocity vector and, accordingly, the intensity vector are measured in an orthogonal coordinate system, for example, in a Cartesian one. Consider the projection of a vector Urω on the OX axis. In this case, the projection will be a function of only the x coordinate, and the Y and Z coordinates will act as parameters and determine the straight line parallel to the OX axis, along which the observation point x changes. The wave numbers k1Ψ2r and k1U2r, accordingly, the solutions of equations (Eqs. (29) and (30)) depend on the values of these parameters. In fact, we turn to the case of one-dimensional propagation of acoustic radiation along the OX axis passing through the point r0X0Y0Z0. Let us choose Ψ0rω and U0rω both in the form of plane waves and propagating along the X axis. We can put the moduli of these plane waves Ψ0ω and U0ω equal to the moduli of the acoustic pressure P0 and the component vx of the particle velocity vector on the OX axis. In this case, the component of the vector Green’s function will be equal to the one-dimensional Green’s function

Gxx1=12ik0expik0xx1 we find the solution corresponding to the first approximation at the point X

Ψxy0z0ω=P0expikx0+xk1Ψ2x1y0z0P02ik0expik0xdx1++xk1Ψ2x1y0z0P02ik0expik0x+2ik0x1dx1=P0expikx0+Ψ1x+Ψ2xE31
Uxxy0z0ω=Vxexpik0x+xk1U2xy0z0Vxexpik0x2ik0dx1++xk1U2x1y0z0Vx2ik0expik0x+2ik0x1dx1=Vxexpik0x+U1x+U2xE32

Here P0expik0x и U0expik0x represent the primary radiation, the second terms are the radiation scattered forward in the region x1x, and the third terms are the radiation scattered back. Solutions (Eqs. (31) and (32)) and the relations (Eq. (21)) allow us, in the first approximation, to find an expression for the projection of the complex intensity vector on the on the OX axis

Ixxy0z0ω=Zpxy0z0Zvxy0z0Zp0Zv0P0Vx++P0U1x+P0U2x+Ψ1xVx++Ψ1xU1x+Ψ1xU2x+Ψ2xVx++Ψ2xU1x+Ψ2xU2xE33

In this expression, the first term describes the complex intensity vector of the primary radiation, the fifth term corresponds to the forward propagating secondary radiation, and the ninth term corresponds to the backscattered radiation. The other terms describe the mutual energy of the primary and scattered radiation. If the field is measured arriving at a point x only from the region, x1x then the dependence of the projection of the complex intensity vector on the OX axis takes the following form:

Ixxy0z0ωexpiΦxy0z0=C0x0y0z0Cxy0z0P0Vx1+i2k0αvx0y0z0αpxy0z0++14k02αvxy0z0αpxy0z0.E34

In this expression αpx=xk1Ψ2x1dx1 and αvx=xk1U2x1dx1. The modulus I0xy0z0ω and phase Фxy0z0ω of the complex of acoustic intensity vector are respectively equal:

I0xω=C0P0Vx4k02Cxy0z04k02+αpxy0z024k02+αvxy0z02,Φxy0z0ω=arctg2k0αvxy0z0αpxy0z04k02+αvxy0z0αpxy0z0.E35

The field intensity vector is

Re12Ixxy0z0ω=12C0Cxy0z0P0Vx1+14k02αvxy0z0αpxy0z0,E36

and for the average field energy density we have the following expression:

εxy0z0ω=Pxy0z0ωPxy0z0ωρxy0z0C2xy0z0==P02ρ0C2xy0z01+14k02αpxy0z0αpxy0z0.E37

From expressions (Eqs. (22) and (23)) it is seen that the inhomogeneous nature of the speed of sound will have a more significant effect on the particle velocity vector than on the acoustic pressure. This makes it possible in principle to create methods for separately measuring the contribution to the acoustic field in an inhomogeneous medium of the density of the medium and the speed of sound in it.

In conclusion, we note that the proposed method makes it possible to analytically and numerically solve the problems of mathematical modeling of a shallow sea, remote sensing of natural media, problems of acoustics of a shallow sea, modeling acoustic and optical metamaterials, etc. Note that for applied problems of acoustics, both fields of the particle velocity vector and the intensity vector in any inhomogeneous medium have a vortex character. Therefore, the algorithms for solving applied problems of ocean and especially shallow sea acoustics, problems of modeling the propagation of acoustic energy in composite media and metamaterials should take into account the vortex component of the vector acoustic field intensity and curvature of the streamlines of the acoustic field.

Advertisement

Acknowledgments

The research is supported by the grant from the Russian Science Foundation (project # 19-12-00323).

References

  1. 1. Caloz Ch., Itoh T. Electromagnetic metamaterials: transmission line theory and microwave applications (the engineering approach). A John Wiley & Sons, Inc., 2006. 352 p
  2. 2. Prather, D. W., S.Shi. Formulation and application of the finite-difference time-domain method for the analysis of axially symmetric diffractive. J. Opt. Soc. Am. A.1999. Vol.16, Iss. 5. P. 1131-1142. Doi:10.1364/JOSAA.16.001131
  3. 3. Taflove A. and Hagness S. C. Computational Electrodynamics: The Finite-Difference Time-Domain Method. Publisher: Artech House; 3 edition ,2005. 1038 pages
  4. 4. C.T. Tai, Diadic Green, Function in Electromagnetic Teory, 2nd ed. Piscataway, NJ, USA IEEEPress,1994
  5. 5. George W. Hanson Dyadic Green’s Function for a Multilayered Planar Medium—A Dyadic Eigenfunction Approach. IEEE Transactions on Antennas and Propagation 2005, 52(12), P.3350 – 3356. DOI:10.1109/TAP.2004.836409
  6. 6. Dzyuba, V.P., Zapolsky, A.M. The G,reen function for the elastic half-space with a curvilinear boundary. Doklady Earth Sciences V. 372, May 2000, Pages 709-711
  7. 7. V.A. Gordienko Vector-phase methods in acoustics. Moscow: Fizmatlit; 2007. 480 p. ISBN978-5-92210864
  8. 8. F. J. Fahy Sound intensity. 2nd ed. London: Elsevier; 1989
  9. 9. Dzyuba V.P. Scalar-vector methods of theoretical acoustics. Vladivostok: Dalnauka; 2006. 192 p. ISBN 5-8044-0559-4
  10. 10. Akulichev V.A., Dzyuba V.P., Gladkov P. V. On Acoustic Tomography Scheme for Reconstruction of Hydrophysical Parameters for Marine Environment. Theoretical and Computational Acoustics. Ed. Er- Chang Shang and Qihi Li. Published by World Scientific; 2002.P.107 -114
  11. 11. Baster K.J. , Lanchle L.C., Mc. ConnelJ.A. Development of a Velocity Gradient Underwater Acoustic Intensity Sensor. JASA .1999.Vol. 106. P. 3178-3188. doi:10.1121/1.428172
  12. 12. Mann J.A., Tichy J., Romano A.J. Instantaneous and time-average energy transfer in acoustic fields. JASA. 1987.Vol. 82. №1. P. 17-30. DOI:10.1121/1.395562
  13. 13. K. J. Taylor Absolute measurement of acoustic particle velocity. The Journal of the Acoustical Society of America. 1976 Vol. 59, 691 DOI:10.1121/1.380896
  14. 14. Nelson, P.A. and Elliott, S.J. (1991) Active control of sound: Academic Press; 1993. 250pp. ISBN: 0125154259
  15. 15. Medwin H.. Sounds in the sea: from ocean acoustics to acoustical oceanography. Cambridge, UK: Cambridge Univ. Press. 2005. 643 pp
  16. 16. David R. Dall’Ostoa and Peter H. Dahl Properties of the acoustic intensity vector field in a shallow water waveguide. The Journal of the Acoustical Society of America 2012. 131, 2023; DOI: 10.1121/1.3682063
  17. 17. Woon Siong Gan New Acoustics Based on Metamaterials . Springer Nature Singapore Pte Ltd. 2018. P. 314. DOI:10.1007/978-981-10-6376-3
  18. 18. Morse, P.N., Ingard, K.U. Theoretical Acoustics. New York: McGraw-Hill; 1987. 949P
  19. 19. Jan-Niklas Hau1, and Björn Müller Acoustic wave propagation in a temporal evolving shear-layer for low-Mach number perturbations. Physics of Fluids 2018; 30, 016105. DOI:10.1063/1.4999044
  20. 20. L. Friedland, G. Marcus, J. S. Wurtele, and P. Michel3, Excitation and control of large amplitude standing ion acoustic waves. Phys. Plasmas 2019; 26, 092109 DOI: 10.1063/1.5122300
  21. 21. Zhu, X., Li, K., Zhang, P. et al. Implementation of dispersion-free slow acoustic wave propagation and phase engineering with helical-structured metamaterials. Nat Commun. 2016; 11731. DOI:10.1038/ncomms11731
  22. 22. N. Anders Petersson & Björn Sjögreen High Order Accurate Finite Difference Modeling of Seismo-Acoustic Wave Propagation in a Moving Atmosphere and a Heterogeneous Earth Model Coupled Across a Realistic Topography.Journal of Scientific Computing 2018; volume 74, pages 290–323. DOI:10.1016/j.jcp.2020.109386
  23. 23. S.MishrabCh.SchwabaJ.ŠukyscMulti- level Monte Carlo finite volume methods for uncertainty quantification of acoustic wave propagation in random heterogeneous layered medium. Journal of Computational Physics. Volume 312, 1 May 2016, Pages 192-217. DOI:10.1016/j.jcp.2016.02.014
  24. 24. Elias Perras1, and Chuanzeng Zhang1 Analysis of acoustic wave propagation in composite laminates via aspectral element method// PAMM • Proc. Appl. Math. Mech. 2019;19:e201900282. DOI:10.1002/pamm.201900282
  25. 25. Yijun Mao, Jiancheng Cai, Yuanyuan Gu and Datong Qi Direct Evaluation of Acoustic-Intensity Vector Field Around an Impedance Scattering Body AIAA JOURNA. 2015; Vol. 53, No. 5, May. DOI:10.2514/1.J053431
  26. 26. Granino Arthur Korn, Theresa M. Korn Mathematical handbook for scientists and engineers, McGraw-Hill, 1968. P. 1130

Written By

Vladimir P. Dzyuba and Roman Romashko

Submitted: July 22nd, 2020 Reviewed: October 30th, 2020 Published: December 4th, 2020