Abstract:High order upwind-biased schemes based on the general explicit difference formulas with arbitrary order accuracy from Li (2005) are developed. The dissipation and dispersion errors of these upwind-biased schemes are assessed by means of the Fourier analysis. The results show that the numerical phase speeds for even order schemes are faster than the real phase speed. But the numerical phase speeds for odd order schemes are slower than the real phase speed. Moreover, the dissipation and dispersion errors of the even order scheme are smaller than the next lower-order odd scheme. For the upwind-biased schemes derived, there is an increase of the dissipation error in the high wave number range. However, the increase of numerical dissipation in the low wave number range is very small. Some numerical dissipation is needed in the high wave number range for some flow problems. The decision that a numerical modeler should always make is how much numerical dissipation is need. The upwind-biased schemes provide some additional alternatives for adjusting numerical dissipation. These schemes are applied to a one-dimensional model in order to test their computational performance. Firstly, consider one-dimensional advection in a constant velocity field. The two types of test problems that are used to evaluate the accuracy of a numerical scheme are the Gaussian function and the square wave function. The first example is the advection of the square wave function. This function reveals a numerical method's capability to handle Gibb's oscillations that arise in the vicinity of discontinuities. The second example is the advection of the Gaussian function. This function provides a demonstration of each scheme's ability to transport well resolved, smoothly varying functions over larger distances. Testing results show that the errors of numerical schemes gradually decrease as the order of the scheme increases and further changes are moderate for higher than sixth order scheme. A high order scheme, although not maintaining a positive definite field, will reduce the problem of negative quantities dramatically. Secondly, consider the inviscid Burgers' equation. The inviscid Burgers' equation is used as the test model because it is the simplest equation that allows scale collapse (shock formation), and because it has analytic solutions. There are many important phenomena in the atmosphere that are associated with sharp gradients. One requirement of a computational method for the scale collapse phenomena would seem to be an ability to accurately represent a very sharp gradient in the computational domain. In addition, any numerical dispersion error due to the presence of a shock should not contaminate the smooth solution away from the shock. Numerical results show that there is a drastic improvement when going from the first to sixth order, while further changes are moderate. The results indicate the usefulness of the upwind-biased scheme in handling the advection near discontinuities and the development of scale collapse regions in atmospheric numerical models. In conclusion, the sixth-order upwind-biased scheme appears to be the best balance between efficiency and accuracy. It appears that the sixth-order upwind-biased scheme is well suited for many atmospheric modeling applications where advection plays a significant role.

Effects Of Differencing Schemes On Simulation Of Dense Gas-Particle Two-Phase Flows

Advances in Cryogenic Engineering: Volumes 35A and B

Difference scheme for two-phase flow

An Efficient Difference Scheme for Gas Dynamics Equations

FLUID SIMULATION USING AN ADAPTIVE SEMI-DISCRETE CENTRAL-UPWIND SCHEME

Numerical Simulation of Hydrodynamics in Downers Using a CFD–DEM Coupled Approach

Eulerian-Lagrangian simulation of hydrodynamics of gas-solid flows in downer

Comparison of Several Upwind Schemes in Transonic Flows

Numerical simulation for two-phase flows using hybrid scheme

Improvement of Unresolved CFD-DEM by Velocity Field Reconstruction on Unstructured Grids

A Comparison and Analysis of High Order Upwind-Biased Schemes

Accuracy and Stability Enhancements in the Incompressible Finite-Volume-particle Method for Multiphase Flow Simulations.

Simulation of two-phase flows at large density ratios and high Reynolds numbers using a discrete unified gas kinetic scheme

A class of high-order finite difference schemes with minimized dispersion and adaptive dissipation for solving compressible flows

Comparison of Finite Difference Schemes and Physical Conservation Features

Numerical Dissipation Effects on Massive Separation Around Tandem Cylinders

Optimized Low-Dissipation and Low-Dispersion Schemes for Compressible Flows

A Comparison of Two-Fluid Model, Dense Discrete Particle Model and CFD-DEM Method for Modeling Impinging Gas–solid Flows

A Two‐layer Mesh Method for Discrete Element Simulation of Gas‐particle Systems with Arbitrarily Polyhedral Mesh

Assessment of Model Formulations in the Discrete Particle Simulation of Gas−Solid Flow

A Critical Comparison of Two-Fluid Model, Discrete Particle Method and Direct Numerical Simulation for Modeling Dense Gas-Solid Flow of Rough Spheres