Return to computing page for the first course APMA0330
Return to computing page for the second course APMA0340
Return to Mathematica tutorial for the first course APMA0330
Return to Mathematica tutorial for the second course APMA0340
Return to the main page for the
first course APMA0330
Return to the main page for the
second course APMA0340
Return to Part VI of the course APMA0340
Introduction to Linear Algebra with Mathematica
The order of the above differential equation is the order of the highest derivative occuring in the equation, The equation is linear if it depends linearly on u and its derivatives. If all highest-order derivatives of u occur linearly with coefficients depending only on independent variables, u, and lower-order derivatives of
u, then the equation is quasilinear. A simple example of a fist-order PDE is
\[
u_t + a(u)\,u_x = 0 ,
\]
where \( u_t \overset{\mathrm def}{=} \partial u/\partial t \) and \( u_x \overset{\mathrm def}{=} \partial u/\partial x . \) When \( a(u) \equiv a \) is a constant, the linear equation ut + aux = 0 is called the transport equation. When \( a(u) \equiv u, \) the quasilinear equation
ut + uux = 0 is called the inviscid Burger's equation, which arises in the study of a one-dimensional stream of particles or fluid having zero viscosity.
where \( a \) is a constant. We can reduce the transport equation to an exact ordinary differential equation by considering x-variable as a function of time t. So we consider the transport equation along some curve x(t). Then we will have
if we require \( {\text d}x/{\text d} t = a. \)
In this case, x = 𝑎 t + x_{0}, where x_{0} is the x-intercept of the curve. Along this curve we have \( \frac{\text d}{{\text d}t}\, u(x(t),t) = 0 , \) that is, u is a constant, which we know from the initial condition to be h(x_{0}). Thus, u(x,t) =
h(x_{0}) = h(x - at). Indeed, ih h is C¹ (= continuously differentiable), then we can check that u(x,t) = h(x - at satisfies the transport equation and prescribed initial condition. This solution corresponds to "tranporting" the initial data h(x) without change along the x-axis at a speed \( {\text d}x/{\text d} t = a. \) The lines x = 𝑎 t + x_{0} are called the characteristic curves for u,sub>t + 𝑎 u_{x} = 0. The reduction of a first order partial differential equation to an ordinary differential equation along its characteristics is called the method of characteristics.
Let us consider the quasilinear equation for a function of two variables x and y:
where the coefficient functions 𝑎, b, and c are continuous in x, y, and u. If u(x,y) is a solution of the above partial differential equation, let us consider the graph z = u(x,y). This surface has normal vector N_{0} = ⟨ -u_{x}(x_{0},y_{0}), -u_{y}(x_{0},y_{0}) ⟩ at the point (x_{0},y_{0},u(x_{0},y_{0})). However, if we let z_{0} = u(x_{0},y_{0}), then the above partial differential equation implies that that the vector
\( {\bf V} = \langle a\left( x_0 , y_0 , z_0 \right) , b\left( x_0 , y_0 , z_0 \right) , c\left( x_0 , y_0 , z_0 \right) \rangle \) is perpendicular to this normal vector, and hence must lie in the tangent plane to the graph of z = u(x,y) at the point z_{0}.
In other words, \( {\bf V}(x,y,z) = \langle a\left( x , y , z \right) , b\left( x , y , z \right) , c\left( x , y , z \right) \rangle \) defines a vector field in ℝ³, to which graphs of solutions must be tangent at each point. Surfaces that are tangent at each point to a vector field in ℝ³ are called integral surfaces of the vector field. Thus, in order to find a solution of the given partial differential equation \( a\,u_x + b\,u_y = c , \) we should try to find integral surfaces. Of course, there may be many integral surfaces of V, so we might try to be more specific and find the integral surface containing a given curve Γ ⊂ ℝ^{3}. So we are led to formulating the following:
The Cauchy problem: :
For a given curve Γ ⊂ ℝ³, can we find a solution u of the first order partial differential equation whose graph contains Γ?
⧫
In special case when Γ is the graph ⟨ x, h(x) ⟩ in the xz-plane of a function h, the Cauchy problem is just an initial value problem with the obvious interpretation of the variable y as "yime".
We constract the integral surfaces by using characteristic curves that are the integral curves of the vector field V. These curves are called the characteristic curves is they satisfy the system of equations
An infinitely differentiable function φ = φ(x, t) that vanish in a neighborhood of the boundary ∂Ω of Ωis called test function on Ω. We say that u = u(x, t) is a weak solution of equation \eqref{Eq.First.1} if
El-Shandwily, M.E., Solutions of Maxwells equations for general nonperiodic waves in lossy media, IEEE Transactions on Electromagnetic Compatibility, 1988, Volume: 30 , Issue: 4, pp. 577--582. doi: 10.1109/15.8774
Gupta, V.G. and Sharma, B., Application of Sumudu Transform in Reaction-Diffusion Systems and Nonlinear Waves, Applied Mathemathical Sciiences, 2007, Vol. 4, No. 9, pp. 435--446.
Return to Mathematica page
Return to the main page (APMA0340)
Return to the Part 1 Matrix Algebra
Return to the Part 2 Linear Systems of Ordinary Differential Equations
Return to the Part 3 Non-linear Systems of Ordinary Differential Equations
Return to the Part 4 Numerical Methods
Return to the Part 5 Fourier Series
Return to the Part 6 Partial Differential Equations
Return to the Part 7 Special Functions