Conference on Computation Physics-2006  (I27)
This presentation is the property of its rightful owner.
Sponsored Links
1 / 51

Xiwei HU, Zhonghe JIANG, Shu ZHANG and Minghai LIU PowerPoint PPT Presentation


  • 57 Views
  • Uploaded on
  • Presentation posted in: General

Conference on Computation Physics-2006 (I27) The propagation of a microwave in an atmospheric pressure plasma layer : 1 and 2 dimensional numerical solutions. Xiwei HU, Zhonghe JIANG, Shu ZHANG and Minghai LIU H uazhong U niversity of S cience & T echnology Wuhan, P. R. China

Download Presentation

Xiwei HU, Zhonghe JIANG, Shu ZHANG and Minghai LIU

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


Xiwei hu zhonghe jiang shu zhang and minghai liu

Conference on Computation Physics-2006 (I27)The propagation of a microwave in an atmospheric pressure plasma layer:1 and 2 dimensional numerical solutions

Xiwei HU, Zhonghe JIANG,

Shu ZHANG and Minghai LIU

Huazhong University of Science &Technology

Wuhan, P. R. China

August 30, 2006


Xiwei hu zhonghe jiang shu zhang and minghai liu

IIntroductionand motivationIIOne dimensional solutionIIITwo dimensional solutionIVConclusions


I introduction and motivation

IIntroduction and motivation


The classical mechanism

The classical mechanism

  • firstly, the EM wave transfer its wave energy to the quiver kinetic energy of plasma electrons through electric field action of waves.

  • Then, the electrons transfer their kinetic energy to the thermal energy of electrons, ions or neutrals in the plasmas through COLLISIONS between electrons or between electrons and other particles.


The electron fluid motion equation

The electron fluid motion equation

  • f0 is the microwave frequency,

  • νee , νei and νe0 is the collision frequency of electron-electron, electron-ion and electron-neutral, respectively.


Xiwei hu zhonghe jiang shu zhang and minghai liu

Pure plasma (produced by strong laser):νe=νee+νei,

Pure magnetized plasma (in magnetic confinement devices, e.g. tokamak): νe=0,

The mixing of plasma and neutral (in ionosphere or in low pressure discharge): νe=νe0.

In all of above cases:νe / f0 << 1

Taking the WKB (or ekonal) approximation

The solution of electron fluid equation is


The appleton formula

The Appleton formula


Xiwei hu zhonghe jiang shu zhang and minghai liu

Whenp=50 – 760 Torr

νe0≈6-466G(109) Hz,

electron density of APP

ne ≈1010– 1012 cm-3,

correspondent cut off frequency

ωc≈2- 20 GHz,

so

νe0 ≥or >>ωc≈2πf0.

f0 : frequency of electromagnetic wave


The goal of our work

The goal of our work

  • Study the propagation behaviors of microwave by solving the coupled wave (Maxwell) equation and electron fluid motion equation directly in time and space domain instead of in frequency and wavevector domain.


Ii one dimensional case

II One dimensional case

II.1 The integral-differential equation

II.2 The numerical method, basic wave form and precision check

II.3 The comparisons with the Appleton formula

II.4 Outline of numerical results


Ii 1 the integral differential equation

II.1The integral-differential equation


The coupled set of equations

The coupled set of equations

  • Begin with the EM wave equation

  • Coupled with the electron fluid motion equation


Xiwei hu zhonghe jiang shu zhang and minghai liu

  • Combinewave and electron motion equations, we have got a integral-differential equation:

  • Obtain numerically the full solutions of EM wave field in space and time domain


Ii 2 the numerical method precision check and basic wave forms

II.2The numerical method, precision check and basic wave forms


Numerical method

Numerical Method

  • Compiler:

    Visual C++ 6.0

  • Algorithm:

    —average implicit difference method for differential part

    —composite Simpson integral method for integral part


Check the precision of the code

Check the precision of the code

  • Compare the numerical phase shift with the analytic result inνe0 =0.

  • The analytic formula for phase shift


Bell like electron density profile

Bell-like electron density profile


Phase shift when e0 0

Phase shift Δφ when νe0 =0


Waveform of e y x n e 0 5 n c d 2 0 e0 0 1 0

Waveform of Ey (x)ne = 0.5 nc,d = 2 λ0, νe0 = 0.1 ω0


Wave forms passed plasma passed vacuum interference phase shift n e 0 5 n c d 2 0 e0 1 0 0

Wave forms: passed plasma, passed vacuum, interference, phase shift.ne = 0.5 nc,d = 2 λ0, νe0 = 1.0 ω0


The reflected plane wave e 2

The reflected plane wave E2


Xiwei hu zhonghe jiang shu zhang and minghai liu

II.3

The comparison with the Appleton formula


Brief summary 1

Brief summary (1)

  • When n0 /nc <1, the reflected wave is weak, the Δφ andT obtained from analytic (Appleton) formula and numerical solutions are agree well.

  • When n0 /nc >1, the wave reflected strongly, the Appleton formula is no longer correct. We have to take the full solutions of time and space to describe the behaviors of a microwave passed through the APP.


Ii 4 outline of numerical results phase shift transmissivity t reflectivity r absorptivity a

II.4 Outline of numerical resultsPhase shift ΔφTransmissivity TReflectivity RAbsorptivity A


Determination

Determination

  • E0—incident electric field of EM wave,

    E1—transmitted electric field,

    E2—reflected electric field

  • Transmissivity:

    T=E1 /E0 , Tdb =-20 lg (T).

  • Reflectivity:

    R=E2 /E0 , Rdb =-20 lg (R).

  • Absorptivity: A=1 - T2 - R2


Three models of n e x n e m x dx n e constant m 1 2 3

The bell-like profile

2. The trapezium profile

3. The linear profile

Three models of ne(x)∫ne{m} (x) dx =Ne=constant, m=1,2,3.


Effects of profiles are not important

Effects of profiles are not important


The phase shift

The phase shift | Δφ |


Briefly summary 2

Briefly summary (2)

1.\Δφ\ increases withn0andd.

2. When νe0 → 0,\Δφ\ → the maximum value in pure (collisionless) plasmas.

3. Then, \Δφ\ decreases withνe0/ω0increasing.

4. When νe0/ω0 >>1, Δφ→0 –the pure neutral gas case.


The transmissivity t db and the absorptivity a reach their maximum at e0 0 1

The transmissivity Tdb and The absorptivity A reach their maximum atνe0/ω0 ≈1


Briefly summary 3

Briefly summary (3)

  • All four quantities Δφ, T, R, A depend on

    --the electron density ne(x),

    --the collision frequency νe0 ,

    --the plasma layer width d.


Is more important than and d

is more important than and d

  • According to the collision damping mechanism, the transferred wave energy is approximately proper to the total number of electrons, which is in the wave passed path.

  • represents the total number of electrons in a volume with unit cross-section and width d when the average linear density of electron is .


T db seems a simple function of the product of n and d

TdB seems a simple function of the product of n and d

  • Let

    TdB (nd)=F(ne , νe)

  • When νe> 1,

    F(ne , νe) = Const.

  • When νe < 1 ,

    F(ne , νe) increases slowly with ne


F n e e

F(ne ,νe )


Iii two dimensional case

III Two dimensional case

III.1 The geometric graph and arithmetic

III.2 Comparison between one and two dimensional results in normal incident case

III.3 Outline of numerical results


Iii 1 geometric graph for fdtd integral differential equations

III.1Geometric graph for FDTDIntegral-differential equations


When microwave obliquely incident into an app layer

When microwave obliquely incident into an APP layer

  • The propagation of wave becomes a problem at least in two dimension space.

  • Then, the incidence angleθand thepolarization (S or P mode) of incident wave will influence the attenuation and phase shift of wave.


The equations in two dimension case

The equations in two dimension case

  • Maxwell equation for the microwave.

  • Electron fluid motion equation for the electrons.


Xiwei hu zhonghe jiang shu zhang and minghai liu

s-polarized p-polarized


Xiwei hu zhonghe jiang shu zhang and minghai liu

Combine Maxwell’s and motion equationsintegral-differential equations

  • S-polarized integral-differential equations:

  • P-polarized integral-differential equations:


Iii 2 comparison between one and two dimensional results in normal incident case

III.2Comparison between one and two dimensional results in normal incident case


Iii 3 the numerical results about the effects of incidence angles and polarizations

III.3The numerical resultsabout the effects of incidence angles and polarizations


The influence of incidence angle

The influence of incidence angle


The effects of the density profile

The effects of the density profile


Iv conclusion

IVConclusion


Xiwei hu zhonghe jiang shu zhang and minghai liu

1. When nmax /nc >1, the Appleton formula should be replayed by the numerical solutions.

2. The larger the microwave incidence angle is, the bigger the absorptivity of microwave is.

3. The absorptivity of P (TE) mode is generally larger than the one of S (TM) mode incidence microwave.


Xiwei hu zhonghe jiang shu zhang and minghai liu

4. The bigger the factor is, the better the absorption of APP layer is.

5. The absorptivity reaches it maximum when .

6.The less the gradient of electron density is, the larger (smaller) the absorptivity (reflectivity) is.


Xiwei hu zhonghe jiang shu zhang and minghai liu

Thanks !


  • Login