A finite volume element method is developed for analyzing unsteady scalar reaction-diffusion problems in two dimensions. The method combines the concepts that are employed in the finite volume and the finite element m...A finite volume element method is developed for analyzing unsteady scalar reaction-diffusion problems in two dimensions. The method combines the concepts that are employed in the finite volume and the finite element method together. The finite volume method is used to discretize the unsteady reaction-diffusion equation, while the finite element method is applied to estimate the gradient quantities at cell faces. Robustness and efficiency of the combined method have been evaluated on uniform rectangular grids by using available numerical solutions of the two-dimensional reaction-diffusion problems. The numerical solutions demonstrate that the combined method is stable and can provide accurate solution without spurious oscillation along the high-gradient boundary layers.展开更多
A new numerical algorithm for telegraph equations with homogeneous boundary con- ditions is proposed. Due to the damping terms in telegraph equations, there is no royal conservation law according to Noether's theorem...A new numerical algorithm for telegraph equations with homogeneous boundary con- ditions is proposed. Due to the damping terms in telegraph equations, there is no royal conservation law according to Noether's theorem. The algorithm origins from the discovery of a transform applied to a telegraph equation, which transforms the telegraph equation to a Klein-Gordon equation. The Symplectic method is then brought in this algorithm to solve the Klein-Gordon equation, which is based on the fact that the Klein-Gordon equation with the homogeneous boundary condition is a perfect Hamiltonian system and the symplectic method works very well for Hamiltonian systems. The transformation itself and the inverse transformation theoretically bring no error to the numerical computation. Therefore the error only comes from the symplectic scheme chosen. The telegraph equation is finally explicitly computed when an explicit symplectic scheme is utilized. A relatively long time result can be expected due to the application of the symplectic method. Mean- while, we present order analysis for both one-dimensional and multi-dimensional cases in the paper. The efficiency of this approach is demonstrated with numerical examples.展开更多
文摘A finite volume element method is developed for analyzing unsteady scalar reaction-diffusion problems in two dimensions. The method combines the concepts that are employed in the finite volume and the finite element method together. The finite volume method is used to discretize the unsteady reaction-diffusion equation, while the finite element method is applied to estimate the gradient quantities at cell faces. Robustness and efficiency of the combined method have been evaluated on uniform rectangular grids by using available numerical solutions of the two-dimensional reaction-diffusion problems. The numerical solutions demonstrate that the combined method is stable and can provide accurate solution without spurious oscillation along the high-gradient boundary layers.
基金Acknowledgments. This work is supported by the Natural Science Foundation of China (NSFC) under grant 11371287 and the International Science and Technology Cooperation Program of China under grant 2010DFA14700.
文摘A new numerical algorithm for telegraph equations with homogeneous boundary con- ditions is proposed. Due to the damping terms in telegraph equations, there is no royal conservation law according to Noether's theorem. The algorithm origins from the discovery of a transform applied to a telegraph equation, which transforms the telegraph equation to a Klein-Gordon equation. The Symplectic method is then brought in this algorithm to solve the Klein-Gordon equation, which is based on the fact that the Klein-Gordon equation with the homogeneous boundary condition is a perfect Hamiltonian system and the symplectic method works very well for Hamiltonian systems. The transformation itself and the inverse transformation theoretically bring no error to the numerical computation. Therefore the error only comes from the symplectic scheme chosen. The telegraph equation is finally explicitly computed when an explicit symplectic scheme is utilized. A relatively long time result can be expected due to the application of the symplectic method. Mean- while, we present order analysis for both one-dimensional and multi-dimensional cases in the paper. The efficiency of this approach is demonstrated with numerical examples.