15.11.2014 Views

MATLAB Mathematics - SERC - Index of

MATLAB Mathematics - SERC - Index of

MATLAB Mathematics - SERC - Index of

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

5 Differential Equations<br />

Note The Differential Equations Examples browser enables you to view the<br />

code for the PDE examples, and also run them. Click on the link to invoke the<br />

browser, or type odeexamples('pde')at the command line.<br />

Example<br />

pdex1<br />

pdex2<br />

pdex3<br />

pdex4<br />

pdex5<br />

Description<br />

Simple PDE that illustrates the straightforward formulation,<br />

computation, and plotting <strong>of</strong> the solution<br />

Problem that involves discontinuities<br />

Problem that requires computing values <strong>of</strong> the partial<br />

derivative<br />

System <strong>of</strong> two PDEs whose solution has boundary layers at<br />

both ends <strong>of</strong> the interval and changes rapidly for small t<br />

System <strong>of</strong> PDEs with step functions as initial conditions<br />

Introduction to PDE Problems<br />

pdepe solves systems <strong>of</strong> parabolic and elliptic PDEs in one spatial variable x<br />

and time t , <strong>of</strong> the form<br />

c⎛x, t, u,<br />

∂u ------ ⎞ ∂u<br />

------ = x – m ------ ∂ x m f ⎛x, t, u,<br />

∂u ------ ⎞<br />

⎝ ∂x⎠<br />

∂t ∂x ⎝<br />

⎛ ⎝ ∂x⎠⎠<br />

⎞ + s ⎛ x , t , u ,<br />

∂u ------ ⎞<br />

⎝ ∂x⎠<br />

(5-3)<br />

The PDEs hold for t 0 ≤ t ≤ t f and a ≤ x≤ b. The interval [ a,<br />

b]<br />

must be finite.<br />

m can be 0, 1, or 2, corresponding to slab, cylindrical, or spherical symmetry,<br />

respectively. If m > 0 , then a ≥ 0 must also hold.<br />

In Equation 5-3, fxtu∂u ( , , , ⁄ ∂x) is a flux term and sxtu∂u ( , , , ⁄ ∂x)<br />

is a source<br />

term. The flux term must depend on ∂u ⁄ ∂x . The coupling <strong>of</strong> the partial<br />

derivatives with respect to time is restricted to multiplication by a diagonal<br />

matrix cxtu∂u ( , , , ⁄ ∂x)<br />

. The diagonal elements <strong>of</strong> this matrix are either<br />

identically zero or positive. An element that is identically zero corresponds to<br />

an elliptic equation and otherwise to a parabolic equation. There must be at<br />

least one parabolic equation. An element <strong>of</strong> c that corresponds to a parabolic<br />

5-90

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!