[论文解读] Homogenization of parabolic equations with a continuum of space and time scales

Houman Owhadi, Lei Zhang|ArXiv.org|Dec 21, 2005

Advanced Mathematical Modeling in Engineering参考文献 37被引用 5

一句话总结

本论文提出了一种针对系数不具尺度分离性或遍历性的线性抛物型方程的均匀化框架，通过调和坐标实现正则性，利用Cordes型条件。仅需求解 $ n $ 次，即可实现时空均匀化，将计算成本降低高达35,000倍，从而高效求解地球物理学和油藏模拟中的多尺度问题。

ABSTRACT

This paper addresses the issue of homogenization of linear divergence form parabolic operators in situations where no ergodicity and no scale separation in time or space are available. Namely, we consider divergence form linear parabolic operators in $Ω\subset \R^n$ with $L^\infty(Ω imes (0,T))$-coefficients. It appears that the inverse operator maps the unit ball of $L^2(Ω imes (0,T))$ into a space of functions which at small (time and space) scales are close in $H^1$-norm to a functional space of dimension $n$. It follows that once one has solved these equations at least $n$-times it is possible to homogenize them both in space and in time, reducing the number of operations counts necessary to obtain further solutions. In practice we show that under a Cordes type condition that the first order time derivatives and second order space derivatives of the solution of these operators with respect to harmonic coordinates are in $L^2$ (instead of $H^{-1}$ with Euclidean coordinates). If the medium is time independent then it is sufficient to solve $n$ times the associated elliptic equation in order to homogenize the parabolic equation.

研究动机与目标

• 解决空间与时间上缺乏尺度分离或遍历性的抛物型PDE的均匀化问题。
• 开发一种方法，对具有 $ L^∞ $ 系数的线性散度型抛物型算子，同时均匀化空间与时间尺度。
• 降低重复求解多尺度抛物型问题的计算成本，尤其适用于油藏模拟与油气开采。
• 在调和坐标下，通过Cordes型条件建立解的正则性，从而实现有效的尺度提升。
• 证明在时间不变情况下，仅需求解 $ n $ 次椭圆问题即可实现抛物型算子的均匀化。

提出的方法

• 通过求解一个边界条件为 $ F = x $、且在 $ \Omega $ 内满足 $ \operatorname{div}(a(x,0)\nabla F) = 0 $ 的抛物型系统，引入调和坐标。
• 定义张量 $ \sigma = {}^{t}\nabla F a \nabla F $，并利用Cordes参数 $ \beta_\sigma $ 衡量各向异性和确保正则性。
• 在 $ \beta_\sigma < 1 $ 且 $ \sigma $ 的迹有界时，证明 $ u \circ F^{-1} \in L^2(0,T; W^{2,2}_D(\Omega)) $。
• 使用层级矩阵方法，以 $ O(N(\ln N)^{n+3}) $ 次运算求解 $ n $ 重椭圆问题，其中 $ N $ 为细尺度自由度。
• 在粗网格上构建均匀化抛物型算子，自由度为 $ N^\alpha $ 或固定值 $ M $，实现高达35,000倍的压缩比。
• 通过求解 $ n $ 次并利用所得结构，高效近似 $ K > n $ 个时间依赖介质的解。

实验结果

研究问题

• RQ1 在空间与时间上均无尺度分离或遍历性的情况下，能否实现抛物型方程的均匀化？
• RQ2 调和坐标变换下的Cordes型条件是否足以保证正则性，从而实现有效的尺度提升？
• RQ3 在时间不变情况下，是否仅通过求解 $ n $ 次关联的椭圆问题即可近似抛物型PDE的解？
• RQ4 在多尺度抛物型问题中同时均匀化空间与时间尺度时，可实现的计算成本压缩因子是多少？
• RQ5 该方法在具有不连续、多尺度系数的随机分形介质上的表现如何？

主要发现

• 对于时间不变系数，求解 $ n $ 次椭圆问题即可实现抛物型方程的均匀化，自由度从 $ N $ 降低至 $ N^\alpha $，其中 $ \alpha = 0.2 $。
• 在数值实验中，当 $ N = 16,641 $ 且 $ M = 9 $ 个粗网格自由度时，该方法将计算复杂度降低了高达35,000倍。
• 在具有随机傅里叶模的时间依赖情况下，仅需67个粗时间步长，而非1,332个细时间步长，实现时间步长减少20倍，总复杂度降低35,000倍。
• 对于 $ \lambda_{\max}/\lambda_{\min} = 160.3 $ 的时间依赖随机分形介质，细时间步长从3,482减少至175，压缩因子达35,000倍。
• 在所有测试案例中，粗网格的 $ L^1 $、$ L^2 $、$ L^\infty $ 和 $ H^1 $ 误差均低于0.007，且 $ M = 9 $ 个节点，表明具有高精度。
• 在某些配置下，粗网格的 $ H^1 $ 误差始终低于细网格，表明均匀化带来了更高的稳定性。