Different solutions of the diffusion equation and its applications

In this report, we solved the advection–diffusion equation under pollutants deposition on the ground surface, taking wind speed and vertical diffusion depend on the vertical height. Also, we estimated a simple diffusion model from point source in an urban atmosphere and the conservative material with downwind was evaluated. Then, we calculated the extreme ground-level concentration as a function of stack height and plume rise in two cases. Comparison between the proposed models and the emission from the Egyptian Atomic Research Reactor at Inshas had been done. Lastly, we discussed the results in this report.


Background
The Gaussian plume equation was semi-analytical solution, assuming that wind velocity and eddy diffusivities were constant. The non-Gaussian models agree well with the observed data by Hinrichsen [17]. Wortmann et al. [35] and Essa et al. [10] calculated a new technique to get the concentration of contaminants in the planetary boundary layer.
Early solutions used only two dimensions [26,28,34]. The Green's function method was studied by Stakgold [30]. Yeh [37] calculated advection-diffusion equation in three dimensions under boundary conditions. Also, these equations have been studied in turbulent models for unbounded domain by [18,38]. Also, [33,32] utilized dispersion modeling. All techniques have been restricted to a single isolated point source located at the origin. The details of several solutions were investigated by Carslaw and Jaeger [7], Sutton [31], Yaglom [36], Pasquill [24], Berlyand [4] and Lin and Hildemann [21]. Also, the diffusion from a point source in an urban atmosphere was estimated by Essa and El-Otaify [9]. Heines and Peters [16], and Bennett [3] investigated the physical effect of model for the dry deposition of pollutants to a rough surface. Recently, Essa et al. [11,12] estimated the solution of the advection-diffusion equation in two dimensions with variable vertical eddy diffusivity and wind speed using Hankel transform. Also, Essa et al. [11,12] found the solution of advection-diffusion equation in three dimensions using Hankel transform.
In the first part: Pasquill and Smith [22,23] introduced the advection-diffusion equation in steady state for a continuous point source as follows: where u(z) and K(z) are the mean horizontal wind velocity and vertical eddy diffusivity as a function of vertical height. Assuming that u(z) = 0 at z = 0, as defined by [29] as follows:

Open Access
Beni-Suef University Journal of Basic and Applied Sciences where φ m (z/L) is the non-dimension wind gradient, L is the Monin-Obukhov length, u * is the friction velocity, and k is the von-Karman's constant equals 0.4. One can get by integrating Eq. (2) w.r.t (z) the following velocity profiles [9]: The vertical diffusion coefficient was taken from Hanna et al. [13,14]as follows: where K o is a constant value. Equation (1) is solved using the boundary conditions as follows: (1)-The mixing layer was assumed not allowing fluid to pass through to the pollutants And where h is the height of the planetary boundary layer (PBL).
(2)-Assuming that there is a deposition at the ground surface as follows: where v d is the deposition velocity.
The solution of Eq. (1) is solved theoretical under different stabilities [10] as follows: Integrating Eq. (1) from zero to h with respect to z, and applying Eqs. (9b) and (10), one gets:

For neutral conditions
Substituting C(x, z) and u (z) by Eqs. (11) and (3)  where F 0 is a constant. Then, the solution in the neutral case was obtained as follows: Taking where x dn is the decay distance of pollutant in neutral condition. Equation (17) becomes: Now, F 0 can be determined by using the formula: where Q is the emission rate which after substitution of Eqs. (3) and (19), one gets: Then, the concentration of pollutants in neutral condition becomes: Similarly, C(x, z) in stable condition takes form [10]: where X ds is given in the stable case as follows: where M = N -(5.2 h 2 /12L). Also, C(x,z) in unstable condition is given [10] as follows: where X du is given in unstable case as follows: where D is taken from [10].
In the second part: We assumed two-dimensional structure with homogeneity in the lateral coordinate, where ū(z) is described in Fig. 1. The ground surface is treated with no deposition of matter exist.
The conservation of mass can be taken from [9]: where u(z) is the average wind velocity, C(z) is the concentration of pollutants, and H is the effective stack height of the plume.

The Effective stack Height (H)
We estimated the height of the plume as follows: where w o is the exit velocity of the pollutants (m/s) and D 1 is the inside stack diameter (m). Then, H equals: The concentration profile will be assumed in form [9]: where C o is the concentration value at the edge of the plume. C(z) is the concentration at vertical height; α 1 and α 2 , etc., are constants.
The fitting in Fig. 2 is obtained by taking the first two terms as follows: If α 1 is a r percentage of the concentration as follows: if r = 0, then After substituting, Eq. (26) becomes using different stabilities as follows:

Neutral case
The integration over "z" yields:

Stable case
After integration, C o takes the form: 6 Unstable case where f(z) stands for: f (z) = 16(z+z 0 ) L . By substituting from the above equation into Eq. (35) and integrating w.r.t "z, " one gets: where was taking from [9]: The time dosage integral of the concentration "D" is defined as: The mass balance condition takes the form: Substituting from Eq. (36) in Eq. (37), one gets: In the third part: The maximum ground-level concentration with two cases was estimated [9]. The concentration at the ground reaches a maximum and thereafter decreases as given by Hans et al. [15] with the formula: σ i (i = y, z): are dispersion parameters of concentration in the lateral and vertical directions (m).
The extreme concentration at ground occurs along the plume centerline (y = 0) at ground (z = 0) [25]). σ y and σ z can be estimated as follows: where a, b, c, d values are taken from Curtiss [8].
The ground-level concentration along centerline was estimated from Eq. (40) taking σ y and σ z from Eq. (41) one gets: The maximum concentration value can be obtained after solving the following equation: The maximum concentration and effective height were obtained as follows:

First: the effective height as a function of x
A-Unstable or neutral conditions: The buoyancy flux parameter, F b, was written [5] as follows: For (h s < 305 m), g is the acceleration (ms −2 ), v s is the exit vertical speed (ms −1 ), r s is the exit radius (m), T s is the gas exit temperature (K), and T a is the ambient temperature (K) at h s [27].
The critical x* is given by: For x ≤ x*, we have: Taking the constant equals 1.6 [6], one gets: This equation is used if T s > T a (2/3 law) [27]. u h can be calculated at 10 m as follows: where the parameter p is given from [20]. Substituting from (47)  The maximum value for χ * has the form: In neutral condition, taking d = 0.95 [8] χ * becomes: In stable stability, there are two methods for maximum concentration at z = 0 as follows: (First)-In stable, where d = 0.67 [8] the maximum value for χ * is: (Second)-In stable (E and F), the stability parameter is written as follows: where �θ �Z = 0.02 K/m for E and �θ �Z = 0.035 K/m for F ( [27].
h is written as follows: H* becomes: (53) Also, χ* has the form:

Second case: H is a constant of downwind distance:
If H is constant, then, ∂H ∂x = 0 , assuming H c = h s + 3w o D/u, where H c is constant. H* becomes: The maximum value for χ * has the form:

First case study
The decay distance and the concentration of pollutant at Inshas, Nuclear Research Center, Egyptian Atomic Energy Authority, were calculating. The stability was found in stable condition by using these data. The value of L was taken from Gifford work [13]. One uses Eq. (4) to calculate the values of u * in terms of the wind velocity where f is the Coriolis parameter. v d equals 0.04 m/s for Iodine [19]. Hence, the decay distance of iodine has been estimated using Eq. (23). The friction velocity and the mixing height with downwind distance are illustrated in Fig. 3. One gets the change of the friction velocity and mixing height from 0.25 to 0.3 m/s and 180 to 190 m, respectively [10]. The values of mixing height proportional with the decay distance are shown in Fig. 4; the normalized concentration (C/Q) of iodine at heights z = 0, 10, 50 m as a function of downwind distance has been evaluated , for h/L > 0 by Eq. (22). One finds that the maximum normalized ground-level concentration equals 18*10 -6 at 100 m and then decreases to reach near minimum value at 5000 m as shown in Fig. 5.

Second case study
C o /Q is applied on the first research reactor at Inshas, Nuclear Research Center, Egyptian Atomic Energy Authority. h s equals 43 m, D equals 1 m, w o equals 4 m/s, and the total ventilation equals 39,965 m 3 /hr. Also, one takes α = -1. Figure 6 shows that a straight-line fit well to these data in three conditions between Co/Q and H.

Third case study
One gets the maximum concentration on the second research reactor at Inshas, Nuclear Research Center, Egyptian Atomic Energy Authority. h s equals 27 m, D  [1]. One finds that the maximum concentrations are 6.2*10 -4 , 4.8* 10 -4 and 2*10 -4 Bq/m 3 at 40 m in unstable, neutral and stable conditions, respectively; when plume rise depends on downwind distance as shown in Fig. 7, these values are less than the maximum system activity (MSA). One finds that the effective height is quadratic with downwind distance in unstable case and the correlation coefficient (R) equals 0.99, (x max ) equals 1000 m, and (H max ) reaches 978 m. In neutral case, the curve is still quadratic and smaller than that in the unstable case, R = 1. Also, in the stable case, the relation is linear, R = 0.984, x max reaches to maximum value, and H max will reach 59 m as shown in Fig. 8.
We find that the maximum concentrations are 2.5* 10 -4 , 2.1* 10 -4 and 2* 10 -4 Bq/m 3 in neutral, unstable and stable conditions, respectively, at 40 m when H is a constant as shown in Fig. 9. Also, the relation between "H" and "x" is linear in unstable condition, H max = 163 m when x max = 1000 m and R = 0.997. In neutral case, the relation is quadratic and R = 1 and H max = 83 m at x max = 1000 m. Also, in stable case the relation is still quadratic and R = 0.998, H max = 58.5 m at x max = 1000 m as shown in Fig. 10. We must take H that depends on x because one gets the values of "H" and the concentration near the true. This is the best result to reach it.

Conclusions
The concentration released from an elevated point source in the presence of both deposition and elevated mixed layer is estimated. The decay distance in the downwind distance is calculated. The maximum normalized concentration equals 18*10 -6 (s/m 3 ) at 100 m and then decreases to reach minimum value at 5000 m at the ground level. The logarithmic law of wind velocity in different conditions is used to get the plume rise, effective stack height and the normalized concentration at the axis of the plume from the reactor release through different stability classes. We find that the ground-level concentration of iodine (I 131 ) agrees with the observed concentrations.
The values of extreme concentration are 6.2*10 -4 , 4.8* 10 -4 and 2*10 -4 Bq/m 3 at 40 m in unstable, neutral and stable conditions, respectively; when the plume rise depends on downwind distance and then decreases to values less than MSA, after that the values are close to zero at 400 m from the stack. In stable case, the maximum values of the ground-level concentration of air