[论文解读] Regularized Stokes Immersed Boundary Problems in Two Dimensions: Well-posedness, Singular Limit, and Error Estimates

Jiajun Tong|arXiv (Cornell University)|Apr 21, 2019

Lattice Boltzmann Simulation Studies参考文献 39被引用 4

一句话总结

本文提出了一种使用光滑化δ函数来平滑奇异力和流体耦合的二维Stokes浸入边界问题的正则化版本。建立了全局适定性，证明了当正则化参数趋于零时，正则化弦动力学收敛到未正则化问题，并推导出误差估计，表明当选择合适的光滑化函数时可实现更高的精度。

ABSTRACT

Inspired by the numerical immersed boundary method, we introduce regularized Stokes immersed boundary problems in two dimensions to describe regularized motion of a 1-D closed elastic string in a 2-D Stokes flow, in which a regularized $δ$-function is used to mollify the flow field and singular forcing. We establish global well-posedness of the regularized problems, and prove that as the regularization parameter diminishes, string dynamics in the regularized problems converge to that in the Stokes immersed boundary problem with no regularization. Viewing the un-regularized problem as a benchmark, we derive error estimates under various norms for the string dynamics. Our rigorous analysis shows that the regularized problems achieve improved accuracy if the regularized $δ$-function is suitably chosen. This may imply potential improvement in the numerical method, which is worth further investigation.

研究动机与目标

• 对使用光滑化δ函数的二维Stokes浸入边界问题的正则化连续形式进行严格分析。
• 建立二维Stokes流中一维弹性弦的正则化系统的全局适定性。
• 证明当正则化参数趋于零时，正则化弦动力学收敛到未正则化问题。
• 推导正则化与未正则化解之间误差估计，其依赖于正则化参数和光滑化函数的选择。
• 基于分析得到的误差界，提出对数值浸入边界方法的潜在改进。

提出的方法

• 使用光滑、紧支集的光滑化函数替代Dirac δ函数，构建正则化的二维Stokes IB问题，以正则化奇异力和速度场。
• 通过消除流体变量，推导弦运动的轮廓动力方程，将弦速度表示为涉及二维Stokes方程基本解的奇异积分。
• 在周期Sobolev空间中使用Hölder型与Sobolev型估计，控制弦构型及其导数的差异。
• 结合特征线法与积分估计，控制正则化与未正则化弦动力学随时间的差异。
• 对速度核差异应用精确的点态估计与$L^2$型估计，以控制弦演化中的误差。
• 推导出依赖于正则化参数$ρ$、光滑化函数特性及解正则性的各种范数（如$L^2$、$Ϫ{H}^1$）下的误差界。

实验结果

研究问题

• RQ1 正则化的二维Stokes浸入边界问题是否对$H^{5/2}$类初始弦构型存在全局解？
• RQ2 当正则化参数$ρ \to 0$时，正则化问题的解是否收敛到未正则化问题的解？
• RQ3 正则化与未正则化弦动力学之间的误差估计如何依赖于光滑化函数的选择？
• RQ4 正则化形式能否实现优于标准数值浸入边界方法的精度？
• RQ5 何种正则化δ函数能最小化弦动力学中的误差？

主要发现

• 对于$H^{5/2}(\bbT)$类初始弦构型，正则化问题在全局范围内是适定的，确保所有时间下解的存在性与唯一性。
• 当正则化参数$ρ \to 0$时，正则化问题中的弦动力学在$L^2$与$Ϫ{H}^1$范数下收敛到未正则化问题的解。
• 推导出误差估计，表明弦动力学差异被界于$C\lambda^{-2}\|Y_{1}^{\prime\prime}\|_{\dot{H}^\theta}^2$与高阶项之间，其中$λ$控制光滑化函数的宽度。
• 误差界对光滑化函数的选择极为敏感：构造更优的光滑化函数（如具有更高阶光滑性者）可提升收敛速率。
• 分析表明，正则化形式可实现优于标准数值方法的精度，提示了方法改进的路径。
• 进一步推导出$Ϫ{H}^1$型范数下的精细误差估计：$\|g_{Y_1} - g_{Y_2}\|_{\dot{H}^1} \leq C\lambda^{-3}(\|Y_1\|_{\dot{H}^{2+\theta}} + \|Y_2\|_{\dot{H}^{2+\theta}})^3\|\delta Y\|_{\dot{H}^{3/2}}$，其精度优于较弱的替代估计。