B.P. Leonard, A.P. Lock and M.K. Macvean
The NIRVANA project is concerned with the development of anonoscillatory, integrally reconstructed,volume‐averaged numerical advectionscheme. The conservative, flux‐based…
Abstract
The NIRVANA project is concerned with the development of a nonoscillatory, integrally reconstructed, volume‐averaged numerical advection scheme. The conservative, flux‐based finite‐volume algorithm is built on an explicit, single‐step, forward‐in‐time update of the cell‐average variable, without restrictions on the size of the time‐step. There are similarities with semi‐Lagrangian schemes; a major difference is the introduction of a discrete integral variable, guaranteeing conservation. The crucial step is the interpolation of this variable, which is used in the calculation of the fluxes; the (analytic) derivative of the interpolant then gives sub‐cell behaviour of the advected variable. In this paper, basic principles are described, using the simplest possible conditions: pure one‐dimensional advection at constant velocity on a uniform grid. Piecewise Nth‐degree polynomial interpolation of the discrete integral variable leads to an Nth‐order advection scheme, in both space and time. Nonoscillatory results correspond to convexity preservation in the integrated variable, leading naturally to a large‐Δt generalisation of the universal limited. More restrictive TVD constraints are also extended to large Δt. Automatic compressive enhancement of step‐like profiles can be achieved without exciting “stair‐casing”. One‐dimensional simulations are shown for a number of different interpolations. In particular, convexity‐limited cubic‐spline and higher‐order polynomial schemes give very sharp, nonoscillatory results at any Courant number, without clipping of extrema. Some practical generalisations are briefly discussed.
Details
Keywords
B.P. LEONARD and SIMIN MOKHTARI
In 1982, Smith and Hutton published comparative results of several different convection‐diffusion schemes applied to a specially devised test problem involving…
Abstract
In 1982, Smith and Hutton published comparative results of several different convection‐diffusion schemes applied to a specially devised test problem involving near‐discontinuities and strong streamline curvature. First‐order methods showed significant artificial diffusion, whereas higher‐order methods gave less smearing but had a tendency to overshoot and oscillate. Perhaps because unphysical oscillations are more obvious than unphysical smearing, the intervening period has seen a rise in popularity of low‐order artificially diffusive schemes, especially in the numerical heat‐transfer industry. This paper presents an alternative strategy of using non‐artificially diffusive higher‐order methods, while maintaining strictly monotonic transitions through the use of simple flux‐limiter constraints. Limited third‐order upwinding is usually found to be the most cost‐effective basic convection scheme. Tighter resolution of discontinuities can be obtained at little additional cost by using automatic adaptive stencil expansion to higher order in local regions, as needed.
Details
Keywords
A.N. Pavlov, S.S. Sazhin, R.P. Fedorenko and M.R. Heikal
Detailed results of numerical calculations of transient, 2D incompressible flow around and in the wake of a square prism at Re = 100, 200 and 500 are presented. An implicit…
Abstract
Detailed results of numerical calculations of transient, 2D incompressible flow around and in the wake of a square prism at Re = 100, 200 and 500 are presented. An implicit finite‐difference operator‐splitting method, a version of the known SIMPLEC‐like method on a staggered grid, is described. Appropriate theoretical results are presented. The method has second‐order accuracy in space, conserving mass, momentum and kinetic energy. A new modification of the multigrid method is employed to solve the elliptic pressure problem. Calculations are performed on a sequence of spatial grids with up to 401 × 321 grid points, at sequentially halved time steps to ensure grid‐independent results. Three types of flow are shown to exist at Re = 500: a steady‐state unstable flow and two which are transient, fully periodic and asymmetric about the centre line but mirror symmetric to each other. Discrete frequency spectra of drag and lift coefficients are presented.
Details
Keywords
Carlos Enrique Torres-Aguilar, Jesús Xamán, Pedro Moreno-Bernal, Iván Hernández-Pérez, Ivett Zavala-Guillén and Irving Osiris Hernández-López
The purpose of this study is to propose a novel relaxation modified factor to accelerate the numerical solution of the radiative transfer equation (RTE) with several…
Abstract
Purpose
The purpose of this study is to propose a novel relaxation modified factor to accelerate the numerical solution of the radiative transfer equation (RTE) with several high-resolution total variation diminishing schemes. The methodology proposed is denoted as the X-factor method.
Design/methodology/approach
The X-factor method was compared with the technique deferred-correction (DC) for the calculations of a two-dimensional cavity with absorting-emiting-scatteting gray media using the discrete ordinates method. Four parameters were considered to evaluate: the absorption coefficient, the emissivity of boundary surface, the scattering albedo and under-relaxation factor.
Findings
The results showed the central processing unit (CPU) time of X-factor method was lower than DC. The reductions of CPU time with the X-factor method were observed from 0.6 to 75.4%.
Originality/value
The superiority of the X-factor method over DC was showed with the reduction of CPU time of the numerical solution of RTE for evaluated cases.
Details
Keywords
V. PENNATI, M. MARELLI and L.M. DE BIASE
In this paper new cubic v‐splines monotonic one‐dimensional profiles are presented, for the finite volume solution of convection‐diffusion problems. By studying the profile in…
Abstract
In this paper new cubic v‐splines monotonic one‐dimensional profiles are presented, for the finite volume solution of convection‐diffusion problems. By studying the profile in normalized variables, some weight functions have been determined for the profile. Being free of the requirement that the volumes be equal, the volume size can be reduced where needed. Numerical properties of the proposed method were formally analysed and are confirmed by numerical examples included here.
Details
Keywords
A. Pascau, C. Pérez and D. Sánchez
A new discretization scheme named NOTABLE (New Option forthe Treatment of Advection in the Boundary Layer Equations) ispresented. Despite its name, this scheme is intended to be…
Abstract
A new discretization scheme named NOTABLE (New Option for the Treatment of Advection in the Boundary Layer Equations) is presented. Despite its name, this scheme is intended to be used in a general transport equation to discretize the convective term. It is formally third‐order accurate in regions of smooth solution and first‐order accurate at grid points having local maxima. Within the finite‐volume formulation it relates the face values to the nodal values via a non‐linear function. This scheme has been compared with well‐known high‐order schemes like QUICK and it has always given more accurate solutions. After assessing the scheme in several unidimensional test cases for which an exact solution is available, its performance in a complex swirling flow is addressed.
Details
Keywords
Carlos Enrique Torres-Aguilar, Pedro Moreno-Bernal, Jesús Xamán, Ivett Zavala Guillen and Irving Osiris Hernández-López
This paper aims to present an evolutionary algorithm (EA) to accelerate the convergence for the radiative transfer equation (RTE) numerical solution using high-order and…
Abstract
Purpose
This paper aims to present an evolutionary algorithm (EA) to accelerate the convergence for the radiative transfer equation (RTE) numerical solution using high-order and high-resolution schemes by the relaxation coefficients optimization.
Design methodology/approach
The objective function minimizes the residual value difference between iterations in each control volume until its difference is lower than the convergence criterion. The EA approach is evaluated in two configurations, a two-dimensional cavity with scattering media and absorbing media.
Findings
Experimental results show the capacity to obtain the numerical solution for both cases on all interpolation schemes tested by the EA approach. The EA approach reduces CPU time for the RTE numerical solution using SUPERBEE, SWEBY and MUSCL schemes until 97% and 135% in scattering and absorbing media cases, respectively. The relaxation coefficients optimized every two numerical solution iterations achieve a significant reduction of the CPU time compared to the deferred correction procedure with fixed relaxation coefficients.
Originality/value
The proposed EA approach for the RTE numerical solution effectively reduces the CPU time compared to the DC procedure with fixed relaxation coefficients.
Details
Keywords
Tian-Yu Wu, Jianfei Zhang, Yanjun Dai, Tao-Feng Cao, Kong Ling and Wen-Quan Tao
To present the detailed implementation processes of the IDEAL algorithm for two-dimensional compressible flows based on Delaunay triangular mesh, and compare the performance of…
Abstract
Purpose
To present the detailed implementation processes of the IDEAL algorithm for two-dimensional compressible flows based on Delaunay triangular mesh, and compare the performance of the SIMPLE and IDEAL algorithms for solving compressible problems. What’s more, the implementation processes of Delaunay mesh generation and derivation of the pressure correction equation are also introduced.
Design/methodology/approach
Programming completely in C++.
Findings
Five compressible examples are used to test the SIMPLE and IDEAL algorithms, and the comparison with measurement data shows good agreement. The IDEAL algorithm has much better performance in both convergence rate and stability over the SIMPLE algorithm.
Originality/value
The detail solution procedure of implementing the IDEAL algorithm for compressible flows based on Delaunay triangular mesh is presented in this work, seemingly first in the literature.
Details
Keywords
A finite difference scheme for convection term discretization, calledBSOU (stands for Bounded Second Order Upwind), is developed and itsperformance is assessed against exact or…
Abstract
A finite difference scheme for convection term discretization, called BSOU (stands for Bounded Second Order Upwind), is developed and its performance is assessed against exact or benchmark solutions in linear and non‐linear cases. It employs a flux blending technique between first order upwind and second order upwind schemes only in those regions of the flow field where spurious oscillations are likely to occur. The blending factors are calculated with the aid of the convection boundedness criterion. In all cases the scheme performed very well, minimizing the numerical diffusion errors. The scheme is transportive, conservative, bounded, stable and accurate enough so as to be suitable for inclusion into a general purpose solution algorithm.
Details
Keywords
E. Jahanbakhsh, R. Panahi and M.S. Seif
This study aims to present compatible computational fluid dynamics procedure for calculation of incompressible three‐dimensional time‐dependent flow with complicated free surface…
Abstract
Purpose
This study aims to present compatible computational fluid dynamics procedure for calculation of incompressible three‐dimensional time‐dependent flow with complicated free surface deformation. A computer software is developed and validated using a variety of academic test cases.
Design/methodology/approach
Two fluids are modeled as a single continuum with a fluid property jump at the interface by solving a scalar transport equation for volume fraction. In conjunction, the conservation equations for mass and momentum are solved using fractional step method. Here, a finite volume discretisation and colocated arrangement are used.
Findings
The developed code results in accurate simulation of interfacial flows, e.g. Rayleigh‐Taylor instability, sloshing and dambreaking problems. All results are in good concordance with experimental data especially when there are two phases with high density ratio.
Research limitations/implications
Turbulence, which has great importance in a wide variety of real world phenomena, is not considered in the present formulation and left for future researches.
Originality/value
Here, an integrated numerical simulation for transient interfacial flows is presented. In this way, the pressure integral term in Navier‐Stokes equation is discretised based on a newly developed interpolation which results in non‐oscillative velocity field especially in free surface.