定常的热传导-对流问题的Galerkin/Petrov最小二乘混合元方法

1.引言 热传导-对流问题是大气动力学中的一个重要的方程,这个方程组也称为强迫耗散的非线性系统方程组,其较Navier-Stokes方程多了一个未知函数温度场,且温度与速度和压力之间存在着复杂的非线性关系.从热动力学可知,任何运动都会产生热量即有温度,而且温度与速度和压力之间必定互相转化,因此对该非线性系统的研究更具有实际意义.[1]先对     

Difference schemes on uniform grids performed by general discrete operators

Our aim is to set the foundations of a discrete vectorial calculus on uniform n-dimensional grids, that can be easily reformulated on general irregular grids. As a key tool we first introduce the notion of tangent space to any grid node. Then we define the concepts of vector field, field of matrices and inner products on the space of grid functions and on the space of vector fields, mimicking the continuous setting. This allows us to obtain the discrete analogous of the basic first order differential operators, gradient and divergence, whose composition define the fundamental second order difference operator. As an application, we show that all difference schemes, with constant coefficients, for first and second order differential operators with constant coefficients can be seen as difference operators of the form for suitable choices of q, and  . In addition, we characterize special properties of the difference scheme, such as consistency, symmetry and positivity in terms of q, and  .     

Stability analysis of abnormal multiplication of plankton using parameter identification technique

This paper presents the mathematical approach for the abnormal multiplication of plankton. An abnormal multiplication can be expressed as an unstable problem and the stability of the system is investigated by introducing eigenvalues of a mathematical equation. The stability of the system can be judged by an eigenvalue based on the Lyapunov's stability theory. In this paper, the Arnoldi‐QR method is used to obtain eigenvalues and eigenvectors of the system. The mode superposition method is employed to create spatial distribution needed to analyse the stability. To obtain the objective eigenvalue, the parameter identification technique is employed. The finite element method is used for the discretization in space. Lake Kasumigaura, which is located in Ibaraki Prefecture in Japan, is selected and actual data in 1975, 1976, 1991 and 2000 are used in order to investigate the stability of the specified lake in Japan. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical analysis of the rotating viscous flow approaching a solid sphere

A numerical simulation is performed to investigate the flow induced by a sphere moving along the axis of a rotating cylindrical container filled with the viscous fluid. Three‐dimensional incompressible Navier–Stokes equations are solved using a finite element method. The objective of this study is to examine the feature of waves generated by the Coriolis force at moderate Rossby numbers and that to what extent the Taylor–Proudman theorem is valid for the viscous rotating flow at small Rossby number and large Reynolds number. Calculations have been undertaken at the Rossby numbers (R_o) of 1 and 0.02 and the Reynolds numbers (R_e) of 200 and 500. When R_o=O(1), inertia waves are exhibited in the rotating flow past a sphere. The effects of the Reynolds number and the ratio of the radius of the sphere and that of the rotating cylinder on the flow structure are examined. When R_o ? 1, as predicted by the Taylor–Proudman theorem for inviscid flow, the so‐called ‘Taylor column’ is also generated in the viscous fluid flow after an evolutionary course of vortical flow structures. The initial evolution and final formation of the ‘Taylor column’ are exhibited. According to the present calculation, it has been verified that major theoretical statement about the rotating flow of the inviscid fluid may still approximately predict the rotating flow structure of the viscous fluid in a certain regime of the Reynolds number. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulation of single‐ and two‐phase flows on sliding unstructured meshes using finite volume method

The paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.     

Fully adaptive multiresolution finite volume schemes for conservation laws

The use of multiresolution decompositions in the context of finite volume schemes for conservation laws was first proposed by A. Harten for the purpose of accelerating the evaluation of numerical fluxes through an adaptive computation. In this approach the solution is still represented at each time step on the finest grid, resulting in an inherent limitation of the potential gain in memory space and computational time. The present paper is concerned with the development and the numerical analysis of fully adaptive multiresolution schemes, in which the solution is represented and computed in a dynamically evolved adaptive grid. A crucial problem is then the accurate computation of the flux without the full knowledge of fine grid cell averages. Several solutions to this problem are proposed, analyzed, and compared in terms of accuracy and complexity.

     

低速二维流动的粒子仿真研究

本文运用信息保存法对低速二维的流动现象进行模拟，考察了低速条件下的有限平板绕流以及微槽道气体流动问题。研究表明：在对低速流动的模拟过程中，运用IP法在能够获得较好的结果的同时，具有比DSMC方法更高的计算效率。     

Contractivity and Analyticity in l p of Some Approximation of the Heat Equation

We consider the lumped mass method with piecewise linear finite elements in two dimensions. When the triangulation is of Delaunay type it is known that the discrete scheme satisfies a maximum principle. In this work we pursue the analysis and prove that the discrete semi-group is l ^p contractive in a sector, which implies smoothing effects and provide some resolvent estimates.     

Efficient discretization of pressure‐correction equations on non‐orthogonal grids

An alternative discretization of pressure‐correction equations within pressure‐correction schemes for the solution of the incompressible Navier–Stokes equations is introduced, which improves the convergence and robustness properties of such schemes for non‐orthogonal grids. As against standard approaches, where the non‐orthogonal terms usually are just neglected, the approach allows for a simplification of the pressure‐correction equation to correspond to 5‐point or 7‐point computational molecules in two or three dimensions, respectively, but still incorporates the effects of non‐orthogonality. As a result a wide range (including rather high values) of underrelaxation factors can be used, resulting in an increased overall performance of the underlying pressure‐correction schemes. Within this context, a second issue of the paper is the investigation of the accuracy to which the pressure‐correction equation should be solved in each pressure‐correction iteration. The scheme is investigated for standard test cases and, in order to show its applicability to practical flow problems, for a more complex configuration of a micro heat exchanger. Copyright © 2003 John Wiley & Sons, Ltd.