A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first ...A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.展开更多
In this paper,two classes of Riesz space fractional partial differential equations including space-fractional and space-time-fractional ones are considered.These two models can be regarded as the generalization of the...In this paper,two classes of Riesz space fractional partial differential equations including space-fractional and space-time-fractional ones are considered.These two models can be regarded as the generalization of the classical wave equation in two space dimensions.Combining with the Crank-Nicolson method in temporal direction,efficient alternating direction implicit Galerkin finite element methods for solving these two fractional models are developed,respectively.The corresponding stability and convergence analysis of the numerical methods are discussed.Numerical results are provided to verify the theoretical analysis.展开更多
Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been...Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been discovered that the higher-order accurate method can give reliable and efficient computational results, as well as better resolution of the complex flow fields with multi-scale structures. Compact finite difference schemes, which feature higher-order accuracy and spectral-like resolution with smaller stencils and easier application of boundary conditions, has attracted more and more interest and attention.展开更多
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.展开更多
An alternating direction implicit (ADI) Galerkin method with moving finite element spaces is formulated for a class of second order hyperbolic equations in two space variables. A priori H 1 error estimate is derived.
Based on rectangular partition and bilinear interpolation,we construct an alternating-direction implicit(ADI)finite volume element method,which combined the merits of finite volume element method and alternating direc...Based on rectangular partition and bilinear interpolation,we construct an alternating-direction implicit(ADI)finite volume element method,which combined the merits of finite volume element method and alternating direction implicit method to solve a viscous wave equation with variable coefficients.This paper presents a general procedure to construct the alternating-direction implicit finite volume element method and gives computational schemes.Optimal error estimate in L2 norm is obtained for the schemes.Compared with the finite volume element method of the same convergence order,our method is more effective in terms of running time with the increasing of the computing scale.Numerical experiments are presented to show the efficiency of our method and numerical results are provided to support our theoretical analysis.展开更多
In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is i...In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is introduced recently.In this paper,a theoretical proof of the stabilityof the three-dimensional(3-D)ADI-FDTD method is presented.It is shown that the 3-D ADI-FDTDmethod is unconditionally stable and free from the CFL condition restraint.展开更多
In this paper, an alternating direction Galerkin finite element method is presented for solving 2D time fractional reaction sub-diffusion equation with nonlinear source term. Firstly, one order implicit-explicit metho...In this paper, an alternating direction Galerkin finite element method is presented for solving 2D time fractional reaction sub-diffusion equation with nonlinear source term. Firstly, one order implicit-explicit method is used for time discretization, then Galerkin finite element method is adopted for spatial discretization and obtain a fully discrete linear system. Secondly, Galerkin alternating direction procedure for the system is derived by adding an extra term. Finally, the stability and convergence of the method are analyzed rigorously. Numerical results confirm the accuracy and efficiency of the proposed method.展开更多
文摘A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.
基金supported by the Guangxi Natural Science Foundation[grant numbers 2018GXNSFBA281020,2018GXNSFAA138121]the Doctoral Starting up Foundation of Guilin University of Technology[grant number GLUTQD2016044].
文摘In this paper,two classes of Riesz space fractional partial differential equations including space-fractional and space-time-fractional ones are considered.These two models can be regarded as the generalization of the classical wave equation in two space dimensions.Combining with the Crank-Nicolson method in temporal direction,efficient alternating direction implicit Galerkin finite element methods for solving these two fractional models are developed,respectively.The corresponding stability and convergence analysis of the numerical methods are discussed.Numerical results are provided to verify the theoretical analysis.
文摘Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been discovered that the higher-order accurate method can give reliable and efficient computational results, as well as better resolution of the complex flow fields with multi-scale structures. Compact finite difference schemes, which feature higher-order accuracy and spectral-like resolution with smaller stencils and easier application of boundary conditions, has attracted more and more interest and attention.
基金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.
基金the National Natural Sciences Foundation of China
文摘An alternating direction implicit (ADI) Galerkin method with moving finite element spaces is formulated for a class of second order hyperbolic equations in two space variables. A priori H 1 error estimate is derived.
基金supported by the National Natural Science Foundation of China grants No.11971241.
文摘Based on rectangular partition and bilinear interpolation,we construct an alternating-direction implicit(ADI)finite volume element method,which combined the merits of finite volume element method and alternating direction implicit method to solve a viscous wave equation with variable coefficients.This paper presents a general procedure to construct the alternating-direction implicit finite volume element method and gives computational schemes.Optimal error estimate in L2 norm is obtained for the schemes.Compared with the finite volume element method of the same convergence order,our method is more effective in terms of running time with the increasing of the computing scale.Numerical experiments are presented to show the efficiency of our method and numerical results are provided to support our theoretical analysis.
基金Supported by the Specialized Research Fund for the Doctoral Program of Higher Education(No.20010614003)
文摘In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is introduced recently.In this paper,a theoretical proof of the stabilityof the three-dimensional(3-D)ADI-FDTD method is presented.It is shown that the 3-D ADI-FDTDmethod is unconditionally stable and free from the CFL condition restraint.
文摘In this paper, an alternating direction Galerkin finite element method is presented for solving 2D time fractional reaction sub-diffusion equation with nonlinear source term. Firstly, one order implicit-explicit method is used for time discretization, then Galerkin finite element method is adopted for spatial discretization and obtain a fully discrete linear system. Secondly, Galerkin alternating direction procedure for the system is derived by adding an extra term. Finally, the stability and convergence of the method are analyzed rigorously. Numerical results confirm the accuracy and efficiency of the proposed method.
