Calculation of impulsively started incompressible viscous flows

The paper's focus is the calculation of unsteady incompressible 2D flows past airfoils. In the framework of the primitive variable Navier–Stokes equations, the initial and boundary conditions must be assigned so as to be compatible, to assure the correct prediction of the flow evolution. This requirement, typical of all incompressible flows, viscous or inviscid, is often violated when modelling the flow past immersed bodies impulsively started from rest. Its fulfillment can however be restored by means of a procedure enforcing compatibility, consisting in a pre‐processing of the initial velocity field, here described in detail. Numerical solutions for an impulsively started multiple airfoil have been obtained using a finite element incremental projection method. The spatial discretization chosen for the velocity and pressure are of different order to satisfy the inf–sup condition and obtain a smooth pressure field. Results are provided to illustrate the effect of employing or not the compatibility procedure, and are found in good agreement with those obtained with a non‐primitive variable solver. In addition, we introduce a post‐processing procedure to evaluate an alternative pressure field which is found to be more accurate than the one resulting from the projection method. This is achieved by considering an appropriate ‘unsplit’ version of the momentum equation, where the velocity solution of the projection method is substituted. Copyright © 2004 John Wiley & Sons, Ltd.     

A fourth‐order compact finite difference scheme for the steady stream function–vorticity formulation of the Navier–Stokes/Boussinesq equations

A fourth‐order compact finite difference scheme on the nine‐point 2D stencil is formulated for solving the steady‐state Navier–Stokes/Boussinesq equations for two‐dimensional, incompressible fluid flow and heat transfer using the stream function–vorticity formulation. The main feature of the new fourth‐order compact scheme is that it allows point‐successive overrelaxation (SOR) or point‐successive underrelaxation iteration for all Rayleigh numbers Ra of physical interest and all Prandtl numbers Pr attempted. Numerical solutions are obtained for the model problem of natural convection in a square cavity with benchmark solutions and compared with some of the accurate results available in the literature. Copyright © 2003 John Wiley & Sons, Ltd.     

GOLDEN AGES OR MAGIC MOMENTS?

ABSTRACT. The paper is an attempt to refute Quinn's ‘the Golden Age of fisheries population models has ended “hypothesis. The approach is to show it is possible to provide simple matrix based theory for use with delay difference equations applied to products of the 0th to 4th powers of length and numbers at length (Alias” proto‐moments“). It compares these to models just based upon numbers and biomass. It also shows examples (of comparatively little moment) indicating how using multiple proto‐moments could improve the ability of delay difference equations to handle the size dependent processes of maturity, selection, natural mortality rates and predation rates and thus improving the biological reality of these models.     

Temperature factor effect on the structure of the separated flow within a supersonic gas stream

The effect of the temperature factor, that is, the ratio of the body temperature to the freestream stagnation temperature, on the structure of the separated flow formed in the presence of a concave corner in a supersonic stream is studied. The strong influence of the temperature factor on the separation zone length and the flow-generated aerodynamic characteristics is established. It is shown that for fairly large deflection angles this flow cannot be described by free interaction, or triple deck, theory.     

Parametric Centers for Trigonometric Abel Equations

This article is devoted to one-dimensional perturbative theory on R × S ¹. There is a recursive formula for the successive obstructions to parametric center at any order of the perturbation parameter. The first obstruction is studied by means of complex analysis techniques. This extends to the trigonometric case what was done previously for the polynomial case (Israel J. Math. 142, 273–283, 2004). This article is dedicated to Professor Zhang Zhi-Fen on the occasion of her 80th Birthday     

Dispersion analysis of the least‐squares finite‐element shallow‐water system

The frequency or dispersion relation for the least‐squares mixed formulation of the shallow‐water equations is analysed. We consider the use of different approximation spaces corresponding to co‐located and staggered meshes, respectively. The study includes the effect of Coriolis, and the dispersion properties are compared analytically and graphically with those of the mixed Galerkin formulation. Numerical solutions of a test problem to simulate slow Rossby modes illustrate the theoretical results. Copyright © 2003 John Wiley & Sons, Ltd.     

The LEM exponential integrator for advection–diffusion–reaction equations

We implement a second-order exponential integrator for semidiscretized advection–diffusion–reaction equations, obtained by coupling exponential-like Euler and Midpoint integrators, and computing the relevant matrix exponentials by polynomial interpolation at Leja points. Numerical tests on 2D models discretized in space by finite differences or finite elements, show that the Leja–Euler–Midpoint (LEM) exponential integrator can be up to 5 times faster than a classical second-order implicit solver.     

粘性／非粘性方程有限元和谱方法混合逼近

讨论有粘性和无粘性Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程有限元和谱方法混合逼近问题，弱交面匹配关系被应用。收敛估计显示有限元精度和谱精度可在各自的子区域中得到。     

“解方程器”软件包(EQNS-SOLVER)

作者提出了对于任意多项式方程组的一种带有普遍性的求解方法。多项式方程组可能没有解,或者有无穷多个解,或者有有限个解,解是实数也可以是复数。在有限个解的情况下,所创建的方法能够求出所有的解。作者完成了一个软料包,称为《解方程器》,曾在微型计算机上用FORTRAN语言完成计算过程。本文是《解方程器》软料包的说明,内容包括基本原理、方法、程序结构、应用举例以及使用的步骤。     

在对称群Sm中解循环方程组

本文研究了在对称群S_m中解循环方程组:x_mX_(m+1)…x_(ω+k-1)σ_m,其中σ_m∈S_m,ω=1,2,…,s,确定了这个循环方程组有解的充分必要条件和解的个数,而且给出了求解过程。