走向现代数学学术报告 - 高志明研究员（No. 658）

题目：The positivity-preserving finite volume scheme with fixed stencils for radiation diffusion problems on general polyhedral meshes

报告人：高志明 研究员（北京应用物理与计算数学研究所）

时间：2023年11月10日，19:30

腾讯会议ID：774106541

摘要：Diffusion process is very important in many coupled systems for real life applications such as radiation hydrodynamics, plasma physics, reservoir modelling and so on. Finite volume discretizations of such systems are very popular among the scientists and engineers due to their nice properties such as local conservation. In many situations, one of the significant features of the finite volume scheme is that the discrete solution should be non-negative. The severely distorted meshes (often occur in the simulation of inertial confinement fusion) and highly anisotropic media (in reservoir modelling) may result in the violation of positivity-preserving property. Many scientists have made great efforts on the construction of positivity-preserving finite volume schemes of diffusion equation that approximately have a second-order accuracy on severely distorted grids in case that the diffusion tensor is taken to be highly anisotropic, at times heterogeneous, and/or discontinuous.

In this talk, two kinds of nonlinear cell-centered positivity-preserving finite volume schemes are proposed for radiation diffusion problems on general three-dimensional polyhedral meshes. Firstly, the one-sided flux on the cell-faces is discretized using the fixed stencil of all vertices. Then by using the nonlinear two-point flux approximation, the cell-centered discretization scheme is obtained by eliminating the vertex auxiliary unknowns. On this basis, a new explicit weighted second-order vertex interpolation algorithm for arbitrary polyhedral meshes is designed to eliminate the vertex auxiliary unknowns in the scheme. In addition, an improved Anderson acceleration algorithm is adopted for the nonlinear iteration. Finally, some benchmark examples are given to verify the convergence and positivity-preserving property of the two discretization schemes.

报告人个人简介：高志明，北京应用物理与计算数学研究所研究员，博士生导师，中国工程物理研究院第八届科学技术委员会计算数学专业组委员。2008年毕业于西安交通大学计算数学系，博士论文获全国优秀博士论文提名。长期从事极端条件下辐射扩散方程数值方法研究、多物理耦合数值模拟软件研制及应用工作。在 SIAM J. Sci. Comput.，J.Comput.Phys., Comput Methods Appl. Mech. Engrg.等杂志发表论文60余篇，主持国家自然科学基金、JKW创新基金、中物院院长基金多项。