Medical, Pharma, Engineering, Science, Technology and Business

^{1}Mathematics and Theoretical Physics, NRC, Atomic Energy Authority, Cairo, Egypt

^{2}Department of Mathematics, Faculty of Science, Zagazig University, Egypt

- *Corresponding Author:
- Essa KSM

Mathematics and Theoretical Physics

NRC, Atomic Energy Authority, Cairo, Egypt

**Tel:**0020244717553

**E-mail:**mohamedksm56@yahoo.com

**Received date:** November 16, 2015 **Accepted date:** December 04, 2015 **Published date:** December 10, 2015

**Citation:** Essa KSM, Marrouf AA, El-Otaify MS, Mohamed AS, Ismail G (2015) New Technique for Solving the Advection-diffusion Equation in Three Dimensions using Laplace and Fourier Transforms. J Appl Computat Math 4:272. doi:10.4172/2168-9679.1000272

**Copyright:** © 2015 Essa KSM, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Applied & Computational Mathematics

A steady-state three-dimensional mathematical model for the dispersion of pollutants from a continuously emitting ground point source in moderated winds is formulated by considering the eddy diffusivity as a power law profile of vertical height. The advection along the mean wind and the diffusion in crosswind and vertical directions was accounted. The closed form analytical solution of the proposed problem has obtained using the methods of Laplace and Fourier transforms. The analytical model is compared with data collected from nine experiments conducted at Inshas, Cairo (Egypt). The model shows a best agreement between observed and calculated concentration.

Advection-diffusion equation; Laplace transform; Fourier transform; Bessel function

Environmental problems caused by the huge development and the big progress in industrial, which cause's a lot of pollutions. The transport of these pollutants can be adequately described by the advection–diffusion equation. In the last few years, there has been increased research interest in searching for analytical solutions for the advection–diffusion equation (ADE). Therefore, it is possible to construct a theoretical model for the dispersion from a continuous point source from an Eulerian perspective, given adequate boundary and initial conditions and the knowledge of the mean wind velocity field and of the concentration **turbulent fluxes**. The exact solution of the linear advection–dispersion (or diffusion) transport equation for both transient and steady-state regimes has been obtained [1].

The two- and three-dimensional advection–diffusion equation with spatially variable velocity and diffusion coefficients has been provided analytically [2]. A mathematical treatment has been proposed for the ground level concentration of pollutant from the continuously emitted point source [3]. Essa and El-Otaify studied a mathematical model for hermitized atmospheric dispersion (self adjoined by itself) in low winds with eddy diffusivities as linear functions of the downwind distance [4]. More recently the generalized analytical model describing the crosswind-integrated concentrations is presented [5]. Also, an analytical scheme is described to solve the resulting two dimensional steady-state advection–diffusion equation for horizontal wind speed as a generalized function of vertical height above the ground and eddy diffusivity as a function of both downwind distance from the source and vertical height.

On the other hand, the literature presents several methods to analytically solve the partial differential equations governing transport phenomena [6-10]. For example, the method of separation-of-variables is one of the oldest and most widely used techniques. Similarly, the classical Green’s function method can be applied to problems with source terms and inhomogeneous boundary conditions on finite, Semiinfinite, and infinite regions [10,11]. Integral transform techniques, such as the **Laplace** and Fourier transform methods, employ a mathematical operator that produces a new function by integrating the product of an existing function and a kernel function between suitable limits.

In this study, we obtained a **mathematical model** for dispersion of air pollutants in moderated winds by taking into account the diffusion in vertical height direction and advection along the mean wind. The eddy diffusivity is assumed to be power law in the vertical length. We provided analytical solutions to the advection–diffusion equation for three-dimensional with the physically relevant boundary conditions. The moderate data collected during the convective conditions. From nine experiments conducted at Inshas site, Cairo-Egypt [4], which used to investigate the analytical solution.

The dispersion of pollutants in the atmosphere is governed by the basic atmospheric diffusion equation. Under the assumption of incompressible flow, atmospheric diffusion equation based on the Gradient transport theory can be written in the rectangular coordinate system as:

(1)

where C(x, y, z) is the mean concentration of a pollutant (Bq/m^{3}), (μg/m^{3}) and (ppm); in which t is the time, S and R are the source and removal terms, respectively; (u, v, w) and (k_{x}, k_{y}, k_{z}) are the components of wind and diffusivity vectors in x, y and z directions, respectively, in an Eulerian frame of reference.

The following assumptions are made in order to simplify equation (1):

1) Steady -state conditions are considered, i.e., ∂C/∂t = 0

2) We are going to study Eq. (1) in case, when the components of wind (v, w) tends to zero

3) Source and removal (physical/chemical) pollutants are ignored so that S=0 and R=0

4) Under the moderate to strong winds, the transport due advection dominates over that due to longitudinal diffusion:

With the above assumptions, equation (1) reduces to:

(2)

Under the following boundary conditions:

(3)

(4)

(5)

(6)

where δ (…) is Dirac's delta function, and h is the mixing height. The wind speed u and eddy diffusivity k_{y}, k_{z} is expressed as a functions of power law of z as:

(7)

Where *α ,β ,m,n* are turbulence parameters and depend on **atmospheric stability**.

Eq. (2) can solve analytically as follows:

Transform the variable x to s by applying the Laplace transform on Eq. (2) to become

(8)

Again transform the variable y to λ by applying the Fourier transform on Eq. (8) to become

(9)

Eq. (9) simplified to the form

(10)

Now we will solve the homogeneous equation of Eq. (10) which takes the form

(11)

Transform then Eq. (11) becomes

(12)

Again transform then Eq. (12) become

(13)

where

But Eq. (13) is modified Bessel equation which has solution [12].

(14)

(15)

where A and B are constant

Now the general solution of the non-homogeneous Eq. (10) takes the form

(16)

where A_{*} and B_{*} are constants

Apply the boundary condition Eq. (3) on Eq. (16) which become

(17)

Apply the boundary condition Eq. (5) on Eq. (16) which gives

(18)

Substitute B_{*} Eq. (16) in Eq. (17) which gives:

(19)

Apply inverse Fourier and inverse Laplace respectively on Eq. (19) we get:

(20)

The diffusion data for the estimating were gathered during ^{135}I **isotope tracer**nine experiments in moderate wind with unstable conditions at Inshas, Cairo. During each run, the tracer was released from source has height 43 m for twenty four hours working, where the air samples were collected during half hour at a height 0.7 m. We collected air samples from 92 m to 184 m around the source in AEA, Egypt. The study area is flat, dominated by sandy soil with poor vegetation cover. The air samples collected were analyzed in Radiation Protection Department, NRC, AEA, Cairo, Egypt using a high volume air sampler with 220 V/50 Hz bias [13]. Meteorological data have been provided by the measurements done at 10 m and 60 m.

For the concentration computations, we require the knowledge of wind speed, wind direction, source strength, the dispersion parameters, mixing height and the **vertical scale velocity**. Wind speeds are greater than 3 m/s most of the time even at 10 m level. Further the variation wind direction with time is also visible. Thus in the present study, we have adopted dispersion parameters for urban terrain which are based on power law functions. The analytical expressions depend upon downwind distance, vertical distance and atmospheric stability. The atmospheric stability has been calculated from Monin-Obukhov length scale (1/L) [14] based on friction velocity, temperature, and surface heat flux.

The concentration is computed using data collected at vertical distance of a 30 m multi-level micrometeorological tower. In all a test runs were conducted for the purpose of computation. The concentration at a receptor can be computed in the following way:

Applying formula Eq. (21) which contains eddy diffusivities as function with power law at y = 0.0 for half hourly averaging.

As an illustration, results computed from these approaches are shown in **Table 1**, for nine typical tests conducted at Ins has site, Cairo-Egypt [4]. This table shows that the observed and predicted concentrations for ^{135}I using Eq. (20) with power law of eddy diffusivities and the wind speed are very near to each other of ^{135}I.

**Figure 1** shows the variation of predicted and observed concentration of ^{135}I with the downwind distance. One gets very good agreement between observed and predicted concentration.

**Figure 2** shows that the predicted concentrations which are estimated from Eq. (20) are a factor of two with the observed concentration.

Now, the statistical method is presented and comparison among analytical, statically and observed results will be offered [13]. The following standard statistical performance measures that characterize the agreement between prediction (C_{p} = C_{pred}) and observations (C_{o}=C_{obs}):

1. Normalized mean square error (NMSE), It is an estimator of the overall deviations between predicted and observed concentrations. Smaller values of NMSE indicate a better model performance. It is defined as:

2. Fractional bias (FB): It provides information on the tendency of the model to overestimate or underestimate the observed concentrations. The values of FB lie between -2 and +2 andit has a value of zero for an ideal model. It is expressed as:

3. Correlation coefficient (R): It describes the degree of association between predicted and observed concentrations and is given by:

4. Fraction within a factor of two (FAC2) is defined as:

FAC2 = fraction of the data for which

0.5 ≤ (C_{p}/C_{o}) ≤ 2

Where σp and σo are the standard deviations of C_{p} and C_{o} respectively. Here the over bars indicate the average over all measurements (Nm). A perfect model would have the following idealized performance: NMSE = FB = 0 and COR = FAC2 = 1.0

From the statistical method of **Table 2**, we find that the predicted concentrations for 135I lie inside factor of 2 with observed data. Regarding to NMSE, FB and COR the predicted concentrations for ^{135}I are better with observed data.

In this paper, a steady-state three-dimensional mathematical model for the dispersion of pollutants from a continuously emitting ground point source in moderated winds is formulated. Besides advection along the mean wind, the model takes into account the diffusion in crosswind and vertical directions. The eddy diffusivity and the wind speed are assumed to be power law in the vertical height z.

The closed form analytical solution of the proposed problem has obtained using the methods of Laplace and **Fourier transforms**.

In general, the present model is compared with data collected from nine experiments conducted at Inshas, Cairo (Egypt). One gets the predicted concentrations are in a best agreement with the corresponding observation. Moreover, the Statistical results here are in agreement with the analytical results.

This work has been completed by the support of Egyptian Atomic Energy Authority, and the authors thank for this support. The First author’s thank is extended to all members of Mathematics and theoretical Physics for providing the experimental data of ^{135}I.

- Guerrero JSP, Pimentel LCG, Skaggs TH, van Genuchten MT (2009) Analytical solution of the advection-diffusion transport equation using a change-of-variable and integral transform technique. International Journal of Heat and Mass Transfer 52: 3297-3304.
- Zoppou C, Knight JH (1999) Analytical solution of a spatially variable coefficient advection diffusion equation in up to three dimensions. Applied Mathematical Modeling 23: 667-685.
- Embaby M,Mayhoub AB,Essa KSM,Etman S (2002) Maximum ground level concentration of air pollutant.Atmosfera 15: 185-191.
- Essa KSM, El-Otaify MS(2007) Mathematical model for hermitized atmosphericdispersion in low winds with eddy diffusivities linear functions downwind distance. Meteorology and Atmospheric Physics 96: 265-275.
- Kumar P,Sharan M (2014) An analytical model for dispersion of pollutants from a continuous source in the atmospheric boundary layer.Proc R Soc A pp:1-24.
- Courant D, Hilbert D (1953)Methods of Mathematical Physics.Wiley Interscience Publications, New York, USA.
- Morse PM, Feshbach H (1953) Methods of Theoretical Physics. McGraw-Hill, New York, USA.
- Carslaw HS, Jaeger JC (1959) Conduction of Heat in Solids(2ndedn.).Oxford University press, Oxford.
- Sneddon IN (1972) The Use of Integral Transforms.McGraw-Hill.
- Ozisik MN (1980) Heat Conduction.Wiley, New York.
- Leij FJ, v an Genuchten MT (2000) Analytical modeling of non-aqueous phase liquid dissolution with Green’s functions. Transport in Porous Media 38: 141-166.
- Irving J,Mullineux (1959) Mathematics in Physics and Engineering. Pure and applied physics.
- Essa KSM, Mubarak F,Khadra SA (2005) Comparison of Some Sigma Schemes for Estimation of Air Pollutant Dispersion in Moderate and Low Winds.Atmospheric Science Letter 6:90-96.
- Donald G (1972)Relation among Stability Parameters in the Surface Layer. Boundary LayerMeteorology 3: 47-58.

Select your language of interest to view the total content in your interested language

- Adomian Decomposition Method
- Algebraic Geometry
- Analytical Geometry
- Applied Mathematics
- Axioms
- Balance Law
- Behaviometrics
- Big Data Analytics
- Binary and Non-normal Continuous Data
- Binomial Regression
- Biometrics
- Biostatistics methods
- Clinical Trail
- Complex Analysis
- Computational Model
- Convection Diffusion Equations
- Cross-Covariance and Cross-Correlation
- Differential Equations
- Differential Transform Method
- Fourier Analysis
- Fuzzy Boundary Value
- Fuzzy Environments
- Fuzzy Quasi-Metric Space
- Genetic Linkage
- Hamilton Mechanics
- Hypothesis Testing
- Integrated Analysis
- Integration
- Large-scale Survey Data
- Matrix
- Microarray Studies
- Mixed Initial-boundary Value
- Molecular Modelling
- Multivariate-Normal Model
- Noether's theorem
- Non rigid Image Registration
- Nonlinear Differential Equations
- Number Theory
- Numerical Solutions
- Physical Mathematics
- Quantum Mechanics
- Quantum electrodynamics
- Quasilinear Hyperbolic Systems
- Regressions
- Riemannian Geometry
- Robust Method
- Semi Analytical-Solution
- Sensitivity Analysis
- Smooth Complexities
- Soft biometrics
- Spatial Gaussian Markov Random Fields
- Statistical Methods
- Theoretical Physics
- Theory of Mathematical Modeling
- Three Dimensional Steady State
- Topology
- mirror symmetry
- vector bundle

- Total views:
**7707** - [From(publication date):

December-2015 - Dec 03, 2016] - Breakdown by view type
- HTML page views :
**7655** - PDF downloads :
**52**

OMICS International Journals

- Clinical Research
- Radiology
- Toxicology
- Nephrology
- Haematology
- Ophthalmology
- Anesthesiology
- Gasteroenterology

- Nursing
- Health Care
- Alternative Medicine
- Medicine
- Oncology
- Neurology
- Diabetes & Endocrinology
- Psychiatry

OMICS International Conferences 2016-17