Considers the extension of a new class of higher‐order compact methods to nonuniform grids and examines the effect of pollution that arises with differencing the associated metric…
Abstract
Considers the extension of a new class of higher‐order compact methods to nonuniform grids and examines the effect of pollution that arises with differencing the associated metric coefficients. Numerical studies for the standard model convection diffusion equation in 1D and 2D are carried out to validate the convergence behaviour and demonstrate the high‐order accuracy.
Details
Keywords
Nagesh Babu Balam and Akhilesh Gupta
Modelling accurately the transient behaviour of natural convection flow in enclosures been a challenging task because of a variety of numerical errors which have limited achieving…
Abstract
Purpose
Modelling accurately the transient behaviour of natural convection flow in enclosures been a challenging task because of a variety of numerical errors which have limited achieving the higher order temporal accuracy. A fourth-order accurate finite difference method in both space and time is proposed to overcome these numerical errors and accurately model the transient behaviour of natural convection flow in enclosures using vorticity–streamfunction formulation.
Design/methodology/approach
Fourth-order wide stencil formula with appropriate one-sided difference extrapolation technique near the boundary is used for spatial discretisation, and classical fourth-order Runge–Kutta scheme is applied for transient term discretisation. The proposed method is applied on two transient case studies, i.e. convection–diffusion of a Gaussian Pulse and Taylor Vortex flow having analytical solution.
Findings
Error magnitude comparison and rate of convergence analysis of the proposed method with these analytical solutions establish fourth-order accuracy and prove the ability of the proposed method to truly capture the transient behaviour of incompressible flow. Also, to test the transient natural convection flow behaviour, the algorithm is tested on differentially heated square cavity at high Rayleigh number in the range of 103-108, followed by studying the transient periodic behaviour in a differentially heated vertical cavity of aspect ratio 8:1. An excellent comparison is obtained with standard benchmark results.
Research limitations/implications
The developed method is applied on 2D enclosures; however, the present methodology can be extended to 3D enclosures using velocity–vorticity formulations which shall be explored in future.
Originality/value
The proposed methodology to achieve fourth-order accurate transient simulation of natural convection flows is novel, to the best of the authors’ knowledge. Stable fourth-order vorticity boundary conditions are derived for boundary and external boundary regions. The selected case studies for comparison demonstrate not only the fourth-order accuracy but also the considerable reduction in error magnitude by increasing the temporal accuracy. Also, this study provides novel benchmark results at five different locations within the differentially heated vertical cavity of aspect ratio 8:1 for future comparison studies.
Details
Keywords
Abstract
Purpose
The purpose of this paper is to develop a Vortex-In-Cell (VIC) method with the semi-Lagrangian scheme and apply it to the high-Re lid-driven cavity flow.
Design/methodology/approach
The VIC method is developed for simulating high Reynolds number incompressible flow. A semi-Lagrangian scheme is incorporated in the convection term to produce unconditional stability, which gets rid of the constraint of the convection Courant-Friedrichs-Lewy (CFL) condition; the adaptive time step is used to maintain the numerical stability of the diffusion term; and the velocity boundary condition is readily converted to the vorticity formulation to suit discontinuous boundary treatment. The VIC simulation results are compared with those produced by other gird methods reported in open literature studies.
Findings
The lid-driven cavity flow is simulated from Re = 100 to 100,000. Similar vortex birth mechanisms are exhibited though, but distinct flow characteristics are revealed. At Re = 100 to 7,500, the cavity flow is confirmed steady. At Re = 10,000, 15,000 and 20,000, the cavity flow is periodical with a primary vortex held spatially at the center. In particular, at Re = 100,000 highly turbulent characteristics is first revealed and an analogous primary vortex is formed but in motion rather than stationary, which is caused by the considerable flow separation at all the boundaries.
Originality/value
In the lid-driven cavity, the flow becomes extremely complex and highly turbulent at Re = 100,000, and the analogous primary vortex structure is observed. Boundary layer separation is observed at all walls, producing small vortices and causing the displacement of the analogous primary vortex. Such a finding original and has not yet been reported by other investigators. It may provide a basis for conducting in-depth studies of the lid-driven cavity flow.
Details
Keywords
Ramesh Chand Mittal, Sudhir Kumar and Ram Jiwari
The purpose of this study is to extend the cubic B-spline quasi-interpolation (CBSQI) method via Kronecker product for solving 2D unsteady advection-diffusion equation. The CBSQI…
Abstract
Purpose
The purpose of this study is to extend the cubic B-spline quasi-interpolation (CBSQI) method via Kronecker product for solving 2D unsteady advection-diffusion equation. The CBSQI method has been used for solving 1D problems in literature so far. This study seeks to use the idea of a Kronecker product to extend the method for 2D problems.
Design/methodology/approach
In this work, a CBSQI is used to approximate the spatial partial derivatives of the dependent variable. The idea of the Kronecker product is used to extend the method for 2D problems. This produces the system of ordinary differential equations (ODE) with initial conditions. The obtained system of ODE is solved by strong stability preserving the Runge–Kutta method (SSP-RK-43).
Findings
It is found that solutions obtained by the proposed method are in good agreement with the analytical solution. Further, the results are also compared with available numerical results in the literature, and a reasonable degree of compliance is observed.
Originality/value
To the best of the authors’ knowledge, the CBSQI method is used for the first time for solving 2D problems and can be extended for higher-dimensional problems.
Details
Keywords
The purpose of this paper is to determine both analytically and numerically the existence of smooth, cusped and sharp shock wave solutions to a one-dimensional model of…
Abstract
Purpose
The purpose of this paper is to determine both analytically and numerically the existence of smooth, cusped and sharp shock wave solutions to a one-dimensional model of microfluidic droplet ensembles, water flow in unsaturated flows, infiltration, etc., as functions of the powers of the convection and diffusion fluxes and upstream boundary condition; to study numerically the evolution of the wave for two different initial conditions; and to assess the accuracy of several finite difference methods for the solution of the degenerate, nonlinear, advection--diffusion equation that governs the model.
Design/methodology/approach
The theory of ordinary differential equations and several explicit, finite difference methods that use first- and second-order, accurate upwind, central and compact discretizations for the convection terms are used to determine the analytical solution for steadily propagating waves and the evolution of the wave fronts from hyperbolic tangent and piecewise linear initial conditions to steadily propagating waves, respectively. The amplitude and phase errors of the semi-discrete schemes are determined analytically and the accuracy of the discrete methods is assessed.
Findings
For non-zero upstream boundary conditions, it has been found both analytically and numerically that the shock wave is smooth and its steepness increases as the power of the diffusion term is increased and as the upstream boundary value is decreased. For zero upstream boundary conditions, smooth, cusped and sharp shock waves may be encountered depending on the powers of the convection and diffusion terms. For a linear diffusion flux, the shock wave is smooth, whereas, for a quadratic diffusion flux, the wave exhibits a cusped front whose left spatial derivative decreases as the power of the convection term is increased. For higher nonlinear diffusion fluxes, a sharp shock wave is observed. The wave speed decreases as the powers of both the convection and the diffusion terms are increased. The evolution of the solution from hyperbolic tangent and piecewise linear initial conditions shows that the wave back adapts rapidly to its final steady value, whereas the wave front takes much longer, especially for piecewise linear initial conditions, but the steady wave profile and speed are independent of the initial conditions. It is also shown that discretization of the nonlinear diffusion flux plays a more important role in the accuracy of first- and second-order upwind discretizations of the convection term than either a conservative or a non-conservative discretization of the latter. Second-order upwind and compact discretizations of the convection terms are shown to exhibit oscillations at the foot of the wave’s front where the solution is nil but its left spatial derivative is largest. The results obtained with a conservative, centered second--order accurate finite difference method are found to be in good agreement with those of the second-order accurate, central-upwind Kurganov--Tadmor method which is a non-oscillatory high-resolution shock-capturing procedure, but differ greatly from those obtained with a non-conservative, centered, second-order accurate scheme, where the gradients are largest.
Originality/value
A new, one-dimensional model for microfluidic droplet transport, water flow in unsaturated flows, infiltration, etc., that includes high-order convection fluxes and degenerate diffusion, is proposed and studied both analytically and numerically. Its smooth, cusped and sharp shock wave solutions have been determined analytically as functions of the powers of the nonlinear convection and diffusion fluxes and the boundary conditions. These solutions are used to assess the accuracy of several finite difference methods that use different orders of accuracy in space, and different discretizations of the convection and diffusion fluxes, and can be used to assess the accuracy of other numerical procedures for one-dimensional, degenerate, convection--diffusion equations.
Details
Keywords
This imposing volume is the third in order of publication of the Princeton series. Its eight hundred pages comprise sections on: Fundamentals of Thermodynamics by Rossini;…
Abstract
This imposing volume is the third in order of publication of the Princeton series. Its eight hundred pages comprise sections on: Fundamentals of Thermodynamics by Rossini; Fundamental Physics of Gases by Herzfeld, Griffing, Hirschfelder, Curtiss, Bird and Spotz; Thermodynamic Properties of Real Gases and their Mixtures by Beattie; Transport Properties of Gases and Gaseous Mixtures by Hirschfelder, Curtiss, Bird and Spotz; Critical Phenomena by Rice; Properties of Liquids and Liquid Solutions by Richardson and Brinkley; Properties of Solids and Solid Solutions by Ewald; Relaxation Phenomena in Gases by Herzfeld; Gases at Low Densities by Estermann; and Thermodynamics of Irreversible Processes by Curtiss.
This paper aims to develop a novel numerical method based on bi-cubic B-spline functions and alternating direction (ADI) scheme to study numerical solutions of advection diffusion…
Abstract
Purpose
This paper aims to develop a novel numerical method based on bi-cubic B-spline functions and alternating direction (ADI) scheme to study numerical solutions of advection diffusion equation. The method captures important properties in the advection of fluids very efficiently. C.P.U. time has been shown to be very less as compared with other numerical schemes. Problems of great practical importance have been simulated through the proposed numerical scheme to test the efficiency and applicability of method.
Design/methodology/approach
A bi-cubic B-spline ADI method has been proposed to capture many complex properties in the advection of fluids.
Findings
Bi-cubic B-spline ADI technique to investigate numerical solutions of partial differential equations has been studied. Presented numerical procedure has been applied to important two-dimensional advection diffusion equations. Computed results are efficient and reliable, have been depicted by graphs and several contour forms and confirm the accuracy of the applied technique. Stability analysis has been performed by von Neumann method and the proposed method is shown to satisfy stability criteria unconditionally. In future, the authors aim to extend this study by applying more complex partial differential equations. Though the structure of the method seems to be little complex, the method has the advantage of using small processing time. Consequently, the method may be used to find solutions at higher time levels also.
Originality/value
ADI technique has never been applied with bi-cubic B-spline functions for numerical solutions of partial differential equations.
Details
Keywords
Vitus Mwinteribo Tabie, Chong Li, Wang Saifu, Jianwei Li and Xiaojing Xu
This paper aims to present a broad review of near-a titanium alloys for high-temperature applications.
Abstract
Purpose
This paper aims to present a broad review of near-a titanium alloys for high-temperature applications.
Design/methodology/approach
Following a brief introduction of titanium (Ti) alloys, this paper considers the near-α group of Ti alloys, which are the most popular high-temperature Ti alloys developed for a high-temperature application, particularly in compressor disc and blades in aero-engines. The paper is relied on literature within the past decade to discuss phase stability and microstructural effect of alloying elements, plastic deformation and reinforcements used in the development of these alloys.
Findings
The near-a Ti alloys show high potential for high-temperature applications, and many researchers have explored the incorporation of TiC, TiB SiC, Y2O3, La2O3 and Al2O3 reinforcements for improved mechanical properties. Rolling, extrusion, forging and some severe plastic deformation (SPD) techniques, as well as heat treatment methods, have also been explored extensively. There is, however, a paucity of information on SiC, Y2O3 and carbon nanotube reinforcements and their combinations for improved mechanical properties. Information on some SPD techniques such as cyclic extrusion compression, multiaxial compression/forging and repeated corrugation and straightening for this class of alloys is also limited.
Originality/value
This paper provides a topical, technical insight into developments in near-a Ti alloys using literature from within the past decade. It also outlines the future developments of this class of Ti alloys.