## The Dynamic Behavior of the Electrically Charged Cloud of the Ice Crystals

Department of Structural Mechanic, Technical University of Lodz, Lodz, Poland### Abstract

The paper includes the derivation of the equation of the two-dimensional, dynamic behavior of electrically charged cloud of ice crystals. A large crystal rotation angles and a continuous distribution of charges on the surface of the crystals are included in deliberations. Finally, possible solutions of model equation are discussed and compared with solutions available in the literature. The resulting model can be used as a mechanical basis for optic models of the atmospheric phenomenon called the “miracle of the sun”.

### At a glance: Figures

**Keywords:** the miracle of the Sun, finite difference method, non-linear vibrations, mathematical modeling

*Applied Mathematics and Physics*, 2014 2 (1),
pp 19-26.

DOI: 10.12691/amp-2-1-6

Received November 30, 2013; Revised January 18, 2014; Accepted February 07, 2014

**Copyright**© 2013 Science and Education Publishing. All Rights Reserved.

### Cite this article:

- Wirowski, Artur. "The Dynamic Behavior of the Electrically Charged Cloud of the Ice Crystals."
*Applied Mathematics and Physics*2.1 (2014): 19-26.

- Wirowski, A. (2014). The Dynamic Behavior of the Electrically Charged Cloud of the Ice Crystals.
*Applied Mathematics and Physics*,*2*(1), 19-26.

- Wirowski, Artur. "The Dynamic Behavior of the Electrically Charged Cloud of the Ice Crystals."
*Applied Mathematics and Physics*2, no. 1 (2014): 19-26.

Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |

### 1. Introduction

**1.1. The Motivation to Take the Research Topic**

The topic of the vibrations of the electrically charged cloud of the ice crystals, have been already the subject of the author’s research in the work ^{[1]}. The non-linear one-dimensional model of rotational vibration has been presented in that paper. It has been proposed as the mechanical basis of dynamic, atmospheric optical phenomenon commonly called the “miracle of the sun”. The first simply model of this phenomenon was presented in ^{[2]}.

However, one-dimensional model is not sufficient to a comprehensive consideration of phenomenon called the “miracle of the sun”, because it allows for only the analysis of the reflection of light from a single crystal of ice. The real cloud which is a three-dimensional body, is causing multiple reflection of sunlight. Hence the need to define a two-dimensional mechanical model, which in the future could be used to better and more precise modeling of the phenomenon called the “miracle of the sun”, and thus for better understand it.

**1.2. Aim and the Methodology of the Paper**

The aim of the paper is derivation of the differential equation which describes dynamic behavior of the cloud of electrically charged ice crystals. The considerations are limited to plate crystals, that assumes that the the horizontal dimension of the crystals is significantly larger than their thickness. A large crystal rotation angles and a continuous distribution of charges on the surface of the crystals are included in deliberations.

The methodology of the derivations of the model equation is analogous to that used in ^{[1]}. However, unlike to that considerations, four adjacent crystals distributed in two-dimensional space is taken into account. The starting point for the considerations are the electrostatic forces acting on a single crystal from adjacent. Next we calculate the rotating moments acting on the crystal. The resulting very complex equation is simplify to a short form that is usefull for further considerations. Ultimately by using the finite difference method formulas, we transform equation to the continuous form.

### 2. The Modelling

**2.1. The Assumptions**

Let us introduce the following modeling assumptions (Figure 1):

**Fig**

**ure**

**1.**The schematic representation of the modeling assumptions. The individual regular hexagonal plate crystals are shown in enlargement

- A homogeneous cloud of identical, electrically charged ice crystals is given,

- The crystals are evenly distributed in the cloud both vertically and horizontally,

- The crystals have the shape of regular hexagonal plates, i.e. their height is negligibly small in comparison to other dimensions,

- The gravitational force acting on every crystal is balanced by the force of air resistance acting on a crystal falling down uniformly,

- The air resistance forces influencing the crystal rotation are linearly proportional to the linear velocity in the rotational motion,

- The electric charge of the crystal is uniformly distributed on its surface,

- We assume that the synchronization occurs between the crystals, i.e. the difference between the rotation of adjacent crystals is negligible.

**2.2. The Equation of the Rotational Motion of an Individual Crystal**

(1) |

where:

- the moment of inertia for the axis of crystal’s rotation,

- the total resultant moment acting on the crystal,

- the angular acceleration in the rotational motion,

- the density of ice.

Due to the continuous charge distribution on the crystals’ surfaces, in the following modeling we consider the moments derived from the interaction between the appropriate surfaces of adjacent crystals. Then we can expand the resultant moment as (Figure 2):

(2) |

where: is a conductive vector of an infinitely small segment of the surface of the analyzed crystal, are the forces of the electrostatic interaction between the infinitely small segments of the surface of adjacent crystals respectively and , is the force of air resistance acting on an infinitely small segment of the crystal’s surface. The coordinates , are local coordinates, running respectively on the surfaces of the center and adjacent crystals and adopting the values from -½ a to ½ a for the extreme edges of the crystal and 0 in the midst of each surface (Figure 2). is the external forcing vibrations moment.

**Fig**

**ure**

**2.**The electrostatic forces between the infinitely small segments of the surfaces of the analyzed crystals

The resistance force, in accordance to the assumptions, is proportional to the linear velocity of segment , forces are Coulomb forces and therefore:

(3) |

where:

R is the coefficient of air resistance,

are vectors of beginnings and ends at the points P_{} and O respectively (Figure 2), which are the points of the analyzed segments of crystals’ surfaces and .

are the differences between the angles of crystals’ rotations respectively for *****th* and central crystal.

is the charge of the analyzed segment of the crystal’s surface,

is the relative dielectric constant of air,

is the dielectric constant of vacuum.

By substituting these relations into equation (2) and expanding the vector product we obtain:

(4) |

The equation (4) may be subject to certain simplifications. Let us note, that due to the uniform charge distribution on the surface of each crystal, we have:

(5) |

Due to the assumption of the synchronization of the crystals, which means that the adjacent crystals are rotated by almost the same angles, we have:

(6) |

Due to these relations, the equation (4) can be simplified after some transformations to the form:

(7) |

**2.3. The Derivation of the Discrete Equation of the Rotational Motion of the Crystal**

Although the solution of double integrals (7) exists in the general case, it is very complicated. In view of the fact that it is also a discrete equation of the rotation angle, which subsequently will be subject of conversion into a continuous form, the use of the exact solutions of the equation (7) for this purpose is impossible. To solve the given problem we can find the approximate solution of double integrals (7) with a reasonable accuracy. We use the assumption of the synchronization of crystals, i.e. the fact that the differences between the rotation angles of adjacent crystals are relatively small. In this case, we may assume that the moment induced by the crystals’ rotation by identical angles is independent of the moment induced by the components of the differences of the appropriate angles. Small differences between the angles of rotation crystals 1-4 in two directions are divided on the symmetric and antisymmetric parts. The result are a 5 component of the moment acting on the central crystal: principal moment of rotation deriving from identical angles of rotation and four components (symmetric-symmetric, antisymmetric-symmetric, symmetric-antisymmetric and antisymmetric-antisymmetric) deriving from different of angles of rotation of crystals 1-4. (Figure 3):

**Fig**

**ure**

**3.**The assumption of superposition of the moments caused by crystals’ rotation by identical angles and moments caused by the symmetric and antisymmetric components of the differences of the appropriate angles, which is described by (9)

Let us denote the integrals from equation (7) by :

(8) |

We assume, that the following superposition of the moments is true:

(9) |

Let us calculate the moment . In this case we have:

Both integrals (7) give the same solution for the couple of the crystals andand the second couple and due to the symmetry of the task (Figure 3), we get:

(10) |

It can be proved that:

(11) |

due to the symmetry of the moments induced by the electrostatic forces from corresponding pairs of crystals, as it is shown in Figure 3. For the remaining we have:

After the transformations and the substitutions, which are analogous to the ones performed in the calculation of *K*_{0}_{,} we obtain:

(12) |

We can calculate the integrals (10) and (12) in an accurate manner. By denoting respectively (in symbols the upper operator corresponds to , the lower one to )

we get:

(13) |

and by denoting:

We get (in (14) summation convention is used):

(14) |

Finally, the equation of the crystals’ rotation takes the form:

(15) |

**2.4. The Approximation of the Discrete Equation of the Motion**

The obtained equation (15) is a very complicated formula of the parameters φ_{0} and Δ Let us notice that in order to be able to transform it to a continuous model using the finite difference method formulas, the coefficients *K*_{0} and *K*_{AA} must depend on the parameters φ_{0}_{ }and Δ in the simplest possible way.

Hence, the series of graphs illustrating the correlation between the coefficients *K*_{0} and *K*_{AA} and the parameters φ_{0}_{ }and Δ have been performed. The Figure 4 and Figure 5 have been made for the selected allowable parameters of the model:

In the case of *K*_{0} we get good results in the range of angles after using a polynomial approximation of the 3rd order. It is worth noting that the polynomial approximating *K*_{0} has the form (Figure 4):

(16) |

The remaining coefficients of the 3-rd order development are equal nearly to zero due to the specific shape of the graph *K*_{0} and they can be ignored, regardless of the combination of the other parameters of the model.

**Fig**

**ure**

**4.**The dependency of the coefficients

*K*

_{0}from the discrete parameter

_{0}for the selected combinations of the parameters

*a*and

We apply the linear approximation of the coefficient *K*_{AA}, which is the exact approximation in the range of angles. .

(17) |

**Fig**

**ure**

**5**. The dependency of the coefficients

*K*

_{AA}from the discrete parameter for the selected combinations of the parameters

*a*and

We can also note that the approximation coefficients ,and are positive, regardless of the parameters of the model, and the following inequality is satisfied:

(18) |

As a result of the polynomial approximation of the coefficients *K*_{0} and *K*_{AA }the following equation is obtained:

(19) |

where

The equation (19) has coefficients, which are independent of the discrete parameters φ_{0} and Δ and it can be used to construct the continuous model of the vibrations of the cloud of the ice crystals.

**2.5. The Continuous Model of the Vibrating Cloud**

Next we transform the discrete equation (19) to a continuous form in order to build the continuous model of the vibrating cloud. Let us rotate the coordinate system by an angle of ∏/4 to obtain a regular finite difference method grid, which is shown in Figure 6:

**Fig**

**ure**

**6**. The adopted finite difference method grid with the rotated coordinate system x’ and y’

We will use the finite difference method formulas ^{[12]}:

(20) |

Let us note, that (8) (Figure 3). By substituting the relations (20) into the equation (19) and after appropriate transforming and grouping the components, we get:

(21) |

where

Then we go back with the turnover of the coordinate system by an angle of ∏/4 and we transform the coordinate system according to patterns:

(22) |

The corresponding derivatives are as follows:

(23) |

(24) |

By substituting (23) and (24) into (21) we obtain the final equation in x and y coordinates:

(25) |

The equation (25) describes the non-linear vibrations of the cloud of the electrically charged ice crystals in the continuous form in two dimentions. It takes into account the continuous distribution of the electric charge on the surface of the crystals, their rotation by large angles, the external forcing and the air resistance.

**2.6. The Discussion of the Possible Solutions**

The equation (25) is the second degree differential equation with respect to time and space. It is a generalization in the two-dimensional space equation called in the literature as the equation φ^{4} ^{[3, 4]} presented in ^{[1]}. It has analytical solutions only for some very specific boundary and initial conditions, and specific parameters of the equation. General analytical solutions of the equation (25) does not exist, so it should be solve numerically.

However, some solutions of special cases of simplified problems can suggest us, how they might look like potential solutions of equation (25). In ^{[1]} it has been shown the selected solution of this equation with the simplification of assuming a total lack of wave propagation in the spatial directions.

Then the equation (25) is converted into a simpler form which is called Duffing equation ^{[5]}. It has been shown, that the oscillations described by this equation may be of a harmonic, quasi-harmonic or chaotic, in dependency on the parameters of the equation and the initial-boundary conditions.

In the case of the equation (25) it can be expected that the disturbance of any kind, will be transmitted in the directions x and y, depending on the initial-boundary conditions, external torque which inducement vibration, damping parameters, and the model parameters corresponding to the propagation of the vibration between the crystals. We will have to deal in some areas with chaotic vibrations in other with harmonics, and these areas will be vary in time and space. The computer simulations will allow us for the analysis and visualization of the solutions of equation (25). It will be the subject of further research of the author.

**2.7. The Comparison with Data from Motivation of Paper**

As previously mentioned, in the general case, the equation (25) can be solved only numerically. However, in certain very specific cases, we can find analytical solutions of this equation and compare them with the solutions derived from linear and one-dimensional model proposed in ^{[2]} and use them to model the phenomenon commonly called “miracle of the sun”. Below we present a very simple example of the use of equation (25) for modeling this phenomenon.

Suppose we are given the cloud of ice crystals satisfying the assumptions of the model with dimensions *L*_{x} to* **L*_{y}. For simplicity, let’s assume that we are dealing only with the issue of sinusoidal forced and linear vibration without damped, ie:

(26) |

The assumption of linearity of vibration is equivalent to the assumption of small deflection angles, smaller than 20-30°. (see Figure 4). Then the equation (25) simplifies greatly to the form:

(27) |

Suppose we are given the following boundary conditions:

(28) |

and the initial conditions:

(29) |

We get the following solution of equation (27) with boundary conditions (28) and initial conditions (29):

(30) |

Due to the nature of forcing vibration and the boundary-initial conditions, the solution does not depend on the model parameters, but they affect only on the amplitude of vibration. Thus, we can simplify solution (30) to the form:

(31) |

where *A*φ_{}_{}**is a constant amplitude dependent from parameters of the model and amplitude of forced vibration.

When we have a simple form of field of rotation of crystals (31), we can take advantage of the optical conditions given in ^{[1]}, which defining the relative brightening or darkening of the sun. Figure 7 shows an example resulting obtain from formula (31) in the spatial coordinates, that is, the distribution of crystals in the analyzed cloud. In addition, to illustrate the idea of the phenomenon “miracle of the sun”, the situation of four observers is considered. The sun’s rays fall from the right side on the vibrating cloud crystals. At time *t* shown in Figure 7, as a result of reflections from the ice crystals, a part of the rays is dispersed, a part is concentrated and others rays passed without encountering any a crystal. Therefore, depending on the position of the observer in relative to the crystal cloud, the visual effect will be very different. At the time *t* shown in the figure:

- the observer A will see a darkening of the image of the sun, because the cloud will disperse part of the sunlight that would reach him,

- the observer B is observed brightening the image of the sun, because the cloud will focus part of the rays,

- the observer C will see additional area of brightening in the sky above the real sun,

- the observer D will not notice any change in the appearance of the sun.

**Fig**

**ure**

**7**

**.**An example of the use of solutions of equation (31) for the analysis of the phenomenon called “miracle of the sun”. For the simulation we used 66 crystals and 48 sun rays. The crystal system corresponds to the moment

*and amplitude 20°. The figure shows the situation of the four observers and their observations. By using the two dimensional model, we can explained the fact that their observations may be extremely different*

As the solution obtain from the formula (31) is dependent on time, the situation which is shown on Figure 7 is dynamically changing. A similar dynamics, time and spatial diversity of observation was characteristic of the witnesses descriptions of the phenomenon called miracle of the sun in 1917 ^{[6]}. The previous attempts of the scientific explanation of the phenomenon of the “miracle of the sun” ^{[7, 8]} did not explain such a large spatial variability at all.

As we can see, from the two-dimensional solutions of equation (25) we can provide much greater opportunities for modeling the phenomenon of the “miracle of the sun” than the one-dimensional solution shown in ^{[2]}.

**2.8. The Other Applications of the Model**

The proposed model can be also used to modeling other physical phenomena. From the mechanics point of view is not an essential if we are dealing crystals of ice or other tiles made from different material. Just as the forces of interaction between them may not necessarily be associated with electrostatic forces. We can imagine that the role of the interaction forces act for example magnetic forces.

This model could be used for example to model the mechanics of the Kerr phenomena ^{[9]}, or Cotton-Mouton phenomena ^{[10]}. Similarly, it could be used to model the behavior of some of the liquid crystals in the dynamically changing electric field ^{[11]}.

### 3. The Conclusions

In this paper, the differential equation which describe the mechanical vibrations of the electrically charged cloud of plated ice crystals was obtained. It determines a mathematical description of the dynamic behavior of this media. The derived equation is a nonlinear, the second degree differential equation with respect to time and space. Although its general analytical solutions do not exist, it could be used to find numerical solutions of the problem of dynamic behavior of electrically charged cloud of ice crystals. That solutions could be used for further modeling of atmospheric phenomena such as optical phenomenon called “miracle of the sun”.

### References

[1] | A. Wirowski, The Non-linear Modeling of the Rotational Vibrations of the Electrically Charged Cloud of the Ice Crystals, Open Journal of Mathematical Modeling, Vol. 1, No. 2 (2013), p. 46-57. | ||

In article | CrossRef | ||

[2] | A. Wirowski, Modelling of the phenomenon known as “the miracle of the Sun” as the reflection of light from ice crystals oscillating synchronously, Journal of Modern Physics, Vol. 3 No. 3, 2012, pp. 282-289. | ||

In article | CrossRef | ||

[3] | A.A. Soliman, H.A. Abdo, New Exact Solutions of Non-linear Variants of the RLW, the PHI-four and Boussinesq Equations Based on Modified Extended Direct Algebraic Method, International Journal of Non-linear Science, Vol. 7 (2009), No. 3. | ||

In article | |||

[4] | A. Bekir, New Exact Travelling Wave Solutions for Regularized Long-wave, Phi-Four and Drinfeld-Sokolov Equations, International Journal of Non-linear Science, Vol. 6 (2008), No. 1. | ||

In article | |||

[5] | I. Kovacic, M.J. Brennan, The Duffing Equation, Non-linear Oscillators and their Behaviour, John Wiley & Sons, 2011. | ||

In article | CrossRef | ||

[6] | J. De Marchi, The Immaculate Heart, Farrar, Straus and Young, New York, 1952. | ||

In article | |||

[7] | S. Campbell, The Miracle of the Sun at Fatima, Journal of Meteorology, Vol. 14, Np. 142, 1989. | ||

In article | |||

[8] | P. Simons, Weather Secrets of Miracle at Fatima, The Times, 17 February 2005. | ||

In article | |||

[9] | M. Faraday, J. Kerr, P. Zeeman, The Effects of a Magnetic Field on Radiation, American Book Company, 1900. | ||

In article | |||

[10] | P. Daveze, B. Payet, L. Delaunay, Using dynamic Cotton and Mouton effect to study Brownian relaxation and magnetization of ferrofluids, Proc. SPIE 3098, Optical Inspection and Micromeasurements II, 456 (July 4, 1997). | ||

In article | CrossRef | ||

[11] | R. Vepa, Dynamics of Smart Structures, Wiley 2010. | ||

In article | CrossRef | ||

[12] | J.W.Thomas, Numerical Partial Differental Equations-Finite Difference Methods, Vol. 1, Springer, 1995. | ||

In article | CrossRef | ||