Lect 2 PDF
Lect 2 PDF
-u
qq
=h
1
(u
, u
q
, u, , q)
u
q
=h
2
(u
, u
q
, u, , q)
Parabolic PDE:
u
=h
3
(u
, u
q
, u, , q)
or u
qq
=h
4
(u
, u
q
, u, , q)
Elliptic PDE:
u
+u
qq
=h
5
(u
, u
q
, u, , q)
We are assured of a nonsingular mapping, provided that
the Jacobian of the transformation is
nonzero (Taylor, 1955).
Any real nonsingular transformation does not change the
type of PDE.
Hyperbolic PDEs
Two distinct families of characteristics
exist.
Second-order wave equation: u
tt
=c
2
u
xx
-<x<+, u(x,0)=f(x), u
t
(x,0)=g(x)
Solution
Transformation to characteristic coordinates
permits simple integration of the wave equation:
u
q
=0 where =x+ct, q=x-ct.
First term: propagation of the initial data along
the characteristics
Second term: effect of the data within the closed
interval at t=0
Right running characteristic: +(1/c)
Left running characteristic: -(1/c)
Fundamental Property of
Hyperbolic PDE
Limited domain of dependence is bounded
by the characteristics that pass through
the point (x
0
, t
0
).
Any disturbance that occurs outside of this
interval can never influence the solution at
(x
0
, t
0
).
Parabolic PDEs
Parabolic PDEs associated with diffusion
processes.
Solution of parabolic equation at time t
1
depends
upon the entire physical domain (tst
1
), including
any side BCs.
Rayleigh problem: 2-D flat plate impulsively accelerated to
a velocity U from rest.
Similarity solution
(Hansen, 1964)
Growth of this layer is controlled by v and t.
Velocity change in the layer is induced by diffusion of the
plate velocity into the initially undisturbed fluid.
Error function
Elliptic PDEs
The discriminant is negative (b
2
-4ac<0) and
the characteristic differential equation has
no real solution.
It means that there is no real characteristics
exist and it is meaningless to find them.
The solution of elliptic PDE depends upon
conditions on the entire boundary of the
closed domain.
Well-Posed Problem
A unique solution to PDEs cannot be
obtained if the initial data or boundary
conditions is improper used.
A problem involving a PDE to be well-
posed, the solution to the problem must
exist, must be unique, and must depend
continuously upon the initial or boundary
data (Hadamard, 1952).
Different Types of BCs
Dirichlet BCs [u|
B
=f] + PDEs|
D
Dirichlet Problem
Neumann BCs [cu/cn|
B
=f] +PDEs|
D
Neumann Problem
Robin!s BCs: [(cu/cn+ku)|
B
=f] + PDEs|
D
Robin!s Problem
u
Systems of Eqns.
In applying numerical methods to physical
problems, systems of eqns. are frequently
encountered.
The process is governed by a higher-order PDE,
the PDE can usually be converted to a system of
1st order equations.
Wave equation: u
tt
=c
2
u
xx
cv/ct=ccw/cx, cw/ct=ccv/cx where v=cu/ct,
w=ccu/cx
Laplace equation: u
xx
+u
yy
=0
cu/cx=cv/cy, cu/cy=-cv/cx
System of 1st Order Eqns.
Consider curves C on which all but the highest
order derivatives are specified and inquire about
conditions that will indicate that the highest
derivatives are not uniquely determined.
We let a parameter t vary along curves C and use
the Chain rule to write the following Eqns.
Writing the four equations in matrix form:
A unique solution for the first derivatives of u(x, y)
and v(x, y) does not exist if the determinant of
the coefficient matrix is zero.
Classification of 1st order system:
D>0: Hyperbolic system, D=0: Parabolic system,
D<0: Elliptic system.
In 1st order system, if the roots of the
characteristic eqns. are contain both real
and complex parts, then the system is
mixed and may exhibit hyperbolic, parabolic,
and elliptic behavior.
The classification of systems of 2nd order
PDEs is very complex.
It is difficult to determine the mathematical
behavior of these systems except for simple
cases.
Other differential equations of interest
2
nd
equations: Wave equation, Heat
equation, Laplaces equation.
A number of other very important
equations should be mentioned, since tyey
govern common physical phenomena or
they are used as simple models for more
complex problems. In many cases, exact
analytical solutions for these equations
exist.