Search results

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.

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.

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.

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.

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.

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%.

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.

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.

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.

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.

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.

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.

The proposed EA approach for the RTE numerical solution effectively reduces the CPU time compared to the DC procedure with fixed relaxation coefficients.

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.

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.

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.

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.

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.

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.

Based on the normalized variable diagram, the weakness of the Gaskell and Lau's convective boundedness criterion (GL‐CBC) is revealed by numerical example. By careful consideration of the smoothness of the normalized variable variation pattern, more rigorous constraints on the interface value interpolation are found. A new CBC is thus proposed, whose feasibility and correctness are demonstrated by the inspection of ten existing bounded schemes and a numerical example.