Search results

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.

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.

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 10³-10⁸, 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.

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.

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.

The purpose of this paper is to analyse the force convection laminar boundary layer flow on irregular boundary in diverging channel with the presence of magnetic field effects. Effects of various fluid parameters such as suction/injection, viscous dissipation, magnetic parameter and heat source/sink on velocity and temperature profiles are numerically analyzed. Moreover, numerically investigated on skin-friction and heat transfer coefficients when suction/injection occur.

The governing coupled partial differential equations are transformed to dimensionless form using non-similarity transformations. The non-dimensional partial differential equations are linearized by quasi-linearization technique and solved by varga's algorithm with numerical finite difference scheme on a non-uniform mesh.

The computation results are presented in terms of temperature, heat transfer and skin friction coefficients; these are useful for determining surface heat requirements. It was found that, in finite difference scheme for non-uniform mesh with quasi-linearization technique method gives smoothness of solution compared to finite difference scheme for uniform mesh, and this evidence is graphically represented in Figure 2.

The impacts of viscous dissipation (Ec) and magnetic parameter (Ha) on temperature profiles, skin friction and heat transfer are analyzed, which determine the heat generation/absorption to ensure the MHD flow of the laminar boundary layer on irregular boundary over a diverging channel.

The purpose of this paper is to develop an efficient numerical algorithm for the self‐consistent solution of Schrodinger and Poisson equations in one‐dimensional systems. The goal is to compute the charge‐control and capacitance‐voltage characteristics of quantum wire transistors.

The paper presents a numerical formulation employing a non‐uniform finite difference discretization scheme, in which the wavefunctions and electronic energy levels are obtained by solving the Schrödinger equation through the split‐operator method while a relaxation method in the FTCS scheme (“Forward Time Centered Space”) is used to solve the two‐dimensional Poisson equation.

The numerical model is validated by taking previously published results as a benchmark and then applying them to yield the charge‐control characteristics and the capacitance‐voltage relationship for a split‐gate quantum wire device.

The paper helps to fulfill the need for C‐V models of quantum wire device. To do so, the authors implemented a straightforward calculation method for the two‐dimensional electronic carrier density n(x,y). The formulation reduces the computational procedure to a much simpler problem, similar to the one‐dimensional quantization case, significantly diminishing running time.

This paper aims to propose a new adaptive numerical method to find more accurate numerical solution for the heat source optimal control problem (OCP).

The main aim of this paper is to present an adaptive collocation approach based on the interpolating wavelets to solve an OCP for finding optimal heat source, in a two-dimensional domain. This problem arises when the domain is heated by microwaves or by electromagnetic induction.

This paper shows that combination of interpolating wavelet basis and finite difference method makes an accurate structure to design adaptive algorithm for such problems which usually have non-smooth solution.

The proposed numerical technique is flexible for different OCP governed by a partial differential equation with box constraint over the control or the state function.

The space fractional PDEs (SFPDEs) play an important role in the fractional calculus field. Proposing a high-order, stable and flexible numerical procedure for solving SFPDEs is the main aim of most researchers. This paper devotes to developing a novel spectral algorithm to solve the FitzHugh–Nagumo models with space fractional derivatives.

The fractional derivative is defined based upon the Riesz derivative. First, a second-order finite difference formulation is used to approximate the time derivative. Then, the Jacobi spectral collocation method is employed to discrete the spatial variables. On the other hand, authors assume that the approximate solution is a linear combination of special polynomials which are obtained from the Jacobi polynomials, and also there exists Riesz fractional derivative based on the Jacobi polynomials. Also, a reduced order plan, such as proper orthogonal decomposition (POD) method, has been utilized.

A fast high-order numerical method to decrease the elapsed CPU time has been constructed for solving systems of space fractional PDEs.

The spectral collocation method is combined with the POD idea to solve the system of space-fractional PDEs. The numerical results are acceptable and efficient for the main mathematical model.

The centrifugal instability mechanism of boundary layers over concave surfaces is responsible for the development of quasi-periodic, counter-rotating vortices aligned in a streamwise direction known as Görtler vortices. By distorting the boundary layer structure in both the spanwise and the wall-normal directions, Görtler vortices may modify heat transfer rates. The purpose of this study is to conduct spatial numerical simulation experiments based on a vorticity–velocity formulation of the incompressible Navier–Stokes system of equations to quantify the role of the transition in the heat transfer process.

Experiments are conducted using an in-house, parallel, message-passing code. Compact finite difference approximations and a spectral method are used to approximate spatial derivatives. A fourth-order Runge–Kutta method is adopted for time integration. The Poisson equation is solved using a geometric multigrid method.

Results show that the numerical method can capture the physics of transitional flows over concave geometries. They also show that the heat transfer rates in the late stages of the transition may be greater than those for either laminar or turbulent ones.

The numerical method can be considered as a robust alternative to investigate heat transfer properties in transitional boundary layer flows over concave surfaces.

Variable density flows play an important role in many technological devices and natural phenomena. The purpose of this paper is to develop a robust and accurate method for low Mach number flows with large density and temperature variations.

Low Mach number approximation approach is used in the paper combined with a predictor-corrector method and accurate compact scheme of fourth and sixth order. A novel algorithm is formulated for the projection method in which the boundary conditions for the pressure are implemented in such a way that the continuity equation is fulfilled everywhere in the computational domain, including the boundary nodes.

It is shown that proposed implementation of the boundary conditions considerably improves a solution accuracy. Assessment of the accuracy was performed based on the constant density Burggraf flow and for two benchmark cases for the natural convection problems: steady flow in a square cavity and unsteady flow in a tall cavity. In all the cases the results agree very well with exemplary solutions.

A staggered or half-staggered grid arrangement is usually used for the projection method for both constant and low Mach number flows. The staggered approach ensures stability and strong pressure-velocity coupling. In the paper a high-order compact method has been implemented in the framework of low Mach number approximation on collocated meshes. The resulting algorithm is accurate, robust for large density variations and is almost free from the pressure oscillations.

Mixed recirculatory flow in the annuli of stationary and rotating horizontal cylinders were studied numerically. A set of distorted ‘false transient’ parameters were introduced to speed up the steady state solution of the unsteady vorticity, energy and stream function—vorticity equations. The inner cylinder of the annuli is assumed heated and rotating at Reynolds numbers that exclude the effects of centrifugal acceleration and three‐dimensional Taylor vortices. The Prandtl number considered is in the range of 0.01 to 1.0 and Rayleigh number in the range of 10² to 10⁶. Radius ratios of the cylinders considered are 1.25, 2.5 and 5.0. For a radius ratio of 2.5, inner cylinder rotation in the Reynolds number range of 0 to 1120 was considered. Vertical eccentricities in the range of ±2/3 were studied for cases of the rotating inner cylinder. Numerical experiments show that the mean Nusselt number increases with Rayleigh number for both cases of concentric and eccentric stationary inner cylinder. At a Prandtl number of order 1.0 with a fixed Rayleigh number, when the inner cylinder is made to rotate, the mean Nusselt number decreases throughout the flow. At lower Prandtl number of the order 0.1 to 0.01, the mean Nusselt number remained fairly constant with respect to the rotational Reynolds number.

Large‐eddy simulation (LES) of a turbulent channel flow is performed using different subfilter‐scale (SFS) models and test filter functions. The SFS models used are the dynamic Smagorinsky model (DSM) and the dynamic mixed model (DMM). The DMM is a linear combination between the scale‐similarity model and the DSM. The test filter functions investigated are the sharp cut‐off (in spectral space) and smooth filter that is commutative up to fourth‐order. The filters are applied either in the homogeneous directions or in all three spatial directions. The governing equations are discretized using a fourth‐order energy‐conserving finite‐difference scheme. The influence from the test filter function and the SFS model on the LES results are investigated and the effect of two‐dimensional versus three‐dimensional test filtering are investigated. The study shows that the combination of SFS model and filter function highly influences the computational results; even the effect on the zeroth‐order moment is large.

The purpose of this paper is to study the effect of non‐uniform double slot suction (injection) into a steady laminar boundary layer flow over a yawed cylinder when fluid properties such as viscosity and Prandtl number are inverse linear functions of temperature. Non‐similar solutions have been obtained from the starting point of the streamwise co‐ordinate to the exact point of separation.

The governing equations are tackled by the implicit finite difference scheme in combination with the quasi‐linearization technique. Quasi‐linear technique can be viewed as a generalization of the Newton‐Raphson approximation technique in functional space. An iterative sequence of linear equations is carefully constructed to approximate the nonlinear equations for achieving quadratic convergence and monotonicity. The quadratic convergence and monotonicity are unique characteristics of the quasilinear implicit finite difference scheme, which makes this scheme superior to built‐in iteration of upwind or finite amplitude techniques.

The results indicate that the separation can be delayed by non‐uniform double slot suction and also by moving the slot downstream. However, the effect of non‐uniform double slot injection is just the opposite. Yaw angle has very little affect on the location of the point of separation.

This analysis is useful in understanding many boundary layer problems of practical importance for undersea applications, for example, in suppressing recirculating bubbles and controlling transition and/or separation of the boundary layer over submerged bodies.