Texts on physics, maths and programming (Несколько текстов для зачёта), страница 31
Описание файла
Файл "Texts on physics, maths and programming" внутри архива находится в папке "3". Документ из архива "Несколько текстов для зачёта", который расположен в категории "". Всё это находится в предмете "английский язык" из 5 семестр, которые можно найти в файловом архиве МГТУ им. Н.Э.Баумана. Не смотря на прямую связь этого архива с МГТУ им. Н.Э.Баумана, его также можно найти и в других разделах. Архив можно найти в разделе "остальное", в предмете "английский язык" в общих файлах.
Онлайн просмотр документа "Texts on physics, maths and programming"
Текст 31 страницы из документа "Texts on physics, maths and programming"
In conclusion, the properties of the solutions of Eq. (7) ( (11)) discussed above actually reveal the source of the problems encountered using numerical methods. Suppose the trajectory starts at a point (x0, v0) in R near (0, 0). Then this trajectory converges to the separatrix Γ very rapidly, in fact, the distance between θt(x0, v0) and Γ decreases with time as e−ct. Once θt is very close to Γ (which is ), it moves towards R near (π/β, 0) as eβkt/c. Hence, the time required to get very close to (π/β, 0) is roughly eβkt/c=π/β t=c/βk log (π/β), thus it follows from the orders of magnitude of the constants (c, k, β) that the time required is roughly O(102). So with a time step of about h = 0.1, the number of steps required to come close to seeing the asymptotic nature of the solution as t → ∞ is O(103). Using this many steps is not only computationally expensive; it also makes the contamination of the approximate numerical data by round off errors almost inevitable. Moreover, one can see from the intrinsic nature of the solutions described above that the basic problems with numerical solutions cannot be ameliorated through the use of standard rescaling techniques, nor can they completely circumvented using the usual stiff equation solvers. One infers from all of this that we shall need to employ other methods to determine useful approximate trajectories of Eq. (7) that stretch all the way from an initial point (x0, v0) to a very small neighborhood of the sink (π/β, 0) and that these methods must be of an asymptotic nature for t → ∞.
4. Equation for particle trajectories
The equation of motion of fine particles in a fluid with an acoustic sound field is given from Eq. (7). In the literature, the inertia term x″ is neglected because it is small and the singular perturbation solution shown in Eq. (8) is used.
Using Eq. (8) forces the initial velocity to be k/c sin βx0 (it is usually assumed = 0) and leads to errors in the solution that for very small time t, can be as large as k/c. The solution obtained by using the approximation equation (8) with initial value x(0) = x0 can easily be obtained:
| (29) |
By separation of variables:
| (29) |
| (30) |
| (31) |
To obtain a more accurate approximation than the singular perturbation solution, let the assumed solution of Eq. (7) be of the form:
x=u+xs. | (32) |
Substituting Eq. (32) in Eq. (7), one obtains:
| (33) |
Assuming that u is small, it follows from Eq. (30) that the terms in brackets above are small, and it is reasonable to approximate u by solving the following initial value problem:
| (34) |
Here it is assumed that x′(0) = v(0) = 0 in Eq. (7). The solution of Eq. (34) can be readily found by elementary means to be:
| (35) |
| (36) |
where,
| (37) |
| (38) |
and r1 >>>r2. It is therefore plausible that the following approximate solutions (t) are quite accurate for all t:
| (39) |
| (40) |
In fact, it can be proved rigorously that
| (41) |
and
| (42) |
Note also then that:
| (43) |
Therefore, particle trajectories in the acoustic medium can be calculated using Eq. (39) with the given values of constants (c, k, β). In graph Fig. 6, the vertical axis represents the time in seconds for the sound field treatment. The values on the horizontal axis represent the position of particles at corresponding times. The distances are in micrometers and for convenience taken at λ/16 of wavelength intervals from the first pressure anti-node at the face of the transducer. All the particles move toward the pressure node to which the particles are supposed to move. Particles with 2 μm diameter take about 20 s to reach the node. Fig. 5 shows where the sinks (pressure nodes) repeat at every 2π/β distance. Thus the behavior of the particle movement shown in Fig. 6 using the trajectory equation represents the actual behavior of the system.
Fig. 6. Particle trajectories.
5. Concentration equation
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 vs = k/c sinβxs is used, one can simply substitute v(x, t) = vs(x, t) = vs(x) = k/c sinβxs in Eq. (44), with xs a function of t and x0, 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) = x0, 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 x0 in terms of x and t be denoted by x0 = ψ (x, t). This cannot be done in closed form, but the following can readily be shown, by solving Eq. (30) for x0, 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 x0 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 r1 r2, 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 vs(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):
n(x,0)=n0(const). | (53) |
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(101). 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 (n0) 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.