Fourth-order alternating direction implicit compact finite difference schemes for two-dimensional Schrödinger equations

2011 ◽  
Vol 61 (4) ◽  
pp. 593-614 ◽  
Author(s):  
Zhen Gao ◽  
Shusen Xie
2001 ◽  
Vol 7 (3) ◽  
pp. 283-297 ◽  
Author(s):  
Mehdi Dehghan

Two different finite difference schemes for solving the two-dimensional parabolic inverse problem with temperature overspecification are considered. These schemes are developed for indentifying the control parameter which produces, at any given time, a desired temperature distribution at a given point in the spatial domain. The numerical methods discussed, are based on the (3,3) alternating direction implicit (ADI) finite difference scheme and the (3,9) alternating direction implicit formula. These schemes are unconditionally stable. The basis of analysis of the finite difference equation considered here is the modified equivalent partial differential equation approach, developed from the 1974 work of Warming and Hyett [17]. This allows direct and simple comparison of the errors associated with the equations as well as providing a means to develop more accurate finite difference schemes. These schemes use less central processor times than the fully implicit schemes for two-dimensional diffusion with temperature overspecification. The alternating direction implicit schemes developed in this report use more CPU times than the fully explicit finite difference schemes, but their unconditional stability is significant. The results of numerical experiments are presented, and accuracy and the Central Processor (CPU) times needed for each of the methods are discussed. We also give error estimates in the maximum norm for each of these methods.


Mathematics ◽  
2021 ◽  
Vol 9 (19) ◽  
pp. 2508
Author(s):  
Jesús Amo-Navarro ◽  
Ricardo Vinuesa ◽  
J. Alberto Conejero ◽  
Sergio Hoyas

In fluid mechanics, the bi-Laplacian operator with Neumann homogeneous boundary conditions emerges when transforming the Navier–Stokes equations to the vorticity–velocity formulation. In the case of problems with a periodic direction, the problem can be transformed into multiple, independent, two-dimensional fourth-order elliptic problems. An efficient method to solve these two-dimensional bi-Laplacian operators with Neumann homogeneus boundary conditions was designed and validated using 2D compact finite difference schemes. The solution is formulated as a linear combination of auxiliary solutions, as many as the number of points on the boundary, a method that was prohibitive some years ago due to the large memory requirements to store all these auxiliary functions. The validation has been made for different field configurations, grid sizes, and stencils of the numerical scheme, showing its potential to tackle high gradient fields as those that can be found in turbulent flows.


2010 ◽  
Vol 2010 ◽  
pp. 1-17 ◽  
Author(s):  
Yaw Kyei ◽  
John Paul Roop ◽  
Guoqing Tang

We derive a family of sixth-order compact finite-difference schemes for the three-dimensional Poisson's equation. As opposed to other research regarding higher-order compact difference schemes, our approach includes consideration of the discretization of the source function on a compact finite-difference stencil. The schemes derived approximate the solution to Poisson's equation on a compact stencil, and thus the schemes can be easily implemented and resulting linear systems are solved in a high-performance computing environment. The resulting discretization is a one-parameter family of finite-difference schemes which may be further optimized for accuracy and stability. Computational experiments are implemented which illustrate the theoretically demonstrated truncation errors.


Sign in / Sign up

Export Citation Format

Share Document