Abstract
The scattering of waves in an inhomogeneous infinite string with one change in the density of the string is well known. In this paper, we study the case where there are two discontinuities in the density of the string. It turns out that we can write the solution for this case as well. The form of the solution will be given in finite sums of reflected and transmitted waves over finite time. The finite sums become converging infinite series for infinite time.
1. Introduction
The scattering of waves in an infinite string due to a single change in the density of the string has been studied in graduate text books such as [1, 2]. In [1] the author notes that in practice a scattering of this sort occurs in physical systems such as submarine cables or telephone lines when a join in the system creates a transmitted as well as a reflected wave. Both of these waves can be computed in terms of the incoming signal. For practical purposes it is desirable to suppress the reflected waves by attaching, at the join, a damping mechanism or a point mass so that the reflected waves do not interfere with the incoming signals. All three cases of no suppression, suppression by damping, and suppression by point mass have been studied in [1] by assigning appropriate boundary conditions at the join. In this paper we study the case where the string has more than one change in its density and no suppression occurs at the joins. We write the form of the solution when there are only two changes in the density. We also show what the reflected and transmitted waves due to an incident wave on both sides of a join at an arbitrary point in the string look like. This will allow us to theoretically write the form of the reflected and transmitted waves for higher number of changes. But because of the repetitive bouncing back and forth of the waves between each interface writing a general formula for the solution of the wave problem with more than two discontinuities will be impossible. Just imagine that all the waves inside one pair of joins get transmitted to the neighboring pair and in turn bounce against the joins and create waves of their own. Writing a formula to account for all of these waves will require a horrendous amount of work if at all possible, and we are not going to tackle it here.
To fix the ideas, consider two semi-infinite strings with densities and , , joined together at the origin . The speed of the waves traveling along the two pieces is and , respectively. An incoming wave , such that , will be scattered to partly reflected and partly outgoing at . Suppose the resulting wave satisfies that where Under the geometric continuity (the string is continuous at ) and dynamical continuity (the transverse force is continuous at ) conditions, and the initial conditions, the solution is well known to be (see [1, 2]) Now, consider the above setup extended to a string made of three pieces in the intervals , , and for some constant . Denote the densities and their corresponding wave speeds by , , and , , , respectively. We will study the behavior of an incoming wave , where , , at the interfaces and show how, under similar conditions as above, similar formulas for the scattered waves can be written. We note that each outgoing wave from an interface will be incident on the next, and the reflected waves at the next interface will be incident on the previous. Each incident wave will, in turn, be scattered into reflected and outgoing waves at the interfaces. We will write the solution to the wave problem as sums of these scattered waves.
2. The Problem
Consider the problem Assume that a wave , where , , is incoming. At the wave scatters as follows: The transmitted wave becomes incident at the point . Let us denote We will show here, for the sake of completeness, that also scatters at and creates new waves.
Theorem 1. The incoming wave given by (2.3) scatters at , for , in the form
Proof. First note that when , and , , by the definition of in (2.3). Consider the wave function , which we assume as on . The solution to the wave equation where is of the form We use the geometric and dynamical continuity conditions, and the initial conditions, to find , , , and . From continuity of the wave function at , we have Using dynamical condition (2.8) at and differentiating equation (2.10) with respect to , we have the system Solving the system (2.11) for and and integrating, we have From initial conditions (2.9) at , The right-hand sides of (2.14) provide a system that can be solved for and : At the time and , we have Differentiating the right-hand side of (2.17) and solving the system obtained from (2.17) and (2.18), namely, we will have By (2.15), when the arguments of and are negative the functions are equal. This means that Notice that by (2.21) if the argument of satisfies , then . For the argument of in (2.12) and (2.13) is larger than . Therefore, and (2.12) and (2.13) become, From (2.22) and (2.25) we obtain Equation (2.26) determines as follows: Denote the argument of in (2.27) by , then Due to (2.28), Finally, putting (2.23) and (2.27) in (2.12) determines : A similar change of independent variable in (2.30) yields Plugging , , , and determined by (2.22), (2.23), (2.29), and (2.31), respectively, in (2.7) yields in (2.4).
By an argument similar to the one in Theorem 1, we can show the following. The reflected wave at , namely in (2.4), becomes incident at the point . Let us denote Then, the wave scatters at , for , as follows: We note here that is the wave that moves to the left with , . It is not difficult to check that satisfies and .
In this manner, the forms of outgoing waves through the interfaces and and the ones bouncing back and forth between the two are determined. In order to write the solution to the problem where, subject to for , , , positive constants, and the incoming wave for some , we need sums of the waves in each interval , , and . One way to write such solution is through the use of composition of the arguments of the scattered waves. For this purpose we introduce the following functions: Let us write the functions and introduced in (2.3) and (2.32), respectively, in terms of the variable : From in (2.37) and (2.38)-(2.39), can be written in terms of as follows: In terms of the functions in (2.37), (2.38) and (2.40), (2.4) and (2.33) in the time intervals and will, respectively, become Now, we begin to write the solution of the wave problem (2.34)-(2.36) in each interval , , and using (2.2), (2.41), and (2.42). Consider the following interpretations of (2.2), (2.41), and (2.42). First we have the following wave on the interval : When this wave hits the join at 0 from the left, by (2.2) the reflected wave gets a coefficient of the form and an argument change to , that is, the reflected wave corresponding to (2.43) is is as given in (2.37). The transmitted wave due to (2.43) on the other hand gets a coefficient of and argument change to , that is, the transmitted wave traveling to the right at 0 is Now let us see what changes the signal (2.45) undergo when it hits the join at on the left. From (2.41) part of the signal gets reflected to the left acquiring the coefficient and argument change to . The transmitted wave traveling right acquires the coefficient and the argument change . Lastly, by (2.42) a signal that hits the join at 0 on the right its transmitted part acquires the coefficient and the change of argument . Its reflected part acquires the coefficient and the argument change of .
Note that the signals inside the interval keep scattering, but the ones that are transmitted outside this interval continue traveling to the right or left forever. So, now let us look at a few more wave signals that are produced inside . Look at the wave in (2.42). As it travels to the right it reaches the join as an incident wave on the left. By the argument above its reflected component picks up a coefficient of and an argument change of . The transmitted component picks up a coefficient of and an argument change . These two wave are, respectively, Now the wave (2.46) is incident on the join 0 from the right. The argument above shows that the transmitted part must pick up the coefficient and the change of argument , and the reflected part picks up the coefficient and the argument change of . These two new waves are given below, respectively, The waves (2.47) and (2.48) will be moving to the right and left away from the interval indefinitely. But the wave (2.49) will be incident on the join on the left, and the process of reflection and transmission repeats as before. It is not too difficult now to decipher the general pattern of the reflected and transmitted waves in and the outgoing waves outside of this interval. Here they are. The waves outgoing in are of the form The waves in moving left are The waves in moving right are Finally the ones moving right for are Now, we recall that given in (2.2) is only good for the time interval , and after , the transmitted part hits the join at and splits. Therefore, in this time interval we have the wave defined as We note that is the original incoming signal incident on 0 plus the reflected and transmitted waves. The transmitted wave has not reached the join yet, and so there is no wave beyond the point . Therefore in (2.54) satisfies the wave problem (2.34)–(2.36) in this time interval. We now bring in another wave in the interval . We define it by The function in (2.55) represents the reflected and transmitted waves at , caused by the incident signal . We see this incident wave in the definition of in (2.54). We argue that the sum of the wave in (2.54) and (2.55), , over the interval satisfies the wave problem (2.34)–(2.36). The fact that they satisfy the wave (2.34) over each space interval is clear. The boundary conditions at 0 follow from the fact that is the solution to the wave problem (1.1)–(1.4) and the function = for . The boundary conditions at are satisfied, because the sum of the parts for and is the the solution of the wave problem (2.5), (2.8)–(2.9), with given by (2.3) there. The initial conditions are satisfied because satisfies them, and due to the argument of staying positive when .
Now we consider wave over the interval . It is defined by This is the wave due to the scattering of , at the join 0. This incident wave can be seen in the definition of in (2.55). We claim that the sum , where are given by (2.54)–(2.56), satisfies the wave problem (2.34)–(2.36) over the interval . We note that when added to the incident wave is the wave that was given in (2.33) except for the notation. This can be seen by noting that in (2.33) is defined in terms of in (2.32) and is given in terms of in (2.3). We constructed so that it satisfies the boundary conditions in (2.36), as well as the wave equation (2.34) for . The sum contains for . In light of the previous argument about and the way contributes to the sum, it is not difficult to see why satisfies the boundary conditions in (2.36) and the wave equation in (2.34) in the interval . The fact that at time all terms involving except are zero in shows that the initial conditions of (2.30) are also satisfied.
If one waits other units of time, another scattering happens at from the wave and then again at 0 from the reflected wave of . This process continues forever, and new waves appear every units of time. The functions can be defined as before such that their sum will satisfy the wave problem (2.34)–(2.36) over longer and longer time intervals. In the statement below we write the solution in the form of a finite series whose upper limit depends on the number of elapsed. In doing so we use the general forms of the waves in (2.50)–(2.53).
We summarize these results in the following theorem.
Theorem 2. Let , , be an incoming wave. Then the solution to the problem (2.34)–(2.36) in the time interval , is given by Furthermore, assume that and its derivatives of up to the order two are bounded in . Then, the solution to (2.34)–(2.36) over the time interval is the limit as of the above finite time solution.
Proof. Note that for as described in (2.54), for , where is given in (2.55), and for , where is defined in (2.56). We showed that in all these cases these sums satisfy the wave problem (2.34)–(2.36). For , we will be adding reflected and transmitted waves at each interface 0 or which, when added to their corresponding incident waves, satisfy the boundary conditions in (2.36). We also argued that the initial conditions are satisfied because for close to 0 all terms involving except will be zero. So solves the wave problem in .
It remains to show the limiting case when . For that, let us just focus on one piece of the function in the interval :
First notice that by definitions in (2.37) the argument of in the summation (2.58) can be simplified to
The coefficients of satisfy . Let
then . Now with being bounded, the partial sums of the series (2.58) are bounded above by the partial sums of a convergent geometric series in powers of . Therefore it is absolutely and uniformly convergent. On the other hand, the term-by-term differentiation of (2.58) in terms of results in an extra coefficient and no extra coefficient in terms of , as can be seen by (2.59). Since and are also bounded, the resulting series also converge absolutely and uniformly to the derivatives of the limit of the series (2.58). A similar argument can be applied to the series for in the intervals and .
We have shown the form of the scattering of an incident wave on the right- and left-hand side of the join 0. But we have only shown the scattering of an incident wave on the left of . If one wanted to study a higher number of interfaces, then it would be necessary to know what happens when an incoming wave hits an interface at an arbitrary point on the right. Here, and for the completeness of the argument, we point out that the transmitted and reflected waves at from an incoming wave from the right would be as follows: Theoretically, this should enable us to extend the result of Theorem 2 to a higher number of discontinuities in the density of the string. The difficulty would be the ability to keep track of all incoming and outgoing waves as well as the ones that bounce back and forth between the interfaces. Since there are so many such waves, with even three interfaces, a general solution for more than two interfaces is impossible to write down. However, for a specific, relatively short length of time the solution can be found, when the discontinuities are few, say three. In this case the form of the solution will also depend on the distance between the discontinuities. Then, one can ask up to how many joins or for what length of time will the scattered waves be tractable. Another interesting question would lie in the area of inverse scattering. Since by our experience the solution contains the location of a discontinuity, in the absence of such knowledge will it be possible to find the coordinate of the interface from the form of the scattered waves.