微极Navier-Stokes方程的一种二阶时间步算法.pdf

Advances in Applied Mathematics 应用数学进展应用数学进展,2024,13(4),1362-1367 Published Online April 2024 in Hans.https:/www.hanspub.org/journal/aam https:/doi.org/10.12677/aam.2024.134127 文章引用文章引用:王思全.微极 Navier-Stokes 方程的一种二阶时间步算法J.应用数学进展,2024,13(4):1362-1367.DOI:10.12677/aam.2024.134127 微极微极Navier-Stokes方程的一种二阶时间步算法方程的一种二阶时间步算法 王思全王思全 云南师范大学数学学院，云南 昆明 收稿日期：2024年3月17日；录用日期：2024年4月11日；发布日期：2024年4月18日 摘摘 要要 针对微极针对微极Navier-Stokes方程方程(MNSE)，我们提出了一种新的二阶时间步算法。对，我们提出了一种新的二阶时间步算法。对MNSE中的非线性项进中的非线性项进行了线性化处理，并在线速度、压力和角速度的离散解中加入了“曲率稳定”项，旨在改进常用的“速行了线性化处理，并在线速度、压力和角速度的离散解中加入了“曲率稳定”项，旨在改进常用的“速率稳定”。该方法不仅克服了阻力产生的数值不稳定，并且在不增加计算复杂性的情况下将解的精度率稳定”。该方法不仅克服了阻力产生的数值不稳定，并且在不增加计算复杂性的情况下将解的精度从一阶提高到二阶。然后我们给出了算法的无条件稳定性。最后，通过数值实验验证了预测的从一阶提高到二阶。然后我们给出了算法的无条件稳定性。最后，通过数值实验验证了预测的收敛速收敛速度。度。关键词关键词 微极微极Navier-Stokes，二阶格式二阶格式，稳定性稳定性 A Second-Order Time-Step Algorithm for the Micropolar Navier-Stokes Equations Siquan Wang School of Mathematics,Yunnan Normal University,Kunming Yunnan Received:Mar.17th,2024;accepted:Apr.11th,2024;published:Apr.18th,2024 Abstract We propose a new second-order time-step algorithm for the micropolar Navier-Stokes equations(MNSE).The nonlinear term in MNSE is linearized,and the term“curvature stabilization”is added to the discrete solutions of online velocity,pressure and angular velocity,which aims to improve the commonly used“rate stabilization”.This method not only overcomes the numerical instability caused by resistance,but also improves the accuracy of the solution from the first order to the second order without increasing the computational complexity.And then we give the unconditional stability of the algorithm.Finally,the predicted convergence rate is verified by numerical experi-ments.王思全 DOI:10.12677/aam.2024.134127 1363 应用数学进展 Keywords Micropolar Navier-Stokes,Second-Order Scheme,Stability Copyright 2024 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 微极 Navier-Stokes 方程(MNSE)是一种重要的物理模型，该方程一般由 Navier-Stokes 方程耦合角动量方程而成，可以描述具有微观结构的流体运动现象，如材料质点具有平移和旋转自由度时的不可压缩流体的运动等。针对具有微观结构的流体模型的新型高性能的数值算法研究具有重要的理论和现实价值。1964 年，Eringen 提出了一种通用的微极流体理论来描述这类流体的运动，同时考虑了内部子结构的特性。文1证明了磁微极流体运动强解的存在性和唯一性。事实上，大量的研究人员使用数值方法和技术来讨论这类问题。例如：Nochetto 等人2提出了线速度和角速度解耦的 MNSE 一阶半隐式全离散有限元方法；Salgado 研究了微极 Navier-Stokes 方程(MNSE)的全离散分数阶时间步长投影有限元法3；文4针对Navier-Stokes 方程提出了一种新的、二阶精度的数值正则化方法；文5提出并分析了微极 Navier-Stokes方程(MNSE)的两种全离散投影数值格式；文6将 SAV 方法与压力投影方法相结合，提出了求解微极Navier-Stokes 方程的无条件能量稳定数值格式；文7提出了不可压缩微极性流体流动控制方程的二阶Gauge-Uzawa 格式；文8针对二维/三维非定常不可压缩热微极流体方程，提出了三种分裂有限元格式，分别为标准分裂格式、一致分裂格式。文9将对流项的标量辅助变量(SAV)方法和耦合项的隐显(IMEX)处理方法相结合，提出了求解微极 Navier-Stokes 方程的解耦、线性和无条件能量稳定数值格式。主要研究内容是针对微极 Navier-Stokes 方程提出一种新的、二阶无条件稳定的时间步方法，给出该格式的无条件稳定性。最后用数值实验验证该格式的收敛速度。2.符号和符号和预备知识预备知识 微极 Navier-Stokes 方程可以描述如下：对于有界正则域()23ddor=，给定最终时间T，找到(:0,:0,ddRpdR，u和(:0,ddRw满足：()()()()()02042trrtaddarrsspsjccjcccss+=+=+=+uuuuwfuwwuwwwug (1)边界条件为：()()()(,0,0,0,x tx tx tT=uw，初始条件为：()()00,0,.,0 xxx=uu ww。这里的,pu和w分别是线速度、压力和角速度。f和g的函数分别表示平滑的外部施加的力和力矩。我们假设所有材料常数0,raj s s c c和dc的运动粘度均为常数且为正数。我们设 0120,radadsss ccccccc=+=+=+，并且12,0c c。特别的，当2d=时，我们令()()12,0,0,0,u uw=uw，那么有 210 0,0uuwwxyyx=，uw Open AccessOpen Access王思全 DOI:10.12677/aam.2024.134127 1364 应用数学进展 我们令()pL空间中的范数表示为pL，()pkW空间中的范数表示为pkW。特别的，()2dL空间中的范数表示为，()2dL空间中的范数表示为。接下来，我们定义一些 Sobolev 空间：()()()()()1102201:0:0:0dddXHHQLqLqdxYHn=vvzz 此外，无散函数的空间由下式给出：():,0,VXqqQ=vv 对于在整个时间间隔()0,T上定义的函数(),x tu，我们定义以下范数：()()()()()10,;0,;00ess sup,1pTppLT XLT XXXt Tttdtp=uuuu 令1V表示V的对偶空间，1V的范数由下式给出()10s,upv V=f vfv 我们定义离散无散空间：():,0,hhhhhhhVXqqQ=vv 通常的显式偏斜对称三线性形式定义为：()()()*11,22b=u v wuv wuw v 我们可以将(1)的弱形式表示为：找到(:0,:0,TX pTQu和(:0,TXw，在时间(0,tT上满足：()()()()()()()()()()()()()()*0*12,2,0,4,2,trtrrsbpsqjcjbcss+=+=+=+u vuvu u vvw vf vuw zwzu w zwzw zu zg z 3.算法及其稳定性算法及其稳定性 提出了一种二阶全离散时间步进算法，该算法将非线性项的线性处理与离散弱解在时间上的稳定项相结合，并给出了该方法的无条件稳定性。3.1.二阶时间步算法二阶时间步算法 我们定义最终时间T、时间步数N、时间步TtN=和,1,2.ntn t nN=。1ntt+=处的插值由下 式给出：()1112:1nnnnaaaa+=+王思全 DOI:10.12677/aam.2024.134127 1365 应用数学进展 1ntt+=处的外推由下式给出：()11:2nnnaaa+=接下来我们给出二阶全离散时间步算法。算法算法 3.1 给定11,nnnnhhhhuuww和1,nnhhpp，找到()()111,0,1,1nnnhhhhhhpXQXnN+=uw满足：()()()()()()()()()()()()()()()()()()11011111111111112134,2,2,034,2,4nnnhhhhnhhnhnhhnhhnrnhhhnhnnnhhhhnhnhhnhhnhhrstbpsqjjbtccs+=+=+uuuvuvuuvvwvfvuwwwzuwzwzwz()()()()()111,2,nhhnrnhhhs+=+gwzuzz (2)3.2.稳定性分析稳定性分析 为了分析稳定性，我们通过n nF引入对称矩阵4FII=+，以及对称矩阵22nnG如下：5242224IIGII+=+我们定义F范数和G分别为：()2,n nFGFG=uuuuuuu wwww 定理定理 3.2 算法 3.1 无条件稳定并且满足以下界限：()()()()()()()1222221112111221111122112200001221111112212333332142NNNhhnhnhnhnNnnnnnnhhhhhhFFnhhhhhhGGNnnnNtjsccjNjjN tsc+=+=+=+fguwuwwuuuwwwuwuwuw 4.数值实验数值实验 在这一部分中，我们使用了两个数值实验来验证所提出的方法的理论结果。我们用 BDF2Reg 测试了王思全 DOI:10.12677/aam.2024.134127 1366 应用数学进展 算法 3.1 的收敛速度。我们使用 FreeFem+完成以下实验。在这个实验中，我们的目的是验证理论结果的收敛速度。我们考虑域20,1=，设置其他参数01radjssccc=。我们假设已知的函数()()()()()()()()()4323213243221210223cos,10 232cos,10 2121 cos,xxxyyytux y tuxxxyyytp x y txytw x y tuu+=+=+u 作为真解。对于线速度、压力和角速度，我们使用混合有限元212PPP。误差估计定义如下：()()1/222,11,Nnnhhntt=vvvvvu w 我们通过改变网格宽度h和时间步长t来测试二阶全离0.1,/16,0.1TtT=散时间步进算法3.1 的精度。对于表 1 中的空间收敛阶数，我们取最终时间和不同的网格宽度。对于表 2 中的时间收敛顺序，我们采用最终时间3,128,0.1ThT=和变化的时间步长。我们可以看到，空间收敛和时间收敛的顺序都为 2。因此，该实验验证了理论的正确性。Table 1.Order of spatial convergence 表表 1.空间收敛阶数 h 2,1huu 速率 2,1hww 速率 1/2 0.0150004 0.0086915 1/4 0.00402464 1.89807 0.00239977 1.85671 1/8 0.00102998 1.96624 0.00061566 1.96269 1/16 0.000259289 1.98999 0.000154974 1.99011 1/32 6.49757e-005 1.99659 3.88132e005 1.99741 Table 2.Order of time convergence 表表 2.时间收敛阶数 t 2,1huu 速率 2,1hww 速率 T/4 0.0198007 0.0197803 T/8 0.00569487 1.79782 0.00568793 1.79809 T/16 0.00144155 1.98204 0.00143825 1.98359 T/32 0.000370322 1.96078 0.00036441 1.98068 5.结语结语 本文提出了微极 Navier-Stokes 方程的一种新的二阶全离散时间步算法，我们对非线性项进行线性化处理，在时间上增加离散弱解的稳定项，在空间上增加协调有限元。然后给出了该算法的无条件稳定性，并通过数值试验验证了所提出算法的理论结果。本文提出的算法将解的精度从一阶提高到了二阶，并且可以轻松应用于许多其他模型，是研究这些线性化方法和 CFD 背景下新的稳定性的初步尝试。稳定项如何单独影响稳定性仍需要更多研究。一般来说，需要更广泛的数值实验，特别是对于存在间歇性和大时间变化的问题。王思全 DOI:10.12677/aam.2024.134127 1367 应用数学进展 参考文献参考文献 1 John,V.,Knobloch,P.and Novo,J.(2018)Finite Elements for Scalar Convection-Dominated Equations and Incom-pressible Flow Problems:A Never Ending Story?Computing and Visualization in Science,19,47-63.https:/doi.org/10.1007/s00791-018-0290-5 2 Nochetto,R.H.,et al.(2014)The Micropolar Navier-Stokes Equations:A Priori Error Analysis.Mathematical Models&Methods in Applied Sciences,24,1237-1264.https:/doi.org/10.1142/S0218202514500018 3 Salgado,A.J.(2015)Convergence Analysis of Fractional Time-Stepping Techniques for Incompressible Fluids with Microstructure.Journal of Scientific Computing,64,216-233.https:/doi.org/10.1007/s10915-014-9926-x 4 Jiang,N.,Mohebujjaman,M.,Rebholz,L.and Trenchea,C.(2016)An Optimally Accurate Discrete Regularization for Second Order Time Stepping Methods for Navier-Stokes Equations.Computer Methods in Applied Mechanics and En-gineering,310,388-405.https:/doi.org/10.1016/j.cma.2016.07.017 5 Yang,Y.-B.and Jiang,Y.-L.(2017)Analysis of Two Decoupled Time-Stepping Finite Element-Methods for Incom-pressible Fluids with Microstructure.International Journal of Computer Mathematics,95,686-709.https:/doi.org/10.1080/00207160.2017.1294688 6 Shen,J.and Zheng,N.(2022)Efficient and Unconditional Energy Stable Schemes for the Micropolar Navier-Stokes Equations.CSIAM Transactions on Applied Mathematics,3,57-81.https:/doi.org/10.4208/csiam-am.SO-2021-0008 7 Slayi,S.,Arwadi,T.E.and Dib,S.(2021)Stabilized Gauge Uzawa Scheme for an Incompressible Micropolar Fluid Flow.Applied Numerical Mathematics,167,45-72.https:/doi.org/10.1016/j.apnum.2021.04.003 8 Ren,Y.and Liu,D.(2023)Research on Three Kinds of Splitting Finite Element Schemes for 2D/3D Unsteady Incom-pressible Thermomicropolar Fluid Equations.International Journal for Numerical Methods in Fluids,95,1148-1173.https:/doi.org/10.1002/fld.5188 9 Zhang,X.-D.,Long,X.-N.(2023)Unconditional Stability and Error Analysis of an Euler IMEX-SAV Scheme for the Micropolar Navier-Stokes Equations.Applied Numerical Mathematics,192,214-240.https:/doi.org/10.1016/j.apnum.2023.05.027