Outline Introduction neering Composite Materials Design with - - PDF document

outline
SMART_READER_LITE
LIVE PREVIEW

Outline Introduction neering Composite Materials Design with - - PDF document

12/4/2010 Cullen College of Engineering University of Houston Design/Modeling for Periodic neering trical & Computer Engin N Nano Structures for EMC/EMI St t f EMC/EMI Ji Chen Elect Department of Electrical and Computer Engineering


slide-1
SLIDE 1

12/4/2010 1

University of Houston Cullen College of Engineering neering

Design/Modeling for Periodic N St t f EMC/EMI

trical & Computer Engin

Ji Chen

Nano Structures for EMC/EMI

Elect

1

Department of Electrical and Computer Engineering University of Houston Houston, TX, 77204

University of Houston Cullen College of Engineering neering

Outline

  • Introduction
  • Composite Materials Design with Numerical

trical & Computer Engin

  • Composite Materials Design with Numerical

Mixing-Law

  • FDTD design of Nano-scale FSS
  • Stochastic analysis

Elect

2

  • Conclusions
slide-2
SLIDE 2

12/4/2010 2

University of Houston Cullen College of Engineering neering

Introduction

Electromagnetic Compatibility /Electromagnetic Interference trical & Computer Engin Elect

3

http://en.wikipedia.org/wiki/Electromagnetic_compatibility

University of Houston Cullen College of Engineering neering Shielding trical & Computer Engin Elect

4

NASA JSC Nanotube

http://research.jsc.nasa.gov/BiennialResearchReport/PDF/Eng-8.pdf http://www.ursi.org/Proceedings/ProcGA05/pdf/E01.3(01681).pdf

slide-3
SLIDE 3

12/4/2010 3

University of Houston Cullen College of Engineering neering

Shielding Materials with Periodic Structures

trical & Computer Engin

650nm Pitch 450nm length 100nm Arm s

Elect

5

FSS

University of Houston Cullen College of Engineering neering Composite Materials with Numerical Mixing-Law trical & Computer Engin Maxwell-Garnett equation Elect

6

q

∑ ∑

= =

+ − − + − + =

N k e k e k k N k e k e k k e

f f

1 1 e eff

2 1 2 3 ε ε ε ε ε ε ε ε ε ε ε

L L L L + + + − − + + + + + − + + − = + − ) 2 ( / ) ( ) ( 2 ) 2 ( ) 2 ( / ) ( ) 2 ( ) ( 2

1 2 3 1 3 2 1 2 1 1 1 3 1 3 2 1 1 1 eff eff

ε ε ε ε ε ε ε ε ε ε ε ε ε ε ε ε ε ε ε ε a a a a f

e e e e e e e e

slide-4
SLIDE 4

12/4/2010 4

University of Houston Cullen College of Engineering neering Composite Materials with Numerical Mixing-Law trical & Computer Engin

Periodic Composite Material 3D Periodical Structure

Elect

7

Unit Element

University of Houston Cullen College of Engineering neering

FDM FDM

Zmax V PBC V trical & Computer Engin Xmax Ymax Voltage PBC PBC PBC Voltage Voltage Voltage PBC: periodic boundary condition Elect

8

PBC: periodic boundary condition

Unit Element

slide-5
SLIDE 5

12/4/2010 5

University of Houston Cullen College of Engineering neering

Effect of Effect of I Inclusion nclusion E Electrical lectrical P Properties roperties

trical & Computer Engin Elect

9 _ _ _ _

=1cm, 1, 3, 1

r host rx inclusion ry inclusion rz inclusion

ε ε ε ε = = = = l

University of Houston Cullen College of Engineering neering

Effect of Effect of I Inclusion nclusion E Electrical lectrical P Properties roperties

trical & Computer Engin Elect

10 _ _ _ _

=1cm, 1, 3 , 1

r host rx inclusion ry inclusion rz inclusion

ε ε ε ε = = = = l

_ _ _ _

=1cm, 1, 30, 3

r host rx inclusion ry inclusion rz inclusion

ε ε ε ε = = = = l 30 3

slide-6
SLIDE 6

12/4/2010 6

University of Houston Cullen College of Engineering neering trical & Computer Engin Elect

11

University of Houston Cullen College of Engineering neering

Electrical Properties Versus Frequency Electrical Properties Versus Frequency

trical & Computer Engin Elect

12

slide-7
SLIDE 7

12/4/2010 7

University of Houston Cullen College of Engineering neering

Cubic Inclusion Cubic Inclusion

trical & Computer Engin Vf=0.764 Elect

13

University of Houston Cullen College of Engineering neering

Numerical Results Versus Frequency Numerical Results Versus Frequency

trical & Computer Engin

(1 )

eff incl f host f

V V ε ε ε = × + × −

Elect

14

slide-8
SLIDE 8

12/4/2010 8

University of Houston Cullen College of Engineering neering Optimization for Multiphase Mixture trical & Computer Engin volume fractions 0.1. espr: 4.42 and 4.04 Elect

15

Conductivity : 0.00715 S/m and 1.369 S/m University of Houston Cullen College of Engineering neering

FDTD Modeling

trical & Computer Engin Elect

16

slide-9
SLIDE 9

12/4/2010 9

University of Houston Cullen College of Engineering neering

Challenge in the Modeling of IR FSSs Challenge in the Modeling of IR FSSs

  • PEC assumption is not valid any

more

  • The metal is highly frequency-

dependent no

Lorentz-Drude Model (gold) ( ) ( ) ( )

f b

ε ω ε ω ε ω +

trical & Computer Engin

dependent now

  • It has both negative permittivity

and conductive loss

  • But all of the tradition microwave

designs are based on this assumption.

  • FDTD Modeling for Periodic

Structures

( ) ( ) ( ) ( )

( )

2 2 2 2 1

= 1

f r r r k p i p i i i

f j j ε ω ε ω ε ω ω ω ω ω ω ω

=

= + ⎛ ⎞ ⎛ ⎞ Ω ⎜ ⎟ − + ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ − Γ − + Γ ⎝ ⎠ ⎝ ⎠

Elect

17

University of Houston Cullen College of Engineering neering

Finite-Sized Electromagnetic Source

?

PBC PBC

?

z z k0 kx kz

“Brute Force”

trical & Computer Engin

Plane Wave Incidence Finite Size Source Incidence Plane wave incidence

a a x x ( ) ( ) ( ) ( )

, ,

x

jk a

E x a H x a E x H x e− ⎡ ⎤ ⎡ ⎤ + + = ⎣ ⎦ ⎣ ⎦

…. …. …. ….

Elect

18

  • Periodic boundary condition (PBC) can be applied for above

equation in both frequency and time domains

Finite size source incidence

Assumption is no longer valid

  • In time domain, “Brute Force” FDTD simulation is needed
slide-10
SLIDE 10

12/4/2010 10

University of Houston Cullen College of Engineering neering

ASM-FDTD Method

z z

( , , ) th unit cell

n tot x

na y t n + E

trical & Computer Engin

….

a x

….

a x

…. ….

Spectral domain transformation of the finite size source The field in the 0th unit cell excited by this source can be represented as ( ) ( ) ( ) ( )

, , ,

x

jk na i i x n

J x y k J x y x x na y y e δ δ

∞ − =−∞

= − − −

% ( ) ( )

/ /

, 2

a i i x x a

a J x y J k dk

π π

π − =

%

/

( ) ( )

a

a k dk

π

E E

n=0 n=1 n=2 n=-2 n=-1

Elect

19

  • The 0th unit cell spectral domain solution can be obtained by

FDTD simulation when following PBC is applied

  • The field in nth cell is found by

/

( , , ) ( , , ) 2

tot tot x x a

x y t k y t dk

π

π − =

E E

( , , ) (0, , )

x

jk a tot tot

a y t y t e− = E E

0 (

, , )

tot x

k y t E

/ /

( , , ) ( , , ) 2

x

a jk na n tot tot x x a

a x na y t k y t e dk

π π

π

− −

+ =

E E

University of Houston Cullen College of Engineering neering

Numerical Examples

trical & Computer Engin

dipole source 45 mm above the structure

Elect

20

dipole source 45 mm above the structure dipole strip 12 mm by 3 mm periodicity 15 mm in both directions field sampled 90 mm under the dipole source

slide-11
SLIDE 11

12/4/2010 11

University of Houston Cullen College of Engineering neering

Different Model Impacts On The Final Result Different Model Impacts On The Final Result

trical & Computer Engin Elect

21

Simulated transmission intensity of an Au cross slot array on the quartz substrate (a=650nm, L=500nm, W=110nm, thickness=300nm and substrate dielectric constant is 2.1316

University of Houston Cullen College of Engineering neering

Three Practical Patterns I: Standard Case Three Practical Patterns I: Standard Case

trical & Computer Engin Elect

22

Simulated transmission intensity of an Au cross slot array on the quartz substrate (a=660nm, L=510nm, W=160nm, thickness=220nm and substrate dielectric constant is 2.1316

slide-12
SLIDE 12

12/4/2010 12

University of Houston Cullen College of Engineering neering

Five Practical Patterns II: Faced Centered Case Five Practical Patterns II: Faced Centered Case

trical & Computer Engin Elect

23

Simulated transmission intensity of an Au cross slot array on the quartz substrate (a=660nm, L=500nm, W=170nm, thickness=220nm and substrate dielectric constant is 2.1316

University of Houston Cullen College of Engineering neering

Variation quantification Composite mixtures have inherent randomness

trical & Computer Engin

The homogenized EM property needs to be evaluated Sources of randomness includes

  • - material electrical properties

Elect

24

  • - component geometry
slide-13
SLIDE 13

12/4/2010 13

University of Houston Cullen College of Engineering neering

Variation quantification techniques Monte-Carlo (MC) : Simple to implement, computationally

expensive

trical & Computer Engin

expensive

Perturbation: Limited to small fluctuation Stochastic collocation method (SCM): Can handle large

fluctuations highly efficient transforms the stochastic analysis into a

Elect

25

fluctuations, highly efficient, transforms the stochastic analysis into a series of deterministic simulations

University of Houston Cullen College of Engineering neering

Suppose we have N random parameters

– we use the abbreviation

Introduction to SCM

1

{ }N

n n

ξ

=1

{ }N

n n

ξ

=1

{ }N

n n

ξ

=

1

{ }N

n n

ξ

=

1 2

{ , ,..., }

N

ξ ξ ξ ξ = r

1 2

{ , ,..., }

N

ξ ξ ξ ξ = r

( )

( )

n n

ρ ξ

trical & Computer Engin – the parameters could be distributed according to a joint PDF – each could be distributed independently according to its probability density function (PDF)

1

{ }N

n n

ξ

=

( )

ρ ξ r

n

ξ

( )

n n

ρ ξ

( )

1

( )

N n n n

ρ ξ ρ ξ

=

=∏ r

Elect

26

Realization = a output from the deterministic simulation tool for a specific choice of

( ) f ξ r

1

{ }N

n n

ξ ξ

=

= r

slide-14
SLIDE 14

12/4/2010 14

University of Houston Cullen College of Engineering neering

Introduction to SCM (cont.)

One may be interested in statistics of outputs

– average or expected value trical & Computer Engin – variance

[ ]

( ) ( ) E f f d ξ ρ ξ ξ

Γ

=∫ r r r [ ]

( )

2 2 2 2

[ ] [ ] ( ) ( ) ( ) ( ) Var f E f E f f d f d ξ ρ ξ ξ ξ ρ ξ ξ = −

∫ ∫

r r r r r r ( ( )) ( ) G f d ξ ρ ξ ξ

Γ

r r r

Elect

27

– higher moments

( )

( ) ( ) ( ) ( ) f d f d ξ ρ ξ ξ ξ ρ ξ ξ

Γ Γ

= −

∫ ∫

University of Houston Cullen College of Engineering neering

Integrals of the type

Introduction to SCM (cont.)

( ( )) ( ) G f d ξ ρ ξ ξ

Γ

r r r

trical & Computer Engin

cannot, in general, be evaluated exactly Thus, these integrals are approximated using a quadrature rule

1

( ( )) ( ) ( ) (( ( )))

Q q q q q

G f d G f ξ ρ ξ ξ ω ρ ξ ξ

Γ =

=∑

r r r r r

Elect

28

for some choice of quadrature weights and quadrature points

q

1

{ }Q

q q

ω

= 1

{ }Q

q q

ξ

=

r

slide-15
SLIDE 15

12/4/2010 15

University of Houston Cullen College of Engineering neering

Introduction to SCM (cont.)

To use such a rule, one needs to know the simulation output

at each of the quadrature points

  • - for this purpose, one can build a polynomial approximation

( ) f ξ r

{ }

1 Q q q

ξ

=

r

( ( )) G f ξ r % ( ( )) G f ξ r trical & Computer Engin and then evaluate that approximation at the quadrature points

  • - the simplest means of doing this is to use the set of Lagrange

interpolation polynomials corresponding to the sample points

( ) { }

1

sample

N s s

L ξ

=

r

( ) ( )

( )

( ) ( )

sample

N s s

G f G f L ξ ξ ξ = ∑ r r %

{ }

1

sample

N s s

ξ

=

r

1

( ( )) ( ) ( ) (( ( )))

Q q q q q

G f d G f ξ ρ ξ ξ ω ρ ξ ξ

Γ =

=

∫ ∑

r r r r r

Elect

29

  • Then we have the approximation for the mean:

( ) ( )

( )

1

( ) ( )

s s s

G f G f L ξ ξ ξ

=

[ ]

( )

1 1

( ) ( ) ( ) ( )

sample

N Q s q s q q s q

E f f d f L ξ ρ ξ ξ ξ ω ρ ξ ξ

Γ = =

= = ∑

∑ ∫

r r r r r r ( )

1 1 1

( ( )) ( )

sample

q N Q s q s q q s q

G f L ξ ω ρ ξ ξ

= = =

= ∑

r r r

University of Houston Cullen College of Engineering neering

i

' ε

_

= 1.0 m,

1, = =

r host host x y z

μ

ε σ = =

l l l

Numerical example

trical & Computer Engin

_

[ ]:15% [ 8, [

] ] 1

r inclusion

inclusion

E VF E E ε σ =

=

Mean values:

Goal: Evaluating the global sensitivity of effective permittivity

Elect

30

g g y p y due to variation in certain mixing parameters

Varying parameters: inclusion relative permittivity , inclusion conductivity

and volume fraction; all of which are assumed to be Gaussian variables.

slide-16
SLIDE 16

12/4/2010 16

University of Houston Cullen College of Engineering neering

Convergence Test

  • Operating Frequency: 1e10 Hz
  • Variable: inclusion permittivity, has a variance around its mean value

30% ±

trical & Computer Engin Elect

31

mean Standard deviation University of Houston Cullen College of Engineering neering

The Effect of Inclusion Relative Permittivity Variation

trical & Computer Engin

_

8

r inclusion

ε =

Elect

32

This case involves a variance of 30% around the mean value

_

8

r inclusion

ε =

_ host c r host

j σ ε ε ωε = −

slide-17
SLIDE 17

12/4/2010 17

University of Houston Cullen College of Engineering neering

The Effect of Inclusion Conductivity Variation

trical & Computer Engin

σ

Elect

33

This case involves a variance of 30% around the mean value

1

inclusion

σ

=

_ inclusion c r inclusion

j σ ε ε ωε = −

University of Houston Cullen College of Engineering neering

The Effect of Volume Fraction Variation

trical & Computer Engin Elect

34

This case involves a variance of 30% around the mean value %

15% VF =

slide-18
SLIDE 18

12/4/2010 18

University of Houston Cullen College of Engineering neering

The Effect of Volume Fraction Variation

trical & Computer Engin Elect

35

This case involves a variance of 30% around the mean value %

15% VF =

University of Houston Cullen College of Engineering neering

Conclusion

  • New FDTD methods for periodic structures

trical & Computer Engin

  • Applications include

– FSS – Meta-materials – Nano-scale devices

  • Multi layer periodic structures at different

Elect

36

  • Multi-layer periodic structures at different

periodicities