不可压缩流动的并行数值方法

不可压缩流动的数值模拟是计算流体力学的重要组成部分. 基于有限元离散方法, 本文设计了不可压缩Navier-Stokes (N-S)方程支配流的若干并行数值算法. 这些并行算法可归为两大类: 一类是基于两重网格离散方法, 首先在粗网格上求解非线性的N-S方程, 然后在细网格的子区域上并行求解线性化的残差方程, 以校正粗网格的解; 另一类是基于新型完全重叠型区域分解技巧, 每台处理器用一局部加密的全局多尺度网格计算所负责子区域的局部有限元解. 这些并行算法实现简单, 通信需求少, 具有良好的并行性能, 能获得与标准有限元方法相同收敛阶的有限元解. 理论分析和数值试验验证了并行算法的高效性     

服从OldroydB型微分模型的粘弹性流体问题的数值解法

本文对服从OldroydB型微分模型的粘弹性流体问题给出了一种数值逼近算法．该算法对压力方程采用标准混合有限元方法,对速度方程采用并行非重叠区域分解方法和特征线法．这种并行算法在子区域上用Galerkin方法,通过积分平均方法显式地给出内边界的数值流．在本文最后还给出了该算法的最优L^2。一误差估计．     

钝体分离旋涡流动的区域分解、杂交数值模拟——Ⅰ.理论方法及其应用

为克服涡旋法不能精确预计物体附近小尺度流动结构的理论缺陷,减少高Reynolds数流动N-S方程差分解的困难,本文提出一种区域分解、杂交耦合N-S方程有限差分解及涡旋法的新的数值模型和理论方法.将流场分解为内外两区,在靠近物体表面、范围为O(R)的内区进行N-S方程有限差分解,外区作Lagrange-Euler涡旋法解,建立了分区流动的联结、耦合条件,给出了杂交耦合求解的数值计算方法.用本方法作了Re=10²,10³的圆柱绕流计算,考察了区域交界面位置变化时解的稳定性.与全场N-S方程解及实验结果的比较表明本文方法能精确预计流动分离及近场流动的详细结构,并可有效地计算流动的总体特性,且比全场N-S方程解显著节省机时和计算量.     

外部绕流的Navier-Stokes方程的边界层方程和维数分裂法

采用基于物体表面二维曲面的半测地坐标系(S-coordinate)建立了一个新的外部绕流边界层方程(boundary layer equations,BLE).BLE是一个关于物体的未知法向粘性应力张量和压力的非线性偏微分方程,其解的存在性得到了证明.此外,通过在二维流形上应用若干个2D-3C偏微分方程组来近似Navier-Stokes方程,获得了三维Navier-Stokes方程的维数分裂法.最后,对球和椭球的外部绕流问题给出了算例.     

透平机械内部三维可压缩Navier-Stokes方程的流面方法和维数分裂算法

在这篇文章中,运用经典的张量分析方法,把流动区域用-个二维流形序列分割成一系列流层之并,推得在流层内半测地坐标之下的Navier-Stokes方程,在流形的法线方向应用向后Euler差分,推导了两维流形上的可压缩Navier-Stokes方程,和流函数满足的方程.在这个基础上,提出了一种维数分裂法的新算法.这种方法不同于区域分解法.对于三维问题,在区域分解法中我们必须在每个子区域上仍解三维问题,但是在这种新方法中,只需要在每个子区域上求解二维问题,不过是几个二维流形上的NS方程.文中还给出了-个透平机械内部流动的数值计算实例.     

环形扩压通道内有旋绕流动的研究

本文计算了环形截面的扩压通道内带进气旋绕的流动.在小横向流假定下.用三维边界层积分方程法求解内外壁面附近的流动.通过对子午面上与流线子午投影准正交方向的速度梯度方程和流量不变方程的迭代求解得出边界层外的势流场.计算与实验结果基本符合.本研究可用于分析环形扩压器内带进气予旋的流动.     

基于对称化原理的对称区域分裂法

本文利用对称化原理,讨论了一种只需在子区域上计算两个完全独立子问题就可得到原问题解的对称区域分裂法,并用此方法求解线性算子方程和线性透射问题.此方法可作为并行算法在MIMD计算机上使用.     

定常Stokes方程一种基于完全区域分解的有限元并行算法

基于完全区域分解技巧,提出了一种求解定常Stokes方程的有限元并行算法．该算法中,所有子问题都是定义在整个求解区域上,但绝大部分自由度来自其所负责的子区域,从而使得算法稍加修改现有的串行程序即可实现相应的并行计算,实现简单,通信需求少．数值结果验证了算法的高效性．     

不可压流体的边界层问题

研究三维有界区域在边界上有流动的不可压流体的边界层问题,导出了Navier-Stokes方程区域内部的近似方程(Euler方程和线性化的Euler方程)和边界附近近似的方程(零阶边界层方程与一阶边界层方程),证明了这种近似的合理性.     

耗散MKdV方程的指数吸引子

讨论了无界区域R~1上的MKdV方程,运用带权空间构造一类紧算子和算子分解的方法,得到该方程在H~2（R~1）上指数吸引子的存在性.     

Dimension splitting method for the three dimensional rotating Navier-Stokes equations

In this paper, we propose a dimensional splitting method for the three dimensional (3D) rotating Navier-Stokes equations. Assume that the domain is a channel bounded by two surfaces and is decomposed by a series of surfaces ■i into several sub-domains, which are called the layers of the flow. Every interface i between two sub-domains shares the same geometry. After establishing a semi-geodesic coordinate (S-coordinate) system based on ■i , Navier-Stoke equations in this coordinate can be expressed as the sum of two operators, of which one is called the membrane operator defined on the tangent space on ■i , another one is called the bending operator taking value in the normal space on ■i . Then the derivatives of velocity with respect to the normal direction of the surface are approximated by the Euler central difference, and an approximate form of Navier-Stokes equations on the surface ■i is obtained, which is called the two-dimensional three-component (2D-3C) Navier-Stokes equations on a two dimensional manifold. Solving these equations by alternate iteration, an approximate solution to the original 3D Navier-Stokes equations is obtained. In addition, the proof of the existence of solutions to 2D-3C Navier-Stokes equations is provided, and some approximate methods for solving 2D-3C Navier-Stokes equations are presented.     

A novel variational method for 3D viscous flow in flow channel of turbomachines based on differential geometry

ABSTRACT

This paper presents a novel variational method for treating three-dimensional rotational Navier-Stokes equations in flow channel of turbomachines. The proposed method establishes a new semi-geodesic coordinate system on the central surface of blades. From the perspective of differential geometry, the system under concern is split into a set of membrane operator equations on two-dimensional manifolds and bending operator equations along hub circle. The third variable of the new coordinate system is approximated by the central difference scheme. We derive a new formulation of two-dimensional Navier-Stokes equations with three components on the manifolds in the variational sense. The well-posedness of the proposed variational formulation is rigorously justified.     

Time stepping for vectorial operator splitting

We present a fully implicit finite difference method for the unsteady incompressible Navier-Stokes equations. It is based on the one-step θ-method for discretization in time and a special coordinate splitting (called vectorial operator splitting) for efficiently solving the nonlinear stationary problems for the solution at each new time level. The resulting system is solved in a fully coupled approach that does not require a boundary condition for the pressure. A staggered arrangement of velocity and pressure on a structured Cartesian grid combined with the fully implicit treatment of the boundary conditions helps us to preserve the properties of the differential operators and thus leads to excellent stability of the overall algorithm. The convergence properties of the method are confirmed via numerical experiments.     

Rational Chebyshev collocation method for the similarity solution of two dimensional stagnation point flow

In this study, we propose an efficient and accurate numerical technique that is called the rational Chebyshev collocation (RCC) method to solve the two dimensional flow of a viscous fluid in the vicinity of a stagnation point named Hiemenz flow. The Navier-Stokes equations governing the flow, are reduced to a third-order ordinary differential equation of a boundary value problem with a semi-infinite domain by using similarity transformation. The rational Chebyshev method reduces this nonlinear ordinary differential equation to a system of algebraic equations. This technique is a powerful type of the collocation methods for solving the boundary value problems over a semi-infinite interval without truncating it to a finite domain. We also present the comparison of this work with others and show that the present method is more accurate and efficient.     

Global strong solutions of Navier-Stokes equations with interface boundary in three-dimensional thin domains

In this article, we study the spectrum of the Stokes operator in a 3D two layer domain with interface, obtain the asymptotic estimates on the spectrum of the Stokes operator as thickness ε goes to zero. Based on the spectral decomposition of the Stokes operator, a new average-like operator is introduced and applied to the study of Navier-Stokes equation in the two layer thin domains under interface boundary condition. We prove the global existence of strong solutions to the 3D Navier-Stokes equations when the initial data and external forces are in large sets as the thickness of the domain is small. This article is a continuation of our study on the Stokes operator under Navier friction boundary condition. Due to the viscosity distinction between the two layers, the Stokes operator displays radically different spectral structure from that under Navier friction boundary condition, then causes great difficulty to the analysis.     

Navier-Stokes equations in 3D thin domains with Navier friction boundary condition

In this article we study the 3D Navier-Stokes equations with Navier friction boundary condition in thin domains. We prove the global existence of strong solutions to the 3D Navier-Stokes equations when the initial data and external forces are in large sets as the thickness of the domain is small. We generalize the techniques developed to study the 3D Navier-Stokes equations in thin domains, see [G. Raugel, G. Sell, Navier-Stokes equations on thin 3D domains I: Global attractors and global regularity of solutions, J. Amer. Math. Soc. 6 (1993) 503-568; G. Raugel, G. Sell, Navier-Stokes equations on thin 3D domains II: Global regularity of spatially periodic conditions, in: Nonlinear Partial Differential Equations and Their Application, College de France Seminar, vol. XI, Longman, Harlow, 1994, pp. 205-247; R. Temam, M. Ziane, Navier-Stokes equations in three-dimensional thin domains with various boundary conditions, Adv. Differential Equations 1 (1996) 499-546; R. Temam, M. Ziane, Navier-Stokes equations in thin spherical shells, in: Optimization Methods in Partial Differential Equations, in: Contemp. Math., vol. 209, Amer. Math. Soc., Providence, RI, 1996, pp. 281-314], to the Navier friction boundary condition by introducing a new average operator Mε in the thin direction according to the spectral decomposition of the Stokes operator Aε. Our analysis hinges on the refined investigation of the eigenvalue problem corresponding to the Stokes operator Aε with Navier friction boundary condition.     

Boundary Shape Control of the Navier-Stokes Equations and Applications

In this paper, the geometrical design for the blade's surface in an impeller or for the profile of an aircraft, is modeled from the mathematical point of view by a boundary shape control problem for the Navier-Stokes equations. The objective function is the sum of a global dissipative function and the power of the fluid. The control variables are the geometry of the boundary and the state equations are the Navier-Stokes equations. The Euler-Lagrange equations of the optimal control problem are derived, which are an elliptic boundary value system of fourth order, coupled with the Navier-Stokes equations. The authors also prove the existence of the solution of the optimal control problem, the existence of the solution of the Navier-Stokes equations with mixed boundary conditions, the weak continuity of the solution of the Navier-Stokes equations with respect to the geometry shape of the blade's surface and the existence of solutions of the equations for the Gateaux derivative of the solution of the Navier-Stokes equations with respect to the geometry of the boundary.     

Couette-Taylor流的谱Galerkin逼近

利用谱方法对轴对称的旋转圆柱间的Couette-Taulor流进行数值模拟．首先给出Navier-Stokes方程流函数形式,利用Couette流把边界条件齐次化．其次给出Stokes算子的特征函数的解析表达式,证明其正交性,并对特征值进行估计．最后利用Stokes算子的特征函数作为逼近子空间的基函数,给出谱Galerkin逼近方程的表达式．证明了Navier-Stokes方程非奇异解的谱Galerkin逼近的存在性、唯一性和收敛性,给出了解谱Galerkin逼近的误差估计,并展示了数值计算结果．     

A High-Order Kernel-Free Boundary Integral Method for Incompressible Flow Equations in Two Space Dimensions

This paper presents a fourth-order kernel-free boundary integral method for the time-dependent, incompressible Stokes and Navier-Stokes equations defined on irregular bounded domains. By the stream function-vorticity formulation, the incompressible flow equations are interpreted as vorticity evolution equations. Time discretization methods for the evolution equations lead to a modified Helmholtz equation for the vorticity, or alternatively, a modified biharmonic equation for the stream function with two clamped boundary conditions. The resulting fourth-order elliptic boundary value problem is solved by a fourth-order kernel-free boundary integral method, with which integrals in the reformulated boundary integral equation are evaluated by solving corresponding equivalent interface problems, regardless of the exact expression of the involved Green's function. To solve the unsteady Stokes equations, a four-stage composite backward differential formula of the same order accuracy is employed for time integration. For the Navier-Stokes equations, a three-stage third-order semi-implicit Runge-Kutta method is utilized to guarantee the global numerical solution has at least third-order convergence rate. Numerical results for the unsteady Stokes equations and the Navier-Stokes equations are presented to validate efficiency and accuracy of the proposed method.     

On the rotating Navier-Stokes equations with mixed boundary conditions

The stationary and nonstationary rotating Navier-Stokes equations with mixed boundary conditions are investigated in this paper. The existence and uniqueness of the solutions are obtained by the Galerkin approximation method. Next, θ-scheme of operator splitting algorithm is applied to rotating Navier-Stokes equations and two subproblems are derived. Finally, the computational algorithms for these subproblems are provided.