Stabilization Methods for Spectral Element Computations of Incompressible Flows

A stabilization method for the spectral element computation of incompressible flow problems is investigated. It is based on a filtering procedure which consists in filtering the velocity field by a spectral vanishing Helmholtz-type operator at each time step. Relationship between this filtering procedure and SVV-stabilization method, introduced recently in [JCP, 2004, 196(2), p680], is established. A number of numerical examples are presented to show the accuracy and stabilization capability of the method.     

Numerical Study of a Modified Time-Stepping θ-Scheme for Incompressible Flow Simulations

In [Turek (1996). Int. J. Numer. Meth. Fluids 22, 987–1011], we had performed numerical comparisons for different time stepping schemes for the incompressible Navier–Stokes equations. In this paper, we present the numerical analysis in the context of the Navier–Stokes equations for a modified time-stepping θ-scheme which has been recently proposed by Glowinski [Glowinski (2003). In: Ciarlet, P. G., and Lions, J. L. (eds.), Handbook of Numerical Analysis, Vol. IX, North-Holland, Amsterdam, pp. 3–1176]. Like the well-known classical Fractional-Step-θ-scheme which had been introduced by Glowinski [Glowinski (1985). In Murman, E. M. and Abarbanel, S. S. (eds.), Progress and Supercomputing in Computational Fluid Dynamics, Birkh?user, Boston MA; Bristeau et al. (1987). Comput. Phys. Rep. 6, 73–187], too, and which is still one of the most popular time stepping schemes, with or without operator splitting techniques, this new scheme consists of 3 substeps with nonequidistant substepping to build one macro time step. However, in contrast to the Fractional-Step-θ-scheme, the second substep can be formulated as an extrapolation step for previously computed data only, and the two remaining substeps look like a Backward Euler step so that no expensive operator evaluations for the right hand side vector with older solutions, as for instance in the Crank–Nicolson scheme, have to be performed. This modified scheme is implicit, strongly A-stable and second order accurate, too, which promises some advantageous behavior, particularly in implicit CFD simulations for the nonstationary Navier–Stokes equations. Representative numerical results, based on the software package FEATFLOW [Turek (2000). FEATFLOW Finite element software for the incompressible Navier–Stokes equations: User Manual, Release 1.2, University of Dortmund] are obtained for typical flow problems with benchmark character which provide a fair rating of the solution schemes, particularly in long time simulations.Dedicated to David Gottlieb on the occasion of his 60th anniversary     

Implicit–Explicit Schemes for BGK Kinetic Equations

In this work a new class of numerical methods for the BGK model of kinetic equations is presented. In principle, schemes of any order of accuracy in both space and time can be constructed with this technique. The methods proposed are based on an explicit–implicit time discretization. In particular the convective terms are treated explicitly, while the source terms are implicit. In this fashion even problems with infinite stiffness can be integrated with relatively large time steps. The conservation properties of the schemes are investigated. Numerical results are shown for schemes of order 1, 2 and 5 in space, and up to third-order accurate in time.     

Viscoelastic State for the Solutions of Spherically Symmetric Viscoelasticity Problems

The solutions of the spherically symmetric, linear, isothermal, and transient viscoelasticity problems via reciprocity theorem have been investigated for a specific material. The integral form of stress–strain relations has been used. The Laplace transform of a viscoelastic state, which is necessary for the integral equation arising as a result of reciprocity theorem, has been derived. This integral equation has been solved by Laplace transform. A sample problem has been solved to test the presented formulation. A numerical application of the analytic solution of this problem has been given.     

一种复合型数值方法的改进与算法实现

文中对原有的复合方法进行了改进,提出了一种新的复合型数值方法,这种方法解决了原方法求解含源项的问题时存在数值震荡的缺点。为了验证算法的有效性,对计算流体力学中一阶浅水波问题的两个典型的算例进行了数值模拟,在解的光滑性、锐利性等几个方面做了分析,结果表明改进后的方法既消除了数值震荡现象,又能准确地逼近真解,求出锐利的间断,是一种行之有效的数值方法。     

Experimental and numerical investigation of fracture in double-edge notched steel plates

Double-edge notched (DENT) steel plates were pulled until complete fracture and several experimental observations were made (using profilometry and scanning electron microscopy). The essential work of fracture (EWF) model was found to be well verified. Numerical simulations—up to the maximum load only—of some experiments were performed using the finite element method (FEM), and incorporating geometric and material non-linearities (large deformation elasto-plasticity). Some experimental measurements were compared with the corresponding numerical computations and excellent agreement was found.     

用于非均匀介质重建的选代算法收敛性的研究

研究轴对称二维非均匀介质重建迭代算法的收敛性问题。首先借助于玻恩近似将非线性积分方程线性化，然后，利用玻恩迭代法和变形玻恩迭代法，由散射场数据对轴对称二维非均匀介质进行重建。通过几个例子研究了影响迭代算法收敛性的几种因素，如送代算法本身、积分离散化后的网格划分、正则化方法中的正则化参数选取等。     

泥沙输移矢量式及其应用

本文基于Parker等人泥沙输移矢量式的分析和研究，增加了上举力的作用，并考虑了泥沙颗粒受力非平衡条件，采用牛顿力学的矢量分析法，推导了泥沙输移的通用矢量式。经对公式简化，且用Ikeda和笔者试验资料的验证，证明计算值与试验值吻合较好。     

RADIATION FIELD FROM A RECTANGULAR PLAQUE SOURCE FOR RADIATION PROCESSING

Evaluation of the Hubbell integral for an omnidirectional detector response to a rectangular plaque gamma emitting radiation source using a different numerical computation method on a personal computer was made. For a plaque consisting of several strip gamma emitters, exposure rates at different positions above the source were calculated using 'exact' evaluation and an approximate method. Optimization method is used to find source activity of each strip, so as to acquire desirable dose uniformity of the irradiated products. In uses of the Hubbell integral, by introducing appropriate programming, computation time can be saved.     

Affine Dependence Classification for Communications Minimization

This paper introduces results on placement and communications minimization for systems of affine recurrence equations. We show how to classify the dependences according to the number and nature of communications they may result in. We give both communication-free conditions and conditions for an efficient use of broadcast or neighbor-to-neighbor communication primitives. Since the dependences of a problem can generally not be all communication-free, we finally introduce a heuristic to globally minimize the communications based on the classification of dependences.