A Wavelet-Galerkin method is proposed to solve the singular perturbation problem with boundary layers numerically. Because there are boundary layers in the solution of the singular perturbation problem, the approximat...A Wavelet-Galerkin method is proposed to solve the singular perturbation problem with boundary layers numerically. Because there are boundary layers in the solution of the singular perturbation problem, the approximation spaces with different scale wavelets and boundary bases are chosen. In addition, the computation of the inner integrals is transformed to an eigenvalue problem. Therefore, a high accuracy method with reasonable computation is obtained. On the other hand, there is an explicit diagonal preconditioning which makes the condition number of the stiff matrix become bounded by a constant. The error estimate of the Wavelet-Galerkin method and the analysis of the computation complexity are given. The numerical examples show that the method is feasible and effective for solving the singular perturbation problem with boundary layers numerically.展开更多
In this paper, we use Daubechies scaling functions as test functions for the Galerkin method, and discuss Wavelet-Galerkin solutions for the Hamilton-Jacobi equations. It can be proved that the schemes are TVD schemes...In this paper, we use Daubechies scaling functions as test functions for the Galerkin method, and discuss Wavelet-Galerkin solutions for the Hamilton-Jacobi equations. It can be proved that the schemes are TVD schemes. Numerical tests indicate that the schemes are suitable for the Hamilton-Jacobi equations. Furthermore, they have high-order accuracy in smooth regions and good resolution of singularities.展开更多
基金Doctoral Program Foundation ofHigher Education of China
文摘A Wavelet-Galerkin method is proposed to solve the singular perturbation problem with boundary layers numerically. Because there are boundary layers in the solution of the singular perturbation problem, the approximation spaces with different scale wavelets and boundary bases are chosen. In addition, the computation of the inner integrals is transformed to an eigenvalue problem. Therefore, a high accuracy method with reasonable computation is obtained. On the other hand, there is an explicit diagonal preconditioning which makes the condition number of the stiff matrix become bounded by a constant. The error estimate of the Wavelet-Galerkin method and the analysis of the computation complexity are given. The numerical examples show that the method is feasible and effective for solving the singular perturbation problem with boundary layers numerically.
基金the National Natural Science Foundation of China(No.10571178)
文摘In this paper, we use Daubechies scaling functions as test functions for the Galerkin method, and discuss Wavelet-Galerkin solutions for the Hamilton-Jacobi equations. It can be proved that the schemes are TVD schemes. Numerical tests indicate that the schemes are suitable for the Hamilton-Jacobi equations. Furthermore, they have high-order accuracy in smooth regions and good resolution of singularities.
