Vibrating String — Question 5

PDF ↗

Question 5

A unit string with fixed ends, unit density and unit tension experiences uniform viscous drag: utt+2βut=uxx,β=2π,u(x,0)=sin⁡(πx)+sin⁡(2πx)+sin⁡(3πx),ut(x,0)=0.u_{tt}+2\beta u_t=u_{xx},\qquad \beta=2\pi,\qquad u(x,0)=\sin(\pi x)+\sin(2\pi x)+\sin(3\pi x),\quad u_t(x,0)=0.

Tasks

  1. Derive the modal equations and classify each of the three active modes as over-, critically or underdamped.

  2. Construct the full solution, satisfying both initial conditions for every mode.

  3. Derive the energy dissipation law for E=12∫01(ut2+ux2)dxE=\tfrac 12\int_0^1(u_t^2+u_x^2)\,dx. Explain why energy is nonincreasing even when a modal displacement oscillates.

  4. Find the leading long-time displacement profile and its exact exponential decay rate. More generally, state when all positive string modes are underdamped.

Original worksheet page 1: question and worked solution for 9-8-005
Show solutionHide solution

Question 5 – Solution

Strategy. The same drag rate competes with different natural frequencies, so one string can exhibit all three damping regimes at once.

Step 1: Classify the modal roots. With u=∑qn(t)sin⁡(nπx)u=\sum q_n(t)\sin(n\pi x), each coefficient satisfies qn″+2βqn′+(nπ)2qn=0q_n''+2\beta q_n'+(n\pi)^2q_n=0. For β=2π\beta=2\pi, the discriminant sign is the sign of 4−n24-n^2. Thus mode one is overdamped, mode two is critical, and mode three is underdamped.

Step 2: Impose both initial conditions. The required coefficient functions, all with qn(0)=1q_n(0)=1 and qn′(0)=0q_n'(0)=0, are q1=e−2πt[cosh(3πt)+23sinh(3πt)],q2=e−2πt(1+2πt),q3=e−2πt[cos(5πt)+25sin(5πt)].\begin{aligned} q_1&=e^{-2\pi t}\left[\cosh(\sqrt 3\pi t)+\frac 2{\sqrt 3}\sinh(\sqrt 3\pi t)\right],\\ q_2&=e^{-2\pi t}(1+2\pi t),\\ q_3&=e^{-2\pi t}\left[\cos(\sqrt 5\pi t)+\frac 2{\sqrt 5}\sin(\sqrt 5\pi t)\right]. \end{aligned} Then u=∑n=13qn(t)sin⁡(nπx)\boxed{u=\sum_{n=1}^3q_n(t)\sin(n\pi x)} satisfies the damped PDE and both endpoints. The added sinh, linear and sine terms are necessary to make the initial velocities zero.

Step 3: Account for dissipated energy. Differentiate EE and integrate uxuxtu_xu_{xt} by parts: E′=[uxut]01+∫01ut(utt−uxx)dx=−2β∫01ut2dx≤0.E'=[u_xu_t]_0^1+\int_0^1u_t(u_{tt}-u_{xx})\,dx =\boxed{-2\beta\int_0^1u_t^2\,dx\le 0}. The boundary term is zero because fixed ends have zero velocity. Oscillation changes how energy is partitioned; drag still removes energy whenever the velocity is not identically zero.

Step 4: Identify the slow surviving mode. Expanding q1q_1 into two exponentials, its slower coefficient is C=(1+2/3)/2C=(1+2/\sqrt 3)/2 and its slower decay rate is γ=(2−3)π\gamma=(2-\sqrt 3)\pi. The other two modes decay as e−2πte^{-2\pi t} times bounded or linear factors. Therefore eγtu(x,t)→Csin⁡(πx)uniformly for 0≤x≤1.\boxed{e^{\gamma t}u(x,t)\longrightarrow C\sin(\pi x) \quad\text{uniformly for }0\le x\le 1.} For a nonnegative general drag parameter, mode nn is underdamped when β<nπ\beta<n\pi; all positive modes are underdamped exactly when 0≤β<π0\le\beta<\pi. The threshold compares drag with the lowest frequency, not with the highest one.

Original worksheet page 2: question and worked solution for 9-8-005

Original worksheet layout. Use Enlarge or open the PDF for a closer view.