Euler Equations — Question 10

PDF ↗

Question 10

For x≥1x\ge 1, consider a finite interval of forcing in an Euler equation: x2y″−xy′+y=g(x),y(1)=y′(1)=0,g(x)={x,1≤x<e,0,x>e.x^2y''-xy'+y=g(x),\qquad y(1)=y'(1)=0,\qquad g(x)=\begin{cases}x,&1\le x<e,\\0,&x>e.\end{cases} Seek a solution that is C1C^1 across x=ex=e and solves the equation on each side. No classical second derivative at the switching point is assumed.

Tasks

  1. Set t=ln⁡xt=\ln x and Y=etZY=e^tZ. Derive a general definite-integral solution for zero initial data, valid for piecewise continuous forcing gg.

  2. Apply it to the specified forcing and give explicit formulas on 1≤x≤e1\le x\le e and x≥ex\ge e.

  3. Verify continuity of yy and y′y' at ee. Compute the one-sided second derivatives and explain why a C2C^2 solution at that point is impossible.

  4. Determine the large-xx behavior after the forcing stops. Explain why the solution does not return to zero just because the right-hand side becomes zero.

Original worksheet page 1: question and worked solution for 6-4-010
Show solutionHide solution

Question 10 – Solution

Strategy. Remove the repeated exponent, integrate in logarithmic time, and match the finite-forcing transition.

Step 1: Derive the integral response. With D=d/dtD=d/dt, the transformed equation is (D−1)2Y=g(et)(D-1)^2Y=g(e^t), so Z″=e−tg(et)Z''=e^{-t}g(e^t). The data at t=0t=0 give Z(0)=Z′(0)=0Z(0)=Z'(0)=0. Thus Z(t)=∫0t(t−s)e−sg(es)ds,y(x)=x∫1xln⁡(x/r)g(r)r2dr.Z(t)=\int_0^t(t-s)e^{-s}g(e^s)\,ds,\qquad \boxed{y(x)=x\int_1^x\ln(x/r)\frac{g(r)}{r^2}\,dr.} The second formula follows from r=esr=e^s, including ds=dr/rds=dr/r.

Step 2: Evaluate the finite forcing interval. Here Z″=1Z''=1 for 0<t<10<t<1 and Z″=0Z''=0 for t>1t>1. Integrating with the zero initial data gives Z=t2/2Z=t^2/2 up to 11, then Z=t−1/2Z=t-1/2. Therefore y(x)={12x(ln⁡x)2,1≤x≤e,x(ln⁡x−12),x≥e.\boxed{y(x)=\begin{cases} \tfrac 12x(\ln x)^2,&1\le x\le e,\\ x(\ln x-\tfrac 12),&x\ge e. \end{cases}} Both formulas have the same value at ee, and the first satisfies both initial conditions at 11.

Step 3: Check the join and the derivative jump. For Y=etZY=e^tZ, y′=Z+Z′y'=Z+Z' and y″=(Z′+Z″)/xy''=(Z'+Z'')/x. At t=1t=1, both sides have Z=1/2Z=1/2 and Z′=1Z'=1, so y(e)=e/2,y′(e)=3/2,y″(e−)=2/e,y″(e+)=1/e.y(e)=e/2,\qquad y'(e)=3/2,\qquad y''(e^-)=2/e,\quad y''(e^+)=1/e. The finite jump in forcing produces a jump −1/e-1/e in y″y'', while y,y′y,y' remain continuous. The unequal one-sided limits rule out C2C^2 regularity, irrespective of any value assigned to g(e)g(e).

Step 4: Interpret the surviving homogeneous motion. After ee, the solution is the nonzero homogeneous combination xln⁡x−x/2x\ln x-x/2. Thus y(x)/(xln⁡x)→1\boxed{y(x)/(x\ln x)\to 1} as x→∞x\to\infty. The completed forcing interval leaves nonzero value and slope at the switch. Those become initial data for the homogeneous equation; zero subsequent forcing does not reset them to zero.

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

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