Monday, March 20, 2017

electricity - General integral to find resistance



My question is: is there a simple and truly general equation for the resistance between two electrical equipotential surfaces?. Obviously, if so, what is it, and if not, why? It would be very difficult to solve, granted, but I just want to see a calculus equation that is fully descriptive. I have two frameworks under which this could be entertained, I'll write those out and then explain the motivation.


To start with, we need propose that the volume separating the two surfaces has a volumetric resistivity, $\rho$ in units of $(\Omega m)$.



We can limit the discussion to a defined volume, then the surfaces reside in that volume or on the surface of it. This volume may have a constant resistivity $\rho$ while everywhere outside the volume is completely electrically insulating.



An alternative to the above approach that might make the task more or less difficult would be to replace a constant resistivity with a spatial dependence $\rho(\vec{r})$ and no longer require a boundary condition. In that case we only have 3 mathematical inputs to the problem, which is the resistivity defined for all $\vec{r}$ and a definition of the two surfaces, $S_1$ and $S_2$.



The basic algebraic formulation that I find insufficient is:


$$R = \rho \frac{\ell}{A}$$


Where $l$ is the length of the restive material that is any shape which has translation symmetry over that length, and $A$ is the cross-sectional area. Obviously, this is a rather simple equation that won't apply to more complicated geometry. Even more sophisticated academic sources seem to give equations that fall short of what I'm asking. For example:



$$R = \rho \int_0^l \frac{1}{A(x)} dx$$


I think it's obvious that an equation such as this is built upon a myriad of assumptions. For a thought experiment, imagine that the area starts out as very small and then pans out to very large quickly. Well, accounting for the larger area in the above sense underestimates the resistance, because the charge has to diffuse out perpendicular to the average direction of flow as well as parallel to it.


I have some reasons to suspect this might actually be rather difficult. A big reason is that all the approaches I'm familiar with require the flow paths to be established beforehand, which can't be done for what I'm asking. So maybe this will result in two interconnected calculus equations.



I had an interest in Squishy Circuits, and it occurred to me that I can't quickly and simply write down the equation for resistance between two points. The unique thing about Squishy Circuits is that it calls for two types of dough, one that conducts and one that is mostly insulating. However, the recipes aren't perfect and because of that, the young children who play with these circuits regularly encounter the limits of conductor and insulator definitions. If you make your conductor dough too long and/or too thin, you will encounter dimming of the light you connect with it. Similarly, a thin insulator layer will lead to a lot of leakage current which also dims the light.



Answer



I will do the case where the material is homogeneous and isotropic, $\rho = \sigma^{-1}$ is a constant proportional to the identity matrix. We are interested in the steady state, where none of our variables depend on time.


We have $\nabla \times E = 0$ from Faradays law and, $\nabla\cdot J = 0$ from the equation of continuity, where $J$ is the current density. Ohm's law tell us that $J=\sigma E$. Taking the divergence of Ohm's Law we get $\nabla \cdot E = 0$. Therefore in steady state $E$ must be a divergence-free, curl-free function.


This means that the potential $\phi$ ($\nabla\phi =E$) obeys Laplace's equation,


$\nabla^2 \phi = 0$.



To solve this we need the appropriate boundary conditions which are as follows. At the boundary where your resistor connects to a lead the potential $\phi$ must be the same value as the potential on the lead. At the boundary of your resistor not connected to a lead you can't have current flowing out, so the appropriate condition is $\nabla\phi\cdot n = 0$, where $n$ is the surface normal. This is sufficient to determine a single solution to Laplace's equation.


Laplace's equation is a very nice and friendly equation and there is a lot of material on numeric and analytic solutions available, although the mixed boundary conditions will be annoying.


Once you have your solution to Laplace's equation you need the total current. To get that, pick any cross sectional surface $S$ of your resistor, and using $J = \sigma E = \sigma\nabla\phi$, we get the total current I is equal to


$I = \sigma\int dS\cdot\nabla\phi$.


Then you can use $R = V/I$ to get an effective resistance. The case with a non homogeneous or non isotropic material is similar, you just end up with a different equation from Laplace's equation, which may be a little more annoying to solve.


I can't imagine a dough based system will need this level of precision though :). For anything dough based probably just assuming its a cylinder and using $R = \rho L /A$ will get you close enough no? I always thought of dough physics as really being an order of magnitude sort of game, like astronomy.


No comments:

Post a Comment

classical mechanics - Moment of a force about a given axis (Torque) - Scalar or vectorial?

I am studying Statics and saw that: The moment of a force about a given axis (or Torque) is defined by the equation: $M_X = (\vec r \times \...