High Order Weighted Essentially Non-oscillatory Research Articles

Computation of three-dimensional incompressible flows using high-order weighted essentially non-oscillatory finite-difference lattice Boltzmann method

In the present work, an accurate and robust solution methodology based on the high-order weighted essentially non-oscillatory (WENO) finite-difference lattice Boltzmann method (LBM) in the three-dimensional generalized curvilinear coordinates is presented and applied for simulating the three-dimensional incompressible flows over complicated configurations with curved boundaries. Here, the incompressible form of the lattice Boltzmann equation in three dimensions is considered and the discretization of the spatial derivative terms is performed with the fifth-order WENO finite-difference method and the temporal derivative term is discretized with the fourth-order Runge–Kutta scheme to ensure the accuracy and stability of the solution method for both the steady and unsteady problems. The three-dimensional lattice Boltzmann equation applied here is based on a nineteen discrete velocity model for transforming the microscopic properties to the macroscopic ones. To assess the accuracy and robustness of the present three-dimensional high-order finite-difference LBM solver, different incompressible flow benchmarks and practical test cases are studied that are the cavity flow, the Beltrami flow, the flow in the curved ducts of rectangular cross sections, and the flow over a sphere for different flow conditions. The decay of the homogeneous isotropic turbulence is also computed to examine the suitability of the present solution method to be applied as the direct numerical simulation of turbulent flows. It is demonstrated that the solution methodology presented based on the high-order WENO finite-difference LBM in the three-dimensional generalized curvilinear coordinate can be used for accurately and effectively computing the three-dimensional practical incompressible flow problems.

Numerical study of liquid jet and shock wave induced by two-bubble collapse in open field

The shock waves and high-speed jets induced by the violent collapse of bubbles can cause serious damage to nearby structures. To fully understand the shock wave and jet characteristics induced by bubble-bubble interaction, we developed a diffuse-interface method for the simulation of compressible two-phase flow, which combined interface compression technique with interface sharpening technique to keep the sharpness of the shock wave and phase interface. The thermodynamically consistent five-equation model was improved using interface compression technique to maintain a constant thickness of the phase interface. To further suppress numerical dissipation, the high-order scheme WENO (Weighted Essentially Non-Oscillatory) and interface sharpening function THINC (Tangent of Hyperbola for INterface Capturing) were used for the flux reconstruction at the grid boundary to ensure that the total variation at the cell boundary was minimized (Boundary Variation Diminishing, BVD), thus the spatial reconstruction scheme, named WENO-THINC-BVD, was developed. Moreover, we extended the present method to a block-structured adaptive mesh refinement framework to improve grid resolution and save computing resources. Numerical results of several benchmark tests fully demonstrated the advantages of the present method in simulating complex compressible flows containing shock-shock and shock-interface interactions. Based on the present high-fidelity numerical methods, we investigated the shock wave and liquid jet induced by the two-bubble interaction in an open space. A phase diagram of the liquid jet propagation direction as a function of the initial bubble pressure and the inter-bubble distance was summarized. Finally, according to the bubble-bubble interaction process, an empirical formula for the collapse strength and another for the attenuation of shock waves were developed.

A new family of downwind-limited, scale-invariant WENO schemes with optimal accuracy

Despite improvement over decades, there are still at least two key issues not fully solved in Weighted Essentially Non-oscillatory (WENO) schemes. Firstly, most WENO schemes hardly satisfy the scale-invariant property, which undermines their robustness and accuracy when simulating flows involving different scales such as reactive multi-component flows and plasma. Secondly, most of existing Essentially Non-oscillatory (ENO) constraints are not sufficient and necessary conditions to achieve optimal accuracy. Therefore, in this work, we design a new family of WENO schemes based on a three-cell stencil to address these two issues. The new WENO non-linear weights are designed as bounded asymptotic functions which are dependent on cell Normalised Volume Integrated Average (NVIA) values and a shape parameter θ. Thus, the new WENO weighting functions are scale-independent. The shape parameter θ controls the errors between the nonlinear weight and the ideal weight. As the value of θ increases, the new weighting functions become more accurate but also become more likely to produce numerical oscillations. Thus, we propose a downwind limiting method as a supplementary constraint of ENO property. The maximum value of θ which corresponds to the optimal accuracy of a weighting function can then be determined. The accuracy analysis via the standard Taylor expansion at critical points reveals that the proposed weighting functions outperform classical WENO schemes. The accuracy, scale-invariant property and ENO property of proposed schemes are further validated through benchmark tests, which shows the superior performance of the new schemes in comparison with existing three-cell-based WENO schemes. Thus, this work proposes a new direction to construct scale-invariant WENO weighting functions with optimal accuracy, which can potentially be extended to higher-order WENO schemes.

A class of high-order improved fast weighted essentially non-oscillatory schemes for achieving optimal order at any critical points

The weighted essentially non-oscillatory (WENO) scheme is one of the most popular shock-capturing schemes, and constructing a more efficient and higher-order WENO scheme has always been an intention of optimization design. In the general WENO reconstruction framework, the smoothness indicator plays an important role in identifying whether the sub-stencils are in discontinuous or smooth regions. However, the classical smoothness indicator is the most expensive one in the whole reconstruction algorithm, and its computational complexity increases sharply with the improvement of the accuracy order. Therefore, a class of efficient and superior WENO schemes called improved fast WENO (IFWENO) are proposed based on the fast WENO (FWENO). To improve efficiency, the smoothness indicator of the IFWENO scheme is simplified from the traditional version, and the nonlinear weight calculation method is modified. The parameter ε is carefully designed to obtain the superior property that the accuracy of the spatial derivatives will not degrade at any order critical point in smooth regions. The reason for the instability occurring in the high-order FWENO is revealed, and the parameter p is likewise specifically selected to improve robustness at discontinuities. The excellent multi-scale resolution of the proposed IFWENO scheme is proven by theoretical analyses and numerical experiments. Through several typical examples, the consistently high accuracy and efficiency of the designed scheme in both smooth and discontinuous regions are verified.

Efficient implementation of high-order WENO schemes with sharing function for solving Euler equations

Due to the high-order accuracy and essentially non-oscillatory (ENO) property, the weighted ENO (WENO) schemes have a wide range of successful applications. The component-wise reconstruction WENO (CP WENO) scheme for fluxes or variables has simple formulations but it may produce numerical oscillations near discontinuities when solving the Euler equations. Although the characteristic-wise reconstruction WENO (CH WENO) scheme can reduce such oscillations, it involves too many characteristic projection operations. In this paper, first, we introduced a sharing function to indicate the discontinuities in the Euler equations and then constructed new adaptive characteristic-wise WENO (Ada-WENO) scheme and common-weights WENO (Co-WENO) scheme with this function. Several one and two dimensional problems are used to test the performances of Ada-WENO and Co-WENO. Numerical results show that, Ada-WENO can reduce the computational cost of CH WENO while maintaining its oscillation-free property since it only switches from CP WENO to CH WENO near discontinuities, and Co-WENO can reduce the cost and oscillations of CP WENO, but it may still generate few oscillations.

A modified high-order symmetrical WENO scheme for hyperbolic conservation laws

Abstract This paper designs a modified weighted essentially non-oscillatory (WENO) scheme for solving hyperbolic conservation laws. Using the switching principle based on inflection points, the new scheme automatically adapts between linear upwind and WENO schemes. If there is at least one inflection point in the largest stencil available for reconstruction, a symmetrical WENO (SWENO) scheme is considered for the reconstruction of the numerical flux; otherwise the numerical flux is directly approximated by the reconstruction polynomial. By comparing the new scheme introduced in this paper with the classical WENO scheme and another improved scheme that has been proposed recently D. Chai, G. Xi, Z. Sun, Z. Wangand Z. Huang,Comput. Fluids, vol. 170, pp. 176–186, 2018), we can point out the robustness and better efficiency of this scheme. To examine and explain the features of the new scheme, a number of examples such as Euler equations have been considered.

On high-order numerical schemes for viscous relativistic hydrodynamics through the Kelvin–Helmholtz instability

ABSTRACT This work assesses the dissipative properties of high-order numerical methods for relativistic hydrodynamics. A causal theory of physical dissipation is included within a finite volume high-resolution shock-capturing framework based on the Israel–Stewart theory to study high-order WENO (weighted-essentially non-oscillatory) schemes for simulating the relativistic Kelvin–Helmholtz instability. We provide an estimation of the numerical dissipation of high-order schemes based on results obtained both with and without physically resolved dissipation and determine an empirical relationship between the numerical dissipation and the grid resolution. We consider the appearance of secondary flow features within the evolution of the Kelvin–Helmholtz instability and determine that they are numerical artifacts — this is partly based on arguments presented in terms of a frame-dependent form of the relativistic Reynolds number. There is a potential advantage of using high-order schemes in terms of their accuracy and computational cost on coarser grid resolutions when directly compared to low-order schemes on a fine grid in the presence of physical viscosity. It is possible to find reasonable agreement between numerical results that employ lower-order schemes using a finer grid resolution and results that employ higher order schemes at a coarser grid resolution when sufficient viscosity is present. Overall, the present analysis gives an insight into the numerical dissipation of high-order shock-wave capturing schemes which can be relevant to computational studies of astrophysical phenomena in the relativistic regime. The results presented herein are problem and scheme-dependent and serve to highlight the different roles of numerical and physical dissipation.

Preconditioned Jacobian-free Newton–Krylov fully implicit high order WENO schemes and flux limiter methods for two-phase flow models

Motivated by the high-resolution properties of high-order Weighted Essentially Non-Oscillatory (WENO) and flux limiter (FL) for steep-gradient problems and the robust convergence of Jacobian-free Newton-Krylov (JFNK) methods for nonlinear systems, the preconditioned JFNK fully implicit high-order WENO and FL schemes are proposed to solve the transient two-phase two-fluid models. Specially, the second-order fully-implicit BDF2 is used for the temporal operator and then the third-order WENO schemes and various flux limiters can be adopted to discrete the spatial operator. For the sake of the generalization of the finite-difference-based preconditioning acceleration methods and the excellent convergence to solve the complicated and various operational conditions, the random vector instead of the initial condition is skillfully chosen as the solving variables to obtain better sparsity pattern or more positions of non-zero elements in this paper. Finally, the WENO_JFNK and FL_JFNK codes are developed and then the two-phase steep-gradient problem, phase appearance/disappearance problem, U-tube problem and linear advection problem are tested to analyze the convergence, computational cost and efficiency in detailed. Numerical results show that WENO_JFNK and FL_JFNK can significantly reduce numerical diffusion and obtain better solutions than traditional methods. WENO_JFNK gives more stable and accurate solutions than FL_JFNK for the test problems and the proposed finite-difference-based preconditioning acceleration methods based on the random vector can significantly improve the convergence speed and efficiency.

Fast Sparse Grid Simulations of Fifth Order WENO Scheme for High Dimensional Hyperbolic PDEs

The weighted essentially non-oscillatory (WENO) schemes, especially the fifth order WENO schemes, are a popular class of high order accurate numerical methods for solving hyperbolic partial differential equations (PDEs). However when the spatial dimensions are high, the number of spatial grid points increases significantly. It leads to large amount of operations and computational costs in the numerical simulations by using nonlinear high order accuracy WENO schemes such as a fifth order WENO scheme. How to achieve fast simulations by high order WENO methods for high spatial dimension hyperbolic PDEs is a challenging and important question. In the literature, sparse-grid technique has been developed as a very efficient approximation tool for high dimensional problems. In a recent work [Lu, Chen and Zhang, Pure and Applied Mathematics Quarterly, 14 (2018) 57-86], a third order finite difference WENO method with sparse-grid combination technique was designed to solve multidimensional hyperbolic equations including both linear advection equations and nonlinear Burgers’ equations. Numerical experiments showed that WENO computations on sparse grids achieved comparable third order accuracy in smooth regions of the solutions and nonlinear stability as that for computations on regular single grids. In application problems, higher than third order WENO schemes are often preferred in order to efficiently resolve the complex solution structures. In this paper, we extend the approach to higher order WENO simulations specifically the fifth order WENO scheme. A fifth order WENO interpolation is applied in the prolongation part of the sparse-grid combination technique to deal with discontinuous solutions. Benchmark problems are first solved to show that significant CPU times are saved while both fifth order accuracy and stability of the WENO scheme are preserved for simulations on sparse grids. The fifth order sparse grid WENO method is then applied to kinetic problems modeled by high dimensional Vlasov based PDEs to further demonstrate large savings of computational costs by comparing with simulations on regular single grids. Several open problems are discussed at last.

Very high order WENO schemes using efficient smoothness indicators

Very high order weighted essentially non-oscillatory (WENO) schemes can obtain more accurate results than those relatively low order WENO schemes for many problems. However, as the order increases, the candidate sub-stencils of WENO schemes will contain more node points. Then, the classical smoothness indicators become very complex and need much more floating point operations. Although they could be written in a more succinct form which is a sum of perfect squares, only a little improvement of the computational efficiency is obtained. Furthermore, for the power parameter p in the nonlinear weights of the previous (2r-1)th order WENO schemes, only the range of choices is given and the optimal value could be problem dependent. In this article, another class of very high order WENO schemes is constructed by using more efficient smoothness indicators. Compared with very high order WENO schemes using classical smoothness indicators, the new ones are more time efficient, and less error prone for programming. Furthermore, they have the same approximate dispersion relation (ADR) as their underlying linear schemes. And a more important advantage is that they have both good behavior near critical points and good stability near discontinuities. For Euler equations, a predetermined order reduction technique including two inequalities is proposed to deal with the problems caused by interactions between discontinuities of different characteristic fields. One inequality is constructed using this kind of smoothness indicators of pressure to suppress the non-physical oscillations. And the other is proposed according to the governing equation of the pressure to prevent the emergence of negative pressure. Several test problems are presented to demonstrate the good properties of the proposed WENO schemes.

Direct numerical simulation of supersonic flow and acoustics over a compression ramp

We present direct numerical simulations of the shock wave boundary layer interaction (SBLI) at Mach number 2.9 over a 24° ramp. We study both the numerical accuracy and flow physics. Two classes of spatial reconstruction schemes are employed: the monotonic upstream-centered scheme for conservation laws and the Weighted Essentially Non-Oscillatory (WENO) scheme, of accuracy ranging from 2nd- to 11th-order. Using the canonical Taylor–Green vortex test-case, a simple and computationally inexpensive rescaling of the candidate stencil values—within the context of the high-order WENO scheme—is proposed for reducing the numerical dissipation, particularly in under-resolved simulations. For the compression ramp case, higher-order schemes are shown to capture the size of the SBLI separation zone more accurately, a consequence of resolving much finer turbulence structures. For second- and fifth-order schemes, the energy of the unresolved small scale turbulence shifts the cascade of the turbulence kinetic energy (TKE) spectrum, thus resulting in more energetic large scale turbulent structures. Consequently, the λ-shock foot shifts further downstream, leading to a smaller separation bubble size. Nonetheless, other statistical quantities, such as the turbulence anisotropy invariant map and the turbulence kinetic energy budget terms, show little dependence on the type and order of the spatial reconstruction scheme. Finally, using the more accurate ninth-order WENO results, it is reasoned that the interaction of the λ-shock with the post-shock relaxation region drives the low-frequency oscillation of the λ-shock.

High order well-balanced finite difference WENO interpolation-based schemes for shallow water equations

A numerical framework of the generalized form of high order well-balanced finite difference weighted essentially non-oscillatory (WENO) interpolation-based schemes is proposed for the shallow water equations. It demonstrates more flexible construction process than the classical WENO reconstruction-based schemes. The weighted compact nonlinear schemes and finite difference alternative WENO schemes are two specific cases. To maintain the exact C-property, the splitting technique for the source term in the finite difference scheme [Xing and Shu, J. Comput. Phys. 208 (2005)] and the reconstruction technique in the finite volume WENO scheme [Xing and Shu, J. Comput. Phys. 214 (2006)] are adopted. The proposed scheme can be proved mathematically to maintain the exact C-property and demonstrates numerically that it is well-balanced by construction for the stationary water surface. Moreover, the local characteristic projections are employed to further mitigate the Gibbs oscillations. The proposed generic high order WENO schemes not only achieve high order accuracy but also capture the high gradients/shock waves essentially non-oscillatory. Meanwhile, the small perturbation problems can be resolved well on a coarse grid.

A new method towards high-order weno schemes on structured and unstructured grids

A new method is presented to develop high-order weighted essentially non-oscillatory (WENO) finite-volume schemes on structured and unstructured grids. The one-dimensional reconstruction used in the WENO-ZQ scheme is extended by introducing the concept of phantom points, which are constructed with the least-squares approximation to form a series of high-order WENO schemes. Numerical validations are performed on both regular structured grids and irregular unstructured grids. Analyses of the spectral property and numerical accuracy of the new schemes indicate attractive merits, including less dissipation and dispersion errors and low computational costs, which can be easily implemented into existing structured and unstructured grid CFD solvers.

A new type of third-order finite volume multi-resolution WENO schemes on tetrahedral meshes

In this continuing paper of Zhu and Shu (2018) [62], Zhu and Shu (2019) [63], we design a new third-order finite volume multi-resolution weighted essentially non-oscillatory (WENO) scheme for solving hyperbolic conservation laws on tetrahedral meshes. We only use the information defined on a hierarchy of nested central spatial stencils without introducing any equivalent multi-resolution representation. Comparing with classical third-order finite volume WENO schemes Zhang and Shu (2009) [56] on tetrahedral meshes, the crucial advantages of such new multi-resolution WENO schemes are their simplicity and compactness with the application of only three unequal-sized central stencils for reconstructing unequal degree polynomials in the WENO type spatial procedures, their easy choice of arbitrary positive linear weights without considering the topology of the tetrahedral meshes, their optimal order of accuracy in smooth regions, and their suppression of spurious oscillations near strong discontinuities. The linear weights of such new multi-resolution WENO scheme can be any positive numbers on the condition that their sum is one. This is the first time that a series of unequal-sized hierarchical central spatial stencils are used in designing high-order finite volume WENO scheme on tetrahedral meshes. By performing such new spatial reconstruction procedures and adopting a third-order TVD Runge-Kutta method for time discretization, the occupied memory is decreased and the computing efficiency is increased. This new third-order finite volume multi-resolution WENO scheme is suitable for large scale engineering applications and could maintain good convergence property for steady-state problems on tetrahedral meshes. Benchmark examples are computed to demonstrate the robustness and good performance of these new finite volume WENO schemes.

Convergence to Steady-State Solutions of the New Type of High-Order Multi-resolution WENO Schemes: a Numerical Study

A new type of high-order multi-resolution weighted essentially non-oscillatory (WENO) schemes (Zhu and Shu in J Comput Phys, 375: 659–683, 2018) is applied to solve for steady-state problems on structured meshes. Since the classical WENO schemes (Jiang and Shu in J Comput Phys, 126: 202–228, 1996) might suffer from slight post-shock oscillations (which are responsible for the residue to hang at a truncation error level), this new type of high-order finite-difference and finite-volume multi-resolution WENO schemes is applied to control the slight post-shock oscillations and push the residue to settle down to machine zero in steady-state simulations. This new type of multi-resolution WENO schemes uses the same large stencils as that of the same order classical WENO schemes, could obtain fifth-order, seventh-order, and ninth-order in smooth regions, and could gradually degrade to first-order so as to suppress spurious oscillations near strong discontinuities. The linear weights of such new multi-resolution WENO schemes can be any positive numbers on the condition that their sum is one. This is the first time that a series of unequal-sized hierarchical central spatial stencils are used in designing high-order finite-difference and finite-volume WENO schemes for solving steady-state problems. In comparison with the classical fifth-order finite-difference and finite-volume WENO schemes, the residue of these new high-order multi-resolution WENO schemes can converge to a tiny number close to machine zero for some benchmark steady-state problems.

A simple, high-order and compact WENO limiter for RKDG method

In this paper, a new limiter using weighted essentially non-oscillatory (WENO) methodology is investigated for the Runge–Kutta discontinuous Galerkin (RKDG) methods for solving hyperbolic conservation laws. The idea is to use the high-order DG solution polynomial itself in the target cell and the linear polynomials which are reconstructed by the cell averages of solution in the target cell and its neighboring cells to reconstruct a new high-order polynomial in a manner of WENO methodology. Since only the linear polynomials need to be prepared for reconstruction, this limiter is very simple and compact with a stencil including only the target cell and its immediate neighboring cells. Numerical examples of various problems show that the new limiting procedure can simultaneously achieve uniform high-order accuracy and sharp, non-oscillatory shock transitions.

Completed repeated Richardson extrapolation for compressible fluid flows

Abstract Richardson extrapolation is a powerful approach for reducing the spatial discretization errors and increasing, in this way, the accuracy of the computed solution obtained by using many numerical methods for solving different scientific and engineering problems. This approach has been used in a variety of computational fluid dynamics problems to reduce the numerical error, but it has been mainly restricted to the computation of incompressible fluid flows and on grids with coincident nodes. The purpose of this work is to present a completed repeated Richardson extrapolation (CRRE) procedure for a more generic type of grids not necessarily with coincident nodes, and test it on compressible fluid flows. Three tests are performed for one- and quasi-one-dimensional Euler equations, i.e. (i) the Rayleigh flow, (ii) the isentropic flow, and (iii) the adiabatic flow through a nozzle. The last test presents a normal shock wave. In order to build a simple CRRE solver, these problems are solved using a first-order upwind-type finite difference method as the base scheme. The normal shock wave problem is also solved using a high-order weighted essentially non-oscillatory (WENO) scheme to compare it with the CRRE procedure. The procedure we have proposed can increase the achieved accuracy and significantly decrease the magnitude of the spatial error in all three tests. Its performance can be best demonstrated in the Rayleigh flow test, where the spatial discretization error is reduced by seven orders of magnitude and the achieved accuracy is increased from 0.998 to 6.62 on the grid with 10240 nodes. Similar performance can be observed for the isentropic flow, for which the spatial discretization error is reduced by nine orders of magnitude and the achieved accuracy is again increased from 0.996 to 6.73 on the grid with 10240 nodes. Finally, in the adiabatic flow with a normal shock wave, the procedure can reduce the spatial discretization error both upstream and downstream of the shock. We do remark, however, that the more expensive high-order WENO scheme has errors of lower magnitude upstream of the shock and has a sharper shock transition for this shocked test case.

Bounded and compact weighted essentially nonoscillatory limiters for discontinuous Galerkin schemes: Triangular elements

Two new classes of compact weighted essentially nonoscillatory (WENO) polynomial limiters are presented for second-, third-, fourth-, and fifth-order discontinuous Galerkin (DG) schemes on irregular simplex elements. The presented WENO-DG procedures are extensions of the high-order WENO finite-volume and finite-difference schemes of Zhu and Shu (2017) [25], (2019) [26] to high-order unstructured DG schemes. A compact positivity preserving limiter is applied to the solutions to ensure pressure and density remain within physical ranges at all time. It is then verified that the bounded WENO-DG maintains the formal order of accuracy of the underlying DG schemes in the smooth regions. The performance of the proposed WENO-DG is also demonstrated with inviscid test cases including the classical Riemann problems, shock-turbulence interaction, scramjet, blunt body flows, and the double Mach Reflection problems.

A new type of multi-resolution WENO schemes with increasingly higher order of accuracy on triangular meshes

In this paper, we continue our work in [46] and propose a new type of high-order finite volume multi-resolution weighted essentially non-oscillatory (WENO) schemes to solve hyperbolic conservation laws on triangular meshes. Although termed “multi-resolution WENO schemes”, we only use the information defined on a hierarchy of nested central spatial stencils and do not introduce any equivalent multi-resolution representation. We construct new third-order, fourth-order, and fifth-order WENO schemes using three or four unequal-sized central spatial stencils, different from the classical WENO procedure using equal-sized biased/central spatial stencils for the spatial reconstruction. The new WENO schemes could obtain the optimal order of accuracy in smooth regions, and could degrade gradually to first-order of accuracy so as to suppress spurious oscillations near strong discontinuities. This is the first time that only a series of unequal-sized hierarchical central spatial stencils are used in designing arbitrary high-order finite volume WENO schemes on triangular meshes. The main advantages of these schemes are their compactness, robustness, and their ability to maintain good convergence property for steady-state computation. The linear weights of such WENO schemes can be any positive numbers on the condition that they sum to one. Extensive numerical results are provided to illustrate the good performance of these new finite volume WENO schemes.

Adaptive mapping for high order WENO methods

In this paper, a novel mapping approach through the use of adaptive mapping functions is introduced for high order weighted essentially non-oscillatory (WENO) methods. The new class of adaptive mapping functions are designed to adjust themselves to the solution based on a simple parameter calculated using the smoothness indicators that are readily available during computation. It is shown that this adaptive nature allows the resultant mapped WENO scheme to maintain sub-stencil weights close to the optimal weights in smooth regions without amplifying the weights of non-smooth stencils containing discontinuities. Therefore, adaptive mapping achieves enhanced accuracy in smooth regions and is more resistant against spurious oscillations near discontinuities. Taylor series analysis of the seventh order finite volume WENO scheme has been performed to demonstrate the loss of accuracy of the original WENO method near critical points. The convergence rates of the seventh order finite volume WENO scheme with adaptive mapping have been shown through a simple numerical example. Excellent results have been obtained for one-dimensional linear advection cases especially over long output times. Improved results have also been obtained for one- and two-dimensional Euler equation test cases.