Vibrating String — Question 10

PDF ↗

Question 10

A unit fixed-end string with unit tension and density is released from rest with u(x,0)=x(1−x)u(x,0)=x(1-x). Let uNu_N retain only sine modes with index n≤Nn\le N, where N≥1N\ge 1 is an integer, from the exact finite-energy solution. Define the error energy Eerr(t)=12∫01[(ut−uN,t)2+(ux−uN,x)2]dx.E_{\mathrm{err}}(t)=\frac 12\int_0^1\bigl[(u_t-u_{N,t})^2+(u_x-u_{N,x})^2\bigr]\,dx.

Tasks

  1. Derive the Fourier solution and its total energy. State the regularity caveat at the initial boundary corners.

  2. Derive an exact series for EerrE_{\mathrm{err}} and prove the all-time estimate Eerr≤16/(3π4N3)E_{\mathrm{err}}\le 16/(3\pi^4N^3).

  3. Prove the all-time displacement estimate ∥u−uN∥∞≤4/(π3N2)\|u-u_N\|_\infty\le 4/(\pi^3N^2). Find a sufficient integer NN for error at most 10−310^{-3} everywhere.

  4. Find a sufficient integer NN for relative energy error at most 10−410^{-4}. Explain why neither estimate depends on elapsed time and why wave evolution supplies no heat-like smoothing factor.

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

Question 10 – Solution

Strategy. Orthogonality gives a conserved tail energy, while absolute Fourier tails give a different uniform displacement certificate.

Step 1: Construct the solution and normalize energy. Two integrations by parts give bn=2∫01x(1−x)sin⁡(nπx)dx=4(1−(−1)n)(nπ)3.b_n=2\int_0^1x(1-x)\sin(n\pi x)\,dx =\frac{4(1-(-1)^n)}{(n\pi)^3}. Thus u(x,t)=8π3∑n≥1noddsin⁡(nπx)cos⁡(nπt)n3,E=12∫01(1−2x)2dx=16.\begin{gathered} \boxed{u(x,t)=\frac 8{\pi^3}\sum_{\substack{n\ge 1\\n\ {\mathrm{odd}}}} \frac{\sin(n\pi x)\cos(n\pi t)}{n^3},\qquad E=\frac 12\int_0^1(1-2x)^2\,dx=\frac 16.} \end{gathered} Displacement and first-derivative series converge uniformly, and the PDE has a finite-energy interpretation. A globally C2C^2 corner solution is impossible: fixed endpoints would have utt=0u_{tt}=0 there, while the initial second spatial derivative is −2-2.

Step 2: Bound the conserved energy tail. Each omitted mode contributes (nπ)2bn2/4(n\pi)^2b_n^2/4, independent of time. Orthogonality therefore gives Eerr=16π4∑n>Nnodd1n4≤16π4∑n>N1n4≤163π4N3.\begin{gathered} E_{\mathrm{err}}=\frac{16}{\pi^4} \sum_{\substack{n>N\\n\ {\mathrm{odd}}}}\frac 1{n^4} \le\frac{16}{\pi^4}\sum_{n>N}\frac 1{n^4} \le\boxed{\frac{16}{3\pi^4N^3}}. \end{gathered} The last inequality is the decreasing-function integral bound ∑n>Nn−4≤∫N∞s−4ds\sum_{n>N}n^{-4}\le\int_N^\infty s^{-4}\,ds.

Step 3: Certify displacement independently. Using |sin⁡|,|cos⁡|≤1|\sin|,|\cos|\le 1 and the same tail comparison, ∥u−uN∥∞≤8π3∑n>N1n3≤4π3N2.\boxed{\|u-u_N\|_\infty\le\frac 8{\pi^3}\sum_{n>N}\frac 1{n^3} \le\frac 4{\pi^3N^2}.} This holds simultaneously for every xx and every tt. For tolerance δ=10−3\delta=10^{-3}, it suffices to take N≥4/(π3δ)≈11.36N\ge\sqrt{4/(\pi^3\delta)}\approx 11.36, so N=12\boxed{N=12} suffices. The bound includes even indices for convenience; it is sufficient, not a claim of the smallest possible retained set.

Step 4: Choose an energy-based cutoff. Since E=1/6E=1/6, Eerr/E≤32/(π4N3)E_{\mathrm{err}}/E\le 32/(\pi^4N^3). For tolerance η=10−4\eta=10^{-4}, it suffices that N≥(32π4η)1/3≈14.87,N=15.N\ge\left(\frac{32}{\pi^4\eta}\right)^{1/3}\approx 14.87, \qquad \boxed{N=15}. Each wave mode preserves its energy and merely oscillates. There is no factor decaying with n2tn^2t as in diffusion, so waiting does not reduce the truncation’s energy error. These two cutoffs measure different errors and need not agree.

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

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