Eigenvalues and Eigenfunctions — Question 10

PDF ↗

Question 10

A two-material interval has stiffness p(x)={1,0<x<1/2,4,1/2<x<1,−(py′)′=λy,y(0)=y(1)=0.p(x)=\begin{cases}1,&0<x<1/2,\\4,&1/2<x<1,\end{cases} \qquad -(py')'=\lambda y,\qquad y(0)=y(1)=0. Require yy and the flux py′py' to be continuous at x=1/2x=1/2. Solutions are C2C^2 on each side. Ordinary slopes need not be continuous across the interface.

Tasks

  1. Use integration by parts on the two subintervals to show that every real eigenvalue is positive. Explain cancellation of the interface terms.

  2. Set k=λk=\sqrt\lambda and write the endpoint-compatible sine on each side. Derive the two interface equations for their amplitudes and a determinant condition.

  3. With t=k/4t=k/4, factor the determinant and characterize all positive eigenvalues. Treat sin⁡t=0\sin t=0 separately to avoid losing modes when solving for the amplitude ratio.

  4. Identify the smallest eigenvalue and give its eigenfunction with left amplitude one. Verify continuity and the slope ratio at the interface, and sketch the mode with its physical corner.

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

Question 10 – Solution

Strategy. Match displacement and stiffness-weighted slope, then preserve every branch of the determinant equation.

Step 1: Exclude nonpositive eigenvalues. Adding the two integrations by parts gives λ∫01y2dx=∫01p(y′)2dx\lambda\int_0^1y^2dx=\int_0^1p(y')^2dx: endpoint values vanish and the interface terms cancel because both yy and py′py' match. The right side is positive for a nonzero admissible function; equality would make it constant on both sides, hence identically zero. Therefore λ>0\lambda>0.

Step 2: Form the interface system. The endpoint-compatible pieces are Asin⁡(kx)A\sin(kx) on the left and Bsin⁡(k(1−x)/2)B\sin(k(1-x)/2) on the right. With t=k/4t=k/4, matching gives Asin⁡2t=Bsin⁡t,Acos⁡2t=−2Bcos⁡t.A\sin 2t=B\sin t,\qquad A\cos 2t=-2B\cos t. The second equation uses flux, not ordinary derivative. The determinant vanishes exactly when 2sin⁡2tcos⁡t+cos⁡2tsin⁡t=02\sin 2t\cos t+\cos 2t\sin t=0.

Step 3: Preserve both spectral branches. Factoring yields sin⁡t(6cos⁡2t−1)=0\sin t(6\cos^2t-1)=0. Thus the complete spectrum is λ=16t2,t>0,t=mπ(m≥1)orcos⁡2t=1/6.\boxed{\lambda=16t^2,\qquad t>0,\quad t=m\pi\ (m\geq 1)\ \text{or}\ \cos^2t=1/6.} At t=mπt=m\pi, the first matching equation is zero and the second gives B=(−1)m+1A/2B=(-1)^{m+1}A/2. On the other branch, sin⁡t≠0\sin t\ne 0 and B=2Acos⁡tB=2A\cos t. Each branch gives a nonzero solution and a one-dimensional eigenspace. The apparent root t=0t=0 is excluded by Step 1.

Step 4: Construct the first mode and its corner. The smallest positive root is t0=arccos⁡(1/6)∈(0,π/2)t_0=\arccos(1/\sqrt 6)\in(0,\pi/2). Hence λ1=16t02\boxed{\lambda_1=16t_0^2} and, taking A=1A=1, y={sin⁡(4t0x),0≤x≤1/2,2/3sin⁡(2t0(1−x)),1/2≤x≤1.\boxed{y=\begin{cases}\sin(4t_0x),&0\leq x\leq 1/2,\\ \sqrt{2/3}\sin(2t_0(1-x)),&1/2\leq x\leq 1.\end{cases}} At the join both values equal 5/3\sqrt 5/3. The one-sided slopes are −8t0/3-8t_0/3 and −2t0/3-2t_0/3, so y′(1/2−)=4y′(1/2+)y'(1/2-)=4y'(1/2+) as required by flux continuity. The slope change is physical; smoothing it would violate the model.

See the diagram in the original worksheet below.

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

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