# Growth and Form

Is There Anything Interesting to Say About Dimensions in Physics?
1,103
Full-Text Views:
134
Citations (Scopus):
0
Citations (Crossref):
0

## 1. INTRODUCTION

### 1.1. Dimensional Analysis and Reduction of Complexity

Any quantitative science relies on the necessity of introducing specific quantities to be defined through a measure, implying a protocol, tools and units. The geometry, if intended in its etymology γεωµητρια (measurement of land) is such an example. The measure itself may be intended through different procedures, the most natural of which is the comparison with a reference quantity. Regarding geometry, the length is what is intended to be specified in quantitative terms.

We define in abstract terms the geometrical dimension [L], which is a primitive concept, as stressed below. Such a concept is richer than it may appear and has evolved through the concept of dimension illustrated below.

The point has dimensions zero [L0], the line 1 [L1], the area 2 [L2] and the volume 3 [L3]. The reference quantity can be a unit sub-segment, a unit square a unit cube (see Fig. 1).

What we have described so far is very elementary. According to this picture a length is just the sum of a certain number (3 in Fig. 1) of the chosen reference measure, say 3 s, if s is the reference measure. Accordingly surface and volume scale as 32, 33 respectively.

The essence of Physics is the experiment, which lies on the concept of measurement. A prerequisite for the process of measure is the development of a protocol, according to which one attributes a gauge to characteristic quantities, thus getting a numerical meaning for pure concepts, often referred as “fundamental” or “primary” dimensions (PD), like length [L], time [T] and mass [M], if we limit ourselves to mechanical problems.

Let us remain, for the moment, within the land of pure abstraction and think in terms of dimensions, without referring to any specific system of units. Any physical quantity like energy [E], velocity [U]1, force [F]. . . can be expressed as “dimensional products” of the aforementioned PD and for example, we find

[E]=[ML2T2][U]=[LT1][F]=[MLT2] (1)
where the exponent defines how a “fundamental” dimension enters in the specification of another physical quantity. The velocity is indeed just the ratio between space (length) and time, the dimensional identification of the force follows from the second Newton law and E reflects the dimensional content of the kinetic energy.

The question arises on what fundamental does mean. In philosophical terms we may refer to length, time and matter as “pre-categories” through which “reality” can be organized within a deductive context.

We simply rephrase the point of view of the Kantian philosophy to state that L, T, M are the tools which are “triggered” by the experience to formulate “a priori analytical propositions” and transform perceptions into a “phenomenon”.

In principle perception in physics can be organized not only through what we have defined PD, but through other “means” as, for example, E, F, U. The idea we like to convey is that the meaning of “fundamental” is just a matter of convention, as shown below in an appropriate mathematical context.

In general, any mechanical quantity will be specified by the dimensional equation

[Q]=[LαTβMγ] (2)

Where α, β, γ are not necessarily positive and integers.

We provide an idea of how the process of assigning dimensions proceeds. An example is provided by the definition of the viscosity of a fluid, which can be viewed as its resistance to flow. The concept is illustrated in Fig. 2, where we have reported a liquid within two plates (with surfaces A), the upper of which moves and drags lamina of flow, whose velocity increases with the height (vertical direction in Fig. 2). We can accordingly state that the liquid resistive force is

FV=μAΔuΔy (3)

The dimensional content of the coefficient μ is readily obtained

μ=FVA(ΔuΔy)1[μ]=[ML1T1]

Before proceeding further, we make a step forward in the process of abstraction and define the “space of dimensions”, spanned by the unit vectors

l_(1,0,0),t_(0,1,0),m_(0,0,1) (4)

In this space the quantity Q is written as

q_(α,β,γ)=αl_+βt_+γm_ (5)

It is therefore worth noting that, keeping three linearly independent triples of primary units

q_1=α1l_+β1t_+γ1m_q_2=α2l_+β2t_+γ2m_q_3=α3l_+β3t_+γ3m_ (6)

They are invertible

(l_t_m_)=(α1β1γ1α2β2γ2α3β3γ3)1(q_1q_2q_3) (7)
because the determinant of the associated matrix is not vanishing.

Accordingly, the statement of fundamental and derived dimensions is, at least in mathematical terms, only matter of a choice, as discussed more accurately below and in the forthcoming sections.

For example, we find that by keeping energy, force and velocity as “fundamental”, the following dimensional correspondence with time, mass and length can be established

t_=e_f_v_[T]=[EF1V1]m_=e_2v_[M]=[EV2]l_=e_f_[L]=[EF1] (8)

There is however nothing deep in what we have noted, which, said in these terms, is just a reshuffling of symbols, the forthcoming sub-section contains a more substantive step further.

### 1.2. Dimensional Analysis and the Principle of Homogeneous Dimensions

The topics and concepts, we have discussed so far, can be organized into what is recognized as dimensional analysis (DA), which provides a powerful mean to get a preliminary insight into the “scaling” of physical quantities, within a given phenomenological context.

It is worth to adding, as important premise, the Principle of Dimension Homogeneity, which, albeit tacitly assumed, is now explicitly worded as follows

All additive terms in a dimensional equation have the same dimensions

This statement clarifies the obscure wording (dimensional product), used before. We accordingly note that power is the ratio between work and time, work is the same as energy (at least in terms of dimensions), therefore we find

[P]=[ET1]=[ML2T2T1]=[ML2T3]

A paradigmatic example of application of DA is associated with the drag forces, emerging as consequence of the viscosity and we use dimensionality argument to guess an analytical form for the force itself. By analytical we mean how the various physical quantities defining the problem are embedded to get a force.

A preliminary analysis indicates that all the quantities of interest can be arranged as in the table shown below, where we have included velocity (U) Surface (S), density (ρ) and viscosity (μ) along with the relevant PD content.

USρμFL12311T10012M00111 (9)

On the right column we have added the force and the relevant dimensions to be matched.

We have four physical quantities and 3 primary dimensions. This gives us 4 choices to aggregate U, S, ρ, µ to end up with a force, one of which is shown below:

[F]1=[UαρβSγ]l_2t_+m_=(α3β+2γ)l_+βm_αt_α=2,β=1,γ=1 (10)

We can accordingly conclude that the force may be cast as

F112ρv2S (11)

The factor 12 is just due to the fact that 12ρv2 is precisely the kinetic pressure of the fluid.

Eq. (11) is dimensionally correct, but we cannot exclude any proportionality coefficient, which can be expressed as function of a non-dimensional quantity emerging from the combination of the variables entering the definition of our problem.

To clarify this point and attempt a solution, we consider different arrangements of the reference quantities and assume that

[F]2=[μαρβSγ]l_2t_+m_=(α+β)m_+(2γα3β)l_αt_α=2,β=1,γ=0F3=[μαρβUγ]α=2,β=1,γ=0F4=[μαSβUγ]α=1,β=1/2,γ=1 (12)

F2=F3μ2ρF4μvS (13)

Together with dimensional, non-dimensional quantities play a pivotal role in the analysis of physical phenomena. A non-dimensional quantity naturally emerges from Eq. (13) as the ratio between forces, i.e.

F4F2,3ρvSμ (14)

This is a non-trivial quantity and is recognized as the Reynolds number:

Re=ρvSμ (15)

The other ratios yield

F1F2,312Re2,F1F412Re (16)

We have not developed any theory to specify how Re should be appended to the definition of force. We can however assume that Eq. (11) is true except that it should be multiplied by a not specified function of the Reynolds number, i.e.

F=cd2ρv2Scd=f(Re) (17)
where cd is the drag coefficient.

The practical outcome of this discussion is that if we are interested to carry out an experiment to specify cd we are not obliged to perform four groups of different measures, but we can just study the behavior of 2F/(ρ v2 S) vs. Re.

From the discussion in this introductory section, Dimensional analysis emerges as a working tool useful to analyze a physical problem and getting scale relationships between wisely chosen groups of variables. A significant byproduct of the procedure is the identification of key quantities useful to reduce the complexity of the problem under study. This is certainly true, but we believe it is much more, as we will try to convey in the forthcoming sections.

## 2. THE BUCKINGHAM THEOREM AND ITS CONSEQUENCES

### 2.1. Ship Design and the Froude Number

In the previous section we have mixed up and anticipated different points, which we are going to develop in a more systematic way. We have mentioned the procedure of combining quantities like surfaces, density, viscosity and velocity to derive either the drag force and a non-dimensional number inferred as a byproduct of this melting criterion. The previous ideas trace back to the seminal works of Rayleigh [1,2] and Buckingham [3,4] which will be discussed below.

The Rayleigh method, even though brilliant, amounts to grouping a number of physical quantities, fix primer dimensions, combine the exponents and match them to another specific quantity. The procedure strongly relies on the principle of dimensional homogeneity, on experience and physical intuition as well.

To illustrate the method, an example is better than any comment. We ask therefore the following question: in a given physical phenomenon the velocity depends on the pressure [p] ≡ [ML−1T−2], and density [ρ] = [ML−3], can we use DA to guess the dependence of velocity on the previously quoted quantities?

To answer the question, we require that

[U]=[pαργ] (18)

Which, embedded into the dimensional vector space, yields the relation

l_t_=(α+γ)m_2αt_(α+3γ)l_ (19)

The principle of identity between vectors (i.e. the dimension homogeneity condition), imposes that

α+γ=02α=1α+3γ=1 (20)
yielding α=12,γ=12. We can therefore conclude that such a dependence occurs through pressure and density, according to the relation
v=pρ (21)

We have “discovered” an obvious identity, which becomes only slightly less trivial, if we look at the problem from a different perspective, namely by asking whether DA can help us to infer the dependence of the velocity of propagation of a pressure wave through a liquid, in terms of its elasticity (K) and density. The use of the outlined procedure yields v=Kρ, which is the same as Eq. (21), if we note that K (associated with the Young modulus of the liquid) has the dimensions of a pressure.

The next example is the guessing of a dimensionless number (the Froude number) of central importance in the design of ships.

The analysis of the motion of ships during their traveling at sea, at lakes or at rivers is an extremely difficult problem. It should be therefore noted while traveling the boat-wave interaction usually determines the formation of steep diverging bow and stern waves, possibly breaking. Furthermore the ship geometry and the velocity are responsible for a complex phenomenology, as, for example, the generation of vorticity in the wave-breaking region, the appearance of free-surface scars, and the entrainment of air, to name a few [5].

To better appreciate the physical nature of the previously quoted phenomena and the relationship with the Reynolds number, we go back to fundamental principles.

A body in a fluid may be subject to different forces, including those of viscous FD and gravitational Fg nature. Whatever forces we consider, the second Newton laws imposes that their effect on the body will be a change of its momentum, we can accordingly write (FR denotes any the other force acting on the body)

FD+Fg+FR=ΔpΔt (22)

Accordingly we can write

ND+Ng+NR=κNqFqΔpΔt (23)
where Nq (for q = D, g, R) are characteristic non-dimensional numbers. To make a comparison between the first two, we use the theorem of kinetic energy2 and set
ΔpΔt=ΔTΔLT=12Mv2 (24)

Thus finding, from Eq. (11)

NDρv2SMv2L1=ρvLLMLT1L1=ρvSμ=ReNgMgMv2L1=gLv2=(NF)2NF=vgLFroudeNumber (25)
and NR is a number depending on the specific force FR, as discussed later in the article. We have so far just combined dimensions obtaining non-dimensional quantity, which, if not embedded within a proper physical, dynamical and design context, are hardly appreciated. Hence, talking about NF independently of gravity wave dynamics and ship sailing does not reveal much of its importance [6]. These points will be touched later at the end of the article.

The discussion we have developed so far may appear simplistic, we will see at the end of the paper how the nondimensionalization of the Navier-Stokes equation leads to analogous conclusions.

### 2.2. The Buckingham Π Theorem

According to the previous discussion it has emerged that, in correspondence of a certain number of dimensional variables, there exist a number of non-dimensional quantities. It is therefore quite natural to specify how many of these numbers can be found in correspondence of the list of dimensional terms.

The answer is contained in an important theorem of dimensional analysis, proposed by Buckingham [3,4] and asserting that

“Given a list of n dimensional quantities with k primary units, the number of corresponding non dimensional products is n – k”

Or, in alternative form:

“A dimensionally homogeneous equation, involving n physical variables, all expressible in terms of k independent fundamental physical quantities, can be expressed in terms of p = n – k dimensionless parameters”

A preliminary example is certainly useful to clarify what the previous statements mean. We consider therefore the following problem, given an equation

F=f(w,h,μ,ρ,U) (26)

Accounting for the dependence of F with the dimension of a force on

[w]=[h]=[L],[μ]=[ML1T1],[ρ]=[ML3],[U]=[LT1] (27)

Namely a total of 6 variables. Because we have chosen M, L, T as primary dimensions (see below), we should be able of working out 3 non-dimensional parameters, characterizing the same equation. The above general formulation is sufficient to work out a solution, in any case if like a less abstract formulation, we can reformulate it in the following way [7]:

a thin rectangular plate having a width w and a height h is located so that it is normal to a moving stream of fluid. We can guess that the drag force, exerted by the fluid on the plate, depends on the plate dimensions (w, h), on the fluid viscosity and density (µ, ρ respectively) and on the velocity of the fluid approaching the plate. We can determine a suitable set of non-dimensional quantities (Π′s, see below) to study this problem experimentally.

The dimensional matrix, characterizing the problem is drawn below

FwhμρUL111131T200101M100110 (28)

Since we are looking for non-dimensional quantities, we proceed as follows.

The first step is the definition of a number of “repeating variables”3 to be embedded with the others to get the desired non-dimensional terms. The number of repeated variables are (usually) the same as the primary dimensions, we choose (w, U, ρ) and set (the symbol Π stands for product)

Π1=Fwα1Uβ1ργ1(1+γ1)m_+(1+α1+β13γ1)l_(β1+2)t_Π2=hwα2Uβ2ργ2γ2m_+(1+α2+β23γ2)l_β2t_Π3=μwα3Uβ3ργ3(1+γ3)m_+(1+α3+β33γ3)l_(1+β3)t_ (29)

The Buckingham Π′s are required to be null vectors Π ≡ (0, 0, 0) in the dimension vector space, therefore equating to zero the coefficients of (l_,t_,m_), we find

(α1β1γ1α2β2γ2α3β3γ3)=(221100111) (30)
and the relevant list of non-dimensional quantities is given below
Π1=Fw2U2ρ,Π2=hw,Π3=μwUρ (31)

Which, even though in a more organized way, yields the same information we have drawn in the previous parts of the article. It is furthermore evident that Eq. (26) can be organized in another functional dependence involving the three “numbers”, which can be written as

Π1=Φ(Π2,Π3) (32)

The “proof” of the theorem is hidden in the example itself.

Without developing any formal discussion, we note that, once we have established a number of r repeating variables out of n, we can use them as pivot to exploit the Rayleigh method to get, in correspondence of any of the k = nr remaining variables, a quantity with the inverse dimension of one corresponding to the set k. In other words, following the Rouché-Capelli theorem, k is the solutions subspace dimension of a homogeneous linear system with n variables and m equations, being r = min(rank(A),m)4. We have found indeed that w−2U−2ρ−1, defining Π1 in Eq. (31), has the dimensions of the inverse of a force.

There are subtleties in the procedure requiring much care and experience. Some obvious prescriptions are that the choice of the n variables should include quantities which can be kept constant, and all the variables be independent. In the previous example, we have considered separately the geometrical dimensions (w, h) of the plate, they, even though with the same dimensions, are independent each other. Selecting S = w h and w would be inappropriate, because not independent.

The choice of keeping separated w, h has been motivated by the fact that if, for example, wh we can expect significant variation in the drag than the case wh, for which the variable can be reduced by choosing, as independent, S only. The assumption of MLT as fundamental dimensions is not compulsory, we can use another set of units, provided that we preserve the homogeneity condition and even get some computational advantage.

The following example [7] can be useful to understand how things are going.

An open, cylindrical tank having a diameter D is supported around its bottom circumference and is filled to a depth h with a liquid having a specific weight γ. The vertical deflection, δ, of the center of the bottom is a function of D, h, d, γ, and E, where d is the thickness of the bottom and E is the modulus of elasticity of the bottom material.

In mathematical terms we have:

δ=f(D,h,d,γ,E)γ=ρg (33)

Furthermore, E has the same dimensions of the pressure. If we choose F, L, T as primary units, we get

δDhdγEL111132T000000F000011 (34)

Inspecting the table, we note that time, in this system of units, does not play any role. We can therefore choose two repeating variables only, for example D, γ and eventually find 4 (6 – 2) non-dimensional Π, namely

Π1=δD,Π2=hD,Π3=dD,Π4=EDγ (35)
and
δD=Φ(hd,dD,EDγ) (36)

A further possibility to be considered is the choice, for the same problem, of different repeating variables. We consider indeed the case of the pressure dropper unit length, along a pipe of diameter D, through which is flowing a liquid with viscosity μ, density ρ and velocity U. We are accordingly looking for a dependence of the type

P=f(D,ρ,μ,U) (37)
where [P′] = F/L3. The repeating variables may be for example D, U, ρ yielding
PDρU2=Φ1(ρUDμ) (38)
or
PD2μU=Φ2(ρUDμ) (39)

The two results are clearly equivalent, provided that

Φ2(Re)=ReΦ1(Re) (40)
and both useful to test the validity of the so-called Darcy-Weisbach equation [8], according to which
P=12fDρv2D (41)
with fD being the friction coefficient, which, for laminar fluid motion, is just given by
fD64Re (42)

### 2.3. Star Oscillations and Water Droplets

Before closing this section, we like to present a few examples from fields not strictly associated with hydraulics, to this aim we consider the example of star oscillations. There is no doubt that sun and other stars are subject to small oscillations. We can use the DA, developed so far, to infer how the oscillation frequency depends on the characteristic density and radius of the star itself.

The gravitational force is the most significant force acting on the star, and the driving quantity of this force is the gravitational constant G, it is however not difficult to argue that the problem can be reduced to the following equation [8]:

ω=Ω(ρ,R,G) (43)

Here R is the radius of the star and G is the gravitational constant with the following dimensions

[G]=[M1L3T2]=[ρ1T2] (44)

The answer is therefore straightforward, noting indeed that [ω] = [T–1], we can conclude that

ω=κGρ (45)
where κ is a numerical constant.

The conclusion is, within certain limits surprising, the star oscillation depends only on its density and on the gravitational constant, but not on its radius. The “bewilderment” is however ungrounded and the relationship of Eq. (45) with the Gauss theorem5 is easily understood.

The discussion can be extended to small spherical water droplets. We may use an analogous argument to infer the oscillation frequency in terms of the surface tension τ = [MT–2], density and radius, in this case too, not much algebra is involved and we end up with [9]

ωτρR3 (46)

The tacit assumption in deriving the identity (46) was that the gravity force can be neglected with respect to those due to the surface tension. This statement can be expressed through the condition

GρτρR3 (47)
which in turn yields an upper limit on the water droplet radius to neglect gravity force effects, namely
R(τρ2G)1/3 (48)
and taking into account the water characteristic parameters ρ = 103 kg/m3, τ = 0.07 N · m−1 and G6.674081011m3kgs2, we end up with R ≪ 10 m (not exactly a droplet).

### 2.4. Surface Tension and Capillary Waves

As a byproduct of this discussion, we consider the gravity and capillary wave oscillation frequencies.

We remind that a gravity wave is a wave formed at the boundary of two media, when the gravity tends to restore the equilibrium. If we denote by λ, k=2πλ its wavelength and wave vector respectively, the associated frequency is readily obtained as

ωGgk (49)
with g = 9.81 ms−2.

The capillary waves are short wavelength water waves, in which the equilibrium is restored by forces associated with the surface tension. We obtain after replacing e.g. R3 in Eq. (47) with k−3:

ωGτk3ρ (50)

Equating Eq. (49) and Eq. (50), we find for the wave vector

k=ρgτ (51)
which fixes the upper limit of the capillary wave wavelength (about 1 cm in the case of water).

Going back to the Froude number, we have now a physical argument to justify the relevant introduction. The wave group V G(g) and phase V G(p) velocities associated to the dispersion relation in Eq. (49), write

V G(p)=ωGkgkV G(g)=ωGk12gk (52)

In shallow water, with depth d, the two velocities coincide and are provided by

VG=gd (53)

The Froude number arises therefore after confronting the wave gravity velocity with that of the flow velocity.

A few more comments are however in order, the problem is sometimes a little more complicated than our smart assumptions. To get a more general and reasonable picture of gravity waves, we make the assumption that they depend on the parameters listed in the table below (encompassing water depth d, wave height h, surface tension and gravity acceleration constant g) along with the relevant dimensions [10]

dhρτgkωL1130110T0002201M0011000 (54)

The Buckingham theorem yields the following functional dependence

ω2gk=Φ(hk,dk,ρgτk2) (55)

The non-dimensional quantity Bo=ρgτk2 is the Bond number, which suggests that when the wavelengths are long (Bolarge), the influence of the surface tension can be neglected. If furthermore d k << 1 (shallow water), expanding Ф at lowest order in d k we obtain the dispersion relation (53).

As final remark we stress that the group velocity associated with the dispersion relation (55) is written as

V G(g)=12gk(S+kSΦ)S=Φ(hk,dk,ρgτk2)Φ=dΦdk (56)

In this section we have given a general idea of how further we can go by mixing physical intuition, DA, common sense . . . in understanding physical phenomenology with relatively easy math. The forthcoming section regarding the use of DA in differential equations will further corroborate what has been obtained so far.

## 3. DIMENSIONAL ANALYSIS, DIFFERENTIAL EQUATIONS AND HOMOGENEOUS FUNCTIONS

We make, in the following, a step forward towards the understanding of the importance of a non-superficial view to dimension analysis in physical problems, by discussing the relevant importance within different context including quantum mechanics, general relativity, mathematics, …

Let us start with an important observation framing dimensional analysis within the context of scale invariance. We reformulate the principle of dimensional homogeneity as the following condition. Consider a dimensional functional dependence of the type reported in Eq. (29) (e.g. Π1=Fwα1Uβ1ργ1) and note that

(ww+UU+ρρ)Π1=(α1+β1+γ1)Π1 (57)

This statement is less trivial than it might appear, since it implicitly contains the dimensional homogeneity. Namely if we change by an amount B each of the variables w, U, ρ the invariance is ensured if F changes by B(α1+β1+γ1)

Accordingly, any function f (x1, x2, ..., xn) is dimensionally homogeneous if

Bnf(x1,x2,...,xn)=f(Bx1,Bx2,...,Bxn) (58)

Which, after deriving with respect to B, we obtain

nBn1f(x1,x2,...,xn)=s=1nx s'(x s')f(x 1',x 2',...,x n')x s'=Bxs (59)
after keeping B = 1, Eq. (59) yields the so called Euler equation for homogeneous functions
nf(x1,x2,...,xn)=s=1nxsf(x1,x2,...,xn)xs (60)

Let us now replace the rhs with a derivative with respect to a parameter η and transform it into the evolution equation

ϕ(x1,x2,...,xn;η)η=s=1nxsϕ(x1,x2,...,xn;η)xsϕ(x1,x2,...,xn;0)=f(x1,x2,...,xn) (61)

The solution of the previous Cauchy problem is

ϕ(x1,x2,...,xn;η)=eηs=1nxsxsf(x1,x2,...,xn) (62)

The notation (x1, x2, ..., xn; η) means that all the x-variables are implicit functions of η.

The explicit solution of Eq. (62) is obtained as a consequence of the identity [11,12]

eλxxf(x)=f(eλx) (63)
in the form
ϕ(x1,x2,...,xn,η)=f(eηx1,eηx2,...,eηxn) (64)

If we assume that the initial function of our problem is simply a monomial, namely

f(x1,x2,...,xn)=s=1nx sks (65)
the solution of Eq. (64) can eventually be written as
ϕ(x1,x2,...,xn;η)=eη(k1+...+kn)s=1nx sks (66)

The Buckingham theorem can now be formulated by using the more general framework of homogeneous functions.

We first remind the fundamental dogma of physics, let’s say the zero-th principle, which has been the leitmotif of the previous discussion. It states that the physical laws are independent of the chosen system of units. This amounts to saying that physical laws, to make any sense, must be homogeneous in all dimensions.

We consider therefore the following equation accounting for a physical law, namely

x0=f(x1,x2,...,xn) (67)
with xi being physical quantities in dimensional terms as
xi=[Lαi,LMαi,MTαi,T] (68)

Therefore, any change λ associated with one of the dimensions (say L), produces the transformation

x i'=λαi,L[Lαi,LMαi,MTαi,T] (69)

We get therefore

x 0'=f(x 1',x 2',...,x n') (70)

By keeping the derivative with respect to λ of both sides of Eq. (70) we obtain

α0,LλαL1x0=s=1nαs,Lxsxsf(x1,...,xn) (71)

Repeating the same procedure for the same change in the other dimensions and keeping λ = 0, we obtain

s=1nαs,Lf(x1,...,xn)ln(xs)=0s=1nαs,Tf(x1,...,xn)(lnxs)=0s=1nαs,Mf(x1,...,xn)(lnxs)=0(lnxs)=xsxs (72)

This means that we have obtained 3 independent equations fixing the dimensional variables; the remaining n − 3 are non-dimensional. It is evident that, if the number of units is larger, say m, the number of non-dimensional variables is nm. See the forthcoming section for further comments.

## 4. BUCKINGHAM THEOREM AND VIBRATIONS IN CRYSTALS

### 4.1. The Buckingham Theorem and Exponential Matrix Transformations

In the introductory section of the article, we have mentioned a kind of space of physical dimensions, realizing an abstract geometrical structure in which the corresponding vectors are specified as linear combinations of the fundamental “units” l_,t_,m_. The combination coefficients correspond to the exponents of the ordinary dimensional equations.

The transition from a triple of fundamental units to another is fixed by an ordinary vector transformation, therefore the transition from L, T, M to energy, viscosity and surface tension is specified by

e_=2l_2t_+m_μ_=l_t_+m_σ_=2t_+m_ (73)

In the following we use the notation (borrowed from [13,14]) to put Eq. (73) and their analog in the form of a generalized exponential equation, which reads

(E,M,Σ)=(L,T,M)AˆA=(221111021) (74)

We have exploited a notation in which a row vector is raised to a matrix. The operation leading to the explicit form of the transformation is obtained by the use of the following rule

[X1,X2]=[x1,x2](ξ1,1ξ1,2ξ2,1ξ2,2)=[x 1ξ1,1x 2ξ1,2,x 1ξ2,1x 2ξ2,2] (75)

The previous rule can easily be extended the higher vector/matrix dimensions, therefore, regarding Eq. (74), we find

(E,M,Σ)=(ML2T2,ML1T1,MT2) (76)

The inverse of the transformation is ensured, as easily checked, by

(L,M,T)=(E,M,Σ)Aˆ1Aˆ1=(1/201/21/213/2122) (77)

Within this context the search for the Π′s of the Buckingham method corresponds to the search of an exponential matrix transformation.

For example, the formalism applied to define the Π′s of Eq. (31) gives the transformation

[Π1,Π2,Π3]=[F,w,h,U,ρ,μ]LˆΠ1=Fw2U2ρ,Π2=hw,Π3=μwUρLˆ=(120210011000010111) (78)

Which according to the previously stipulated rules, yields6

[Π1,Π2,Π3]=[Fw2U2ρ,hw,μwUρ] (79)

It is worth stressing that this method does not add anything to DA and Buckingham, but it is simply a little detour, allowing a more advanced mathematical view to an apparently naïve problem. It is evident that the formal and conceptual apparatus associated with dimensional analysis (in physics and not only) are by no means elementary and that the relevant understanding requires multidisciplinary capabilities, ranging from applied Physics to quantum field theory. We have gone through different topics with the idea of conveying the underlying conceptual subtleties and the practical usefulness. This concluding section of the article is an attempt to raise further elements underscoring the two aspects of the relevant studies and what might be the possible future developments in this field.

### 4.2. Einstein’s Estimation of the Vibration Frequencies of Atoms

In order to better underscore what we mentioned at the end of the previous subsection, we like to resort an example paradigmatic for its importance and simplicity, regarding the use of DA by Einstein [15] to determine the characteristic frequency of vibration of atoms embedded in a crystal (see also [16,17]). If we make a list of the physical quantities involved into the relevant phenomenology, we note that it is determined by those defined in the following table

where MA is the mass of the atom, dA the inter-atomic distance, namely the distance between neighbor atoms, and χ is associated with the forces restoring the atom displacement from the original position. It is measured in terms of the coefficient of compressibility, having the dimensions of the inverse of the elasticity modulus.

The problem of expressing νA in terms of the remaining three independent variables follows from the principle of dimensional homogeneity and from the dimension matching techniques developed so far, the relevant conclusion amounts therefore to the equation

νA=CdAMAχ (81)

According to Einstein himself (as also reported in [16,17]) the difference between exact and DA treatment (of any physical problem) is a “mere” numerical factor, achievable through the analysis of the experimental data. In the case of this specific example Einstein found C = 1/3. No other Buckingham Π′s other than C are available since the total number of variables is 4 and the repeating counterparts 3. As already stressed, the proportionality constant C is not achievable using dimensional analysis only.

We can however make a step further towards the understanding of its value by introducing further variables allowing the comparison with tabulated physical constants (see [16,17]). The procedure we follow is that of including in our list of fundamental dimensions also the mol (the SI unit for the amount of a substance).

The introduction of this further degree of freedom allows the redefinition of atomic mass and of inter-atomic distance in terms of the molar mass Mm, molar volume vm, the Avogadro number NA and density ρ, which yield

MA=MmNA,dA=(vmNA)13=[MmNAρ]13 (82)

Thus rewriting Eq. (81) as

νA=CχNAMmρ3 (83)

What is new? A further row, including the MOL unit7 and a further column for the universal constant NA.

The corresponding dimension table is given below

MmρχNAνL03100M11100T00201MOL10010 (84)

Still we have only one Π which is now expressed as

Π=νAχMmρNA3 (85)

Using the experimental data for iron [16,17] given below

IronMmχρν55.81036101278606.75 (86)
plugged into Eq. (85) yields 0.33, consistent with the Einstein prediction.

We can proceed, even further, by including in our set of fundamental variables the Thermodynamic Temperature ϑ (unit 0K) the key element is the use of the melting temperature ϑs linked to χ by the relation

χ=RNAϑs (87)
where R is the Gas constant R = 8.31446261815324 JK−1Mol−1.

The associated dimensional matrix can therefore be written as reported below

ϑsρMaNARνL030020M011010T000021ϑ100010MOL001110 (88)

By repeating the same procedure as before and noting that the number of dimensionless parameters is 6 – 5 = 1, we end up with

Π=νϑ s1/2ρ1/3M a5/6N A1/3R1/2 (89)

The use of the experimental data (for Gold) reported below

GoldMmϑsρν197.21031347193102.79 (90)
yields C = 3.

The conclusion of this discussion is that going deeper with the dimensional analysis yields the possibility of entering deeply into the structure of the formulae expressing physical quantities. This is a further element of study to be pursued, which finds an echo in a statement attributed to Enrico Fermi [18]:

Those who really understand the nature of a given phenomenon must obtain the basic laws from dimensional considerations.

## 5. GENERAL RELATIVITY AND QUANTUM MECHANICS

### 5.1. Mercury Perihelion Precessions and the Pioneer Anomaly

We like now to comment on non-dimensional quantity often appearing in General Relativity (GR) and provided by GMRc2, (clight velocity) it rules crucial phenomena like

1. 1.

The deflection of light by a massive object

2. 2.

Gravitational redshift

3. 3.

The precession of Mercury perihelion

and also the global GPS (Global Position System) positioning.

We have already noted that (see Eq. (44))

[G]=[ρ1T2] (91)
which yields
[T2]=[GML3][D0]=[GMU2L] (92)
with U the velocity and D any physical dimension.

We have established that if we combine gravitational constant, mass velocity squared and length we get a dimensionless quantity.

Let us now consider the force exerted by a mass on a ray of light at grazing incidence on a spherical body of mass M and radius R (see Fig. 3). The effect of the deflection can be viewed as a kick induced by the gravitational force exerted by the body on the light photons. If we use for the “photon mass” the quantity mγ=ωc2 we obtain that the Energy variation of the photons before the deflection (Ei) and after the escape (E0) at very large distance (dR) is

ΔE=EiE0GMRc2ω (93)

It should be noted that the energy decreases after the deflection.

The associated fractional energy change ΔEE, directly linked to the deflection angle, is therefore

|ΔEE||Δpp|θGMRc2 (94)

We have no idea about the proportionality factor, the classical mechanics yields 2, GR requires the solution of the geodesic equations, not exactly an easy task [1922]. The conclusion of the computation is a factor 4.

It is evident that the deflection is proportional to mass; in our solar system the largest mass is the sun, using therefore the relevant parameters we find θ ≅ 2 · 10−6 rad. The experiment [20] measured 8 · 10−6 rad in agreement with the Einstein prediction.

The gravitational redshift [1921], namely the phenomenon of loose of energy of a photon moving through a gravitational field, is a naïve consequence of the previous discussion. We have indeed

ΔEE=ΔωωΔ(1λ)(1λ)GMRc2 (95)
which eventually yields
λ0λiλi=GMRc2λ0=(1+GMRc2)λi (96)
namely the escaped light has larger wavelength, hence lower frequency, as
ν0=(1+GMRc2)1νiν0(1GMRc2)νi (97)

Even though not immediately recognized, the previous relations have close relevance to the problem of Mercury perihelion precession, which owing to its closer proximity to the sun is more significantly affected by GR corrections. This means that a massive object modifies either the spatial and time structure of the space time metric

dt=(1GMc2R)dtdx=(1+GMc2R)dx (98)

The change of metrics determines a change in the GR Kepler equation, which determines an emerging new term correcting the Newtonian result [2023]. According to this correction, the revolution angle around the sun executed in one year by the planet is [1923]

Δϑ2π(1+3GMc2R) (99)

The advance per year of the perihelion can be seen in Fig. 4.

These tiny effects, ruled by the dimensionless constant GMc2R, have an impact on our daily life. In GPS the time measured by the clocks on the satellites runs faster than that measured at ground. In the case of satellites at distance from the earth 3RT the relative time variation with respect to ground is

ΔTTGMTRTc23RT4RT51010 (100)

This means that the absolute error accumulated in one day (86,400s) is ∆T ≅ 4.3 · 10−5 s, and the associated error (if not corrected) in measuring length (∆L = cT) would be about 14 km [19].

Before closing this topic, we like to mention the “heterodox” speculations recently published [24] in which GR predictions, like Mercury perihelion precession, have been ascribed to a structure of the space time viewed as a “dilatant-fluid” which, among other things, produces a kind of viscosity determining the deceleration of spaceships during their motion. The speculations in [24] have been motivated by the so-called Pioneer anomaly [25], consisting in a much shorter traveled distance with respect to that predicted by the computation. The deceleration inferred from the trajectory data has been estimated to be −8.74 · 10−10 m · s−2, a value which has been justified including in the equations of motion the effect of the recoil of the photons emitted via black-body radiation by the spacecraft. The alternative explanation put forward in [26] is that of assuming a drag force which has some resemblance with the already quoted Stokes force including the viscosity. As it is well known the Stokes force is written as

FD=6πr(vη) (101)
where r is the radius of the moving object, v its velocity and η the viscosity. According to the prescription of [24] the corresponding drag force of the dilatant-vacuum is restated by replacing the term in parenthesis on the right-hand side as
(vη)D=[γ1]κγ=11β2,β=vc (102)
where κ is a unit with the dimensions kg s−2, necessary to get the dimension of a force. The “new” force can eventually be written as
FD,S=6πrvηSηs=Dκv (103)
and ηs has the dimensions of P a · s. In the non-relativistic limit the above equation reads
FD,S3πr(κβ2) (104)

Putting this force in the classical equation for the planet motion and replacing the orbital velocity with

v2=GMaβ2=GMac2 (105)
we find that in one planetary period the angle is
Δϑ2π(1+6[1(GMac2)]1211e2)2π(1+31e2GMac2) (106)
with a, e being the orbit major axis and eccentricity respectively.

The results in Eq. (106) are equivalent to those in Eq. (99) and have been obtained within the context of classical mechanics, with special relativity corrections included in the dilatant part of the Stokes term. In both formulations the leading corrections are due to the dimensionless parameter GM/(Rc2).

Why have we reported this example?

Not to deny the GR to give credit to the (doubtful) introduction of an ad hoc force causing decelerating effects, rather to corroborate the importance of non-dimensional parameters which, once properly chosen, may provide a useful hint even though inserted in not a properly rigorous physical and mathematical context (see also our previous derivations of gravitational redshift and light deflection).

### 5.2. Quantum Mechanics and Dimensional Analysis

We have tried to convey the importance of thinking in terms of dimensions and what may be the impact on different aspects of Physical discipline. We must add a warning, which is absolutely necessary to be extremely cautious in using the tool we have described. The use of DA requires experience to guess an appropriate list of physical quantities. The wrong list would lead to meaningless results. An example, in absence of the Planck constant, the attempt of determining the size of atoms resulted in a complete failure [27,28].

The idea of introducing Planck’s constant as a fundamental quantity in atomic physics was disregarded for some time (a decade since the Planck paper on Black-Body radiation [29]). It had not yet the status of a fundamental constant, it being associated with thermal phenomena. It was, as noted in [28], a constant to be explained rather than to be used for explanations. It was later shown that charge mass and h can be combined to get a length (see below) but this observation did not provide any further progress, since nobody had indicated how the combination of this quantity had to come in the analysis of atomic structure.

The breakthrough came with the Bohr paper [30] and this is well known, but there were earlier suggestions and this story is worth to be mentioned. To the authors’ knowledge the first towards the search for a connection between the quantum of action and the atoms was Haas in his study of the stability of Thomson’s atom [31].

The reasoning of Haas was as follows. His main assumption was that one electron on the surface of an atom is attracted by the +Ne charge of the nucleus and repelled by the other remaining (N–1) electrons. If the charges are symmetrically distributed the electrons would experience the attraction of a single positive charge. If the centrifugal force is compared to the Coulomb attractive force one can get the speed of the electron around the nucleus, namely

The electron rotation period around the nucleus is therefore

T=2πr[14πɛ0e2mer]1/2 (108)

The next and final step by Haas was the assumption that the electron potential energy multiplied by the revolution time is just the Planck constant ħ (a quantity with the dimensions of an action, see below) namely

=2π[e24πɛ0mer]1/2 (109)
which, if inverted, yields the correct size of the atomic radius.

However, the intent of Haas was that of specifying the Planck constant (the mysterious quantity at the time) in terms of the atom parameters instead of the contrary. This attempt by Haas was not taken in due consideration and was strongly criticized for mixing up thermodynamics and mechanics. The importance of the idea was later recognized.

We can now proceed in terms of DA and see how the problem should be correctly afforded. The atomic size is dominated by forces of electromagnetic nature, we list accordingly physical quantities, r (Atomic Radius), h (Planck constant), ɛ0 (electric permittivity), me (electron mass) and e (electron charge). We are now obliged to enlarge our system of fundamental units, by including the charge whose dimensions are denoted by [Q]. We can therefore set

r=eam ebɛ 0chd (110)

Recalling that h has the dimensions of an action [ET] and [ε0] = [Q2F−1L−2], we have five physical quantities and four repeating variables and therefore we end up with (κ is a constant)

r=κɛ0h2mee2 (111)

Which introducing the fine structure constant

α=14πɛ0e2c=h2π (112)
can be written as
r=κ4παλeλeelectronComptonwavelength (113)

We have mentioned the constant α, a number which has created a great deal of speculation in the past. An account of this debate is outside of the scope of this paper; the interested reader is therefore referred to [32] where a survey is presented.

A final example we want to discuss is the derivation of the physical dimension of the magnetic charges (assuming that they exist) [33].

We remind that the force acting on electric and magnetic charges at rest in SI units are (μ0 is the vacuum magnetic permeability)

Fe=qeEFm=qmBμ0 (114)

Multiplying both sides of the previous equations, we get

FeFM=qeqmEBμ0 (115)

Recalling that Bμ0=HMagneticIntensity, from the dimensional point of view we find

[F2]=[qeqm][EH] (116)

Noting furthermore that [EH] has the same dimensions of the Poynting vector, namely Power/Surface, we end up with

[qeqm]=[F2L2][FLT1]=[FLT] (117)

From the dimensional point of view the product of an electric and a magnetic charge is an angular momentum. This would open further speculations on Dirac monopoles and charge quantization [33], which in unit of the reduced Planck constant can be written as

qeqm2π=n (118)

The relevant technical details go well beyond the scope of this article and will not be included in the present discussion.

In this paper, we have put together different examples from diverse branches of Physics. There are problems which we did not even touch, for example, regarding the nondimensionalization of differential equations describing a physical phenomenon.

The Navier-Stokes equation accounts for the momentum conservation of an incompressible fluid. It can be cast in the form

which is a transcription, in partial differential notation, of Eq. (22). The equation is written in dimensional form and the goal is that of writing it in terms of non-dimensional variables. It can be readily done and we find
Nstt*V*+(V**)V*=Neu*P*+1N F2g+1Re*2V* (120)
where we have defined
NstStroudhalNumber=fLVNeuEulerNumber=PP0ρV2 (121)
and
t*=ft,x*=xL,V*=V|V|g*=g|g|,*=L,P*=PPP0Pfcharacteristicfrequency (122)

The procedure of nondimensionalization has revealed two other quantities of crucial importance in the dynamics of fluids, namely the Stroudhal [34] and Euler [35] numbers.

Without entering into a discussion of the relevant importance, we note that the second represents the ratio between pressure and kinetic pressure. It is related to the cavitation number, which quantifies the potential of a flow to cavitate, namely the ability of a fluid to create bubbles, which after collapsing create a shock wave damaging machineries of a ship.

The Stroudhal number is important whenever non-stationary (oscillating) fluids are considered. It represents indeed a measure of the ratio of the inertial forces to those associated with the unsteadiness of the flow. If it turns out to be a small quantity with respect to the others entering the normalized form of the Navier-Stokes equation, the problem can be considered stationary.

Before concluding, we like to go back the incipit of the article where we have mentioned the concept of dimension in Euclidean geometry, in which we have defined dimensions of integer order (0, 1, 2, 3). The Fractal geometry admits non-integer dimensions. This has opened new avenues in the understanding of the interplay between (physical/geometrical) dimensions and scaling laws. This aspect of the problem will be the topic of a forthcoming investigation.

We hope that we have provided at least two feelings:

• The problem of physical dimensions is an issue of fundamental importance in Physics, rich of conceptual aspects and providing unsuspected links between apparently uncorrelated fields.

• Our paper, albeit long, has scratched the surface of a huge field of research.

Regarding this last point, we have ignored important aspects of the discussion like those associated with the renormalization group, DA in quantum field theory and dimensional collapse. These topics deserve a deeper understanding of the problems and are worth to be treated in a separate discussion.

## Footnotes

U instead of V to avoid confusion with volume.

We remind that the theorem asserts that the variation of the kinetic energy of a body is equivalent to the work done on that body.

This choice is just matter of convenience.

Being min(rank(A),m) the minimum rank of the A matrix written in Eq. (28).

The differential form of the Gauss’s law for gravity assumes the expression ∇2ϕ = 4πGρ, being ϕ a scalar potential.

The matrix Lˆ is not square and does not have an inverse, therefore the first of Eq. (75) is apparently not invertible. The problem can be overcome and the interested reader can find an appropriate discussion in [10].

The MOL in the SI denotes the unit measuring the amount of a substance.

## REFERENCES

R Butterfield, Dimensional Analysis for Geotechnical Engineers, Journal of Geotechnique, Vol. Vol. 49, No. 3, 1999, pp. 357-366.
PW Bridgman, Dimensional Analysis, Yale University Press, New Haven, 1931.
E Buckingham, Model Experiments and the Forms of Empirical Equations, American Society of Mechanical Engineers, Transactions, Vol. Vol. 37, 1915, pp. 263-292.
A Olivieri, F Pistani, R Wilson, L Benedetti, F La Gala, EF Campana, and F Stern, Froude Number and Scale Effects and Froude Number 0.35 Wave Elevations and Mean-Velocity Measurements for Bow and Shoulder Wave Breaking of Surface Combatant DTMB 5415. IIHR Report n. 441, Hydroscience & Engineering College of Engineering, The University of Iowa, Iowa City, Iowa 52242-1585, USA, 2004.
CB Barrass, Ship Design and Performance for Masters and Mates, Elsevier, Burlington (MA), 2004.
MJ Moran, HN Shapiro, BR Munson, and DP DeWitt, Introduction to Thermal Systems Engineering: Thermodynamics, Fluid Mechanics, and Heat Transfer, John Wiley & Sons, New York (NY), 2003.
J Dunkel, Lecture notes, MIT edu., 2014.
H Hanche-Olsen, Buckingham’s Pi-Theorem.
CG Callan, Broken Scale Invariance in Scalar Field Theory, Physical Review D, Vol. Vol. 2, 1970, pp. 1541.
A Einstein, Elementare Betrachtungen über die Thermische Molekularbewegung in Festen Körpern, Annalen der Physic, Vol. Vol. 35, 1911, pp. 679-694. Deutschland
B Carneiro, Análise Dimensional e Teoria da Semelhança e dos Modelos Físicos, Editoria UFRJ, Rio de Janeiro, 1996, pp. 258.
E Fermi, The Power of Dimensional Thinking, Quantum, Vol. Vol. 2, No. 34, 1992.
M Berry, Principles of Cosmology and Gravitation, Cambridge University Press, London, 1976.
G Dattoli, M Migliorati, and K Zhukovsky, An Elementary Account of Relativistic Cosmology, La Rivista del Nuovo Cimento, Vol. Vol. 10, 2006, pp. 1-85.
S Weinberg, Gravitation and Cosmology: Principles and Applications of The General Theory of Relativity, Wiley & Sons, Inc., New York, 1972, pp. 188.
R Gilmore, Dimensional Analysis, Physics Department, Drexel University, Philadelphia (PA 19104), 2004.
P Stehle, Order, Chaos, Order: The Transition from Classical to Quantum Physics, Oxford University Press, New York, Oxford, 1994.
J Mehra and H Rechenberg, The Historical Development of Quantum Theory. Vol. I, Part 1, Springer-Verlag, New York, 1982.
N Bohr, On the Constitution of Atoms and Molecules, Philosophical Magazine. Series 6., Vol. Vol. 26, 1913, pp. 857-875.
AH Haas, Über eine neue Theoretische Methode zur Bestimmung des elektrischen Elementarquantums des Wasserstoffatoms, Physikallische Zeitschrift, Vol. Vol. 11, 1910, pp. 510.
G Dattoli, The Fine Structure Constant and Numerical Alchemy, 2010. arXiv:1009.1711 [physics.gen-ph].
JD Jackson, Classical Electrodynamics, 3rd edition, Wiley & Sons, Inc., New York-London, 1963.
FM White, Fluid Mechanics, 4th edition, McGraw Hill, New York, 1999. (ISBN 978-0-07-116848-9).
GK Batchelor, An Introduction to Fluid Dynamics, Cambridge University Press, Cambridge, 1967. (ISBN 0-521-09817-3).

ris
TY  - JOUR
AU  - G. Dattoli
AU  - E. Di Palma
AU  - E. Sabia
PY  - 2022
DA  - 2022/05/06
TI  - Is There Anything Interesting to Say About Dimensions in Physics?
JO  - Growth and Form
SP  - 3
EP  - 25
VL  - 3
IS  - 1-2
SN  - 2589-8426
UR  - https://doi.org/10.55060/j.gandf.220605.001
DO  - https://doi.org/10.55060/j.gandf.220605.001
ID  - Dattoli2022
ER  -

enw
bib