Concentration equation

Date: 2015-10-07; view: 430.

The approximate solution of Eq. (39), obtained earlier will be used to obtain an approximate value of the particle concentrations. Since the particles are neither created nor disappear during the irradiation of sound, the equation of conservation of particles must hold:

(44)

In order to use this equation, one must find an expression for the velocity field v = v(x, t) as a function of just x and t. Now if the singular perturbation approximation v v_s = k/c sinβx_s is used, one can simply substitute v(x, t) = v_s(x, t) = v_s(x) = k/c sinβx_s in Eq. (44), with x_s a function of t and x₀, and get the result in the literature. A more accurate solution is obtained by employing the approximate formula of Eq. (39):

(45)

for the solution of Eq. (7) subject to the initial conditions x(0) = x₀, x′(0) = 0. In order to do this one must find the velocity of the particle at distance x from the pressure anti-node at time t > 0, so that one can find the approximation for v as a function of x and t.

To use the approximation equation (45) for the determination of the particle concentration, one must extract from it an approximation for v(x, t) to be substituted into Eq. (44), hence one must first solve:

(46)

Let the solution of the above equation for x₀ in terms of x and t be denoted by x₀ = ψ (x, t). This cannot be done in closed form, but the following can readily be shown, by solving Eq. (30) for x₀, to yield a good approximate solution for all 0 < x < π/β, t > 0:

(47)

From the above equation one can obtain a good approximation for the velocity field to be used in Eq. (44) by replacing x₀ with ψ in Eq. (40):

(48)

where n is the velocity at a given x and t. Substituting the value of ψ from Eq. (47) into Eq. (48), and using the fact that r₁ r₂, yields the approximation:

(49)

As e^−βkt/c 1 for all except very large t, Eq. (49) yields:

(50)

an even simpler approximation for v:

Hence,

(51)

Observe that converges very rapidly to v_s(x) = k/csinβx with increasing t. Substituting for v in Eq. (44), one obtains:

(52)

which must be solved subject to the initial condition (Cauchy problem):

Therefore, the following approximation formula for the concentration is applicable for any time t:

(54)

Observe that this is very close to the formula:

(55)

obtained using the singular perturbation approximation, except for t very close to zero. For 0 < t 1, the difference between Eq. (54) and the above (55) formula can be O(10¹). It is worth noting in closing that it can be shown by a long and tedious-but straightforward-computation that the error in using Eq. (54) is 2βk/c, for all 0 x π/β, t 0.

The ratio of the initial concentration (n₀) to final concentration (n) at that time (t), is taken as 1 when time = 0. Using Eq. (54) the particle concentrations with increasing t values are given on graph in Fig. 7 represent the concentration at different times; note how the concentrations are highest at the sink π/β (Fig. 5). Thus one sees that the approximate concentration closely follows the behavior of the actual system.

Fig. 7. Particle concentration.

The approximate solutions for particle trajectories and concentration will be used to compare the experimental observations with expected theoretical values. There is no closed form solution for Eq. (7), and direct application of a Runge–Kutta solver led to rather surprising results that were completely unsatisfactory. Stiff ODE solvers also appeared somewhat inadequate, so a viable approach was to use an approximate solution. If the experimental results can be compared to Fig. 6 and Fig. 7, then it can be stated that the mathematical model derived has solutions that are in good agreement with the actual behavior of the system and can quite possibly be generalized for any type of medium or particle.