Summary of Separation of Variables — Question 4

PDF ↗

Question 4

On 0<x<10<x<1, compare Ht=HxxH_t=H_{xx} and Wtt=WxxW_{tt}=W_{xx} with the moving boundary values H(0,t)=W(0,t)=t,H(1,t)=W(1,t)=2t.H(0,t)=W(0,t)=t,\qquad H(1,t)=W(1,t)=2t. Take H(x,0)=W(x,0)=0H(x,0)=W(x,0)=0 and Wt(x,0)=1+xW_t(x,0)=1+x.

Tasks

  1. For general boundary values A(t),B(t)A(t),B(t), derive the transformed PDE and initial data after subtracting ℓ(x,t)=(1−x)A(t)+xB(t)\ell(x,t)=(1-x)A(t)+xB(t) in each problem.

  2. For the specified data, decide whether the affine lifting itself solves each original PDE.

  3. Construct the complete heat solution using a stationary correction to its transformed forcing and a sine expansion. Derive the required coefficients.

  4. Verify the reconstructed original data and find the long-time heat offset relative to the moving line. Contrast it with the wave field and state the heat release-corner regularity caveat.

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

Question 4 – Solution

Strategy. Homogenizing boundary values creates a source involving the first or second time derivative of the lifting, depending on the PDE.

Step 1: Transform the PDE and every initial datum. Since ℓxx=0\ell_{xx}=0, the zero-endpoint fields v=H−ℓv=H-\ell and z=W−ℓz=W-\ell obey vt−vxx=−ℓt,ztt−zxx=−ℓtt.v_t-v_{xx}=-\ell_t,\qquad z_{tt}-z_{xx}=-\ell_{tt}. Their initial displacements are the respective original profiles minus ℓ(x,0)\ell(x,0); the wave initial velocity is the original velocity minus ℓt(x,0)\ell_t(x,0). Transforming only the endpoints would omit essential source and initial terms.

Step 2: Test the actual moving line. Here ℓ=t(1+x)\ell=t(1+x), so ℓt=1+x\ell_t=1+x but ℓtt=ℓxx=0\ell_{tt}=\ell_{xx}=0. It fails the heat PDE but solves the wave PDE. The wave’s transformed displacement and velocity are both zero, giving W(x,t)=t(1+x).\boxed{W(x,t)=t(1+x).} The heat correction has vt−vxx=−(1+x)v_t-v_{xx}=-(1+x) and v(x,0)=0v(x,0)=0.

Step 3: Solve the forced heat correction. Seek P″=1+xP''=1+x with P(0)=P(1)=0P(0)=P(1)=0. Integration gives P(x)=x22+x36−2x3=x(x−1)(x+4)6.P(x)=\frac{x^2}{2}+\frac{x^3}{6}-\frac{2x}{3} =\frac{x(x-1)(x+4)}6. The sine coefficients of 1+x1+x are an=2[1−2(−1)n]/(nπ)a_n=2[1-2(-1)^n]/(n\pi). Twice integrating P″P'' against a sine gives pn=−an/(nπ)2p_n=-a_n/(n\pi)^2. Thus H=t(1+x)+P(x)+∑n≥12[1−2(−1)n](nπ)3e−(nπ)2tsin⁡(nπx).\boxed{H=t(1+x)+P(x)+ \sum_{n\ge 1}\frac{2[1-2(-1)^n]}{(n\pi)^3} e^{-(n\pi)^2t}\sin(n\pi x).}

Step 4: Verify and interpret the offset. At both ends, PP and every sine vanish. At t=0t=0 the summable sine coefficients reconstruct −P-P uniformly, giving zero initial temperature. For t>0t>0, Gaussian damping permits differentiated convergence; substitution verifies Ht=HxxH_t=H_{xx}. The initial flat field has Hxx=0H_{xx}=0 while the boundary time derivatives are one and two, so joint classical corner smoothness is not asserted. Finally H−ℓ→P<0H-\ell\to P<0 on 0<x<10<x<1, uniformly on the closed interval, whereas W−ℓ=0W-\ell=0 at all times. An identical boundary lifting can have different physical consequences in the two PDEs.

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

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