解Stokes特征值问题的一种两水平稳定化有限元方法

基于局部Gauss积分,研究了解Stokes特征值问题的一种两水平稳定化有限元方法．该方法涉及在网格步长为H的粗网格上解一个Stokes特征值问题,在网格步长为h=O(H2)的细网格上解一个Stokes问题．这样使其能够仍旧保持最优的逼近精度,求得的解和一般的稳定化有限元解具有相同的收敛阶,即直接在网格步长为h的细网格上解一个Stokes特征值问题．因此,该方法能够节省大量的计算时间．数值试验验证了理论结果．     

Acceleration of two-grid stabilized mixed finite element method for the Stokes eigenvalue problem

This paper provides an accelerated two-grid stabilized mixed finite element scheme for the Stokes eigenvalue problem based on the pressure projection. With the scheme, the solution of the Stokes eigenvalue problem on a fine grid is reduced to the solution of the Stokes eigenvalue problem on a much coarser grid and the solution of a linear algebraic system on the fine grid. By solving a slightly different linear problem on the fine grid, the new algorithm significantly improves the theoretical error estimate which allows a much coarser mesh to achieve the same asymptotic convergence rate. Finally, numerical experiments are shown to verify the high efficiency and the theoretical results of the new method.     

Meshless solution of two-dimensional incompressible flow problems using the radial basis integral equation method

The two-dimensional incompressible fluid flow problems governed by the velocity–vorticity formulation of the Navier–Stokes equations were solved using the radial basis integral (RBIE) equation method. The RBIE is a meshless method based on the multi-domain boundary element method with overlapping subdomains. It solves at each node for the potential and its spatial derivatives. This feature of the RBIE is advantageous in solving the velocity–vorticity formulation of the Navier–Stokes equations since the calculated velocity gradients can be used to compute the vorticity that is prescribed as a boundary condition to the vorticity transport equation. The accuracy of the numerical solution was examined by solving the test problem with known analytical solution. Two benchmark problems, i.e. the lid driven cavity flow and the thermally driven cavity flow were also solved. The numerical results obtained using the RBIE showed very good agreement with the benchmark solutions.     

Integrability and Self-Similarity in Transient Stimulated Raman Scattering

Summary. The phenomenon of stimulated Raman scattering (SRS) can be described by three coupled PDEs which define the pump electric field, the Stokes electric field, and the material excitation as functions of distance and time. In the transient limit these equations are integrable, i.e., they admit a Lax pair formulation. Here we study this transient limit. The relevant physical problem can be formulated as an initial-boundary value (IBV) problem where both independent variables are on a finite domain. A general method for solving IBV problems for integrable equations has been introduced recently. Using this method we show that the solution of the equations describing the transient SRS can be obtained by solving a certain linear integral equation. It is interesting that this equation is identical to the linear integral equation characterizing the solution of an IBV problem of the sine-Gordon equation in light-cone coordinates. This integral equation can be solved uniquely in terms of the values of the pump and Stokes fields at the entry of the Raman cell. The asymptotic analysis of this solution reveals that the long-distance behavior of the system is dominated by the underlying self-similar solution which satisfies a particular case of the third Painlevé transcendent. This result is consistent with both numerical simulations and experimental observations. We also discuss briefly the effect of frequency mismatch between the pump and the Stokes electric fields. Received December 10, 1996; second revision received October 10, 1997; final revision received January 20, 1998     

二维矩形域内Stokes流问题的辛解析和数值方法

给出了一种新的解析求解二维矩形域中的Stokes流动问题的方法——辛体系方法（Hamilton体系方法）．在辛体系下,基本问题归结为本征值和本征解的问题．由于辛本征解之间存在辛正交共轭关系,问题的解和边界条件均可以由本征解描述和表示．利用辛本征解空间的完备性,建立一套封闭的求解问题方法．研究结果表明零本征值本征解描述了基本流动,而非零本征值本征解则表示问题的局部效应．数值结果给出了几种有代表性的流动情况,显示了该求解方法对求解许多问题的有效性．同时,这种方法也为研究其他问题提供了一条思路．     

A variational conjugate gradient method for determining the fluid velocity of a slow viscous flow

The problem considered is that of determining the fluid velocity for linear hydrostatics Stokes flow of slow viscous fluids from measured velocity and fluid stress force on a part of the boundary of a bounded domain. A variational conjugate gradient iterative procedure is proposed based on solving a series of mixed well-posed boundary value problems for the Stokes operator and its adjoint. In order to stabilize the Cauchy problem, the iterations are ceased according to an optimal order discrepancy principle stopping criterion. Numerical results obtained using the boundary element method confirm that the procedure produces a convergent and stable numerical solution.     

Canonical Integral Equations of Stokes Problem

The canonical boundary reduction, suggested by Feng Kang, also can be applied to the bidimensional steady Stokes problem. In this paper we first give the representation formula for the solution of the Stokes problem via two complex variable functions. Then by means of complex analysis and the Fourier analysis, we find the expressions of the Poisson integral formulas and the canonical integral equations in three typical domains. From these results the canonical boundary element method for solving the Stokes problem can be developed.     

Numerical analysis of the stochastic Stokes equations of Wick type

We propose a finite element method for the numerical solution of the stochastic Stokes equations of the Wick type. We give existence and uniqueness results for the continuous problem and its approximation. Optimal error estimates are derived and algorithmic aspects of the method are discussed. Our method will reduce the problem of solving stochastic Stokes equations to solving a set of deterministic ones. Moreover, one can reconstruct particular realizations of the solution directly from Wiener chaos expansions once the coefficients are available. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

求解定常不可压Stokes方程的两层罚函数方法

借助于两套有限元网格空间提出了一种求解定常不可压Stokes方程的两层罚函数方法.该方法只需要求解粗网格空间上的Stokes方程和细网格空间上的两个易于求解的罚参数方程（离散后的线性方程组具有相同的对称正定系数矩阵）.收敛性分析表明粗网格空间相对于细网格空间可以选择很小，并且罚参数的选取只与粗网格步长和问题的正则性有关.因此罚参数不必选择很小仍能够得到最优解.最后通过数值算例验证了上述理论结果，并且数值对比可知两层罚函数方法对于求解定常不可压Stokes方程具有很好的效果.     

Stokes方程的一个新的非协调四边形单元格式

对于Stokes方程给出了一个新的非协调四边形单元格式.新单元具有构造简单,自由度较少等优势.特别指出的是,该单元在矩形网格下,还是一个Locking-free元,可用于平面弹性问题.尽管该单元不含协调部分,其相容误差估计较困难,通过采用新的技巧和方法得到了最优误差估计.     

An iterative solution method for the Stokes problem with variable viscosity

An iterative method for efficient solution of the Stokes problem with a variable viscosity is considered. A preconditioner for the Shur complement is constructed taking into account the variable viscosity. The efficiency analysis is given. An application of the preconditioner for solving one problem of the mantle convection modeling is considered.     

Two‐level Newton iterative method for the 2D/3D steady Navier‐Stokes equations

A combination method of the Newton iteration and two‐level finite element algorithm is applied for solving numerically the steady Navier‐Stokes equations under the strong uniqueness condition. This algorithm is motivated by applying the m Newton iterations for solving the Navier‐Stokes problem on a coarse grid and computing the Stokes problem on a fine grid. Then, the uniform stability and convergence with respect to ν of the two‐level Newton iterative solution are analyzed for the large m and small H and h << H. Finally, some numerical tests are made to demonstrate the effectiveness of the method. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

A multilevel correction method for Stokes eigenvalue problems and its applications

In this paper, a new multilevel correction scheme is proposed to solve Stokes eigenvalue problems by the finite element method. This new scheme contains a series of correction steps, and the accuracy of eigenpair approximation can be improved after each step. In each correction step, we only need to solve a Stokes problem on the corresponding fine finite element space and a Stokes eigenvalue problem on the coarsest finite element space. This correction scheme can improve the efficiency of solving Stokes eigenvalue problems by the finite element method. As applications of this multilevel correction method, a multigrid method and an adaptive finite element technique are introduced for Stokes eigenvalue problems. Some numerical results are given to validate our schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

The two-grid interpolating element free Galerkin (TG-IEFG) method for solving Rosenau-regularized long wave (RRLW) equation with error analysis

The two-grid method is a technique to solve the linear system of algebraic equations for reducing the computational cost. In this study, the two-grid procedure has been combined with the EFG method for solving nonlinear partial differential equations. The two-grid FEM has been introduced in various forms. The well-known two-grid FEM is a three-step method that has been proposed by Bajpai and Nataraj (Comput. Math. Appl. 2014;68:2277–2291) that the new proposed scheme is an ecient procedure for solving important nonlinear partial differential equations such as Navier–Stokes equation. By applying shape functions of IMLS approximation in the EFG method, a new technique that is called interpolating EFG (IEFG) can be obtained. In the current investigation, we combine the two-grid algorithm with the IEFG method for solving the nonlinear Rosenau-regularized long-wave (RRLW) equation. In other hand, we demonstrate that solutions of steps 1, 2, and 3 exist and are unique and also we achieve an error estimate for them. Moreover, three test problems in one- and two-dimensional cases are given which support accuracy and efficiency of the proposed scheme.     

Numerical solution to a stokes interface problem

A numerical algorithm is described for solving the generalized eigenvalue problem arising in the study of the spectrum of a preconditioned operator in the pressure equation derived from a Stokes interface problem. The algorithm is implemented for two finite element schemes. It is tested for a problem with an analytical solution and is applied to spectrum computations in the case of a piecewise constant viscosity. A large number of numerical experiments are analyzed, and recommendations are given for solving the Stokes interface problem in practice.     

A simultaneous bus route design and frequency setting problem for Tin Shui Wai,Hong Kong

A bus network design problem for Tin Shui Wai, a suburban residential area in Hong Kong, is investigated, which considers the bus services from the origins inside this suburban area to the destinations in the urban areas. The problem aims to improve the existing bus services by reducing the number of transfers and the total travel time of the users. This has been achieved by the proposed integrated solution method which can solve the route design and frequency setting problems simultaneously. In the proposed solution method, a genetic algorithm, which tackles the route design problem, is hybridized with a neighborhood search heuristic, which tackles the frequency setting problem. A new solution representation scheme and specific genetic operators are developed so that the genetic algorithm can search all possible route structures, rather than selecting routes from the predefined set. To avoid premature convergence, a diversity control mechanism is incorporated in the solution method based on a new definition of hamming distance. To illustrate the robustness and quality of solutions obtained, computational experiments are performed based on 1000 perturbed demand matrices. The t-test results show that the design obtained by the proposed solution method is robust under demand uncertainty, and the design is better than both the current design and the design obtained by solving the route design problem and the frequency setting problem sequentially. Compared with the current bus network design, the proposed method can generate a design which can simultaneously reduce the number of transfers and total travel time at least by 20.9% and 22.7% respectively. Numerical studies are also performed to illustrate the effectiveness of the diversity control mechanism introduced and the effects of weights on the two objective values.     

A novel recurrent nonlinear neural network for solving quadratic programming problems

This paper presents a new neural network for solving quadratic programming problems. The new model has a simple form, furthermore it has a good convergence rate with a less number calculation operation than the old models. It converges very fast to exact solution of the dual problem and by substituting in a formulation, the optimal solution of the original problem is obtained. Neural network model with one of numerical method is solved. Finally, simple numerical examples are provided for more illustration.     

Mathematical modeling of the vertical magnetic dipole field in a two-dimensional nonhomogeneous medium

We consider a quasi-three-dimensional problem of remote marine sounding by a high-power stationary source located on land. A transition from the three-dimensional problem to a family of parametric two-dimensional problems is performed. The solution of the remote marine sounding problem is obtained with high accuracy after solving about 20 two-dimensional problems. The integral equations are solved by the modified integral current method, which has proved highly efficient for field computations inside a strongly conducting anomaly. The electric field amplitude is observed to increase with depth. The width of the coastal current channel is estimated by analyzing the vertical magnetic field component.     

A self-consistent field method applied to the dynamics of viscous suspensions

A method for the approximate solution of the problem of many bodies of spherical form in a viscous fluid is developed in the Stokes approximation. Using a purely hydrodynamic approach, based on the use of the concept of a self-consistent field, the classical boundary value problem is reduced to a formal procedure for solving a linear system of algebraic equations in the tensor coefficients, which occur in the solution obtained for the velocity field and pressure of the liquid. A procedure for the approximate solution of this system of equations is constructed for the case of dilute suspensions, when the ratio of the size of the dispersed particles to the characteristic distance between them is a small parameter. Finally, the initial boundary value problem is reduced to solving a recurrent system of equations, in which each subsequent approximation for all the required quantities depends solely on the previous approximations. The system of recurrent equations obtained can be solved analytically in any specified approximation with respect to a small parameter. It is shown that this system of equations contains in itself all possible physical formulations of the problems, and, within the frameworks of the mathematical procedure constructed, they are distinguished solely by a set of specified and required functions. The practical possibilities of the method are in no way limited by the number of dispersed particles in the fluid.     

A Dual-Primal FETI method for incompressible Stokes equations

In this paper, a dual-primal FETI method is developed for incompressible Stokes equations approximated by mixed finite elements with discontinuous pressures. The domain of the problem is decomposed into nonoverlapping subdomains, and the continuity of the velocity across the subdomain interface is enforced by introducing Lagrange multipliers. By a Schur complement procedure, the solution of an indefinite Stokes problem is reduced to solving a symmetric positive definite problem for the dual variables, i.e., the Lagrange multipliers. This dual problem is solved by the conjugate gradient method with a Dirichlet preconditioner. In each iteration step, both subdomain problems and a coarse level problem are solved by a direct method. It is proved that the condition number of this preconditioned dual problem is independent of the number of subdomains and bounded from above by the square of the product of the inverse of the inf-sup constant of the discrete problem and the logarithm of the number of unknowns in the individual subdomains. Numerical experiments demonstrate the scalability of this new method. This work is based on a doctoral dissertation completed at Courant Institute of Mathematical Sciences, New York University. This work was supported in part by the National Science Foundation under Grants NSF-CCR-9732208, and in part by the U.S. Department of Energy under contract DE-FG02-92ER25127.