
{"id":391,"date":"2016-11-03T10:47:03","date_gmt":"2016-11-03T10:47:03","guid":{"rendered":"https:\/\/jhubisz.expressions.syr.edu\/phy360\/?p=391"},"modified":"2016-11-04T01:49:49","modified_gmt":"2016-11-04T01:49:49","slug":"week-9","status":"publish","type":"post","link":"https:\/\/jhubisz.expressions.syr.edu\/phy360\/2016\/11\/03\/week-9\/","title":{"rendered":"Week 9"},"content":{"rendered":"<p>During this week, we derived the wave equation by taking the continuum limit of the discrete system (where there were a finite number of masses attached to a string). \u00a0We started from the equation of motion for a single mass, which we derived last week:<\/p>\n<img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cddot%7By%7D_p+%3D+-+%5Cfrac%7BT%7D%7Bl+m%7D+%5Cleft%5B+2y_p+-+y_%7Bp-1%7D+-+y_%7Bp%2B1%7D+%5Cright%5D&#038;bg=T&#038;fg=000000&#038;s=0' alt='\\ddot{y}_p = - \\frac{T}{l m} \\left[ 2y_p - y_{p-1} - y_{p+1} \\right]' title='\\ddot{y}_p = - \\frac{T}{l m} \\left[ 2y_p - y_{p-1} - y_{p+1} \\right]' class='latex' \/>\n<p>which we can we write as<\/p>\n<img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cddot%7By%7D_p+%3D+-+%5Cfrac%7BT+l%7D%7Bm%7D+%5Cfrac%7B2+y_p+-+y_%7Bp-1%7D-y_%7Bp%2B1%7D%7D%7Bl%5E2%7D&#038;bg=T&#038;fg=000000&#038;s=0' alt='\\ddot{y}_p = - \\frac{T l}{m} \\frac{2 y_p - y_{p-1}-y_{p+1}}{l^2}' title='\\ddot{y}_p = - \\frac{T l}{m} \\frac{2 y_p - y_{p-1}-y_{p+1}}{l^2}' class='latex' \/>\n<p>then we use the fact that the second derivative of a function can be approximated by<\/p>\n<img src='https:\/\/s0.wp.com\/latex.php?latex=f%27%27%28x%29+%5Capprox+%5Cfrac%7Bf%28x%2B%5Cepsilon%29+%2B+f%28x-%5Cepsilon%29+-+2+f%28x%29%7D%7B%5Cepsilon%5E2%7D&#038;bg=T&#038;fg=000000&#038;s=0' alt='f&#039;&#039;(x) \\approx \\frac{f(x+\\epsilon) + f(x-\\epsilon) - 2 f(x)}{\\epsilon^2}' title='f&#039;&#039;(x) \\approx \\frac{f(x+\\epsilon) + f(x-\\epsilon) - 2 f(x)}{\\epsilon^2}' class='latex' \/>\n<p>this allows us to re-write the equation of motion as<\/p>\n<img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cfrac%7B%5Cpartial%5E2+y%7D%7B%5Cpartial+t%5E2%7D+%3D+%5Cfrac%7BT%7D%7B%5Cmu%7D+%5Cfrac%7B%5Cpartial%5E2+y%7D%7B%5Cpartial+x%5E2%7D&#038;bg=T&#038;fg=000000&#038;s=0' alt='\\frac{\\partial^2 y}{\\partial t^2} = \\frac{T}{\\mu} \\frac{\\partial^2 y}{\\partial x^2}' title='\\frac{\\partial^2 y}{\\partial t^2} = \\frac{T}{\\mu} \\frac{\\partial^2 y}{\\partial x^2}' class='latex' \/>\n<p>or, in terms of a velocity:<\/p>\n<img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cfrac%7B%5Cpartial%5E2+y%7D%7B%5Cpartial+x%5E2%7D+%3D+%5Cfrac%7B1%7D%7Bv%5E2%7D+%5Cfrac%7B%5Cpartial%5E2+y%7D%7B%5Cpartial+t%5E2%7D&#038;bg=T&#038;fg=000000&#038;s=0' alt='\\frac{\\partial^2 y}{\\partial x^2} = \\frac{1}{v^2} \\frac{\\partial^2 y}{\\partial t^2}' title='\\frac{\\partial^2 y}{\\partial x^2} = \\frac{1}{v^2} \\frac{\\partial^2 y}{\\partial t^2}' class='latex' \/>\n<p>So we now have the\u00a0<strong>wave equation!<\/strong> \u00a0This is our new favorite equation, and it is important to note that we obtained it from the behavior of a large number of coupled simple harmonic oscillators &#8211; that means that there is nothing intrinsically new here, but we have found a convenient formalism for dealing with systems with a large (essentially infinite) numbers of degrees of freedom.<\/p>\n<p>Let us think about a &#8220;case study&#8221; of a string fixed at both ends (x = 0 and x = L).<\/p>\n<p>We can now go about solving for the normal modes of oscillation. \u00a0It is not difficult to check that the following is a solution:<\/p>\n<img src='https:\/\/s0.wp.com\/latex.php?latex=y_n+%28x%2Ct%29+%3D+A_n+%5Ccos+%28%5Comega_n+t+-+%5Cdelta_n%29+%5Csin%5Cleft%28%5Cfrac%7B2+%5Cpi+x%7D%7B%5Clambda_n%7D%5Cright%29&#038;bg=T&#038;fg=000000&#038;s=0' alt='y_n (x,t) = A_n \\cos (\\omega_n t - \\delta_n) \\sin\\left(\\frac{2 \\pi x}{\\lambda_n}\\right)' title='y_n (x,t) = A_n \\cos (\\omega_n t - \\delta_n) \\sin\\left(\\frac{2 \\pi x}{\\lambda_n}\\right)' class='latex' \/>\n<p>Plugging this into the wave equation, we see that the wave equation is solved if <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cfrac%7B%5Comega_n%7D%7Bv%7D+%3D+%5Cfrac%7B2+%5Cpi%7D%7B%5Clambda_n%7D&#038;bg=T&#038;fg=000000&#038;s=0' alt='\\frac{\\omega_n}{v} = \\frac{2 \\pi}{\\lambda_n}' title='\\frac{\\omega_n}{v} = \\frac{2 \\pi}{\\lambda_n}' class='latex' \/><\/p>\n<p>Further more, in order for <img src='https:\/\/s0.wp.com\/latex.php?latex=y_n&#038;bg=T&#038;fg=000000&#038;s=0' alt='y_n' title='y_n' class='latex' \/> to be zero at x=L, we need <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cfrac%7B2+%5Cpi+L%7D%7B%5Clambda_n%7D+%3D+n+%5Cpi&#038;bg=T&#038;fg=000000&#038;s=0' alt='\\frac{2 \\pi L}{\\lambda_n} = n \\pi' title='\\frac{2 \\pi L}{\\lambda_n} = n \\pi' class='latex' \/>, or <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Clambda_n+%3D+%5Cfrac%7B2+L%7D%7Bn%7D&#038;bg=T&#038;fg=000000&#038;s=0' alt='\\lambda_n = \\frac{2 L}{n}' title='\\lambda_n = \\frac{2 L}{n}' class='latex' \/> where n is any positive\u00a0integer.<\/p>\n<p>Thus, for a normal mode solution to the wave equation with the boundary conditions we have chosen (y(x=0) = y(x=L) = 0), we have:<\/p>\n<img src='https:\/\/s0.wp.com\/latex.php?latex=%5Clambda_n+%3D+%5Cfrac%7B2+L%7D%7Bn%7D+%5Ctext%7B+and+%7D+%5Comega_n+%3D+%5Cfrac%7Bn+%5Cpi+v%7D%7BL%7D&#038;bg=T&#038;fg=000000&#038;s=0' alt='\\lambda_n = \\frac{2 L}{n} \\text{ and } \\omega_n = \\frac{n \\pi v}{L}' title='\\lambda_n = \\frac{2 L}{n} \\text{ and } \\omega_n = \\frac{n \\pi v}{L}' class='latex' \/>\n<p>Now since the wave equation is linear, any superposition of the normal mode solutions is also a solution, and we can write any motion of the string as<\/p>\n<img src='https:\/\/s0.wp.com\/latex.php?latex=y%28x%2Ct%29+%3D+%5Csum_%7Bn%3D1%7D%5E%7B%5Cinfty%7D+A_n+%5Ccos%5Cleft%5B+%5Cfrac%7Bn+%5Cpi+v+t%7D%7BL%7D+-+%5Cdelta_n+%5Cright%5D+%5Csin+%5Cleft%5B+%5Cfrac%7Bn+%5Cpi+x%7D%7BL%7D+%5Cright%5D&#038;bg=T&#038;fg=000000&#038;s=0' alt='y(x,t) = \\sum_{n=1}^{\\infty} A_n \\cos\\left[ \\frac{n \\pi v t}{L} - \\delta_n \\right] \\sin \\left[ \\frac{n \\pi x}{L} \\right]' title='y(x,t) = \\sum_{n=1}^{\\infty} A_n \\cos\\left[ \\frac{n \\pi v t}{L} - \\delta_n \\right] \\sin \\left[ \\frac{n \\pi x}{L} \\right]' class='latex' \/>\n<p>(Of course always keep in mind that this only works sufficiently close to equilibrium, where the equation is linear &#8211; for large displacements\/very jagged string configurations\/etc, non-linear terms will be important, and will not be amenable to this type of analysis).<\/p>\n<p>A good question that we will address next week is how to get out the unknown coefficients <img src='https:\/\/s0.wp.com\/latex.php?latex=A_n&#038;bg=T&#038;fg=000000&#038;s=0' alt='A_n' title='A_n' class='latex' \/> and <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdelta_n&#038;bg=T&#038;fg=000000&#038;s=0' alt='\\delta_n' title='\\delta_n' class='latex' \/>. \u00a0This process is referred to as <strong>Fourier analysis<\/strong>, in which the sum over normal modes is referred to as a <strong>Fourier Series\u00a0<\/strong>expansion.<\/p>\n<p>As a final note for this week, we note that we can derive the wave equation in a similar way for systems with spatial extent in more than one dimension. \u00a0A string is 1 D, but we could consider a sheet\/membrane, or the surface of water, as examples of media in which 2D waves can travel. \u00a0Also, sound waves typically are transmitted in 3D spaces (like the classroom, or a shower that you might find yourself singing in).<\/p>\n<p>In 2D, with transverse displacements from equilibrium characterized by a function <img src='https:\/\/s0.wp.com\/latex.php?latex=z%28x%2Cy%2Ct%29&#038;bg=T&#038;fg=000000&#038;s=0' alt='z(x,y,t)' title='z(x,y,t)' class='latex' \/> the wave equation is:<\/p>\n<img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cfrac%7B%5Cpartial%5E2+z%7D%7B%5Cpartial+x%5E2%7D+%2B+%5Cfrac%7B%5Cpartial%5E2+z%7D%7B%5Cpartial+y%5E2%7D+%3D+%5Cfrac%7B1%7D%7Bv%5E2%7D+%5Cfrac%7B%5Cpartial%5E2+z%7D%7B%5Cpartial+t%5E2%7D&#038;bg=T&#038;fg=000000&#038;s=0' alt='\\frac{\\partial^2 z}{\\partial x^2} + \\frac{\\partial^2 z}{\\partial y^2} = \\frac{1}{v^2} \\frac{\\partial^2 z}{\\partial t^2}' title='\\frac{\\partial^2 z}{\\partial x^2} + \\frac{\\partial^2 z}{\\partial y^2} = \\frac{1}{v^2} \\frac{\\partial^2 z}{\\partial t^2}' class='latex' \/>\n<p>In 3D, where a sound wave is characterized by differences in the density of the ambient gas <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cxi+%28x%2Cy%2Cz%2Ct%29&#038;bg=T&#038;fg=000000&#038;s=0' alt='\\xi (x,y,z,t)' title='\\xi (x,y,z,t)' class='latex' \/>, we have<\/p>\n<img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cfrac%7B%5Cpartial%5E2+%5Cxi%7D%7B%5Cpartial+x%5E2%7D+%2B+%5Cfrac%7B%5Cpartial%5E2+%5Cxi%7D%7B%5Cpartial+y%5E2%7D+%2B+%5Cfrac%7B%5Cpartial%5E2+%5Cxi%7D%7B%5Cpartial+z%5E2%7D+%3D+%5Cfrac%7B1%7D%7Bv%5E2%7D+%5Cfrac%7B%5Cpartial+y%7D%7B%5Cpartial+t%5E2%7D&#038;bg=T&#038;fg=000000&#038;s=0' alt='\\frac{\\partial^2 \\xi}{\\partial x^2} + \\frac{\\partial^2 \\xi}{\\partial y^2} + \\frac{\\partial^2 \\xi}{\\partial z^2} = \\frac{1}{v^2} \\frac{\\partial y}{\\partial t^2}' title='\\frac{\\partial^2 \\xi}{\\partial x^2} + \\frac{\\partial^2 \\xi}{\\partial y^2} + \\frac{\\partial^2 \\xi}{\\partial z^2} = \\frac{1}{v^2} \\frac{\\partial y}{\\partial t^2}' class='latex' \/>\n<p>In more than 1D, the patterns of nodes of the normal modes can become quite beautiful, and due to possible degeneracies (i.e. when a 2D medium is a square rather than a rectangle or some other shape), can have unexpected shapes. \u00a0Our DEMO of the square plate with salt (or looked at with the strobe light) showed you some of this richness. \u00a0Here is a picture of a guitar back when it is set into motion associated with one of its normal modes:<\/p>\n<p><a href=\"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-content\/uploads\/sites\/5\/2016\/11\/schmily-150x150.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-405\" src=\"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-content\/uploads\/sites\/5\/2016\/11\/schmily-150x150.jpg\" alt=\"schmily-150x150\" width=\"150\" height=\"150\" \/><\/a><\/p>\n<p>Normal modes are beautiful things, and in our engineering feats, it is of great value to understand the role that they play so that we can either eliminate them when they are problematic, or accentuate them when they are of value.<\/p>\n<p>This weeks Mathematica notebook is here:<\/p>\n<p><a href=\"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-content\/uploads\/sites\/5\/2016\/11\/Week-9.nb\">https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-content\/uploads\/sites\/5\/2016\/11\/Week-9.nb<\/a><\/p>\n<p>In it, you should explore the transition from the discretuum to the continuum, and also seek to understand the structure of the normal modes (i.e. their shape, their relative frequencies, and the behavior of superpositions of normal modes).<\/p>\n","protected":false},"excerpt":{"rendered":"<p>During this week, we derived the wave equation by taking the continuum limit of the discrete system (where there were a finite number of masses attached to a string). \u00a0We started from the equation of motion for a single mass, &hellip; <a href=\"https:\/\/jhubisz.expressions.syr.edu\/phy360\/2016\/11\/03\/week-9\/\">Continue reading <span class=\"meta-nav\">&rarr;<\/span><\/a><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1],"tags":[],"_links":{"self":[{"href":"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-json\/wp\/v2\/posts\/391"}],"collection":[{"href":"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-json\/wp\/v2\/comments?post=391"}],"version-history":[{"count":14,"href":"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-json\/wp\/v2\/posts\/391\/revisions"}],"predecessor-version":[{"id":408,"href":"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-json\/wp\/v2\/posts\/391\/revisions\/408"}],"wp:attachment":[{"href":"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-json\/wp\/v2\/media?parent=391"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-json\/wp\/v2\/categories?post=391"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/jhubisz.expressions.syr.edu\/phy360\/wp-json\/wp\/v2\/tags?post=391"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}