Figure 1 depicts the problem at hand. As shown, the blue marked circle is the boundary of the spherical charge distribution. It contains a net charge, +Q which is evenly distributed within the sphere of radius R, e.g. R = 5 units. Within this sphere there is a void colored pink of radius r; r < R. The bubble may move within the sphere on an arbitrary path. Point A is an exterior point of interest. We wish to calculate the electric field, E at A. Since at point A each sphere acts as a point charge with the concentrated charge at their respective centers, first we get the field of the large sphere, and then we treat the void as a uniform charge distribution with a net negative, −Q', charge. Since both spheres contain charges with even distributions this yields the charge within the bubble to, Q' =

Q(r/R)³. The blue arrow connects the origin to A, and r 1 shown with the black arrow connects the origin to the center of the void, and the pink arrow connects the center of the void to A, i.e. Δ r = r A − r 1 . Hence the electric field at A becomes [2] [3],

E A = k Q r A r A 3 + k Q ′ Δ r | Δ r | 3 , (1)

As shown, the E A contains two vector fields. The field from the large sphere, R is along the extension of the blue arrow, while the field from the pink is along the opposite direction of the pink arrow. A combination of these two vectors is given by (1), where k = 1 / ( 4 π ϵ 0 ) = 9 × 10 9 , SI units. Breaking the vectors of interest in horizontal and vertical components gives,

r A = r A cos [ θ ] i ^ + r A sin [ θ ] j ^ , (2)

Δ r = ( r A cos [ θ ] i ^ + r A sin [ θ ] j ^ ) − ( r 1 cos [ ϕ ] i ^ + r 1 sin [ ϕ ] j ^ ) , (3)

where θ and ϕ, are the angles between the r A and r 1 with the horizontal axis. This gives,

| Δ r | 3 = ( r A 2 + r 1 2 − 2 r A r 1 cos [ θ − ϕ ] ) 3 2 , (4)

Substituting (4) in (1) and simplifications yield the components of E A ,

E A x k Q = 1 r A 2 cos [ θ ] − ( r R ) 3 1 | Δ r | 3 ( r A cos [ θ ] − r 1 cos [ ϕ ] ) , (5)

E A y k Q = 1 r A 2 sin [ θ ] − ( r R ) 3 1 | Δ r | 3 ( r A sin [ θ ] − r 1 sin [ ϕ ] ) , (6)

The magnitude of E A and its angular orientation are, respectively

E A = E A x 2 + E A y 2 , (7)

Θ = arctan [ E A y E A x ] , (8)

And since one of our objectives is to analyze the impact of the mobile void, and its relevant impact on the E , we set the value of kQ = 1. Our study falls into one of three categories.

Case 1. The bubble has the size half of the size of the radius of charge distribution, e.g. r = 1/2R. The simplest scenario is when point A is along the extension of the radial vector and the center of the void and the center of the spherical charge are aligned. This is shown in Figure 2.

In this situation as shown in Figure 2 angles, θ and ϕ overlap. The E_A becomes,

E A = 1 r A 2 − 1 8 1 ( r A − R 2 ) 2 , (9)

Here we consider a static condition where the void is stationary. The field at A is time-independent and is given by (9). The extension of this case is a scenario where the void is smaller than the case discussed. The formulation is straight-forward, the situation is depicted in Figure 3.

Case 2. This case deals with either a displaced A that comes about by moving the observation point A around the initial position and/or preferably by circulating the void around the center of the charge distribution, as shown in Figure 4. Either one of these scenarios has an implication, the said vectors become misaligned making θ ≠ ϕ . Here, (5) and (6) become active.

In Figure 4, we show a situation where r < 1/2R. Of course the size of the bubble may be adjusted at will. For the given situation and chosen angles θ and ϕ, (5) and (6) give the components of the E_A. Instead of giving a set of numeric values for field components because in the forthcoming case 3, here we would rather postpone the numeric calculation.

Case 3a. The impact of the mobile bubble

Here we consider a situation where the void is mobile. We have set the formulation of the problem such that any kind of path could be handled. For the sake of demonstration only here we consider a situation where the bubble circulates at a constant pace around the center of the charge distribution. To achieve this goal we replace the static orientation angle, θ with a time-dependent entity, i.e. θ → ω t . Where ω is the angular speed, ω = 2πft with the frequency f at our disposal. The trajectory of the orbital motion is shown in Figure 5.

The far-right figure is a snapshot of the animation. If the animation could have been activated one could see the circulating bubble around the center tracing the shown pink circle. We notice the powerful CAS’s usefulness. Without it, these would have been a long computational procedure.

Case 3b. The impact of the oscillating-circulating mobile bubble

As pointed the bubble may move on an arbitrary path. To make the situation just a tad challenging here we consider a situation where the bubble with its constant radius not only circulates about the center, O but simultaneously oscillates along the position vector r₁. This is implemented by replacing the constant ϕ with ϕ → 2 π   f ϕ t . Here f_ϕ is the frequency of the oscillating void along its radial direction. We noticed this situation somewhat resembles the Lissajous [4] phenomenon. Here we demonstrate two different examples. First we set the {θ, ϕ} = (2π 1 t, 2π 3 t} and {θ, ϕ} = {2π 1 t, 2π 5 t}, respectively. For the former case, the ratio θ/ϕ = 1/3 and for the latter θ/ϕ = 1/5. These ratios control the number of peddles in the carnation shown by the pink paths. Figure 6 is associated with the first and Figure 7 with the second ratio, respectively.

In Figure 6 and Figure 7, the far left drawings have overlapped curves, i.e. the blue-green appears green.

We begin this project thinking about how to handle the impact of a void in a chosen charge distribution. In the course of analyzing the problem not only we addressed the issue but ended up expanding the problem from a static situation to a dynamic one. Our analysis handles the impact of a mobile bubble within a charge distribution. A few cases are presented. At some point, the mobility of the void is related to that of the Lissajous curves. The presented analysis has the potential to be extended even more, e.g. multiple bubbles with arbitrary pathways. The kinematics of the mobile voids is totally under our control. Mathematica can achieve the set goals. Although the presented article focuses only on the static aspect of the proposed scenario, the author enjoyed watching the animations. The learned lesson is priceless. The interested reader may find [5] [6] resourceful for making the reported graphs, calculations, and animations.

The author 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. (2022) Impact of a Mobile Void within a Charge Distribution on the Electric Field. American Journal of Computational Mathematics, 12, 349-354. https://doi.org/10.4236/ajcm.2022.124025