Question #78896

Obtain a solution of the wave equation
D2u(x,t)/dt2 = 16*d2u(x,t)/ dx2
For 0<= x<= π and t>0, and the following boundary and initial conditions:
U(0,t)=u(π,t)=0,
U(x,0)=x(π-x) and du(x,0)/dt =0

Expert's answer

ANSWER on Question #78896 – Math – Differential Equations

QUESTION

Obtain a solution of the wave equation


2u(x,t)t2=162u(x,t)x2\frac {\partial^ {2} u (x , t)}{\partial t ^ {2}} = 16 \cdot \frac {\partial^ {2} u (x , t)}{\partial x ^ {2}}


for 0xπ0 \leq x \leq \pi and t>0t > 0 and the following boundary and initial conditions:


{u(0,t)=0u(π,t)=0boundary conditions\left\{ \begin{array}{l} u (0, t) = 0 \\ u (\pi , t) = 0 \end{array} - \text{boundary conditions} \right.{u(x,0)=x(πx)u(x,0)t=0initial conditions\left\{ \begin{array}{c} u (x, 0) = x (\pi - x) \\ \frac {\partial u (x , 0)}{\partial t} = 0 \end{array} - \text{initial conditions} \right.

SOLUTION

0 STEP: separation of variables.

Let


u(x,t)=X(x)T(t){2u(x,t)t2=2t2(X(x)T(t))=X(x)d2(T(t))dt2=X(x)T(t)2u(x,t)x2=2x2(X(x)T(t))=T(t)d2(X(x))dx2=X(x)T(t)u (x, t) = X (x) T (t) \rightarrow \left\{ \begin{array}{l} \frac {\partial^ {2} u (x , t)}{\partial t ^ {2}} = \frac {\partial^ {2}}{\partial t ^ {2}} \big (X (x) T (t) \big) = X (x) \cdot \frac {d ^ {2} \big (T (t) \big)}{d t ^ {2}} = X (x) \cdot T ^ {\prime \prime} (t) \\ \frac {\partial^ {2} u (x , t)}{\partial x ^ {2}} = \frac {\partial^ {2}}{\partial x ^ {2}} \big (X (x) T (t) \big) = T (t) \cdot \frac {d ^ {2} \big (X (x) \big)}{d x ^ {2}} = X ^ {\prime \prime} (x) \cdot T (t) \end{array} \right.


Boundary conditions:


{u(0,t)=0u(π,t)=0{u(0,t)=X(0)T(t)=0,t>0u(π,t)=X(π)T(t)=0,t>0{X(0)=0}{X(π)=0\left\{ \begin{array}{l} u (0, t) = 0 \\ u (\pi , t) = 0 \end{array} \right. \to \left\{ \begin{array}{l} u (0, t) = X (0) T (t) = 0, \forall t > 0 \\ u (\pi , t) = X (\pi) T (t) = 0, \forall t > 0 \end{array} \right. \to \boxed {\{X (0) = 0 \}} \\ \left\{ \begin{array}{l} X (\pi) = 0 \end{array} \right.


Then,


2u(x,t)t2=162u(x,t)x2X(x)T(t)=16X(x)T(t)×116X(x)T(t)\frac {\partial^ {2} u (x , t)}{\partial t ^ {2}} = 16 \cdot \frac {\partial^ {2} u (x , t)}{\partial x ^ {2}} \rightarrow X (x) \cdot T ^ {\prime \prime} (t) = 16 \cdot X ^ {\prime \prime} (x) \cdot T (t) | \times \frac {1}{16X(x)T(t)} \rightarrowX(x)T(t)16X(x)T(t)=16X(x)T(t)16X(x)T(t)116T(t)T(t)=X(x)X(x)=λ\frac {X (x) \cdot T ^ {\prime \prime} (t)}{16X(x)T(t)} = \frac {16 \cdot X ^ {\prime \prime} (x) \cdot T (t)}{16X(x)T(t)} \rightarrow \boxed {\frac {1}{16} \cdot \frac {T ^ {\prime \prime} (t)}{T (t)} = \frac {X ^ {\prime \prime} (x)}{X (x)} = - \lambda}


1 STEP: We solve the Sturm-Liouville problem.

(More information: https://en.wikipedia.org/wiki/Sturm%E2%80%93Liouville_theory)

In our case,


{X(x)X(x)=λX(0)=0X(π)=0\left\{ \begin{array}{l} \frac {X ^ {\prime \prime} (x)}{X (x)} = - \lambda \\ X (0) = 0 \\ X (\pi) = 0 \end{array} \right.X(x)X(x)=λX(x)=λX(x)X(x)+λX(x)=0\frac {X ^ {\prime \prime} (x)}{X (x)} = - \lambda \rightarrow X ^ {\prime \prime} (x) = - \lambda X (x) \rightarrow X ^ {\prime \prime} (x) + \lambda X (x) = 0


Let us find the solutions of the given equation in the form


X(x)=ekxX(x)=k2ekxX (x) = e ^ {k x} \rightarrow X ^ {\prime \prime} (x) = k ^ {2} \cdot e ^ {k x}


Then,


X(x)+λX(x)=0k2ekx+λekx=0ekx(k2+λ)=0k2=λX ^ {\prime \prime} (x) + \lambda X (x) = 0 \rightarrow k ^ {2} \cdot e ^ {k x} + \lambda e ^ {k x} = 0 \rightarrow e ^ {k x} (k ^ {2} + \lambda) = 0 \rightarrow k ^ {2} = - \lambdak2=λ[k1=λ=iλk2=λ=iλk ^ {2} = - \lambda \rightarrow \left[\begin{array}{l}k _ {1} = \sqrt {- \lambda} = i \sqrt {\lambda}\\k _ {2} = - \sqrt {- \lambda} = - i \sqrt {\lambda}\end{array}\right.


Then,


X(x)=C1eiλx+C2eiλxA1cos(λx)+A2sin(λx)X (x) = C _ {1} e ^ {i \sqrt {\lambda} x} + C _ {2} e ^ {- i \sqrt {\lambda} x} \equiv A _ {1} \cos (\sqrt {\lambda} x) + A _ {2} \sin (\sqrt {\lambda} x)X(x)=A1cos(λx)+A2sin(λx)\boxed {X (x) = A _ {1} \cos (\sqrt {\lambda} x) + A _ {2} \sin (\sqrt {\lambda} x)}X(0)=0=A1cos(λ0)+A2sin(λ0)=A1cos(0)+A2sin(0)=A11+A20X (0) = 0 = A _ {1} \cos (\sqrt {\lambda} \cdot 0) + A _ {2} \sin (\sqrt {\lambda} \cdot 0) = A _ {1} \cos (0) + A _ {2} \sin (0) = A _ {1} \cdot 1 + A _ {2} \cdot 0 \rightarrowA1=0\boxed {A _ {1} = 0}X(π)=0=A2sin(λπ)sin(λπ)=0λπ=πn,n=1,2,3,X (\pi) = 0 = A _ {2} \sin (\sqrt {\lambda} \pi) \rightarrow \sin (\sqrt {\lambda} \pi) = 0 \rightarrow \sqrt {\lambda} \pi = \pi n, n = 1, 2, 3, \dotsλn=n2,n=1,2,3,\boxed {\lambda_ {n} = n ^ {2}, n = 1, 2, 3, \dots}


Conclusion,


{Xn(x)=Asin(nx)λn=n2n=1,2,3,\left\{ \begin{array}{c} X _ {n} (x) = A \cdot \sin (n x) \\ \lambda_ {n} = n ^ {2} \\ n = 1, 2, 3, \ldots \end{array} \right.


2 STEP: Finding the general solution.


116T(t)T(t)=X(x)X(x)=λn116T(t)T(t)=n2T(t)=16n2T(t)\frac {1}{16} \cdot \frac {T ^ {\prime \prime} (t)}{T (t)} = \frac {X ^ {\prime \prime} (x)}{X (x)} = - \lambda_ {n} \rightarrow \frac {1}{16} \cdot \frac {T ^ {\prime \prime} (t)}{T (t)} = - n ^ {2} \rightarrow T ^ {\prime \prime} (t) = - 16 n ^ {2} T (t) \rightarrowT(t)+16n2T(t)=0T ^ {\prime \prime} (t) + 16 n ^ {2} T (t) = 0


Let us find the solutions of the given equation in the form


T(t)=ektT(t)=k2ektT (t) = e ^ {k t} \rightarrow T ^ {\prime \prime} (t) = k ^ {2} \cdot e ^ {k t}


Then,


T(t)+16n2T(t)=0k2ekt+16n2ekt=0ekt(k2+16n2)=0k2=16n2T ^ {\prime \prime} (t) + 16 n ^ {2} T (t) = 0 \rightarrow k ^ {2} \cdot e ^ {k t} + 16 n ^ {2} e ^ {k t} = 0 \rightarrow e ^ {k t} (k ^ {2} + 16 n ^ {2}) = 0 \rightarrow k ^ {2} = - 16 n ^ {2}k2=16n2[k1=16n2=4ink2=16n2=4ink ^ {2} = - 16 n ^ {2} \rightarrow \left[\begin{array}{c}k _ {1} = \sqrt {- 16 n ^ {2}} = 4 i n\\k _ {2} = - \sqrt {- 16 n ^ {2}} = - 4 i n\end{array}\right.


Then,


Tn(x)=C1e4int+C2e4intA1cos(4nt)+A2sin(4nt)T _ {n} (x) = C _ {1} e ^ {4 i n t} + C _ {2} e ^ {- 4 i n t} \equiv A _ {1} \cos (4 n t) + A _ {2} \sin (4 n t)Tn(x)=An(1)cos(4nt)+An(1)sin(4nt)\boxed {T _ {n} (x) = A _ {n} ^ {(1)} \cos (4 n t) + A _ {n} ^ {(1)} \sin (4 n t)}


Then,


un(x,t)=Xn(x)Tn(t)=(Asin(nx))(An(1)cos(4nt)+An(1)sin(4nt))u _ {n} (x, t) = X _ {n} (x) \cdot T _ {n} (t) = (A \cdot \sin (n x)) \cdot \left(A _ {n} ^ {(1)} \cos (4 n t) + A _ {n} ^ {(1)} \sin (4 n t)\right) \rightarrowun(x,t)=(Bn(1)cos(4nt)+Bn(2)sin(4nt))sin(nx)particularsolution\boxed {u _ {n} (x, t) = \left(B _ {n} ^ {(1)} \cos (4 n t) + B _ {n} ^ {(2)} \sin (4 n t)\right) \sin (n x) - p a r t i c u l a r s o l u t i o n}


Conclusion,


u(x,t)=n=1un(x,t)general solutionu (x, t) = \sum_ {n = 1} ^ {\infty} u _ {n} (x, t) - \text{general solution}u(x,t)=n=1(Bn(1)cos(4nt)+Bn(2)sin(4nt))sin(nx)\boxed {u (x, t) = \sum_ {n = 1} ^ {\infty} \left(B _ {n} ^ {(1)} \cos (4 n t) + B _ {n} ^ {(2)} \sin (4 n t)\right) \sin (n x)}


3 STEP: Determine the coefficients Bn(1)B_{n}^{(1)} and Bn(2)B_{n}^{(2)}.

To do this, you must use the initial conditions.


u(x,0)t=0t(n=1(Bn(1)cos(4nt)+Bn(2)sin(4nt))sin(nx))t=0==(n=1(4nBn(1)sin(4nt)+4nBn(2)cos(4nt))sin(nx))t=0==n=1(4nBn(1)sin(4n0)+4nBn(2)cos(4n0))sin(nx)==n=1(4nBn(1)0+4nBn(2)1)sin(nx)=n=14nBn(2)sin(nx)=0\begin{array}{l} \frac {\partial u (x , 0)}{\partial t} = 0 \rightarrow \\ \frac {\partial}{\partial t} \left(\sum_ {n = 1} ^ {\infty} \left(B _ {n} ^ {(1)} \cos (4 n t) + B _ {n} ^ {(2)} \sin (4 n t)\right) \sin (n x)\right) _ {t = 0} = \\ = \left(\sum_ {n = 1} ^ {\infty} \left(- 4 n B _ {n} ^ {(1)} \sin (4 n t) + 4 n B _ {n} ^ {(2)} \cos (4 n t)\right) \sin (n x)\right) _ {t = 0} = \\ = \sum_ {n = 1} ^ {\infty} \left(- 4 n B _ {n} ^ {(1)} \sin (4 n \cdot 0) + 4 n B _ {n} ^ {(2)} \cos (4 n \cdot 0)\right) \sin (n x) = \\ = \sum_ {n = 1} ^ {\infty} \left(- 4 n B _ {n} ^ {(1)} \cdot 0 + 4 n B _ {n} ^ {(2)} \cdot 1\right) \sin (n x) = \sum_ {n = 1} ^ {\infty} 4 n B _ {n} ^ {(2)} \sin (n x) = 0 \\ \end{array}


Then,


4nBn(2)=0Bn(2)=04 n B _ {n} ^ {(2)} = 0 \rightarrow \boxed {B _ {n} ^ {(2)} = 0}


Conclusion,


u(x,t)=n=1Bn(1)cos(4nt)sin(nx)\boxed {u (x, t) = \sum_ {n = 1} ^ {\infty} B _ {n} ^ {(1)} \cos (4 n t) \sin (n x)}u(x,0)=x(πx)u(x,0)=(n=1Bn(1)cos(4nt)sin(nx))t=0=n=1Bn(1)cos(4n0)sin(nx)==n=1Bn(1)1sin(nx)n=1Bn(1)sin(nx)=x(πx)\begin{array}{l} u(x, 0) = x(\pi - x) \rightarrow \\ u(x, 0) = \left(\sum_{n=1}^{\infty} B_n^{(1)} \cos(4nt) \sin(nx)\right)_{t=0} = \sum_{n=1}^{\infty} B_n^{(1)} \cos(4n \cdot 0) \sin(nx) = \\ = \sum_{n=1}^{\infty} B_n^{(1)} \cdot 1 \cdot \sin(nx) \rightarrow \sum_{n=1}^{\infty} B_n^{(1)} \sin(nx) = x(\pi - x) \end{array}


As we know


0πsin(mx)sin(nx)dx={π2,n=m0,nm\int_{0}^{\pi} \sin(mx) \cdot \sin(nx) \, dx = \begin{cases} \dfrac{\pi}{2}, & n = m \\ 0, & n \neq m \end{cases}


In our case,


0π×n=1Bn(1)sin(nx)=x(πx)×sin(mx)dxBm(1)=0πx(πx)sin(mx)dx=0π(πxx2)sin(mx)dxBm(1)=0ππxsin(mx)dx0πx2sin(mx)dx=I1I2\begin{aligned} \int_{0}^{\pi} \times \left| \sum_{n=1}^{\infty} B_n^{(1)} \sin(nx) = x(\pi - x) \right| & \times \sin(mx) \, dx \\ B_m^{(1)} = \int_{0}^{\pi} x(\pi - x) \sin(mx) \, dx &= \int_{0}^{\pi} (\pi x - x^2) \sin(mx) \, dx \rightarrow \\ B_m^{(1)} = \int_{0}^{\pi} \pi x \sin(mx) \, dx &- \int_{0}^{\pi} x^2 \sin(mx) \, dx = I_1 - I_2 \end{aligned}I1=0ππxsin(mx)dx=π0πxsin(mx)dxdv=[u=xdu=dxdv=sin(mx)dxv=cos(mx)m]==π(xcos(mx)m0π0πcos(mx)dxm==π(πcos(πm)m(0cos(m0)m)+1m0πcos(mx)dx)==π(π(1)mm+1msin(mx)m0π)=π(π(1)mm+sin(mπ)m2sin(m0)m2)==π(π(1)mm+00)=π2(1)mm\begin{array}{l} I_{1} = \int_{0}^{\pi} \pi x \sin(mx) \, dx = \pi \cdot \int_{0}^{\pi} \underbrace{x} \cdot \underbrace{\sin(mx) \, dx}_{dv} = \begin{bmatrix} u = x \to du = dx \\ dv = \sin(mx) \, dx \\ v = \frac{-\cos(mx)}{m} \end{bmatrix} = \\ = \pi \cdot \left(-\frac{x \cdot \cos(mx)}{m}\right|_{0}^{\pi} - \int_{0}^{\pi} \frac{-\cos(mx) \, dx}{m} = \\ = \pi \cdot \left(-\frac{\pi \cdot \cos(\pi m)}{m} - \left(-\frac{0 \cdot \cos(m \cdot 0)}{m}\right) + \frac{1}{m} \int_{0}^{\pi} \cos(mx) \, dx\right) = \\ = \pi \cdot \left(-\frac{\pi \cdot (-1)^{m}}{m} + \frac{1}{m} \cdot \left.\frac{\sin(mx)}{m}\right|_{0}^{\pi}\right) = \pi \cdot \left(-\frac{\pi \cdot (-1)^{m}}{m} + \frac{\sin(m \cdot \pi)}{m^{2}} - \frac{\sin(m \cdot 0)}{m^{2}}\right) = \\ = \pi \cdot \left(-\frac{\pi \cdot (-1)^{m}}{m} + 0 - 0\right) = -\frac{\pi^{2} \cdot (-1)^{m}}{m} \end{array}


Conclusion,


I1=π2(1)mm\boxed{I_{1} = -\frac{\pi^{2} \cdot (-1)^{m}}{m}}I2=0πx2sin(mx)dxdv=[u=x2du=2xdxdv=sin(mx)dxv=cos(mx)m]==x2cos(mx)m0π0π2xcos(mx)dxm==π2cos(mπ)m(02cos(m0)m)+2m0πxcos(mx)dx==π2(1)mm+2m0πxcos(mx)dxdv=[u=xdu=dxdv=cos(mx)dxv=sin(mx)m]==π2(1)mm+2m(xsin(mx)m0π0πsin(mx)mdx=\begin{array}{l} I_{2} = \int_{0}^{\pi} \underbrace{x^{2}} \cdot \underbrace{\sin(mx) \, dx}_{dv} = \begin{bmatrix} u = x^{2} \to du = 2x \, dx \\ dv = \sin(mx) \, dx \\ v = \frac{-\cos(mx)}{m} \end{bmatrix} = \\ = -\frac{x^{2} \cdot \cos(mx)}{m}\Bigg|_{0}^{\pi} - \int_{0}^{\pi} \frac{-2x \cos(mx) \, dx}{m} = \\ = -\frac{\pi^{2} \cdot \cos(m\pi)}{m} - \left(-\frac{0^{2} \cdot \cos(m \cdot 0)}{m}\right) + \frac{2}{m} \cdot \int_{0}^{\pi} x \cos(mx) \, dx = \\ = -\frac{\pi^{2} \cdot (-1)^{m}}{m} + \frac{2}{m} \cdot \int_{0}^{\pi} \underbrace{x} \cdot \underbrace{\cos(mx) \, dx}_{dv} = \begin{bmatrix} u = x \to du = dx \\ dv = \cos(mx) \, dx \\ v = \frac{\sin(mx)}{m} \end{bmatrix} = \\ = -\frac{\pi^{2} \cdot (-1)^{m}}{m} + \frac{2}{m} \cdot \left(\frac{x \cdot \sin(mx)}{m}\right|_{0}^{\pi} - \int_{0}^{\pi} \frac{\sin(mx)}{m} \, dx = \\ \end{array}=π2(1)mm+2m(πsin(mπ)m0sin(m0)m1m0πsin(mx)dx)==π2(1)mm+2m(001m(cos(mx)m0π)==π2(1)mm+2m3(cos(mπ)cos(m0))==π2(1)mm+2m3((1)m1)\begin{array}{l} = - \frac {\pi^ {2} \cdot (- 1) ^ {m}}{m} + \frac {2}{m} \cdot \left(\frac {\pi \cdot \sin (m \pi)}{m} - \frac {0 \cdot \sin (m \cdot 0)}{m} - \frac {1}{m} \int_ {0} ^ {\pi} \sin (m x) d x\right) = \\ = - \frac {\pi^ {2} \cdot (- 1) ^ {m}}{m} + \frac {2}{m} \cdot \left(0 - 0 - \frac {1}{m} \cdot \left(- \frac {\cos (m x)}{m} \right| _ {0} ^ {\pi}\right) = \\ = - \frac {\pi^ {2} \cdot (- 1) ^ {m}}{m} + \frac {2}{m ^ {3}} \cdot (\cos (m \pi) - \cos (m \cdot 0)) = \\ = - \frac {\pi^ {2} \cdot (- 1) ^ {m}}{m} + \frac {2}{m ^ {3}} \cdot ((- 1) ^ {m} - 1) \\ \end{array}I2=π2(1)mm+2m3((1)m1)\boxed {I _ {2} = - \frac {\pi^ {2} \cdot (- 1) ^ {m}}{m} + \frac {2}{m ^ {3}} \cdot ((- 1) ^ {m} - 1)}


Then,


Bm(1)=I1I2=π2(1)mm(π2(1)mm+2m3((1)m1))B _ {m} ^ {(1)} = I _ {1} - I _ {2} = - \frac {\pi^ {2} \cdot (- 1) ^ {m}}{m} - \left(- \frac {\pi^ {2} \cdot (- 1) ^ {m}}{m} + \frac {2}{m ^ {3}} \cdot ((- 1) ^ {m} - 1)\right)\rightarrowBm(1)=π2(1)mm+π2(1)mm2m3((1)m1)B _ {m} ^ {(1)} = - \frac {\pi^ {2} \cdot (- 1) ^ {m}}{m} + \frac {\pi^ {2} \cdot (- 1) ^ {m}}{m} - \frac {2}{m ^ {3}} \cdot ((- 1) ^ {m} - 1) \rightarrowBm(1)=2((1)m1)m3\boxed {B _ {m} ^ {(1)} = \frac {- 2 \cdot ((- 1) ^ {m} - 1)}{m ^ {3}}}


As we know


{(1)m1=0,m=2k,k=1,2,3,4,(1)m1=2,m=2k1,k=1,2,3,4,\left\{ \begin{array}{l} (- 1) ^ {m} - 1 = 0, m = 2 k, k = 1, 2, 3, 4, \dots \\ (- 1) ^ {m} - 1 = - 2, m = 2 k - 1, k = 1, 2, 3, 4, \dots \end{array} \right. \to{Bm(1)=0,m=2k,k=1,2,3,4,Bm(1)=4m3,m=2k1,k=1,2,3,4,\left\{ \begin{array}{c} B _ {m} ^ {(1)} = 0, m = 2 k, k = 1, 2, 3, 4, \dots \\ B _ {m} ^ {(1)} = \frac {4}{m ^ {3}}, m = 2 k - 1, k = 1, 2, 3, 4, \dots \end{array} \right.


Conclusion,


u(x,t)=n=1Bn(1)cos(4nt)sin(nx)=n=1(2((1)n1)n3)cos(4nt)sin(nx)u (x, t) = \sum_ {n = 1} ^ {\infty} B _ {n} ^ {(1)} \cos (4 n t) \sin (n x) = \sum_ {n = 1} ^ {\infty} \left(\frac {- 2 \cdot ((- 1) ^ {n} - 1)}{n ^ {3}}\right) \cos (4 n t) \sin (n x)u(x,t)=k=1(4(2k1)3)cos(4(2k1)t)sin((2k1)x)\boxed {u (x, t) = \sum_ {k = 1} ^ {\infty} \left(\frac {4}{(2 k - 1) ^ {3}}\right) \cos (4 (2 k - 1) t) \sin ((2 k - 1) x)}


ANSWER:


u(x,t)=k=1(4(2k1)3)cos(4(2k1)t)sin((2k1)x)u(x,t) = \sum_{k=1}^{\infty} \left(\frac{4}{(2k-1)^3}\right) \cos(4(2k-1)t) \sin((2k-1)x)


Answer provided by https://www.AssignmentExpert.com

Need a fast expert's response?

Submit order

and get a quick answer at the best price

for any assignment or question with DETAILED EXPLANATIONS!

LATEST TUTORIALS
APPROVED BY CLIENTS