Search results

The flux vector splitting (FVS) schemes are known for their higher resistance to shock instabilities and carbuncle phenomena in high-speed flow computations, which are generally accompanied by relatively large numerical diffusion. However, it is desirable to control the numerical diffusion of FVS schemes inside the boundary layer for improved accuracy in viscous flow computations. This study aims to develop a new methodology for controlling the numerical diffusion of FVS schemes for viscous flow computations with the help of a recently developed boundary layer sensor.

The governing equations are solved using a cell-centered finite volume approach and Euler time integration. The gradients in the viscous fluxes are evaluated by applying the Green’s theorem. For the inviscid fluxes, a new approach is introduced, where the original upwind formulation of an FVS scheme is first cast into an equivalent central discretization along with a numerical diffusion term. Subsequently, the numerical diffusion is scaled down by using a novel scaling function that operates based on a boundary layer sensor. The effectiveness of the approach is demonstrated by applying the same on van Leer’s FVS and AUSM schemes. The resulting schemes are named as Diffusion-Regulated van Leer’s FVS-Viscous (DRvLFV) and Diffusion-Regulated AUSM-Viscous (DRAUSMV) schemes.

The numerical tests show that the DRvLFV scheme shows significant improvement over its parent scheme in resolving the skin friction and wall heat flux profiles. The DRAUSMV scheme is also found marginally more accurate than its parent scheme. However, stability requirements limit the scaling down of only the numerical diffusion term corresponding to the acoustic part of the AUSM scheme.

To the best of the authors’ knowledge, this is the first successful attempt to regulate the numerical diffusion of FVS schemes inside boundary layers by applying a novel scaling function to their artificial viscosity forms. The new methodology can reduce the erroneous smearing of boundary layers by FVS schemes in high-speed flow applications.

This paper introduces a novel algorithm for solving the two‐dimensional Euler and Navier‐Stokes compressible equations using a one‐step effective flux vector‐splitting implicit method. The new approach makes a contribution by deriving a simple and yet effective implicit scheme which has the features of an exact factorization and avoids the solving of block‐diagonal system of equations. This results in a significant improvement in computational efficiency as compared to the standard Beam‐Warming and Steger implicit factored schemes. The current work has advantageous characteristics in the creation of higher order numerical implicit terms. The scheme is stable if we could select the correct values of the scalars (λ±ξ and λ±η) for the respective split flux‐vectors (F± and G±) along the ξ− and η−directions. A simple solving procedure is suggested with the discussion of the implicit boundary conditions, stability analysis, time‐step length and convergence criteria. This method is spatially second‐order accurate, fully conservative and implemented with general co‐ordinate transformations for treating complex geometries. Also, the scheme shows a good convergence rate and acceptable accuracy in capturing the shock waves. Results calculated from the program developed include transonic flows through convergence‐divergence nozzle and turbine cascade. Comparisons with other well‐documented experimental data are presented and their agreements are very promising. The extension of the algorithm to 3D simulation is straightforward and under way.

The purpose of this paper is to present two low diffusive convective-pressure flux split finite volume algorithms for solving incompressible flows in artificial compressibility framework.

The present method follows the framework similar to advection upwind splitting method of Liou and Steffen for compressible flows which is used by Vierendeels et al. to solve incompressible flow equations. Instead of discretizing the total inviscid flux using upwind scheme, the inviscid flux is first split into convective and pressure parts, and then discretized the two parts differently. The convective part is discretized using upwind method and the pressure part using central differencing. Since the Vierendeels type scheme may not be able to capture the divergence free velocity field due to the presence of artificial dissipation term, a strategy to progressively withdraw the dissipation with time step is proposed here that can ascertain the divergence free velocity condition to the level of residual error. This approach helps in reducing the amount of numerical dissipation due to upwind discretization, which is evident from the numerical test examples.

Upwind treatment of only the convective part of the inviscid flux terms, instead of the whole inviscid flux term, leads to more accurate solutions even at relatively coarse grids, which is substantiated by numerical test examples.

The method is presently applicable to Cartesian grid.

Although similar formulation is reported by Vierendeels et al., no detailed study of the accuracy is presented. Discretization and solution reconstructions used in the present approach differ from the approach reported by Vierendeels et al. A modification to Vierendeels type scheme is proposed that can help in achieving divergence free velocity condition. Finally the efficacy of the present approach to produce very accurate solutions even on coarse grids is successfully demonstrated using a few benchmark problems.

To develop a high‐order compact finite‐difference method for solving flow problems containing shock waves.

A numerical algorithm based on high‐order compact finite‐difference schemes is developed for solving Navier‐Stokes equations in two‐dimensional space. The convective flux terms are discretized by using advection upstream splitting method (AUSM). The developed method is then used to compute some example laminar flow problems. The problems considered have a range of Mach number that corresponds to subsonic incompressible flow to hypersonic compressible flows that contain shock waves and shock/boundary‐layer interaction.

The paper shows that the AUSM flux splitting and high‐order compact finite‐difference methods can be used accurately and robustly in resolving shear layers and capturing shock waves. The highly diffusive nature of conventional flux splitting especially on coarse grids makes them inaccurate for boundary layers even with high‐order discretization.

This paper presents a high‐order numerical method that can accurately and robustly capture shock waves without deteriorating oscillations and resolve boundary layers and shock/boundary layer interaction.

A class of flux‐splitting explicit second‐order finite difference schemes is set up. An ‘optimal’ scheme is defined for 1‐D flows and applied to 2‐D flows with CFL being able to reach 2. The results obtained show that this ‘optimal’ scheme is well adapted to the unsteady flows.

The purpose of this paper is to present a new hybrid Euler flux fonction for use in a finite-volume Euler/Navier-Stokes code and adapted to compressible flow problems.

The proposed scheme, called AUFSRR can be devised by combining the AUFS solver and the Roe solver, based on a rotated Riemann solver approach (Sun and Takayama, 2003; Ren, 2003). The upwind direction is determined by the velocity-difference vector and idea is to apply the AUFS solver in the direction normal to shocks to suppress carbuncle and the Roe solver across shear layers to avoid an excessive amount of dissipation. The resulting flux functions can be implemented in a very simple manner, in the form of the Roe solver with modified wave speeds, so that converting an existing AUFS flux code into the new fluxes is an extremely simple task.

The proposed flux functions require about 18 per cent more CPU time than the Roe flux. Accuracy, efficiency and other essential features of AUFSRR scheme are evaluated by analyzing shock propagation behaviours for both the steady and unsteady compressible flows. This is demonstrated by several test cases (1D and 2D) with standard finite-volume Euler code, by comparing results with existing methods.

The hybrid Euler flux function is used in a finite-volume Euler/Navier-Stokes code and adapted to compressible flow problems.

The AUFSRR scheme is devised by combining the AUFS solver and the Roe solver, based on a rotated Riemann solver approach.

A high‐order compact upwind algorithm is developed for solving Navier‐Stokes equations in two‐space dimensions. The method is based on advection upstream splitting method and fourth‐order compact finite‐difference schemes. The convection flux terms of the Navier‐Stokes equations are discretized by a compact cell‐centered differencing scheme while the diffusion flux terms are discretized by a central fourth‐order compact scheme. The midpoint values of the flux functions required by the cell‐centered compact scheme are determined by a fourth‐order MUSCL approach. For steady‐state solutions; first‐order implicit time integration, with LU decomposition, is employed. Computed results for a laminar flow past a flat plate and the problem of shock‐wave boundary layer interaction are presented.

An implicit method for the solution of transonic flows modelled by the time‐dependent Euler equations is presented. The method is characterized by a robust linearization for first‐ and second‐order versions of Roe's flux‐difference splitting scheme, an implicit treatment of the boundary conditions and the implementation of an adaptive grid strategy for global efficiency. The performance of the method is investigated for the GAMM test circular‐arc bump configuration and for the RAE 2822 aerofoil.

The purpose of this study is to present and demonstrate a numerical method for solving chemically reacting flows. These are important for energy conversion devices, which rely on chemical reactions as their operational mechanism, with heat generated from the combustion of the fuel, often gases, being converted to work.

The numerical study of such flows requires the set of Navier-Stokes equations to be extended to include multiple species and the chemical reactions between them. The numerical method implemented in this study also accounts for changes in the material properties because of temperature variations and the process to handle steep spatial fronts and stiff source terms without incurring any numerical instabilities. An all-speed numerical framework is used through simple low-dissipation advection upwind splitting (SLAU) convective scheme, and it has been extended in a multi-component species framework on the in-house density-based flow solver. The capability of solving turbulent combustion is also implemented using the Eddy Dissipation Concept (EDC) framework and the recent k-kl turbulence model.

The numerical implementation has been demonstrated for several stiff problems in laminar and turbulent combustion. The laminar combustion results are compared from the corresponding results from the Cantera library, and the turbulent combustion computations are found to be consistent with the experimental results.

This paper has extended the single gas density-based framework to handle multi-component gaseous mixtures. This paper has demonstrated the capability of the numerical framework for solving non-reacting/reacting laminar and turbulent flow problems. The all-speed SLAU convective scheme has been extended in the multi-component species framework, and the turbulent model k-kl is used for turbulent combustion, which has not been done previously. While the former method provides the capability of solving for low-speed flows using the density-based method, the later is a length-scale-based method that includes scale-adaptive simulation characteristics in the turbulence modeling. The SLAU scheme has proven to work well for unsteady flows while the k-kL model works well in non-stationary turbulent flows. As both these flow features are commonly found in industrially important reacting flows, the convection scheme and the turbulence model together will enhance the numerical predictions of such flows.

Two numerical methods, based on high order finite volume formulations and upwind schemes, are used to compute the two‐ and three‐dimensional flow field in a transonic nozzle. The influence of numerical diffusivity, boundary treatment and mesh structure is explored for inviscid and turbulent configurations. First order computations provide significantly different inviscid results. However, high order methods lead to similar solutions. An explanation of the error generated through the shockwave is proposed in this case. The two‐dimensional interaction of the shock with the thin turbulent boundary layer developing on the bump wall is also presented. Good agreement between both approaches is obtained considering the rapid thickening of the boundary layer due to the shock. Furthermore, the downstream velocity recovery is almost identical. Only slight discrepancies occur in the main flow, near the outer edge of the boundary layer. These seem to be related to the way the turbulence model deals with the free stream turbulence. Finally, preliminary three‐dimensional unstructured turbulent results are presented and discussed.