Partial Differential Equations/The transport equation

From testwiki
Revision as of 00:12, 22 July 2015 by imported>Mathmensch (Exercises)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Template:Navigation

In the first chapter, we had already seen the one-dimensional transport equation. In this chapter we will see that we can quite easily generalise the solution method and the uniqueness proof we used there to multiple dimensions. Let d. The inhomogenous d-dimensional transport equation looks like this:

(t,x)×d:tu(t,x)𝐯xu(t,x)=f(t,x)

, where f:×d is a function and 𝐯d is a vector.

Solution

The following definition will become a useful shorthand notation in many occasions. Since we can use it right from the beginning of this chapter, we start with it.

Template:TextBox

Before we prove a solution formula for the transport equation, we need a theorem from analysis which will play a crucial role in the proof of the solution formula.

Template:TextBox We will omit the proof.

Template:TextBox Note that, as in chapter 1, that there are many solutions, one for each continuously differentiable g in existence.

Proof:

1.

We show that u is sufficiently often differentiable. From the chain rule follows that g(x+𝐯t) is continuously differentiable in all the directions t,x1,,xd. The existence of

xn0tf(s,x+𝐯(ts))ds,n{1,,d}

follows from the Leibniz integral rule (see exercise 1). The expression

t0tf(s,x+𝐯(ts))ds

we will later in this proof show to be equal to

f(t,x)+𝐯x0tf(s,x+𝐯(ts))ds,

which exists because

x0tf(s,x+𝐯(ts))ds

just consists of the derivatives

xn0tf(s,x+𝐯(ts))ds,n{1,,d}

2.

We show that

(t,x)×d:tu(t,x)𝐯xu(t,x)=f(t,x)

in three substeps.

2.1

We show that

tg(x+𝐯t)𝐯xg(x+𝐯t)=0(*)

This is left to the reader as an exercise in the application of the multi-dimensional chain rule (see exercise 2).

2.2

We show that

t0tf(s,x+𝐯(ts))ds𝐯x0tf(s,x+𝐯(ts))ds=f(t,x)(**)

We choose

F(t,x):=0tf(s,x𝐯s)ds

so that we have

F(t,x+𝐯t)=0tf(s,x+𝐯(ts))ds

By the multi-dimensional chain rule, we obtain

ddtF(t,x+𝐯t)=(tF(t,x+𝐯t)x1F(t,x+𝐯t)xdF(t,x+𝐯t))(1𝐯)=tF(t,x+𝐯t)+𝐯xF(t,x+𝐯t)

But on the one hand, we have by the fundamental theorem of calculus, that tF(t,x)=f(t,x𝐯t) and therefore

tF(t,x+𝐯t)=f(t,x)

and on the other hand

xnF(t,x+𝐯t)=xn0tf(s,x+𝐯(ts))ds

, seeing that the differential quotient of the definition of xn is equal for both sides. And since on the third hand

ddtF(t,x+𝐯t)=t0tf(s,x+𝐯(ts))ds

, the second part of the second part of the proof is finished.

2.3

We add (*) and (**) together, use the linearity of derivatives and see that the equation is satisfied.

Initial value problem

Template:TextBox

Proof:

Quite easily, u(0,x)=g(x+𝐯0)+00f(s,x+𝐯(ts))ds=g(x). Therefore, and due to theorem 2.3, u is a solution to the initial value problem of the transport equation. So we proceed to show uniqueness.

Assume that v is an arbitrary other solution. We show that v=u, thereby excluding the possibility of a different solution.

We define w:=uv. Then

(t,x)×d:tw(t,x)𝐯xw(t,x)=(tu(t,x)𝐯xu(t,x))(tv(t,x)𝐯xv(t,x))=f(t,x)f(t,x)=0(*)xd:w(0,x)=u(0,x)v(0,x)=g(x)g(x)=0(**)

Analogous to the proof of uniqueness of solutions for the one-dimensional homogenous initial value problem of the transport equation in the first chapter, we define for arbitrary (t,x)×d,

μ(t,x)(ξ):=w(tξ,x+𝐯ξ)

Using the multi-dimensional chain rule, we calculate μ(t,x)(ξ):

μ(t,x)(ξ):=ddξw(tξ,x+𝐯ξ) by defs. of the  symbol and μ=(tw(tξ,x+𝐯ξ)x1w(tξ,x+𝐯ξ)xdw(tξ,x+𝐯ξ))(1𝐯)chain rule=tw(tξ,x+𝐯ξ)+𝐯xw(tξ,x+𝐯ξ)=0(*)

Therefore, for all (t,x)×d μ(t,x)(ξ) is constant, and thus

(t,x)×d:w(t,x)=μ(t,x)(0)=μ(t,x)(t)=w(0,x+𝐯t)=(**)0

, which shows that w=uv=0 and thus u=v.

Exercises

  1. Let f𝒞1(×d) and 𝐯d. Using Leibniz' integral rule, show that for all n{1,,d} the derivative

    xn0tf(s,x+𝐯(ts))ds

    is equal to

    0txnf(s,x+𝐯(ts))ds

    and therefore exists.

  2. Let g𝒞1(d) and 𝐯d. Calculate tg(x+𝐯t).
  3. Find the unique solution to the initial value problem

    {(t,x)×3:tu(t,x)(234)xu(t,x)=t5+x16+x27+x38x3:u(0,x)=x19+x210+x311.

Sources

Template:Navigation