Solution to Partial Differential Equations: Methods and Applications (Robert McOwen) Section 2.3
- For complex-valued functions u and v on Ω, let ⟨u,v⟩=∫Ωuˉvdx, where ˉv denotes the complex conjugate of v. If the coefficients aα(x) in (32a) are complex-valued functions, define the adjoint L∗ so that ⟨Lu,v⟩=⟨u,L∗v⟩ for all u∈Cm(Ω), v∈Cm0(Ω).
- For (complex) vector-valued functions →u(x)=(u1(x),…,uN(x)) on Ω, let ⟨→u,→v⟩=∫Ω(u1ˉv1+⋯+uNˉvN)dx. If the coefficients aα(x) in (32a) are N×N matrix-valued functions, then (32a) defines a system of mth-order operators. Define the adjoint L∗ so that ⟨L→u,→v⟩=⟨→u,L∗→v⟩ for all →u∈Cm(Ω,RN), →v∈Cm0(Ω,RN).
- From (32a), we find that
⟨Lu,v⟩=∫Ω(Lu)(ˉv)dx=∫Ω(∑|α|≤maα(x)Dαu)ˉvdx=∑|α|≤m∫Ω(Dαu)(aα(x)ˉv)dx=∑|α|≤m(−1)|α|∫Ω(u)(Dα(aα(x)ˉv))dx=∫Ω(u)[∑|α|≤m(−1)|α|Dα(aα(x)ˉv)]dx:=⟨u,L∗v⟩,
where the adjoint operator L∗ is given byL∗v:=¯∑|α|≤m(−1)|α|Dα(aα(x)ˉv)=∑|α|≤m(−1)|α|Dα(¯aα(x)v)for v∈Cm0(Ω).
- Suppose aα(x)=[aijα(x)]1≤i,j≤N is a N×N matrix-valued function. Then we have
⟨L→u,→v⟩=∫Ω[(N∑j=1∑|α|≤ma1jα(x)Dαuj)ˉv1+⋯+(N∑j=1∑|α|≤maNjα(x)Dαuj)ˉvN]dx=∫Ω(N∑j=1∑|α|≤ma1jα(x)Dαuj)ˉv1dx+⋯+∫Ω(N∑j=1∑|α|≤maNjα(x)Dαuj)ˉvNdx=N∑j=1∑|α|≤m[∫Ω(Dαuj)(a1jα(x)ˉv1)dx+⋯+∫Ω(Dαuj)(aNjα(x)ˉvN)dx]=N∑j=1∑|α|≤m(−1)|α|[∫Ω(uj)Dα(a1jα(x)ˉv1)dx+⋯+∫Ω(uj)Dα(aNjα(x)ˉvN)dx]=N∑j=1∫Ω(uj)∑|α|≤m(−1)|α|Dα(a1jα(x)ˉv1+⋯+aNjα(x)ˉvN)dx=⟨→u,L∗→v⟩,
where the adjoint operator L∗ is given byL∗→v=¯∑|α|≤m(−1)|α|Dα(aTα(x)ˉv)=∑|α|≤m(−1)|α|Dα(¯aTα(x)v)for v∈Cm0(Ω,RN).
Here AT denotes the tranpose of the matrix A. - Let Lu=uμη as in Example 1. Derive the transmission condition (44) along the characteristic μ=0.
- Consider the first-order equation ut+cux=0.
- If f∈C(R), show that u(x,t)=f(x−ct) is a weak solution.
- Can you find any discontinuous weak solutions?
- Is there a transmission condition for a weak solution with jump discontinuity along the characteristic x=ct?
- To show u(x,t)=f(x−ct) is a weak solution, it suffices to verify that ∫R2f(x−ct)(L′v)(x,t)dxdt=0 for v∈C10(R2), where L′=−∂t−c∂x is the adjoint operator L=∂t+c∂x. Let s=x−ct, t=t and ˜v(s,t)=v(s+ct,t) for s,t∈R2. Clearly, ˜v∈C10(R2) and satisfies ˜vt(s,t)=vt(s+ct,t)+cvx(s+ct,t). Thus, the improper double integral becomes
∫R2f(x−ct)(L′v)(x,t)dxdt=−∫∞−∞∫∞−∞f(x−ct)(vt(x,t)+cvx(x,t))dxdt=−∫∞−∞∫∞−∞f(s)˜vt(s,t)dsdt=−∫∞−∞f(s)∫∞−∞˜vt(s,t)dtds=−∫∞−∞f(s)⋅0ds=0for v∈C10(R2).
Here the Jacobian is 1. Therefore, u(x,t)=f(x−ct) is a weak solution for ut+cux=0. - When f is discontinuous but Lebesgue integrable, u(x,t)=f(x−ct) is still a weak solution because the computation in part (a) also holds true.
- Suppose u is a weak solution with jump discontinuity along the characteristic x=ct. Let ˜u(s,t)=u(s+ct,t). Then ˜u is a weak solution of ˜L˜u:=˜ut(s,t)=0 with jump discontinuity along the line s=0, where ˜L=∂t. Clearly, the adjoint operator ˜L′ of ˜L is −∂t. Let R2±={(s,t)∈R2:±s>0,t∈R} and ˜u±=˜u|R2±. Then for any test function v∈C10(R2), we have
0=∫R2˜u(s,t)˜L′v(s,t)dsdt=∫R2˜u(s,t)⋅(−∂tv(s,t))dsdt=−∫∞0∫∞−∞˜u+(s,t)vt(s,t)dtds−∫0−∞∫∞−∞˜u−(s,t)vt(s,t)dtds=∫∞0[−˜u+(s,t)v(s,t)|t=∞t=−∞+∫∞−∞˜u+t(s,t)v(s,t)dt]ds+∫0−∞[−˜u−(s,t)v(s,t)|t=∞t=−∞+∫∞−∞˜u−t(s,t)v(s,t)dt]ds=∫∞0∫∞−∞˜u+t(s,t)v(s,t)dtds+∫0−∞∫∞−∞˜u−t(s,t)v(s,t)dtds=0.
Thus, there is no transmission condition for a weak solution with jump discontinuity along the characteristic x=ct. - If f∈L1loc(Ω), define ⟨Ff,v⟩≡∫Ωf(x)v(x)dx. Show that Ff is a distribution in Ω.
- If ξ1,…,ξk are points in Ω and a1,…,ak are real numbers, show that F≡a1δξ1+⋯+akδξk is a distribution in Ω.
- For Ω=Rn, find an infinite sequence {ξk}∞k=1 for which f≡a1δξ1+⋯ is a distribution. (This gives an example of a distribution that cannot be realized as a distributional derivative of an integrable function.)
- Let a,b∈R and v,w∈D(Ω). It is easy to verify that
F(av+bw)=⟨F,av+bw⟩=⟨k∑i=1aiδξi,av+bw⟩=∫Ωk∑i=1aiδξi(x)[av(x)+bw(x)]dx=k∑i=1ai∫Ωδξi(x)[av(x)+bw(x)]dx=k∑i=1ai[av(ξi)+bw(ξi)]=ak∑i=1aiv(ξi)+bk∑i=1aiw(ξi)=ak∑i=1ai∫Ωδξi(x)v(x)dx+bk∑i=1ai∫Ωδξi(x)w(x)dx=a∫Ωk∑i=1aiδξi(x)[v(x)]dx+b∫Ωk∑i=1aiδξi(x)[w(x)]dx=a⟨k∑i=1aiδξi,v⟩+ba⟨k∑i=1aiδξi,w⟩=a⟨F,v⟩+b⟨F,w⟩=aF(v)+bF(w),
which shows F is a linear mapping on D(Ω). Moreover, let {vj}∞j=1 be a sequence of functions in D(Ω) with supp(vj)⊆K for some compact set K⊊Ω (independent of j) and Dαvj→0 uniformly in K. Since K∩{ξ1,…,ξk} is finite, we may write it as {ξj1,…,ξjr} for r∈N. Set A=max1≤i≤r|aji|. Then for any ε>0, there exists N∈N such that |vj(x)|≤ε/(rA) for j>N and x∈K. For j>N, we have|F(vj)−0|=|⟨F,vj⟩|=|∫ΩF(x)vj(x)dx|=|∫Ωk∑i=1aiδξi(x)vj(x)dx|≤r∑i=1|aji||∫Kδξji(x)vj(x)dx|=r∑i=1|aji||vj(ξji)|≤r∑i=1A⋅εrA=ε.
This implies limj→∞F(vj)=0. Therefore, F is continuous on D(Ω) and F is a distribution in Ω. - Choose ξj=je1:=(j,0,…,0) for j∈N, and then we prove F=∑∞i=1aiδξi is a distribution in Rn. Let a,b∈R and v,w∈D(Rn). Then there exists a compact subset K⊊Rn such that supp(v)⊆K and supp(w)⊆K. Since K is bounded, K∩{ξj}∞j=1 is finite and K∩{ξj}∞j=1={ξj1,…,ξjr} for r∈N. Then one may check that
F(av+bw)=⟨F,av+bw⟩=⟨∞∑j=1ajδξj,av+bw⟩=∫Rn∞∑j=1ajδξj(x)[av(x)+bw(x)]dx=∫Kr∑i=1ajiδξji(x)[av(x)+bw(x)]dx=r∑i=1aji∫Kδξji(x)[av(x)+bw(x)]dx=r∑i=1aji[av(ξji)+bw(ξji)]=ar∑i=1ajiv(ξji)+br∑i=1ajiw(ξji)=ar∑i=1aji∫Kδξji(x)v(x)dx+br∑i=1aji∫Kδξji(x)w(x)dx=a∫Kr∑i=1ajiδξji(x)v(x)dx+b∫Kr∑i=1ajiδξji(x)w(x)dx=a∫Rn∞∑j=1ajδξj(x)v(x)dx+b∫Rn∞∑j=1ajδξj(x)w(x)dx=a⟨∞∑j=1ajδξj,v⟩+b⟨∞∑j=1ajδξj,w⟩=a⟨F,v⟩+b⟨F,w⟩=aF(v)+bF(w),
which shows F is a linear mapping on D(Rn). Moreover, let {vk}∞k=1 be a sequence of functions in D(Rn) with supp(vk)⊆K for some compact set K⊊Rn (independent of j) and Dαvk→0 uniformly in K. Since K is bounded, K∩{ξj}∞j=1 is finite and K∩{ξj}∞j=1={ξji}ri=1 for r∈N. Set A=max1≤i≤r|aji|. Then for any ε>0, there exists N∈N such that |vk(x)|<ε/(rA) for k>N and x∈K. For j>N, we have|F(vk)−0|=|⟨F,vk⟩|=|∫RnF(x)vk(x)dx|=|∫Rn∞∑j=1ajδξj(x)vk(x)dx|=|∫Kr∑i=1ajiδξji(x)vk(x)dx|=|r∑i=1ajivk(ξji)|≤r∑i=1|aji||vk(ξji)|≤r∑i=1A⋅εrA=ε.
This implies limk→∞F(vk)=0. Therefore, F is continuous on D(Rn) and F is a distribution in Rn. - Prove the Lemma on p. 67 when f∈C(Rn) has a compact support.
- Find an example to show that the result fails to hold if g∈L1loc(Rn) is replaced by g∈D′(Rn).
- Suppose that K=supp(f) is a compact subset of Rn. Let x∈Rn and y∈ˉB1(x) arbitrarily and define h(z)=f(x−z)−f(y−z) for z∈Rn. Clearly, h is continuous function compact support, and
supp(h)⊆(x−K)∪(ˉB1(x)−K):=˜K,
wherex−K={x−z∈Rn:z∈K} and ˉB1(x)−K={y−z∈Rn:y∈ˉB1(x),z∈K}.
On the other hand, f is uniformly continuous in Rn. Then for any ε>0, there exists δ=δ(ε)>0 such that |f(x−z)−f(y−z)|<ε(∫˜K|g(z)|dz)−1 whenever |x−y|<δ and z∈Rn. Thus, we find that|(f∗g)(x)−(f∗g)(y)|=|∫Rn[f(x−z)−f(y−z)]g(z)dz|=|∫˜K[f(x−z)−f(y−z)]g(z)dz|≤∫˜K|f(x−z)−f(y−z)||g(z)|dz<ε.
This shows limy→x(f∗g)(y)=(f∗g)(x), which means f∗g∈C(Rn). - For n=1, let g=δ′∈D′(R) and f(x)=|x| for x∈R. Then one may check that
(f∗g)(x)=∫Rf(x−y)g(y)dy=∫R|x−y|δ′(y)dy=−∫Rddy|x−y|δ(y)dy=∫Rsgn(x−y)δ(y)dy=sgn(x),
which is not continuous on R. - Let Γ be a hypersurface (such as a sphere) in Rn, let a(z) be a continuous function of z∈Γ, and let dz denote the surface measure on Γ.
- Show that ⟨F,v⟩=∫Γa(z)v(z)dz is a distribution in Rn. (It is natural to denote this distribution by F=aδΓ.)
- Suppose we can choose a unit normal ν along Γ. Formulate a definition for the "conormal distribution" a∂νδΓ in Rn.
- Note that Γ is a hypersurface in Rn so Γ is a closed set in Rn. For c1,c2∈R and v1,v2∈D(Rn), we have
⟨F,c1v1+c2v2⟩=∫Γa(z)[c1v1(z)+c2v2(z)]dz=c1∫Γa(z)v1(z)dz+c2∫Γa(z)v2(z)dz=c1⟨F,v1⟩+c2⟨F,v2⟩,
which shows F is a linear mapping on D(Rn). Here we have used the fact that v1,v2∈C∞(Γ)⊊D(Rn). Moreoever, let {vj}∞j=1 be a sequence of functions in D(Rn) with supp(vj)⊆K for some compact set K⊊Rn (independent of j) and Dαvj→0 uniformly in K. Since a is continuous on the compact set Γ∩K, the maximum of |a| over Γ∩K is finite and can be denoted by A:=maxz∈Γ∩Ka(z). For any ε>0, there exists N∈N such that|vj(x)|≤ε[A∫Γ∩Kdz]−1 for j>N and x∈K.
For j>N, we have|⟨F,vj⟩−0|=|∫Γa(z)vj(z)dz|=|∫Γ∩Ka(z)vj(z)dz|≤∫Γ∩K|a(z)||vj(z)|dz≤∫Γ∩KA⋅ε[A∫Γ∩Kdz]−1dz=ε.
This implies limj→∞F(vj)=0. Therefore, F is continuous on D(Rn) and F is a distribution in Rn. - Suppose a∈C1(Γ). Then we define the conformal distribution a∂νδΓ by
⟨a∂νδΓ,v⟩=−∫Γ∇(a(z)v(z))⋅ν(z)dzfor v∈D(Rn).
For c1,c2∈R and v1,v2∈D(Rn), we have⟨a∂νδΓ,c1v1+c2v2⟩=−∫Γ∇(a(z)[c1v1(z)+c2v2(z)])⋅ν(z)dz=−c1∫Γ∇(a(z)v1(z))⋅ν(z)dz−c2∫Γ∇(a(z)v2(z))⋅ν(z)dz=c1⟨a∂νδΓ,v1⟩+c2⟨a∂νδΓ,v2⟩,
which shows ⟨a∂νδΓ,⋅⟩ is a linear mapping on D(Rn). Here we have used the fact that v1,v2∈C∞(Γ)⊊D(Rn). Moreover, let {vj}∞j=1 be a sequence of functions in D(Rn) with supp(vj)⊆K for some compact set K⊊Rn (independent of j) and Dαvj→0 uniformly in K. Since a is continuously differentiable on the compact set Γ∩K, ‖a‖C1(Γ∩K):=maxz∈Γ∩K|a(z)|+maxz∈Γ∩K|∇a(z)| is finite. For any ε>0, there exists N∈N such thatvj(x)≤ε2[‖a‖C1(Γ∩K)∫Γ∩Kdz]−1,|∇vj(x)|≤ε2[‖a‖C1(Γ∩K)∫Γ∩Kdz]−1for j>N and x∈K.
For j>N, we have|⟨a∂νδΓ,vj⟩|=|∫Γ∇(a(z)vj(z))⋅ν(z)dz|=|∫Γ[vj(z)∇a(z)+a(z)∇vj(z)]⋅ν(z)dz|=|∫Γ∩K[vj(z)∇a(z)+a(z)∇vj(z)]⋅ν(z)dz|≤∫Γ∩K(|vj(z)||∇a(z)|+|a(z)||∇vj(z)|)dz≤∫Γ∩K‖a‖C1(Γ∩K)(|vj(z)|+|∇vj(z)|)dz≤‖a‖C1(Γ∩K)∫Γ∩Kε[‖a‖C1(Γ∩K)∫Γ∩Kdz]−1dz=ε.
This implies limj→∞⟨a∂νδΓ,vj⟩=0. Therefore, ⟨a∂νδΓ,⋅⟩ is continuous on D(Rn) and F is a distribution in Rn. - Let
fn(x)={n2for −1n<x<1n;0for |x|≥1n.
Show that fn(x)→δ(x) as distributions on R. - If fn(x) and f(x) are integrable functions such that for any compact set K⊂Ω we have ∫K|fn(x)−f(x)|dx→0 as n→∞, then fn→f as distributions.
- Let a∈R, a≠0.
- Find a fundamental solution for L=d/dx−a on R (i.e., solve dF/dx−aF=δ).
- Show that a fundamental solution for L=d2/dx2−a2=(d/dx+a)(d/dx−a) on R is given by
F(x)={a−1sinhaxif x>0;0if x<0.
- Multiplying equation by e−ax, it gives
ddx(e−axF)=e−axdFdx−ae−axF=e−axδ(x)=δ(x),
which impliese−axF(x)=H(x)={1if x>0;0if x<0.
Hence F(x)=eaxH(x)={eaxif x>0;0if x<0.
Note that the adjoint operator L′ of L is −d/dx−a. It is easy to verify that⟨LF,v⟩=⟨F,L′v⟩=∫RF(x)(−dvdx−av)dx=−∫∞0eaxdvdxdx−a∫∞0eaxv(x)dx=−eaxv(x)|∞0+a∫∞0eaxv(x)dx−a∫∞0eaxv(x)dx=v(0)=⟨δ,v⟩,
which shows LF=δ. - Note that L is self-adjoint, i.e., L=L′. Then for v∈D(R), it is easy to verify that
⟨LF,v⟩=⟨F,L′v⟩=∫RF(x)(d2vdx2−a2v(x))dx=∫∞0a−1sinh(ax)(d2vdx2−a2v(x))dx=a−1sinh(ax)dvdx|∞0−∫∞0cosh(ax)dvdxdx−∫∞0asinh(ax)v(x)dx=−cosh(ax)v(x)|∞0+∫∞0asinh(ax)v(x)dx−∫∞0asinh(ax)v(x)dx=v(0)=⟨δ,v⟩,
which implies LF=δ. - Regularity of u defined by (60):
- If f(x)∈L1(R) has a compact support, show that (60) defines a continuous weak solution of u″=f(x).
- If, in addition to (a), f(x) is a bounded function on R, show u∈C1(R), and u′(x)=12(∫x−∞f(y)dy−∫∞xf(y)dy).
- If, in addition to (a), f(x) is continuous, show that u is in fact a classical solution: u∈C2(R) with u″=f.
- Since f∈L1(R) has a compact support, f∈L1loc(R). On the other hand, g(x)=12|x|∈C(R). Thus, by Lemma at Page 67, u(x)=(g∗f)(x)=12∫R|x−y|f(y)dy is continuous on R.
- From (60), we have
u(x)=12∫x−∞(x−y)f(y)dy+12∫∞x(y−x)f(y)dy.
A direct computation givesu(x+h)−u(x)=12∫x+h−∞(x+h−y)f(y)dy−12∫x−∞(x−y)f(y)dy+12∫∞x+h(y−x−h)f(y)dy−12∫∞x(y−x)f(y)dy=h2∫x+h−∞f(y)dy+12∫x+hx(x−y)f(y)dy−h2∫∞x+hf(y)dy−12∫x+hx(y−x)f(y)dy.
Since f is bounded, there exists M>0 such that |f(y)|≤M for y∈R. Then it is easy to observe that|∫x+hx±(x−y)f(y)dy|≤|∫x+hx(y−x)Mdy|=Mh22,
which implies limh→01h∫x+hx±(x−y)f(y)dy=0 by the squeeze theorem. Therefore, by the definition, we findu′(x)=limh→0u(x+h)−u(x)h=12∫x+h−∞f(y)dy−12∫∞x+hf(y)dy+12h∫x+hx(x−y)f(y)dy−12h∫x+hx(y−x)f(y)dy=12∫x−∞f(y)dy−12∫∞xf(y)dy.
Here we have used the fact that functions ∫∞xf(y)dy and ∫x−∞f(y)dy are continuous in x because f∈L1(R) is bounded with compact support. This also shows u′∈C(R), and hence u∈C1(R). - When continuous function f∈L1(R) has a compact support, f is bounded on R. Thus, by (b), u∈C1(R). Moreover, by the fundamental theorem of calculus, functions ∫x−∞f(y)dy and ∫∞xf(y)dy are differentiable in x, and then u′ is also differentiable. By the fundamental theorem of calculus, we arrive at u″(x)=12f(x)−12⋅(−f(x))=f(x) for x∈R. Since f∈C(R), u″∈C(R) and then u∈C2(R).
- Use the fundamental solution (59) to solve the initial u″=f(x) for x>0 with u(0)=u0 and u′(0)=u′0, where f∈C∞([0,∞)) and f=O(|x|−2−ϵ) as |x|→∞.
- Do the same for the boundary value problem u″=f(x) for 0<x<ℓ with u(0)=0=u(ℓ).
- Extend ˜f by ˜f(x)=f(x) if x>0; 0 if x≤0. Then we consider the equation ˜u″(x)=˜f(x) for x∈R. Plugging (59) into (56), we obtain the particular solution
up(x)=(F∗f)(x)=∫RF(x−y)˜f(y)dy=12∫R|x−y|˜f(y)dy=12∫∞0|x−y|f(y)dyfor x∈R.
To fulfill the initial condition, we need to find a function uh satisfyingu″h(x)=0,uh(0)=u0−up(0)=u0−12∫∞0yf(y)dy,u′h(0)=u′0−u′p(0)=u′0+12∫∞0f(y)dy.
Here we have used the result of Exercise 11(b). It is easy to obtainuh(x)=uh(0)+u′h(0)x=u0−12∫∞0yf(y)dy+(u′0+12∫∞0f(y)dy)xfor x∈R.
Thus, we find˜u(x)=uh(x)+up(x)=u0−12∫∞0yf(y)dy+(u′0+12∫∞0f(y)dy)x+12∫∞0|x−y|f(y)dyfor x∈R.
In particular, the solution is u(x)=˜u|x≥0(x) for x≥0.
It is remark that the solution u can be represented asu(x)=u(0)+∫x0u′(t)dt=u0+∫x0[u′(0)+∫t0u″(s)ds]dt=u0+u′0x+∫x0∫t0f(s)dsdt=u0+u′0x+∫x0∫xsf(s)dtds=u0+u′0x+∫x0(x−s)f(s)dsfor x≥0..
- Extend ˜f by ˜f(x)=f(x) if 0<x<ℓ; 0 if x≤0 or x≥ℓ. Then we consider the equation ˜u″(x)=˜f(x) for x∈R. Plugging (59) into (56), we obtain the particular solution
up(x)=(F∗f)(x)=∫RF(x−y)˜f(y)=12∫R|x−y|˜f(y)dy=12∫ℓ0|x−y|f(y)dyfor x∈R.
To fulfill the boundary condition, we need to find a function uh satisfyingu″h(x)=0,uh(0)=−up(0)=−12∫ℓ0yf(y)dy,uh(ℓ)=−up(ℓ)=−12∫ℓ0(ℓ−y)f(y)dy.
It is easy to obtainuh(x)=−(12ℓ∫ℓ0(ℓ−y)f(y)dy)x−(12ℓ∫ℓ0yf(y)dy)(ℓ−x)=xℓ∫ℓ0(2y−ℓ)f(y)dy−12∫ℓ0yf(y)dyfor x∈R.
Thus, we find˜u(x)=uh(x)+up(x)=xℓ∫ℓ0(2y−ℓ)f(y)dy−12∫ℓ0yf(y)dy+12∫ℓ0|x−y|f(y)dyfor x∈R.
In particular, the solution is u(x)=˜u|0≤x≤ℓ(x) for 0≤x≤ℓ.
Note that u(ℓ)=0=u(0), which implies ∫ℓ0u′(t)dt=0, and hence0=∫ℓ0[u′(0)+∫t0u″(s)ds]dt=ℓu′(0)+∫ℓ0∫t0f(s)dsdt.
Thus, u′(0)=−1ℓ∫ℓ0∫t0f(s)dsdt. Therefore, the solution can be represented asu(x)=u(0)+∫x0u′(t)dt=∫x0[u′(0)+∫t0u″(s)ds]dt=−(1ℓ∫ℓ0∫t0f(s)dsdt)x+∫x0∫t0f(s)dsdt=−xℓ∫ℓ0(ℓ−s)f(s)ds+∫x0(x−s)f(s)dsfor 0≤x≤ℓ.
- Using δ(μ,η)=δ(μ)δ(η), show that each of the following functions is a fundamental solution of L=∂2/∂μ∂η:
F1(μ,η)=H(μ)H(η),F2(μ,η)=−H(μ)H(−η),F3(μ,η)=−H(−μ)H(η),F4(μ,η)=H(−μ)H(−η).
- Use part (a) with the change of variables μ=x+ct, η=x−ct to obtain four distinct fundamental solutions for the one-dimensional wave operator, each having support in one of the wedges determined by the lines x=±ct.
- Using δ(μ,η)=δ(μ)δ(η), show that each of the following functions is a fundamental solution of L=∂2/∂μ∂η:
- We firstly note that L is self-adjoint, i.e., L=L′. For v∈D(R2), it is easy to verify that
⟨LF1,v⟩=⟨F1,L′v⟩=∫R2H(μ)H(η)vμη(μ,η)dμdη=∫∞0∫∞0vμη(μ,η)dμdη=∫∞0vη(μ,η)|μ=∞μ=0dη=−∫∞0vη(0,η)dη=−v(0,η)|∞η=0=v(0,0)=⟨δ,v⟩,⟨LF2,v⟩=⟨F2,L′v⟩=−∫R2H(μ)H(−η)vμη(μ,η)dμdη=−∫0−∞∫∞0vμη(μ,η)dμdη=−∫0−∞vη(μ,η)|μ=∞μ=0dη=∫0−∞vη(0,η)dη=v(0,η)|η=0η=−∞=v(0,0)=⟨δ,v⟩,⟨LF3,v⟩=⟨F3,L′v⟩=−∫R2H(−μ)H(η)vμη(μ,η)dμdη=−∫∞0∫0−∞vμη(μ,η)dμdη=−∫∞0vη(μ,η)|μ=0μ=−∞dη=−∫∞0vη(0,η)dη=−v(0,η)|η=∞η=0=v(0,0)=⟨δ,v⟩⟨LF4,v⟩=⟨F4,L′v⟩=∫R2H(−μ)H(−η)vμη(μ,η)dμdη==∫0−∞∫0−∞vμη(μ,η)dμdη=∫0−∞vη(μ,η)|μ=0μ=−∞dη=∫0−∞vη(0,η)dη=v(0,η)|η=0η=−∞=v(0,0)=⟨δ,v⟩.
This shows Fi is a fundamental solution of LFi=δ for i=1,2,3,4. - By (a) with the change of variables μ=x+ct and η=x−ct, the fundamental solutions for −4c2˜uμη(μ,η)=utt(x,t)−c2uxx(x,t)=δ are
U1(x,t)=−(2c)−1H(x+ct)H(x−ct),U2(x,t)=(2c)−1H(x+ct)H(−x+ct),U3(x,t)=(2c)−1H(−x−ct)H(x−ct),U4(x,t)=−(2c)−1H(−x−ct)H(−x+ct).
- If K(x,t) is a distribution in Rn depending on the parameter t, then we say K depends continuously on t if, for every v∈C∞0(Rn), f(t)≡⟨K(⋅,t),v⟩=∫K(x,t)v(x)dx defines a continuous function of t. Use this idea to define K is continuously differentiable in t.
- According to Section 1.2, a weak solution of uy+(G(u))x=0 in a rectangle Ω=(a,b)×(c,d) satisfies
ddy∫x2x1u(x,y)dx+G(u(x,y))|x=x2x=x1=0
for every [x1,x2]⊂(a,b) and c<y<d. Show that this implies∫Ωuvy+G(u)vxdxdy=0for every v∈C10(Ω)
(i.e., the notion of weak solution of Section 1.2 is consistent with that of Section 2.3). - The mth-order operator (32a) is elliptic at x if its principal symbol (32c) has no nonzero characteristics [i.e., σL(x,ξ)≠0 for all ξ∈Rn∖{0}]. In this case, show that m must be an even integer.
Solution
Solution
Suppose u is a weak solution of Lu=0, i.e., ∫Ωuμηvdμdη=∫Ω(Lu)vdμdη=0 for v∈C20(Ω). Recall that the operator L is self-adjoint, i.e., L=L′. From (40), we have0=∫Ωu(L′v)dμdη=∫Ω+u+(L′v)dμdη+∫Ω−u−(L′v)dμdη=∫Ω+u+vμηdμdη+∫Ω−u−vμηdμdη=∫∞0∫∞−∞u+vμηdηdμ+∫0−∞∫∞−∞u−vμηdηdμ=−∫∞0∫∞−∞u+ηvμdηdμ−∫0−∞∫∞−∞u−ηvμdηdμ=−∫∞−∞∫∞0u+ηvμdμdη−∫∞−∞∫0−∞u−ηvμdμdη=−∫∞−∞[u+ηv|μ=∞μ=0−∫∞0u+μηvdμ]dη−∫∞−∞[u−ηv|μ=0μ=−∞−∫0−∞u−μηvdμ]dη=∫∞−∞u+η(0,η)v(0,η)dη−∫∞−∞u−η(0,η)v(0,η)dη=∫∞−∞[u+η(0,η)−u−η(0,η)]v(0,η)dη.
Thus, the transmission condition is∫∞−∞[u+η(0,η)−u−η(0,η)]v(0,η)dη=0for v∈C20(Ω).
Solution
Solution
Let a,b∈R and v,w∈D(Ω). It is easy to verify thatFf(av+bw):=⟨Ff,av+bw⟩=∫Ωf(x)[av(x)+bw(x)]dx=a∫Ωf(x)v(x)dx+b∫Ωf(x)w(x)dx=a⟨Ff,v⟩+b⟨Ff,w⟩=aFf(v)+bFf(w),
which shows Ff is a linear mapping on D(Ω). Moreover, let {vj}∞j=1 be a sequence of functions in D(Ω) with supp(vj)⊆K for some compact set K⊊Ω (independent of j) and Dαvj→0 uniformly in K. Then for any given ε>0, there exists a N∈N such that |vj(x)|<ε(∫K|f(x)|dx)−1 for j>N and x∈K. For j>N, we have|Ff(vj)−0|=|⟨Ff,vj⟩|=|∫Ωf(x)vj(x)dx|=|∫Kf(x)vj(x)dx|≤∫K|f(x)||vj(x)|dx≤∫K|f(x)|⋅ε(∫K|f(x)|dx)−1dx=ε.
This implies limj→∞Ff(vj)=0. Therefore, Ff is continuous on D(Ω) and Ff is a distribution in Ω.Solution
Solution
Solution
Solution
For v∈C∞0(R), we find that∫Rfn(x)v(x)dx=n2∫1/n−1/nv(x)dx.
Since v is continuous on R, for any ε>0, there exists δ>0 such that |v(x)−v(0)|<ε whenever |x|<δ. Then for n>1/δ, we have|∫Rfn(x)v(x)dx−v(0)|=n2|∫1/n−1/n[v(x)−v(0)]dx|≤n2∫1/n−1/n|v(x)−v(0)|dx<n2∫1/n−1/nεdx=ε.
This shows limn→∞∫Rfn(x)v(x)dx=v(0)=∫Rδ(x)v(x)dx, i.e., fn(x)→δ(x) as distributions on R.Solution
Let v∈C∞0(R) and K=supp(v). Then we observe that|⟨fn,v⟩−⟨f,v⟩|=|⟨fn−f,v⟩|=|∫Ω(fn(x)−f(x))v(x)dx|=|∫K(fn(x)−f(x))v(x)dx|≤∫K|fn(x)−f(x)||v(x)|dx≤‖v‖C0(K)∫K|fn(x)−f(x)|dx→0as n→∞.
This means limn→∞⟨fn,v⟩=⟨f,v⟩ for v∈D(Ω), i.e., fn→f as distributions on Ω.Solution
Solution
Solution
Solution
Solution
K is continuously differentiable in t if for every v∈C∞0(Rn), f(t)=∫RnK(x,t)v(x)dx defines a continuously differentiable function of t.Solution
Let xi=a+(b−a)i/m for i=0,1,…,m and yj=c+(d−c)j/n for j=1,…,n. Set Ωij=[xi−1,xi]×[yj−1,yj] for 1≤i≤n and 1≤j≤m. Then {Ωij}1≤i≤n,1≤j≤m forms a partition of Ω. Put Δx=(b−a)/m and Δy=(d−c)/n. Then using the summation by parts, we havem∑i=1n∑j=1[u(xi,yj)v(xi,yj)−v(xi,yj−1)Δy+G(u(xi,yj))v(xi,yj)−v(xi−1,yj)Δx]ΔyΔx=m∑i=1[u(xi,yn)v(xi,yn)−u(xi,y1)v(xi,y0)Δy−n∑j=2v(xi,yj−1)u(xi,yj)−u(xi,yj−1)Δy]ΔyΔx+n∑j=1[G(u(xm,yj))v(xm,yj)−G(u(x1,yj))v(x0,yj)Δx−m∑i=2v(xi−1,yj)G(u(xi,yj))−G(u(xi−1,yj))Δx]ΔyΔx=−m∑i=1n∑j=2[v(xi,yj−1)u(xi,yj)−u(xi,yj−1)Δy]ΔyΔx−m∑i=2n∑j=1[v(xi−1,yj)G(u(xi,yj))−G(u(xi−1,yj))Δx]ΔyΔx=−m−1∑i=1n−1∑j=1{v(xi,yj)[u(xi,yj+1)−u(xi,yj)Δy+G(u(xi+1,yj))−G(u(xi,yj))Δx]}ΔyΔx
Here we have used the fact that v∈C∞0(Ω) so v(x0,yj)=v(xm,vj)=v(xi,y0)=v(xi,yn)=0 for 1≤i≤m and 1≤j≤n. Note thatG(u(xi+1,yj))−G(u(xi,yj))=−∫xi+1xiuy(x,yj)dxfor 1≤i≤m−1, 1≤j≤n−1.
Then the double sum becomes−m−1∑i=1n−1∑j=1{v(xi,yj)[u(xi,yj+1)−u(xi,yj)Δy−1Δx∫xi+1xiuy(x,yj)dx]}ΔyΔx→∫ba∫dcv(x,y)[uy(x,y)−uy(x,y)]dxdy=0as n,m→∞.
This shows ∫Ωuvy+G(u)vxdxdy=0 for v∈C10(Ω).Warning: This problem is quite difficult so I cannot promise this solution is correct. Please provide some references or theorem that I may use to make this solution better.
Solution
Suppose by contradiction that m is an odd integer. Then σL(x,ξ)=∑|α|=maα(x)ξα is a polynomial of odd degree. Since aα(x) is not all zero, we may assume ae1(x)=a(m,0,…,0)(x)≠0 by rearranging the indices. Thus, we haveσL(x,ξ)=ame1(x)ξm1+∑|α|=m,α≠me1aα(x)ξα.
Since m is odd and aα(x) are real numbers, by the fundamental theorem of algebra, there exists ξ1=ξ1(ξ2,…,ξm,x) for any (ξ2,…,ξm)≠0. This means the principal symbol has nonzero characteristics, which leads a contradiction. Hence m must be an even integer.
沒有留言:
張貼留言