COMPACT DIFFERENCE SCHEMES FOR THE DIFFUSION AND SCHRODINGER EQUATIONS. APPROXIMATION,STABILITY, CONVERGENCE,EFFECTIVENESS, MONOTONY

Various compact difference schemes （both old and new, explicit and implicit, one-level and two-level）, which approximate the diffusion equation and SchrSdinger equation with periodical boundary conditions are constructed by means of the general approach. The results of numerical experiments for various initial data and right hand side are presented. We evaluate the real order of their convergence, as well as their stability, effectiveness, and various kinds of monotony. The optimal Courant number depends on the number of grid knots and on the smoothness of solutions. The competition of various schemes should be organized for the fixed number of arithmetic operations, which are necessary for numerical integration of a given Cauchy problem. This approach to the construction of compact schemes can be developed for numerical solution of various problems of mathematical physics.     

A compact 9 point stencil based on integrated RBFs for the convection–diffusion equation

In this paper, a high-order compact stencil for solving the convection–diffusion equation in two dimensions is proposed. The convection and diffusion terms are both approximated by means of radial basis functions (RBFs) that are constructed over 3×3$3\times 3$ rectangular stencils. Salient features here are that (i) integration is employed to construct local RBF approximations; and (ii) through the constants of integration, values of the convection–diffusion equation at several selected nodes on the stencil are also enforced. Numerical results indicate that (i) the inclusion of the governing equation into the stencil leads to a significant improvement in accuracy; (ii) when the convection dominates, accurate solutions are obtained at a regime of the RBF width which makes the RBFs peaked; and (iii) high levels of accuracy are achieved using relatively coarse grids.     

非定常不可压缩Navier-Stokes方程的旋度形式压力校正格式分析

将求解不可压缩流动的旋度形式压力校正格式从Stokes方程延拓到非定常不可压缩Navier-Stokes方程.在第一步需要求解一个线性的对流-扩散方程,在第二步求解一个Stokes问题.首先给出格式的稳定性分析,然后采用P_2—P_1元分别使用标准形式的压力校正格式和旋度形式的压力校正格式进行了数值模拟,模拟结果表明,对于速度的L~2,H~1误差,两种格式几乎一样,对于压力的L~2误差,旋度形式的压力校正格式略有改进.     

Numerical investigation of hypersonic step-flows

Hypersonic aerospace vehicles are exposed to extreme flight conditions with heavy contour loads during their mission. Especially at ridges and sharp corners, the wall heat flux and pressure may cause serious damage to the body. Sometimes, the surface material cannot resist the high loading and fails completely. In this work the laminar hypersonic flow over forward and backward facing steps is investigated by CFD techniques and the results are compared with experimental data. The selected flow conditions correspond to cold hypersonic flow according to the availability of experimental data. The Navier-Stokes equations in the high temperature gas approximation of a thermally perfect gas in local equilibrium serve as the model for the physical problem. A multiblock finite-volume method is used to discretize consistently all spatial derivatives appearing in the balance equations. A second order in space Godunov-type method is utilized for the non-diffusive part of the governing equations whereas centered differences are used for the diffusive part. Time integration is performed by a second order implicit scheme. In each time step, the resulting nonlinear system of equations is solved by Newton's method employing a relaxation scheme based on conjugate gradients for the linear equation system. The results obtained permit a close insight into the physics of the flow problems under consideration and by this provide valuable information for construction concepts of hypersonic vehicles. Besides a careful comparison of the numerical results with experimental data, numerical aspects like the grid influence are addressed. Received 9 November 1998 / Accepted 2 December 1999     

Existence and stability of periodic solutions of high-order Hopfield neural networks with impulses and delays

By using the continuation theorem of coincidence degree theory and constructing suitable Lyapunov functions, the global exponential stability and periodicity are investigated for a class of delayed high-order Hopfield neural networks (HHNNs) with impulses, which are new and complement previously known results. Finally, an example with numerical simulation is given to show the effectiveness of the proposed method and results. The numerical simulation shows that our models can occur in many forms of complexities including periodic oscillation and the Gui chaotic strange attractor.     

High‐order finite difference schemes for incompressible flows

This paper presents a new high‐order approach to the numerical solution of the incompressible Stokes and Navier–Stokes equations. The class of schemes developed is based upon a velocity–pressure–pressure gradient formulation, which allows: (i) high‐order finite difference stencils to be applied on non‐staggered grids; (ii) high‐order pressure gradient approximations to be made using standard Padé schemes, and (iii) a variety of boundary conditions to be incorporated in a natural manner. Results are presented in detail for a selection of two‐dimensional steady‐state test problems, using the fourth‐order scheme to demonstrate the accuracy and the robustness of the proposed methods. Furthermore, extensions to higher orders and time‐dependent problems are illustrated, whereas the extension to three‐dimensional problems is also discussed. Copyright © 2010 John Wiley & Sons, Ltd.     

Short Communication: Application of WENO (Weighted Essentially Non-oscillatory) Approach to Navier–Stokes Computations

The direct implementation of the essentially non-oscillatory schemes for flow simulation over complex geometries sometimes results in insufficiently robust numerical algorithms. In order to overcome this difficulty, it is suggested to use the weighted essentially non-oscillatory approach for multidimensional Navier–Stokes computations. The results indicate a significant improvement in accuracy and robustness, especially for low Mach and high supersonic flows.     

Evaluation of linear and nonlinear convection schemes on multidimensional non‐orthogonal grids with applications to KVLCC2 tanker

Almost all evaluations of convection schemes reported in the literature are conducted using simple problems on uniform orthogonal grids; thus, having limited contribution when solving industrial computational fluid dynamics (CFD), where the grids are usually non‐orthogonal with distortions. Herein, several convection schemes are assessed in uniform and distorted non‐orthogonal grids with emphasis on industrial applications. Linear and nonlinear (TVD) convection schemes are assessed on analytical benchmarks in both uniform and distorted grids. To evaluate the performance of the schemes, four error metrics are used: dissipation, phase and L₁ errors, and the schemes' effective order of accuracy. Qualitative and quantitative deterioration of these error metrics as a function of the grid distortion metrics are investigated, and rigorous verifications are performed. Recommendations for effective use of the convection schemes based on the range of grid aspect ratio (AR), expansion ratio (ER) and skewness (Q) are included. A ship hydrodynamics case is studied, involving a Reynolds averaged Navier–Stokes simulation of a bare‐hull KVLCC2 tanker using linear and nonlinear convection schemes coupled with isotropic and anisotropic Reynolds‐stress (ARS) turbulence models using CFDShip‐Iowa v4. Predictions of local velocities and turbulent quantities from the midships to the nominal wake plane are compared with experimental fluid dynamics (EFD), and rigorous verification and validation analyses for integral forces and moments are performed for 0° and 12° drift angles. Best predictions are observed when coupling a second‐order TVD scheme with the anisotropic turbulence model. Further improvements are observed in terms of prediction of the vortical structures for 30° drift when using TVD2S‐ARS coupled with DES. Copyright © 2009 John Wiley & Sons, Ltd.