Rotor wake capture improvement based on high-order spatially accurate schemes and chimera grids

A high-order upwind scheme has been developed to capture the vortex wake of a helicopter rotor in the hover based on chimera grids. In this paper, an improved fifth-order weighted essentially non-oscillatory (WENO) scheme is adopted to interpolate the higher-order left and right states across a cell interface with the Roe Riemann solver updating inviscid flux, and is compared with the monotone upwind scheme for scalar conservation laws (MUSCL). For profitably capturing the wake and enforcing the period boun...     

Numerical simulation for 2D double-diffusive convection (DDC) in rectangular enclosures based on a high resolution upwind compact streamfunction model I: numerical method and code validation

Applied Mathematics and Mechanics - A high resolution upwind compact streamfunction numerical algorithm for two-dimensional (2D) double-diffusive convection (DDC) is developed. The unsteady...     

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

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

EFFICIENT COMPUTATIONS USING UPWIND BIASED SCHEMES

A new multiblock unfactored implicit upwind scheme for inviscid two-dimensional flow calculations is presented. Spatial discretization is carried out by means of an upwind first-order method; an original extension to higher accuracy is also presented. The integration algorithm is constructed in a ‘δ’ form that provides a direct derivation of the scheme and leads to an efficient computational method. Fast solutions of the linear systems arising at each time step are obtained by means of the bi-conjugate gradient stabilized technique. The computational results for super/hypersonic steady state flows illustrate the efficiency and accuracy of the algorithm.     

Multigrid solutions of the Euler and Navier–Stokes equations on unstructured grids

A computationally efficient multigrid algorithm for upwind edge‐based finite element schemes is developed for the solution of the two‐dimensional Euler and Navier–Stokes equations on unstructured triangular grids. The basic smoother is based upon a Galerkin approximation employing an edge‐based formulation with the explicit addition of an upwind‐type local extremum diminishing (LED) method. An explicit time stepping method is used to advance the solution towards the steady state. Fully unstructured grids are employed to increase the flexibility of the proposed algorithm. A full approximation storage (FAS) algorithm is used as the basic multigrid acceleration procedure. Copyright © 1999 John Wiley & Sons, Ltd.     

Finite element analysis of the viscous flow in a vaned radial diffuser

A new finite element method was used to analyze an experimental model of a radial vaned diffuser. The new method includes a streamline upwind formulation for the advection terms in the governing equations. The streamline upwind significantly reduces numerical diffusion while maintaining the stability of the conventional upwind formulation. The new finite element method also incorporates an iterative equal-order, velocity-pressure solution method based on the well-known SIMPLER algorithm. The results of the analysis are compared to flow visualization studies of the experimental model. The flow separation point for the four blade diffuser was predicted to occur at 19, 6% of the blade length from the leading edge. The experimentally determined value was 23% of the blade length. For the eight blade diffuser model, separation was predicted to occur at 43% of the blade length from the leading edge, as compared to the experimentally observed value of 50% of the blade length. With this performance comparison, the proposed finite element method has been demonstrated to be reliable for predicting complex fluid flows.     

An upwind finite element scheme for high-Reynolds-number flows

A new upwind finite element scheme for the incompressible Navier-Stokes equations at high Reynolds number is presented. The idea of the upwind technique is based on the choice of upwind and downwind points. This scheme can approximate the convection term to third-order accuracy when these points are located at suitable positions. From the practical viewpoint of computation, the algorithm of the pressure Poisson equation procedure is adopted in the framework of the finite element method. Numerical results of flow problems in a cavity and past a circular cylinder show excellent dependence of the solutions on the Reynolds number. The influence of rounding errors causing Karman vortex shedding is also discussed in the latter problem.     

ALE迎风有限元法研究进展

综述了ALE描述方法及ALE描述下的运动学关系, 介绍了流体力学中运动边界的ALE追踪方法; 对在解决带自由液面复杂流动问题中与ALE技术密切相关的网格生成和自动更新技术给出了简要概述并讨论了3种不同的网格生成与自动更新方法, 讨论了用经典的有限元方法模拟带自由液面液体流动问题时非物理振荡产生的原因并分析了两种消除非物理振荡的途径, 重点讨论了迎风格式并详述了迎风有限元的发展, 评述了几种将ALE描述和迎风格式相结合的ALE迎风有限元法的最新进展, 并展望了该领域进一步研究的方向.      

An adaptive multigrid finite-volume scheme for incompressible Navier–Stokes equations

An algorithm for the solutions of the two-dimensional incompressible Navier–Stokes equations is presented. The algorithm can be used to compute both steady-state and time-dependent flow problems. It is based on an artificial compressibility method and uses higher-order upwind finite-volume techniques for the convective terms and a second-order finite-volume technique for the viscous terms. Three upwind schemes for discretizing convective terms are proposed here. An interesting result is that the solutions computed by one of them is not sensitive to the value of the artificial compressibility parameter. A second-order, two-step Runge–Kutta integration coupling with an implicit residual smoothing and with a multigrid method is used for achieving fast convergence for both steady- and unsteady-state problems. The numerical results agree well with experimental and other numerical data. A comparison with an analytically exact solution is performed to verify the space and time accuracy of the algorithm.     

Assessment of different reconstruction techniques for implementing the NVSF schemes on unstructured meshes

Three new far‐upwind reconstruction techniques, New‐Technique 1, 2, and 3, are proposed in this paper, which localize the normalized variable and space formulation (NVSF) schemes and facilitate the implementation of standard bounded high‐resolution differencing schemes on arbitrary unstructured meshes. By theoretical analysis, it is concluded that the three new techniques overcome two inherent drawbacks of the original technique found in the literature. Eleven classic high‐resolution NVSF schemes developed in the past decades are selected to evaluate performances of the three new techniques relative to the original technique. Under the circumstances of arbitrary unstructured meshes, stretched meshes, and uniform triangular meshes, for each NVSF scheme, the accuracies and convergence properties, when implementing the four aforementioned far‐upwind reconstruction techniques respectively, are assessed by the pure convection of several scalar profiles. The numerical results clearly show that New‐Technique‐2 leads to a better performance in terms of overall accuracy and convergence behavior for the 11 NVSF schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

High‐order upwind compact scheme and simulation of turbulent premixed V‐flame

A high‐order accurate upwind compact difference scheme with an optimal control coefficient is developed to track the flame front of a premixed V‐flame. In multi‐dimensional problems, dispersion effect appears in the form of anisotropy. By means of Fourier analysis of the operators, anisotropic effects of the upwind compact difference schemes are analysed. Based on a level set algorithm with the effect of exothermicity and baroclinicity, the flame front is tracked. The high‐order accurate upwind compact scheme is employed to approximate the level set equation. In order to suppress numerical oscillations, the group velocity control technique is used and the upwind compact difference scheme is combined with the random vortex method to simulate the turbulent premixed V‐flame. Distributions of velocities and flame brush thickness are obtained by this technique and found to be comparable with experimental measurement. Copyright © 2005 John Wiley & Sons, Ltd.     

UNIFIED COMPUTATION OF FLOW WITH COMPRESSIBLE AND INCOMPRESSIBLE FLUID BASED ON ROE''''S SCHEME

A unified numerical scheme for the solutions of the compressible and incompressible Navier-Stokes equations is investigated based on a time-derivative preconditioning algorithm. The primitive variables are pressure, velocities and temperature. The time integration scheme is used in conjunction with a finite volume discretization. The preconditioning is coupled with a high order implicit upwind scheme based on the definition of a Roe's type matrix. Computational capabilities are demonstrated through computations of high Mach number, middle Mach number, very low Mach number, and incompressible flow. It has also been demonstrated that the discontinuous surface in flow field can be captured for the implementation Roe's scheme.     

Fourier analysis of a class of upwind schemes in shallow water systems for gravity and Rossby waves

A Fourier analysis has been performed for a class of upwind finite volume schemes, including the study of phase speed, group velocity, damping and dispersion. In the first part, pure gravity waves are investigated. As expected, most upwind schemes lead to a significant damping, but they exhibit a better phase behavior than most centered schemes. In the second part, the Coriolis parameter is considered and the Rossby modes are studied. In this case, all selected upwind schemes lead to a severe damping. The numerical results are also compared with those obtained by using a slope limiter approach. It is concluded that most upwind schemes with or without slope limiters present poor results for an accurate calculation of the Rossby modes. Copyright © 2008 John Wiley & Sons, Ltd.     

Finite volume multigrid method of the planar contraction flow of a viscoelastic fluid

This paper reports on a numerical algorithm for the steady flow of viscoelastic fluid. The conservative and constitutive equations are solved using the finite volume method (FVM) with a hybrid scheme for the velocities and first‐order upwind approximation for the viscoelastic stress. A non‐uniform staggered grid system is used. The iterative SIMPLE algorithm is employed to relax the coupled momentum and continuity equations. The non‐linear algebraic equations over the flow domain are solved iteratively by the symmetrical coupled Gauss–Seidel (SCGS) method. In both, the full approximation storage (FAS) multigrid algorithm is used. An Oldroyd‐B fluid model was selected for the calculation. Results are reported for planar 4:1 abrupt contraction at various Weissenberg numbers. The solutions are found to be stable and smooth. The solutions show that at high Weissenberg number the domain must be long enough. The convergence of the method has been verified with grid refinement. All the calculations have been performed on a PC equipped with a Pentium III processor at 550 MHz. Copyright © 2001 John Wiley & Sons, Ltd.     

A note on two upwind strategies for RBF‐based grid‐free schemes to solve steady convection–diffusion equations

In this paper, two radial basis function (RBF)‐based local grid‐free upwind schemes have been discussed for convection–diffusion equations. The schemes have been validated over some convection–diffusion problems with sharp boundary layers. It is found that one of the upwind schemes realizes the boundary layers more accurately than the rest. Comparisons with the analytical solutions demonstrate that the local RBF grid‐free upwind schemes based on the exact velocity direction are stable and produce accurate results on domains discretized even with scattered distribution of nodal points. Copyright © 2009 John Wiley & Sons, Ltd.     

A study on the numerical stability of the two-fluid model near ill-posedness

The two-fluid model is widely used in studying gas–liquid flow inside pipelines because it can qualitatively predict the flow field at low computational cost. However, the two-fluid model becomes ill-posed when the slip velocity exceeds a critical value, and computations can be quite unstable before the flow reaches the ill-posed condition. In this work, computational stability of various convection schemes together with the Euler implicit method for the time derivatives in conjunction with the two-fluid model is analyzed. A pressure correction algorithm for the two-fluid model is carefully implemented to minimize its effect on numerical stability. von Neumann stability analysis shows that the central difference scheme is more accurate and more stable than the 1st-order upwind, 2nd-order upwind, and QUICK schemes. The 2nd-order upwind scheme is much more susceptible to instability than the 1st-order upwind scheme and is inaccurate for short waves. Excellent agreement is obtained between the predicted and computed growth rates of harmonic disturbances. The instability associated with the two-fluid model discretized system of equations is related to but quantitatively different from the instability associated with ill-posedness of the two-fluid model. When the computation becomes unstable due to the ill-posedness, the machine roundoff errors from a selected range of short wavelengths, which scale with the grid size, are amplified rapidly to render the computation of any targeted long wavelength variation useless. For the viscous two-fluid model with wall friction and interfacial drag, a small-amplitude long wavelength disturbance grows due to viscous Kelvin–Helmholtz instability without triggering the grid scale short waves when the system remains well posed. Under such a condition, central difference is found to be the most accurate discretization scheme among those investigated.     

Insights on a sign‐preserving numerical method for the advection–diffusion equation

In this paper we explore theoretically and numerically the application of the advection transport algorithm introduced by Smolarkiewicz to the one‐dimensional unsteady advection–diffusion equation. The scheme consists of a sequence of upwind iterations, where the initial iteration is the first‐order accurate upwind scheme, while the subsequent iterations are designed to compensate for the truncation error of the preceding step. Two versions of the method are discussed. One, the classical version of the method, regards the second‐order terms of the truncation error and the other considers additionally the third‐order terms. Stability and convergence are discussed and the theoretical considerations are illustrated through numerical tests. The numerical tests will also indicate in which situations are advantageous to use the numerical methods presented. Copyright © 2008 John Wiley & Sons, Ltd.     

VISCOELASTIC MODELLING OF ENTRANCE FLOW USING MULTIMODE LEONOV MODEL

A simulation of planar 2D flow of a viscoelastic fluid employing the Leonov constitutive equation has been presented. Triangular finite elements with lower-order interpolations have been employed for velocity and pressure as well as the extra stress tensor arising from the constitutive equation. A generalized Lesaint–Raviart method has been used for an upwind discretization of the material derivative of the extra stress tensor in the constitutive equation. The upwind scheme has been further strengthened in our code by also introducing a non-consistent streamline upwind Petrov–Galerkin method to modify the weighting function of the material derivative term in the variational form of the constitutive equation. A variational equation for configurational incompressibility of the Leonov model has also been satisfied explicitly. The corresponding software has been used to simulate planar 2D entrance flow for a 4:1 abrupt contraction up to a Deborah number of 670 (Weissenberg number of 6·71) for a rubber compound using a three-mode Leonov model. The predicted entrance loss is found to be in good agreement with experimental results from the literature. Corresponding comparisons for a commercial-grade polystyrene, however, indicate that the predicted entrance loss is low by a factor of about four, indicating a need for further investigation. © 1997 by John Wiley & Sons, Ltd.     

Augmented upwind numerical schemes for the groundwater transport advection‐dispersion equation with local operators

When solute transport is advection‐dominated, the advection‐dispersion equation approximates to a hyperbolic‐type partial differential equation, and finite difference and finite element numerical approximation methods become prone to artificial oscillations. The upwind scheme serves to correct these responses to produce a more realistic solution. The upwind scheme is reviewed and then applied to the advection‐dispersion equation with local operators for the first‐order upwinding numerical approximation scheme. The traditional explicit and implicit schemes, as well as the Crank‐Nicolson scheme, are developed and analyzed for numerical stability to form a comparison base. Two new numerical approximation schemes are then proposed, namely, upwind–Crank‐Nicolson scheme, where only for the advection term is applied, and weighted upwind‐downwind scheme. These newly developed schemes are analyzed for numerical stability and compared to the traditional schemes. It was found that an upwind–Crank‐Nicolson scheme is appropriate if the Crank‐Nicolson scheme is only applied to the advection term of the advection‐dispersion equation. Furthermore, the proposed explicit weighted upwind‐downwind finite difference numerical scheme is an improvement on the traditional explicit first‐order upwind scheme, whereas the implicit weighted first‐order upwind‐downwind finite difference numerical scheme is stable under all assumptions when the appropriate weighting factor (θ) is assigned.     

A high‐order discontinuous Galerkin method for extension problems

We present a novel technique for solving extension problems such as the extension velocity, by reformulating the problem into an elliptic differential equation. We introduce a novel discretization using an upwind flux without any additional stabilization. This leads to a triangular matrix structure, which can be solved using a marching algorithm and high‐order accuracy, even in the presence of singularities.