间断Galerkin有限元和有限体积混合计算方法研究

通过局部坐标变换而建立的非正交单元间断Galerkin(DG)有限元计算方法计算精度高， 但计算量大、内存需求大；而非结构网格有限体积方法虽然准确计算热流的问题目 前还没有完全解决，但其具有计算速度快和内存需求小的优点. 该研究是将有 限元和有限体积方法的优点结合，发展有限元和有限体积的混合方法. 在物面 附近黏性占主导作用的区域内采用有限元方法进行计算，在远离物面的区域采用快速的有限 体积方法进行计算，在有限元和有限体积方法结合处要保证通量守恒. 通过算例说明有 限元和有限体积混合方法既能保证黏性区域的热流计算精度和流场结构的分辨率，又能 降低内存需求和提高计算效率，使有限元方法应用于复杂外形(实际工程问题)的计 算成为可能.     

THE APPLICATION OF A FINITE DIFFERENCE/FINITE VOLUME ALGORITHM TO SOLVE MAXWELL'S EQUATIONS IN THE TIME DOMAIN

Abstract

A finite volume/finite difference method based on Ni's multigrid formulation is introduced for the solution of Maxwell's equations. The scheme is presented for the cases of transverse magnetic scattering from two-dimensional circular and square cylinders, as well as from NACA 0012 airfoil. The codes are validated against the traditional Method of Moments, which is analogous to a panel method in CFD. The circular cylinder scattering is compared to the analytical series solution for better understanding how the roles of numerical dispersion and dissipation errors affect the solution. The reflecting boundary conditions are modeled by the idea of inducing fields inside the conductor and a method of modeling the singularities that arise at a sharp corner is presented. Absorbing boundary conditions are modeled by integrating along the characteristic compatibility equations in the direction of the outgoing wave.     

Whitham—Broer—Kaup浅水波方程新的多孤子解

Whitham-Broer-Kaup（简记WBK）方程具有重要的意义,至今人们只给出了它的单孤子解,本利用齐次平衡法并借助数学给出它新的多孤子解,并作为特例得到一类变式Boussinesq方程的多孤子解。     

COMPARISON OF H AND P FINITE ELEMENT APPROXIMATIONS OF THE SHALLOW WATER EQUATIONS

A p-type finite element scheme is introduced for the three-dimensional shallow water equations with a harmonic expansion in time. The wave continuity equation formulation is used which decouples the problem into a Helmholtz equation for surface elevation and a momentum equation for horizontal velocity. An exploration of the applicability of p methods to this form of the shallow water problem is presented, with a consideration of the problem of continuity errors. The convergence rates and relative computational efficiency between h- and p- type methods are compared with the use of three test cases representing various degrees of difficulty. A channel test case establishes convergence rates, a continental shelf test case examines a problem with accuracy difficulties at the shelf break, and a field-scale test case examines problems with highly irregular grids. For the irregular grids, adaptive h combined with uniform p refinement was necessary to retain high convergence rates. © 1997 John Wiley & Sons, Ltd.     

A HYBRID SPECTRAL/FINITE VOLUME METHOD FOR VISCOUS FLOWS

Abstract

This paper presents a hybrid spectral/finite volume method for steady-state compressible viscous flows. The method is evaluated for accuracy via test cases for various Mach numbers. The domain is divided into a viscous region and an inviscid region. The viscous region uses the full Navier-Stokes equations, while the inviscid region employs the Euler equations. A high order Chebyshev collocation spectral method is developed for the viscous region to resolve boundary layers. This method avoids the dense grids needed by finite-volume methods to resolve the viscous areas. A low order finite-volume method based on a Lax-Wendroff type scheme is employed for the inviscid region. A special interface formulation is developed for coupling the spectral with the finite-volume method. Comparisons with analytic results as well as convergence histories are presented.     

AN ADAPTIVE UNSTRUCTURED TRI-TREE ITERATIVE SOLVER FOR MIXED FINITE ELEMENT FORMULATION OF THE STOKES EQUATIONS

An iterative adaptive equation solver for solving the implicit Stokes equations simultaneously with tri-tree grid generation is developed. The tri-tree grid generator builds a hierarchical grid structure which is mapped to a finite element grid at each hierarchical level. For each hierarchical finite element grid the Stokes equations are solved. The approximate solution at each level is projected onto the next finer grid and used as a start vector for the iterative equation solver at the finer level. When the finest grid is reached, the equation solver is iterated until a tolerated solution is reached. In order to reduce the overall work, the element matrices are integrated analytically beforehand. The efficiency and behaviour of the present adaptive method are compared with those of the previously developed iterative equation solver which is preconditioned by incomplete LU factorization with coupled node fill-in. The efficiency of the incomplete coupled node fill-in preconditioner is shown to be largely dependent on the global node numbering. The preconditioner is therefore tested for the natural node ordering of the tri-tree grid generator and for different ways of sorting the nodes.     

AN ENTROPY VARIABLE FORMULATION AND APPLICATIONS FOR THE TWO-DIMENSIONAL SHALLOW WATER EQUATIONS

A new symmetric formulation of the two-dimensional shallow water equations and a streamline upwind Petrov–Galerkin (SUPG) scheme are developed and tested. The symmetric formulation is constructed by means of a transformation of dependent variables derived from the relation for the total energy of the water column. This symmetric form is well suited to the SUPG approach as seen in analogous treatments of gas dynamics problems based on entropy variables. Particulars related to the construction of the upwind test functions and an appropriate discontinuity-capturing operator are included. A formal extension to the viscous, dissipative problem and a stability analysis are also presented. Numerical results for shallow water flow in a channel with (a) a step transition, (b) a curved wall transition and (c) a straight wall transition are compared with experimental and other computational results from the literature.     

ON BREAKING WAVES AND WAVE-CURRENT INTERACTION IN SHALLOW WATER: A 2DH FINITE ELEMENT MODEL

A two-dimensional (horizontal plane) coastal and estuarine region model, capable of predicting the combined effects of gravity surface shallow- water waves (shoaling, refraction, diffraction, reflection and breaking), and steady currents, is described and numerical results are compared with those obtained experimentally. Two series of observations within a wave flume and a combined wave-current facility were developed. In the first case, the wave was generated via a hinged paddle located within a deepened section at one end of the channel, as, in the second case, the wave propagating with or against the current was generated by a plunger-type wavemaker; the re-circulating current was introduced via one passing tank connected to a centrifugal pump. Several comparisons for a number of 1D situations and one 2D horizontal plane case are presented.     

变深度浅水域中非定常船波

以Green—Naghdi(G—N)方程为基础，采用波动方程／有限元法计算船舶经过变深度浅水域时非定常波浪特性．把运动船舶对水面的扰动作为移动压强直接加在Green-Naghdi方程里，以描述运动船体和水面的相互作用．以Series60 CB＝0．6船为算例，给出自由面坡高，波浪阻力在船舶经过一个水下凸包时变化规律，并与浅水方程的结果进行了比较．计算结果表明，当船舶经过凸包时，波浪阻力先增加，后减少，并逐渐趋于正常．同时发现，当船速小于临界速度时(Fr＝√gh＜1．0)，G—N方程给出的船后尾波比浅水方程的结果明显，波浪阻力也比浅水方程的结果有所提高，频率散射必须考虑．当船速大于临界速度时(Fr＝√gh＞1．0)，G—N方程的计算结果与浅水方程差别不大，频率散射的影响可以忽略．     

MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTATION (Ⅰ)-THE CONTINUOUS-TIME CASE

An initial-boundary value problem for shallow equation system consisting of water dynamics equations, silt transport equation, the equation of bottom topography change, and of some boundary and initial conditions is studied, the existence of its generalized solution and semidiscrete mixed finite element (MFE) solution was discussed, and the error estimates of the semidiscrete MFE solution was derived. The error estimates are optimal.     

固体力学应力分析中的有限域法

本文给出了用有限域法进行了和学应力分析的一般原理，并与目前广泛采用的有限元不做了比较。文中的研究表明：有限域法与有限元法极为相似，前者使用单位位移加权，后者从虚位移原理出发；两种方法的实施过程也有共同点，它们都进行网络部分及逼近。     

一维孔隙率浅水方程及其数值离散

通过孔隙率方法来描述挡水物对过水能力的影响建立了一维孔隙率浅水方程. 采用有限体积方法和Roe格式的近似Riemann解建立了孔隙率浅水方程的离散模式. 对底坡和孔隙率源项采用特性方向分解的方法进行处理,使模型精确满足C(Conservative)特性,增加了模型的稳定性. 通过算例模拟证明了模型可以对河道中的挡水物作用进行模拟,且计算结果表明模型具有和谐、稳定、分辨率高等优点.      

BOUNDARY CONDITIONS IN SHALLOW WATER MODELS—AN ALTERNATIVE IMPLEMENTATION FOR FINITE ELEMENT CODES

Finite element solution of the shallow water wave equations has found increasing use by researchers and practitioners in the modelling of oceans and coastal areas. Wave equation models, most of which use equal-orderC⁰ interpolants for both the velocity and the surface elevation, do not introduce spurious oscillation modes, hence avoiding the need for artificial or numerical damping. An important question for both primitive equation and wave equation models is the interpretation of boundary conditions. Analysis of the characteristics of the governing equations shows that for most geophysical flows a single condition at each boundary is sufficient, yet there is not a consensus in the literature as to what that boundary condition must be or how it should be implemented in a finite element code. Traditionally (partly because of limited data), surface elevation is specified at open ocean boundaries while the normal flux is specified as zero at land boundaries. In most finite element wave equation models both of these boundary conditions are implemented as essential conditions. Our recent work focuses on alternative ways to numerically implement normal flow boundary conditions with an eye towards improving the mass-conserving properties of wave equation models. A unique finite element formulation using generalized functions demonstrates that boundary conditions should be implemented by treating normal fluxes as natural conditions with the flux interpreted as external to the computational domain. Results from extensive numerical experiments show that the scheme does conserve mass for all parameter values. Furthermore, convergence studies demonstrate that the algorithm is consistent, as residual errors at the boundary diminish as the grid is refined.     

对流扩散方程的摄动有限体积(PFV)方法及讨论

在有限体积(FV)方法的重构近似中,引入数值摄动处理,即把界面数值通量摄动展开成网格间距的幂级数,并利用积分方程自身的性质求出幂级数的系数,同时获得高精度迎风和中心型摄动有限体积(PFV)格式.对标量输运方程给出积分近似为二阶、重构近似为二、三和四阶迎风和中心型PFV格式,这些PFV格式的结构形式及使用基点数与一阶迎风格式完全一致,迎风PFV格式满足对流有界准则;二阶和四阶中心PFV格式对网格Peclet数的任意值均为正型格式,比常用的二阶中心格式优越.用一维标量输运和方腔流动算例说明PFV格式的优良性能,并把PFV方法与性质相近的摄动有限差分(PFD)方法及相关的高精度方法作了对比分析.     

A NEW FLUX-CONSERVING NEWTON SCHEME FOR THE NAVIER-STOKES EQUATIONS

SUMMARY

A new numerical method is developed for the two-dimensional, steady Navier-Stokes equations. Using local polynomial expansions to represent the discrete primitive variables on each cell, we construct a scheme which has the following properties: First, the local discrete primitive variables are functional solutions of both the integral and differential forms of the Navier-Stokes equations. Second, fluxes are balanced across cell interfaces using exact functional expressions (to the order of accuracy of the local expansions). No interpolation, flux models, or flux limiters are required. Third, local and global conservation of mass, momentum, and energy are explicitly provided for. Finally, the discrete primitive variables and their derivatives are treated in a unified and consistent manner. All are treated as unknowns to be solved together for simulating the local and global flux conservation.

A general third-order formulation for the steady, compressible Navier-Stokes equations is presented. As a special case, the formulation is applied to incompressible flow, and a Newton's method scheme is developed for the solution of laminar channel flow. H is shown that, at Reynolds numbers of 100, 1000, and 2000, the developing channel flow boundary layer can be accurately resolved using as few as six to ten cells per channel width.     

用格子Boltzmann方程模拟浅水波问题

提出了用格子Ｂｏｌｔｚｍａｎｎ方程（ＬＢＥ）模拟浅水波问题的方法．通过无粘气体运动方程与浅水波方程的比较，确定了ＬＢＥ方法中平衡态的形式，使宏观方程与浅水波方程一致．计算了二维浅水波的一个问题，数值结果与精确解做了比较．     

A NON–LINEAR ADAPTED TRI–TREE MULTIGRID SOLVER FOR FINITE ELEMENT FORMULATION OF THE NAVIER–STOKES EQUATIONS

An iterative adaptive equation multigrid solver for solving the implicit Navier–Stokes equations simultaneously with tri-tree grid generation is developed. The tri-tree grid generator builds a hierarchical grid structur e which is mapped to a finite element grid at each hierarchical level. For each hierarchical finite element multigrid the Navier–Stokes equations are solved approximately. The solution at each level is projected onto the next finer grid and used as a start vector for the iterative equation solver at the finer level. When the finest grid is reached, the equation solver is iterated until a tolerated solution is reached. The iterative multigrid equation solver is preconditioned by incomplete LU factorization with coupled node fill-in. The non-linear Navier–Stokes equations are linearized by both the Newton method and grid adaption. The efficiency and behaviour of the present adaptive method are compared with those of the previously developed iterative equation solver which is preconditioned by incomplete LU factorization with coupled node fill-in.     

轴对称高速碰撞问题的拉格朗日无结构网格有限体积法的并行计算

给出计算轴对称高速碰撞问题的拉格朗日无结构三角形网格有限体积法的并行格式，并给出以小巨型机ＡｌｉａｎｔＦＸ／４０为目标计算机的算例数值模拟结果和效率分析     

近海水域海洋流体动力学方程的数值模拟

近海水域海洋流体动力学方程数值模拟的研究,最早可追溯到1919年Defant的工作,50年代末开始得到蓬勃发展。60年代和70年代主要发展二维全流模式。现在对于二维全流模式数值模拟的研究应该说是相当完善的,不仅有系统的理论,提出了许多有效的方法,而且广泛地应用于诸如潮汐环流、风暴潮等的工程和预报的实际计算。三维模式是70年代开始发展的,近10年来发展很快。三维模式比较起二维模式来,其数值模拟自然要困难许多。这引起了研究者广泛的兴趣。现在三维模式的数值模拟,虽然提出的方法很多,也已有了不少成果,但研究工作仍是方兴未艾。      

SINGLE-LAYER FINITE ELEMENT APPROXIMATION OF 3D NAVIER–STOKES EQUATIONS WITH FREE BOUNDARIES FOR SHALLOW WATER FLOWS

In order to simulate flows in the shallow water limit, the full incompressible Navier–Stokes equations with free boundaries are solved using a single layer of finite elements. This implies a polynomial approximation of the velocity profile in the vertical direction, which in turn distorts the wave speed. This fact is verified by numerical results: the wave speed depends on the vertical discretization. When at least two layers of finite elements are used, the boundary layer at the bottom can be simulated and the correct solution for the shallow water limit is recovered. Then this algorithm is applied to the prediction of Tsunami event.