• Support PF! Buy your school textbooks, materials and every day products Here!

Cylindrical capacitor with varying dielectric

  • Thread starter brainslush
  • Start date
  • #1
26
0

Homework Statement



Consider a long cylindrical coaxial capacitor with an inner conductor of radius a, and outer conductor of radius b, and a dielectric with a relative electric permittivity or dielectric ε(r), varying with the cylindrical radius. The capacitor is charged to the voltage V. Choose the radial dependence of ε(r) such that the energy density in the capacitor is constant.
Calculate the electric field inside the capacitor.

Homework Equations


Energy density equation

The Attempt at a Solution


from the energy density eq. which is defined by [itex]u = \frac{1}{2} \textbf{E} \circ \textbf{D}[/itex] we get that
1/ε(r) [itex]\propto \left| \textbf{E} \right|^{2}[/itex]

[itex]\Phi(b) - \Phi(a) = V[/itex]

Honestly I'm not quite sure which ansatz one can use. I was trying to solve it with the Gauss's law but didn't find a satisfying solution.
 

Answers and Replies

  • #2
TSny
Homework Helper
Gold Member
12,537
2,956
You might try using the basic equation [itex]\cdot[/itex]D = [itex]\cdot[/itex](εE) = [itex]\rho[/itex] along with the energy condition ε(r)E2 = constant to get a differential equation for ε(r).

ρ represents the free charge density which would be zero between the plates.
 
  • #3
gabbagabbahey
Homework Helper
Gold Member
5,002
6
I think Gauss' Law for dielectrics is probably the easiest way to go here. You say you did not get a satisfying answer....can you post your attempt so we can see where you may be going wrong?
 
  • #4
TSny
Homework Helper
Gold Member
12,537
2,956
brainslush,

Take gabbagabbahey's advice! My suggestion will get you there, but not as easily.:redface:
 
  • #5
26
0
[itex]\nabla \circ \textbf{D} = \epsilon (\nabla \circ \textbf{E}) + (\nabla \circ \epsilon) \textbf{E} = 0[/itex]

Now this gives the diferential eq. (we only care about the radial direction)

[itex]\frac{\frac{\partial}{\partial \rho} \epsilon (r)}{\epsilon (r)} = -\frac{\frac{\partial}{\partial r} E(r)}{E(r)}[/itex]

which then resolves into

[itex]\epsilon = \frac{c}{\left|E\right|}[/itex]

Which somewhat makes no sense to me, since this would imply that E is constant and therefore epsilon is constant. I probably messed something up.
 
  • #6
ehild
Homework Helper
15,497
1,878
[itex]\nabla \circ \textbf{D} = \epsilon (\nabla \circ \textbf{E}) + (\nabla \circ \epsilon) \textbf{E} = 0[/itex]

Now this gives the diferential eq. (we only care about the radial direction)

[itex]\frac{\frac{\partial}{\partial \rho} \epsilon (r)}{\epsilon (r)} = -\frac{\frac{\partial}{\partial r} E(r)}{E(r)}[/itex]
Your equation is wrong. Check out divergence in cylindrical coordinates.
http://en.wikipedia.org/wiki/Divergence

Anyway, it would be easier to use Gauss' Theorem for the displacement vector D. Edit: in integral form...

ehild
 
Last edited:
  • #7
TSny
Homework Helper
Gold Member
12,537
2,956
ε
[itex]\nabla \circ \textbf{D} = \epsilon (\nabla \circ \textbf{E}) + (\nabla \circ \epsilon) \textbf{E} = 0[/itex]

Now this gives the diferential eq. (we only care about the radial direction)

[itex]\frac{\frac{\partial}{\partial \rho} \epsilon (r)}{\epsilon (r)} = -\frac{\frac{\partial}{\partial r} E(r)}{E(r)}[/itex]

which then resolves into

[itex]\epsilon = \frac{c}{\left|E\right|}[/itex]

Which somewhat makes no sense to me, since this would imply that E is constant and therefore epsilon is constant. I probably messed something up.
Sorry for leading you off onto this path! If you really want to do it this way, then as echild pointed out you need to express the divergence in cylindrical coordinates. [(By the way, in the last term of your first equation you should have a gradiant of ε rather than a divergence (which wouldn't make sense).] Anyway, it is easier to first solve the differential equation ∇[itex]\cdot[/itex]D = 0 for D in the region between the plates and then relate E to D via the dielectric constant. You can then discover the functional form of ε by the requirement of constant energy density. This avoids having to find a differential equation for ε(r), although it's not particularly difficult to do so.

Of course ∇[itex]\cdot[/itex]D = ρ is just the differential form of Gauss' law for D. The integral form of Gauss' law for D will allow you to easily find the form of D between the plates, as others have suggested.
 
  • #8
gabbagabbahey
Homework Helper
Gold Member
5,002
6
Just to clarify, my suggestion was to use Gauss' Law in integral form. Using it in differential form is also very easy provided you are comfortable with taking the divergence of a function in cylindrical coordinates - if you just apply the divergence directly to the product [itex]\mathbf{D}=\epsilon(r)E(r)\hat{r}[/itex] without expanding it using the product rule, the differential equation you get will be trivial to solve.
 
  • #9
26
0
Of course...

[itex]\nabla \circ (\epsilon \textbf{E}) = (\nabla \epsilon) \circ \textbf{E} + \epsilon (\nabla \circ \textbf{E}) \Rightarrow [/itex]

[itex]\frac{\frac{\partial}{\partial \rho} \epsilon (\rho)}{\epsilon (\rho)} = -\frac{\frac{\partial}{\partial \rho} (\rho E(\rho))}{\rho E(\rho)} = -\frac{1}{\rho}-\frac{\frac{\partial}{\partial \rho} E(\rho)}{E(\rho)} \Rightarrow[/itex]

[itex]\epsilon = \frac{c_{0}}{\rho \left| E(\rho) \right|} [/itex] by using the condition [itex]\epsilon \propto \frac{1}{\left| \textbf{E} \right|^{2}}[/itex] we get

[itex]\epsilon = \left(\frac{c_{1}}{\rho}\right)^{2}[/itex]

So we can say that [itex] E \propto \rho[/itex]

then we can use [itex]\int_{a}^{b} E d\rho = V = \int_{a}^{b} c \rho d\rho[/itex]

to get [itex]c = \frac{2V}{b^{2}-a^{2}}[/itex]

[itex]E = \frac{2V \rho}{a^{2}-b^{2}}[/itex]
 
  • #10
TSny
Homework Helper
Gold Member
12,537
2,956
Of course...

[itex]\nabla \circ (\epsilon \textbf{E}) = (\nabla \epsilon) \circ \textbf{E} + \epsilon (\nabla \circ \textbf{E}) \Rightarrow [/itex]

[itex]\frac{\frac{\partial}{\partial \rho} \epsilon (\rho)}{\epsilon (\rho)} = -\frac{\frac{\partial}{\partial \rho} (\rho E(\rho))}{\rho E(\rho)} = -\frac{1}{\rho}-\frac{\frac{\partial}{\partial \rho} E(\rho)}{E(\rho)} \Rightarrow[/itex]

[itex]\epsilon = \frac{c_{0}}{\rho \left| E(\rho) \right|} [/itex] by using the condition [itex]\epsilon \propto \frac{1}{\left| \textbf{E} \right|^{2}}[/itex] we get

[itex]\epsilon = \left(\frac{c_{1}}{\rho}\right)^{2}[/itex]

So we can say that [itex] E \propto \rho[/itex]

then we can use [itex]\int_{a}^{b} E d\rho = V = \int_{a}^{b} c \rho d\rho[/itex]

to get [itex]c = \frac{2V}{b^{2}-a^{2}}[/itex]

[itex]E = \frac{2V \rho}{a^{2}-b^{2}}[/itex]
Looks good. Note that you could also write ∇[itex]\cdot[/itex](εE) = 0 as [itex]\frac{1}{ρ}\frac{d(ρεE)}{dρ}[/itex] = 0 which immediately gives ρεE = co
 
  • #11
26
0
Thanks for the help
 
  • #12
gabbagabbahey
Homework Helper
Gold Member
5,002
6
[itex]\int_{a}^{b} E d\rho = V = \int_{a}^{b} c \rho d\rho[/itex]
Careful with your signs :wink:
 

Related Threads on Cylindrical capacitor with varying dielectric

Replies
0
Views
2K
  • Last Post
Replies
10
Views
1K
  • Last Post
Replies
5
Views
8K
  • Last Post
Replies
1
Views
2K
  • Last Post
Replies
1
Views
2K
Replies
1
Views
3K
Replies
1
Views
807
Replies
1
Views
12K
Replies
2
Views
3K
Replies
7
Views
3K
Top