本文我们考虑一维对流扩散方程
解的逐点衰减估计.此方程可模拟溶质中溶液浓度的衰减情况, 在水利工程, 环境工程及化工, 冶金, 航空等领域都有很重要的应用.方程中$ x \in \mathbb{R} $代表空间变量, $ t $是时间, $ D $是扩散常系数, $ u $是流速, 因为$ \frac{{\partial u}}{{\partial x}} = 0, $所以$ u $只与时间有关, 和空间变量没关系. $ c $可代表溶液浓度.本文中我们研究$ c $随着时间和空间的衰减变化情况.
数学家们对对流扩散方程的研究由来已久. Escobedo和Zuazua在文[1]中研究了方程$ {u_t} - \Delta u = a \cdot ( {{{\left| u \right|}^{q - 1}}u} ) $的柯西问题.由Banach空间不动点理论他们证明了当$ q > 1 $时, 此方程存在唯一经典解$ u \in ( {[ {0,\infty } );{L^1}( {{R^n}} )})$.在文[2]中, Kirane和Qafsaoui考虑了
的解, 此处$ b $是正的有界函数, 且
其中$ {a_{\alpha \beta }} $是正的常系数, $ L_0^ * $是$ {L_0} $的伴随算子.他们得到了解的整体存在性.徐红梅和马慧玲[3]研究了式(1.1) 解的整体存在性.她们通过构造一个Banach空间柯西序列的方法, 得到了此方程解的整体存在性和某些衰减估计.还有很多计算工作者也研究了对流扩散方程的数值解, 这可参看文献[4-6].
本文中, 用$C$表示一般常数, $ {W^{m,p}}( \mathbb{R} ),m \in {Z_ + }, p \in [ {1,\infty } ) $表示Sobolev空间, 模定义为$ {\left\| f \right\|_{{w^{m,p}}}}: = \sum\limits_{k = 0}^m {{{\left\| {\partial _x^kf} \right\|}_{{L_P}}}} $, 特别的$ {H^S}: = {W^{s,2}} $. $ m,n,N $为非负整数. $ F( f ) = \mathop f\limits^ \wedge = \displaystyle\int {{e^{ - ix\xi }}f( x )dx} $表示函数$ f $的傅立叶变换.本文中所有卷积都是关于空间变量$ x $的.本文是在[3]中解整体存在基础上得到的解的逐点衰减估计, 所以先列出文献[3]的结论.
定理1.1 当$ D > u,l \geqslant 1,{\left\| {{c_0}} \right\|_{{H^l}}} = E $充分小, 方程(1.1) 有整体解$ c( {x,t} ) $存在, 且$ {\left\| c \right\|_{{H^l}}} \leqslant CE $.
在定理1.1的基础上, 得到本文结论.
定理1.2 当$ D > u,l \geqslant 1,E = \max ({\left\| {{c_0}} \right\|_{{H^l}}},{\left\| {{c_0}} \right\|_{{L_1}}}) $充分小, $ {c_0} $有紧支集且$ {\left| {{c_0}} \right|_{}} \leqslant C{(1 + {\left| x \right|^2})^{ - 1}} $, 则式(1.1) 的解$ c(x,t) $满足
为更好的理解(1.2) 式, 令$ x = kut, $则
当$ k = 1, $即沿特征线方向$ c(x,t) $以热核的速度衰减, 否则以$ {t^{ - \frac{{1 + n}}{2} - 1}} $的速度衰减.所以可得
注1 方程(1.1) 是典型的双曲-抛物耦合方程, 定理1.2既反映了方程(1.1) 的双曲性质, 即解沿着特征线传播, 也反映了方程的抛物性质, 即解有与热核算子相同的衰减速度.
注2 溶质中溶液浓度在流体流动的方向最大, 这也与物理现象相符.
本文安排如下, 在第二章中, 由duhamel原理给出方程(1.1) 解的表达式, 并给出一些不等式, 解的逐点衰减估计将在第三章给出.
对方程$ \frac{{\partial c}}{{\partial t}} + u\frac{{\partial c}}{{\partial x}} = D{c_{xx}} + {c_{xt}} $两边关于变量$ x $作傅立叶变换, 考虑到$ u $只与时间有关, 和空间变量没关系, 得到
由常微分方程求解得
令
由duhamel原理知方程(1.1) 的解$ c( {x,t} ) $可表示为
为得到$ c(x,t) $的衰减估计, 需要对$ G,H $做详细分析.由(2.2) 式知, 当$ \xi $有界时, $ |\hat G( {\xi ,t} )| \leqslant {e^{ - b{\xi ^2}t}},|\hat H( {\xi ,t} )| \leqslant {e^{ - b{\xi ^2}t}},b \geqslant 0 $.知道$ {F^{ - 1}}({e^{ - b{\xi ^2}t}}) = {(2\pi )^{ - 1}}{t^{ - \frac{1}{2}}}{e^{ - \frac{{{x^2}}}{t}}} $, 当然这里对$ G,H $不可能得到如$ {e^{ - \frac{{{x^2}}}{t}}} $这么好的估计, 但可用$ {B_N}( {x,t} ) = {(1 + \frac{{{x^2}}}{{1 + {t^{}}}})^{ - N}} $来近似代替$ {e^{ - \frac{{{x^2}}}{t}}} $, 对这部分处理结果见引理2.1.当$ \xi $充分大时, 有$ \left| {\hat G} \right| \leqslant {e^{ - bt}},\left| {\hat H} \right| \leqslant {e^{ - bt}},b \geqslant 0 $, 所以这一部分$ G,H $衰减是没有问题的, 但$ \hat G,\hat H $的$ {L_1} $模不存在, 所以在处理高频部分的卷积时, 需要分析$ G,H $的结构, 对这部分处理结果见引理2.2.由以上分析, 需要对频谱分情况讨论, 所以作光滑截断函数$ {\chi _1}(\xi ) = \left\{ \begin{gathered} 1,\quad \left| \xi \right| \leqslant 1 \hfill \\ 0,\quad \left| \xi \right| \geqslant 2 \hfill \\ \end{gathered} \right., $ $ {\chi _2}(\xi ) = \left\{ \begin{gathered} 1\quad \left| \xi \right| \geqslant 2 \hfill \\ 0\quad \left| \xi \right| \leqslant 1 \hfill \\ \end{gathered} \right. $且$ {\chi _1}(\xi ) + {\chi _2}(\xi ) = 1 $.令$ {\hat G_i} = {\chi _i}\hat G,{\hat H_i} = {\chi _i}\hat H,i = 1,2 $.
引理2.1 当$ D > u $时, 对任意的正整数N, 存在仅依赖于N和n的常数$ {C_{N,n}} $, 有
证 当$ \left| \xi \right| $充分小时, 由Taylor展开, 有
所以$ {\hat G_1}( {\xi ,t} ){e^{iu\xi t}} = {e^{( { - ( {D - u} ){\xi ^2} + O( {{\xi ^3}} )} )t}} $, 于是
由文献[7]的引理1得
所以
同理可得
引理2.2 存在正常数$ b $和函数$ g_1^i( x ) $, $ g_2^i( x ) $, $ ( {i = 1,2} ) $, 有
此处$ \delta ( x ) $是Dirac函数, 且
证 当$ |\xi | $充分大, 由泰勒展开得
于是
所以当$ t $充分大时, 存在正常数$ b $有
其中$ m \geqslant 1 $.由文献[7]的引理3.2, 得到结论.
有了解的表达式和对格林函数的分析, 下面给出解的逐点衰减估计.
引理3.1 若$ {c_0} $有紧支集, $ E = \max ({\left\| {{c_0}} \right\|_{{H^l}}},{\left\| {{c_0}} \right\|_{{L_1}}}) $充分小, 当t充分大时, 有
证 由文献[8]的引理2.4及引理2.1有
引理3.2 当$ l \geqslant 1 $, $ E = \max ({\left\| {{c_0}} \right\|_{{H^l}}},{\left\| {{c_0}} \right\|_{{L_1}}}) $, 有
证 由(2.2) 式,
取$ m = 2N $充分大, 则
由引理2.2,
当$ l \geqslant 1, $由Sobolev嵌入不等式$ ||{c_0}|{|_{{L^\infty }}} \leqslant ||{c_0}|{|_{{H^l}}} $和(3.2) 式得到
由(3.1)、(3.3) 式得
引理得证.
在对非线性部分进行估计时, 需先有下面准备工作.
令$ M(t) = \mathop{\sup}\limits_{\begin{subarray}{l} 0 \leqslant s \leqslant t\\ \left| n \right| \leqslant l \end{subarray}} {( {1 + s} )^{\frac{{1 + n}}{2}}}\left| {\partial _x^nc(x,s)} \right|(1 + \frac{{{{\left| {x - us} \right|}^2}}}{{1 + s}}), $则
引理3.3 当$ \left| n \right| \leqslant l $, 有
证 由引理2.1, (3.5) 式, 得
为证明引理, 仅需证明
这可由文[9]的引理3.5得到.
引理3.4 当$ \left| n \right| \leqslant l $, 有
证 由引理2.2和(3.5) 式得
由文[9]的引理3.4, 引理3.6得
由(3.7) 式,
由文[9]的引理3.6, 得
由(3.8), (3.9), (3.10) 式得到引理.
由引理3.1, 引理3.2, 引理3.3, 引理3.4和(2.3) 式, 得到
由$ M(t) $的定义得到$ M( t ) \leqslant CE + C{M^2}( t ) $.因为$ E $充分小, 所以$ M( t ) $有界.由(3.4) 式, 得到本文结论, 即定理1.2.