[论文解读] Div-Curl Problems and Stream Functions in 3D Lipschitz Domains

Matthias Kirchhart, Erick Schulz|arXiv (Cornell University)|Jan 1, 2020

Advanced Numerical Methods in Computational Mathematics参考文献 1被引用 2

一句话总结

本文提出了一种在三维利普希茨区域中对速度场的旋度问题进行适定构造的方法，以求解散度-旋度问题，即从给定的涡度 F = curl U 恢复一个散度为零的速度场 U ∈ L²(Ω)。通过向量势 A ∈ H¹(Ω) 实现该解，使得 U = curl A，且在适当的泛函分析条件下证明了其存在性与唯一性。

ABSTRACT

We consider the problem of recovering the divergence-free velocity field $\mathbf{U}\in\mathbf{L}^2(\Omega)$ of a given vorticity $\mathbf{F}=\mathrm{curl}\,\mathbf{U}$ on a bounded Lipschitz domain $\Omega\subset \mathbb{R}^3$. To that end, we solve the 'div-curl problem' for a given $\mathbf{F}\in\bigl[ \mathbf{H}_0(\mathrm{curl};\Omega)\bigr]'$. The solution is given in terms of a vector potential (or stream function) $\mathbf{A}\in\mathbf{H}^1(\Omega)$ such that $\mathbf{U}=\mathrm{curl}\,{\mathbf{A}}$. After discussing existence and uniqueness of solutions and associated vector potentials, we propose a well-posed construction for the stream function. A numerical example of the construction is presented at the end.

研究动机与目标

• 解决在有界三维利普希茨区域 Ω 中，从其涡度 F = curl U 重构一个散度为零的速度场 U ∈ L²(Ω) 的逆问题。
• 为 F ∈ [H₀(curl; Ω)]′ 的散度-旋度问题建立解的存在性与唯一性。
• 构造一个适定的向量势 A ∈ H¹(Ω)，使得 U = curl A，确保解具有稳定性和物理意义。
• 提供一种可计算的框架，以数值上可行的方式计算流函数 A，并通过一个数值示例加以支持。

提出的方法

• 在对偶空间 [H₀(curl; Ω)]′ 中表述散度-旋度问题，以处理广义涡度场 F。
• 利用索伯列夫空间中的泛函分析工具，证明对于给定的 F，解 U ∈ L²(Ω) 的存在性与唯一性。
• 构造一个向量势 A ∈ H¹(Ω)，使得 U = curl A，从而保证 U 的散度为零。
• 为 A 建立一个适定的变分形式，保证在 H¹(Ω) 范数下的稳定性与收敛性。
• 采用伽辽金型方法或有限元法进行数值计算，并通过一个数值示例加以验证。

实验结果

研究问题

• RQ1 能否在三维利普希茨区域中，唯一地从其涡度 F = curl U 重构出一个散度为零的速度场 U ∈ L²(Ω)？
• RQ2 在什么条件下，可保证对于给定的 F ∈ [H₀(curl; Ω)]′，存在唯一一个向量势 A ∈ H¹(Ω)，使得 U = curl A？
• RQ3 如何在不规则的三维区域中实现流函数 A 的稳定且适定的构造？
• RQ4 对偶空间 [H₀(curl; Ω)]′ 在广义涡度的散度-旋度问题表述中起什么作用？
• RQ5 如何以保证收敛性与精度的方式数值计算流函数 A？

主要发现

• 建立了 A ∈ H¹(Ω) 的适定流函数构造，确保对于给定的 F ∈ [H₀(curl; Ω)]′，U = curl A 是 L²(Ω) 中的唯一解。
• 利用利普希茨区域上索伯列夫空间中的泛函分析方法，证明了解 U 的存在性与唯一性。
• 向量势 A 满足所需的正则性与边界条件，确保 U 为散度为零且属于 L²(Ω)。
• 该方法提供了一个稳定且可数值实现的框架，通过一个数值示例展示了构造过程。
• 该方法将经典流函数理论扩展至边界规则性较低的三维区域（如利普希茨区域），扩大了其在流体力学与偏微分方程理论中的适用范围。