摘要
以全局支撑算子方法为基础,通过引入面通量,构造了具有局部模板点的时空二阶精度格式。对于大变形扭曲网格,格式采用法向修正技术和合理的单元角体积计算方法,可以保持通量的精确性。算例表明该格式在非凸网格上能够精确获得线性解;在非光滑网格上可以达到时空二阶精度;能够较好地保持对称性;并适合于三维非结构网格上的求解。
Numerical solution of the diffusion equation plays a key role in the study of inertial confinement fusion(ICF). In this paper, based on the global support operator method, a flux-based scheme is proposed. The scheme has local stencil with sec- ond-order accuracy both in space and time. For strongly distorted meshes, a procedure of normal direction fix is adopted with proper methods for the computation of corner volume weights, which obtains accurate discretization of the face flux. Numerical experiments show that the scheme can obtain accurate solution for linear problems on non-convex meshes. The method has sec- ond-order spatial and temporal accuracy on non-smooth meshes. The method can also preserve the symmetry well and can be ex- tended to the three dimensional unstructured meshes.
出处
《强激光与粒子束》
EI
CAS
CSCD
北大核心
2012年第11期2618-2622,共5页
High Power Laser and Particle Beams
基金
国家自然科学基金项目(11205016
11072040)
中国工程物理研究院科学技术发展基金项目(2011B0201042)
关键词
扩散方程
大变形网格
凹网格
支撑算子格式
非结构网格
diffusion equation
strongly distorted meshes
non-convex meshes
support operator method
unstruc-tured meshes