Electric Sandwich :: Try Open Physics Challenge !

Given two identical Solid uniformly charged (ρ\rho , RR) hemisphere , initially system in equilibrium Find time period of small vertical oscillation of charged particle qq and mass mm.

Details :

d=2m d=2m

R=1mR=1m

*hemisphere are fixed

Take It as challenge , and try to solve it. And then Share your all approaches with us!

#Calculus #Physics #ElectricityAndMagnetism #Challenge #Time_period

Note by Deepanshu Gupta
6 years, 2 months ago

No vote yet
1 vote

  Easy Math Editor

This discussion board is a place to discuss our Daily Challenges and the math and science related to those challenges. Explanations are more than just a solution — they should explain the steps and thinking strategies that you used to obtain the solution. Comments should further the discussion of math and science.

When posting on Brilliant:

  • Use the emojis to react to an explanation, whether you're congratulating a job well done , or just really confused .
  • Ask specific questions about the challenge or the steps in somebody's explanation. Well-posed questions can add a lot to the discussion, but posting "I don't understand!" doesn't help anyone.
  • Try to contribute something new to the discussion, whether it is an extension, generalization or other idea related to the challenge.
  • Stay on topic — we're all here to learn more about math and science, not to hear about your favorite get-rich-quick scheme or current world events.

MarkdownAppears as
*italics* or _italics_ italics
**bold** or __bold__ bold

- bulleted
- list

  • bulleted
  • list

1. numbered
2. list

  1. numbered
  2. list
Note: you must add a full line of space before and after lists for them to show up correctly
paragraph 1

paragraph 2

paragraph 1

paragraph 2

[example link](https://brilliant.org)example link
> This is a quote
This is a quote
    # I indented these lines
    # 4 spaces, and now they show
    # up as a code block.

    print "hello world"
# I indented these lines
# 4 spaces, and now they show
# up as a code block.

print "hello world"
MathAppears as
Remember to wrap math in \( ... \) or \[ ... \] to ensure proper formatting.
2 \times 3 2×3 2 \times 3
2^{34} 234 2^{34}
a_{i-1} ai1 a_{i-1}
\frac{2}{3} 23 \frac{2}{3}
\sqrt{2} 2 \sqrt{2}
\sum_{i=1}^3 i=13 \sum_{i=1}^3
\sin \theta sinθ \sin \theta
\boxed{123} 123 \boxed{123}

Comments

We will proceed by finding the potential at the axis of a solid hemisphere at a distance of aa from its center. We can do this by first finding potential due to a disk at it's axis and then integrating to get a sphere. I am simply giving formula here, and leaving the integration as an exercise to the reader.

V(a)V(a) is the potential at distance aa when R=1R=1 :

ϕ=ρ2ε0V(a)=ϕ[a23aa2+13a2+13a+13a+12]\phi =\frac { \rho }{ 2{ \varepsilon }_{ 0 } } \\ V(a)=\phi \left[ \frac { { a }^{ 2 } }{ 3 } -\frac { a\sqrt { { a }^{ 2 }+1 } }{ 3 } -\frac { \sqrt { { a }^{ 2 }+1 } }{ 3a } +\frac { 1 }{ 3a } +\frac { 1 }{ 2 } \right]

Now we know that the motion of this particle can be approximated to S.H.M. We will use the fact that total energy remains constant:

K.E of particle+ electrostatic P.E of particle= constant\text{K.E of particle+ electrostatic P.E of particle= constant}

mv22+qV(a+x)+qV(ax)=K\Rightarrow \frac { m{ v }^{ 2 } }{ 2 } +qV(a+x)+qV(a-x)=K

Since R.H.S is constant w.r.t time, we will differentiate to get:

mvdvdt+qdV(a+x)dx×dxdt+qdV(ax)dx×dxdt=0mv\frac { dv }{ dt } +q\frac { dV(a+x) }{ dx } \times \frac { dx }{ dt } +q\frac { dV(a-x) }{ dx } \times \frac { dx }{ dt } =0

mdvdt=qdV(a+x)dxqdV(ax)dx\Rightarrow m\frac { dv }{ dt } =-q\frac { dV(a+x) }{ dx }-q\frac { dV(a-x) }{ dx }

We can now find the R.H.S of above equation by differentiating the V(a)V(a) function.

dV(a+x)dx=ϕ3[2(a+x)2+1+(a+x)2+1(a+x)21(a+x)2+2(a+x)]\frac { dV(a+x) }{ dx } =\frac { \phi }{ 3 } \left[ -2\sqrt { { (a+x) }^{ 2 }+1 } +\frac { \sqrt { { (a }+x)^{ 2 }+1 } }{ { (a+x) }^{ 2 } } -\frac { 1 }{ { (a+x) }^{ 2 } } +2(a+x) \right]

In the above equation, we use the binomial approximation on each and every term (even the terms under the square root) and then neglect all terms where power of xx is more than 11. I'm skipping the simplification, and giving the final answer at a=1a=1:

dV(a+x)dx=ϕ3[12+(852)2x]\frac { dV(a+x) }{ dx } =\frac { \phi }{ 3 } \left[ 1-\sqrt { 2 } +\frac { (8-5\sqrt { 2 } ) }{ 2 } x \right]

Hence: dV(ax)dx=ϕ3[12(852)2x]\frac { dV(a-x) }{ dx } =\quad -\frac { \phi }{ 3 } \left[ 1-\sqrt { 2 } -\frac { (8-5\sqrt { 2 } ) }{ 2 } x \right]

Therefore:

mdvdt=qdV(a+x)dxqdV(ax)dx=ϕ3(852)qx\Rightarrow m\frac { dv }{ dt } =-q\frac { dV(a+x) }{ dx }-q\frac { dV(a-x) }{ dx } =-\frac { \phi }{ 3 } (8-5\sqrt { 2 } )qx

mdvdt=ρq6ε0(852)qx m\frac { dv }{ dt }=\frac { -\rho q }{ 6{ \varepsilon }_{ 0 } } (8-5\sqrt { 2 } )qx

Finally, we have the expression of time period as:

2π6mε0ρq(852)2\pi \sqrt { \frac { 6m{ \varepsilon }_{ 0 } }{ \rho q(8-5\sqrt { 2 } ) } }

Raghav Vaidyanathan - 6 years, 2 months ago

Log in to reply

Log in to reply

I will tell you what actually happened when I tried to solve this problem, I made the integral for electric field and in the integral itself I just added a little xx with rr and after series of approximations calculated the result.

But that doesn't seem elegant, hence I thought of differentiation method.

Ronak Agarwal - 6 years, 2 months ago

Log in to reply

@Ronak Agarwal That is acceptable, that is basically what you would do with differentiation as well,

as zero order term of field is 0 (since it is equillibrium point), and 1st order is dominant and restoring, and finally second order and higher orders are neglected. So, approximation here is just as rigorious as differentiation

@Raghav Vaidyanathan - nice, atleast you confirmed the answer, and yours seems simpler though there are some complex derivatives to be done, and maybe we should expect it to be bit complicated since this is not some mathematically simple problem anyway :)

Mvs Saketh - 6 years, 2 months ago

@Ronak Agarwal Okay. I actually found a lot of similarity between the method above and yours. After posting this, I tried using Newton-Leibniz rule in order to calculate the derivative of V(a)V(a) directly instead of getting V(a)V(a) and then differentiating. It does work. But I'm quite unsatisfied with the fact that this problem takes a while to solve, no matter what method one uses.

Raghav Vaidyanathan - 6 years, 2 months ago

Nice Raghav ! Now I realise finding field alway's is not good idea ... atleast for me , who know's very li'l calculus ..! You did it nicely bro ! Thanks for sharing !

Deepanshu Gupta - 6 years, 2 months ago

BullShit BullShit

Firstly I am using a simple result that, electric field due to a disc of radius rr having a surface charge density σ \sigma at on point on the axis at a distance zz from the centre of the disc is given by :

E=σ2ε0(1zz2+r2) \displaystyle E = \dfrac { \sigma }{ 2{ \varepsilon }_{ 0 } } (1-\frac { z }{ \sqrt { { z }^{ 2 }+{ r }^{ 2 } } } )

In my diagram I have taken a angular elemental disc whose radius is Rsin(θ)R\sin(\theta) and width of disc is given by =Rsin(θ)dθR\sin(\theta)d\theta, also to note that here σ=ρ×(Width)\sigma = \rho \times (Width)

Now electric field due to that element at a point that is at a distance rr from the centre of the hemisphere is given by :

dE=ρRsinθ2ε0(1(r+Rcosθ)(r+Rcosθ)2+(Rsinθ)2)dθ\displaystyle dE = \frac { \rho R\sin { \theta } }{ 2{ \varepsilon }_{ 0 } } (1-\frac { (r+R\cos { \theta ) } }{ \sqrt { ({ r+R\cos { \theta } ) }^{ 2 }+{ (R\sin { \theta } ) }^{ 2 } } } )d\theta

Integrating it our we have :

E(r)=ρR2ε00π/2(sinθ(r+Rcosθ)sinθ(r+Rcosθ)2+(Rsinθ)2)dθ\displaystyle E(r) = \frac { \rho R }{ 2{ \varepsilon }_{ 0 } } \int _{ 0 }^{ \pi /2 }{ (\sin { \theta } -\frac { (r+R\cos { \theta ) } \sin { \theta } }{ \sqrt { ({ r+R\cos { \theta } ) }^{ 2 }+{ (R\sin { \theta } ) }^{ 2 } } } )d\theta }

Now we know that for small xx we have :

E(r+x)=E(r)+xE(r) E(r+x)=E(r)+xE'(r)

and

E(rx)=E(r)xE(r) E(r-x)=E(r)-xE'(r)

We have net electric field due to both hemisphere's when particle is at a distance xx from it's mean position is given by :

E=2xE(r) E = -2xE'(r)

Now by newton leibnitz rule we have :

E(r)=ρ2ε00π/2(R3sin3θ((r+Rcosθ)2+(Rsinθ)2)3/2)dθ\displaystyle E'(r) = \frac { \rho }{ 2{ \varepsilon }_{ 0 } } \int _{ 0 }^{ \pi /2 }{ (\frac { { R }^{ 3 }\sin ^{ 3 }{ \theta } }{ { (({ r+R\cos { \theta } ) }^{ 2 }+{ (R\sin { \theta } ) }^{ 2 }) }^{ 3/2 } } )d\theta }

Evaluting the derivative at RR , (since the mean position is at a distance R from the centre of the two hemispheres.)

E(R)=ρ2ε00π/2(sin3θ2)dθ=ρ12ε0(852) \displaystyle E'(R) = \frac { \rho }{ 2{ \varepsilon }_{ 0 } } \int _{ 0 }^{ \pi /2 }{ (\sin ^{ 3 }{ \frac { \theta }{ 2 } } )d\theta } = \frac { \rho }{ 12{ \varepsilon }_{ 0 } } (8-5\sqrt { 2 } )

Hence acceleration of the particle is given by :

a=ρq6mε0(852)x \displaystyle a = -\dfrac { \rho q }{ 6m{ \varepsilon }_{ 0 } } (8-5\sqrt { 2 } )x

Hence the time period is given by :

T=2π6mε0ρq(852) \displaystyle T = 2\pi \sqrt { \dfrac { 6m{ \varepsilon }_{ 0 } }{ \rho q(8-5\sqrt { 2 } ) } }

I have jumped many steps it seems hence if you have any confusion ask in the comments.

Ronak Agarwal - 6 years, 2 months ago

Log in to reply

i see you have not used wolfram, Nice +1

Mvs Saketh - 6 years, 2 months ago

Log in to reply

did you solved it ? @Mvs Saketh

Deepanshu Gupta - 6 years, 2 months ago

Log in to reply

@Deepanshu Gupta i got to the expression of field for hemisphere and then thought it can be solved so i stopped

Mvs Saketh - 6 years, 2 months ago

Log in to reply

@Mvs Saketh Please post how much you reach .... Ronak and me discuss this earliear that this is very difficult to solve , if we use conventional method ... Roonak approach is very cool ! That's why I posted this as a note !

Deepanshu Gupta - 6 years, 2 months ago

Log in to reply

@Deepanshu Gupta My approach is no different from Ronaks, i too use Newton Leibinz method to avoid Dauting Integrals or wolfram where needed, i dont think i have anything more to add to him , But i shall add if i do get some idea :)

Mvs Saketh - 6 years, 2 months ago

Log in to reply

@Mvs Saketh Please Tell what should we do after this ? Please suggest , if want to do in different manner as Ronak awesomely did,....

E(r)=ρ2ε0π2cosθ(r+sinθ)cosθ(r)2+1+2rcosθdθ\displaystyle{E(r)=\cfrac { \rho }{ 2\varepsilon } { \int _{ 0 }^{ \cfrac { \pi }{ 2 } }{ \cos { \theta } -\cfrac { (r+\sin { \theta } )\cos { \theta } }{ \sqrt { { (r) }^{ 2 }+1+2r\cos { \theta } } } d\theta } }}

I had puted the value R=1 .

Deepanshu Gupta - 6 years, 2 months ago

Log in to reply

@Deepanshu Gupta Well, you can use WolfRam, though i know not if it can do it,

But,i dont think i can solve it, (not in a simple manner surely)

Mvs Saketh - 6 years, 2 months ago

Log in to reply

@Mvs Saketh I think there is an easier way to do this problem rather than finding field. Why not find potential energy? Then we can use total energy equation and differentiate it to get SHM equation.

Raghav Vaidyanathan - 6 years, 2 months ago

@Deepanshu Gupta I think there is an easier way to do this problem rather than finding field. Why not find potential energy? Then we can use total energy equation and differentiate it to get SHM equation.

Raghav Vaidyanathan - 6 years, 2 months ago

Log in to reply

@Raghav Vaidyanathan Have you done it with that ....?

Deepanshu Gupta - 6 years, 2 months ago

@Raghav Vaidyanathan yes you can, but whether it is really easy is not has to be tried, we do not yet know if the integral for potential of hemisphere is simpler (though it usually is the case that potential is simpler than field) (have you tried that?)

Mvs Saketh - 6 years, 2 months ago

Log in to reply

@Mvs Saketh It is simpler, but I think I am going wrong in some the integration somewhere. I don't think I'm getting the boundary conditions right when I substitute them in my equation

Raghav Vaidyanathan - 6 years, 2 months ago

Log in to reply

@Raghav Vaidyanathan You can post important expression here ..... The purpose of this note is exatly that ... to discuss various Possiblities ! @Raghav Vaidyanathan Please post what'sever you get ... may be we together complete it ... You can skips unnecceserry steps for posting here to save ur time ..!

Deepanshu Gupta - 6 years, 2 months ago

@Mvs Saketh Yes @Mvs Saketh @Deepanshu Gupta @Ronak Agarwal

I have tried it using the potential method and have gotten the right answer. I will post an overview soon.

Raghav Vaidyanathan - 6 years, 2 months ago

Log in to reply

@Raghav Vaidyanathan we are waiting !

Deepanshu Gupta - 6 years, 2 months ago

I think , u are decide to not to reply me ... am i right ? :) ... I have recently mentioned you 2-3 times but you are not responding ...

Deepanshu Gupta - 6 years, 2 months ago

Log in to reply

@Deepanshu Gupta Woah bro, Chill not at all, i was solving your cylinder problem and hence did not reply ,

Also, i do not get notifications if some one mentions me ,

Mvs Saketh - 6 years, 2 months ago

Log in to reply

@Mvs Saketh mention and disput's are are comes on Emails only ... But okay leave that , I actually replied to your comment on My problem Fearless boy..... and where I asked something please reply there ...!

Deepanshu Gupta - 6 years, 2 months ago
×

Problem Loading...

Note Loading...

Set Loading...