Least-squares finite-element method for shallow-water equations with source terms

Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.     

Vorticity–velocity formulation of the 3D Navier–Stokes equations in cylindrical co‐ordinates

A finite difference method is presented for solving the 3D Navier–Stokes equations in vorticity–velocity form. The method involves solving the vorticity transport equations in ‘curl‐form’ along with a set of Cauchy–Riemann type equations for the velocity. The equations are formulated in cylindrical co‐ordinates and discretized using a staggered grid arrangement. The discretized Cauchy–Riemann type equations are overdetermined and their solution is accomplished by employing a conjugate gradient method on the normal equations. The vorticity transport equations are solved in time using a semi‐implicit Crank–Nicolson/Adams–Bashforth scheme combined with a second‐order accurate spatial discretization scheme. Special emphasis is put on the treatment of the polar singularity. Numerical results of axisymmetric as well as non‐axisymmetric flows in a pipe and in a closed cylinder are presented. Comparison with measurements are carried out for the axisymmetric flow cases. Copyright © 2003 John Wiley & Sons, Ltd.     

Large‐eddy simulation with complex 2‐D geometries using a parallel finite‐element/spectral algorithm

A parallel stabilized finite‐element/spectral formulation is presented for incompressible large‐eddy simulation with complex 2‐D geometries. A unique discretization scheme is developed consisting of a streamline‐upwind Petrov–Galerkin/Pressure‐Stabilized Petrov–Galerkin (SUPG/PSPG) finite‐element discretization in the 2‐D plane with a collocated spectral/pseudospectral discretization in the out‐of‐plane direction. This formulation provides an efficient approach for solving 3‐D flows over arbitrary 2‐D geometries. Utilizing this discretization and through explicit temporal treatment of the non‐linear terms, the system of equations for each Fourier mode is decoupled within each time step. A novel parallelization approach is then taken, where the computational work is partitioned in Fourier space. A validation of the algorithm is presented via comparison of results for flow past a circular cylinder with published values for Re=195, 300, and 3900. Copyright © 2003 John Wiley & Sons, Ltd.     

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

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

A multigrid pseudo-spectral method for incompressible Navier–Stokes flows

A pseudo-spectral solver with multigrid acceleration for the numerical prediction of incompressible non-isothermal flows is presented. The spatial discretization is based on a Chebyshev collocation method on Gauss–Lobatto points and for the discretization in time the second-order backward differencing scheme (BDF2) is employed. The multigrid method is invoked at the level of algebraic system solving within a pressure-correction method. The approach combines the high accuracy of spectral methods with efficient solver properties of multigrid methods. The capabilities of the proposed scheme are illustrated by a buoyancy driven cavity flow as a standard benchmark case. To cite this article: K. Krastev, M. Schäfer, C. R. Mecanique 333 (2005).     

Three‐dimensional modeling of non‐hydrostatic free‐surface flows on unstructured grids

A three‐dimensional numerical model is presented for the simulation of unsteady non‐hydrostatic shallow water flows on unstructured grids using the finite volume method. The free surface variations are modeled by a characteristics‐based scheme, which simulates sub‐critical and super‐critical flows. Three‐dimensional velocity components are considered in a collocated arrangement with a σ‐coordinate system. A special treatment of the pressure term is developed to avoid the water surface oscillations. Convective and diffusive terms are approximated explicitly, and an implicit discretization is used for the pressure term to ensure exact mass conservation. The unstructured grid in the horizontal direction and the σ coordinate in the vertical direction facilitate the use of the model in complicated geometries. Solution of the non‐hydrostatic equations enables the model to simulate short‐period waves and vertically circulating flows. Copyright © 2015 John Wiley & Sons, Ltd.     

Time‐accurate intermediate boundary conditions for large eddy simulations based on projection methods

Projection methods are among the most adopted procedures for solving the Navier–Stokes equations system for incompressible flows. In order to simplify the numerical procedures, the pressure–velocity de‐coupling is often obtained by adopting a fractional time‐step method. In a specific formulation, suitable for the incompressible flows equations, it is based on a formal decomposition of the momentum equation, which is related to the Helmholtz–Hodge Decomposition theorem of a vector field in a finite domain. Owing to the continuity constraint also in large eddy simulation of turbulence, as happens for laminar solutions, the filtered pressure characterizes itself only as a Lagrange multiplier, not a thermodynamic state variable. The paper illustrates the implications of adopting such procedures when the decoupling is performed onto the filtered equations system. This task is particularly complicated by the discretization of the time integral of the sub‐grid scale tensor. A new proposal for developing time‐accurate and congruent intermediate boundary conditions is addressed. Several tests for periodic and non‐periodic channel flows are presented. This study follows and completes the previous ones reported in (Int. J. Numer. Methods Fluids 2003; 42, 43 ). Copyright © 2005 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.     

A unified approach to direct and inverse boundary layer solutions

A unified approach is presented for solving the two-dimensional incompressible boundary layer equations. Solutions are obtained for direct and inverse options using the same equation formulation by a simple interchange of boundary conditions. A modified form of the mechul function scheme obtains inverse solutions with specification of transformed wall shear, skin friction coefficient or displacement thickness distributions. Direct solutions may be obtained without altering the block tridiagonal structure of the system by simply requiring no corrections on the streamwise pressure gradient parameter. Fourth-order spline discretization approximates normal derivatives with two- and three-point backward differences approximating streamwise derivatives, yielding a fully implicit solution method. The resulting spline/finite difference equations are solved by Newton-Raphson iteration together with partial pivoting. The results of the study demonstrate the importance of proper linearization of all equations. The successful use of spline discretization is also tied to the use of strong two-point boundary conditions at the wall for cases involving reversed flow. Numerical solutions are presented for several non-similar flows and compared with published results.     

A critical evaluation of seven discretization schemes for convection–diffusion equations

A comparative study of seven discretization schemes for the equations describing convection-diffusion transport phenomena is presented. The (differencing) schemes considered are the conventional central- and upwind-difference schemes, together with the Leonard,¹ Leonard upwind¹ and Leonard super upwind difference¹ schemes. Also tested are the so called locally exact difference scheme² and the quadratic-upstream difference scheme.^3,4 In multidimensional problems errors arise from ‘false-diffusion’ and function approximations. It is asserted that false diffusion is essentially a multidimensional source of error. No mesh constraints are associated with errors in function approximation and discretization. Hence errors associated with discretization only may be investigated via one-dimensional problems. Thus, although the above schemes have been tested for one- and two-dimensional flows with sources, only the former are presented here. For 1D flows, the Leonard super upwind difference scheme and the locally exact scheme are shown to be far superior in accuracy to the others at all Peclet numbers and for most source distributions, for the test cases considered. Furthermore, the latter is shown to be considerably cheaper in computational terms than the former. The stability of the schemes and their CPU time requirements are also discussed.     

Finite element implementation of two‐equation and algebraic stress turbulence models for steady incompressible flows

The main purpose of this paper is to describe a finite element formulation for solving the equations for k and ε of the classical k–ε turbulence model, or any other two‐equation model. The finite element discretization is based on the SUPG method together with a discontinuity capturing technique to deal with sharp internal and boundary layers. The iterative strategy consists of several nested loops, the outermost being the linearization of the Navier–Stokes equations. The basic k–ε model is used for the implementation of an algebraic stress model that is able to account for the effects of rotation. Some numerical examples are presented in order to show the performance of the proposed scheme for simulating directly steady flows, without the need of reaching the steady state through a transient evolution. Copyright © 1999 John Wiley & Sons, Ltd.     

Development of an Artificial Compressibility Methodology with Implicit LU-SGS Method

A numerical method has been developed to solve the steady and unsteady incompressible Navier-Stokes equations in a two-dimensional, curvilinear coordinate system. The solution procedure is based on the method of artificial compressibility and uses a third-order flux-difference splitting upwind differencing scheme for convective terms and second-order center difference for viscous terms. A time-accurate scheme for unsteady incompressible flows is achieved by using an implicit real time discretization and a dual-time approach, which introduces pseudo-unsteady terms into both the mass conservation equation and momentum equations. An efficient fully implicit algorithm LU-SGS, which was originally derived for the compressible Eulur and Navier-Stokes equations by Jameson and Toon [1], is developed for the pseudo-compressibility formulation of the two dimensional incompressible Navier-Stokes equations for both steady and unsteady flows. A variety of computed results are presented to validate the present scheme. Numerical solutions for steady flow in a square lid-driven cavity and over a backward facing step and for unsteady flow in a square driven cavity with an oscillating lid and in a circular tube with a smooth expansion are respectively presented and compared with experimental data or other numerical results.     

Implicit preconditioned high‐order compact scheme for the simulation of the three‐dimensional incompressible Navier–Stokes equations with pseudo‐compressibility method

This paper combines the pseudo‐compressibility procedure, the preconditioning technique for accelerating the time marching for stiff hyperbolic equations, and high‐order accurate central compact scheme to establish the code for efficiently and accurately solving incompressible flows numerically based on the finite difference discretization. The spatial scheme consists of the sixth‐order compact scheme and 10th‐order numerical filter operator for guaranteeing computational stability. The preconditioned pseudo‐compressible Navier–Stokes equations are marched temporally using the implicit lower–upper symmetric Gauss–Seidel time integration method, and the time accuracy is improved by the dual‐time step method for the unsteady problems. The efficiency and reliability of the present procedure are demonstrated by applications to Taylor decaying vortices phenomena, double periodic shear layer rolling‐up problem, laminar flow over a flat plate, low Reynolds number unsteady flow around a circular cylinder at Re = 200, high Reynolds number turbulence flow past the S809 airfoil, and the three‐dimensional flows through two 90°curved ducts of square and circular cross sections, respectively. It is found that the numerical results of the present algorithm are in good agreement with theoretical solutions or experimental data. Copyright © 2011 John Wiley & Sons, Ltd.     

TRANSIENT SIMULATION OF INTERNAL SEPARATED FLOWS USING AN INTELLIGENT HIGHER-ORDER SPATIAL DISCRETIZATION SCHEME

This paper summarizes the method-of-lines (MOL) solution of the Navier–Stokes equations for an impulsively started incompressible laminar flow in a circular pipe with a sudden expansion. An intelligent higher-order spatial discretization scheme, which chooses upwind or downwind discretization in a zone-of-dependence manner when flow reversal occurs, was developed for separated flows. Stability characteristics of a linear advective–diffusive equation were examined to depict the necessity of such a scheme in the case of flow reversals. The proposed code was applied to predict the time development of an impulsively started flow in a pipe with a sudden expansion. Predictions were found to show the expected trends for both unsteady and steady states. This paper demonstrates the ease with which the Navier–Stokes equations can be solved in an accurate manner using sophisticated numerical algorithms for the solution of ordinary differential equations (ODEs). Solutions of the Navier–Stokes equations in primitive variables formulation by using the MOL and intelligent higher-order spatial discretization scheme are not available to date. © 1997 by John Wiley & Sons, Ltd.     

Computation of non-equilibrium hypersonic flow with artificially upstream flux vector splitting (AUFS) schemes

The AUFS scheme has been presented for solving the Euler equations [Sun, M., Takayama, K., 2003. An artificially upstream flux vector splitting scheme for the Euler equations. Journal of Computational Physics, 189, 305–329]. An extension of this high resolution scheme-based on upwind numerical methods has been developed to calculate a two-dimensional hypersonic viscous flowfield in thermochemical non-equilibrium. The time-dependent Navier–Stokes governing equations are computed by using a multi-block finite volume technique on a structured mesh. The convective fluxes at the interfaces are evaluated using a flux vector splitting (FVS) method with a second-order reconstruction of the interface values and the viscous terms are discretised by second-order central differences. A better evaluation of aerodynamic parameters are obtained with this AUFS scheme and they are also compared to those obtained by previous works. The freestream flow conditions of these computations correspond to high-enthalpy flows with a Mach number range between 6.4 and 25.9. The obtained numerical results indicate that the AUFS scheme is accurate, robust, and efficient for the calculation of hypersonic flow.     

An implicit meshless method for application in computational fluid dynamics

An implicit meshless scheme is developed for solving the Euler equations, as well as the laminar and Reynolds‐averaged Navier–Stokes equations. Spatial derivatives are approximated using a least squares method on clouds of points. The system of equations is linearised, and solved implicitly using approximate, analytical Jacobian matrices and a preconditioned Krylov subspace iterative method. The details of the spatial discretisation, linear solver and construction of the Jacobian matrix are discussed; and results that demonstrate the performance of the scheme are presented for steady and unsteady two dimensional fluid flows. Copyright © 2012 John Wiley & Sons, Ltd.     

Interval Arithmetic and Static Interval Finite Element Method

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｅａｎａｌｙｓｉｓａｎｄｄｅｓｉｇｎｏｆｓｔｒｕｃｔｕｒｅｓ,ｓｏｍｅｕｎａｖｏｉｄａｂｌｅｕｎｃｅｒｔａｉｎｔｉｅｓ ,ｓｕｃｈａｓｔｈａｔｏｆｍａｔｅｒｉａｌａｎｄｇｅｏｍｅｔｒｉｃａｌｐｒｏｐｅｒｔｉｅｓ,ｌｏａｄｓ ,ａｎｄｓｏｏｎ ,ｓｈｏｕｌｄｂｅｒｅａｓｏｎａｂｌｙｔａｋｅｎｉｎｔｏａｃｃｏｕｎｔ.Ｉｎｔｈｅｐａｓｔｄｅｃａｄｅｓ,ｔｈｅｓｅｕｎｃｅｒｔａｉｎｔｉｅｓｗｅｒｅｍｏｓｔｌｙｔｒｅａｔｅｄｗｉｔｈｐｒｏｂａｂｉｌｉｔｙｔｈｅｏｒｙｏｒｒａｎｄｏｍｐ…     

A PRESSURE-CORRECTION METHOD FOR THE SOLUTION OF INCOMPRESSIBLE VISCOUS FLOWS ON UNSTRUCTURED GRIDS

An unstructured grid, finite volume method is presented for the solution of two-dimensional viscous, incompressible flow. The method is based on the pressure-correction concept implemented on a semi-staggered grid. The computational procedure can handle cells of arbitrary shape, although solutions presented herein have been obtained only with meshes of triangular and quadrilateral cells. The discretization of the momentum equations is effected on dual cells surrounding the vertices of primary cells, while the pressure-correction equation applies to the primary-cell centroids and represents the conservation of mass across the primary cells. A special interpolation scheme s used to suppress pressure and velocity oscillations in cases where the semi-staggered arrangement does not ensure a sufficiently strong coupling between pressure and velocity to avoid such oscillations. Computational results presented for several viscous flows are shown to be in good agreement with analytical and experimental data reported in the open literature.     

亚、跨、超音速及不可压流动的数值分析方法的研究

为了对亚、跨、超音速及不可压无粘流动进行数值模拟,将LU－SGS方法与预处理方法结合,给出了PLU－SGS方法。方程离散基于有限体积法,采用高阶精度AUSMPW格式。方程求解采用了特征边界条件。通过典型算例的数值试验对比分析,表明PLU－SGS方法可以有效地对亚、跨、超音速及不可压流动进行数值模拟,并具有较高的计算精度和收敛速度。     

可压缩多介质流动的间断有限元方法

采用间断有限元方法、LS方法和通量装配技术相结合,建立了一种计算可压缩多介质流动的有效 方法。计算中以光滑Heavside函数构造流体比热比和重新初始化方程中的符号距离函数,并采用通量装配 技术抑制界面附近的非物理振荡。为解决可压缩多介质流动提供一种新的手段。