<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AJCM</journal-id><journal-title-group><journal-title>American Journal of Computational Mathematics</journal-title></journal-title-group><issn pub-type="epub">2161-1203</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ajcm.2012.23024</article-id><article-id pub-id-type="publisher-id">AJCM-23170</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Numerical Modeling of the Measure of Global Environmental Needs with Applications Laser-LIDAR
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>erhat</surname><given-names>Mohammedi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Belgacem</surname><given-names>Zergui</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>H.</surname><given-names>Soubari</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Said</surname><given-names>Bensaada</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Physics, Université Mohamed Khider de Biskra, Biskra, Algeria</addr-line></aff><aff id="aff2"><addr-line>Photonics Systems Laboratory, Université de Strasbourg, Strasbourg, France</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>farwane@yahoo.fr(EM)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>28</day><month>09</month><year>2012</year></pub-date><volume>02</volume><issue>03</issue><fpage>194</fpage><lpage>198</lpage><history><date date-type="received"><day>March</day>	<month>28,</month>	<year>2012</year></date><date date-type="rev-recd"><day>June</day>	<month>12,</month>	<year>2012</year>	</date><date date-type="accepted"><day>June</day>	<month>20,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The functions of Bessel are used extensively in the various problems of the science and the technology. A laser offers practical remote sensing technologies for measuring environmental changes on both global and local scales. We describe a computer model that was developed to simulate the performance of three-dimensional (3D) laser radars (lidars). The principle of the problem consists in interpreting information on the absorption of the laser impulse in a spectral line assigned to the chemical body that one studied. Our purpose is to estimate the vertical variation of extinction and atmospheric transmission due to aerosol particles near the air-geographical surface interface. The feasibility and effectiveness of the proposed method is demonstrated by computer simulation.
 
</p></abstract><kwd-group><kwd>Transmission; Lasers; Optical Propagation; Atmospheric Effects</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The measure air pollution over large European cities with lidar, mobile differential absorption lidar (DIAL). Regardless of the measures aiming to reduce the broadcasts of atmospheric pollutants, it appears indispensable all at once to improve information and to pursue the efforts of research. Systems have been developed to produce 2D and 3D maps of concentrations of nitrous oxide, nitrogen dioxide, sulfuric dioxide, and ozone. The city, concentration gradients are high and rapidly changing periods of active traffic flow (<xref ref-type="fig" rid="fig1">Figure 1</xref>). In absorption of light in the spectral lines of the specter often provides precious information on the properties physics of matter [1-3]. This is how the displacements of the lines (effect Doppler) inform on the speed of the movement directed of matter, and the width of the line, on the temperature and the density of this one. One uses the laser radiance extensively to determine the content of the atmosphere in different chemical bodies and sprays, and more especially to detect concentrations petty of the sparkling impurities distributed in the atmosphere. In this study, a three-dimensional simulation MATLAB program for multi-pollutants dispersion from an industrial stack has been presented.</p></sec><sec id="s2"><title>2. Formalism of Method</title><p>The most efficient method is probably the one of the absorption compared, that implies the use of the laser radars (lidar). It consists has send in the atmosphere of the impulses lasers to two neighboring frequencies <img src="3-1100095\3d493e0e-083e-42a6-b785-f2f71ab6b127.jpg" /> and<img src="3-1100095\b6a3bebf-415f-4e2c-8c86-a2efa122d9e3.jpg" />, nearly confounds itself with the center of the line of absorption of the studied body, i.e. <img src="3-1100095\90776680-7d97-41c8-917f-197616aa983c.jpg" />and <img src="3-1100095\0c00e0ba-9d5f-4a4c-8c6d-f66df6c46a25.jpg" /> the other is located out of this line. All processes of interaction of the radiance with matter for the neighboring frequencies <img src="3-1100095\225a3f0a-d4fa-48bd-8b2a-5ea7dec00ba9.jpg" /> and <img src="3-1100095\1bbd9bca-f1c5-4651-a6d2-712541fc7ecb.jpg" /> are appreciably equal [4-6]. Are indeed <img src="3-1100095\2714f09a-ef53-428f-9fa6-2e42ebf1c728.jpg" /> the contour of the line absorption, i.e. the spectral intensity of the laser impulse; <img src="3-1100095\5a8c3174-8e5a-4abf-8c7b-99f302cbc65c.jpg" />the contour of the line of absorption of the studied body; <img src="3-1100095\c09fadd1-3a47-4b82-9304-5db69ebdd435.jpg" />the spectral coefficient of absorption for the other processes of interaction of the radiation with matter; The power of the laser radiation captured in the hypothesis of homogeneity of the atmosphere [7-9], will express itself then by:</p><disp-formula id="scirp.23170-formula73060"><label>(1)</label><graphic position="anchor" xlink:href="3-1100095\d77defda-8008-4917-b81a-e0a3d537c325.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="3-1100095\2a935644-7bf4-49f3-a505-8b2b72ba12fc.jpg" /> is the concentration (%) mass (of the component considered in the atmosphere; m: the absorbing matter mass crossed by the laser impulse m = LSρ<sub>0</sub>, L: the distance the impulse between emitter and the receiving S: area of the surface of the receiving antenna and ρ<sub>0</sub> the density of the atmosphere (<xref ref-type="fig" rid="fig2">Figure 2</xref>)</p><disp-formula id="scirp.23170-formula73061"><label>(2)</label><graphic position="anchor" xlink:href="3-1100095\ba234891-7b16-4596-b3d8-b6f3b47c1a4e.jpg"  xlink:type="simple"/></disp-formula><p>we deal with a simplified model of a coupled transport and Bessel (Figures 3-5) equations.</p><disp-formula id="scirp.23170-formula73062"><label>(3)</label><graphic position="anchor" xlink:href="3-1100095\3e828798-e454-4edd-a987-70f71883ef73.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23170-formula73063"><label>(4)</label><graphic position="anchor" xlink:href="3-1100095\c541bbdd-1c72-40f9-b7c3-b5dbe757db7b.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23170-formula73064"><label>(5)</label><graphic position="anchor" xlink:href="3-1100095\38105f7a-371f-40e1-8235-3fab2ca31be7.jpg"  xlink:type="simple"/></disp-formula><p>We can interpret the function transmittance as a wave that has to fulfill the homogeneous wave equation.</p><disp-formula id="scirp.23170-formula73065"><label>(6)</label><graphic position="anchor" xlink:href="3-1100095\d3bfd5c8-0374-4155-84c2-2db4bf6697e8.jpg"  xlink:type="simple"/></disp-formula><p><img src="3-1100095\713149d3-e56d-4b48-8634-efb8e8753be2.jpg" /></p><disp-formula id="scirp.23170-formula73066"><label>(7)</label><graphic position="anchor" xlink:href="3-1100095\bbbc66c4-c271-4290-81ac-10d863f67e04.jpg"  xlink:type="simple"/></disp-formula><p>When one suppose the following conditions</p><p><img src="3-1100095\4be1c58b-690b-4a64-81e6-472a014dab45.jpg" />, <img src="3-1100095\8d47c099-3e1b-4178-a314-1308c09094af.jpg" />, <img src="3-1100095\1368de34-63f3-480d-b143-5ea13ba61a86.jpg" /></p><p><img src="3-1100095\27c20b9f-8287-4982-a988-f0cc32e529de.jpg" />, and <img src="3-1100095\709b820f-98d8-4481-8219-f091a08a8357.jpg" /></p><p><img src="3-1100095\4d8d8da0-30de-4e6a-ba48-c24238ad6ecd.jpg" /></p><p>One will be able to write the transmittance:</p><disp-formula id="scirp.23170-formula73067"><label>, (8)</label><graphic position="anchor" xlink:href="3-1100095\506b267e-c53c-4d8e-92c8-c13077ce15d3.jpg"  xlink:type="simple"/></disp-formula><p><img src="3-1100095\0c296920-4cbb-4abb-aa71-ce7378c6920f.jpg" /></p><disp-formula id="scirp.23170-formula73068"><label>(9)</label><graphic position="anchor" xlink:href="3-1100095\7173bba7-a093-4141-af3e-3a6ffc81d7b3.jpg"  xlink:type="simple"/></disp-formula><p>For the left boundaries we have do discretize the following equation:</p><disp-formula id="scirp.23170-formula73069"><label>(10)</label><graphic position="anchor" xlink:href="3-1100095\c30adac7-5e66-4b71-9800-e4079d20efb3.jpg"  xlink:type="simple"/></disp-formula><p>Noticing that ρ &lt; 1 and that for value fixed of z and one has<img src="3-1100095\933d683e-0801-4507-80a5-ac8524f9cc16.jpg" />, <img src="3-1100095\ef025a9c-ce55-426f-a3a9-2cfe14d261c1.jpg" /></p><p>The gotten formulas permit to calculate the function transmission T for values arbitrary of<img src="3-1100095\5d8a5aa3-f6d8-4c98-ac58-49225cc8954e.jpg" />, <img src="3-1100095\4135a2f4-2067-4d1f-8a17-f19fd8dcf92d.jpg" />, <img src="3-1100095\ddee88f2-7b23-4bce-9cf7-c57fbe862e4a.jpg" />, <img src="3-1100095\d93fed57-0b12-4782-a731-5a16efdf68df.jpg" />, <img src="3-1100095\8b0452cc-b016-4b58-8ecd-2d3d89eb57a3.jpg" />, <img src="3-1100095\3cfa03b5-6039-43c8-8dd5-b90170335874.jpg" />, as well as to examine some cases different limits. Let’s suppose for example that<img src="3-1100095\e8ba19f7-afaf-4880-80cb-2ef9d0beeaf3.jpg" />, i.e. that the frequency of the signal confounds itself with the center of the absorption line [10-12]. One then:</p><p><img src="3-1100095\11dd2161-0b48-49ac-8cd6-f3a330b9eec7.jpg" />If<img src="3-1100095\03d2264e-9d17-405d-b6bb-822663178583.jpg" />, or <img src="3-1100095\1c1125bc-b2b8-4ce6-a075-0646b004acfc.jpg" /></p><p>One has then according to the Equation (10)</p><p><img src="3-1100095\c1fcbc5f-8ab7-4f3b-851b-cc9d9d224f17.jpg" /></p><p>In particular, if<img src="3-1100095\a0afd512-b9fe-4022-a8d8-742849d20728.jpg" />, one has <img src="3-1100095\0d3022f8-5c04-4e57-96f0-e3a69e631156.jpg" /></p><p>The contribution [Piazzola &amp; all] to coefficient K in inverse kilometers by aerosol particulates can written as</p><disp-formula id="scirp.23170-formula73070"><label>(11)</label><graphic position="anchor" xlink:href="3-1100095\51739425-08ed-4648-989e-5203d0327667.jpg"  xlink:type="simple"/></disp-formula><p>We calculate the atmospheric transmission using the expression</p><disp-formula id="scirp.23170-formula73071"><label>(12)</label><graphic position="anchor" xlink:href="3-1100095\03f268be-d1c7-4fc0-bd05-7e6fa5d1765a.jpg"  xlink:type="simple"/></disp-formula><p>The relative variation of the atmospheric transmission is given by</p><disp-formula id="scirp.23170-formula73072"><label>(13)</label><graphic position="anchor" xlink:href="3-1100095\83b84b88-bcd1-4392-b50c-fe395920a7bf.jpg"  xlink:type="simple"/></disp-formula></sec><sec id="s3"><title>3. Simulation Distribution-Diffusion</title><p>Before you begin to format your paper, for the 3D dimensional diffusion equation we apply a second order finite difference scheme in space and a higher order discretization scheme in time. We deal with higher order time-discretization methods. Therefore [13-16] we</p><p>propose the Runge-Kutta as adapted time-discretization methods to reach higher order results. For the time-discretization we use the following higher order discretization methods (Figures 6-8).</p><p>We deal with the following semi-discretized partial differential equations; such equations are used in each iterative splitting step:</p><disp-formula id="scirp.23170-formula73073"><label>(14)</label><graphic position="anchor" xlink:href="3-1100095\af5b1c6c-0a6a-407c-9f78-25f6437b15a6.jpg"  xlink:type="simple"/></disp-formula><p>where ν is the operator that we implicit solve in the equation and is the explicit operator, with a previous solution, e.g. last iterative solution. One supposes that the total flux <img src="3-1100095\0ce96d52-4650-42e4-b99e-97360e963e9d.jpg" /> one generalizes the problem (11) under the form: &#160;</p><p><img src="3-1100095\58016f94-2c33-461d-b1b0-ea28697f47e9.jpg" /><img src="3-1100095\31bb92e5-796b-4aed-9ac0-4e611c8fea7a.jpg" /></p><p>one approached the calculation with the method of Galerkin non consolidated.</p><p>For flexible specification of these model properties, a great number of variables and functions are available. with <img src="3-1100095\7f2f488d-6609-43ca-981c-9ddd1fcf7c07.jpg" /> = 10 m<sup>2</sup>/s, β = −0.3 m/s. In the simulations id given by Quartaroni &amp; All [17,18], one took a concentration u0 of 1 particle by m<sup>3</sup> a height of 10 km the number of global.</p><p>Peclet is therefore<img src="3-1100095\ca56b0de-1344-424b-a5a7-8b649b042abe.jpg" />, with <img src="3-1100095\71ee9c7d-e8f0-49b0-9c96-b2469a194517.jpg" /> = 10 m<sup>2</sup>/s, β = −0.3 m/s. Several examples on the effects of meteorological parameters (i.e., wind velocity, ambient air temperature, atmospheric stability and surface roughness) on pollutant dispersion were illustrated using the program.</p></sec><sec id="s4"><title>4. Conclusions</title><p>The aim of this paper was to assess the vertical variations in atmospheric transmission calculated from vertical aerosol concentration profiles recorded in stratospheric</p><p>area. The analysis of the data revealed a negative transmission gradient between 0.3 and 0.4 m height during winds of marine origin lower than 8 m/s, whereas a positive gradient occurs between 0.5 to 0.6 to 0.9 m height during high-wind-speed periods (<img src="3-1100095\8a863f04-84cc-4444-800e-4681c02e04ec.jpg" />&gt; 8 m/s), these distribution particles gradients induce a relative difference in atmospheric transmission between the two sample heights, which is maximal in the 3 - 5 μm band for winds lower 8 m/s, compared to the 8 - 12 μm band.</p><p>However, comparisons with in situ data show the need for realistic source function and to model the specific situation of the geographical location. The present paper provides a simple approach to the complex problem of the aerosol dynamics during severally seasons. A realistic simulation would take into account a set of the aqueous or gaseous chemistry equations which we have ignored here. This tool is under development in our laboratory, and includes specialized modules: dust/radiation interaction, wet scavenging, an improved dust scheme, and a more detailed emission map for the carbonaceous particles. Finally, the pollutants, emitted once in the atmosphere, can be transported on long distances and to cause some damages in regions relatively faraway of the places of emission.</p></sec><sec id="s5"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.23170-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">H. Nakatsuka, S. Asaka, H. Itoh, K. Ikeda, and M. Matsuoka, Observation of bifurcation to chaos in an all-optical bistable system, Phys. Rev. Lett., 50(2), 109-112, 1983. doi:10.1103/PhysRevLett.50.109</mixed-citation></ref><ref id="scirp.23170-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple"> K. Tamura, E.P. Ippen, H.A. Haus, and L.E. Nelson, 77-fs pulse generation from a stretched-pulse mode locked all fiber ring laser, Opt. Lett., 18(13), 1080-1082,1993. doi:10.1364/OL.18.001080</mixed-citation></ref><ref id="scirp.23170-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">G. Yandong, S. Ping, and T. Dingyuan, 298a??fs passively mode locked ring fiber soliton laser, Microw. Opt. Tech. Lett., 32(5), 320-333, 2002.  
doi:10.1002/mop.10170</mixed-citation></ref><ref id="scirp.23170-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">K.K. Gupta, N. Onodera, and M. Hyodo, Technique to generate equal amplitude, higher-order optical pulses in rational harmonically modelocked fiber ring lasers, Electron. Lett., 37(15), 948-950, 2001.  
doi:10.1049/el:20010584</mixed-citation></ref><ref id="scirp.23170-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Y. Shi, M. Sejka, and O. Poulsen, A unidirectional Er3+-doped fiber ring laser without isolator, IEEE Photon. Technol. Lett., 7(3), 290-292, 1995.  
doi:10.1109/68.372749</mixed-citation></ref><ref id="scirp.23170-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">J.M. Oh and D. Lee, Strong optical bistability in a simple L-band tunable erbium doped fiber ring laser, J. Quantum Electron., 40(4), 374-377, 2004.  
doi:10.1109/JQE.2004.824699</mixed-citation></ref><ref id="scirp.23170-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Q. Mao and J.W.Y. Lit, L-band fiber laser with wide tuning range bases on dual wave length optical bistability in linear  overlapping grating cavities, IEEE J. Quantum Electron., 39(10), 1252-1259, 2003.  
doi:10.1109/JQE.2003.817239</mixed-citation></ref><ref id="scirp.23170-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">L. Luo, T.J. Tee, and P.L. Chu, Bistability of erbium doped fiber laser, Opt. Commun., 146, 151-157, 1998.  
doi:10.1016/S0030-4018(97)00502-6</mixed-citation></ref><ref id="scirp.23170-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">P.P. Banerjee, Nonlinear Optics—Theory, Numerical Modeling, and Applications, Marcel Dekker Inc., New York, 2004.</mixed-citation></ref><ref id="scirp.23170-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">P. Meystre, On the use of the mean-field theory in optical bistability, Opt. Commun., 26(2), 277-280, 1978.  
doi:10.1016/0030-4018(78)90071-8</mixed-citation></ref><ref id="scirp.23170-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">L.A. Orozco, H.J. Kimble, A.T. Rosenberger, L.A. Lugiato, M.L. Asquini, M. Brambilla, and L.M. Narducci, Single-mode instability in optical bistability, Phys. Rev.A, 39(3), 1235-1252, 1989. doi:10.1103/PhysRevA.39.1235</mixed-citation></ref><ref id="scirp.23170-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">K. Ikeda, Multiple-valued stationary state and its instability of the transmitted light by a ring cavity, Opt. Commun., 30(2), 257-261, 1979.  
doi:10.1016/0030-4018(79)90090-7</mixed-citation></ref><ref id="scirp.23170-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">M.J. Ogorzalek, Chaos and complexity in nonlinear electronic circuits, World Sci. Ser. Nonlinear Sci., Ser. A, (22), 1997.</mixed-citation></ref><ref id="scirp.23170-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">M.P. Kennedy, Three steps to chaos—Part I: Evolution, IEEE Trans. Circuits Syst. I, Fund. Theory Appl., 40(10), 640-656, 1993.</mixed-citation></ref><ref id="scirp.23170-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">N.J. Doran and D. Wood, Nonlinear-optical loop mirror, Opt. Lett., 13(1), 56-58, 1988.  
doi:10.1364/OL.13.000056</mixed-citation></ref><ref id="scirp.23170-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">S. Zhang, Y. Liu, D. Lenstr, M.T. Hill, H. Ju, G.D. Khoe, and H.J.S. Dorren, Ringlaser optical flip-flop memory with single active element, IEEE J. Sel. Top. Quant. Electron., 10(5), 1093-1100, 2004.</mixed-citation></ref><ref id="scirp.23170-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">M.E. Fermann, F. Haberl, M. Hofer, and H. Hochreiter, Nonlinear amplifying loop mirror, Opt. Lett., 15(13), 752-754, 1990. doi:10.1364/OL.15.000752</mixed-citation></ref><ref id="scirp.23170-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Alfio Quarteroni, Ricardo Sacco , Fausto Saleri, Méthodes numériques,Algorithmes, analyse et applications, Edition Springer 2008.</mixed-citation></ref></ref-list></back></article>