_{1}

Two electrically charged rings of different sizes are assembled along their common vertical symmetry axis through their centers. The bottom ring is secured on a horizontal support while the top one is loose. For a set of practical values characterizing the charged rings we envision a scenario where the mutual electric repulsion between the rings and the weight of the top ring results in stable nonlinear oscillations. To quantify the characteristics of the oscillations, we utilize a Computer Algebra System specifically
*Mathematica* [1]. We accompany the analysis with a simulation for a comprehensive visual understanding.

Reviewing and instructing off the popular undergraduate [

The project as briefly outlined in the abstract concerns with the oscillations of a charged ring that is placed horizontally floating above another charged ring clamped on a horizontal table along their common vertical symmetry axis. Charging both rings positively (negatively) creates repulsive electrostatic force. The floating ring being massive gets pull down by its weight. Adjusting the relevant physical parameters makes it feasible to balance these two forces making the latter motionless yielding to the equilibrium along the vertical axis, and therefore, intuitively speaking dropping or pushing the floating charged ring above or below the equilibrium in the absentee of non-conservative agents such as air friction results in stable oscillations. One would also speculate because of the nature of the charged ring-ring interaction, the oscillation would be nonlinear. To quantify the oscillations conducive to the relevant quantities such as: frequency, period, and amplitude, one needs to establish the equation of motion. The latter requires a methodology, meaning either the mutual interactive force between the rings are to be considered or from energy point of view of their electrostatic potential energy. The former requires calculation of the electric field of the bottom ring at off axis region where the floating ring is at, where the latter calls for evaluation of associated electrostatic potential. Contemplating these two we favored the latter. Information cited in [

As shown in _{1} on the lower ring to a similarly characterized differential charged length on the upper ring, p_{2}. Intuitively speaking any physical quantity such as electrostatic potential energy comes about from summing a distance dependent expression between these two points, p 1 p 2 ¯ . And because the charges are distributed on each ring, the latter procedure for a chosen separation distance between the rings requires double integration over azimuthal angles of each ring. Assuming analytically this goal is achievable the output of the integration being a function of separation distance is to be differentiated with respect to the separation distance result the interactive force between the rings.

The distance p 1 p 2 ¯ is

p 1 p 2 ¯ = R 2 + r 2 + z 2 − 2 R r cos [ α − β ] , (1)

here, R, and r are the radii of the lower and upper circles: α, and β are their azimuthal angles, respectively. The physical quantity of interest is the electrostatic potential energy,

U ( z ) = 1 ( 2 π ) 2 k q 1 q 2 ∫ 0 2 π ∫ 0 2 π 1 p 1 p 2 ¯ d α d β . (2)

here, q_{1} and q_{2} are the charges of the lower and the upper rings, respectively. The k is the electrostatic coupling constant, k = 1 4 π ε 0 in metric units it is 8.99 × 10^{9} Nm^{2}/C^{2}.

For a chosen separation height z, integrating (2) applying Mathematica [

1) Select values for radii of the hoops and other relevant parameters.

2) Select a range for height z, 0 ≤ z ≤ ∞ .

3) For the list of chosen z’s subject to step 2 integrate (2). The output is a list of paired {z, integrated (2)}.

4) Plot the paired list of step 3 and fit the data with a physically meaningful z-dependent analytic function.

5) Utilize the fitted function in step 5, apply the fundamental relationship, F → = − k ^ [ d d z U ( z ) ] conducive the interactive force.

6) Apply result of step 5, set up the equation of motion.

7) Solve the equation of motion. Utilize the solution plot the z(t).

Foreseeing the ultimate goal is to form the equation of motion of the upper ring, among the relevant physical parameters the values includes the mass of the floating ring m. The radii of the rings, their respective charges, and the constants such as k and g are included as well; units are metric.

values = { k → 9 × 10 9 , q 1 → 2 × 10 − 6 , q 2 → 2 × 10 − 6 , R → 20 × 10 − 2 , r → 10 × 10 − 2 , m → 1 × 10 − 3 , g → 9.8 } ;

Compatible with the radii of the rings we consider z ∞ → 10 . According to step 3, the numeric double integration of (2) yields,

table Integrations =

It is assuring to note in the limit of r → 0 i.e. by shrinking the small ring into a point-like charge the integrand and consequently the integration converges to the well-known result of potential energy of a charged ring and a point-like charge.

The look of

model z [ z _ ] = a + b z + c e − d z ;

fit z = FindFit [ table Integrations , model z [ z ] , { a , b , c , d } , z ] ;

plotfit z = Plot [ model z [ z ] / . fit z , { z , 0 , 10 } , PlotStyle → Red , PlotRange → All ] ;

Its plot along with the fitted paired data is displayed in

Show [ list Plot table Integrations , plotfit z ]

As shown in

Force z = − D [ model z [ z ] / . fit z , { z , 1 } ] ;

Plot [ Force z , { z , 0 , 5 } , AxesLabel → { “ z ( m ) ” , “ F ( N ) ” } , P lotStyle → Black , GridLines → Automatic ]

As displayed, for small separation distances the force is strong and as the separation distance increases the force is weaker. It also shows the force is not a linear function of separation distance. Hence, nonlinear oscillations should be expected. More on this in forthcoming paragraphs.

Here we search for the separation equilibrium distance. Knowing this distance assists to determine how far up (down) the upper ring needs to be dropped (dropped) results oscillations about the equilibrium. To this end we solve the static equation setting the net force F_{electrostatic} − mg = 0.

Solve [ Force z = = m g / . values , z ] { z → 1 .64658 }

And then according to step 6 we form the equation of motion. Noting the force according to

Equation Of Motion = z ' ' [ t ] − ( Force z / . z → z [ t ] ) ( 1 m / . values ) + g / . values ;

solEquation Of Motion = NDSolve [ { Equation Of Motion = = 0. , z [ 0 ] = = 4. , z ' [ 0 ] = = 0 } , z [ t ] , { t , 0 , 10 } ] ;

We confirm our physical intuition, i.e. the upper ring dropped from a reasonable height above the equilibrium results oscillations.

plot ring Oscillations = Plot [ Evaluate [ z [ t ] / . solEquation Of Motion ] , { t , 0 , 10 } , AxesLabel → { “ t ( s ) ” , “ z ( m ) ” } , PlotStyle → Black , GridLines → Automatic , PlotRange → { 0 , 5 } ]

We also craft a code simulating the oscillating ring. Simulation of the oscillations along with the display of the z(t) helps to form a visual understanding of the problem at hand, this is shown in screen shot of the situation in

As shown in

Critically looking over and instructing off the material published in most of the undergraduate and graduate popular physics as well as mathematical physics texts reveal the incompleteness of most of the presented topics. Although online resources are complementary nonetheless the topic of the interest discussed in this exploratory investigation has not been reported in scientific literature. As described in the abstract and the body of the article technically speaking evaluation of the electric field created by a charged loop in the region of space off the symmetry axis brings about mathematical challenges. Interacting charge distributions on two different spatially separated entities calls for a double integration. Despite the powerful computational algorithms built in Mathematica and Maple their symbolic integrated double integration is complicated and almost useless. In this article we have devised a technical approach overcoming this issue. Accordingly, we have been successful evaluating the double integration numerically converting the output to an analytic continuous useful function. Interested readers may find our newly proposed approach applicable and useful for similar scenarios. Mathematica codes are imbedded in the article so that users familiar with this software [

The author gracefully acknowledges the John T. and Page S. Smith Professorship funds for completing and publishing this work.

The author declares no conflicts of interest regarding the publication of this paper.

Sarafian, H. (2020) Nonlinear Oscillations of a Pair of Charged Rings. American Journal of Computational Mathematics, 10, 571-577. https://doi.org/10.4236/ajcm.2020.104032