1
$\begingroup$

I have the following setup for the linear heat equation on an infinite rod (see below). This is a pretty standard setup and the initial condition is $e^{-x^2}$. This initial condition is not special, it just provides a way to ensure that the code works.

Clear[x, t];
With[{u = u[t, x]}, eq = D[u, t] == k D[D[u, x], x];
 ic = u == Exp[-x^2] /. t -> 0;
 bc = D[u, x] == 0 /. x -> 0;]
asol = DSolveValue[{eq, ic, bc}, u, {t, x}, Assumptions -> {k > 0}];
asol[t, x]

What I want to do is test different ideas and so I want to modify this linear heat equation into the nonlinear heat equation or another PDE that is similar but is different than a linear PDE.

The nonlinear heat equation is defined as the following

\begin{equation} \frac{\partial u}{\partial t} = \frac{\partial}{\partial x}\left[g(u) \frac{\partial u}{\partial x}\right] \end{equation}

and further literature on the nonlinear heat equation can be found here.

If we assume that $g(u) = u$ then I would expect that I could modify the setup as:

Clear[x, t];
With[{u = u[t, x]}, eq = D[u, t] == k D[ u D[u, x], x];
 ic = u == Exp[-x^2] /. t -> 0;
 bc = D[u, x] == 0 /. x -> 0;]
asol = DSolveValue[{eq, ic, bc}, u, {t, x}, Assumptions -> {k > 0}];
asol[t, x]

However, Mathematica simply regurgitates the code and does not produce an output.

I also wish to have a different type of setup such as:

\begin{equation} \frac{\partial u}{\partial t} = \frac{\partial}{\partial x}\left[h(x) \frac{\partial u}{\partial x}\right] \end{equation}

and modify the first block of code for this as well.

If anyone can help me understand why my latter setup doesn't work and can show me how to set it up correctly, I would be very thankful.

Also this is not for a class. It is for my own personal understanding.

$\endgroup$
1
  • 2
    $\begingroup$ The chance to obtain an analytic solution for a nonlinear PDE is basically 0. Try NDSolve. Keep in mind that you have to supply appropriate boundary conditions. However, I am not 100% whether the newest version of Mathematica can handle this nonlinearity. Best, you consult this and this. $\endgroup$ Jan 8, 2020 at 15:11

1 Answer 1

4
$\begingroup$

You can find examples of this for the numerical solver NDSolve in the

FEMDocumentation/tutorial/NonlinearFiniteElementVerificationTests

Specifically the following three tests are taken from the source you give:

There are many more tests and many of those have analytical solutions.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.