A moving mesh finite difference method for equilibrium radiation diffusion equations
- Department of Mathematics, College of Science, China University of Mining and Technology, Xuzhou, Jiangsu 221116 (China)
- Department of Mathematics, University of Kansas, Lawrence, KS 66045 (United States)
- School of Mathematical Sciences and Fujian Provincial Key Laboratory of Mathematical Modeling and High-Performance Scientific Computing, Xiamen University, Xiamen, Fujian 361005 (China)
An efficient moving mesh finite difference method is developed for the numerical solution of equilibrium radiation diffusion equations in two dimensions. The method is based on the moving mesh partial differential equation approach and moves the mesh continuously in time using a system of meshing partial differential equations. The mesh adaptation is controlled through a Hessian-based monitor function and the so-called equidistribution and alignment principles. Several challenging issues in the numerical solution are addressed. Particularly, the radiation diffusion coefficient depends on the energy density highly nonlinearly. This nonlinearity is treated using a predictor–corrector and lagged diffusion strategy. Moreover, the nonnegativity of the energy density is maintained using a cutoff method which has been known in literature to retain the accuracy and convergence order of finite difference approximation for parabolic equations. Numerical examples with multi-material, multiple spot concentration situations are presented. Numerical results show that the method works well for radiation diffusion equations and can produce numerical solutions of good accuracy. It is also shown that a two-level mesh movement strategy can significantly improve the efficiency of the computation.
- OSTI ID:
- 22465664
- Journal Information:
- Journal of Computational Physics, Vol. 298; Other Information: Copyright (c) 2015 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
A massively-parallel, unstructured overset method to simulate moving bodies in turbulent flows
Parallel heterogeneous mesh refinement for multidimensional convection-diffusion equations using an Euler-Lagrange method