Capacitance method of moments - PowerPoint PPT Presentation

slide1 n.
Download
Skip this Video
Loading SlideShow in 5 Seconds..
Capacitance method of moments PowerPoint Presentation
Download Presentation
Capacitance method of moments

play fullscreen
1 / 56
Capacitance method of moments
269 Views
Download Presentation
connor-byers
Download Presentation

Capacitance method of moments

- - - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript

  1. Capacitance method of moments

  2. Q a Capacitance of hollow sphere

  3. Capacitance of parallel plates

  4. Capacitance of a circular disk r R   Voltage at origin   = 0 a 2 - 

  5. Capacitance of a circular disk r  a R  2 -  Charge density is uniform s0

  6. Capacitance of a circular disk r  a R  2 -  Charge density decreases s0(1-r/a)

  7. This technique was developed in a Ph.D. thesis by Ken Mei in 1962. His adviser was Jean van Bladel.

  8. Method of moments • a known charge causes a potential • what is the potential? • a measured or known potential is caused by an unknown charge • what was that charge? • used in em to find current distributions on antennas, scattering objects such as cubes, sleds, planes, and missles • PhD theses, faculty publications

  9. V(a) V(b) Q2 Q1 2 equations 2 unknowns

  10. V(a) V(b) Q2 Q1

  11. MATLAB!!

  12. >> A = [1/4 1/5 ; 1/5 1/4] ; >> V = [1 ; 1] ; >> Q = V’ / A % ‘transpose Q = 2.2222 2.2222 >> QQ = A \ V; QQ = 2.2222 2.2222

  13. complications??? r1 r2 • The potential @ each charge sphere is specified as V1 or V2 • Locations are with respect to an origin.

  14. solution to singularity Assume that the potential is uniform within the sphere and it is equal to the value at its edge r = a.

  15. a = 1

  16. 1 1 1 1 1 1

  17. Q = • 4.69 2.56 4.69 -4.69 -2.56 -4.69 • >>

  18. equal charge

  19. equal charge

  20. a b -a -a a -b -b b Capacitor C = Q / V rj ri

  21. b -b -b b Capacitor C = Q / V Harrington – p. 27 Dwight 200.01 & 731.2

  22. b -b -b b Capacitor C = Q / V

  23. MATLAB • clear; clf; • N=9; • az=37.5; • el=5;

  24. %identify subareas • for m = 1:2 • for h = 1:N • for k = 1:N • a(m - 1)*N*N + (h - 1)*N + k, :)=[m, h, k]; • end • end • end

  25. %calculate matrix elements • for h=1:2*N*N • for k=1:2*N*N • aa=norm(a(h, :) - a(k, :)); • if aa==0 • b(h, k) = 2 * sqrt(pi); • else • b(h, k) =1 / aa; • end • end • end

  26. %set voltages on plates • for h=1:N*N • V(h) = +1/2; • V(h+N*N) = -1/2; • end

  27. %calculate charges • Q = V*inv(b); • %top plate • QA(1:N*N) = Q(1:N*N); • %bottom plate • QB(1:N*N)=Q(N*N+1:2*N*N);

  28. %plot • [x,y]=meshgrid(1:N); • for i=1:N • za(1:N,i)=QA((i-1)*N+1:(i-1)*N+N)'; • zb(1:N,i)=QB((i-1)*N+1:(i-1)*N+N)'; • end • mesh(x,y,za) • hold on • mesh(x,y,zb)

  29. QT=0; for j=1:N*N QT=QT+Q(j); end QT

  30. What have you done? You have learned a technique, to accurately calculate the capacitance of a parallel plate capacitor! Big deal!

  31. Capacitance of a circular disk r  a R  2 -  Charge density is uniform s0

  32. Capacitance of a circular disk r  a R  2 -  Charge density is noniform

  33. C = 3.9420 C = 3.2367 C = 2.7026 C = 3.6332

  34. capacitance of a unit cube

  35. capacitance of a unit cube

  36. capacitance of a unit cube Hwang & Mascogni – “Electrical capacitance of a unit cube” – Journal of Applied Physics 3798-3802 (2004).

  37. 2a C a

  38. note the charge density at the corner

  39. a h Sava V. Savov February 1, 2004

  40. 1 2 3 x. • E = s /2 • V = x s /2 • V1 = 2 [0/2 s1 + 1 s2 + 2/2 s3 ] • V2 = 2 [1/2 s1 + 0 s2 + 1/2 s3] • V3 = 2 [2/2 s1 + 1 s2 + 0/2 s3]

  41. 1 2 3 • V1 0 1 2/2 s1 • V2 = 2 1/2 0 1/2 s2 • V3 2/2 1  0 s3 • pn junction • double layer • VLSI

  42. pn junction • linear quadratic • V[-2d] - 0.36 - 2 • V[-d] - 0.18 - 3/2 • V[0] = 0 or 0 • V[d] + 0.18 + 3/2 • V[2d] + 0.36 + 2

  43. matrix • 0 1 2 3 4/2 • 1/2 0 1 2 3/2 • 2/2 1 0 1 2/2 • 3/2 2 1 0 1/2 • 42 3 2 1 0

  44. MATLAB • [V] = [matrix] [Q] V V r r x x

  45. Calculate vector “V” • clear • clf • m=input('What is the size of the matrix? ...'); • for i=1:m • v(i)=NaN; • end • v

  46. Calculate matrix “a” • for i=1:m • for j=1:m • a(i,j)=NaN; • end • end

  47. for i=1:m • for j=i:m • if i==j • a(i,j)=0; • b(i)=i; • elseif i<j & j<m a(i,j)=(j-b(i)); elseif i<j & j==m • a(i,j)=(j-b(i))/2; • end • end • end