留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

关于带自由表面的Navier-Stokes方程的一点注记

上一篇

下一篇

田清霞, 朱朝生. 关于带自由表面的Navier-Stokes方程的一点注记[J]. 西南师范大学学报(自然科学版), 2020, 45(3): 21-26. doi: 10.13718/j.cnki.xsxb.2020.03.004
引用本文: 田清霞, 朱朝生. 关于带自由表面的Navier-Stokes方程的一点注记[J]. 西南师范大学学报(自然科学版), 2020, 45(3): 21-26. doi: 10.13718/j.cnki.xsxb.2020.03.004
Qing-xia TIAN, Chao-sheng ZHU. A Note on Navier-Stokes Equation with Free Surface[J]. Journal of Southwest China Normal University(Natural Science Edition), 2020, 45(3): 21-26. doi: 10.13718/j.cnki.xsxb.2020.03.004
Citation: Qing-xia TIAN, Chao-sheng ZHU. A Note on Navier-Stokes Equation with Free Surface[J]. Journal of Southwest China Normal University(Natural Science Edition), 2020, 45(3): 21-26. doi: 10.13718/j.cnki.xsxb.2020.03.004

关于带自由表面的Navier-Stokes方程的一点注记

  • 基金项目: 国家自然科学基金项目(11571283)
详细信息
    作者简介:

    田清霞(1994-), 女, 硕士研究生, 主要从事偏微分方程现代理论研究 .

    通讯作者: 朱朝生, 副教授
  • 中图分类号: O175.29

A Note on Navier-Stokes Equation with Free Surface

  • 摘要: 利用文中的一些定义和复合函数的求导方法以及Newton-Leibniz公式,对η函数(一个关键性函数)的一阶和二阶偏导进行了详细推导.
  • 加载中
  • [1] 薛菊峰, 尚月强.非定常不可压Navier-Stokes方程基于欧拉格式的两水平变分多尺度方法[J].西南大学学报(自然科学版), 2018, 40(9):84-90. doi: http://d.old.wanfangdata.com.cn/Periodical/xnnydxxb201809013
    [2] 周海军, 高真圣.Navier-Stokes-Poisson方程的两个注记[J].贵州师范大学学报(自然科学版), 2016, 34(4):54-57. doi: 10.3969/j.issn.1004-5570.2016.04.010
    [3] MASMOUDI N, ROUSSET F.Uniform Regularity and Vanishing Viscosity Limit for the Free Surface Navier-Stokes Equations[J].Archive for Rational Mechanics and Analysis, 2017, 223(1):301-417. doi: 10.1007/s00205-016-1036-5
    [4] TAREK M E, DONGHYUN L.Uniform Regularity for Free-Boundary Navier-Stokes Equations with Surface Tension[J].arXiv: Analysis of PDEs, 2014.
    [5] LEE D.Uniform Estimate of Viscous Free-Boundary Magnetohydrodynamics with Zero Vacuum Magnetic Field[J].SIAM Journal on Mathematical Analysis, 2017, 49(4):2710-2789. doi: 10.1137/16M1089794
    [6] 吴复洲.旋量层与带自由表面的Navier-Stokes方程组的粘性消失极限[D].北京: 清华大学, 2016.http://cdmd.cnki.com.cn/Article/CDMD-10003-1018811772.htm
    [7] doi: http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ0229213859/ MEI Y, WANG Y, XIN Z P.Uniform Regularity for the Free Surface Compressible Navier-Stokes Equations with or Without Surface Tension[J].Mathematical Models and Methods in Applied Sciences, 2018, 28(2):259-336.
    [8] doi: http://cn.bing.com/academic/profile?id=d277cc49558e0f21ecd2f9d37ee91842&encoded=0&v=paper_preview&mkt=zh-cn GUO Y, TICE I.Local Well-Posedness of the Viscous Surface Wave Problem Without Surface Tension[J].Analysis & PDE, 2013, 6(2):287-369.
    [9] BEALE J T.The Initial Value Problem for the Navier-Stokes Equations with a Free Surface[J].Communications on Pure and Applied Mathematics, 1981, 34(3):359-392. doi: 10.1002/cpa.3160340305
    [10] 华东师范大学数学系编.数学分析(下册)[M].4版.北京: 高等教育出版社, 2010: 185-188.
  • 加载中
计量
  • 文章访问数:  605
  • HTML全文浏览数:  605
  • PDF下载数:  93
  • 施引文献:  0
出版历程
  • 收稿日期:  2019-01-03
  • 刊出日期:  2020-03-20

关于带自由表面的Navier-Stokes方程的一点注记

    通讯作者: 朱朝生, 副教授
    作者简介: 田清霞(1994-), 女, 硕士研究生, 主要从事偏微分方程现代理论研究
  • 西南大学数学与统计学院, 重庆 400715
基金项目:  国家自然科学基金项目(11571283)

摘要: 利用文中的一些定义和复合函数的求导方法以及Newton-Leibniz公式,对η函数(一个关键性函数)的一阶和二阶偏导进行了详细推导.

English Abstract

  • 近年来,Navier-Stokes方程的相关研究已经引起了广泛的关注[1-2].文献[3]考虑了不可压缩粘性流体在高Reynolds数下的运动方程(不可压缩Navier-Stokes方程):

    其中流体域Ωt={x$\mathbb{R} $3x3 < h(tx1x2)}.文中通过一个微分型同胚映射Φ(t,·)将运动域上的问题简化到固定域S上,

    其中

    这里产生的η函数是h的延拓,即$\hat{\eta }(\xi, z)=\chi (z\xi)+\hat{h}(\xi)$,其中$\hat{\centerdot }$表示关于变量y的Fourier变换,χ是一个光滑的紧支集函数,使得在B(0,1)上χ=1.另外常数A>0,这使得zφ>0.文献[4-9]都研究了这个η函数.文献[3]中多次对φ进行求导,因此不可避免地必须对η函数进行求导,而η函数的求导过程较为繁琐,并且文献[3]并没有给出具体的计算过程,但是必须有η函数的求导结果才有后续的证明结果.因此本文专门讨论了η函数的一阶偏导和二阶偏导.

  • $F(t)=\int_{{{x}^{2}}+{{y}^{2}}\le {{t}^{2}}}{f}(x, y)\text{d}x\text{d}y$,首先计算F′(t),F″(t).为此作变换$\left\{ \begin{array}{*{35}{l}} x=r\cos \theta \\ y=r\sin \theta \\ \end{array}, 0\le r\le |t|, 0\le \theta \le2\text{ }\!\!\pi\!\!\text{ } \right.$,则有

    一般地,我们定义$G(r)=\int_{0}^{2\text{ }\!\!\pi\!\!\text{ }}{f}(r\cos \theta, r\sin \theta)r\text{d}\theta $.

    1) 当t>0时,rt,则$F(t)=\int_{0}^{t}{G}(r)\text{d}r$,所以

    2) 类似地,当t < 0时,r≤-t,则$F(t)=\int_{0}^{-t}{G}(r)\text{d}r$,所以

    3) 因为当t=0时,F(0)=0,所以当t→0+时,

    类似地,当t→0_时,易得F′_(0)=0.

  • 下面引入Fourier变换和Fourier逆变换的定义:

    由文献[4-6]知

    其中:y=(y1y2),ξ=(ξ1ξ2).这里$\hat{\centerdot }$表示关于变量y的Fourier变换,R表示圆域的半径,χ是一个光滑的紧支集函数,使得在B(0,1)上χ=1.

    下面讨论Ψ(tyz).由于1≤|zξ|≤R,所以1≤z2ξ12+z2ξ22R2,则$\frac{1}{{{z}^{2}}}\le \xi _{1}^{2}+\xi _{2}^{2}\le \frac{{{R}^{2}}}{{{z}^{2}}}$.因此作变换$\left\{ \begin{array}{*{35}{l}} {{\xi }_{1}}=r\cos \theta \\ {{\xi }_{2}}=r\sin \theta \\ \end{array}, \frac{1}{|z|}\le r\le \frac{R}{|z|}, 0\le \theta \le 2\text{ }\!\!\pi\!\!\text{ } \right.$.为了方便记:

    则有

    其中$f(t, {\mathit{\boldsymbol{y}}}, z, r)=\int_{0}^{2\text{ }\!\!\pi\!\!\text{ }}{\chi }(zr\cos \theta, zr\sin \theta)\hat{h}(t, r\cos \theta, r\sin \theta){{\text{e}}^{2\text{i }\!\!\pi\!\!\text{ }r\left({{y}_{1}}\cos \theta +{{y}_{2}}\sin \theta \right)}}r\text{d}\theta $.

    命题1   当z>0时,

      当z>0时,$\frac{1}{z}\le r\le \frac{R}{z}$,则$\mathit{\Psi} (t, {\mathit{\boldsymbol{y}}}, z)=\int_{\frac{1}{z}}^{\frac{R}{z}}{f}(t, {\mathit{\boldsymbol{y}}}, z, r)\text{d}r$.所以,由文献[10]的Newton-Leibniz公式可得:

    进一步求导得:

    由此容易得到(8)式.同理可得方程(9)和(10).下面计算Ψzz(tyz),

    由此容易得到(11)式.

    下面讨论当z < 0时的情形,为了方便记:

    命题2   当z < 0时,

      当z < 0时,$-\frac{1}{z}\le r\le -\frac{R}{z}$,则$\mathit{\Psi} (t, {\mathit{\boldsymbol{y}}}, z)=\int_{-\frac{1}{z}}^{-\frac{R}{z}}{f}(t, {\mathit{\boldsymbol{y}}}, z, r)\text{d}r$.所以易得

    进一步,采用命题1的方法进行计算,很容易得到方程(13)-(16).

    其次,讨论Φ(tyz).由于|zξ|≤1,所以z2ξ12+z2ξ22≤1,则$\xi _{1}^{2}+\xi _{2}^{2}\le \frac{1}{{{z}^{2}}}=\frac{1}{|z{{|}^{2}}}$.因此作变换$\left\{ \begin{array}{*{35}{l}} {{\xi }_{1}}=r\cos \theta \\ {{\xi }_{2}}=r\sin \theta \\ \end{array}, 0\le r\le \frac{1}{|z|}, 0\le \theta \le 2\text{ }\!\!\pi\!\!\text{ } \right.$,则有

    其中$G(t, {\mathit{\boldsymbol{y}}}, r)=\int_{0}^{2\text{ }\!\!\pi\!\!\text{ }}{{\hat{h}}}(t, r\cos \theta, r\sin \theta){{\text{e}}^{2\text{i }\!\!\pi\!\!\text{ }\ r\left({{y}_{1}}\cos \theta +{{y}_{2}}\sin \theta \right)}}r\text{d}\theta $.当z→0时,$\mathit{\Phi} (t, {\mathit{\boldsymbol{y}}}, z)=\int_{{{\mathbb{R}}^{2}}}{{\hat{h}}}(t, {\mathit{\boldsymbol{ }}}\!\!\xi\!\!\text{ }){{\text{e}}^{2\text{i }\pi \text{ }\mathbf{y}\times \text{ }\xi \text{ }}}\text{d}\xi \text{ }$,容易计算各阶偏导数.下面分别讨论z>0和z < 0的情形.为了方便,记$\hat{h}(t, z\cos \theta, z\sin \theta):=\hat{h}, \ {{\text{e}}^{2\text{i }\!\!\pi\!\!\text{ }\ z\left({{y}_{1}}\cos \theta +{{y}_{2}}\sin \theta \right)}}:={\tilde{\rm e}}$.

    命题3   1)当z < 0时,

    2) 当z < 0时,

      当z>0时,$\mathit{\Phi} (t, {\mathit{\boldsymbol{y}}}, z)=\int_{0}^{\frac{1}{z}}{G}(t, {\mathit{\boldsymbol{y}}}, r)\text{d}r$,则$\mathit{\Phi} _{z}^{\prime }(t, {\mathit{\boldsymbol{y}}}, z)=-\frac{1}{{{z}^{2}}}G(t, {\mathit{\boldsymbol{y}}}, z)$.又因为

    所以$\Phi _{z}^{\prime }(t, {\mathit{\boldsymbol{y}}}, z)=-\frac{1}{z}\int_{0}^{2\text{ }\!\!\pi\!\!\text{ }}{{\hat{h}}}{\tilde{\rm e}d}\theta $.进一步求导可得方程(19)-(22).

    z < 0时,$\mathit{\Phi} (t, {\mathit{\boldsymbol{y}}}, z)=\int_{0}^{-\frac{1}{z}}{G}(t, {\mathit{\boldsymbol{y}}}, r)\text{d}r$,则$\mathit{\Phi} _{z}^{\prime }(t, {\mathit{\boldsymbol{y}}}, z)=\frac{1}{{{z}^{2}}}G(t, {\mathit{\boldsymbol{y}}}, z)$.同理可得方程(23)-(26).

参考文献 (10)

目录

/

返回文章
返回