Multilevel preconditioning on the refined interface and optimal boundary solvers for the Laplace equation

In this paper we propose and analyze some strategies to construct asymptotically optimal algorithms for solving boundary reductions of the Laplace equation in the interior and exterior of a polygon. The interior Dirichlet or Neumann problems are, in fact, equivalent to a direct treatment of the Dirichlet-Neumann mapping or its inverse, i.e., the Poincaré-Steklov (PS) operator. To construct a fast algorithm for the treatment of the discrete PS operator in the case of polygons composed of rectangles and regular right triangles, we apply the Bramble-Pasciak-Xu (BPX) multilevel preconditioner to the equivalent interface problem in theH ^1/2-setting. Furthermore, a fast matrix-vector multiplication algorithm is based on the frequency cutting techniques applied to the local Schur complements associated with the rectangular substructures specifying the nonmatching decomposition of a given polygon. The proposed compression scheme to compute the action of the discrete interior PS operator is shown to have a complexity of the orderO(N log ^q N),q [2, 3], with memory needsO(N log² N), whereN is the number of degrees of freedom on the polygonal boundary under consideration. In the case of exterior problems we propose a modification of the standard direct BEM whose implementation is reduced to the wavelet approximation applied to either single layer or hypersingular harmonic potentials and, in addition, to the matrix-vector multiplication for the discrete interior PS operator.     

Augmented high order finite volume element method for elliptic PDEs in non-smooth domains: Convergence study

The accuracy of a finite element numerical approximation of the solution of a partial differential equation can be spoiled significantly by singularities. This phenomenon is especially critical for high order methods. In this paper, we show that, if the PDE is linear and the singular basis functions are homogeneous solutions of the PDE, the augmentation of the trial function space for the Finite Volume Element Method (FVEM) can be done significantly simpler than for the Finite Element Method. When the trial function space is augmented for the FVEM, all the entries in the matrix originating from the singular basis functions in the discrete form of the PDE are zero, and the singular basis functions only appear in the boundary conditions. That is to say, there is no need to integrate the singular basis functions over the elements and the sparsity of the matrix is preserved without special care. FVEM numerical convergence studies on two-dimensional triangular grids are presented using basis functions of arbitrary high order, confirming the same order of convergence for singular solutions as for smooth solutions.     

Turbulent flow computations on 3D unstructured grids

An edge-based finite element method is presented for the simulation of compressible turbulent flows on unstructured tetrahedral grids. A two equation k–ω turbulence model is employed and the standard Galerkin approach is used for spatial discretisation. Stabilisation of the resulting procedure is achieved by the addition of an appropriate diffusion. An explicit multistage time-stepping scheme is used to advance the solution in time to steady state. The performance of the algorithm is demonstrated for the simulation of a high Reynolds number transonic separated flow over a wing.     

Resolution of computational aeroacoustics problems on unstructured grids with a higher-order finite volume scheme

Computational fluid dynamics (CFD) has become increasingly used in the industry for the simulation of flows. Nevertheless, the complex configurations of real engineering problems make the application of very accurate methods that only work on structured grids difficult. From this point of view, the development of higher-order methods for unstructured grids is desirable. The finite volume method can be used with unstructured grids, but unfortunately it is difficult to achieve an order of accuracy higher than two, and the common approach is a simple extension of the one-dimensional case. The increase of the order of accuracy in finite volume methods on general unstructured grids has been limited due to the difficulty in the evaluation of field derivatives. This problem is overcome with the application of the Moving Least Squares (MLS) technique on a finite volume framework. In this work we present the application of this method (FV-MLS) to the solution of aeroacoustic problems.     

Object-oriented modelling and simulation of heat exchangers with finite element methods

The paper describes the derivation of finite-element models of one-dimensional fluid flows with heat transfer in pipes, using the Galerkin/least-squares approach. The models are first derived for one-phase flows, and then extended to homogeneous two-phase flows. The resulting equations have then been embedded in the context of object-oriented system modelling; this allows one to combine the fluid flow model with a model for other phenomena such as heat transfer, as well as with models of other discrete components such as pumps or valves, to obtain complex models of heat exchangers. The models are then validated by simulating a typical heat exchanger plant.     

On the solution of elliptic partial differential equations on regions with corners II: Detailed analysis

In this paper we investigate the solution of boundary value problems on polygonal domains for elliptic partial differential equations. Previously, we observed that when the boundary value problems are formulated as boundary integral equations of classical potential theory, the solutions are representable by series of elementary functions, to arbitrary order, for all but finitely many values of the angles. Here, we extend this observation to all values of the angles. We show that the solutions near corners are representable, to arbitrary order, by linear combinations of certain non-integer powers and non-integer powers multiplied by logarithms.     

非嵌套有限元的W循环多重网格方法

本文讨论了下述情形：１非嵌套网格；２曲边有限元；３非协调元；４拟协调元；５有限元的型函数有特殊性质，都能导致非嵌套的有限元空间．对一个包括上述情形的问题给出了非嵌套有限元的Ｗ循环多重网格方法，并证明了它的收敛性。     

The Finite Element Method on the Sierpinski Gasket

For certain classes of fractal differential equations on the Sierpinski gasket, built using the Kigami Laplacian, we describe how to approximate solutions using the finite element method based on piecewise harmonic or piecewise biharmonic splines. We give theoretical error estimates, and compare these with experimental data obtained using a computer implementation of the method (available at the web site http://mathlab.cit.cornell.edu/\sim gibbons). We also explain some interesting structure concerning the spectrum of the Laplacian that became apparent from the experimental data. March 29, 2000. Date revised: March 6, 2001. Date accepted: March 21, 2001.