Simulation of 2D external viscous flows by means of a domain decomposition method using an influence matrix technique

Two-dimensional external viscous flows are numerically approximated by means of a domain decomposition method which combines a vortex method and a finite difference method. The vortex method is used in the flow region which is dominated by convective effects, whereas the finite difference method is used in the flow region where viscous diffusion effects are dominant. An influence matrix technique combined with the uniformity condition of the pressure is used to enforce the tangential velocity boundary condition. Comparisons between numerical and experimental data show that the method is well adapted for simulating two-dimensional flows.     

Computations of a laminar backward‐facing step flow at Re=800 with a spectral domain decomposition method

The two‐dimensional laminar incompressible flow over a backward‐facing step is computed using a spectral domain decomposition approach. A minimum number of subdomains (two) is used; high resolution being achieved by increasing the order of the basis Chebyshev polynomial. Results for the case of a Reynolds number of 800 are presented and compared in detail with benchmark computations. Stable accurate steady flow solutions were obtained using substantially fewer nodes than in previously reported simulations. In addition, the problem of outflow boundary conditions was examined on a shortened domain. Because of their more global nature, spectral methods are particularly sensitive to imposed boundary conditions, which may be exploited in examining the effect of artificial (non‐physical) outflow boundary conditions. Two widely used set of conditions were tested: pseudo stress‐free conditions and zero normal gradient conditions. Contrary to previous results using the finite volume approach, the latter is found to yield a qualitatively erroneous yet stable flow‐field. Copyright © 1999 John Wiley & Sons, Ltd.     

不可压粘流N－S方程的边界积分解法

对原变量的Ｎ－Ｓ方程进行一阶时间离散,采用共轭梯度法解除压强－速度的耦合．对所得的一系列Ｌａｐｌａｃｅ方程、Ｐｏｓｓｉｏｎ方程和Ｈｅｌｍｈｏｔｚ方程均进行边界积分法求解,首次得到了粘性Ｎ－Ｓ方程的边界积分表示式．圆柱的定常、非定常尾迹计算结果表明了本文方法的有效性．     

Flame computations by a Chebyshev multidomain method

A Chebyshev collocation method is proposed for the computation of laminar flame propagation in a two-dimensional gaseous medium. The method is based on a domain decomposition technique associated with co-ordinate transforms to map the infinite physical subdomains into finite computational ones. The influence matrix method is used to handle the patching conditions at the interfaces. This technique is particularly efficient since at each time step only matrix products have to be performed. The method is tested first on an elliptic model problem; it is then applied to laminar flame computations, including calculations of cellular instabilities of flame fronts.     

不可压粘流N－S方程的边界积分解法

对原变量的Ｎ－Ｓ方程进行一阶时间离散，采用共轭梯度法解除压强－速度的耦合．对所得的一系列Ｌａｐｌａｃｅ方程、Ｐｏｓｓｉｏｎ方程和Ｈｅｌｍｈｏｔｚ方程均进行边界积分法求解，首次得到了粘性Ｎ－Ｓ方程的边界积分表示式．圆柱的定常、非定常尾迹计算结果表明了本文方法的有效性．     

Aerodynamic and aeroacoustic analysis of a pulsating internal flow by a multidomain weak collocation spectral method

A numerical methodology is presented for the coupled aerodynamic and aeroacoustic analysis of time dependent laminar viscous flows in general two-dimensional geometries. The overall procedure is constituted firstly by the solution of the incompressible Navier-Stokes equations and secondly by the solution of a linear evolution equation of second order in space and time which originates from Lighthill's acoustic analogy. A semi-implicit projection method is utilized for the time integration of the incompressible Navier-Stokes system, while a choice between a second order implicit Newmark scheme and a fourth order explicit Runge-Kutta-Nyström method is offered for the temporal discretization of the wave equation. The multidomain weak Legendre collocation spectral method on quadrilateral subdomain topologies is employed for the spatial approximation of both the fluid dynamic and the acoustic problems. A specific pulsating internal flow inside a plane constricted channel is selected as a representative application for the assessment of the capabilities of the proposed discrete algorithm. Numerous results are presented and discussed, so as to thoroughly demonstrate the behavior of the numerical procedure.     

Inviscid instabilities in two-dimensional spatially periodic flows

The instability of two-dimensional vortex arrays with hexagonal symmetry is investigated with respect to inviscid perturbations. A numerical treatment in the framework of Floquet theory provides critical parameters for instability onset and the spatial structure of unstable modes which are in agreement with experiments in electromagnetically driven flows conducting fluids. The stability of point vortex lattices is briefly discussed.     

Chebyshev collocation method and multidomain decomposition for the incompressible Navier-Stokes equations

The two-dimensional incompressible Navier-Stokes equations in primitive variables have been solved by a pseudospectral Chebyshev method using a semi-implicit fractional step scheme. The latter has been adapted to the particular features of spectral collocation methods to develop the monodomain algorithm. In particular, pressure and velocity collocated on the same nodes are sought in a polynomial space of the same order; the cascade of scalar elliptic problems arising after the spatial collocation is solved using finite difference preconditioning. With the present procedure spurious pressure modes do not pollute the pressure field. As a natural development of the present work a multidomain extent was devised and tested. The original domain is divided into a union of patching sub-rectangles. Each scalar problem obtained after spatial collocation is solved by iterating by subdomains. For steady problems a C¹ solution is recovered at the interfaces upon convergence, ensuring a spectrally accurate solution. A number of test cases have been solved to validate the algorithm in both its single-block and multidomain configurations. The preliminary results achieved indicate that collocation methods in multidomain configurations might become a viable alternative to the spectral element technique for accurate flow prediction.     

配置点谱方法-人工压缩法(SCM-ACM)求解不可压缩流体流动

开发了配置点谱方法SCM（spectral collocation method）与人工压缩法ACM（artificial compressibility method）相结合的方法SCM-ACM,用于求解不可压缩粘性流动问题。选取典型的方腔顶盖驱动流为研究测试对象,首先建立人工压缩格式的控制方程,其次采用SCM离散控制方程的空间偏微分项,推导出矩阵形式的代数方程,最后测试了SCM-ACM代码的有效性。结果显示,SCM-ACM能够有效求解不可压缩流动问题,并继承了谱方法的指数收敛特性,且具有ACM求解过程简单及易于实施的特点。     

Multidomain pseudospectral methods for nonlinear convection-diffusion equations

Multidomain pseudospectral approximations to nonlinear convection-diffusion equations are considered. The schemes are formulated with the Legendre-Galerkin method, but the nonlinear term is collocated at the Legendre/Chebyshev-Gauss-Lobatto points inside each subinterval. Appropriate base functions are introduced so that the matrix of the system is sparse, and the method can be implemented efficiently and in parallel. The stability and the optimal rate of convergence of the methods are proved. Numerical results are given for both the single domain and the multidomain methods to make a comparison.     

不规则区域平面弹性问题的正则区域重心插值配点法

将不规则区域嵌入到规则的矩形区域,在矩形区域上将弹性平面问题的控制方程采用重心Lagrange插值离散,得到控制方程矩阵形式的离散表达式。在边界节点上利用重心插值离散边界条件,规则区域采用置换法施加边界条件,不规则区域采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,采用最小二乘法进行求解,得到整个规则区域上的位移数值解。利用重心插值计算得到不规则区域内任意节点的位移值,计算精度可到10^-14以上。数值算例验证了所建立方法的有效性和计算精度。     

配置点谱方法-人工压缩法(SCM-ACM)求解同心圆筒内流体流动

发展了配置点谱方法SCM(Spectral collocation method)和人工压缩法ACM(Artificial compressibility method)相结合的SCM-ACM数值方法,计算了柱坐标系下稳态不可压缩流动N-S方程组。选取典型的同心圆筒间旋转流动Taylor-Couette流作为测试对象,首先,采用人工压缩法获得人工压缩格式的非稳态可压缩流动控制方程;再将控制方程中的空间偏微分项用配置点谱方法进行离散,得到矩阵形式的代数方程;编写了SCM-ACM求解不可压缩流动问题的程序;最后,通过与公开发表的Taylor-Couette流的计算结果对比,验证了求解程序的有效性。结果证明,本文发展的SCM-ACM数值方法能够用于求解圆筒内不可压缩流体流动问题,该方法既保留了谱方法指数收敛的特性,也具有ACM形式简单和易于实施的特点。本文发展的SCM-ACM数值方法为求解柱坐标下不可压缩流体流动问题提供了一种新的选择。     

Simulation of electrically driven jet using Chebyshev collocation method

The model of electrically driven jet is governed by a series of quasi 1D dimensionless partial differential equations (PDEs). Following the method of lines, the Chebyshev collocation method is employed to discretize the PDEs and obtain a system of differential-algebraic equations (DAEs). By differentiating constrains in DAEs twice, the system is transformed into a set of ordinary differential equations (ODEs) with invariants. Then the implicit differential equations solver “ddaskr” is used to solve the ODEs and post-stabilization is executed at the end of each step. Results show the distributions of radius, linear charge density, stretching ratio and also the horizontal velocity at a time point. Meanwhile, the spiral and expanding projections to X-Y plane of the jet centerline suggest the occurring of bending instability.     

不可压缩机翼绕流的有限谱法计算

结合有限谱QUICK格式求解不可压缩粘性流问题。这一格式用于模拟不同攻角下的NACA1200机翼绕流问题。利用体积力,提出了将流场速度从0加速到来流速度的方法。区别于传统的压力梯度为零的边界条件,推导出一个更精确的压力边界条件。为使速度散度保持为零,在泊松方程中给速度散度一个特殊的处理。这一成果说明了有限谱法不但具有很高的精度,而且能灵活地和其他格式一起构造出新的格式,从而成功地应用到复杂流场不可压缩流动的数值计算中。     

虚拟区域法在流固耦合问题中的应用

将多相流领域内的虚拟区域法引入到流固耦合问题的分析中,将固体视为应变率为零的虚拟流体,对流体和虚拟流体均以速度和压强作为基本变量,采用Navier-Stokes方程作为控制方程,同时求解流体域和虚拟流体域, 得到整个计算域的流场分布,应用分布式拉格朗日乘子法在虚拟流体域上施加刚体约束, 以保持虚拟流体的刚体外形和运动形式,最终建立一种流固耦合模型及其数值求解方法. 通过对粒子流问题和流固耦合问题进行数值模拟,验证了此模型的正确性和求解大变形/运动流固耦合问题的有效性.