Next Article in Journal
Eulerian–Eulerian RSTM-PDF Modeling of Turbulent Particulate Flow
Previous Article in Journal
From Ion Fluxes in Living Cells to Metabolic Power Considerations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mylar Balloon and Associated Geometro-Mechanical Moments

by
Vasyl Kovalchuk
1,*,
Vladimir I. Pulov
2 and
Ivaïlo M. Mladenov
3,4
1
Institute of Fundamental Technological Research, Polish Academy of Sciences, 5B, Pawińskiego Str., 02-106 Warsaw, Poland
2
Department of Mathematics and Physics, Technical University of Varna, Studentska Str. 1, 9010 Varna, Bulgaria
3
Institute for Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, Tsarigradsko Chaussee 72, 1784 Sofia, Bulgaria
4
Institute of Mechanics, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Bl. 4, 1113 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(12), 2646; https://doi.org/10.3390/math11122646
Submission received: 10 May 2023 / Revised: 5 June 2023 / Accepted: 8 June 2023 / Published: 9 June 2023

Abstract

:
Starting with identifications of the very fundamental geometric characteristics of a Mylar balloon such as the profile curve, height, volume, arclength, surface area, crimping factor, etc., using the geometrical moments I n ( x ) and I n , we present explicit formulas for them and those of the mechanical moments of both solid and hollow balloons of arbitrary order. This is achieved by relying on the recursive relationships among elliptic integrals and the final results are expressed via the fundamental mathematical constants such as π , lemniscate constant ω ˜ , and Gauss’s constant G. An interesting periodicity modulo 4 was detected and accounted for in the final formulas for the moments. The principal results are illustrated by two tables, a few graphics, and some direct relationships with other fundamental areas in mathematics, physics and geometry are pointed out.

1. Introduction

Before focusing on the actual contents of the paper the reader might be interested to learn something about the Mylar balloon itself. According to Webster’s New World Dictionary, Mylar is an Americanism meaning a trademark for a polyester made in extremely thin sheets of great tensile strength. This term has been introduced by Paulsen [1] in the form of a variational problem under a constraint and with an idea to replace the so-called e-balloon explored intensively in the 70s by NASA (cf., e.g., [2]). More detail about the present day situation with the zero-pressure balloons can be found in [3].
According to Paulsen the Mylar balloon is constructed by taking two identical circular disks of Mylar, sewing them together along their boundaries, and then inflating the resulting object with either air or helium. Surprisingly enough, these balloons turn out not to be spherical as one might expect based on the well-known fact that the sphere possesses the maximal volume for a given surface area. Respectively, this purely experimental fact suggests the following mathematical problem: given a circular Mylar balloon of deflated radius a, what will be the shape of the balloon when it is fully inflated (i.e., inflated to maximize the volume subject to the constraint that the Mylar does not stretch)?
It should be also noted that the Japanese engineer Kawaguchi [4] had arrived at the same figure when looking for a shallowest surface of revolution which has no circumferential stress. This is achieved exactly when the meridional k μ and parallel k π principal curvatures obey everywhere the condition
k μ = 2 k π .
Some years later Mladenov and Oprea [5] proved that this linear-type relationship between curvatures specifies the rotational surfaces uniquely.
The rest of the paper is organized as follows. In Section 2 the geometrical moments I n ( x ) and I n are defined and it is shown how the main characteristics of the Mylar balloon can be expressed via them. In Section 3 the mechanical moments J n and J n of the solid and hollow Mylar balloons are studied and it is shown that in both considered cases they can be connected to the geometrical moments I n . Section 4 is devoted to investigation of the recurrent relations (modulo 4) for the sequence of the geometrical moments I n ( x ) and I n , while in Section 5 the corresponding residue classes of the mechanical moments J n and J n for the solid and hollow Mylar balloons are shown and the first twenty values of the numerical coefficients μ n and μ n are calculated explicitly via Gauss’s constant G with the aid of the Wolfram computer language using the symbolic computation program Mathematica® [6]. Finally, Section 6 concludes the paper and points out a possible direction for the related future work. The main new results of the presented research can be found in Section 2, Section 3, Section 4 and Section 5. The paper deals with the fundamental problem which appears at least in Celestial Mechanics, Probability Theory and Statistics and could help significantly in other areas where one needs to account for corrections keeping some control over them.

2. Balloon’s Characteristics Expressed via Geometrical Moments

The Mylar balloon (see Figure 1) can be defined as a surface of revolution that is obtained when we maximize the volume V (no stretching of the Mylar foil is allowed) for a given (fixed) arclength L a (where a is the radius of the deflated Mylar balloon, i.e., of the Mylar discs) of its profile curve (i.e., the directrix that generates the surface) z = z ( x ) for x changing from x = 0 to x = r , where r is the radius of the inflated Mylar balloon, and obviously we have r < a (see, e.g., [1,7,8]).
It can be shown that the constraint equation for the profile curve can be written as
z ( x ) = x 2 r 4 x 4 , i . e . , z ( x ) = x r t 2 d t r 4 t 4 x [ 0 , r ] .
Therefore, the height, volume, and arclength of the Mylar balloon can be found, respectively, as
τ = 2 z ( 0 ) = 2 0 r t 2 d t r 4 t 4 , V = 2 π 0 r t 4 d t r 4 t 4 , L = r 2 0 r d t r 4 t 4 a .
Similarly, we can calculate the surface area of the inflated Mylar balloon by integrating over all the circular strips of the surface of revolution with the length 2 π x , where x is the distance to the axis of rotation, and the infinitesimal width d s , where d s 2 = d x 2 + d z 2 , i.e.,
S inflated = 2 π x d s = 2 π x d x 2 + d z 2 = 4 π x = 0 r x 1 + ( z ( x ) ) 2 d x
where we can substitute the expression for the derivative z ( x ) given in (2) and obtain
S inflated = 2 π r 2 0 r d x 2 r 4 x 4 = 2 π r 2 arcsin x 2 r 2 0 r = π 2 r 2 .
We can see that the above surface area clearly differs from the surface area of the two sewn-together circular Mylar discs of radius a, i.e., S deflated = 2 π a 2 , describing effective shrinking of the inflated Mylar balloon’s surface compared to the deflated one
S deflated S inflated = 2 π a r 2 1.09422 .
We can also introduce the local measure of the above-described shrinking that is defined as the ratio of the surface area of a small patch on the deflated Mylar balloon to the surface area of the corresponding patch on the inflated Mylar balloon. This local measure is called the crimping factor and is defined through the relation (see, e.g., [1])
C ( x ) = r 2 x 0 x d t r 4 t 4 ·
It can be shown that the crimping factor is minimal for x = 0 (see Figure 2–right), i.e.,
C min = C ( 0 ) = lim x 0 r 2 x 0 x d t r 4 t 4 = 0 0 L’ H o ^ spital’s rule = lim x 0 r 2 r 4 x 4 = 1
and maximal for x = r (again see Figure 2–right), i.e.,
C max = C ( r ) = r 0 r d t r 4 t 4 = L r a r 1.31103 .
Therefore, we can see that all the above relations clearly indicate that the Mylar balloon can be actually characterized by the set of geometrical moments that are defined as
I n ( x ) : = 0 x t n d t r 4 t 4 , I n : = I n ( r ) , I n ( 0 ) = 0 , n = 0 , 1 , 2 ,
Then we can write down the above formulas in the alternative forms, i.e., z ( x ) = I 2 I 2 ( x ) , τ = 2 I 2 , V = 2 π I 4 , L = r 2 I 0 a , S inflated = π 2 r 2 = 4 π r 2 I 1 , and C ( x ) = ( r 2 / x ) I 0 ( x ) .
Finally, the explicit parametrization ( u [ 1 , 1 ] ) of the inflated Mylar balloon’s profile curve (directrix) can be obtained from (2) as (see, e.g., [9])
x ( u ) = r 1 u 2 1 + u 2
z ( u ) = r 2 2 E arcsin 2 u 1 + u 2 , 1 2 F arcsin 2 u 1 + u 2 , 1 2
where F ( φ , k ) and E ( φ , k ) denote the incomplete elliptic integrals of the first and second kind with the Jacobian amplitude φ and elliptic modulus k (see, e.g., [10]). From (11) and (12) we can also obtain the explicit dependency z ( x ) that is given as (see Figure 1)
z ( x ) = r 2 2 E arccos x r , 1 2 F arccos x r , 1 2 .

3. Mechanical Moments of Solid and Hollow Mylar Balloons

It is well known that in order to determine the attraction between gravitating bodies one needs to know the inertia moments of arbitrary order. In most of the astronomical applications only a few of them are used. However, if the difference in masses is significant and the space in which the bodies are considered is confined, the higher moments are essential as well. This holds definitely in the case of the motion around the oblate primary [11]. Strangely enough, the higher moments appear in probability and statistics (cf. [12]) where at least the first six of them have specific names.
In our setting specifically, we suppose that either the volume or the surface area of the Mylar balloon is filled with the physical material of the total mass m that has the homogeneous distribution of the mass (i.e., the constant density). Thus, in the first case we will obtain the solid Mylar balloon (with the constant volumetric mass density ρ = m / V ), whereas in the second case it will be the hollow Mylar balloon (with the constant surface mass density ρ = m / S inflated ).
For such two (solid and hollow) physical realizations of the Mylar balloon we can define the corresponding mechanical moments of mass calculated with respect to the vertical axis of symmetry of the described surface of revolution as
J n = V x n d μ V = ρ V x n d V , J n = S inflated x n d μ S = ρ S inflated x n d S , n = 0 , 1 , 2 ,
where x is the distance from the vertical axis of symmetry to the infinitesimal mass d μ that is positioned either in the volume or on the surface of the Mylar balloon, μ V and μ S are the corresponding measures of the volumetric and surface mass distributions, and the circle indicates that the respective quantities are associated with the hollow Mylar balloon.
Hence, the zeroth (monopole) and second (quadrupole) mechanical moments of mass will correspond to the total mass (the same by construction in both cases) and two polar moments of inertia for the solid and hollow Mylar balloons respectively, i.e.,
J 0 = J 0 = m , J 2 = M = ρ V x 2 d V , J 2 = M = ρ S inflated x 2 d S .

3.1. Mechanical Moments of Solid Mylar Balloon

Using cylindrical coordinates ( x , z , θ ) , where θ [ 0 , 2 π ] , and taking into account the symmetry of the Mylar balloon with respect to the equatorial ( X O Y ) plane (see Figure 1) the mechanical moments J n given by the first expression in (14) can be rewritten as
J n = 2 ρ x = 0 r z = 0 z ( x ) θ = 0 2 π x n + 1 d θ d z d x = 4 π ρ x = 0 r x n + 1 z ( x ) d x = 4 π ρ x = 0 r t = x r x n + 1 t 2 r 4 t 4 d t d x
where we have used the integral representation (2) of the expression z ( x ) . After changing the order of integration (from ( x , t ) [ 0 , r ] × [ x , r ] to ( t , x ) [ 0 , r ] × [ 0 , t ] ) we obtain
J n = 4 π ρ t = 0 r x = 0 t t 2 x n + 1 r 4 t 4 d x d t = 4 π ρ n + 2 t = 0 r t n + 4 d t r 4 t 4 = 4 π ρ n + 2 I n + 4
i.e., the mechanical moments J n can be expressed through the geometrical moments I n + 4 given by (10). It also turns out to be convenient to represent J n in the form
J n = μ n m r n , m = ρ V = 2 π ρ I 4 = J 0 , μ n = 2 n + 2 · I n + 4 I 4 · 1 r n
where μ n are the numerical coefficients that will be explicitly evaluated in Section 5.

3.2. Mechanical Moments of Hollow Mylar Balloon

Similarly to the solid case from the previous section, for the hollow Mylar balloon we can rewrite the mechanical moments J n given by the second expression in (14) as
J n = 2 ρ x = 0 r θ = 0 2 π x n + 1 1 + ( z ( x ) ) 2 d x d θ = 4 π r 2 ρ x = 0 r x n + 1 d x r 4 x 4 = 4 π r 2 ρ I n + 1
where we have used the expression for the derivative z ( x ) given in (2).
As we can see, for the hollow Mylar balloon the mechanical moments J n can be expressed through the geometrical moments I n + 1 and it is convenient to represent them as
J n = μ n m r n , m = ρ S inflated = π 2 r 2 ρ = 4 π r 2 ρ I 1 = J 0 , μ n = 4 π I n + 1 r n = I n + 1 I 1 · 1 r n
where the numerical coefficients μ n will be again explicitly evaluated in Section 5.

4. Recursive Evaluation of Geometrical Moments I n ( x ) and I n

In order to obtain the recursive formula for the geometrical moments I n ( x ) given by (10) we can use the expression (cf., e.g., [13] Formula (250.01)).
I n + 4 ( x ) = n + 1 n + 3 I n ( x ) r 4 x n + 1 n + 3 r 4 x 4 , n = 0 , 1 , 2 ,
or we can reobtain it directly using integration by parts and the differential relation
t 3 d t r 4 t 4 = 1 2 d r 4 t 4 .
Then the first integral expression in (10) taken for n + 4 can be rewritten in the form
I n + 4 ( x ) = 1 2 0 x t n + 1 d r 4 t 4 = t n + 1 2 r 4 t 4 0 x + n + 1 2 0 x t n r 4 t 4 d t
where after interpreting the last integral as r 4 I n ( x ) I n + 4 ( x ) we reobtain (21) for n 0 .
Hence, using the explicit representations of the first four geometrical moments I n ( x ) (for the dependencies, I n I n ( x ) , see Figure 2–left; for the values of I n , see (30)–(38)), i.e.,
I 0 ( x ) = 1 2 K 1 2 F arccos x r , 1 2 · 1 r
I 1 ( x ) = 1 2 arcsin x 2 r 2
I 2 ( x ) = 1 2 2 E 1 2 K 1 2 2 E arccos x r , 1 2 + F arccos x r , 1 2 · r
I 3 ( x ) = 1 2 1 1 x 4 r 4 · r 2
and the expression (21) we can obtain all the geometrical moments I n ( x ) for n 0 , e.g., the next two are given as
I 4 ( x ) = 1 3 1 2 K 1 2 F arccos x r , 1 2 x r 1 x 4 r 4 · r 3
I 5 ( x ) = 1 4 arcsin x 2 r 2 x 2 r 2 1 x 4 r 4 · r 4
where F ( φ , k ) and E ( φ , k ) denote the incomplete elliptic integrals of the first and second kind with the Jacobian amplitude φ and elliptic modulus k, whereas K ( k ) = F ( π / 2 , k ) and E ( k ) = E ( π / 2 , k ) are the complete elliptic integrals of the first and second kind (see [10]).
The corresponding values of the first four geometrical moments I n are given as
I 0 = 1 2 K 1 2 · 1 r 1.31101 · 1 r
I 1 = π 4 0.785398
I 2 = 1 2 2 E 1 2 K 1 2 · r 0.59907 · r
I 3 = 1 2 · r 2 = 0.5 · r 2
whereas the rest of them for all values n 0 can be calculated using the recursive relation
I n + 4 = n + 1 n + 3 I n r 4
e.g., the next four geometrical moments are given as
I 4 = 1 3 I 0 r 4 = 1 3 2 K 1 2 · r 3 0.43701 · r 3
I 5 = 1 2 I 1 r 4 = π 8 · r 4 0.392699 · r 4
I 6 = 3 5 I 2 r 4 = 3 5 2 2 E 1 2 K 1 2 · r 5 0.359442 · r 5
I 7 = 2 3 I 3 r 4 = 1 3 · r 6 0.333333 · r 6 .
By applying (34) iteratively for n p ( mod 4 ) , p = { 0 , 1 , 2 , 3 } , i.e., for each of the four residue classes n | p , it can be readily obtained that
I n + 4 | p = ( n + 1 ) ! ! ! ! ( n + 3 ) ! ! ! ! I p r n + 4 p
where the shortcut notation n ! ! ! ! stands for the quadruple factorial, i.e., a product of positive integers defined by the formula n ! ! ! ! = n ( n 4 ) ( n 8 ) ( n 12 ) It can be easily seen that, when n is an even number, i.e., when n = 2 k with k 0 , then n ! ! ! ! = 2 k k ! ! , where k ! ! stands for the double factorial, i.e., a product of positive integers defined by the formula k ! ! = k ( k 2 ) ( k 4 ) ( k 6 )
Finally, let us present some exemplary calculations of the characteristics of the Mylar balloon expressed through the geometrical moments I n ( x ) presented in Section 2, i.e., the profile curve z ( x ) and the crimping factor C ( x ) are given as (see Figure 2–left for z ( x ) = I 2 I 2 ( x ) , right for C ( x ) )
z ( x ) = I 2 I 2 ( x ) = r 2 2 E arccos x r , 1 2 F arccos x r , 1 2
C ( x ) = I 0 ( x ) x · r 2 = r x 2 K 1 2 F arccos x r , 1 2
(the first expression is identical to the parametrization (13) of the Mylar balloon), whereas the arclength L, equivalently, deflated radius a, height τ , and volume V are calculated as
L a = I 0 r 2 = 1 2 K 1 2 · r 1.31101 · r
τ = 2 I 2 = 2 2 E 1 2 K 1 2 · r 1.19814 · r
V = 2 π I 4 = 2 π 3 I 0 r 4 = 2 π 3 K 1 2 · r 3 2.74581 · r 3 .

5. Residue Classes Modulo 4 of Mechanical Moments

By making use of the recursive Formula (39) for the geometrical moments I n for n 0 , we can arrive at an explicit representation of the mechanical moments J n given by (18) and J n given by (20) classifying them into four non-intersecting residue classes J n | p as
J n | p = μ n | p m r n , J n | p = μ n | p m r n , n p ( mod 4 ) , p = { 0 , 1 , 2 , 3 } .
In the next two sections we will discuss in detail the recursive calculation of the numerical coefficients μ n | p and μ n | p for the cases of the solid and hollow Mylar balloons. Most of the calculations were made by the aid of the Wolfram computer language using the Computer Algebra System Mathematica® [6].

5.1. Calculation of Numerical Coefficients μ n for Solid Mylar Balloon

For the solid Mylar balloon the numerical coefficients μ n | p given by (18) using the recursive Formula (39) can be calculated as
μ n | p = 2 n + 2 · I n + 4 | p I 4 · 1 r n = 6 n + 2 · I n + 4 | p I 0 · 1 r n + 4 = 6 n + 2 · ( n + 1 ) ! ! ! ! ( n + 3 ) ! ! ! ! · I p I 0 · 1 r p
or substituting directly the values of I p , p = { 0 , 1 , 2 , 3 } , given by (30)–(33) we obtain
μ 4 k + 0 = 3 2 k + 1 · ( 4 k + 1 ) ! ! ! ! ( 4 k + 3 ) ! ! ! ! · I 0 I 0 = 3 2 k + 1 · 1 3 · 5 7 · 9 11 4 k + 1 4 k + 3 · 1
μ 4 k + 1 = 6 4 k + 3 · ( 2 k + 1 ) ! ! ( 2 k + 2 ) ! ! · I 1 I 0 · 1 r = 6 4 k + 3 · 1 2 · 3 4 · 5 6 2 k + 1 2 k + 2 · π 2 2 K 1 / 2
μ 4 k + 2 = 3 2 k + 2 · ( 4 k + 3 ) ! ! ! ! ( 4 k + 5 ) ! ! ! ! · I 2 I 0 · 1 r 2 = 3 2 k + 2 · 3 5 · 7 9 · 11 13 4 k + 3 4 k + 5 · π 2 K 1 / 2 2
μ 4 k + 3 = 6 4 k + 5 · ( 2 k + 2 ) ! ! ( 2 k + 3 ) ! ! · I 3 I 0 · 1 r 3 = 6 4 k + 5 · 2 3 · 4 5 · 6 7 2 k + 2 2 k + 3 · 1 2 K 1 / 2
where n = 4 k + p with k 0 and we have used in (49) the relation for the complete elliptic integrals of the first and second kind, i.e.,
2 E 1 2 K 1 2 K 1 2 = π 2 ·
Equivalently, by rewriting (47)–(50) in a more concise form, we obtain
μ n | p = 6 n + 2 · ( n + 1 ) ! ! ! ! ( n + 3 ) ! ! ! ! · 1 , π 2 2 K 1 / 2 , π 2 K 1 / 2 2 , 1 2 K 1 / 2
where the numerical multipliers in the curly brackets are applied in the order of their occurrence (i.e., the number 1 is for μ n | 0 , and so on), in correspondence with the complete set of residues p = { 0 , 1 , 2 , 3 } modulo 4. The elements μ n | p in each of the classes with p = { 0 , 1 , 2 , 3 } are obtained for n running through the values of n p ( mod 4 ) .
As becomes clear from the above formulas, the only rational numerical coefficients are those of the zeroth class, i.e.,
μ n | 0 = 6 n + 2 · 1 3 · 5 7 · 9 11 · 13 15 n + 1 n + 3 , n 0 ( mod 4 ) .
The numerical coefficients in the other classes, i.e., μ n | 1 , μ n | 2 , and μ n | 3 , are irrational transcendental numbers.
The complete elliptic integral of the first kind K ( 1 / 2 ) 1 . 85407 appearing in the Formula (52) for the numerical coefficients μ n | p , or respectively in (45) for the mechanical moments J n | p , can be related to a variety of other special functions and transcendental constants, e.g., the gamma function Γ ( 1 / 4 ) 3 . 62561 , Gauss’s constant
G : = 2 π 0 1 d t 1 t 4 0.834627
or the lemniscate constant ω ˜ = π G 2 . 62206 . The corresponding relations, i.e.,
K 1 2 = π G 2 = Γ ( 1 / 4 ) 2 4 π = ω ˜ 2
allow us to immediately obtain, in addition to (52), the other three equivalent representations for the numerical coefficients μ n | p , n p ( mod 4 ) , p = { 0 , 1 , 2 , 3 } , i.e.,
μ n | p = 6 n + 2 · ( n + 1 ) ! ! ! ! ( n + 3 ) ! ! ! ! · 1 , 2 π 3 / 2 Γ ( 1 / 4 ) 2 , 8 π 2 Γ ( 1 / 4 ) 4 , 2 2 π Γ ( 1 / 4 ) 2
μ n | p = 6 n + 2 · ( n + 1 ) ! ! ! ! ( n + 3 ) ! ! ! ! · 1 , 1 2 G , 1 π G 2 , 1 π G
μ n | p = 6 n + 2 · ( n + 1 ) ! ! ! ! ( n + 3 ) ! ! ! ! · 1 , π 2 ω ˜ , π ω ˜ 2 , 1 ω ˜ .
Therefore, substituting the numerical values of the corresponding special functions and transcendental constants into the expressions (52), (56)–(58) we obtain
μ n | p 6 n + 2 · ( n + 1 ) ! ! ! ! ( n + 3 ) ! ! ! ! · 1 , 0.59907 , 0.456947 , 0.38138 .
The values of the first twenty numerical coefficients μ n | p divided into four residue classes modulo 4 expressed via Gauss’s constant G are presented in Table 1.

5.2. Calculation of Numerical Coefficients μ n for Hollow Mylar Balloon

For the hollow Mylar balloon the numerical coefficients μ n | p given by (20) can be calculated for n = 0 , 1 , 2 as
μ 0 = I 1 I 1 = 1 , μ 1 = I 2 I 1 · 1 r = 2 π G , μ 2 = I 3 I 1 · 1 r 2 = 2 π
whereas for n > 2 we have
μ n | p = I n + 1 | p I 1 · 1 r n = ( n 2 ) ! ! ! ! n ! ! ! ! · 1 , 2 π G , 2 π , 2 G , n p ( mod 4 ) , p = { 0 , 1 , 2 , 3 } .
More specifically, substituting directly the values of I p , p = { 1 , 2 , 3 , 4 } , given by (31)–(35), similarly to (47)–(50), we obtain
μ 4 k + 0 = ( 2 k 1 ) ! ! ( 2 k ) ! ! · I 1 I 1 = 1 2 · 3 4 · 5 6 2 k 1 2 k · 1 , k 1
μ 4 k + 1 = ( 4 k 1 ) ! ! ! ! ( 4 k + 1 ) ! ! ! ! · I 2 I 1 · 1 r = 3 5 · 7 9 · 11 13 4 k 1 4 k + 1 · 2 π G , k 1
μ 4 k + 2 = ( 2 k ) ! ! ( 2 k + 1 ) ! ! · I 3 I 1 · 1 r 2 = 2 3 · 4 5 · 6 7 2 k 2 k + 1 · 2 π , k 1
μ 4 k + 3 = ( 4 k + 1 ) ! ! ! ! ( 4 k + 3 ) ! ! ! ! · I 4 I 1 · 1 r 3 = 1 3 · 5 7 · 9 11 4 k + 1 4 k + 3 · 2 G , k 0
where n = 4 k + p . The values of the first twenty numerical coefficients μ n | p divided into four residue classes modulo 4 expressed via Gauss’s constant G are presented in Table 2.

6. Concluding Remarks

As can be easily seen, by comparing the corresponding expressions, the numerical coefficients μ n and μ n for the solid and hollow Mylar balloons are related to each other through the formula
μ n | p = 6 n + 2 · I n + 4 | p I 0 · 1 r n + 4 = 6 n + 2 · I 1 I 0 · 1 r · I n + 4 | p I 1 · 1 r n + 3 = 3 n + 2 · 1 G · μ n + 3 | p
where we have
n p ( mod 4 ) , p + 3 p ( mod 4 ) , p = { 0 , 1 , 2 , 3 } , p = { 3 , 0 , 1 , 2 } .
From the above expressions it follows that the pair of sequences of numerical coefficients { μ 0 , μ 1 , , μ n , } and { μ 3 , μ 4 , μ n + 3 , } as well as the respective pair of sequences of mechanical moments { J 0 , J 1 , J n , } and { J 3 , J 4 J n + 3 , } can be regarded as pairs of sequences related by cyclically shifted set of residue classes, as is illustrated by the diagrams (n takes the values from the respective class n | p , n p ( mod 4 ) )
p = { 0 , 1 , 2 , 3 } p = { 3 , 0 , 1 , 2 }
{ μ n | 0 , μ n | 1 , μ n | 2 , μ n | 3 } 3 n + 2 · 1 G · μ n + 3 | 3 , μ n + 3 | 0 , μ n + 3 | 1 , μ n + 3 | 2
or, equivalently, by the infinite sequence of diagrams
{ μ 0 | 0 , μ 1 | 1 , μ 2 | 2 , μ 3 | 3 } 3 2 G μ 3 | 3 , 1 G μ 4 | 0 , 3 4 G μ 5 | 1 , 3 5 G μ 6 | 2
{ μ 4 | 0 , μ 5 | 1 , μ 6 | 2 , μ 7 | 3 } 3 6 G μ 7 | 3 , 3 7 G μ 8 | 0 , 3 8 G μ 9 | 1 , 3 9 G μ 10 | 2
{ μ 8 | 0 , μ 9 | 1 , μ 10 | 2 , μ 11 | 3 } 3 10 G μ 11 | 3 , 3 11 G μ 12 | 0 , 3 12 G μ 13 | 1 , 3 13 G μ 14 | 2
and so on. As can be seen from above diagrams, the corresponding elements μ n | p and μ n + 3 | p in each pair of the shifted classes (for n 0 ) differ by specific numerical factors (transcendental numbers) arranged in a sequence as
3 G 1 2 , 1 3 , 1 4 , , 1 n + 2 , 1.79721 , 1.19814 , 0.898605 , , 3.59442 n + 2 ,
where G 0 . 834627 . The presence of the three fundamental mathematical constants π , ω ˜ , and G, as well as the combinatorial character of the above expressions suggests unambiguously that we are dealing here with deep nontrivial relationships in nature.
Let us also state that this work was inspired by the observation based on the longstanding experience which the authors have with Mylar balloons that leads them to recognition of the fundamental fact that all its geometrical characteristics are related to the moments specified in (10). This observation has been immediately expanded to cover mechanical moments in both solid and hollow cases. While the first few moments (in both geometrical and mechanical settings) have direct interpretations as mass, surface area, volume, profile curve, and so on, this is not so obvious in the case of the higher moments. They appear, e.g., in the description of the Newtonian attraction of solid bodies in the non-relativistic theory of gravitation and in the classical potential theory (cf., e.g., [14,15]).
One should also notice the resemblance of the recurrence relationships (21) and the Legendre polynomials which definitely deserves a further study as well.

Author Contributions

Conceptualization, V.K., V.I.P. and I.M.M.; methodology, V.K., V.I.P. and I.M.M.; formal analysis, V.K., V.I.P. and I.M.M.; investigation, V.K., V.I.P. and I.M.M.; writing—original draft preparation, V.K.; writing—review and editing, V.K., V.I.P. and I.M.M.; visualization, V.K., V.I.P. and I.M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This paper presents part of the results obtained during the realization of the joint research project #17 under the title “Classical and quantum models of deformable bodies on curved surfaces”. The project was accomplished within the realm of the scientific cooperation between the Bulgarian and Polish Academies of Sciences for the period of 2022–2023. The authors are thankful to the above-mentioned Institutions for their support of the exchange visits conducted in both directions. The authors are also thankful to the reviewers for their valuable remarks and suggestions helping to improve the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Paulsen, W. What is the shape of a Mylar balloon? Am. Math. Mon. 1994, 101, 953–958. [Google Scholar] [CrossRef]
  2. Smalley, J. Development of the E-Balloon; Technical Report AFCRL-70-0543; National Center for Atmospheric Research: Boulder, CO, USA, 1970. [Google Scholar]
  3. Tang, J.; Pu, S.; Yu, P.; Xie, W.; Li, Y.; Hu, B. Research on trajectory prediction of a high-altitude zero-pressure balloon system to assist rapid recovery. Aerospace 2022, 9, 622. [Google Scholar] [CrossRef]
  4. Kawaguchi, M. The shallowest possible pneumatic forms. Bull. Int. Assoc. Shell Struct. 1977, 18, 3–11. [Google Scholar]
  5. Mladenov, I.; Oprea, J. The Mylar balloon revisited. Am. Math. Mon. 2003, 110, 761–784. [Google Scholar] [CrossRef]
  6. Wolfram Language. Available online: https://en.wikipedia.org/wiki/Wolfram_Language (accessed on 5 June 2023).
  7. Mladenov, I. On the geometry of the Mylar balloon. Comptes Rendus L’acad. Bulg. Sci. 2001, 54, 39–44. [Google Scholar]
  8. Pulov, V.; Hadzhilazova, M.; Mladenov, I. The Mylar balloon: An alternative description. Geom. Integr. Quant. 2015, 16, 256–269. [Google Scholar]
  9. Kovalchuk, V.; Mladenov, I.M. Classical motions of infinitesimal rotators on Mylar balloons. Math. Methods Appl. Sci. 2020, 43, 9874–9887. [Google Scholar] [CrossRef]
  10. Gradstein, I.S.; Ryzhik, I.M. Tables of Integrals, Series, and Products, 7th ed.; Jeffrey, A., Zwillinger, D., Eds.; Academic Press: Oxford, UK, 2007. [Google Scholar]
  11. Celletti, A. Stability and Chaos in Cellestial Mechanics; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  12. Westfall, P. Kurtosis as Peakedness. Am. Stat. 2014, 68, 191–195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Byrd, P.; Friedman, M. Handbook of Elliptic Integrals for Engineers and Scientists, 2nd ed.; Springer: New York, NY, USA, 1971. [Google Scholar]
  14. Ramsey, A. Newtonian Attraction; Cambridge University Press: Cambridge, UK, 1961. [Google Scholar]
  15. Sterne, T. An Introduction to Celestial Mechanics; Interscience: New York, NY, USA, 1960. [Google Scholar]
Figure 1. 3D view of the inflated Mylar balloon (left) and its cross-section through the symmetry axis O Z in the plane X O Z (right). N and S are the North and South Poles, E denotes the Equator of this surface of revolution, and O is the center of the Mylar balloon. The profile curve is defined through the relation z = z ( x ) , whereas τ = 2 z ( 0 ) is the thickness of the Mylar balloon and a, r are its deflated and inflated radii, respectively. Alternatively, the variable u describes the parametrization of the Mylar balloon, i.e., x = x ( u ) and z = z ( u ) with u [ 1 , 1 ] .
Figure 1. 3D view of the inflated Mylar balloon (left) and its cross-section through the symmetry axis O Z in the plane X O Z (right). N and S are the North and South Poles, E denotes the Equator of this surface of revolution, and O is the center of the Mylar balloon. The profile curve is defined through the relation z = z ( x ) , whereas τ = 2 z ( 0 ) is the thickness of the Mylar balloon and a, r are its deflated and inflated radii, respectively. Alternatively, the variable u describes the parametrization of the Mylar balloon, i.e., x = x ( u ) and z = z ( u ) with u [ 1 , 1 ] .
Mathematics 11 02646 g001
Figure 2. Dependencies of I n I n ( x ) , n = 1 , , 10 (left, top to bottom), including the profile curve z ( x ) = I 2 I 2 ( x ) , and the crimping factor C ( x ) (right), for the Mylar balloon with radius r = 1 .
Figure 2. Dependencies of I n I n ( x ) , n = 1 , , 10 (left, top to bottom), including the profile curve z ( x ) = I 2 I 2 ( x ) , and the crimping factor C ( x ) (right), for the Mylar balloon with radius r = 1 .
Mathematics 11 02646 g002
Table 1. First twenty mechanical moments J n | p (solid case) divided into four residue classes modulo 4 expressed via Gauss’s constant G: J n | p = μ n | p m r n , m = 1 3 π 2 r 3 ρ G , n p ( mod 4 ) , p = { 0 , 1 , 2 , 3 } (cf. (57)).
Table 1. First twenty mechanical moments J n | p (solid case) divided into four residue classes modulo 4 expressed via Gauss’s constant G: J n | p = μ n | p m r n , m = 1 3 π 2 r 3 ρ G , n p ( mod 4 ) , p = { 0 , 1 , 2 , 3 } (cf. (57)).
p = 0 p = 1 p = 2 p = 3
n μ n | 0 n μ n | 1 n μ n | 2 n    μ n | 3
011 1 2 G 2 9 10 π G 2 3   4 5 π G
4 5 21 5 9 56 G 6 7 20 π G 2 7  16 45 π G
8 9 77 9 15 176 G 10 77 390 π G 2 11  96 455 π G
12 39 539 13 7 128 G 14 231 1768 π G 2 15 256 1785 π G
16 221 4389 17 189 4864 G 18 209 2210 π G 2 19 512 4851 π G
Table 2. First twenty mechanical moments J n | p (hollow case) divided into four residue classes modulo 4 expressed via Gauss’s constant G: J n | p = μ n | p m r n , m = π 2 r 2 ρ , n p ( mod 4 ) , p = { 0 , 1 , 2 , 3 } (cf. (62)–(65)).
Table 2. First twenty mechanical moments J n | p (hollow case) divided into four residue classes modulo 4 expressed via Gauss’s constant G: J n | p = μ n | p m r n , m = π 2 r 2 ρ , n p ( mod 4 ) , p = { 0 , 1 , 2 , 3 } (cf. (62)–(65)).
p = 0 p = 1 p = 2 p = 3
n μ n | 0 n μ n | 1 n μ n | 2 n    μ n | 3
011 2 π G 2 2 π 3   2 G 3
4 1 2 5 6 5 π G 6 4 3 π 7  10 G 21
8 3 8 9 14 15 π G 10 16 15 π 11  30 G 77
12 5 16 13 154 195 π G 14 32 35 π 15 26 G 77
16 35 128 17 154 221 π G 18 256 315 π 19 442 G 1463
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kovalchuk, V.; Pulov, V.I.; Mladenov, I.M. Mylar Balloon and Associated Geometro-Mechanical Moments. Mathematics 2023, 11, 2646. https://doi.org/10.3390/math11122646

AMA Style

Kovalchuk V, Pulov VI, Mladenov IM. Mylar Balloon and Associated Geometro-Mechanical Moments. Mathematics. 2023; 11(12):2646. https://doi.org/10.3390/math11122646

Chicago/Turabian Style

Kovalchuk, Vasyl, Vladimir I. Pulov, and Ivaïlo M. Mladenov. 2023. "Mylar Balloon and Associated Geometro-Mechanical Moments" Mathematics 11, no. 12: 2646. https://doi.org/10.3390/math11122646

APA Style

Kovalchuk, V., Pulov, V. I., & Mladenov, I. M. (2023). Mylar Balloon and Associated Geometro-Mechanical Moments. Mathematics, 11(12), 2646. https://doi.org/10.3390/math11122646

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop