An analysis of the time integration algorithms for the finite element solutions of incompressible Navier–Stokes equations based on a stabilised formulation

This work is concerned with the analysis of time integration procedures for the stabilised finite element formulation of unsteady incompressible fluid flows governed by the Navier–Stokes equations. The stabilisation technique is combined with several different implicit time integration procedures including both finite difference and finite element schemes. Particular attention is given to the generalised-α method and the linear discontinuous in time finite element scheme. The time integration schemes are first applied to two model problems, represented by a first order differential equation in time and the one dimensional advection–diffusion equation, and subjected to a detailed mathematical analysis based on the Fourier series expansion. In order to establish the accuracy and efficiency of the time integration schemes for the Navier–Stokes equations, a detailed computational study is performed of two standard numerical examples: unsteady flow around a cylinder and flow across a backward facing step. It is concluded that the semi-discrete generalised-α method provides a viable alternative to the more sophisticated and expensive space–time methods for simulations of unsteady flows of incompressible fluids governed by the Navier–Stokes equations.     

Higher Order Finite Volume Central Schemes for Multi-dimensional Hyperbolic Problems

Different ways of implementing dimension-by-dimension CWENO reconstruction are discussed and the most efficient method is applied to develop a fourth order accurate finite volume central scheme for multi-dimensional hyperbolic problems. Fourth order accuracy and shock capturing nature of the scheme are demonstrated in various nonlinear multi-dimensional problems. In order to show the overall performance of the present central scheme numerical errors and non-oscillatory behavior are compared with existing multi-dimensional CWENO based central schemes for various multi-dimensional problems. Moreover, the benefits of the present fourth order central scheme over third order implementation are shown by comparing the numerical dissipation and computational cost between the two.     

Dispersive and Dissipative Properties of Discontinuous Galerkin Finite Element Methods for the Second-Order Wave Equation

Discontinuous Galerkin finite element methods (DGFEM) offer certain advantages over standard continuous finite element methods when applied to the spatial discretisation of the acoustic wave equation. For instance, the mass matrix has a block diagonal structure which, used in conjunction with an explicit time stepping scheme, gives an extremely economical scheme for time domain simulation. This feature is ubiquitous and extends to other time-dependent wave problems such as Maxwell’s equations. An important consideration in computational wave propagation is the dispersive and dissipative properties of the discretisation scheme in comparison with those of the original system. We investigate these properties for two popular DGFEM schemes: the interior penalty discontinuous Galerkin finite element method applied to the second-order wave equation and a more general family of schemes applied to the corresponding first order system. We show how the analysis of the multi-dimensional case may be reduced to consideration of one-dimensional problems. We derive the dispersion error for various schemes and conjecture on the generalisation to higher order approximation in space     

Multilevel methods for discontinuous Galerkin FEM on locally refined meshes

We review shortly the analytical results on multilevel preconditioning for DG methods for elliptic problems. An algorithm is presented which does not loose efficiency under local refinement. It uses subspace smoothing and does not require integration of additional matrices. This algorithm is applied to the interior penalty method and the local discontinuous Galerkin method. A feasible way of implementing the scheme is presented.     

二维有限元网格自动自适应生成

为使有限元法能更有效地应用于日益复杂的科学和工程计算,提出一种自动自适应网格生成方法．该方法在结合推进波前法和法向偏移法的优势的基础上产生,改进了节点的生成方式,融入了简单灵活的密度控制方式和合理准确的边界离散方式．剖分实例显示,该方法能有效地生成高质量网格,自动化程度高,适应能力强;它已成功地应用于复杂海域下波浪与结构物相互作用的有限元计算,并取得了很好的效果．     

The soft way method and the velocity formulation

In this paper the new approach to dynamic contact problems is described. The velocity formulation was assumed and a new time integration scheme was elaborated. The space-time finite element method used in derivation enables control of the accuracy (order of the error) and stability. Methods for the solution of contact problems were discussed. A discretized approach, prepared for large displacements and large rotations, enabled real engineering problems to be solved in a relatively short time.     

A Discontinuous Galerkin Method for Linear Symmetric Hyperbolic Systems in Inhomogeneous Media

The Discontinuous Galerkin (DG) method provides a powerful tool for approximating hyperbolic problems. Here we derive a new space-time DG method for linear time dependent hyperbolic problems written as a symmetric system (including the wave equation and Maxwell’s equations). The main features of the scheme are that it can handle inhomogeneous media, and can be time-stepped by solving a sequence of small linear systems resulting from applying the method on small collections of space-time elements. We show that the method is stable provided the space-time grid is appropriately constructed (this corresponds to the usual time-step restriction for explicit methods, but applied locally) and give an error analysis of the scheme. We also provide some simple numerical tests of the algorithm applied to the wave equation in two space dimensions (plus time).This revised version was published online in July 2005 with corrected volume and issue numbers.     

Fluid–shell structure interaction analysis by coupled particle and finite element method

This paper presents a coupled particle and finite element method for fluid–shell structure interaction analysis. The Moving Particle Semi-Implicit (MPS) method is used to analyze fluid flow and the MITC4 shell element is used in the FEM analysis of the structure. This paper considers partitioned coupling between the fluid and structural solvers. In order to satisfy compatibility in the employed partitioned coupling scheme, the Neumann–Dirichlet condition is applied to both the fluid and the structure. A symplectic time integration scheme is used to preserve energy when analyzing the shell structure. If the frequencies of the shell analysis are much higher than those of the MPS fluid solver, its time integration scheme is sub-cycled. When the presented coupling scheme was applied to simulate the sloshing phenomenon in an elastic thin shell structure, fluid fragmentation and large structural deformations were observed.     

非光滑区域上椭圆型特征值问题的间断有限元方法应用

本文针对非光滑区域上椭圆特征值特征值问题利用间断有限元方法（DG）近似．利用大量的数值算例发现,DG方法对非光滑区域（凹角,裂缝等问题）上Laplace特征值问题的近似比协调有限元、非协调元（如C—R元）,甚至比有限元校正格式有着更好的效果．     

A Generic Stabilization Approach for Higher Order Discontinuous Galerkin Methods for Convection Dominated Problems

In this paper we present a stabilized Discontinuous Galerkin (DG) method for hyperbolic and convection dominated problems. The presented scheme can be used in several space dimension and with a wide range of grid types. The stabilization method preserves the locality of the DG method and therefore allows to apply the same parallelization techniques used for the underlying DG method. As an example problem we consider the Euler equations of gas dynamics for an ideal gas. We demonstrate the stability and accuracy of our method through the detailed study of several test cases in two space dimension on both unstructured and cartesian grids. We show that our stabilization approach preserves the advantages of the DG method in regions where stabilization is not necessary. Furthermore, we give an outlook to adaptive and parallel calculations in 3d.     

Locally Implicit Time Integration Strategies in a Discontinuous Galerkin Method for Maxwell’s Equations

An attractive feature of discontinuous Galerkin (DG) spatial discretization is the possibility of using locally refined space grids to handle geometrical details. However, locally refined meshes lead to severe stability constraints on explicit integration methods to numerically solve a time-dependent partial differential equation. If the region of refinement is small relative to the computational domain, the time step size restriction can be overcome by blending an implicit and an explicit scheme where only the solution variables living at fine elements are treated implicitly. The downside of this approach is having to solve a linear system per time step. But due to the assumed small region of refinement relative to the computational domain, the overhead will also be small while the solution can be advanced in time with step sizes determined by the coarse elements. In this paper, we present two locally implicit time integration methods for solving the time-domain Maxwell equations spatially discretized with a DG method. Numerical experiments for two-dimensional problems illustrate the theory and the usefulness of the implicit–explicit approaches in presence of local refinements.     

一种自适应鲁棒控制方法及其闭环稳定性分析

研究了含有未建模动态的慢时变系统的自适应镇定问题.考虑的对象具有非最小相 位、含未建模动态和大范围时变参数等不良特性,且存在未知但有界外部扰动.这类对象很难 用时不变鲁棒控制器或传统自适应控制器进行镇定.利用l1优化设计结合参数估计的投影算 法,提出了一种自适应鲁棒控制策略.基于l1优化设计的连续性和投影算法的收敛性,证明了 这种控制策略能够持续适应慢时变对象并且保持闭环系统一致稳定性.鲁棒性分析表明这种 控制策略具有良好的鲁棒镇定性.     

Adaptive Multilevel Correction Method for Finite Element Approximations of Elliptic Optimal Control Problems

In this paper we propose an adaptive multilevel correction scheme to solve optimal control problems discretized with finite element method. Different from the classical adaptive finite element method (AFEM for short) applied to optimal control which requires the solution of the optimization problem on new finite element space after each mesh refinement, with our approach we only need to solve two linear boundary value problems on current refined mesh and an optimization problem on a very low dimensional space. The linear boundary value problems can be solved with well-established multigrid method designed for elliptic equation and the optimization problems are of small scale corresponding to the space built with the coarsest space plus two enriched bases. Our approach can achieve the similar accuracy with standard AFEM but greatly reduces the computational cost. Numerical experiments demonstrate the efficiency of our proposed algorithm.     

基于结构网格的高超声速流动大规模并行数值模拟研究

本文采用MPI消息传递模式自主开发出适用于高超声速流动数值模拟的并行计算软件,该软件以三维Navier-Stokes方程为基本控制方程来求解层流问题,应用基于结构网格的有限体积法对计算域进行离散,采用AUSMPW+格式求解对流通量,利用MUSCL插值方法获得高阶精度,时间格式上采用LU-SGS方法进行时间迭代以加快求解定常流动的收敛过程。在高性能计算机上针对不同高超声速流动进行大规模并行计算的结果表明,所开发的CFD并行计算软件具有较高的并行计算效率,为高超声速飞行器气动力/热的准确预测提供了高效工具。     

An adaptive method for the determination of the order of element for acquiring the desired accuracy in 3D elastostatic FEM analysis

An adaptive method for the determination of the order of element (or element order) was developed for the finite element analysis of 3D elastostatic problems. Here the order of element means the order of polynomial function, which interpolates the displacement distribution in the element. This method was based on acquiring the desired accuracy for each finite element. From the numerical experiments, the relationship ξ=k(1/p)^β was deduced, where ξ is the error of the result of the finite element analysis relative to the exact value, p is the order of element, and k and β are constants. Applying this relationship to the two results of computations with different orders of element, the order of element for the third computation was deduced. A computer program using this adaptive determination method for the order of element was developed and applied to several 3D elastostatic problems of various shapes. The usefulness of the method was illustrated by these application results.     

Sensitivity analysis of shock wave Burgers’ equation via a novel algorithm based on scale-3 Haar wavelets

In this paper, a novel technique is being formulated for the numerical solutions of Shock wave Burgers' equations for planar and non-planar geometry. It is well known that Burgers' equation is sensitive to the perturbations in the diffusion term. Thus we use robustness of wavelets generated by dilation and translation of Haar wavelets on third scale to capture the sensitivity information. The present approach is an improved form of the scale-2 Haar wavelet method. The scheme is based on the forward finite difference scheme for time integration, scale-3 Haar wavelets for space integration and the nonlinearity has been tackled via quasilinearzation technique. Through scale-3 Haar wavelet analysis once the wavelet coefficient is calculated then we can compute the solutions at near the perturbation point. The computation cost of the present scheme is negligible. The proposed method is tested on six test problems to check its computational efficiency where the convergence analysis of scale-3 Haar wavelet method is the proof of our computational arguments.     

Primal Discontinuous Galerkin Methods for Time-Dependent Coupled Surface and Subsurface Flow

This paper introduces and analyzes a numerical method based on discontinuous finite element methods for solving the two-dimensional coupled problem of time-dependent incompressible Navier-Stokes equations with the Darcy equations through Beaver-Joseph-Saffman’s condition on the interface. The proposed method employs Crank-Nicolson discretization in time (which requires one step of a first order scheme namely backward Euler) and primal DG method in space. With the correct assumption on the first time step optimal error estimates are obtained that are high order in space and second order in time.     

An adaptive finite element scheme for transient problems in CFD

An adaptive finite element scheme for transient problems is presented. The classic h-enrichment / coarsening is employed in conjunction with a triangular finite element discretization in two dimensions. A mesh change is performed every n timesteps, depending on the Courant number employed and the number of ‘protective layers’ added ahead of the refined region. In order to simplify the refinement/ coarsening logic and to be as fast as possible, only one level of refinement/coarsening is allowed per mesh change. A high degree of vectorizability has been achieved on the CRAY XMP 12 at NRL. Several examples involving shock-shock interactions and the impact of shocks on structures demonstrate the performance of the method, indicating that considerable savings in CPU time and storage can be realized even for strongly unsteady flows.