Two-level stabilized finite element method for Stokes eigenvalue problem

A two-level stabilized finite element method for the Stokes eigenvalue problem based on the local Gauss integration is considered.This method involves solving a Stokes eigenvalue problem on a coarse mesh with mesh size H and a Stokes problem on a fine mesh with mesh size h = O(H 2),which can still maintain the asymptotically optimal accuracy.It provides an approximate solution with the convergence rate of the same order as the usual stabilized finite element solution,which involves solving a Stokes eigenvalue problem on a fine mesh with mesh size h.Hence,the two-level stabilized finite element method can save a large amount of computational time.Moreover,numerical tests confirm the theoretical results of the present method.     

Parallel finite element computation of incompressible magnetohydrodynamics based on three iterations

Based on local algorithms, some parallel finite element(FE) iterative methods for stationary incompressible magnetohydrodynamics(MHD) are presented. These approaches are on account of two-grid skill include two major phases: find the FE solution by solving the nonlinear system on a globally coarse mesh to seize the low frequency component of the solution, and then locally solve linearized residual subproblems by one of three iterations(Stokes-type, Newton, and Oseen-type) on subdomains with fine...     

不可压Navier-Stokes方程基于误差估算的网格自适应解法

首先导出了广义Stokes方程Petrov—Galerkin有限元数值解的当地事后误差估算公式;以非连续二阶鼓包（bump）函数空间为速度、压强误差的近似空间,该估算基于求解当地单元上的广义Stokes问题。然后,证明了误差估算值与精确误差之间的等价性。最后,将误差估算方法应用于Navier—Stokes环境,以进行不可压粘流计算中的网格自适应处理。数值实验中成功地捕获了多强度物理现象,验证了本文所发展的方法。     

Robust high-order semi-implicit backward difference formulas for Navier-Stokes equations

Taylor-Hood finite elements provide a robust numerical discretization of Navier-Stokes equations (NSEs) with arbitrary high order of accuracy in space. To match the accuracy of the lowest degree Taylor-Hood element, we propose a very efficient time-stepping methods for unsteady flows, which are based on high-order semi-implicit backward difference formulas (SBDF) and the inclusion of grad -div term in the NSE. To mitigate the impact on the numerical accuracy (in time) of the extrapolation of the nonlinear term in SBDF, several variants of nonlinear extrapolation formulas are investigated. The first approach is based on an extrapolation of the nonlinear advection term itself. The second formula uses the extrapolation of the velocity prior to the evaluation of the nonlinear advection term as a whole. The third variant is constructed such that it produces similar error on both velocity and pressure to that with fully implicit backward difference formulas methods at a given order of accuracy. This can be achieved by fixing one-order higher than usually done in the extrapolation formula for the nonlinear advection term, while keeping the same extrapolation formula for the time derivative. The resulting truncation errors (in time) between these formulas are identified using Taylor expansions. These truncation error formulas are shown to properly represent the error seen in numerical tests using a 2D manufactured solution. Lastly, we show the robustness of the proposed semi-implicit methods by solving test cases with high Reynolds numbers using one of the nonlinear extrapolation formulas, namely, the 2D flow past circular cylinder at Re=300 and Re = 1000 and the 2D lid-driven cavity at Re = 50 000 and Re = 100 000. Our numerical solutions are found to be in a good agreement with those obtained in the literature, both qualitatively and quantitatively.     

A postprocessing mixed finite element method for the Navier–Stokes equations

A fully discrete postprocessing mixed finite element scheme is considered for solving the time-dependent Navier–Stokes equations. In the PP method, we only consider a non-linear equation in the coarse-level subspace and a linear problem in the fine-level subspace. The analysis shows that the PP scheme can reach the same accuracy as the standard Galerkin method with a very fine mesh size h by an appropriate choice of H. Numerical examples are provided that confirm both the theoretical analysis and the corresponding improvement in computational efficiency.     

TWO-GRID ERROR ESTIMATES FOR THE STREAM FUNCTION FORM OF NAVIER-STOKES EQUATIONS

ＩｎｔｒｏｄｕｃｔｉｏｎＷｅｃｏｎｓｉｄｅｒｔｗｏ＿ｇｒｉｄｍｅｔｈｏｄｆｏｒｔｈｅｓｔｒｅａｍｆｕｎｃｔｉｏｎｆｏｒｍｏｆｔｈｅｓｔａｔｉｏｎａｒｙＮａｖｉｅｒ＿Ｓｔｏｋｅｓｅｑｕａｔｉｏｎｓ.Ｔｈｅａｄｖａｎｔａｇｅｓｏｆｔｈｅｓｔｒｅａｍｆｕｎｃｔｉｏｎｆｏｒｍａｒｅｔｈａｔｔｈｅｉｎｃｏｍｐｒｅｓｓｉｂｉｌｉｔｙｃｏｎｄｉｔｉｏｎｉｓｓａｔｉｓｆｉｅｄａｕｔｏｍａｔｉｃａｌｌｙａｎｄｔｈｅｐｒｅｓｓｕｒｅｉｓｎｏｔｐｒｅｓｅｎｔｉｎｔｈｅｗｅａｋｆｏｒｍ .Ｔｈｅｍｅｔｈｏｄｉｓｂａｓｅｄ…     

High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics

In this paper, we present a class of high‐order accurate cell‐centered arbitrary Lagrangian–Eulerian (ALE) one‐step ADER weighted essentially non‐oscillatory (WENO) finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two‐dimensional unstructured triangular meshes. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space–time Galerkin predictor allows the schemes to be high order accurate also in time by using an element‐local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained either (i) as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu, or (ii) as a solution of multiple one‐dimensional half‐Riemann problems around a vertex, as suggested by Maire, or (iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional Harten–Lax–van Leer Riemann solver recently proposed by Balsara et al. Once the vertex velocity and thus the new node location have been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level tⁿ with the new ones at the next time level t^{n + 1}. If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unnecessary, as the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high‐order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic magnetohydrodynamic equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration. Copyright © 2014 John Wiley & Sons, Ltd.     

三维液体非线性晃动动力学特性的数值模拟

主要讨论圆筒形贮腔中三维液体非线性晃动问题,将任意的拉格朗日-欧拉（即ArbitraryLagrangian-Eulerian,简称ALE）运动学描述引入到Navies-Stokes方程中,在时间域上采用一种速度和压力的分步计算格式进行时间离散;在空间域上利用Galerkin加权余量法对系统方程进行数值离散;得到了数值计算粘性不可压液体非线性晃动的ALE分步有限元法的计算格式,推导了三维液体自由液面上结点法向矢量的数值计算方法,模拟了圆筒形贮腔（包括带圆环形隔板的圆筒形贮腔）中三维液体的非线性晃动;并得到了一些重要的非线必不知所云 性,通过数值模拟结果与实验结果的比较,证明实了本文方法的可靠性与有效性。     

Adaptive variational multiscale method for the Stokes equations

An adaptive variational multiscale method for the Stokes equations is presented in this paper. We solve the coarse scale problem on the coarse mesh and approximate the fine scale solution by solving a series of local residual equations defined on some local fine grids, which can be implemented in parallel. In addition, we also propose a reliable local a posteriori error estimator and construct an adaptive algorithm based on the corresponding a posterior error estimate. Finally, numerical examples are presented to verify the algorithm.Copyright © 2012 John Wiley & Sons, Ltd.     

Adaptive finite element method for an optimal control problem of Stokes flow with L2‐norm state constraint

An adaptive finite element approximation for an optimal control problem of the Stokes flow with an L²‐norm state constraint is proposed. To produce good adaptive meshes, the a posteriori error estimates are discussed. The equivalent residual‐type a posteriori error estimators of the H ¹‐error of state and L²‐error of control are given, which are suitable to carry out the adaptive multi‐mesh finite element approximation. Some numerical experiments are performed to illustrate the efficiency of the a posteriori estimators. Copyright © 2011 John Wiley & Sons, Ltd.     

Local and parallel finite element algorithms based on two-grid discretization for steady Navier-Stokes equations

Local and parallel finite element algorithms based on two-grid discretization for Navier-Stokes equations in two dimension are presented. Its basis is a coarse finite element space on the global domain and a fine finite element space on the subdomain. The local algorithm consists of finding a solution for a given nonlinear problem in the coarse finite element space and a solution for a linear problem in the fine finite element space, then droping the coarse solution of the region near the boundary. By overlapping domain decomposition, the parallel algorithms are obtained. This paper analyzes the error of these algorithms and gets some error estimates which are better than those of the standard finite element method. The numerical experiments are given too. By analyzing and comparing these results, it is shown that these algorithms are correct and high efficient.     

An hp‐adaptive unstructured finite volume solver for compressible flows

The idea of hp‐adaptation, which has originally been developed for compact schemes (such as finite element methods), suggests an adaptation scheme using a mixture of mesh refinement and order enrichment based on the smoothness of the solution to obtain an accurate solution efficiently. In this paper, we develop an hp‐adaptation framework for unstructured finite volume methods using residual‐based and adjoint‐based error indicators. For the residual‐based error indicator, we use a higher‐order discrete operator to estimate the truncation error, whereas this estimate is weighted by the solution of the discrete adjoint problem for an output of interest to form the adaptation indicator for adjoint‐based adaptations. We perform our adaptation by local subdivision of cells with nonconforming interfaces allowed and local reconstruction of higher‐order polynomials for solution approximations. We present our results for two‐dimensional compressible flow problems including subsonic inviscid, transonic inviscid, and subsonic laminar flow around the NACA 0012 airfoil and also turbulent flow over a flat plate. Our numerical results suggest the efficiency and accuracy advantages of adjoint‐based hp‐adaptations over uniform refinement and also over residual‐based adaptation for flows with and without singularities.     

An h4 accurate vorticity-velocity formulation for calculating flow past a cylinder

A method of solution for the two-dimensional Navier-Stokes equations for incompressible flow past a cylinder is given in which the euquation of continuity is solved by a step-by-step integration procedure at each stage of an interative process. Thus the formulation involves the solution of one first-order and one second-order equation for the velocity components, together with the vorticity transport equation. the equations are solved numerically by h⁴-accurate methods in the case of steady flow past a circular cylinder in the Reynolds number range 10–100. Results are in satisfactory agreement with recent h⁴-accurate calculations. An improved approximation to the boundary conditions at large distance is also considered.     

A simple high‐resolution advection scheme

A simple, robust, mass‐conserving numerical scheme for solving the linear advection equation is described. The scheme can estimate peak solution values accurately even in regions where spatial gradients are high. Such situations present a severe challenge to classical numerical algorithms. Attention is restricted to the case of pure advection in one and two dimensions since this is where past numerical problems have arisen. The authors' scheme is of the Godunov type and is second‐order in space and time. The required cell interface fluxes are obtained by MUSCL interpolation and the exact solution of a degenerate Riemann problem. Second‐order accuracy in time is achieved via a Runge–Kutta predictor–corrector sequence. The scheme is explicit and expressed in finite volume form for ease of implementation on a boundary‐conforming grid. Benchmark test problems in one and two dimensions are used to illustrate the high‐spatial accuracy of the method and its applicability to non‐uniform grids. Copyright © 2007 John Wiley & Sons, Ltd.     

极坐标系下Fourier-Legendre谱元方法与有限差分法数值扩散的比较

提出一种Fourier-Legendre谱元方法用于求解极坐标系下的Navier-Stokes方程,其中极点所在单元的径向采用Gauss-Radau积分点,避免了r=0处的1/r坐标奇异性。时间离散采用时间分裂法,引入数值同位素模型跟踪同位素的输运过程验证数值模拟的精度,分别利用谱元法和有限差分法的迎风差分格式求解匀速和加速坩埚旋转流动中的同位素方程。计算结果表明,有限差分法中的一阶迎风差分格式存在严重的数值假扩散,二阶迎风差分格式的数值结果较精确,增加节点可以有效地缓解数值扩散。然而,谱元法具有以较少节点得到高精度解的优势。     

GENERALIZED FINITE SPECTRAL METHOD FOR 1D BURGERS AND KDV EQUATIONS

A generalized finite spectral method is proposed. The method is of high-order accuracy. To attain high accuracy in time discretization, the fourth-order Adams-Bashforth-Moulton predictor and corrector scheme was used. To avoid numerical oscillations caused by the dispersion term in the KdV equation, two numerical techniques were introduced to improve the numerical stability. The Legendre, Chebyshev and Her-mite polynomials were used as the basis functions. The proposed numerical scheme is validated by applications to the Burgers equation (nonlinear convection- diffusion problem) and KdV equation (single solitary and 2-solitary wave problems), where analytical solutions are available for comparison. Numerical results agree very well with the corresponding analytical solutions in all cases.     

高阶谱元区域分解算法求解定常方腔驱动流

主要利用Jacobian-free的Newton-Krylov方法求解定常不可压缩Navier-Stokes方程,将基于高阶谱元法的区域分解Stokes算法的非定常时间推进步作为Newton迭代的预处理,回避了传统Newton方法Jacobian矩阵的显式装配,节省了程序内存,同时降低了Newton迭代线性系统的条件数,且没有非线性对流项的隐式求解,大大加快了收敛速度。对有分析解的Kovasznay流动的计算结果表明,本高阶谱元法在空间上有指数收敛的谱精度,且对定常解的Newton迭代是二次收敛的。本文模拟了二维方腔顶盖一致速度驱动流,同基准解符合得很好,表明本文方法是准确可靠的。本文还考虑了Re=800时方腔顶盖正弦速度驱动流,除得到已知的一个稳定对称解和一对稳定非对称解外,还获得了一对新的不稳定的非对称解。     

Moving mesh methods for Boussinesq equation

The Boussinesq equation is a challenging problem both analytically and numerically. Owing to the complex dynamic development of small scales and the rapid loss of solution regularity, the Boussinesq equation pushes any numerical strategy to the limit. With uniform meshes, the amount of computational time is too large to enable us to obtain useful numerical approximations. Therefore, developing effective and robust moving mesh methods for these problems becomes necessary. In this work, we develop an efficient moving mesh algorithm for solving the two‐dimensional Boussinesq equation. Our moving mesh algorithm is an extension of Tang and Tang (SIAM J. Numer. Anal. 2003; 41 :487–515) for hyperbolic conservation laws and Zhang and Tang (Commun. Pure Appl. Anal. 2002; 1 :57–73) for convection‐dominated equations. Several numerical fluxes (Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction (2nd edn). Springer: Berlin, 1999; WASCOM 99”: 10th Conference on Waves and Stability in Continuous Media, Porto Ercole, Italy, 1999; 257–264; High‐order Methods for Computational Physics. Springer: Berlin, 1999; 439–582; J. Sci. Comput. 1990; 5 :127–149; SIAM J. Numer. Anal. 2003; 41 :487–515; Commun. Pure Appl. Anal. 2002; 1 :57–73) are also discussed. Numerical results demonstrate the advantage of our moving mesh method in resolving the small structures. Copyright © 2009 John Wiley & Sons, Ltd.     

On the use of anisotropic a posteriori error estimators for the adaptative solution of 3D inviscid compressible flows

This paper describes the use of an a posteriori error estimator to control anisotropic mesh adaptation for computing inviscid compressible flows. The a posteriori error estimator and the coupling strategy with an anisotropic remesher are first introduced. The mesh adaptation is controlled by a single‐parameter tolerance (TOL) in regions where the solution is regular, whereas a condition on the minimal element size h_min is enforced across solution discontinuities. This h_min condition is justified on the basis of an asymptotic analysis. The efficiency of the approach is tested with a supersonic flow over an aircraft. The evolution of a mesh adaptation/flow solution loop is shown, together with the influence of the parameters TOL and h_min. We verify numerically that the effect of varying h_min is concordant with the conclusions of the asymptotic analysis, giving hints on the selection of h_min with respect to TOL. Finally, we check that the results obtained with the a posteriori error estimator are at least as accurate as those obtained with anisotropic a priori error estimators. All the results presented can be obtained using a standard desktop computer, showing the efficiency of these adaptative methods. Copyright © 2008 John Wiley & Sons, Ltd.     

Collocated discrete least‐squares (CDLS) meshless method: Error estimate and adaptive refinement

Meshless methods are new approaches for solving partial differential equations. The main characteristic of all these methods is that they do not require the traditional mesh to construct a numerical formulation. They require node generation instead of mesh generation. In other words, there is no pre‐specified connectivity or relationships among the nodes. This characteristic make these methods powerful. For example, an adaptive process which requires high computational effort in mesh‐dependent methods can be very economically solved with meshless methods. In this paper, a posteriori error estimate and adaptive refinement strategy is developed in conjunction with the collocated discrete least‐squares (CDLS) meshless method. For this, an error estimate is first developed for a CDLS meshless method. The proposed error estimator is shown to be naturally related to the least‐squares functional, providing a suitable posterior measure of the error in the solution. A mesh moving strategy is then used to displace the nodal points such that the errors are evenly distributed in the solution domain. Efficiency and effectiveness of the proposed error estimator and adaptive refinement process are tested against two hyperbolic benchmark problems, one with shocked and the other with low gradient smooth solutions. These experiments show that the proposed adaptive process is capable of producing stable and accurate results for the difficult problems considered. Copyright © 2007 John Wiley & Sons, Ltd.