Search results

1 – 10 of over 26000
Article
Publication date: 4 January 2021

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…

128

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

Engineering Computations, vol. 38 no. 6
Type: Research Article
ISSN: 0264-4401

Keywords

Article
Publication date: 1 March 1987

Antonio Strozzi

A class of exact solutions to the elastohydrodynamic problem to be used as test cases is presented. A numerical solution to the elastohydrodynamic problem according to the…

Abstract

A class of exact solutions to the elastohydrodynamic problem to be used as test cases is presented. A numerical solution to the elastohydrodynamic problem according to the Petrov—Galerkin method is developed. The appearance of spurious numerical undulations in the film profile is examined. A comparison between analytical and numerical results is employed to determine which numerical schemes limit the outcome of numerical oscillations without compromising the solution accuracy.

Details

Engineering Computations, vol. 4 no. 3
Type: Research Article
ISSN: 0264-4401

Article
Publication date: 1 June 1994

Michael M. Grigor’ev

The paper gives the description of boundary element method(BEM) with subdomains for the solution ofconvection—diffusion equations with variable coefficients and Burgers’equations…

Abstract

The paper gives the description of boundary element method (BEM) with subdomains for the solution of convection—diffusion equations with variable coefficients and Burgers’ equations. At first, the whole domain is discretized into K subdomains, in which linearization of equations by representing convective velocity by the sum of constant and variable parts is carried out. Then using fundamental solutions for convection—diffusion linear equations for each subdomain the boundary integral equation (in which the part of the convective term with the constant convective velocity is not included into the pseudo‐body force) is formulated. Only part of the convective term with the variable velocity, which is, as a rule, more than one order less than convective velocity constant part contribution, is left as the pseudo‐source. On the one hand, this does not disturb the numerical BEM—algorithm stability and, on the other hand, this leads to significant improvement in the accuracy of solution. The global matrix, similar to the case of finite element method, has block band structure whereas its width depends only on the numeration order of nodes and subdomains. It is noted, that in comparison with the direct boundary element method the number of global matrix non‐zero elements is not proportional to the square of the number of nodes, but only to the total number of nodal points. This allows us to use the BEM for the solution of problems with very fine space discretization. The proposed BEM with subdomains technique has been used for the numerical solution of one‐dimensional linear steady‐state convective—diffusion problem with variable coefficients and one‐dimensional non‐linear Burgers’ equation for which exact analytical solutions are available. It made it possible to find out the BEM correctness according to both time and space. High precision of the numerical method is noted. The good point of the BEM is the high iteration convergence, which is disturbed neither by high Reynolds numbers nor by the presence of negative velocity zones.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 4 no. 6
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 3 May 2016

Bantwal R. (Rabi) Baliga and Iurii Yuri Lokhmanets

The purpose of this paper is to present outcomes of efforts made over the last 20 years to extend the applicability of the Richardson extrapolation procedure to numerical

Abstract

Purpose

The purpose of this paper is to present outcomes of efforts made over the last 20 years to extend the applicability of the Richardson extrapolation procedure to numerical predictions of multidimensional, steady and unsteady, fluid flow and heat transfer phenomena in regular and irregular calculation domains.

Design/methodology/approach

Pattern-preserving grid-refinement strategies are proposed for mathematically rigorous generalizations of the Richardson extrapolation procedure for numerical predictions of steady fluid flow and heat transfer, using finite volume methods and structured multidimensional Cartesian grids; and control-volume finite element methods and unstructured two-dimensional planar grids, consisting of three-node triangular elements. Mathematically sound extrapolation procedures are also proposed for numerical solutions of unsteady and boundary-layer-type problems. The applicability of such procedures to numerical solutions of problems with curved boundaries and internal interfaces, and also those based on unstructured grids of general quadrilateral, tetrahedral, or hexahedral elements, is discussed.

Findings

Applications to three demonstration problems, with discretizations in the asymptotic regime, showed the following: the apparent orders of accuracy were the same as those of the numerical methods used; and the extrapolated results, measures of error, and a grid convergence index, could be obtained in a smooth and non-oscillatory manner.

Originality/value

Strict or approximate pattern-preserving grid-refinement strategies are used to propose generalized Richardson extrapolation procedures for estimating grid-independent numerical solutions. Such extrapolation procedures play an indispensable role in the verification and validation techniques that are employed to assess the accuracy of numerical predictions which are used for designing, optimizing, virtual prototyping, and certification of thermofluid systems.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 26 no. 3/4
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 3 April 2018

Omar Abu Arqub

The purpose of this study is to introduce the reproducing kernel algorithm for treating classes of time-fractional partial differential equations subject to Robin boundary…

416

Abstract

Purpose

The purpose of this study is to introduce the reproducing kernel algorithm for treating classes of time-fractional partial differential equations subject to Robin boundary conditions with parameters derivative arising in fluid flows, fluid dynamics, groundwater hydrology, conservation of energy, heat conduction and electric circuit.

Design/methodology/approach

The method provides appropriate representation of the solutions in convergent series formula with accurately computable components. This representation is given in the W(Ω) and H(Ω) inner product spaces, while the computation of the required grid points relies on the R(y,s) (x, t) and r(y,s) (x, t) reproducing kernel functions.

Findings

Numerical simulation with different order derivatives degree is done including linear and nonlinear terms that are acquired by interrupting the n-term of the exact solutions. Computational results showed that the proposed algorithm is competitive in terms of the quality of the solutions found and is very valid for solving such time-fractional models.

Research limitations/implications

Future work includes the application of the reproducing kernel algorithm to highly nonlinear time-fractional partial differential equations such as those arising in single and multiphase flows. The results will be published in forthcoming papers.

Practical implications

The study included a description of fundamental reproducing kernel algorithm and the concepts of convergence, and error behavior for the reproducing kernel algorithm solvers. Results obtained by the proposed algorithm are found to outperform in terms of accuracy, generality and applicability.

Social implications

Developing analytical and numerical methods for the solutions of time-fractional partial differential equations is a very important task owing to their practical interest.

Originality/value

This study, for the first time, presents reproducing kernel algorithm for obtaining the numerical solutions of some certain classes of Robin time-fractional partial differential equations. An efficient construction is provided to obtain the numerical solutions for the equations, along with an existence proof of the exact solutions based upon the reproducing kernel theory.

Article
Publication date: 2 November 2010

Mohamed Rady, Eric Arquis, Dominique Gobin and Benoît Goyeau

This paper aims to tackle the problem of thermo‐solutal convection and macrosegregation during ingot solidification of metal alloys. Complex flow structures associated with the…

Abstract

Purpose

This paper aims to tackle the problem of thermo‐solutal convection and macrosegregation during ingot solidification of metal alloys. Complex flow structures associated with the development of channels segregate and sharp gradients in the solutal field call for the implementation of accurate methods for numerical modeling of alloy solidification. In particular, the solute transport equation is convection dominated and requires special non‐oscillarity type high‐order schemes to handle the regions of channels segregates.

Design/methodology/approach

In the present study, a time‐splitting approach has been adopted to separately handle solute advection and diffusion. This splitting technique allows the application of accurate total variation dimensioning (TVD) schemes for solution of solute advection. Applications of second‐order Lax‐Wendroff TVD SUPERBEE and fifth‐order weighted essentially non‐oscillatory (WENO) schemes are described in the present article. Classical numerical solution of solute transport using hybrid and central‐difference schemes are also employed for the purpose of comparisons. Numerical simulations for solidification of Pb‐18%Sn in a two‐dimensional rectangular cavity have been carried out using different numerical schemes.

Findings

Numerical results show the difficulty of obtaining grid‐independent solutions with respect to local details in the region of channels. Grid convergence patterns and numerical uncertainty are found to be dependent on the applied scheme. In general, the first‐order hybrid scheme is diffusive and under predicts the formation of channels. The second‐order central‐difference scheme brings about oscillations with possible non‐physical extremes of solute composition in the region of channel segregates due to sharp gradients in the solutal field. The results obtained using TVD and WENO schemes contain no oscillations and show an excellent capture of channels formation and resolution of the interface between solute‐rich and depleted bands. Different stages of channels formation are followed by analyzing thermo‐solutal convection and macrosegregation at different times during solidification.

Research limitations/implications

Accurate prediction of local variation in the solutal and flow fields in the channels regions requires grid refinement up to scales in the order of microscopic dendrite arm spacing. This imposes limitations in terms of large computational time and applicability of available macroscopic models based on classical volume‐averaging techniques.

Practical implications

The present study is very useful for numerical simulation of macrosegregation during ingot casting of metal alloys.

Originality/value

The paper provides the methodology and application of TVD schemes to predict channel segregates during columnar solidification of metal alloys. It also demonstrates the limitations of classical schemes for simulation of alloy solidification.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 20 no. 8
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 26 August 2014

Anjali Verma, Ram Jiwari and Satish Kumar

The purpose of this paper is to propose a numerical scheme based on forward finite difference, quasi-linearisation process and polynomial differential quadrature method to find…

Abstract

Purpose

The purpose of this paper is to propose a numerical scheme based on forward finite difference, quasi-linearisation process and polynomial differential quadrature method to find the numerical solutions of nonlinear Klein-Gordon equation with Dirichlet and Neumann boundary condition.

Design/methodology/approach

In first step, time derivative is discretised by forward difference method. Then, quasi-linearisation process is used to tackle the non-linearity in the equation. Finally, fully discretisation by differential quadrature method (DQM) leads to a system of linear equations which is solved by Gauss-elimination method.

Findings

The accuracy of the proposed method is demonstrated by several test examples. The numerical results are found to be in good agreement with the exact solutions and the numerical solutions exist in literature. The proposed scheme can be expended for multidimensional problems.

Originality/value

The main advantage of the present scheme is that the scheme gives very accurate and similar results to the exact solutions by choosing less number of grid points. Secondly, the scheme gives better accuracy than (Dehghan and Shokri, 2009; Pekmen and Tezer-Sezgin, 2012) by choosing less number of grid points and big time step length. Also, the scheme can be extended for multidimensional problems.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 24 no. 7
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 16 October 2018

Rajni Rohila and R.C. Mittal

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

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 28 no. 11
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 3 August 2021

Zain ul Abdeen and Mujeeb ur Rehman

The purpose of this paper is to obtain a numerical scheme for finding numerical solutions of linear and nonlinear Hadamard-type fractional differential equations.

Abstract

Purpose

The purpose of this paper is to obtain a numerical scheme for finding numerical solutions of linear and nonlinear Hadamard-type fractional differential equations.

Design/methodology/approach

The aim of this paper is to develop a numerical scheme for numerical solutions of Hadamard-type fractional differential equations. The classical Haar wavelets are modified to align them with Hadamard-type operators. Operational matrices are derived and used to convert differential equations to systems of algebraic equations.

Findings

The upper bound for error is estimated. With the help of quasilinearization, nonlinear problems are converted to sequences of linear problems and operational matrices for modified Haar wavelets are used to get their numerical solution. Several numerical examples are presented to demonstrate the applicability and validity of the proposed method.

Originality/value

The numerical method is purposed for solving Hadamard-type fractional differential equations.

Details

Engineering Computations, vol. 39 no. 3
Type: Research Article
ISSN: 0264-4401

Keywords

Article
Publication date: 20 August 2021

Salam Adel Al-Bayati and Luiz C. Wrobel

The purpose of this paper is to describe an extension of the boundary element method (BEM) and the dual reciprocity boundary element method (DRBEM) formulations developed for one…

Abstract

Purpose

The purpose of this paper is to describe an extension of the boundary element method (BEM) and the dual reciprocity boundary element method (DRBEM) formulations developed for one- and two-dimensional steady-state problems, to analyse transient convection–diffusion problems associated with first-order chemical reaction.

Design/methodology/approach

The mathematical modelling has used a dual reciprocity approximation to transform the domain integrals arising in the transient equation into equivalent boundary integrals. The integral representation formula for the corresponding problem is obtained from the Green’s second identity, using the fundamental solution of the corresponding steady-state equation with constant coefficients. The finite difference method is used to simulate the time evolution procedure for solving the resulting system of equations. Three different radial basis functions have been successfully implemented to increase the accuracy of the solution and improving the rate of convergence.

Findings

The numerical results obtained demonstrate the excellent agreement with the analytical solutions to establish the validity of the proposed approach and to confirm its efficiency.

Originality/value

Finally, the proposed BEM and DRBEM numerical solutions have not displayed any artificial diffusion, oscillatory behaviour or damping of the wave front, as appears in other different numerical methods.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 32 no. 5
Type: Research Article
ISSN: 0961-5539

Keywords

1 – 10 of over 26000