An efficient low‐dissipative WENO filter scheme

This paper presents a simple and efficient procedure developed for tracing discontinuities in flow fields. Numerical experiments are carried out to test the new sensor coupled with the associated nonlinear WENO dissipation filter developed to suppress the numerical dissipation. The tests show that, for a problem containing shocks and vortices, the implementation of the new sensor and an optimized WENO scheme can obtain a stable and highly resolved solution. The numerical experiments demonstrated that the new filter scheme performs efficiently both in parallel and serial running for one‐dimensional inviscid flow problems. Direct numerical simulation of a Mach 5 turbulent boundary layer over a flat plate was carried out to demonstrate the applicability of the scheme to the DNS practices. Copyright © 2011 John Wiley & Sons, Ltd.     

数值耗散对压气机流动分离涡模拟的影响研究

分离涡模拟DES是压气机流动模拟中常用的高保真湍流模式。为了使DES准确解析湍流,数值耗散必须限制在合理范围内。然而,当前的压气机流动DES类研究工作中仍然普遍采用高耗散的迎风格式。本文首先基于DES类方法计算的各向同性衰减湍流结果,定量比较了多种不同数值格式的耗散,证实了高耗散迎风格式严重低估中高波数湍流能量。高阶重构格式可以一定程度上改善该问题,但能量耗散仍然过高。本文在高阶重构的基础上,进一步引入自适应耗散函数修改Riemann求解器,构造了自适应耗散格式。该格式在全波数范围都能准确地预测湍流能谱。将该格式配合DES类方法模拟跨声速离心压气机流动,其预测的压比相比于三阶迎风格式,更加接近实验结果。此外,自适应耗散格式显著提高了中小尺度流动结构的分辨率。分析表明,在使用DES类方法模拟压气机流动时,有必要采用数值耗散较低的离散格式,以准确预测压气机总体性能和流动结构。本文构造的自适应耗散格式是一种良好选择。     

Numerical analysis and construction of limiter of high resolution difference scheme

I.IntroductionAdjustingandgoverningproperlythenumericaldissipationanddispersionarethekeytotheconstructionofhighresolutionandnon-oscillationscheme.SinceVanLcerll]introducednuxlimiterandobtainedhighresolutionandnon-oscillationscheme,choisillgandconstructingsuitablelimiterbecomesanimportantwayofdesigninghighresolutiollalld11on-oscillationscheme,andadjustsnumericaldissipationeffects,dispersioneffectsaswellasgroupvelocityeffects.Roe121'l'l,Chakravarthy'andOSherl4]andetal.providedvariouslimitersan…     

Novel local smoothness indicators for improving the third‐order WENO scheme

The local smoothness indicators play an important role in the performance of a weighted essentially nonoscillatory (WENO) scheme. Due to having only 2 points available on each substencil, the local smoothness indicators calculated by conventional methods make the third‐order WENO scheme too dissipative. In this paper, we propose a different method to calculate the indicators by using all the 3 points on the global stencil of the third‐order WENO scheme. The numerical results demonstrate that the WENO scheme with the new indicators has less dissipation and better resolution than the conventional third‐order WENO scheme of Jiang and Shu for both smooth and discontinuous solutions.     

Manipulation of particles using dielectrophoresis

A numerical scheme based on the distributed Lagrange multiplier method (DLM) is used to study the motion of particles of electrorheological suspensions subjected to non-uniform electric fields. At small Reynolds number, the time taken by the particles to collect at the minimums or maximums of the electric field is primarily determined by a parameter defined to be the ratio of the dielectrophoretic and viscous forces. Simulations show that in non-uniform electric fields the collection time is also influenced by a parameter defined by the ratio of the electrostatic particle–particle interaction and dielectrophoretic forces. The collection time decreases as this parameter decreases because when this parameter is less than one, particles move to the regions of high or low electric field regions individually. However, when this parameter is greater than one, particles regroup into chains which then move toward the electric field maximums or minimums without breaking. It is also shown that when the real part of the Clausius–Mosotti factor (β) is negative the positions of the local minimums of the electric field, and thus also the locations where particles collect, can be modified by changing the electric potential boundary conditions.     

Optimized sixth‐order monotonicity‐preserving scheme by nonlinear spectral analysis

In this paper, sixth‐order monotonicity‐preserving optimized scheme (OMP6) for the numerical solution of conservation laws is developed on the basis of the dispersion and dissipation optimization and monotonicity‐preserving technique. The nonlinear spectral analysis method is developed and is used for the purpose of minimizing the dispersion errors and controlling the dissipation errors. The new scheme (OMP6) is simple in expression and is easy for use in CFD codes. The suitability and accuracy of this new scheme have been tested through a set of one‐dimensional, two‐dimensional, and three‐dimensional tests, including the one‐dimensional Shu–Osher problem, the two‐dimensional double Mach reflection, and the Rayleigh–Taylor instability problem, and the three‐dimensional direct numerical simulation of decaying compressible isotropic turbulence. All numerical tests show that the new scheme has robust shock capturing capability and high resolution for the small‐scale waves due to fewer numerical dispersion and dissipation errors. Moreover, the new scheme has higher computational efficiency than the well‐used WENO schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

Positivity‐preserving,flux‐limited finite‐difference and finite‐element methods for reactive transport

A new class of positivity‐preserving, flux‐limited finite‐difference and Petrov–Galerkin (PG) finite‐element methods are devised for reactive transport problems.The methods are similar to classical TVD flux‐limited schemes with the main difference being that the flux‐limiter constraint is designed to preserve positivity for problems involving diffusion and reaction. In the finite‐element formulation, we also consider the effect of numerical quadrature in the lumped and consistent mass matrix forms on the positivity‐preserving property. Analysis of the latter scheme shows that positivity‐preserving solutions of the resulting difference equations can only be guaranteed if the flux‐limited scheme is both implicit and satisfies an additional lower‐bound condition on time‐step size. We show that this condition also applies to standard Galerkin linear finite‐element approximations to the linear diffusion equation. Numerical experiments are provided to demonstrate the behavior of the methods and confirm the theoretical conditions on time‐step size, mesh spacing, and flux limiting for transport problems with and without nonlinear reaction. Copyright © 2003 John Wiley & Sons, Ltd.     

One‐dimensional shock‐capturing for high‐order discontinuous Galerkin methods

Discontinuous Galerkin methods have emerged in recent years as an alternative for nonlinear conservation equations. In particular, their inherent structure (a numerical flux based on a suitable approximate Riemann solver introduces some stabilization) suggests that they are specially adapted to capture shocks. However, numerical fluxes are not sufficient to stabilize the solution in the presence of shocks. Thus, slope limiter methods, which are extensions of finite volume methods, have been proposed. These techniques require, in practice, mesh adaption to localize the shock structure. This is is more obvious for large elements typical of high‐order approximations. Here, a new approach based on the introduction of artificial diffusion into the original equations is presented. The order is not systematically decreased to one in the presence of the shock, large high‐order elements can be used, and several linear and nonlinear tests demonstrate the efficiency of the proposed methodology. Copyright © 2012 John Wiley & Sons, Ltd.     

A nonlocal convective flux limiter for upwind‐biased finite volume simulations

Finite volume methods on structured and unstructured meshes commonly use second‐order, upwind‐biased linear reconstruction schemes to approximate the convective terms. Limiters are employed to reduce the inherent variable overshoot and undershoot of these schemes in discontinuous regions; however, they also can significantly increase the numerical dissipation of a solution in smooth regions. This paper presents a novel nonlocal, nonmonotonic (NLNM) limiter developed by enforcing cell minima and maxima on dependent variable values projected to cell faces. The minimum and maximum values for a cell are determined primarily through recursive reference to the minimum and maximum values of its upwind neighbors. The new limiter has been implemented into an existing flow solver. Various test cases are presented, which highlight the ability of the NLNM limiter to eliminate nonphysical oscillations while maintaining negligible levels of limiter‐induced dissipation into the solution. Results from the new limiter are compared with those from other limited and unlimited second‐order upwind and first‐order upwind schemes. For the cases examined in the study, the NLNM limiter was found to improve accuracy without significantly decreasing overall simulation efficiency and robustness.Copyright © 2011 John Wiley & Sons, Ltd.     

Efficient high‐resolution relaxation schemes for hyperbolic systems of conservation laws

In this work, we present a total variation diminishing (TVD) scheme in the zero relaxation limit for nonlinear hyperbolic conservation law using flux limiters within the framework of a relaxation system that converts a nonlinear conservation law into a system of linear convection equations with nonlinear source terms. We construct a numerical flux for space discretization of the obtained relaxation system and modify the definition of the smoothness parameter depending on the direction of the flow so that the scheme obeys the physical property of hyperbolicity. The advantages of the proposed scheme are that it can give second‐order accuracy everywhere without introducing oscillations for 1‐D problems (at least with) smooth initial condition. Also, the proposed scheme is more efficient as it works for any non‐zero constant value of the flux limiter ? ? [0, 1], where other TVD schemes fail. The resulting scheme is shown to be TVD in the zero relaxation limit for 1‐D scalar equations. Bound for the limiter function is obtained. Numerical results support the theoretical results. Copyright © 2007 John Wiley & Sons, Ltd.     

The composite finite volume method on unstructured meshes for the two‐dimensional shallow water equations

A composite finite volume method (FVM) is developed on unstructured triangular meshes and tested for the two‐dimensional free‐surface flow equations. The methodology is based on the theory of the remainder effect of finite difference schemes and the property that the numerical dissipation and dispersion of the schemes are compensated by each other in a composite scheme. The composite FVM is formed by global composition of several Lax–Wendroff‐type steps followed by a diffusive Lax–Friedrich‐type step, which filters out the oscillations around shocks typical for the Lax–Wendroff scheme. To test the efficiency and reliability of the present method, five typical problems of discontinuous solutions of two‐dimensional shallow water are solved. The numerical results show that the proposed method, which needs no use of a limiter function, is easy to implement, is accurate, robust and is highly stable. Copyright © 2001 John Wiley & Sons, Ltd.     

An improvement of classical slope limiters for high‐order discontinuous Galerkin method

In this paper, we describe some existing slope limiters (Cockburn and Shu's slope limiter and Hoteit's slope limiter) for the two‐dimensional Runge–Kutta discontinuous Galerkin (RKDG) method on arbitrary unstructured triangular grids. We describe the strategies for detecting discontinuities and for limiting spurious oscillations near such discontinuities, when solving hyperbolic systems of conservation laws by high‐order discontinuous Galerkin methods. The disadvantage of these slope limiters is that they depend on a positive constant, which is, for specific hydraulic problems, difficult to estimate in order to eliminate oscillations near discontinuities without decreasing the high‐order accuracy of the scheme in the smooth regions. We introduce the idea of a simple modification of Cockburn and Shu's slope limiter to avoid the use of this constant number. This modification consists in: slopes are limited so that the solution at the integration points is in the range spanned by the neighboring solution averages. Numerical results are presented for a nonlinear system: the shallow water equations. Four hydraulic problems of discontinuous solutions of two‐dimensional shallow water are presented. The idealized dam break problem, the oblique hydraulic jump problem, flow in a channel with concave bed and the dam break problem in a converging–diverging channel are solved by using the different slope limiters. Numerical comparisons on unstructured meshes show a superior accuracy with the modified slope limiter. Moreover, it does not require the choice of any constant number for the limiter condition. Copyright © 2008 John Wiley & Sons, Ltd.     

A new weighting method for improving the WENO‐Z scheme

The calculation of the weight of each substencil is very important for a weighted essentially nonoscillatory (WENO) scheme to obtain high‐order accuracy in smooth regions and keep the essentially nonoscillatory property near discontinuities. The weighting function introduced in the WENO‐Z scheme provides a straightforward method to analyze the accuracy order in smooth regions. In this paper, we construct a new sixth‐order global smoothness indicator (GSI‐6) and a function about GSI‐6 and the local smoothness indicators (IS_k) to calculate the weights. The analysis and numerical results show that, with the new weights, the scheme satisfies the sufficient condition for the fifth‐order convergence in smooth regions even at critical points. Meanwhile, it can also maintain low dissipation for discontinuous solutions due to relative large weights assigned to discontinuous substencils.     

A robust low‐dissipation AUSM‐family scheme for numerical shock stability on unstructured grids

The simple low‐dissipation advection upwind splitting method (SLAU) scheme is a parameter‐free, low‐dissipation upwind scheme that has been applied in a wide range of aerodynamic numerical simulations. In spite of its successful applications, the SLAU scheme could be showing shock instabilities on unstructured grids, as many other contact resolved upwind schemes. Therefore, a hybrid upwind flux scheme is devised for improving the shock stability of SLAU scheme, without compromising on accuracy and low Mach number performance. Numerical flux function of the hybrid scheme is written in a general form, in which only the scalar dissipation term is different from that of the SLAU scheme. The hybrid dissipation term is defined by using a differentiable multidimensional‐shock‐detection pressure weight function, and the dissipation term of SLAU scheme is combined with that of the Van Leer scheme. Furthermore, the hybrid dissipation term is only applied for the solution of momentum fluxes in numerical flux function. Based on the numerical test results, the hybrid scheme is deemed to be a successful improvement on the shock stability of SLAU scheme, without compromising on the efficiency and accuracy. Copyright © 2016 John Wiley & Sons, Ltd.     

Implementing a high‐order accurate implicit operator scheme for solving steady incompressible viscous flows using artificial compressibility method

This paper uses a fourth‐order compact finite‐difference scheme for solving steady incompressible flows. The high‐order compact method applied is an alternating direction implicit operator scheme, which has been used by Ekaterinaris for computing two‐dimensional compressible flows. Herein, this numerical scheme is efficiently implemented to solve the incompressible Navier–Stokes equations in the primitive variables formulation using the artificial compressibility method. For space discretizing the convective fluxes, fourth‐order centered spatial accuracy of the implicit operators is efficiently obtained by performing compact space differentiation in which the method uses block‐tridiagonal matrix inversions. To stabilize the numerical solution, numerical dissipation terms and/or filters are used. In this study, the high‐order compact implicit operator scheme is also extended for computing three‐dimensional incompressible flows. The accuracy and efficiency of this high‐order compact method are demonstrated for different incompressible flow problems. A sensitivity study is also conducted to evaluate the effects of grid resolution and pseudocompressibility parameter on accuracy and convergence rate of the solution. The effects of filtering and numerical dissipation on the solution are also investigated. Test cases considered herein for validating the results are incompressible flows in a 2‐D backward facing step, a 2‐D cavity and a 3‐D cavity at different flow conditions. Results obtained for these cases are in good agreement with the available numerical and experimental results. The study shows that the scheme is robust, efficient and accurate for solving incompressible flow problems. Copyright © 2010 John Wiley & Sons, Ltd.     

On the properties of discrete adjoints of numerical methods for the advection equation

This paper discusses several aspects related to the consistency of the discrete adjoints of upwind numerical schemes. Both linear (finite differences, finite volumes) and nonlinear (slope and flux‐limited) discretizations of the one‐dimensional advection equation are considered. The analysis is focused on uniform meshes and on explicit numerical schemes. We show that the discrete adjoints may lose consistency near the points where upwinding changes, near inflow boundaries where another numerical scheme is employed, and near the locations where the slope/flux limiter is active in the forward simulation. Numerical results are presented to support the theoretical analysis. Copyright © 2007 John Wiley & Sons, Ltd.     

Invariant submodels of the Westervelt model with dissipation

We study three-dimensional Westervelt model of nonlinear hydroacoustics with dissipation. We received all its invariant submodels. With the help of invariant solutions, we explored some wave processes, specifying their physical meaning. The boundary value problems describing these processes are reduced to the nonlinear integro-differential equations. We established the existence and uniqueness of the solutions of these boundary value problems under some additional conditions. Also we considered the invariant solutions of rank 2 and 3. Mechanical relevance of the obtained solutions is as follows: (1) these solutions describe nonlinear and diffraction effects in ultrasonic fields of a special kind, (2) these solutions can be used as a test solutions in the numerical calculations performed in studies of ultrasonic fields generated by powerful emitters.     

High‐order compact numerical schemes for non‐hydrostatic free surface flows

The development of a numerical scheme for non‐hydrostatic free surface flows is described with the objective of improving the resolution characteristics of existing solution methods. The model uses a high‐order compact finite difference method for spatial discretization on a collocated grid and the standard, explicit, single step, four‐stage, fourth‐order Runge–Kutta method for temporal discretization. The Cartesian coordinate system was used. The model requires the solution of two Poisson equations at each time‐step and tridiagonal matrices for each derivative at each of the four stages in a time‐step. Third‐ and fourth‐order accurate boundaries for the flow variables have been developed including the top non‐hydrostatic pressure boundary. The results demonstrate that numerical dissipation which has been a problem with many similar models that are second‐order accurate is practically eliminated. A high accuracy is obtained for the flow variables including the non‐hydrostatic pressure. The accuracy of the model has been tested in numerical experiments. In all cases where analytical solutions are available, both phase errors and amplitude errors are very small. Copyright © 2006 John Wiley & Sons, Ltd.     

High‐order shock capturing schemes for turbulence calculations

This work investigates high‐order central compact methods for simulating turbulent supersonic flows that include shock waves. Several different types of previously proposed characteristic filters, including total variation diminishing, monotone upstream‐centered scheme for conservation laws, and weighted essentially non‐oscillatory filters, are investigated in this study. Similar to the traditional shock capturing schemes, these filters can eliminate the numerical instability caused by large gradients in flow fields, but they also improve efficiency compared with classical shock‐capturing schemes. Adding the nonlinear dissipation part of a classical shock‐capturing scheme to a central scheme makes the method suitable for incorporation into any existing central‐based high‐order subsonic code. The amount of numerical dissipation to add is sensed by means of the artificial compression method switch. In order to improve the performance of the characteristic filters, we propose a hybrid approach to minimize the dissipation added by the characteristic filter. Through several numerical experiments (including a shock/density wave interaction, a shock/vortex interaction, and a shock/mixing layer interaction) we show that our hybrid approach works better than the original method, and can be used for future turbulent flow simulations that include shocks. Copyright © 2009 John Wiley & Sons, Ltd.     

A finite volume‐based high‐order,Cartesian cut‐cell method for wave propagation

Computational aeroacoustics requires numerical techniques capable of yielding low artificial dispersion and dissipation to preserve the amplitude and the frequency characteristics of the physical processes. Furthermore, for engineering applications, the techniques need to handle irregular geometries associated with realistic configurations. We address these issues by developing an optimized prefactored compact finite volume (OPC‐fv) scheme along with a Cartesian cut‐cell technique. The OPC‐fv scheme seeks to minimize numerical dispersion and dissipation while satisfying the conservation laws. The cut‐cell approach treats irregularly shaped boundaries using divide‐and‐merge procedures for the Cartesian cells while maintaining a desirable level of accuracy. We assess these techniques using several canonical test problems, involving different levels of physical and geometric complexities. Richardson extrapolation is an effective tool for evaluating solutions of no high gradients or discontinuities, and is used to evaluate the performance of the solution technique. It is demonstrated that while the cut‐cell method has a modest effect on the order of accuracy, it is a robust method. The combined OPC‐fv scheme and the Cartesian cut‐cell technique offer good accuracy as well as geometric flexibility. Copyright © 2008 John Wiley & Sons, Ltd.