Message Boards Message Boards

0
|
3472 Views
|
3 Replies
|
1 Total Likes
View groups...
Share
Share this post:

How can I solve this impulsive heat equation please?

Posted 3 years ago

I need to solve the following impulsive heat equation:

$$
\left\{\begin{array}{ll}
\partial_{t} \psi(x,t)-\partial_{xx} \psi(x,t)=0, & (x,t)\in (0,1) \times((0, 2) \backslash\{1\}) \\
\psi(0,t)= \psi(1,t)=0, & t \in  (0, 2) \\
\psi(x, 0)= x (1-x), & x \in (0,1) \\
\psi(x, 1)=\psi\left(x, 1^{-}\right)+4, & x \in (0,1)
\end{array}\right.
$$

$1^{-}$ denotes the limit to the left!

This is the code I tried in Mathematica, but it's not giving the results

                               (* problem *)
homogen = If[x = 1, {f[x, t] == Limit[f[x, t], t -> 1, Direction -> "FromBelow"] + 4}, {D[f[x, t], {t, 1}] - D[f[x, t], {x, 2}] == 0}];
                         (*Initial conditions *)
ic = {f[x, 0] == x*(1 - x)};
                   (* Dirichlet boundary conditions*)
bc = {f[0, t] == 0, f[1, t] == 0};
                          (*solution*)
sol = DSolve[{homogen, ic, bc}, f[x, t], {x, 0, 1}, {t, 0, 2}]
POSTED BY: walid zouhair
3 Replies
Posted 3 years ago

You will have to be more specific than

but it's not working!!

Post the code you tried.

POSTED BY: Rohit Namjoshi
Posted 3 years ago

thank you for your reply, but it's not working!!

POSTED BY: walid zouhair
Posted 3 years ago

If[x = 1 ? Perhaps you meant If[x == 1? Try using Piecewise rather than If.

POSTED BY: Rohit Namjoshi
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract