Modeling of ultrasonic waves propagation in inhomogeneous medium using fibered spaces method (k-space)

19/05/2014 4:03pm

Автор: V. S. Morkun, N. V Morkun, A.V. Pikilnyak

Категории: automatization

Abstract:The method for more accurately parameters estimation of the ultrasonic waves propagating in random heterogeneous media, including solid, liquid and gas phases is described. 


Key words:phased array, ultrasound, pulp, control.

Vladimir Morkun
Vice-Rector for research, Doctor of Science, professor of Computer Science, Automation and Control Systems department
Krivyi Rih National University

Natalia Morkun
PhD, Associate professor of Economic Cybernetics and Project Management Department
Krivyi Rih National University


Andrey Pikilnyak
PhD- student of Computer systems and networks department
Research Assistant of the Computer Science, Automation and Control Systems department

1. Introduction. For the control of the basic technological parameters and mineral beneficiation process control, an important task is to control the parameters of complex heterogeneous mediums, including solid, liquid and gas phases.
The basic relations describing the ultrasonic oscillations and waves in the medium, follow from the equation of medium state, Newtonian equations of motion and the continuity equation [1-3]. The result are the wave-type equations that can be solved with appropriate initial and boundary conditions.
Let’s introduce differential characteristic  to describe the ultrasonic waves radiation field. By we mean the intensity of the ultrasonic wave (with wavelength ), which is defined as radiation power per solid angle unit which passing through a unit area perpendicular to the direction  of the point . Here  is the unit vector defining the direction in space,  is the radius vector defining the position of a given point in space.
The purpose is to study the volume ultrasonic waves propagation in the gas-containing iron ore slurry, defining of ultrasonic field characteristics in a liquid medium containing solid particles and gas bubbles, influence laws of the suspended in a liquid particles fluctuations on the performance of the ultrasonic field.
The presence of solid particles and gas bubbles introduces some features to the process of ultrasonic wave energy attenuation and scattering [4].
Wave scattering on the solid phase particles becomes significant when the wavelength  is commensurate with the size of the particles. Consequently, the total intensity of the ultrasonic wave at a given point equal to the sum of the intensities of waves coming from all scattering centers. Scattering cross sections in this case are additive, so the linear absorption and scattering coefficients can be determined by the formulas
,                                            (1)
where n - is the particle concentration;  and s(l) - total cross sections of the acoustic wave absorption and scattering on the particle.
Total cross sections of absorption and scattering depends not only on the wavelength of the ultrasonic vibration, but also on the particle size r. The main characteristic of the ultrasonic radiation field  must be determined from the kinematic equation. Before write this equation, we introduce the concept of the differential at the corners scattering coefficient
,                         (2)
where  - is differential at the corners energy scattering cross section on the solid phase particle.
The value  represents a part of the energy which scattering by particles in the element of solid angle . Obviously, the total scattering cross section ss associated with differential scattering cross section ratio
,                                (3)
The kinetic equation which is solved by function  can be obtained by considering the energy balance in a volume element of the phase space
where .  - is the ultrasound source radiation density function, which determines the average amount of energy emitted per unit time single phase volume.
Under the phase coordinates means the totality of variables r and W, while elementary phase volume is determined by the product .
Equation (4) can be reduced to an integral equation of the form (5)
where ,  is the Dirac delta function;  - is the free term of the integral equation (5), which determines the intensity of the unscattered ultrasonic wave; .
The solution of equation (5) can be written in the form of a Neumann series [5], which is the expansion of the ultrasonic waves scattering multiplicity solution.
However, to obtain an expression analytically even for the singly scattered radiation is impossible. Therefore it is necessary to apply numerical methods for solving integral equations of the form (5).
Let's consider the method of fiber spaces (k-space) for modeling of ultrasonic wave propagation in inhomogeneous medium using coarse grids, with maintaining the required accuracy [6 -8].
We describe the ultrasonic waves propagation depending on the mass conservation equations, momentum conservation law and the equation of state using the first order dual equations, which can be summarized as follows [8, 9]
,                      (6)
,                               (7)
where - the time and space dependent ultrasound pressure perturbations (x - 3D Cartesian axis (x, y, z)); - is the spatially dependent density; - is the spatial dependent sound speed; - is the velocity of the particle and - is the absorption coefficient which equivalent to the inverse of the relaxation time.
Let's represent the all absorption effects with one relaxation time. From (7), the simplified equation can be written as follows

We differentiate (6) with respect to time and variations in (7), and the final equation can be represented as follows
,            (6’)
,               (7’)
Taking into account the permutations (6')

By substituting this equation in (7'), we obtain
,            (8)
The simplification of the pressure deviation to the density gradient can be represented as follows
,                           (9)
Taking into account (8), eq. (9) can be represented as follows
,    (10)
This is a linear wave equation of ultrasonic wave propagation in the heterogeneous medium with the absorption parameters.
Let's simplify (10) by separating the parameters of the sound velocity  and density  from the second derivatives of pressure taking into account the spatial and temporal variables to solve the problem of ultrasoundpropagation using the fiber space method.
The original equation can be written in the form
,              (11)
The normalized pressure can be represented as follows

By substituting this equation in (11) we obtain

After simplifying

Taking into account further simplifications the equation takes the form

Even more simplification can be obtained by determining the functions q(r, t) and v (r, t)efficient sources, which can be summarized as follows

By simplifying (11) we obtain
,                 (12)
This equation can be easily transformed into the frequency domain by using the three-dimensional spatial Fourier transform as follows
,                            (13)
where F(k,t), Q(k,t) и V(k,t) – three-dimensional spatial Fourier transformation of values ,  and  respectively. Equation (13) satisfies the total wavefield, and is defined as the sum of the incident and scattered field , and the scattered wave field.

For the case of an inhomogeneous medium, we introduce an additional source  and by substituting it into (13) we obtain the following expression
,               (14)

where F - is a spatial Fourier transform.
Let's use the substandard finite difference approach to solve this equation [7]. Discretization of the time derivative gives
,                  (15)
Consider the wave equation on the gray scale for the fiber space method (k-space), which includes the non-linear characteristic of ultrasound, which can be represented as follows [10]:

where  - is the nonlinearity source,  - is the nonlinearity coefficient. The harmonic oscillations equation can be represented as follows
,           (16)
where  – additional source;  - is a spatial Fourier transform.

After the spatial Fourier transformation the equation can be expressed as follows

The introduction of the nonlinearity term in fiber space method makes it easier to calculate the actual relief temperature in heterogeneous large scale models.
Conclusions. To build a model of the ultrasonic field in a randomly inhomogeneous medium, the fiber spaces method (k-space), which increased the accuracy of parameter estimation field is used.

  • Kochura Ye.V. Avtomatizatsiya kontrolya kachestva usredneniya syr'ya pered yego obogashcheniyem, Obogashcheniye rud,no. 6 (1995): 41-43.
  • Morkun, V. Ultrasonic control of random heterogeneous mediums parameters.-Non Destructive Testing, ELSEVIER 92, (1991):989-993.
  • Morkun, V., Goncharov, S., Pikilnyak, A., Krivenko, A. Iron ore benefication processes optimization. ТЕKA. Commision of Motorization and Energetics in Agriculture 12, no.4 (2012): 162-166.
  • Morkun V. S., Morkun N. V,. Pikilnyak A.V. Ultrasonic facilities for the ground materials characteristics control, Metallurgical and Mining Industry, 2014, No2,
  • Morkun V. S., Morkun N. V,. Pikilnyak A.V. Iron ore flotation process control and optimization using high-energy ultrasound, Metallurgical and Mining Industry, 2014, No2.
  •  Fano U., Spencer L. Berger, M. Transfer of gamma radiation. Gosatomizdat, 1963
  • Mast T. D. et al. A k-space method for large-scale models of wave propagation in tissue, IEEE Trans.Ultrason., Ferroelect., Freq.Contr. 48, no.2, (2001): 341-354.
  • Tabei M., Mast T. D., Waag R. C.. A new k-space method for coupled first-order acoustic propagation equations, Journal of the Acoustical Society of America (submitted), (2001).
  • Morkun, V. Ultrasonic Control of Random Heterogeneous Mediums Parameters, AMMTIAC, no. 42(1992):
  • Chung A. H., Hynynen K. H., Cline H. E., Jolesz F. A. Quantification of thermal exposure using proton resonance frequency shift, Proc.SMR 4th Meeting, ISSN 1065-9889 3, 1751. (1996).




"Metaljournal" publishing
[email protected]

Project manager - Volodymyr Grynov

 Tel.:+38 (056) 794 36 74
 Fax: +38 (056) 794 36 75

 Mob: +38 (050) 320 69 72

 ISSN 20760507