Capacitance matrix and integral equation method

Click For Summary
SUMMARY

The discussion focuses on calculating the capacitance matrix using the integral equation method, specifically the method of moments, for two infinitely long and thin parallel strips in vacuum. The potential of a line charge is expressed through a Green's function for a two-dimensional Laplace equation, leading to a set of linear equations represented in matrix form as IN = C-1Q. The main challenge presented is extracting the total capacitance C from the capacitance matrix Cij, which represents the capacitance between segments and self-capacitance, and reconciling this with textbook definitions.

PREREQUISITES
  • Integral equation methods in electromagnetism
  • Green's functions for solving Laplace equations
  • Matrix algebra and linear equations
  • Capacitance concepts and calculations
NEXT STEPS
  • Study the method of moments for capacitance calculations
  • Learn about Green's functions in electrostatics
  • Explore matrix inversion techniques in numerical analysis
  • Research self-capacitance and mutual capacitance in parallel conductors
USEFUL FOR

Electrical engineers, physicists, and researchers involved in electromagnetic theory and capacitance calculations will benefit from this discussion.

Jack the Stri
Messages
39
Reaction score
0
Dear all,

I'm having some trouble calculating the capacitance matrix, as outlined below.

So first of all, the lecture notes I'm using use an integral equation method (method of moments) to determine the capacitance of two infinitely long and thin parallel strips in vacuum, at a distance d of each other. The z-axis is taken along one of both strips, and image theory teaches us that this problem is equivalent to one strip at a distance d/2 of a PEC plane.

The notes start from the potential of a line charge, and use the superposition principle to express the potential of the surface charge density, yielding

\phi(\rho) = \int^{}_{c}\rho_{l}(\rho')G(\rho,\rho')dc'

With the position vectors in bold (rho and rho' being each other's image), \rho_{l} the surface charge density, c the strip width and G the Green's function for a two-dimensional Laplace equation in a homogeneous half-space on top of a PEC plane.

The next step is to impose the boundary condition that the strip surface is equipotential, e.g. at 1. Hence

lim_{\rho\rightarrow c} \phi(\rho) = lim_{\rho\rightarrow c} \int^{}_{c}\rho_{l}(\rho')G(\rho,\rho')dc'=1

This is solved by subdividing the strip into N segments of length \Delta, allowing the above expression to be rewritten as
lim_{\rho\rightarrow c}\sum^{N}_{i=1}\rho_{li}\int^{}_{c_{i}}G(\rho,\rho')dc'_{i}=1

The text then goes on to explain that, when one chooses the centres \rho_{j} as the limit c to approach, the sum forms a set of N linear equations with N unknown coefficients. In matrix form this is written:

I_{N}=C^{-1}Q

With IN an Nx1 column matrix with elements 1, Q the Nx1 column matrix with the elementary line charge densities (or at least, the average of the segment, the way I see it), and C-1 the inverted capacitance matrix.

The values of Cij are then approximated to fill in the matrix.

However.

I understand the steps leading up to this, but in this particular case, the only output for the capacitance is one number, not an NxN matrix (which makes sense otherwise the matrix size would be dependent upon the accuracy of your calculation). The question then is, how do you get the capacitance C from this matrix, and, in the more general case of a capacitance matrix, how do you extract each individual Cij (since these are not the same as the Cij mentioned above)?

I can post some more elaboration on the maths followed if necessary.

Edit: LaTeX doesn't seem to be co-operating with me as far as formatting certain formula parts in bold, sorry...
 
Last edited:
Physics news on Phys.org
Bumpity bump. From what I can make out of it each element i,j of the matrix represents the capacitance between two segments \Delta_{i} and \Delta_{j}, or the self-capacitance. The capacitance value would then be the sum of the capacitances for each individual segment, i.e. C = \sum C_{i} = \sum (C_{ii} - \sum_{j(j \neq i)} |C_{ij}|), or the diagonal minus the absolute values of the rest.

This doesn't match the textbook though. Any expertise around here?
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
976
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 9 ·
Replies
9
Views
4K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 6 ·
Replies
6
Views
2K
Replies
3
Views
3K
  • · Replies 5 ·
Replies
5
Views
1K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
1K