Assessment of some iterative methods for non-symmetric linear systems arising in computational fluid dynamics

Various tests have been carried out in order to compare the performances of several methods used to solve the non-symmetric linear systems of equations arising from implicit discretizations of CFD problems, namely the scalar advection-diffusion equation and the compressible Euler equations. The iterative schemes under consideration belong to three families of algorithms: relaxation (Jacobi and Gauss-Seidel), gradient and Newton methods. Two gradient methods have been selected: a Krylov subspace iteration method (GMRES) and a non-symmetric extension of the conjugate gradient method (CGS). Finally, a quasi-Newton method has also been considered (Broyden). The aim of this paper is to provide indications of which appears to be the most adequate method according to the particular circumstances as well as to discuss the implementation aspects of each scheme.     

Front roughening in three-dimensional imbibition

We investigate the structure and dynamics of the interface between two immiscible liquids in a three-dimensional disordered porous medium. We apply a phase-field model that includes explicitly disorder and discuss both spontaneous and forced imbibition. The structure of the interface is dominated by a length scale ξ_× which arises from liquid conservation. We further show that disorder in the capillary and permeability act on different length scales and give rise to different scalings and structures of the interface properties. We conclude with a range of applications.     

一类数论函数的性质(Ⅱ)

For any given positive integer n≥1,the Euler functionφ(n) is defined to be the number of positive integers not exceeding n,which is relatively prime to n.ω(n) is defined to be the number of different prime divisors of n.In order to know the solvability of the function ofφ(φ(φ(n)))=2~(ω(n)),properties of the number theoretical functionφ(φ(n)) is studied in the paper.     

Lattice Boltzmann方法模拟多孔介质Klinkenberg效应

格子Boltzmann数值模拟方法是研究复杂的多孔介质结构特别是Klinkenberg效应的有效方法之一，对处理复杂边值问题尤其有效，用格子Boltzmann方法研究了气流穿越多孔介质问题，并将数值计算结果与实验结果进行了比较，结果表明格子Boltzmann方法是数值模拟气流穿越多孔介质问题的有效方法之一。     

Particle systems on flows

We start with a stochastic flow of diffeomorphisms of the space. Particles enter the space at random times and places. Each particle is carried by the flow for some random amount of time. We examine the point process formed by the particles at a fixed time, on the evolution of that point process as time varies, and on the equilibrium law of the point process.     

Equivalence of u–p and ζ–ψ formulations of the time-dependent Navier–Stokes equations

This paper deals with the non-stationary incompressible Navier--Stokes equations for two-dimensional flows expressed in terms of the velocity and pressure and of the vorticity and streamfunction. The equivalence of the two formulations is demonstrated, both formally and rigorously, by virtue of a condition of compatibility between the boundary and initial values of the normal component of velocity. This condition is shown to be the only compatibility condition necessary to allow for solutions of a minimal regularity, namely H¹ for the velocity, as in most current numerical schemes relying on spatial discretizations of local type.     

沉浮,俯仰及控制面振动引起的跨音速非定常气动载荷计算

本文将文献[9]提出改进的通量分裂方法,应用于随时间变化的贴体网格中,建立了可用于求解非定常Euler方程的通量分裂方法.该方法是以连续的特征值分离为基础,它具有方法简单,便于推广使用的特点.同时克服了Steger-Warming通量分裂方法存在的问题.对通量分裂后的Euler方程.利用MUSCL型迎风差分建立了具有二阶精度的有限体积方程.文中以NACA64A—10翼型为例,对其在跨音速流场中进行沉浮、俯仰及带有振动控制面引起的非定常气动载荷进行了计算.部分计算结果与相应的实验结果进行了比较,吻合良好     

Analysis of plane and axisymmetric flows of incompressible fluids with the stream tube method: Numerical simulation by trust-region optimization algorithm

New concepts for the study of incompressible plane or axisymmetric flows are analysed by the stream tube method. Flows without eddies and pure vortex flows are considered in a transformed domain where the mapped streamlines are rectilinear or circular. The transformation between the physical domain and the computational domain is an unknown of the problem. In order to solve the non-linear set of relevant equations, we present a new algorithm based on a trust region technique which is effective for non-convex optimization problems. Experimental results show that the new algorithm is more robust compared to the Newton-Raphson method.     

A theorem on integer flows on cartesian products of graphs

It is shown that the Cartesian product of two nontrivial connected graphs admits a nowhere‐zero 4‐flow. If both factors are bipartite, then the product admits a nowhere‐zero 3‐flow. © 2003 Wiley Periodicals, Inc. J Graph Theory 43: 93–98, 2003