Burgers方程初值问题解的有效区域
{ u t + u u x = 0 , t > 0 ; u ∣ t = 0 = f ( x ) \begin{cases} u_t + u u_x = 0, & t > 0; \\ u|_{t=0} = f(x) \end{cases} {ut+uux=0,u∣t=0=f(x)t>0;
并描述在 ( x , t ) (x,t) (x,t) 平面上解有效定义的区域,使用以下初始数据之一:
f ( x ) = tanh ( x ) ; f ( x ) = − tanh ( x ) ; f(x) = \tanh(x); \quad f(x) = -\tanh(x); f(x)=tanh(x);f(x)=−tanh(x);
f ( x ) = { − 1 x < − a , x a − a ≤ x ≤ a , 1 x > a ; f(x) = \begin{cases} -1 & x < -a, \\ \frac{x}{a} & -a \leq x \leq a, \\ 1 & x > a; \end{cases} f(x)=⎩⎪⎨⎪⎧−1ax1x<−a,−a≤x≤a,x>a;
f ( x ) = { 1 a x < − a , − x a − a ≤ x ≤ a , − 1 x > a ; f(x) = \begin{cases} \frac{1}{a} & x < -a, \\ -\frac{x}{a} & -a \leq x \leq a, \\ -1 & x > a; \end{cases} f(x)=⎩⎪⎨⎪⎧a1−ax−1x<−a,−a≤x≤a,x>a;
f ( x ) = { − 1 x < 0 , 1 x > 0 ; f(x) = \begin{cases} -1 & x < 0, \\ 1 & x > 0; \end{cases} f(x)={−11x<0,x>0;
f ( x ) = { sin ( x ) ∣ x ∣ < π , 0 ∣ x ∣ > π , f(x) = \begin{cases} \sin(x) & |x| < \pi, \\ 0 & |x| > \pi, \end{cases} f(x)={sin(x)0∣x∣<π,∣x∣>π,
其中 a > 0 a > 0 a>0 是参数。
问题求解
该问题涉及一阶拟线性偏微分方程(无粘性Burgers方程):
u t + u u x = 0 , t > 0 , u_t + u u_x = 0, \quad t > 0, ut+uux=0,t>0,
初始条件为 u ( x , 0 ) = f ( x ) u(x,0) = f(x) u(x,0)=f(x)。方程的特征线方程为 d x d t = u \frac{dx}{dt} = u dtdx=u,沿特征线 u u u 为常数。解为 u ( x , t ) = f ( x 0 ) u(x,t) = f(x_0) u(x,t)=f(x0),其中 x 0 x_0 x0 满足 x = x 0 + f ( x 0 ) t x = x_0 + f(x_0) t x=x0+f(x0)t。解在特征线相交前有效(即激波形成前),此时解为经典解(光滑或连续)。激波形成时间由 1 + f ′ ( x 0 ) t = 0 1 + f'(x_0) t = 0 1+f′(x0)t=0 的最小正解决定(当 f ′ ( x 0 ) < 0 f'(x_0) < 0 f′(x0)<0 时)。
针对每个初始数据 f ( x ) f(x) f(x),求解并描述解有效定义的区域(即经典解存在的区域)。区域以 t t t 的范围和 x ∈ R x \in \mathbb{R} x∈R 给出,除非特别说明。
1. f ( x ) = tanh ( x ) f(x) = \tanh(x) f(x)=tanh(x)
- f ′ ( x ) = sech 2 ( x ) > 0 f'(x) = \text{sech}^2(x) > 0 f′(x)=sech2(x)>0 对所有 x x x。
- 特征线无相交,解全局光滑。
- 解有效定义的区域: t ≥ 0 t \geq 0 t≥0,所有 x x x.
2. f ( x ) = − tanh ( x ) f(x) = -\tanh(x) f(x)=−tanh(x)
- f ′ ( x ) = − sech 2 ( x ) < 0 f'(x) = -\text{sech}^2(x) < 0 f′(x)=−sech2(x)<0 对所有 x x x。
- 激波形成时间 t = 1 t = 1 t=1(在 x 0 = 0 x_0 = 0 x0=0 处 f ′ ( x 0 ) = − 1 f'(x_0) = -1 f′(x0)=−1, t = − 1 / f ′ ( x 0 ) = 1 t = -1/f'(x_0) = 1 t=−1/f′(x0)=1)。
- 对于 t < 1 t < 1 t<1,解光滑; t = 1 t = 1 t=1 时激波形成。
- 解有效定义的区域: 0 ≤ t < 1 0 \leq t < 1 0≤t<1,所有 x x x.
3. f ( x ) = { − 1 x < − a x a − a ≤ x ≤ a 1 x > a f(x) = \begin{cases} -1 & x < -a \\ \frac{x}{a} & -a \leq x \leq a \\ 1 & x > a \end{cases} f(x)=⎩⎪⎨⎪⎧−1ax1x<−a−a≤x≤ax>a
- f ( x ) f(x) f(x) 连续,且 f ′ ( x ) ≥ 0 f'(x) \geq 0 f′(x)≥0( ∣ x ∣ > a |x| > a ∣x∣>a 时 f ′ = 0 f' = 0 f′=0, ∣ x ∣ < a |x| < a ∣x∣<a 时 f ′ = 1 / a > 0 f' = 1/a > 0 f′=1/a>0)。
- 特征线无相交(在 ∣ x ∣ < a |x| < a ∣x∣<a 区域特征线发散)。
- 解有效定义的区域: t ≥ 0 t \geq 0 t≥0,所有 x x x.
4. f ( x ) = { 1 a x < − a − x a − a ≤ x ≤ a − 1 x > a f(x) = \begin{cases} \frac{1}{a} & x < -a \\ -\frac{x}{a} & -a \leq x \leq a \\ -1 & x > a \end{cases} f(x)=⎩⎪⎨⎪⎧a1−ax−1x<−a−a≤x≤ax>a
- 在 x = − a x = -a x=−a 处不连续(除非 a = 1 a = 1 a=1)。
- 情况分析:
- 若 a > 1 a > 1 a>1:
- 在 x = − a x = -a x=−a 处,左状态 u l = 1 / a < 1 = u r u_l = 1/a < 1 = u_r ul=1/a<1=ur(右状态),形成稀疏波。
- 在 ∣ x ∣ < a |x| < a ∣x∣<a 区域, f ′ ( x ) = − 1 / a < 0 f'(x) = -1/a < 0 f′(x)=−1/a<0,激波在 t = a t = a t=a 形成。
- 解在 0 < t < a 0 < t < a 0<t<a 连续且光滑(稀疏波和线性区域连接)。
- 解有效定义的区域: 0 < t < a 0 < t < a 0<t<a,所有 x x x。
- 若 a < 1 a < 1 a<1:
- 在 x = − a x = -a x=−a 处, u l = 1 / a > 1 = u r u_l = 1/a > 1 = u_r ul=1/a>1=ur,形成激波(立即在 t = 0 + t = 0^+ t=0+ 形成)。
- 无经典解区域。
- 解有效定义的区域:无(经典解不存在)。
- 若 a = 1 a = 1 a=1:
- f ( x ) = { 1 x < − 1 − x − 1 ≤ x ≤ 1 − 1 x > 1 f(x) = \begin{cases} 1 & x < -1 \\ -x & -1 \leq x \leq 1 \\ -1 & x > 1 \end{cases} f(x)=⎩⎪⎨⎪⎧1−x−1x<−1−1≤x≤1x>1,连续。
- 在 ∣ x ∣ < 1 |x| < 1 ∣x∣<1 区域, f ′ ( x ) = − 1 < 0 f'(x) = -1 < 0 f′(x)=−1<0,激波在 t = 1 t = 1 t=1 形成。
- 解有效定义的区域: 0 ≤ t < 1 0 \leq t < 1 0≤t<1,所有 x x x.
- 若 a > 1 a > 1 a>1:
5. f ( x ) = { − 1 x < 0 1 x > 0 f(x) = \begin{cases} -1 & x < 0 \\ 1 & x > 0 \end{cases} f(x)={−11x<0x>0
- 在 x = 0 x = 0 x=0 处不连续, u l = − 1 < 1 = u r u_l = -1 < 1 = u_r ul=−1<1=ur。
- 形成稀疏波(中心在 x = 0 x = 0 x=0),解在 t > 0 t > 0 t>0 连续。
- 无激波形成(特征线发散)。
- 解有效定义的区域: t > 0 t > 0 t>0,所有 x x x(在 t = 0 t = 0 t=0 不经典)。
6. f ( x ) = { sin ( x ) ∣ x ∣ < π 0 ∣ x ∣ > π f(x) = \begin{cases} \sin(x) & |x| < \pi \\ 0 & |x| > \pi \end{cases} f(x)={sin(x)0∣x∣<π∣x∣>π
- f ( x ) f(x) f(x) 连续(在 x = ± π x = \pm \pi x=±π 处 f = 0 f = 0 f=0)。
- f ′ ( x ) = cos ( x ) f'(x) = \cos(x) f′(x)=cos(x) 在 ∣ x ∣ < π |x| < \pi ∣x∣<π,在 π / 2 < ∣ x ∣ < π \pi/2 < |x| < \pi π/2<∣x∣<π 时 f ′ < 0 f' < 0 f′<0。
- 激波形成时间 t = 1 t = 1 t=1(在 x 0 → ± π ∓ x_0 \to \pm \pi^{\mp} x0→±π∓ 处 f ′ ( x 0 ) → − 1 f'(x_0) \to -1 f′(x0)→−1, t = − 1 / f ′ ( x 0 ) → 1 t = -1/f'(x_0) \to 1 t=−1/f′(x0)→1)。
- 对于 t < 1 t < 1 t<1,解光滑。
- 解有效定义的区域: 0 ≤ t < 1 0 \leq t < 1 0≤t<1,所有 x x x.
总结
解有效定义的区域指经典解(连续且光滑)存在的区域。结果总结如下:
初始数据 f ( x ) f(x) f(x) | 解有效定义的区域(在 ( x , t ) (x,t) (x,t) 平面上) |
---|---|
f ( x ) = tanh ( x ) f(x) = \tanh(x) f(x)=tanh(x) | t ≥ 0 t \geq 0 t≥0,所有 x x x |
f ( x ) = − tanh ( x ) f(x) = -\tanh(x) f(x)=−tanh(x) | 0 ≤ t < 1 0 \leq t < 1 0≤t<1,所有 x x x |
f ( x ) = { − 1 x < − a x a − a ≤ x ≤ a 1 x > a f(x) = \begin{cases} -1 & x < -a \\ \frac{x}{a} & -a \leq x \leq a \\ 1 & x > a \end{cases} f(x)=⎩⎪⎨⎪⎧−1ax1x<−a−a≤x≤ax>a | t ≥ 0 t \geq 0 t≥0,所有 x x x |
f ( x ) = { 1 a x < − a − x a − a ≤ x ≤ a − 1 x > a f(x) = \begin{cases} \frac{1}{a} & x < -a \\ -\frac{x}{a} & -a \leq x \leq a \\ -1 & x > a \end{cases} f(x)=⎩⎪⎨⎪⎧a1−ax−1x<−a−a≤x≤ax>a | 若 a > 1 a > 1 a>1: 0 < t < a 0 < t < a 0<t<a,所有 x x x 若 a < 1 a < 1 a<1: 无 若 a = 1 a = 1 a=1: 0 ≤ t < 1 0 \leq t < 1 0≤t<1,所有 x x x |
f ( x ) = { − 1 x < 0 1 x > 0 f(x) = \begin{cases} -1 & x < 0 \\ 1 & x > 0 \end{cases} f(x)={−11x<0x>0 | t > 0 t > 0 t>0,所有 x x x |
$$ f(x) = \begin{cases} \sin(x) & | x |
注:
- 对于不连续的初始数据(如情况4、5),解在 t > 0 t > 0 t>0 可能连续(通过稀疏波),但仅在指定区域内经典。
- 激波形成后,解需引入激波条件,但本问题仅要求经典解区域。
问题6的解答: 初始数据 f ( x ) 解有效定义的区域 tanh ( x ) t ≥ 0 , x ∈ R − tanh ( x ) 0 ≤ t < 1 , x ∈ R { − 1 x < − a x a − a ≤ x ≤ a 1 x > a t ≥ 0 , x ∈ R { 1 a x < − a − x a − a ≤ x ≤ a − 1 x > a { a > 1 : 0 < t < a , x ∈ R a < 1 : 无 a = 1 : 0 ≤ t < 1 , x ∈ R { − 1 x < 0 1 x > 0 t > 0 , x ∈ R { sin ( x ) ∣ x ∣ < π 0 ∣ x ∣ > π 0 ≤ t < 1 , x ∈ R \boxed{ \begin{array}{c} \text{问题6的解答:} \\ \\ \begin{array}{c|c} \text{初始数据 } f(x) & \text{解有效定义的区域} \\ \hline \tanh(x) & t \geq 0, \, x \in \mathbb{R} \\ -\tanh(x) & 0 \leq t < 1, \, x \in \mathbb{R} \\ \begin{cases} -1 & x < -a \\ \dfrac{x}{a} & -a \leq x \leq a \\ 1 & x > a \end{cases} & t \geq 0, \, x \in \mathbb{R} \\ \begin{cases} \dfrac{1}{a} & x < -a \\ -\dfrac{x}{a} & -a \leq x \leq a \\ -1 & x > a \end{cases} & \begin{cases} a > 1: & 0 < t < a, \, x \in \mathbb{R} \\ a < 1: & \text{无} \\ a = 1: & 0 \leq t < 1, \, x \in \mathbb{R} \end{cases} \\ \begin{cases} -1 & x < 0 \\ 1 & x > 0 \end{cases} & t > 0, \, x \in \mathbb{R} \\ \begin{cases} \sin(x) & |x| < \pi \\ 0 & |x| > \pi \end{cases} & 0 \leq t < 1, \, x \in \mathbb{R} \\ \end{array} \end{array} } 问题6的解答:初始数据 f(x)tanh(x)−tanh(x)⎩⎪⎪⎨⎪⎪⎧−1ax1x<−a−a≤x≤ax>a⎩⎪⎪⎨⎪⎪⎧a1−ax−1x<−a−a≤x≤ax>a{−11x<0x>0{sin(x)0∣x∣<π∣x∣>π解有效定义的区域t≥0,x∈R0≤t<1,x∈Rt≥0,x∈R⎩⎪⎨⎪⎧a>1:a<1:a=1:0<t<a,x∈R无0≤t<1,x∈Rt>0,x∈R0≤t<1,x∈R