Question #76480

Obtain a solution of the wave equation
∂^2u(x,t)/∂t^2=16(∂^2u(x,t)/∂x^2)

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 ∂u(x,0)/∂t=0

Expert's answer

ANSWER on Question #76480 – 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}{l}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(4nt) + B_n^{(2)} \sin(4nt)\right) \sin(nx)}


3 STEP: Determine the coefficients B1B_1 and B2B_2.

To do this, you must use the initial conditions.


u(x,0)t=0\frac{\partial u(x,0)}{\partial t} = 0 \rightarrowt(n=1(Bn(1)cos(4nt)+Bn(2)sin(4nt))sin(nx))t=0=\frac{\partial}{\partial t} \left(\sum_{n=1}^{\infty} \left(B_n^{(1)} \cos(4nt) + B_n^{(2)} \sin(4nt)\right) \sin(nx)\right)_{t=0} ==(n=1(4nBn(1)sin(4nt)+4nBn(2)cos(4nt))sin(nx))t=0== \left(\sum_{n=1}^{\infty} \left(-4nB_n^{(1)} \sin(4nt) + 4nB_n^{(2)} \cos(4nt)\right) \sin(nx)\right)_{t=0} ==n=1(4nBn(1)sin(4n0)+4nBn(2)cos(4n0))sin(nx)== \sum_{n=1}^{\infty} \left(-4nB_n^{(1)} \sin(4n \cdot 0) + 4nB_n^{(2)} \cos(4n \cdot 0)\right) \sin(nx) ==n=1(4nBn(1)0+4nBn(2)1)sin(nx)=n=14nBn(2)sin(nx)=0= \sum_{n=1}^{\infty} \left(-4nB_n^{(1)} \cdot 0 + 4nB_n^{(2)} \cdot 1\right) \sin(nx) = \sum_{n=1}^{\infty} 4nB_n^{(2)} \sin(nx) = 0


Then,


4nBn(2)=0Bn(2)=04nB_n^{(2)} = 0 \rightarrow \boxed{B_n^{(2)} = 0}


Conclusion,


u(x,t)=n=1Bn(1)cos(4nt)sin(nx)u(x,t) = \sum_{n=1}^{\infty} B_n^{(1)} \cos(4nt) \sin(nx)u(x,0)=x(πx)u(x, 0) = x(\pi - x) \rightarrowu(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{aligned} 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{aligned}


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)dx\int_{0}^{\pi} \times \left| \sum_{n=1}^{\infty} B_n^{(1)} \sin(nx) = x(\pi - x) \right| \times \sin(mx) \, dxBm(1)=0πx(πx)sin(mx)dx=0π(πxx2)sin(mx)dxB_m^{(1)} = \int_{0}^{\pi} x(\pi - x) \sin(mx) \, dx = \int_{0}^{\pi} (\pi x - x^2) \sin(mx) \, dx \rightarrowBm(1)=0ππxsin(mx)dx0πx2sin(mx)dx=I1I2B_m^{(1)} = \int_{0}^{\pi} \pi x \sin(mx) \, dx - \int_{0}^{\pi} x^2 \sin(mx) \, dx = I_1 - I_2I1=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 \frac{\sin(mx)}{m}\right|_{0}^{\pi} = \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)mmI_{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==π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} I_{2} = \int_{0}^{\pi} \underline{x}^{2} \cdot \underbrace{\sin(mx) \, dx}_{dv} = \left[ \begin{array}{c} u = x^{2} \rightarrow du = 2x \, dx \\ dv = \sin(mx) \, dx \\ v = \frac{-\cos(mx)}{m} \end{array} \right] = \\ = -\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} \underline{x} \cdot \underbrace{\cos(mx) \, dx}_{dv} = \left[ \begin{array}{c} u = x \rightarrow du = dx \\ dv = \cos(mx) \, dx \\ v = \frac{\sin(mx)}{m} \end{array} \right] = \\ = -\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 = \\ = -\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(mx) \, dx\right) = \\ = -\frac{\pi^{2} \cdot (-1)^{m}}{m} + \frac{2}{m} \cdot \left(0 - 0 - \frac{1}{m} \cdot \left(-\frac{\cos(mx)}{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)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}{c} (- 1) ^ {m} - 1 = 0, m = 2 k, k = 1, 2, 3, 4, \ldots \\ (- 1) ^ {m} - 1 = - 2, m = 2 k - 1, k = 1, 2, 3, 4, \ldots \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, \ldots \\ B _ {m} ^ {(1)} = \frac {4}{m ^ {3}}, m = 2 k - 1, k = 1, 2, 3, 4, \ldots \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}{(2 k - 1) ^ {3}}\right) \cos (4 (2 k - 1) t) \sin ((2 k - 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