In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids....In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids. Firstly, the original equation is transformed from the physical domain (with a nonuniform mesh) to the computational domain (with a uniform mesh) by using a coordinate transformation. Then, a fourth order compact difference scheme is proposed to solve the transformed elliptic equation on uniform girds. After that, a multigrid method is employed to solve the linear algebraic system arising from the difference equation. At last, the numerical experiments on some elliptic problems with interior/boundary layers are conducted to show high accuracy and high efficiency of the present method.展开更多
In this paper, a family of high-order compact finite difference methods in combination preconditioned methods are used for solution of the Diffusion-Convection equation. We developed numerical methods by replacing the...In this paper, a family of high-order compact finite difference methods in combination preconditioned methods are used for solution of the Diffusion-Convection equation. We developed numerical methods by replacing the time and space derivatives by compact finite-difference approximations. The system of resulting nonlinear finite difference equations are solved by preconditioned Krylov subspace methods. Numerical results are given to verify the behavior of high-order compact approximations in combination preconditioned methods for stability, convergence. Also, the accuracy and efficiency of the proposed scheme are considered.展开更多
An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in ...An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact （HOC） difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit （ADI） scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.展开更多
In this paper, a novel multisymplectic scheme is proposed for the coupled nonlinear Schrodinger-KdV （CNLS-KdV） equations. The CNLS-KdV equations are rewritten into the multisymplectic Hamiltonian form by introducing...In this paper, a novel multisymplectic scheme is proposed for the coupled nonlinear Schrodinger-KdV （CNLS-KdV） equations. The CNLS-KdV equations are rewritten into the multisymplectic Hamiltonian form by introducing some canonical momenta. To simulate the problem efficiently, the CNLS-KdV equations are approximated by a high order compact method in space which preserves N semi-discrete multisymplectic conservation laws. We then discretize the semi-discrete system by using a symplectic midpoint scheme in time. Thus, a full-discrete multisymplectic scheme is obtained for the CNLS-KdV equations. The conservation laws of the full-discrete scheme are analyzed. Some numerical experiments are presented to further verify the convergence and conservation laws of the new scheme.展开更多
文摘In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids. Firstly, the original equation is transformed from the physical domain (with a nonuniform mesh) to the computational domain (with a uniform mesh) by using a coordinate transformation. Then, a fourth order compact difference scheme is proposed to solve the transformed elliptic equation on uniform girds. After that, a multigrid method is employed to solve the linear algebraic system arising from the difference equation. At last, the numerical experiments on some elliptic problems with interior/boundary layers are conducted to show high accuracy and high efficiency of the present method.
文摘In this paper, a family of high-order compact finite difference methods in combination preconditioned methods are used for solution of the Diffusion-Convection equation. We developed numerical methods by replacing the time and space derivatives by compact finite-difference approximations. The system of resulting nonlinear finite difference equations are solved by preconditioned Krylov subspace methods. Numerical results are given to verify the behavior of high-order compact approximations in combination preconditioned methods for stability, convergence. Also, the accuracy and efficiency of the proposed scheme are considered.
基金supported by the National Natural Science Foundation of China(No.51174236)the National Basic Research Program of China(973 Program)(No.2011CB606306)the Opening Project of State Key Laboratory of Porous Metal Materials(No.PMM-SKL-4-2012)
文摘An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact （HOC） difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit （ADI） scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.
基金This work is supported by the NNSFC （Nos. 11771213, 41504078, 11301234, 11271171）, the National Key Research and Development Project of China （No. 2016YFC0600310）, the Major Projects of Natural Sciences of University in Jiangsu Province of China （No. 15KJA110002） and the Priority Academic Program Development of Jiangsu Higher Education Institutions, the Provincial Natural Science Foundation of Jiangxi （Nos. 20161ACB20006, 20142BCB23009, 20151BAB 201012）.
文摘In this paper, a novel multisymplectic scheme is proposed for the coupled nonlinear Schrodinger-KdV （CNLS-KdV） equations. The CNLS-KdV equations are rewritten into the multisymplectic Hamiltonian form by introducing some canonical momenta. To simulate the problem efficiently, the CNLS-KdV equations are approximated by a high order compact method in space which preserves N semi-discrete multisymplectic conservation laws. We then discretize the semi-discrete system by using a symplectic midpoint scheme in time. Thus, a full-discrete multisymplectic scheme is obtained for the CNLS-KdV equations. The conservation laws of the full-discrete scheme are analyzed. Some numerical experiments are presented to further verify the convergence and conservation laws of the new scheme.
