### Author Topic: HA4-P1  (Read 2189 times)

#### Emily Deibert

• Elder Member
• Posts: 100
• Karma: 0
##### HA4-P1
« on: October 09, 2015, 05:07:44 PM »
Problem 1: http://www.math.toronto.edu/courses/apm346h1/20159/PDE-textbook/Chapter2/S2.6.P.html#problem-2.6.1

(a) Here we have a Dirichlet boundary condition. In the domain ${t>0, x \geqslant ct}$, the solution is given by the D'Alembert formula:

u(x,t) = \frac{1}{2}[g(x+ct) + g(x-ct)] + \frac{1}{2c} \int_{x-ct}^{x+ct} h(y)dy

So plugging in the initial conditions as specified in the problem, we arrive at the solution:

u(x,t) = \phi(x+ct)

In the domain ${0<x<ct}$, the solution is given by:

u(x,t) = \frac{1}{2}g(x+ct) + \frac{1}{2c} \int_{0}^{x+ct} h(y)dy + p\bigg(t - \frac{x}{c}\bigg) - \frac{1}{2}g(ct-x) - \frac{1}{2c} \int_0^{ct-x} h(y)dy

So plugging in the initial and boundary conditions as specified in the problem, we arrive at the solution:

u(x,t) = \phi (x+ct) + \chi\bigg(t-\frac{x}{c}\bigg) - \phi(ct-x)

So the solution to (a) is

u(x,t) =
\begin{cases}
\phi(x+ct) & \text{for } \{t>0, x \geqslant ct\}\\
\phi (x+ct) + \chi\bigg(t-\frac{x}{c}\bigg) - \phi(ct-x) & \text{for } \{0<x<ct\}
\end{cases}

(b) In this case we have a Neumann boundary condition. The solution for the domain ${t>0, x \geqslant ct}$ is the same as in part (a), namely:

u(x,t) = \phi(x+ct)

For the domain ${0<x<ct}$, the solution is given by:

u(x,t) = \frac{1}{2}g(x+ct) + \frac{1}{2c} \int_0^{x+ct} h(y)dy - c \int_0^{t-\frac{x}{c}} q(t')dt' + \frac{1}{2}g(ct-x) + \frac{1}{2c} \int_0^{ct-x} h(y)dy

Plugging in our initial and boundary conditions as specified in the problem, we arrive at the solution:

u(x,t) = \phi (x+ct) - \phi (0) - cX\bigg(t-\frac{x}{c}\bigg) + cX(0) + \phi(ct-x)

where

\int\chi(t')dt' = X(t')

So the solution to (b) is

u(x,t) =
\begin{cases}
\phi(x+ct) & \text{for } \{t>0, x \geqslant ct\}\\
\phi (x+ct) - \phi (0) - cX\bigg(t-\frac{x}{c}\bigg) + cX(0) + \phi(ct-x) & \text{for } \{0<x<ct\}
\end{cases}

« Last Edit: October 10, 2015, 07:23:05 AM by Victor Ivrii »

#### Rong Wei

• Sr. Member
• Posts: 43
• Karma: 0
##### Re: HA4-P1
« Reply #1 on: October 14, 2015, 12:51:25 AM »
My solution is Ï•(x+ct)âˆ’1/câˆ«ctâˆ’x0X(tâˆ’xc)+Ï•(ctâˆ’x)for for {0<x<ct},  Because the boundary condition is Ux|x=0=Ï‡(t), so we should integral Ï‡(t). I plug in the derivative of U respect to x, and I got this answer, may be I'm wrong? I'm not sure,

#### Rong Wei

• Sr. Member
• Posts: 43
• Karma: 0
##### Re: HA4-P1
« Reply #2 on: October 14, 2015, 06:17:43 PM »
Now I understand its the part we discussed should be c integrate Xï¼ˆxï¼‰ from t - x/c to 0