退化抛物型方程的一个初值反演问题

数学物理学报 ›› 2020, Vol. 40 ›› Issue (4): 891-903.

退化抛物型方程的一个初值反演问题

杨柳¹(),邓醉茶^1,^2,*()

¹ 兰州交通大学数理学院兰州 730070
² 兰州交通大学计算机科学与技术国家级教学示范中心兰州 730070

收稿日期:2018-11-29 出版日期:2020-08-26 发布日期:2020-08-20
通讯作者: 邓醉茶 E-mail:l_yang218@163.com;zc_deng78@hotmail.com
作者简介:杨柳, E-mail:l_yang218@163.com
基金资助:
国家自然科学基金(11461039);国家自然科学基金(61663018);国家自然科学基金(11961042);兰州交通大学百名青年优秀人才培养计划和甘肃省自然科学基金(18JR3RA122)

An Inverse Initial Value Problem for Degenerate Parabolic Equations

Liu Yang¹(),Zuicha Deng^1,^2,*()

¹ Department of Mathematics, Lanzhou Jiaotong University, Lanzhou 730070
² Computer Science and Technology Experimental Teaching Center, Lanzhou Jiaotong University, Lanzhou 730070

Received:2018-11-29 Online:2020-08-26 Published:2020-08-20
Contact: Zuicha Deng E-mail:l_yang218@163.com;zc_deng78@hotmail.com
Supported by:
NSFC(11461039);NSFC(61663018);NSFC(11961042);the Foundation of a Hundred Youth Talents Training Program of Lanzhou Jiaotong University and the NSF of Gansu Province(18JR3RA122)

摘要/Abstract

摘要：

研究了一类重构退化抛物型方程初值的反问题.这类问题在应用科学的若干领域有着重要的应用.数值求解该问题的关键是构造相应正问题的高阶差分格式.然而，由于退化边界上的主项系数为零，目前广泛用于求解经典热传导方程的虚拟点法不能应用于该模型.该文提出了一种构造二阶精度差分格式的新方法，并证明了该方法的稳定性和收敛性.为了加快收敛速度，采用共轭梯度法求逆问题的数值解，并对算法的效率和精度进行了数值验证.

关键词: 退化抛物型方程, 反初值问题, 共轭梯度法, 收敛性, 数值结果

Abstract:

This paper investigates an inverse problem of reconstructing the initial value in a degenerate parabolic equation. Problems of this type have important applications in several fields of applied science. The key to numerically solve such problem is to construct highorder difference schemes for corresponding forward problem. However, the dumping point method which is widely-used for numerically solving classical heat conduction equations cannot be applied to degenerate parabolic equations, because the principal coefficients are zero on degenerate boundaries. In this paper, a new but quite simple technique is proposed to construct a difference scheme of second order accuracy, and the stability and convergence of the scheme are proved. In order to accelerate the convergence rate, the conjugate gradient method is adopted to obtain numerical solutions of the inverse problem. Numerical verification on the efficiency and accuracy of the proposed algorithm is also performed.

Key words: Degenerate parabolic equation, Inverse initial problem, Conjugate gradient method, Convergence, Numerical results

中图分类号:

O212.62

杨柳,邓醉茶. 退化抛物型方程的一个初值反演问题[J]. 数学物理学报, 2020, 40(4): 891-903.

Liu Yang,Zuicha Deng. An Inverse Initial Value Problem for Degenerate Parabolic Equations[J]. Acta mathematica scientia,Series A, 2020, 40(4): 891-903.

参考文献

1	Buchot J M , Raymond J P . The linearized Crocco equation. J Math Fluid Mech, 2006, 8 (4): 510- 541
2	Cannarsa P , Martinez P , Vancostenoble J . Null controllability of degenerate heat equations. Adv Differ Equ, 2005, 10, 153- 190
3	Cheng J , Liu J J . A quasi Tikhonov regularization for a two-dimensional backward heat problem by a fundamental solution. Inverse Problems, 2008, 24, 065012
4	Deng Z C , Qian K , Rao X B , et al. An inverse problem of identifying the source coefficient in a degenerate heat equation. Inverse Probl Sci Eng, 2015, 23, 498- 517
5	Engl H W , Hanke M , Neubauer A . Regularization of Inverse Problems. Dordrecht: Kluwer Academic Publishers, 1996
6	Johansson T , Lesnic D . Determination of a spacewise dependent heat source. J Comput Appl Math, 2007, 209, 66- 80
7	Johansson T , Lesnic D . A variational method for identifying a spacewise-dependent heat source. IMA J Appl Math, 2007, 72, 748- 760
8	Johansson T , Lesnic D . A procedure for determining a spacewise dependent heat source and the initial temperature. Appl Anal, 2008, 87, 265- 276
9	Lu J F , Guan Z . Numerical Solution of Partial Differential Equations. Beijing: Tsinghua University Press, 2004
10	Oleinik O A, Radkevi? E V. Second Order Differential Equations with Non-Negative Characteristic Form. Providence, RI: American Mathematical Society, 1973
11	Rao X B , Wang Y X , Qian K , et al. Numerical simulation for an inverse source problem in a degenerate parabolic equation. Appl Math Modelling, 2015, 39, 7537- 7553
12	Sun Z Z . Numerical Solution of Partial Differential Equations. Beijing: Science Press, 2005
13	Tort J , Vancostenoble J . Determination of the insolation function in the nonlinear Sellers climate model. Ann I H Poincaré-AN, 2012, 29, 683- 713
14	Yang L , Deng Z C , Yu J N , Luo G W . Optimization method for the inverse problem of reconstructing the source term in a parabolic equation. Math Comput Simul, 2009, 80, 314- 326
15	Yang L , Dehghan M , Yu J N , Luo G W . Inverse problem of time-dependent heat sources numerical reconstruction. Math Comput Simul, 2011, 80, 1656- 1672
16	Yang L , Deng Z C . An inverse backward problem for degenerate parabolic equations. Numer Meth Part Differ Equ, 2017, 33, 1900- 1923

[1]	刘鹏杰, 邵虎, 李佳宜, 董端端. 一个基于 BFGS 更新的惯性三项共轭梯度投影法[J]. 数学物理学报, 2025, 45(5): 1698-1710.
[2]	王琼琼, 唐嘉. 基于切比雪夫多项式加速求解 PageRank 的类海森伯格算法[J]. 数学物理学报, 2025, 45(4): 1291-1300.
[3]	于冬梅, 刘大熠. 求解广义互补问题的 Levenberg-Marquardt 算法[J]. 数学物理学报, 2025, 45(4): 1311-1326.
[4]	马昌凤, 谢亚君, 卜凡. 求解Stein张量方程的张量格式BCGSTAB算法[J]. 数学物理学报, 2024, 44(6): 1652-1664.
[5]	王俊杰. 空间分数阶 KGS 方程组的辛差分格式[J]. 数学物理学报, 2024, 44(5): 1319-1334.
[6]	李姣芬, 孔鲁源, 宋佳铄, 文娅琼. 特征提取中一类矩阵迹函数极值问题的黎曼优化算法[J]. 数学物理学报, 2024, 44(4): 1012-1036.
[7]	马小军, 陈富, 贾芝福. 基于非 Lipschitz 步长策略的临近分裂可行问题的强收敛性研究[J]. 数学物理学报, 2024, 44(4): 1052-1065.
[8]	简金宝, 代钰, 尹江华. 分裂可行性问题的一个惯性共轭梯度投影法[J]. 数学物理学报, 2024, 44(4): 1066-1079.
[9]	聂佳琳, 龙宪军. 求解非光滑鞍点问题的黄金比率原始对偶算法[J]. 数学物理学报, 2024, 44(4): 1080-1091.
[10]	张潇, 张宏武. 分数阶椭圆方程反边值问题的分数 Tikhonov 正则化方法[J]. 数学物理学报, 2024, 44(4): 978-993.
[11]	蔡宇, 周光辉. 一种 WYL 型谱共轭梯度法的全局收敛性[J]. 数学物理学报, 2024, 44(1): 173-184.
[12]	简金宝,林惠,马国栋. 大规模非凸不可分优化问题的分裂序列二次规划算法[J]. 数学物理学报, 2023, 43(4): 1284-1296.
[13]	于一康,牛晶. 一类椭圆型界面问题的数值算法[J]. 数学物理学报, 2023, 43(3): 883-895.
[14]	刘鹏杰, 吴彦强, 邵枫, 张艳, 邵虎. 两个带重启方向的改进 HS 型共轭梯度法[J]. 数学物理学报, 2023, 43(2): 570-580.
[15]	谢亚君,马昌凤. 源于自由边值离散的弱非线性互补问题的m+1阶收敛性算法[J]. 数学物理学报, 2022, 42(5): 1506-1516.