Partial Differential Equations: Exact Solutions Subject to Boundary Conditions This document gives examples of Fourier series and integral transform (Laplace and Fourier) solutions to problems involving a PDE and boundary and/or initial conditions. It also describes how, for certain problems, pdsolve can automatically adjust the arbitrary functions and constants entering the solution of the partial differential equations (PDEs) such that the boundary conditions (BCs) are satisfied.

All Products Maple MapleSim

Home : Support : Online Help : Applications and Example Worksheets : Differential Equations : examples/pdsolve_boundaryconditions

Partial Differential Equations: Exact Solutions Subject to Boundary Conditions
This document gives examples of Fourier series and integral transform (Laplace and Fourier) solutions to problems involving a PDE and boundary and/or initial conditions. It also describes how, for certain problems, pdsolve can automatically adjust the arbitrary functions and constants entering the solution of the partial differential equations (PDEs) such that the boundary conditions (BCs) are satisfied.

Two Examples Solved Automatically and Step-by-step

Example 1: Adjusting an Arbitrary Function _F1(...)

>	$restart &semi;$ $with (PDEtools) &colon;$

>	$U ≔ diff_table (u (x, t)) &colon;$

A substance convecting through a tube at velocity c (for example, a radioactive decaying chemical dissolved in water) has its density $u (x, t)$ modeled by a convection-decay type of equation

>	${pde}_{1} ≔ U_{t} + c U_{x} = - λ U []$

${pde}_{1} ≔ c (\frac{\partial}{\partial x} u (x, t)) + \frac{\partial}{\partial t} u (x, t) = - λ u (x, t)$

(1)

where the size of $λ$ indicates how fast the substance is decaying. In this case, the boundary condition is a generic function of $x$ indicating the value of the density at $t = 0$ .

>	${bc}_{1} ≔ eval (U [], t = 0) = φ (x)$

${bc}_{1} ≔ u (x, 0) = φ (x)$

(2)

Step-by-step Solution

The following is an exact solution for ${pde}_{1}$ . For now, ${bc}_{1}$ is not considered.

>	${sol}_{1} ≔ pdsolve ({pde}_{1})$

${sol}_{1} ≔ u (x, t) =_F1 (\frac{t c - x}{c}) {&ExponentialE;}^{- \frac{λ x}{c}}$

(3)

Perform the following steps to find the form of the mapping of $_F1$ that matches ${bc}_{1}$ when $t = 0.$

Consider $t = 0.$

>	$e1 ≔ eval (\binom{{sol}_{1}}{} \| \binom{}{t = 0}, {bc}_{1})$

$e1 ≔ φ (x) =_F1 (- \frac{x}{c}) {&ExponentialE;}^{- \frac{λ x}{c}}$

(4)

Introduce a new variable q equal to the argument of _F1.

>	$tr ≔ q = - \frac{x}{c} &semi;$ $itr ≔ isolate (tr, x)$

$tr ≔ q = - \frac{x}{c}$

$itr ≔ x = - q c$

(5)

Evaluate at this value of x and then isolate _F1(q). The following mapping is obtained.

>	$eval (e1, itr)$

$φ (- q c) =_F1 (q) {&ExponentialE;}^{λ q}$

(6)

>	$e2 ≔ isolate (,_F1 (q))$

$e2 ≔_F1 (q) = \frac{φ (- q c)}{{&ExponentialE;}^{λ q}}$

(7)

>	$ans_F ≔_F1 = unapply (rhs (), q)$

$ans_F ≔_F1 = (q \mapsto \frac{φ (- q c)}{{&ExponentialE;}^{λ q}})$

(8)

The solution to the PDE and BC is obtained evaluating the general solution at this value of _F1.

>	${sol}_{1}$

$u (x, t) =_F1 (\frac{t c - x}{c}) {&ExponentialE;}^{- \frac{λ x}{c}}$

(9)

>	${ans}_{1} ≔ simplify (eval (, ans_F), \exp)$

${ans}_{1} ≔ u (x, t) = φ (- t c + x) {&ExponentialE;}^{- λ t}$

(10)

Verify this result using pdetest.

>	$pdetest ({ans}_{1}, [{pde}_{1}, {bc}_{1}])$

$[0, 0]$

(11)

Summary: The system of PDEs and BCs is

>	${sys}_{1} ≔ [{pde}_{1}, {bc}_{1}]$

${sys}_{1} ≔ [c (\frac{\partial}{\partial x} u (x, t)) + \frac{\partial}{\partial t} u (x, t) = - λ u (x, t), u (x, 0) = φ (x)]$

(12)

The following is the solution computed by pdsolve automatically by performing the steps explained in the Step-by-step Solution section above.

>	$pdsolve ({sys}_{1})$

$u (x, t) = φ (- t c + x) {&ExponentialE;}^{- λ t}$

(13)

The approach used in this example serves also to address nonlinear equations of the form $\frac{\partial}{\partial t} u (x, t) + c (\frac{\partial}{\partial x} u (x, t)) = Φ (x, t, u (x, t))$

Example 2: Periodic Boundary Conditions

Consider the heat equation that models the temperature $u (x, t)$

>	${pde}_{2} ≔ U_{t} = κ U_{x, x}$

${pde}_{2} ≔ \frac{\partial}{\partial t} u (x, t) = κ (\frac{\partial^{2}}{\partial x^{2}} u (x, t))$

(14)

where $κ$ is any constant of proportionality.

Assume that the boundary conditions at $x = 0$ and $x = L$ are as follows:

>	${bc}_{2} ≔ u (0, t) = 0, u (L, t) = 0$

${bc}_{2} ≔ u (0, t) = 0, u (L, t) = 0$

(15)

Step-by-step Solution

Similar to Example 1, the following is a solution for ${pde}_{2}$ . For now, ${bc}_{2}$ is not considered.

>	${sol}_{2} ≔ pdsolve ({pde}_{2}, build)$

${sol}_{2} ≔ u (x, t) =_C3 {&ExponentialE;}^{κ {_c}_{1} t}_C1 {&ExponentialE;}^{\sqrt{{_c}_{1}} x} + \frac{_C3 {&ExponentialE;}^{κ {_c}_{1} t}_C2}{{&ExponentialE;}^{\sqrt{{_c}_{1}} x}}$

(16)

This time, pdsolve computes a complete (that is, not the most general) solution using separation of variables. The solution appears in terms of arbitrary constants $_C1,_C2,_C3$ and the separation variable ${_c}_{1}$ .

The integration constants can be redefined to absorb $_C3 &colon;$

>	$zz ≔_C1 = \frac{_C1}{_C3},_C2 = \frac{_C2}{_C3}$

$zz ≔_C1 = \frac{_C1}{_C3},_C2 = \frac{_C2}{_C3}$

(17)

Now, evaluate the solution at $x = 0$ and $x = L$ and impose the given boundary conditions to build a system (e1 and e2, so two equations) to be solved for these three constants.

>	$e1 ≔ eval (eval (\binom{{sol}_{2}}{} \| \binom{}{x = 0}, ({bc}_{2}) [1]), \{zz\})$

$e1 ≔ 0 = {&ExponentialE;}^{κ {_c}_{1} t}_C1 + {&ExponentialE;}^{κ {_c}_{1} t}_C2$

(18)

>	$e2 ≔ eval (eval (\binom{{sol}_{2}}{} \| \binom{}{x = L}, ({bc}_{2}) [2]), \{zz\})$

$e2 ≔ 0 = {&ExponentialE;}^{κ {_c}_{1} t}_C1 {&ExponentialE;}^{\sqrt{{_c}_{1}} L} + \frac{{&ExponentialE;}^{κ {_c}_{1} t}_C2}{{&ExponentialE;}^{\sqrt{{_c}_{1}} L}}$

(19)

Solve for the constants.

>	$_EnvAllSolutions ≔ true$

$_EnvAllSolutions ≔ true$

(20)

>	$solve (\{e1, e2\}, \{_C1,_C2, {_c}_{1}\})$

$\{_C1 = 0,_C2 = 0, {_c}_{1} = {_c}_{1}\}, \{_C1 = -_C2,_C2 =_C2, {_c}_{1} = - \frac{π^{2} {_Z1}^{2}}{L^{2}}\}$

(21)

The first solution above leads to a solution of no interest: $u (x, t) = 0$ . The second solution above is the one we want.

>	${ans}_{2} ≔ eval ({sol}_{2}, [2])$

${ans}_{2} ≔ u (x, t) = -_C3 {&ExponentialE;}^{- \frac{κ π^{2} {_Z1}^{2} t}{L^{2}}}_C2 {&ExponentialE;}^{\sqrt{- \frac{π^{2} {_Z1}^{2}}{L^{2}}} x} + \frac{_C3 {&ExponentialE;}^{- \frac{κ π^{2} {_Z1}^{2} t}{L^{2}}}_C2}{{&ExponentialE;}^{\sqrt{- \frac{π^{2} {_Z1}^{2}}{L^{2}}} x}}$

(22)

Here pdetest verifies the answer for pde[2] and the first boundary condition, but fails with the second one.

>	$pdetest ({ans}_{2}, [{pde}_{2}, {bc}_{2}])$

$[0, 0, -_C3_C2 ({&ExponentialE;}^{2 π \sqrt{- \frac{1}{L^{2}}} L |_Z1|} - 1) {&ExponentialE;}^{- \frac{π_Z1 (L^{3} signum (_Z1) \sqrt{- \frac{1}{L^{2}}} + π_Z1 κ t)}{L^{2}}}]$

(23)

However, the answer computed is correct. The second boundary condition can be constructed by evaluating the answer at $x = L$ . This task can be performed within pdetest.

>	$simplify (\binom{{ans}_{2}}{} \| \binom{}{x = L}) assuming L :: real$

$u (L, t) = 0$

(24)

Summary: The system with periodic boundary conditions is

>	${sys}_{2} ≔ [{pde}_{2}, {bc}_{2}]$

${sys}_{2} ≔ [\frac{\partial}{\partial t} u (x, t) = κ (\frac{\partial^{2}}{\partial x^{2}} u (x, t)), u (0, t) = 0, u (L, t) = 0]$

(25)

The solution computed by pdsolve automatically by performing the steps summarized in the Step-by-step Solution section above:

>	$pdsolve ({sys}_{2})$

$u (x, t) = \sum_{n = 1}^{\infty}_C1 (n) csgn (\frac{1}{L}) \sin (\frac{n π x}{L}) {&ExponentialE;}^{- \frac{π^{2} n^{2} κ t}{L^{2}}}$

(26)

>	$pdetest (, {sys}_{2})$

$[0, 0, 0]$

(27)

Notes on this PDE and Initial Values / Boundary Conditions Problem

Given a problem with PDEs and some BCs, note the following:

•	Depending on how you write a general PDE solution, it becomes possible or nearly impossible to adjust it to match the BCs.

•	Certain particular PDE solutions are frequently easier to adjust than general solutions.

•	Some combinations, classified in the literature, of "certain types of PDEs" with "certain type of BCs" can be solved systematically by writing the PDE solution in a special form.

•	By passing PDEs alone to pdsolve, it is not possible to predict the form of the PDE solution. On the other hand, pdsolve accepts a HINT argument so that solutions of a "certain type" can be requested.

Three Textbook Examples

A bar loses heat across its lateral boundary at a rate proportional to the temperature $u (x, t)$ and its BCs state that the temperature at the two ends of the bar are equal to 1.

This problem is modeled by the following equation.

$restart$

>	${sys}_{3} ≔ [\frac{\partial}{\partial t} u (x, t) = k (\frac{\partial^{2}}{\partial x^{2}} u (x, t)) - a u (x, t), u (0, t) = 1, u (1, t) = 1]$

${sys}_{3} ≔ [\frac{\partial}{\partial t} u (x, t) = k (\frac{\partial^{2}}{\partial x^{2}} u (x, t)) - a u (x, t), u (0, t) = 1, u (1, t) = 1]$

(28)

>	$pdsolve ({sys}_{3})$

$u (x, t) = \frac{{&ExponentialE;}^{\sqrt{\frac{a}{k}} x} + {&ExponentialE;}^{- \sqrt{\frac{a}{k}} (x - 1)}}{{&ExponentialE;}^{\sqrt{\frac{a}{k}}} + 1}$

(29)

>	$pdetest (, {sys}_{3})$

$[0, 0, 0]$

(30)

The following is a textbook example where the solution is computed by constructing an infinite series and then reversed using the boundary conditions to compute a closed formula for the coefficients of the series. The three boundary conditions specify the value of $u (x, t)$ at two points (values of $x$ ) and at $t = 0$ as a generic function $f (x)$ .

>	${sys}_{4} ≔ [\frac{\partial}{\partial t} u (x, t) = k (\frac{\partial^{2}}{\partial x^{2}} u (x, t)), u (0, t) = 0, u (λ, t) = 0, u (x, 0) = f (x)]$

${sys}_{4} ≔ [\frac{\partial}{\partial t} u (x, t) = k (\frac{\partial^{2}}{\partial x^{2}} u (x, t)), u (0, t) = 0, u (λ, t) = 0, u (x, 0) = f (x)]$

(31)

>	$pdsolve ({sys}_{4})$

$u (x, t) = \sum_{n = 1}^{\infty} \frac{2 (\int_{0}^{λ} f (x) \sin (\frac{π n x}{λ}) &DifferentialD; x) \sin (\frac{π n x}{λ}) {&ExponentialE;}^{- \frac{k π^{2} n^{2} t}{λ^{2}}}}{λ}$

(32)

The following is the wave equation with conditions at $t = 0$ for $u (x, t)$ and its time derivative $\frac{\partial}{\partial t} u (x, t)$ , as well as conditions at the boundaries x = 0 and x = 1:

>	${sys}_{5} ≔ [\frac{\partial^{2}}{\partial t^{2}} u (x, t) = \frac{1}{16} (\frac{\partial^{2}}{\partial x^{2}} u (x, t)) - \frac{1}{5} (\frac{\partial}{\partial t} u (x, t)), u (0, t) = 0, u (1, t) = 0, u (x, 0) = x (1 - x), D_{2} (u) (x, 0) = x (1 - x)]$

${sys}_{5} ≔ [\frac{\partial^{2}}{\partial t^{2}} u (x, t) = \frac{\frac{\partial^{2}}{\partial x^{2}} u (x, t)}{16} - \frac{\frac{\partial}{\partial t} u (x, t)}{5}, u (0, t) = 0, u (1, t) = 0, u (x, 0) = x (1 - x), D_{2} (u) (x, 0) = x (1 - x)]$

(33)

>	$pdsolve ({sys}_{5})$

$u (x, t) = \sum_{n = 1}^{\infty} \frac{2 I (((I {(−1)}^{n} - I) \sqrt{25 π^{2} n^{2} - 4} + 22 {(−1)}^{1 + n} + 22) {&ExponentialE;}^{- \frac{(2 + I \sqrt{25 π^{2} n^{2} - 4}) t}{20}} + (I \sqrt{25 π^{2} n^{2} - 4} + 22) ({(−1)}^{n} - 1) {&ExponentialE;}^{\frac{(I \sqrt{25 π^{2} n^{2} - 4} - 2) t}{20}}) \sin (π x n)}{\sqrt{25 π^{2} n^{2} - 4} n^{3} π^{3}}$

(34)

>	$pdetest (, {sys}_{5})$

$[0, 0, 0, (\sum_{n = 1}^{\infty} - \frac{4 ({(−1)}^{n} - 1) \sin (π x n)}{n^{3} π^{3}}) + x^{2} - x, (\sum_{n = 1}^{\infty} - \frac{4 ({(−1)}^{n} - 1) \sin (π x n)}{n^{3} π^{3}}) + x^{2} - x]$

(35)

To verify the conditions at t = 0 which pdetest is unable to check, we use plots.

First, we plot together the solution at t = 0 (1000 terms) and the quadratic expression $x \cdot (1 - x)$ and see that they are identical:

>	$plot ([subs (infinity = 1000, t = 0, rhs ()), x (1 - x)], x = 0 .. 1)$

Then, we plot together the derivative of the solution at t = 0 (1000 terms) and the quadratic expression $x \cdot (1 - x)$ and see that they are identical:

>	$plot ([subs (infinity = 1000, t = 0, \frac{\partial}{\partial t} rhs ()), x (1 - x)], x = 0 .. 1)$

More Examples

$restart$

Set infolevel = 2 to receive information on the method used.

>	$infolevel [pdsolve] ≔ 2$

${infolevel}_{pdsolve} ≔ 2$

(36)

>	${pde}_{1} ≔ \frac{\partial^{2}}{\partial x^{2}} u (x, t) - \frac{\frac{\partial^{2}}{\partial t^{2}} u (x, t)}{c^{2}} = 0 &semi;$ ${iv}_{1} ≔ u (x, 0) = f (x), D_{2} (u) ()$

Maple

Erweiterungen für Maple

MapleSim

Erweiterungen für MapleSim

Systementwicklung

Beratende Dienstleistungen

Produkte zur Online-Ausbildung

Ausbildung

Industrie

Automobil und Luftfahrt

Robotik

Maschinendesign & industrielle Automation

Andere

Lösungen für die Industrie

Kaufen: Einzelheiten und Preise

Kaufen

Institutionelle Studentenlizenzierung

Elite-Wartungsprogramm

Support

Produktschulung

Online Hilfe zu Produkten

Webinare und Veranstaltungen

Veröffentlichungen

Content Hubs

Anwendungsbeispiele

Community

Über Maplesoft

Pressezentrum

User Community

Kontakt

Online Help

All Products Maple MapleSim

Maple

Leistungsfähige intuitive Mathematiksoftware

Erweiterungen für Maple

MapleSim

Modellierung auf Systemebene

Erweiterungen für MapleSim

Systementwicklung

Beratende Dienstleistungen

Produkte zur Online-Ausbildung

Ausbildung

Industrie

Automobil und Luftfahrt

Robotik

Maschinendesign & industrielle Automation

Andere

Lösungen für die Industrie

Kaufen: Einzelheiten und Preise

Kaufen

Institutionelle Studentenlizenzierung

Elite-Wartungsprogramm

Support

Produktschulung

Online Hilfe zu Produkten

Webinare und Veranstaltungen

Veröffentlichungen

Content Hubs

Anwendungsbeispiele

Community

Über Maplesoft

Pressezentrum

User Community

Kontakt

Online Help

All Products Maple MapleSim