## Services on Demand

## Journal

## Article

## Indicators

- Cited by SciELO
- Access statistics

## Related links

- Cited by Google
- Similars in SciELO
- Similars in Google

## Share

## Revista Brasileira de Ensino de Física

##
*Print version* ISSN 1806-1117

### Rev. Bras. Ensino Fís. vol.33 no.4 São Paulo Oct./Dec. 2011

#### http://dx.doi.org/10.1590/S1806-11172011000400006

**ARTIGOS GERAIS**

**Wave propagation in a non-uniform string**

**Propagação de ondas em uma corda não-uniforme**

**Fernando Fuzinatto Dall'Agnol ^{1}**

Centro de Tecnologia da Informação Renato Archer, Campinas, SP, Brazil

**ABSTRACT**

In this article I present the behavior of the reflection and transmission coefficients of a pulse at a joint between two strings with mass densities *μ*_{1}and *μ*_{2}. The joint is made of a string segment with mass density varying linearly from *μ*_{1} to *μ*_{2}. It will be shown that the reflection of the pulse at the joint depends largely on the ratio between the pulse width and the length of the joint. Analogies with other physical systems such as antireflection coatings and tsunamis will be considered briefly.

**Keywords:** wave, pulse, string, non-uniform, inhomogeneous, discontinuous, anti-reflection, anti-reflective, tsunamis, pulse energy.

**RESUMO**

Neste artigo mostro o comportamento dos coeficientes de reflexão e transmissão de um pulso propagando através de uma emenda entre cordas com densidades de massa *μ*_{1} e *μ*_{2}. A emenda é um segmento de corda com densidade de massa variando linearmente desde *μ*_{1} a Será mostrado que a reflexão do pulso na emenda depende sensivelmente da razão entre a largura do pulso e o comprimento da emenda. Discutirei brevemente analogias com outros sistemas físicos como camadas anti-refletoras e tsunamis.

**Palavras-chave:** onda, pulso, corda, não-uniforme, não-homogêneo, descontínuo, anti-refletor, tsunamis, energia de pulso.

**1. Introduction**

The basic concepts of wave mechanics in physics courses and textbooks for undergraduate students are usually illustrated with sine functions of a vibrating string. In these functions, properties such as standing waves, wavelength, period/frequency, phase and amplitude are well defined and they can be easily calculated given the tension and the mass density of the string. The relevance of these introductory chapters on wave mechanics is its applicability in many areas of physics like acoustic, optics, electromagnetism and quantum mechanics. In the textbooks and in the internet one can find various animations of waves in uniform strings [1-5] and in strings with different mass densities [6] joined together. There are also various articles describing the behavior of waves in non-uniform strings, e.g. the effect of non-uniformity on natural tones in musical instruments [7], in quantum waves [8,9] and seismological wave propagation inside the Earth [10]. The more sophisticated the physical model, the more phenomena one can describe and explain. In this work I introduce a specific non-uniformity in the string, being a sophistication that is not covered in the introductory texts, but it allows understanding the working principle of the anti-reflection coatings.

In this article, the non-uniformity consists of two semi-infinite strings with mass densities *μ*_{1} and *μ*_{2} joined with a string segment with mass density varying linearly from *μ*_{1} to *μ*_{2}. The objective is to show the behavior of the reflection of a pulse passing through the string joint. To simulate a pulse is preferred instead of a sinusoidal travelling wave because the former is a limited wave package and it is easier to visualize and compare the amplitudes or the pulse width variations as the reflection drops.

**2. The wave equation**

The propagation of a wave in a string is described by the wave equation

where y(x,t) is the transversal displacement of the wave at position *x* and time *t* and *v* is the velocity in the *x* direction or *group velocity,* to be distinguished from the transversal velocity of the string *¶y*/*¶t.* The *v* depends on the tension and mass density of the string according to

where *B* is the tension in the string, expressed in N, and *μ* is the linear mass density, expressed in kg/m.

**3. Physical system**

Figure 1 shows a representation of non-uniform string fixed at two boundaries far from the joint. Since our analysis will be restricted to the central portion of the string the pulse interaction with the boundaries will be neglected; so, in the pictures hereafter these boundaries will not be indicated. A pulse with a Gaussian profile having width D*x* and initial amplitude *y*_{0} propagates to the right. In *x* = 0 there is a joint segment of length D*L,* where the mass density of the string varies linearly from *μ*_{1} to *μ*_{2}. Later in this article it will be important to consider that the string is inextensible so the pulse is formed by mass accumulation (not by strain) and the tension of the string is provided by a hanging weight (not by elastic forces).

After interacting with the joint the pulse will be partially reflected and partially transmitted. Here I'll present a numerical analysis of the reflection *R* coefficient as a function of D*L* and D*x*. *R* is defined by the ratio between the energy of the reflected pulse and the energy of the incident pulse

The initial wave with the characteristics shown in Fig. 1 can be written as a standard Gaussian function

where *x*_{0} is the peak position and Ax is the width of the Gaussian. The mass density of the string varies according to the function

For the description of *R* it is efficient first to consider the energy of the pulse. Figure 2 shows an example of a pulse propagating to the right in a referential system Q. The pulse interacts with the joint and is partially reflected. The frames drawn enclosing the incident, the reflected and the transmitted pulses are the local referential system named *P _{i} (i* = 0, 1 or 2) that moves together with the pulse. This referential coordinate system facilitates the evaluation of the energy of the pulse as was demonstrated by Juenker [11]. In this analysis it is important that the string in inextensible, so the pulse is made by the action of mass only and not by stretching the string. In the system

*P*, the pulse doesn't move; however, the string is seen moving from the right to the left along the pulse like a train in a Gaussian shaped railroad with local velocity

*v*,which will be denoted by

*v*

_{1}in the string with mass density

*μ*

_{1}and

*v*

_{2}for

*μ*

_{2}(see Fig. 3).

The energy of the propagating pulse can easily be obtained from Juenker's procedure as follows: In *P* any mass element, *dm,* has a velocity *v* parallel to the string. From *P* one can conclude that the velocity, v*, of *dm* in *Q* is the vector sum of its velocity in *P* plus the velocity of *P* in *Q* as indicated in Fig. 3. Using the law of cosines, *v** results in

The kinetic energy of a mass element is

where v* was replaced by the expression in Eq. (5). Note from Fig. 3 that cos*θ* = *dx/ds* so the energy of a mass element is given by

The total energy of the pulse is the integral of Eq. (7)

where D*z* = *x _{b}* -

*x*is the length of the frame that completely encloses the pulse and D

_{a}*s*is the integral in

*ds*which gives the total length of the Gaussian curve in the interval D

*z*. Usually Eqs. (7) and (8) cannot be solved analytically when the pulse interacts with a non-uniform segment. Notice that to evaluate the energy in Eq. (8) it doesn't matter if D

*z*is larger than necessary to enclose the pulse, because the difference

*Ds —*D

*z*will tend to a constant. D

*s —*D

*z*is the string accumulated to create the pulse.

The initial, reflected and transmitted pulse energies, *E*_{0}, *E*_{1} and *E*_{2} respectively, are

The *R* and *T* are given by the ratio

**4. Numerical analysis**

In this section I explain the numerical evaluation of the wave function. The propagation of the wave is calculated using finite differences in the time domain (FDTD). The string is divided in a large number of cells with length *dx*. The solution at a time *t* is evaluated from the solutions of the previous two time intervals *t—dt* and *t—2dt* as shown in Fig. 4.

The differential Eq. (1) is written as an equation of differences [12]

where *i* corresponds to the index of the *i ^{th}* node of the string,

*j*is the

*j*time interval and

^{th}The condition 0 *< r* £1 is necessary for Eq. (14) to converge [12]. The *y*(*i*,*j*) has a direct relation to *y*(*x*,*t*) since each ( *x,t)* has a corresponding cell *(i,j)* according to *x* = *x _{ini}+idx* and

*t*=

*jdt.*Initial conditions are

and the boundary conditions are fixed rims

The energies (9) to (10) are obtained numerically from the expression

where, *i*_{1} and *i*_{2} are the first and last cells of the frame that contains the pulse, *μ*(*ί*) is the position depending mass density. Eq. (18) is not applicable if the pulse is interacting with the joint, because the shape of the pulse is changing. Then, Eq. (5) is not valid because there is one more velocity component that depends on the position in a complicated way. To use Eq. (18) the reflected and transmitted pulses must be far from the rims and from the joint, so the numerical propagation must be carried out until this condition has been met.

**5. Results**

Figure 5 compares two cases in which a pulse passes through a discontinuous joint and a smooth joint segment. It can be seen that the amplitude of the reflected pulse in the smooth joint is smaller and wider, which implies that its energy is lower. However, the transmitted wave is taller, so the total energy is conserved.

Figure 6(a) and (b) show the behavior of *R* in which a pulse is passing the joint. Figure 6 (a) shows *R* for *μ*_{1} = 1 kg/m to *μ _{2}* = 4, 9 and 16 kg/m. Figure 6(b) shows the other way around in which the pulse passes from high density

*μ*1 = 4, 9 and 16 kg/m to low density

*μ*

_{2}= 1 kg/m. In all cases the reflection decreases as the length of the joint increases. It can be seen that for D

*L*= 0 the reflection coefficients in both cases are equal. Figure 7 compares

*R*for two pulse widths. The smaller the width, the smaller is the reflection. In the limit of D

*x*<< D

*L*then

*R —*0 and if D

*x*>> D

*L*then

*R —*constant independent of D

*x*or D

*L*.

The transmission coefficient is simply *T* =1-R, so the curves for *T* are not shown here.

**6. Analogies with other physical systems**

As any analogy this analysis only provides insights to some extent to other physical systems, after which it fails due to their particularities. Nevertheless, this analysis is worthwhile, since it helps comprehending qualitatively the transmission of light in anti-reflection coatings and tsunamis. I'll describe these examples briefly.

**6.1.** **Anti-reflection coatings**

The refraction index variation between two media always causes some light reflection at the interface. It is tempting to make an analogy between the index of refraction and wavelength in optics and the mass density and the pulse width respectively in non-uniform strings. Some anti-reflection coatings deal with a smooth gradient of the refractive index and these coatings or transitions show indeed very good anti-reflection behavior. An interesting example is the moth eye, in which anti-reflection pyramidal nanostructures have grown at the interface. These structures are smaller than the wavelength of visible light, so the light passes as if the pyramids are a continuous layer of varying index [13].

**6.2.** **Tsunamis**

The Earth crust accommodation can cause the seabed to lift, lower or displace to the sides. In any case, a large water column is displaced. Also in this case it is tempting to make a comparison between the depth of the water with the mass density and the water displacement with the pulse on a string. Then the decreasing depth of the seabed is analogous to a string joint with decreasing mass density. When the pulse is transmitted to a low density string its amplitude increases. The same is expected in water, which explains why a small amplitude at the deep see becomes a tsunami at the beach. In this case however, the analogy is weak because it fails in two other noticeable characteristics: the pulse in the string is faster and wider in the lighter string whereas tsunami waves become slower and narrower in shallow waters. Water waves have different boundary conditions and involve mass conservation; so, for the amplitude to increase, more water must be taken from the front, rear and beneath the wave and Eq. (1) does not account for this. It is interesting to see the simulation of a tsunami in Ref. [14] and to compare this with the pulse in a string in the simulation in Ref. [15].

**7. Conclusion**

The simulation of pulses in a non-uniform string leads to the following conclusions:

The reflection coefficient is reduced when the transition between two strings is smooth in terms of mass density. An equivalent phenomenon is observed at the interface of two transparent media where the optical reflection also decreases largely upon a smooth transition of the refractive index.

The reflection coefficient depends on the ratio of the pulse width and the length of the joint segment. The analogy in anti-reflective coatings is that these are designed to be efficient in a certain wavelength band, whereas for large wavelengths the coating is too short to diminish the reflection appreciably.

**Acknowledgement**

I'm thankful Prof. Daniel den Engelsen for his help in preparing the manuscript.

**References**

[1] A. Kandus, F.W. Gutmann and C.M.C. Castilho, Revista Brasileira de Ensino de Física 28, 427 (2006). [ Links ]

[2] D. Gámez e R. Welti, Revista Brasileira de Ensino de Física 26, 139 (2004). [ Links ]

[3] M.G. Souza, T.R. Cicogna e A.J. Chiquito, Revista Brasileira de Ensino de Física 29, 5 (2007). [ Links ]

[4] S.L. Piubelli, H.A. Errobidart, S.T. Gobara e N.C.G. Errobidart, Revista Brasileira de Ensino de Física 32, 1501 (2010). [ Links ]

[5] Typical simulation in http://www.youtube.com/watch?v=Rbuhdo0AZDU&feature=related. [ Links ]

[6] Typical simulations in http://www.youtube.com/watch?v=Zn9loyHbW4U and http://www.youtube.com/ watch?v=0qToBvL7q_g. [ Links ]

[7] P. Chumnantas, C.A. Greated and R. Parks, Journal de Physique IV 4, C5-649 (1994). [ Links ]

[8] O.C. Castellani, M.G. Goncalves e A.J. Santiago, Revista Brasileira de Ensino de Física 15, 34 (1993). [ Links ]

[9] M.A. Cândido Ribeiro, V.C. Franzoni, W.R. Passos, E.C. Silva e A.N.F. Aleixo, Revista Brasileira de Ensino de Física 26, 1 (2004). [ Links ]

[10] http://www.youtube.com/watch?v=yOGoKCK17a4{&}NR=1. [ Links ]

[11] D.W. Juenker, Am. J. Phys. 44, 94 (1976). [ Links ]

[12] M.N.O. Sadiku, *Monte Carlo Methods for Electromagnetics* (CRC Press, Boca Raton, 2009). [ Links ]

[13] http://en.wikipedia.org/wiki/Anti-reflective_coating. [ Links ]

[14] http://www.seed.slb.com/flash/science/ features/earth/livingplanet/tsunami_wave/ en/index.htm?width=770&height=370&popup=true. [ Links ]

[15] http://www.youtube.com/watch?v=yAMbu-8QB_Y . [ Links ]

Recebido em 17/2/2011

Aceito em 13/6/2011

Publicado em 29/11/2011

1 E-mail: fernando.dallagnol@cti.gov.br.