[论文解读] Hyperbolic conservation laws on spacetimes. A finite volume scheme based on differential forms

Philippe G. LeFloch, Baver Okutmuştur|ArXiv.org|Oct 1, 2008

Geometric Analysis and Curvature Flows参考文献 26被引用 14

一句话总结

本文提出了一种基于微分形式的新型有限体积格式，用于时空上的双曲守恒律，将通量表述为n-形式场而非向量场。在全局双曲性条件下，建立了熵解半群的存在性，并证明了在几何自然的L¹型距离下具有收缩性质，将Kruzkov理论推广至曲率流形，且几何假设最小化。

ABSTRACT

We consider nonlinear hyperbolic conservation laws, posed on a differential (n+1)-manifold with boundary referred to as a spacetime, and in which the "flux" is defined as a flux field of n-forms depending on a parameter (the unknown variable). We introduce a formulation of the initial and boundary value problem which is geometric in nature and is more natural than the vector field approach recently developed for Riemannian manifolds. Our main assumption on the manifold and the flux field is a global hyperbolicity condition, which provides a global time-orientation as is standard in Lorentzian geometry and general relativity. Assuming that the manifold admits a foliation by compact slices, we establish the existence of a semi-group of entropy solutions. Moreover, given any two hypersurfaces with one lying in the future of the other, we establish a "contraction" property which compares two entropy solutions, in a (geometrically natural) distance equivalent to the L1 distance. To carry out the proofs, we rely on a new version of the finite volume method, which only requires the knowledge of the given n-volume form structure on the (n+1)-manifold and involves the {\sl total flux} across faces of the elements of the triangulations, only, rather than the product of a numerical flux times the measure of that face.

研究动机与目标

• 在一般微分流形上建立双曲守恒律的几何自然表述，避免依赖黎曼或洛伦兹度量。
• 通过基于微分形式的新有限体积方法，建立时空上守恒律熵解的存在性。
• 在几何内蕴距离（等价于L¹范数）下，证明熵解之间的收缩性质，确保唯一性与稳定性。
• 在不预先假设几何结构（除全局双曲性外）的前提下，将Kruzkov的熵解理论推广至曲率时空。
• 提供一种可数值实现的格式，仅需知道n-体积形式及单元面上的总通量，无需依赖度量的通量计算。

提出的方法

• 将守恒律表述为 $ d(\rho(u)) = 0 $，其中 $ \rho(u) $ 是依赖于未知量 $ u $ 的n-形式场（通量场），实现不依赖度量结构的几何表述。
• 通过在紧致子集边界上对测试形式积分，利用Stokes定理将守恒律表达为总通量消失，从而定义弱解。
• 利用时空的三角剖分构造有限体积格式，通过积分 $ \rho(u) $ 在单元面上的通量，而非使用数值通量乘以面测度。
• 引入一个几何一致的熵不等式，涉及离散解的Young测度极限，确保极限下熵的耗散。
• 利用全局双曲性条件定义时间方向，区分过去与未来，从而保证适定性与收缩估计。
• 通过在熵不等式中取极限，利用Young测度的紧致性与弱*收敛性，证明离散格式的收敛性。

实验结果

研究问题

• RQ1 如何在不假设黎曼或洛伦兹度量的前提下，于一般时空上以几何内蕴方式表述双曲守恒律？
• RQ2 能否构造一种仅依赖n-体积形式及面间总通量的有限体积方法，避免依赖度量的通量计算？
• RQ3 在何种几何条件下，熵解半群存在且在几何自然距离下满足收缩性质？
• RQ4 在无度量假设下，如何仅利用微分形式结构证明有限体积格式的收敛性？
• RQ5 在使用微分形式的前提下，Kruzkov的熵解理论在多大程度上可推广至曲率时空？

主要发现

• 本文在全局双曲性条件下，建立了时空上双曲守恒律熵解半群的存在性。
• 在等价于L¹范数的几何自然距离下，证明了两个熵解之间的收缩性质，确保唯一性。
• 当网格尺寸 $ h \to 0 $ 时，有限体积格式收敛至熵解，收敛性通过Young测度技术与熵不等式极限得以证明。
• 该格式仅需知道n-体积形式及单元面的总通量，适用于无度量假设的一般流形。
• 格式中的误差项 $ A^h, B^h, C^h $ 随 $ h \to 0 $ 而趋于零，且满足估计 $ |A^h| \leq C'' \frac{\tau_{\text{max}} + h}{(\tau_{\text{min}})^{1/2}} $ 与 $ |B^h| \leq C \frac{( au_{\text{max}} + h)^2}{h} $，两者均趋于零。
• 通过微分形式的表述，提供了一个比向量场方法更简洁、更通用的框架，统一并拓展了黎曼与洛伦兹情形下的先前结果。