On conforming finite element methods for the inhomogeneous stationary Navier-Stokes equations

Summary We consider the stationary Navier-Stokes equations, written in terms of the primitive variables, in the case where both the partial differential equations and boundary conditions are inhomogeneous. Under certain conditions on the data, the existence and uniqueness of the solution of a weak formulation of the equations can be guaranteed. A conforming finite element method is presented and optimal estimates for the error of the approximate solution are proved. In addition, the convergence properties of iterative methods for the solution of the discrete nonlinear algebraic systems resulting from the finite element algorithm are given. Numerical examples, using an efficient choice of finite element spaces, are also provided.Supported, in part, by the U.S. Air Force Office of Scientific Research under Grant No. AF-AFOSR-80-0083Supported, in part, by the same agency under Grant No. AF-AFOSR-80-0176-A. Both authors were also partially supported by NASA Contract No. NAS1-15810 while they were in residence at the Institute for Computer Applications in Science and Engineering, NASA Langley Research Center, Hampton, VA 23665, USA     

Minimal dissipation hybrid bicompact schemes for hyperbolic equations

New monotonicity-preserving hybrid schemes are proposed for multidimensional hyperbolic equations. They are convex combinations of high-order accurate central bicompact schemes and upwind schemes of first-order accuracy in time and space. The weighting coefficients in these combinations depend on the local difference between the solutions produced by the high- and low-order accurate schemes at the current space-time point. The bicompact schemes are third-order accurate in time, while having the fourth order of accuracy and the first difference order in space. At every time level, they can be solved by marching in each spatial variable without using spatial splitting. The upwind schemes have minimal dissipation among all monotone schemes constructed on a minimum space-time stencil. The hybrid schemes constructed has been successfully tested as applied to a number of two-dimensional gas dynamics benchmark problems.     

Transformation of three-dimensional regions onto rectangular regions by elliptic systems

Summary A transformation method is developed which may be used to solve various types of boundary value problems on three-dimensional regions with an arbitrary boundary. The implementation of the method is illustrated in the solution of a potential flow problem. All computations are performed on a cubic mesh in a rectangular region.This report was prepared as a result of work performed under NASA Contract No. NAS1-14101 while the first author was in residence at ICASE, NASA Langley Research Center, Hampton, VA 23665, USA     

基于不光滑边界的变系数抛物型方程的高精度紧格式

本文讨论基于不光滑边界的变系数抛物型方程求解的高精度紧格式.首先构造一般变系数抛物型方程的高精度紧格式,并在理论上证明格式具有空间方向四阶精度.然后针对非光滑边界条件,引入局部网格加密技巧在奇异点附近进行不均匀的网格加密.数值实验以期权定价中Black-Scholes偏微分方程的求解为例,验证高精度紧格式用于光滑边界条件的微分方程离散可以达到四阶精度.对于处理非光滑边界条件,网格局部加密技巧能有效的提高数值解精度,使得高精度紧格式用于定价欧式期权可以接近四阶精度.     

High‐order difference schemes for 2D elliptic and parabolic problems with mixed derivatives

We propose a 9‐point fourth‐order finite difference scheme for 2D elliptic problems with a mixed derivative and variable coefficients. The same approach is extended to derive a class of two‐level high‐order compact schemes with weighted time discretization for solving 2D parabolic problems with a mixed derivative. The schemes are fourth‐order accurate in space and second‐ or lower‐order accurate in time depending on the choice of a weighted average parameter μ. Unconditional stability is proved for 0.5 ≤ μ ≤ 1, and numerical experiments supporting our theoretical analysis and confirming the high‐order accuracy of the schemes are presented. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 366–378, 2007     

EFFICIENT LINEAR SCHEMES WITH UNCONDITIONAL ENERGY STABILITY FOR THE PHASE FIELD MODEL OF SOLID-STATE DEWETTING PROBLEMS

In this paper, we study linearly first and second order in time, uniquely solvable and unconditionally energy stable numerical schemes to approximate the phase field model of solid-state dewetting problems based on the novel "scalar auxiliary variable" (SAV) approach, a new developed efficient and accurate method for a large class of gradient flows. The schemes are based on the first order Euler method and the second order backward differential formulas (BDF2) for time discretization, and finite element methods for space discretization. The proposed schemes are proved to be unconditionally stable and the discrete equations are uniquely solvable for all time steps. Various numerical experiments are presented to validate the stability and accuracy of the proposed schemes.     

A high‐order finite difference discretization strategy based on extrapolation for convection diffusion equations

We propose a new high‐order finite difference discretization strategy, which is based on the Richardson extrapolation technique and an operator interpolation scheme, to solve convection diffusion equations. For a particular implementation, we solve a fine grid equation and a coarse grid equation by using a fourth‐order compact difference scheme. Then we combine the two approximate solutions and use the Richardson extrapolation to compute a sixth‐order accuracy coarse grid solution. A sixth‐order accuracy fine grid solution is obtained by interpolating the sixth‐order coarse grid solution using an operator interpolation scheme. Numerical results are presented to demonstrate the accuracy and efficacy of the proposed finite difference discretization strategy, compared to the sixth‐order combined compact difference (CCD) scheme, and the standard fourth‐order compact difference (FOC) scheme. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 18–32, 2004.     

Maximum norm error bounds of ADI and compact ADI methods for solving parabolic equations

Alternating direction implicit (ADI) schemes are computationally efficient and widely utilized for numerical approximation of the multidimensional parabolic equations. By using the discrete energy method, it is shown that the ADI solution is unconditionally convergent with the convergence order of two in the maximum norm. Considering an asymptotic expansion of the difference solution, we obtain a fourth‐order, in both time and space, approximation by one Richardson extrapolation. Extension of our technique to the higher‐order compact ADI schemes also yields the maximum norm error estimate of the discrete solution. And by one extrapolation, we obtain a sixth order accurate approximation when the time step is proportional to the squares of the spatial size. An numerical example is presented to support our theoretical results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

On higher-order semi-explicit symplectic partitioned Runge-Kutta methods for constrained Hamiltonian systems

Summary. In this paper we generalize the class of explicit partitioned Runge-Kutta (PRK) methods for separable Hamiltonian systems to systems with holonomic constraints. For a convenient analysis of such schemes, we first generalize the backward error analysis for systems in to systems on manifolds embedded in . By applying this analysis to constrained PRK methods, we prove that such methods will, in general, suffer from order reduction as well-known for higher-index differential-algebraic equations. However, this order reduction can be avoided by a proper modification of the standard PRK methods. This modification increases the number of projection steps onto the constraint manifold but leaves the number of force evaluations constant. We also give a numerical comparison of several second, fourth, and sixth order methods. Received May 5, 1995 / Revised version received February 7, 1996     

JACOBI PSEUDOSPECTRAL METHOD FOR FOURTH ORDER PROBLEMS

In this paper, we investigate Jacobi pseudospectral method for fourth order problems. We establish some basic results on the Jacobi-Gauss-type interpolations in non-uniformly weighted Sobolev spaces, which serve as important tools in analysis of numerical quadratures, and numerical methods of differential and integral equations. Then we propose Jacobi pseudospectral schemes for several singular problems and multiple-dimensional problems of fourth order. Numerical results demonstrate the spectral accuracy of these schemes, and coincide well with theoretical analysis.     

High‐order energy‐preserving schemes for the improved Boussinesq equation

This article proposes a class of high‐order energy‐preserving schemes for the improved Boussinesq equation. To derive the energy‐preserving schemes, we first discretize the improved Boussinesq equation by Fourier pseudospectral method, which leads to a finite‐dimensional Hamiltonian system. Then, the obtained semidiscrete system is solved by Hamiltonian boundary value methods, which is a newly developed class of energy‐preserving methods. The proposed schemes can reach spectral precision in space, and in time can reach second‐order, fourth‐order, and sixth‐order accuracy, respectively. Moreover, the proposed schemes can conserve the discrete mass and energy to within machine precision. Furthermore, to show the efficiency and accuracy of the proposed methods, the proposed methods are compared with the finite difference methods and the finite volume element method. The results of several numerical experiments are given for the propagation of the single solitary wave, the interaction of two solitary waves and the wave break‐up.     

Efficient orthogonal spline collocation methods for solving linear second order hyperbolic problems on rectangles

Summary. Piecewise Hermite bicubic orthogonal spline collocation Laplace-modified and alternating-direction schemes for the approximate solution of linear second order hyperbolic problems on rectangles are analyzed. The schemes are shown to be unconditionally stable and of optimal order accuracy in the and discrete maximum norms for space and time, respectively. Implementations of the schemes are discussed and numerical results presented which demonstrate the accuracy and rate of convergence using various norms. Received November 7, 1994 / Revised version received April 29, 1996     

A high‐order compact boundary value method for solving one‐dimensional heat equations

We combine fourth‐order boundary value methods (BVMs) for discretizing the temporal variable with fourth‐order compact difference scheme for discretizing the spatial variable to solve one‐dimensional heat equations. This class of new compact difference schemes achieve fourth‐order accuracy in both temporal and spatial variables and are unconditionally stable due to the favorable stability property of BVMs. Numerical results are presented to demonstrate the accuracy and efficiency of the new compact difference scheme, compared to the standard second‐order Crank‐Nicolson scheme. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 846–857, 2003.     

Multiderivative Runge-Kutta processes for two-point boundary value problems

A class of finite difference schemes for the solution of a nonlinear system of first order differential equations with two point boundary conditions which shares properties with Runge-Kutta processes and gap schemes is discussed. The order conditions for the coefficients of these processes, techniques for reducing these order conditions in number and the symmetry conditions are given. A symmetricA-stable eight order process which has second, fourth and sixth orderA-stable processes embedded in it is given as an example.Research supported in part by the United States Air Force under contract AFOSR-89-0383.     

Asynchronous and corrected-asynchronous finite difference solutions of PDEs on MIMD multiprocessors

A major problem in achieving significant speed-up on parallel machines is the overhead involved with synchronizing the concurrent processes. Removing the synchronization constraint has the potential of speeding up the computation, while maintaining greater computation flexibility (e.g. differences in processors speed; differences in the data input to processors). We construct asynchronous (AS) finite difference schemes for the solution of PDEs by removing the synchronization constraint. We analyze the numerical properties of these schemes. Based on the analysis, we develop corrected-asynchronous (CA) finite difference schemes which are specifically constructed for an asynchronous processing. We present asynchronous (AS) and corrected-asynchronous (CA) finite difference schemes for the multi-dimensional heat equation. Although our discussion concentrates on the Euler scheme it should serve only as a sample, as it can be extended to other schemes and other PDEs.These schemes are implemented on the shared-memory multi-userSequent Balance machine. Numerical results for one and two dimensional problems are presented. It is shown experimentally that synchronization penalty can be about 50% of run time: in most cases, the asynchronous scheme runs twice as fast as the parallel synchronous scheme. In general, the efficiency of the parallel schemes increases with processor load, with the time-level, and with the problem dimension. The efficiency of the AS may reach 90% and over, but it provides accurate results only for steady-state values. The CA, on the other hand, is less efficient but provides more accurate results for intermediate (non steady-state) values. The results show the potential of developing asynchronous finite deference schemes for steady-state as well as non steadystate problems.This research was partially supported by a grant from The Basic Research Foundation administrated by The Israel Academy of Sciences and Humanities.A reduced version of the paper was presented at the 4th SIAM Conference on Parallel Processing for Scientific Computing, Dec. 11–13, 1989, Chicago, USA.The work by this author was supported by research grant 337 of the Israeli National Council for Research and Development in the years 1990–1991.This research was supported by the National Aeronautics and Space Administration under NASA Contract No. NASI-18107 while the author was in residence at the Institute for Computer Applications in Sciences and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665, USA.     

High accuracy A.D.I. methods for fourth order parabolic equations with variable coefficients

High accuracy alternating direction implicit (A.D.I.) methods are derived for solving fourth order parabolic equations with variable coefficients in one, two, and three space dimensions. Splittings are discussed and numerical results are presented.     

Explicit-implicit difference scheme for the joint solution of the radiative transfer and energy equations by the splitting method

High-order accurate explicit and implicit conservative predictor-corrector schemes are presented for the radiative transfer and energy equations in the multigroup kinetic approximation solved together by applying the splitting method with respect to physical processes and spatial variables. The original system of integrodifferential equations is split into two subsystems: one of partial differential equations without sources and one of ordinary differential equations (ODE) with sources. The general solution of the ODE system and the energy equation is written in quadratures based on total energy conservation in a cell. A feature of the schemes is that a new approximation is used for the numerical fluxes through the cell interfaces. The fluxes are found along characteristics with the interaction between radiation and matter taken into account. For smooth solutions, the schemes approximating the transfer equations on spatially uniform grids are second-order accurate in time and space. As an example, numerical results for Fleck’s test problems are presented that confirm the increased accuracy and efficiency of the method.     

High-order accurate monotone compact running scheme for multidimensional hyperbolic equations

Monotone absolutely stable conservative difference schemes intended for solving quasilinear multidimensional hyperbolic equations are described. For sufficiently smooth solutions, the schemes are fourth-order accurate in each spatial direction and can be used in a wide range of local Courant numbers. The order of accuracy in time varies from the third for the smooth parts of the solution to the first near discontinuities. This is achieved by choosing special weighting coefficients that depend locally on the solution. The presented schemes are numerically efficient thanks to the simple two-diagonal (or block two-diagonal) structure of the matrix to be inverted. First the schemes are applied to system of nonlinear multidimensional conservation laws. The choice of optimal weighting coefficients for the schemes of variable order of accuracy in time and flux splitting is discussed in detail. The capabilities of the schemes are demonstrated by computing well-known two-dimensional Riemann problems for gasdynamic equations with a complex shock wave structure.     

Flux-splitting schemes for parabolic problems

Splitting with respect to space variables can be used in solving boundary value problems for second-order parabolic equations. Classical alternating direction methods and locally one-dimensional schemes could be examples of this approach. For problems with rapidly varying coefficients, a convenient tool is the use of fluxes (directional derivatives) as independent variables. The original equation is written as a system in which not only the desired solution but also directional derivatives (fluxes) are unknowns. In this paper, locally one-dimensional additional schemes (splitting schemes) for second-order parabolic equations are examined. By writing the original equation in flux variables, certain two-level locally one-dimensional schemes are derived. The unconditional stability of locally one-dimensional flux schemes of the first and second approximation order with respect to time is proved.