对称破坏分歧点的分裂迭代方法

构造了求解对称破坏分歧点的扩充系统，采用分裂分块迭代方法逼近对称破坏分歧点，并对2.Box Brusselator反应模型进行了数值模拟．     

Couette-Taylor流的谱Galerkin逼近

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

二维N-S方程的Fourier非线性Galerkin方法*

本文对周期边界条件Navier-Stokes方程,证明了其Fourier非线性Galerkin逼近解的存在唯一性,同时给出了逼近解的误差估计。     

求解Navier-Stokes方程的非退化转向点的扩充系统的一步牛顿迭代法

设（λ0,u0）是Navier-Stokes方程的非退化转向点，其中λ0=1/Re0,Re0为雷诺数，娄N充分大时，在（λ0，u0）的某个邻域内，谱Galerkin逼近方程存在唯一解（λ0^N,u0^N),(λ0^N,u0^N)为谱Galerkin逼近方程的非退化转向点，且有误差估计|λ0^N-λ0| λN 1^-1/2||u0^N-u0|| |u0^N-u0|≤cλN 1^-1,其中λi,i=2,…为Stokes算子的特征值，求解（λ0^N,u0^N）等价于求解某个扩充系统的非奇异解（u0^N,φ0^N,λ0^N)。我们证明，如果选取初值为（u0^m,φ0^m,λ0^m），其中m为与N相比很小的正整数，则这个扩充系统的线性化方程的解（λm^N,um^N)即可达到（λ0^N,u0^N）的精度。     

加罚N-S方程的有限元非线性Galerkin方法

非线性Galerkin方法是对耗散型非线性发展方程的一种数值解法,其空间变量不象一般Galerkin方法那样在线性空间上离散,而是在非线性流形上离散,所得逼近解在时间变量增大时可以更快地逼近其精确解.精细的理论分析可见[1],[2]等,在有限元逼近基础上将此方法应用到Navier-Stokes方程上的工作可参见[3],[4],这些文章主要针对速度与压力同时求解的混合元情形做了讨论.本文在[4]的基础上对加罚Navier-Stokes方程的一种非线性Galerkin方法的半离散和全离散有限元逼近格式分别进行了误差估     

加罚N-S方程的有限元非线性Galerkin方法

非线性Galerkin方法是对耗散型非线性发展方程的一种数值解法,其空间变量不象一般Galerkin方法那样在线性空间上离散,而是在非线性流形上离散,所得逼近解在时间变量增大时可以更快地逼近其精确解.精细的理论分析可见[1],[2]等,在有限元逼近基础上将此方法应用到Navier-Stokes方程上的工作可参见[3],[4],这些文章主要针对速度与压力同时求解的混合元情形做了讨论.本文在[4]的基础上对加罚Navier-Stokes方程的一种非线性Galerkin方法的半离散和全离散有限元逼近格式分别进行了误差估     

混合边界条件下定常Navier-Stokes方程解的正则性

本文考虑混合边界条件下,三维定常Navier-Stokes方程.利用Galerkin逼近方法和差商方法,证明了弱解和强解的存在性.     

Hopf分歧问题及其数值分析

本文讨论Banach空间中一般发展方程的Hopf分歧问题及其数值逼近,文中说明了连续问题及其逼近形式的Hopf分歧解的存在性,并给出近似解的收敛性和误差估计,推广了C.Bernadi的结论,针对非定常不可压Navier-Stokes方程的Hopf分歧解运用谱方法作了逼近。     

N阶分歧问题的数值分析

构造了求解一类非退化分歧点及相关参数的扩充系统，给出了拟牛顿迭代法并证明了收敛性．     

非定常 N-S 方程的非线性 Galerkin 算法及其误差估计

本文给出了二维非定常N-S方程的三种数值格式,其中空间变量用谱非线性Galerkin算法进行离散,时间变量用有限差分离散,并研究了这些格式数值解的逼近精度．最后,给出了部分数值计算结果．     

Stability and convergence of optimum spectral non‐linear Galerkin methods

Our objective in this article is to present some numerical schemes for the approximation of the 2‐D Navier–Stokes equations with periodic boundary conditions, and to study the stability and convergence of the schemes. Spatial discretization can be performed by either the spectral Galerkin method or the optimum spectral non‐linear Galerkin method; time discretization is done by the Euler scheme and a two‐step scheme. Our results show that under the same convergence rate the optimum spectral non‐linear Galerkin method is superior to the usual Galerkin methods. Finally, numerical example is provided and supports our results. Copyright © 2001 John Wiley & Sons, Ltd.     

The coupling system of Navier–Stokes equations and elastic Navier–Lame equations in a blood vessel

In this paper, the blood flow problem is considered in a blood vessel, and a coupling system of Navier–Stokes equations and linear elastic equations, Navier–Lame equations, in a cylinder with cylindrical elastic shell is given as the governing equations of the problem. We provide two finite element models to simulating the three-dimensional Navier–Stokes equations in the cylinder while the asymptotic expansion method is used to solving the linearly elastic shell equations. Specifically, in order to discrete the Navier–Stokes equations, the dimensional splitting strategy is constructed under the cylinder coordinate system. The spectral method is adopted along the rotation direction while the finite element method is used along the other directions. By using the above strategy, we get a series of two-dimensional-three-components (2D-3C) fluid problems. By introduce the S-coordinate system in E³ and employ the thickness of blood vessel wall as the expanding parameter, the asymptotic expansion method can be established to approximate the solution of the 3D elastic problem. The interface contact conditions can be treated exactly based on the knowledge of tensor analysis. Finally, numerical test shows that our method is reasonable.     

Semi-Galerkin approximation and strong solutions to the equations of the nonhomogeneous asymmetric fluids

This paper analyzes an initial/boundary value problem for a system of equations modelling the nonstationary flow of a nonhomogeneous incompressible asymmetric (polar) fluid. Under conditions similar to those usually imposed to the nonhomogeneous 3D Navier–Stokes equations, using a spectral semi-Galerkin method, we prove the existence of a local in time strong solution. We also prove the uniqueness of the strong solution and some global existence results. Several estimates for the solutions and their approximations are given. These can be used to find useful error bounds of the Galerkin approximations.     

Long- time Convergence of Numerical Approximations for Semilinear Parabolic Equations(Ⅰ )

§ 0 .Introduction　 We consider the numerical approximations of the dynamical systems governed bysemilinear parabolic equations,which are discretized by Galerkin and nonlinear Galerkinmethods in space,and by Runge-Kutta method in time.The numerical approximationson a finite time interval have already been widely studied(see[1 ]— [5] ) .We areconcerned with the long-time convergence and error estimates.This article is composedof three parts. In part ,we provide an abstract framework. In§…     

Effect of constriction height on flow separation in a two-dimensional channel

We studied numerically the effect of the constriction height on viscous flow separation past a two-dimensional channel with locally symmetric constrictions. A numerically stable scheme in primitive variables (velocity and pressure) for the solution of two-dimensional incompressible time-dependent Navier–Stokes equations is employed using finite-difference approximation in staggered grid. The wall shear stresses at different heights of the constriction are computed and presented graphically. It is noticed that the maximum stress and the length of the recirculating region associated with two shear layers of the constriction increase with the increase of the area reduction of the constriction. The critical Reynolds number for symmetry breaking bifurcation for the 50%, 60% and 70% area reduction are obtained numerically. The flow field separates after the symmetry breaking bifurcation and the symmetry of the flow depends on the Reynolds number and the height of the constriction.     

A convergent FEM-DG method for the compressible Navier–Stokes equations

This paper presents a new numerical method for the compressible Navier–Stokes equations governing the flow of an ideal isentropic gas. To approximate the continuity equation, the method utilizes a discontinuous Galerkin discretization on piecewise constants and a basic upwind flux. For the momentum equation, the method is a new combined discontinuous Galerkin and finite element method approximating the velocity in the Crouzeix–Raviart finite element space. While the diffusion operator is discretized in a standard fashion, the convection and time-derivative are discretized using discontinuous Galerkin on the element average velocity and a Lax–Friedrich type flux. Our main result is convergence of the method to a global weak solution as discretization parameters go to zero. The convergence analysis constitutes a numerical version of the existence analysis of Lions and Feireisl.     

Fast spectral solutions of the double-gyre problem in a turbulent flow regime

Several semi-analytical models are considered for a double-gyre problem in a turbulent flow regime for which a reference fully numerical eddy-resolving solution is obtained. The semi-analytical models correspond to solving the depth-averaged Navier–Stokes equations using the spectral Galerkin approach. The robustness of the linear and Smagorinsky eddy-viscosity models for turbulent diffusion approximation is investigated. To capture essential properties of the double-gyre configuration, such as the integral kinetic energy, the integral angular momentum, and the jet mean-flow distribution, an improved semi-analytical model is suggested that is inspired by the idea of scale decomposition between the jet and the surrounding flow.     

定常不可压Navier—Stokes方程迎风格式

本文利用齐次定解条件对定常不可压Navier—Stokes方程的非线性项进行处理,给出了相应的一种迎风Galerkin有限元算法;针对这种迎风Galerkin有限元算法,在迎风参数满足一定条件下,利用其三项式具有的一些很好性质,更简单地证明了该问题解的存在唯一性。     

An adaptive Lagrange‐Galerkin method for the shallow‐water equations on the sphere

The shallow‐water equations on the sphere are compressible Navier‐Stokes equations that describe hydrostatic atmospheric dynamics. An adaptive Lagrange‐Galerkin method is applied on the spherical domain to develop a simplified global atmospheric model. Computational results show the experimental convergence of this numerical method.