Carasso Alfred S
Applied and Computational Mathematics Division, National Institute of Standards and Technology, Gaithersburg, MD 20899.
Inverse Probl Sci Eng. 2019;27(12). doi: 10.1080/17415977.2018.1523905.
This paper constructs an unconditionally stable explicit difference scheme, marching backward in time, that can solve a limited, but important class of time-reversed 2D Burgers' initial value problems. Stability is achieved by applying a compensating smoothing operator at each time step to quench the instability. This leads to a distortion away from the true solution. However, in many interesting cases, the cumulative error is sufficiently small to allow for useful results. Effective smoothing operators based on (-Δ) , with real > 2, can be efficiently synthesized using FFT algorithms, and this may be feasible even in non-rectangular regions. Similar stabilizing techniques were successfully applied in other ill-posed evolution equations. The analysis of numerical stabilty is restricted to a related linear problem. However, extensive numerical experiments indicate that such linear stability results remain valid when the explicit scheme is applied to a significant class of time-reversed nonlinear 2D Burgers' initial value problems. As illustrative examples, the paper uses fictitiously blurred 256 × 256 pixel images, obtained by using sharp images as initial values in well-posed, forward 2D Burgers' equations. Such images are associated with highly irregular underlying intensity data that can seriously challenge ill-posed reconstruction procedures. The stabilized explicit scheme, applied to the time-reversed 2D Burgers' equation, is then used to deblur these images. Examples involving simpler data are also studied. Successful recovery from severely distorted data is shown to be possible, even at high Reynolds numbers.
本文构造了一种无条件稳定的显式差分格式,该格式在时间上向后推进,可用于求解一类有限但重要的二维Burgers方程时间反演初值问题。通过在每个时间步应用补偿平滑算子来抑制不稳定性,从而实现稳定性。这会导致解偏离真实解。然而,在许多有趣的情况下,累积误差足够小,能够得到有用的结果。基于(-Δ)且实数>2的有效平滑算子可以使用快速傅里叶变换(FFT)算法有效地合成,甚至在非矩形区域也是可行的。类似的稳定技术已成功应用于其他不适定演化方程。数值稳定性分析仅限于一个相关的线性问题。然而,大量的数值实验表明,当显式格式应用于一类重要的二维Burgers方程时间反演非线性初值问题时,这种线性稳定性结果仍然有效。作为示例,本文使用了通过将清晰图像作为适定的正向二维Burgers方程的初值得到的虚拟模糊256×256像素图像。此类图像与高度不规则的底层强度数据相关联,这会严重挑战不适定的重建过程。然后,将应用于二维Burgers方程时间反演的稳定显式格式用于对这些图像进行去模糊处理。还研究了涉及更简单数据的示例。结果表明,即使在高雷诺数下,从严重失真的数据中成功恢复也是可能的。