一类带有非线性阻尼项的磁流体动力学方程组的解的整体存在性.pdf

应用数学MATHEMATICA APPLICATA2024,37(1):63-72一类带有非线性阻尼项的磁流体动力学方程组的解的整体存在性李林锐1,洪明理1,郑琳2(1.防灾科技学院基础部,河北 三河 0652012.华北水利水电大学数学与统计学院,河南 郑州 450046)摘要：本文研究在多孔介质意义下的一类带有非线性阻尼项a|u|1u(a 0)的不可压的磁流体动力学方程组的解的整体存在性问题,通过古典的能量方法和Sobolev紧性嵌入方法获得了解的整体存在性,利用Gagliardo-Nirenberg插值不等式和其它的一些重要不等式得到了解的正则性结果,建立了弱解和强解的整体存在性,这些结果在很大程度改善了之前相关文献的结果,揭示了磁流体运动的物理现象,为磁流体动力学的发展提供了必要的理论基础.关键词：磁流体动力学方程组;阻尼项;粘性流;Sobolev紧性嵌入中图分类号：O175.29AMS(2010)主题分类：35Q35;76B03文献标识码：A文章编号：1001-9847(2024)01-0063-101.引言本文将考虑如下带有非线性阻尼项的不可压缩磁流体动力学方程组的初始值问题ut+(u )u (b )b u+a|u|1u+p=0,bt+(u )b (b )u b=0,u=b=0,u(x,0)=u0(x),b(x,0)=b0(x),(1.1)其中u表示流体的速度场,b表示磁场,流体的速度场u和磁场b是二维的自由散度场,p表示压力,0表示流体的粘性系数,0表示磁场的扩散系数,阻尼项系数a 0和指标函数 1,阻尼项一般来自于流体运动对外力的抵抗,可以描述诸如多孔介质的流体的运动,摩擦力和拖拽力以及一些耗散项对流体运动的影响等14.当a=0 时,系统(1.1)退化成常规意义下的磁流体动力学方程组,该系统是流体力学的基本宏观模型之一,已经得到广泛而深入的研究26.系统(1.1)中第一个方程表示动量守恒方程；第二个方程是电磁感应方程,这个方程组描述了电场和磁场的耦合作用.磁流体动力学偏微分方程组(magnetohydrodynamics,简称MHD方程组)描述了导电流体在电磁场中运动的状态,在空气动力学、天体物理、地球物理以及宇宙等离子物理学等领域具有非常重要的应用,与其相关的渐近极限问题因其物理背景的重要性、复杂性,其数学方面的挑战性,吸引了许多知名数学家和物理学家的研究兴趣,同时也取得了很多很好的结果.基本的数学问题比如解的适定性问题、整体正则性和解的爆破性有很多深收稿日期：2022-11-19基金项目：中央高校基本科研业务费专题研究项目(ZY20215116);国家自然科学基金项目(12071192)作者简介：李林锐,女,汉族,河南人,副教授,研究方向:流体力中的偏微分方程.64应用数学2024入的研究和结果,在多孔介质中带有非线性阻尼项的磁流体动力学模型更多的物理背景可参见文7-9,14-15.显然,当方程组中的磁场b=0时,系统(1.1)退化成带有非线性阻尼项的不可压的Navier-Stokes方程.近些年,带有非线性阻尼项的Navier-Stokes方程吸引了很多国内外专家学者的关注510,18.蔡晓静和酒全森5研究了带有阻尼项a|u|1u(a 0)的Navier-Stokes方程弱解和强解的整体存在性,证明了当 72时强解的整体存在性和72 5时强解的唯一性;随后,张族锦等6将的下界降低到3,这个下界3具有重要的关键性作用也在文7中被周勇得到证实.受到这些文献的启发,我们将关于带有阻尼项的Navier-Stokes的解的适定性问题推广到带有阻尼项的磁流体动力学方程组.在先前相关的磁流体动力学方程组的研究工作中,更多的专家学者聚焦于不带有非线性阻尼项或者只带有部分粘性项或者部分磁扩散项的广义的磁流体动力学方程组解的存在性和唯一性证明,其中陈文吉、章志飞和周建丰11研究了在周期性区域上带有部分粘性扩散项的三维磁流体动力学方程组在初始速度充分小和初始磁场接近于背景磁场并且满足Diophantine条件时的解整体适定性;潘荣华、周忆和朱忆12在欧拉坐标系下利用带有时间权函数的能量估计方法研究了在初始磁场接近于平衡状态下并且初始条件具有某种对称性条件下三维磁流体动力学方程组不带有磁扩散项的磁流体动力学方程组的古典解的整体存在性;吴家宏和翟晓平13研究了带有真空的非电阻的可压的三维磁流体动力学方程组的柯西问题下局部强解的存在性,并首次给出了三维磁流体动力学方程组在磁场满足Diophantine条件下磁场接近于背景磁场时的解的整体存在性和稳定性;Titi和Trabelsi15解决了在多媒介质通道中三维磁流体动力学方程组的解的适定性问题.综合以上分析,在本文中我们将研究在多孔介质意义下的一类带有非线性阻尼项的不可压的磁流体动力学方程组的解的整体存在性问题,由于磁流体动力学方程组有磁场的耦合作用,在用能量方法的过程中计算更复杂也更富有挑战性,在计算过程中会出现新的问题、难度更大,面对的困难也将更多.同时,在文15中作者证明了同时包含速度场方程和磁场方程中的带有两个非线性阻尼项下的磁流体动力学方程组的强解的存在性,在本文中我们并不需要额外的磁场方向的阻尼项就获得了解的存在性,此外,也给出了弱解的存在性,这同时也推广了文15中的相关结果.2.整体存在性结果对于带有阻尼项的磁流体动力学方程组,首先给出方程组弱解的定义:定义1如果(u,b,p)满足下列三个条件：1)u L(0,T;L2(R2)L2(0,T;H1(R2)L+1(0,T;L+1(R2),b L(0,T;L2(R2)L2(0,T;H1(R2);2)Cc(R2 0,T),div=0,T0R2u tdxdt+T0R2u dxdt+T0R2(u )u (b )b+a|u|1u dxdt=R2u0(x,0)dx,(2.1)且T0R2b tdxdt+T0R2b dxdt+T0R2(u )b (b )udxdt=R2b0(x,0)dx.(2.2)3)(x,t)R2(0,T),都有 u(x,t)=b(x,t)=0成立.本文主要结果如下:第 1 期李林锐等：一类带有非线性阻尼项的磁流体动力学方程组的解的整体存在性65定理1假设初始条件(u0,b0)L2(R2)L2(R2)且满足divu0=divb0=0,令a,0且 1,则带有非线性阻尼项的磁流体动力学方程组(I)有整体弱解(u(x,t),b(x,t)满足u(x,t)L(0,+);L2(R2)L2Loc(0,+);H1(R2)L+1Loc(0,+);L+1)和b(x,t)L(0,+);L2(R2)L2Loc(0,+);H1(R2).特别地,有limt+(u(t)22+b(t)22)=0.此外,limt+min,+t(u()22+b()22+au()+1+1)d=0.进一步地,也可以得到带有非线性阻尼项的磁流体动力学方程组的强解,结论如下:定理2假设初始条件(u0,b0)H1(R2)H1(R2)且满足divu0=divb0=0,令a,0且 52,则带有非线性阻尼项的磁流体动力学方程组(I)有整体强解(u(x,t),b(x,t)满足u(x,t)C0b(0,+);H1(R2)L2Loc(0,+);H2(R2)L+1Loc(0,+);L+1(R2)和b(x,t)C0b(0,+);H1(R2)L2Loc(0,+);H2(R2).此外,如果初始条件u0 L+1(R2),则有tu(x,t)L2Loc(0,+);L2(R2),u LLoc(0,+);L+1(R2)和tb(x,t)L2Loc(0,+);L2(R2).3.预备知识首先我们给出一些在证明过程中用到的引理:引理18若f,g Hk L,则fgHk CfLgHk+fHkgL,当|k时,有D(fg)fDgL2 C(fHkgL+fLgHk1).如果F是光滑函数并且F(0)=0,则对于任意的f Hk L,有F(f)Hk CkfL(1+fHk).引理216(Agmons不等式)若u H2()H10(),R2,则存在常数C使得uL()Cu12H1()u12H2(),uL()Cu12L2()u12H2().引理317(Gagliardo-Nirenberg不等式)若f(x)是定义在R2上的光滑函数,如果1 q,r,m是自然数,假设实数和自然数j满足1p=j2+(1rm2)+1 q,和jm 1,则有DjfLp C1DmfLrf1Lq+C2fLs,66应用数学2024其中s 0是任意数,并且常数C1和C2仅仅依赖于,m,j和s.4.定理的证明为了证明带有非线性阻尼项的磁流体动力学方程组弱解的整体存在性,首先用Faedo-Galerkin近似方法去获得近似解的局部存在性,这一部分方法很经典,可以参考文8中关于Navier-Stokes方程近似解的相关步骤,在此我们不在赘述.我们仅需要获得关于解的一致先验估计去延拓局部解到整体解,并最终应用Aubin紧性原理由近似解的先验估计获得原始系统整体解的存在性.下面将通过以下几个步骤分别获得解的先验估计从而给出定理的证明.定理1的证明步1(u,b)的L2估计对(1.1)的第一个方程两边同时乘以u(x,t),然后在R2上积分可得12ddtu22+u22+au+1+1=R2b b udx,(4.1)对(1.1)的第二个方程两边同时乘以b(x,t),然后在R2上积分可得12ddtb22+b22=R2b u bdx=R2b b udx.(4.2)将上式(4.1)和(4.2)式相加可得,ddt(u22+b22)+2min,(u22+b22)+2au+1+1 0,(4.3)将上式(4.3)式从0到t上积分可得u22+b22 u022+b022,(4.4)t0(u()22+b()22)d u022+b0222min,(4.5)t0u()+1+1d u022+b0222a.(4.6)在上面最后的两个式子中分别令(4.5)、(4.6)式子中t +,则有+0(u()22+b()22)d u022+b0222min,(4.7)+0u()+1+1d u022+b0222a.(4.8)由上面的能量估计可得u(x,t)L(0,+);L2(R2)L+1(0,+);L+1(R2),(4.9)u(x,t)L2(0,+);L2(R2).(4.10)和b(x,t)L(0,+);L2(R2),b(x,t)L2(0,+);L2(R2).(4.11)进一步地,应用H older不等式和Young不等式,可得u22 au+1+1+M3a1()+11.(4.12)利用Poincar e不等式,(4.2)式可变形为ddt(u(t)22+b(t)22+2min,(u(t22+b(t)22)2M3a1()+11.(4.13)通过Gronwalls引理可得,对于任意的t 0,u(t)22+b(t)22(u022+b022)e2min,t+M3()+11a1min,(1 e2min,t).(4.14)第 1 期李林锐等：一类带有非线性阻尼项的磁流体动力学方程组的解的整体存在性67下面令 (0,1),则有R2|u|2dx=|u|u|2dx+|u|u|2dx 2M3+2|u|(|u|)2dx 2M3+2|u|(|u|)+1dx 2M3+52|u|u|+1dx 2M3+52R2|u|+1dx.(4.15)从而可得2a52u22 292M3+2au+12+12.(4.16)将(4.16)代入(4.1)式并加上(4.2)式可得ddt(u(t)22+b(t)22)+2u22+2b(t)22+2a52u22 2a52M3.(4.17)利用Poincar e不等式,在(4.17)式中扔掉2u22,则有ddt(u(t)22+b(t)22)+2b(t)22+2a52u22 2a52M3.(4.18)由于 (0,1),从而可得ddt(u(t)22+b(t)22)+22min,a(u(t)22+b(t)22)2a52M3.(4.19)由Gronwall不等式,可得u(t)22+b(t)22 e22min,at(u022+b022)+a52min,a(1 e22min,at).(4.20)从而可得limt+sup(u(t)22+b(t)22)aM3min,a52.(4.21)在上面式子中令 0,则有 limt+(u(t)22+b(t)22)=0.对式子(4.2)两边同时从(m,t)上积分可得2min,tm(u()22+b()22)d+2atmu()+1+1d u(s)22+b(s)22.(4.22)在(4.22)中不等式左边是有界的并且是关于时间t是单调的,首先令t ,然后令m ,则可以得到:limmmin,+m(u()22+b()22)d+amu()+1+1=0.(4.23)步2(u,b)的H1估计在(1.1)的第一个方程两边同时乘以u,并在R2上分部积分可得12ddtu22+u22 aR2|u|1u udx=R2(u )u udx R2(b )b udx.(4.24)由于aR2|u|1u udx=a|u|12u22.(4.25)68应用数学2024另一方面,由于 3,利用H older不等式和Young不等式可得|R2(u )u udx|R2|u|u|21|u|121|u|dx u|u|211|u|1212(1)3u22|u|12u412u2(121)2+8u222|u|12u22+223u22+8u22.(4.26)对于(4.24)式右边最后一项利用Cauchy-Schwarz不等式和Young不等式可得|R2(b )budx|2b|b|22+8u222b2+13(1)b26(+1)+5+8u222b+13(1)b4+12b2(1)+12+8u22C+1212b22+6b22+8u22.(4.27)在上面式子中最后两步的估计式中用到引理3中的Gagliardo-Nirenberg不等式和Young不等式.接下来,在(1.1)中第二个方程两边同时乘以b并在R2上积分可得12ddtb22+b22=R2(u )b bdx R2b u bdxR2(u )b bdx+32b|u|22+6b22R2(u)b bdx+32b23+11u26+1+5R2(u )b bdx+Cb23+11u4+12u2(1)(+1)2R2(u )b bdx+C+12+12b+13+11u22+8u22+6b22.(4.28)下面估计(4.28)式子右边的第一项,通过分部积分可得R2(u )b bdx=3l,m,n=1R2ullbm2nbmdx=3l,m,n=1R2nullbmnbmdx 3l,m,n=1R2ulnlbmnbmdx=3l,m,n=1R2nullbmnbmdx 3l,m,n=1R2u nbnbdx=3l,m,n=1R22nullbmbmdx+3l,m,n=1R2nulnlbmbmdx3n=1R2u nb nbdx第 1 期李林锐等：一类带有非线性阻尼项的磁流体动力学方程组的解的整体存在性69=R2(u )b bdx+3l,m,n=1R2nulnlbmbmdx 3n=1R2u nbnb=R2(nu )nbbdx 3n=1R2u nb nbdx.(4.29)综合上式可得|R2(u )b bdx|R2|u|b|b|dx|.(4.30)下面类似于在(4.28)式右边第二项的估计可得|R2(u )b bdx|C12+12b+13+11u22+8u22+6b22.(4.31)联合(4.24)-(4.31)可得ddt(u22+b22)+u22+b22+a|u|12u22 C(a12)23+112+12b+13+11u22C+1212b22.(4.32)由于H1(R2),L3+11,如果b0 L3+11,则由后面的估计式(4.53)可知b(t)LLoc(0,+);L3+11(R2),应用Gronwall不等式,通过上面式子(4.31)可得u,b LLoc(0,+);L2(R2),u,b L2Loc(0,+),L2(R2),(4.33)|u|12u L2Loc(0,+),L2(R2).(4.34)综合上面的能量估计式(4.9)、(4.10)、(4.14)、(4.21)、(4.23)和(4.34)式可知,定理1得到证明.定理2的证明首先对(1.1)的第一个方程两边同时乘以tu并在区域R2上积分可得tu22+2ddtu22+R2(u )u tudx+R2(b )b tudx+2+1ddtu+1+1 0.(4.35)对于(4.35)式左端第二项和第三项分别有估计式|R2(u )u tudx|u|u|22+14tu22,(4.36)|R2(b )b tudx|b|b|22+14tu22.(4.37)将(4.36)、(4.37)式分别代入(4.35)式可得tu22+ddtu22+2+1ddtu+1+1 2u|u|22+2b|b|22.(4.38)两边同时关于时间t积分可得t0tu22d+u22+2+1u+1+1u022+2+1u0+1+1+2t0u|u|22d+2t0b|b|22d.(4.39)利用Young不等式可以估计u|u|22R2|u|2|u|41|u|241dx 112|u|12|u|22+u22.(4.40)70应用数学2024同理可估计b|b|22 b22b2 Cb22b322b122 Cb82b22+b22.(4.41)在上面最后一步应用到引理2中的Agmon不等式.由于上面的估计式(4.32)可得u|u|,b|b|L2Loc(0,+),L2(R2).(4.42)将上式(4.42)代入(4.38)式可得tu L2Loc(0,+),L2(R2).(4.43)下面用同样的方法证明tb L2Loc(0,+),L2(R2).对(1.1)式的第二个方程两边同时乘以tb并在区域R2上积分可得tb22+2ddtb22=R2(u )b tbdx+R2(b )u tbdx u|b|22+14tb22+b|u|22+14tb22.(4.44)对(4.44)式两边从0到t上积分可得t0tb22d+b22 b022+2t0(u()22b()22+b()22u()22)d.(4.45)再次应用Gagliardo-Nirenberg不等式和Young不等式可得u|b|22 u2b22 Cu22u22b22+Cu22b22 u22+Cu42b42+Cu22b22.(4.46)同理可得b|u|22 b23+11u26+1+51b23+11u4+12u2(1)+121b+13+11u22+16u22.(4.47)其中在上面估计式子中需要用到b的L3+11(R2)的估计,下面我们将补充该估计式.在(1.1)的第二个方程两边同时乘以|b|+11b并在R2上同时积分可得:13 1+1ddtb3+113+11 R2b b|b|+51dx=R2(b )u b|b|+51dx=R2b (|b|+51b)udx.(4.48)对于上式中左边第二项,利用Stroock-Varopoulos不等式可得R2b b|b|+51dx 492(1)(+1)(+1)2|b|32+112259|b|32+1122.(4.49)另一方面,对于上式中右边的项,估计式如下:R2b (|b|+51b)udx 13+5+1R2|b|32+11|b|32+11|u|dx23R2|b|32+11|b|32+11|u|dx29|b|32+1122+2u|b|32+112229|b|32+1122+2u2+1|b|32+1122+11第 1 期李林锐等：一类带有非线性阻尼项的磁流体动力学方程组的解的整体存在性7129|b|32+1122+2u2+1|b|32+116+12b3213+1129|b|32+1122+29|b|32+1122+C62u2+22+1b3+113+11.(4.50)又由于u2+22+1 1+u+1+1.(4.51)联合(4.49)、(4.50)和(4.51)式,并代入(4.48)式可得ddtb3+113+11+3|b|32+1122 C62(1+u+1+1)b3+113+11.(4.52)由于先验估计式(4.6)可知u L+1(0,+);L+1(R2),从而有b(t)3+11 b03+11expC62t0(1+u()+1+1)d b03+11expC62(t+u022+b0222a).(4.53)从而可得b(t)LLoc(0,+);L3+11(R2).最终,联合(4.53)和式子(4.32)式,从而可得u|b|,b|u|L2Loc(0,+),L2(R2),代入(4.45)式最终可得tb L2Loc(0,+);L2(R2).根据tu和tb的有界性,同时综合上面的式子(4.34)、(4.39)和(4.43),从而完成了定理2的证明.注由于系统(1.1)为同时具有粘性项、磁耗散项和非线性阻尼项的磁流体动力学方程组,因此现有文献如文1821中定理都不适用于系统(1.1),从而说明本文的定理推广、改进且丰富了现有文献中关于磁流体动力学方程组的结果,尤其是完善了之前相关文献中的一些数值分析结果,为磁流体动力学的发展和持续研究提供了理论依据.参考文献：1 HUANG F M,PAN R H.Convergence rate for compressible Euler equations with damping andvacuumJ.Arch.Ration.Mech.Anal.,2003,166:359-376.2 HSIAO L.Quasilinear Hyperbolic Systems and Dissipative MechanismsM.Singapore:World Scien-tific,1998.3 BRESCH D,DESJARDINS B,LIN C K.On some compressible fluid models:Korteweg,lubrication,and shallow water systemsJ.Comm.Partial Differential Equations.,2003,28(3-4):843-868.4 BRESCH D,DESJARDINS B.Existence of global weak solutions for a 2D viscous shallow waterequations and convergence to the quasi-geostrophic modelJ.Comm.Math.Phys.,2003,238(1-2):211-223.5 CAI X J,JIU Q S,Weak and strong solutions for the incompressible Navier-Stokes equations withdampingJ.Journal of Mathematical Analysis and Applications.,2008,343:799-809.6 ZHANG Z J,WU X L,LU M.On the uniqueness of strong solution to the incompressible Navier-Stokes equations with dampingJ.J.Math.Anal.Appl.,2011,377:414-419.7 ZHOU Y.Regularity and uniqueness for the 3D incompressible Navier-Stokes equations with damp-ingJ.Appl.Math.Lett.,2012,25:1822-1825.8 ANTONTSEV S N,OLIVEIRA H B.The Navier-Stokes problem modified by an absorption termJ.Applicable Analysis,2010,89,12:1805-1825.9 ZHOU Y.Remarks on regularities for the 3D MHD equationsJ.Discrete Contin.Dyn.Syst.,2005,12:881-886.10 LIU Yanlin,ZHANG Ping.Global solutions of 3-D Navier-Stokes system with small unidirectionalderivativeJ.Arch.Ration.Mech.Anal.,2020,235:1405-1444.72应用数学202411 CHEN W,ZHANG Z,ZHOU J.Global well-posedness for the 3-D MHD equations with partialdiffusion in periodic domainJ.Sci China Math,2022,65:309-318.12 PAN R H,ZHOU Y,ZHU Y.Global classical solutions of three dimensional viscous MHD systemwithout magnetic diffusion on periodic boxesJ.Arch.Ration.Mech.Anal.,2018,227:637-662.13 WU J,ZHAI X P.Global small solutions to the 3D compressible viscous Non-resistive MHDsystemEB/OL.(2022-11-11)2022-11-19.https:/arxiv.org/abs/2211.06231.14 ZHANG Z J,WU C P,YAO Z A.Remarks on global regularity for the 3D MHD system withdampingJ.Applied Mathematics and Computation.,2018,333:1-7.15 TITI E S,TRABELSI S.Global well-posedness of a 3D MHD model in porous mediaJ.Journal ofGeometric mechanics,2019,11(4):621-637.16 AGMON S.Lectures on Elliptic Boundary Value ProblemsM.Providence,RI:AMS Chelsea pub-lishing,2010.17 LIONS P L.Mathematical topics in fluid mechanics.Vol 2.Compressible models-Oxford SciencePublicationsM.Oxford Lecture Series in Mathematics and Its Applications,10.New York:TheClarendon Press,Oxford University Press,1998.18 ZHANG Z J,ZHONG D X,GAO S J,et al.Fundamental Serrin type regularity criteria for 3D MHDfluid passing through the porous mediumJ.Applicable Analysis,2017,96:2130-2139.19 LIU Y,ZHANG T.Global weak solutions to a 2D compressible non-resistivity MHD system withnonmonotone pressure law and nonconstant viscosityJ.J.Math.Anal.Appl.,2021,502:125244.20 WU G,ZHANG Y,ZOU W.Optimal time-decay rates for the 3D compressible magnetohydrodynamicflows with discontinuous initial data and large oscillationsJ.J.Lond.Math.Soc.,2021,103:817-845.21 CHEMIN J,MCCORMICK D S,ROBINSON J C,et al.Local existence for the non-resistive MHD-equations in Besov spacesJ.Adv.Math.,2016,286:1-31.The Global Existence for a Class of MagnetohydrodynamicEquations with Nonlinear Damping TermLI Linrui1,HONG Mingli1,ZHENG Lin2(1.Basic Courses Department,Institute of Disaster Prevention,Sanhe 065201,China;2.School of Mathematics and Statistics,North China University of Water Resources andElectric Power,Zhengzhou 450046,China)Abstract:In this paper we consider the global existence results for a new class of magnetohydro-dynamic(MHD)equations with nonlinear damping term a|u|1u(a 0)in porous media,we obtain theglobal weak solutions and strong solution by classical energy method and Sobolev compactness embeddingmethod and the regularized results with the aid of Gagliardo-Nirenberg interpolation inequality and theothers important inequalities.The results of study improve and extend the ones in the previous worksto a large extent,reveal the physical phenomenon of fluid motion,and provide the necessary theoreticalbasis for the development of magnetohydrodynamic.Key words:Magnetohydrodynamic equations;Damping term;Viscous flow;Sobolev compactnessimbedding