Search results
1 – 10 of 610An unconditionally positive definite finite difference scheme termed as UPFD has been derived to approximate a linear advection-diffusion-reaction equation which models…
Abstract
Purpose
An unconditionally positive definite finite difference scheme termed as UPFD has been derived to approximate a linear advection-diffusion-reaction equation which models exponential travelling waves and the coefficients of advection, diffusion and reactive terms have been chosen as one (Chen-Charpentier and Kojouharov, 2013). In this work, the author tests UPFD scheme under some other different regimes of advection, diffusion and reaction. The author considers the case when the coefficient of advection, diffusion and reaction are all equal to one and also cases under which advection or diffusion or reaction is more important. Some errors such as L1 error, dispersion, dissipation errors and relative errors are tabulated. Moreover, the author compares some spectral properties of the method under different regimes. The author obtains the variation of the following quantities with respect to the phase angle: modulus of exact amplification factor, modulus of amplification factor of the scheme and relative phase error.
Design/methodology/approach
Difficulties can arise in stability analysis. It is important to have a full understanding of whether the conditions obtained for stability are sufficient, necessary or necessary and sufficient. The advection-diffusion-reaction is quite similar to the advection-diffusion equation, it has an extra reaction term and therefore obtaining stability of numerical methods discretizing advection-diffusion-reaction equation is not easy as is the case with numerical methods discretizing advection-diffusion equations. To avoid difficulty involved with obtaining region of stability, the author shall consider unconditionally stable finite difference schemes discretizing advection-diffusion-reaction equations.
Findings
The UPFD scheme is unconditionally stable but not unconditionally consistent. The scheme was tested on an advection-diffusion-reaction equation which models exponential travelling waves, and the author computed various errors such as L1 error, dispersion and dissipation errors, relative errors under some different regimes of advection, diffusion and reaction. The scheme works best for very small values of k as k → 0 (for instance, k = 0.00025, 0.0005) and performs satisfactorily at other values of k such as 0.001 for two regimes; a = 1, D = 1, κ = 1 and a = 1, D = 1, κ = 5. When a = 5, D = 1, κ = 1, the scheme performs quite well at k = 0.00025 and satisfactorily at k = 0.0005 but is not efficient at larger values of k. For the diffusive case (a = 1, D = 5, κ = 1), the scheme does not perform well. In general, the author can conclude that the choice of k is very important, as it affects to a great extent the performance of the method.
Originality/value
The UPFD scheme is effective to solve advection-diffusion-reaction problems when advection or reactive regime is dominant and for the case, a = 1, D = 1, κ = 1, especially at low values of k. Moreover, the magnitude of the dispersion and dissipation errors using UPFD are of the same order for all the four regimes considered as seen from Tables 1 to 4. This indicates that if the author is to optimize the temporal step size at a given value of the spatial step size, the optimization function must consist of both the AFM and RPE. Some related work on optimization can be seen in Appadu (2013). Higher-order unconditionally stable schemes can be constructed for the regimes for which UPFD is not efficient enough for instance when advection and diffusion are dominant.
Details
Keywords
This paper aims to present an interpolating element-free Galerkin (IEFG) method for the numerical study of the time-fractional diffusion equation, and then discuss the stability…
Abstract
Purpose
This paper aims to present an interpolating element-free Galerkin (IEFG) method for the numerical study of the time-fractional diffusion equation, and then discuss the stability and convergence of the numerical solutions.
Design/methodology/approach
In the time-fractional diffusion equation, the time fractional derivatives are approximated by L1 method, and the shape functions are constructed by the interpolating moving least-squares (IMLS) method. The final system equations are obtained by using the Galerkin weak form. Because the shape functions have the interpolating property, the unknowns can be solved by the iterative method after imposing the essential boundary condition directly.
Findings
Both theoretical and numerical results show that the IEFG method for the time-fractional diffusion equation has high accuracy. The stability of the fully discrete scheme of the method on the time step is stable unconditionally with a high convergence rate.
Originality/value
This work will provide an interpolating meshless method to study the numerical solutions of the time-fractional diffusion equation using the IEFG method.
Details
Keywords
A.A. Aderogba, M. Chapwanya and J.K. Djoko
For a partial differential equation with a fourth-order derivative such as the Cahn-Hilliard equation, it is always a challenge to design numerical schemes that can handle the…
Abstract
Purpose
For a partial differential equation with a fourth-order derivative such as the Cahn-Hilliard equation, it is always a challenge to design numerical schemes that can handle the restrictive time step introduced by this higher order term. The purpose of this paper is to employ a fractional splitting method to isolate the convective, the nonlinear second-order and the fourth-order differential terms.
Design/methodology/approach
The full equation is then solved by consistent schemes for each differential term independently. In addition to validating the second-order accuracy, the authors will demonstrate the efficiency of the proposed method by validating the dissipation of the Ginzberg-Lindau energy and the coarsening properties of the solution.
Findings
The scheme is second-order accuracy, the authors will demonstrate the efficiency of the proposed method by validating the dissipation of the Ginzberg-Lindau energy and the coarsening properties of the solution.
Originality/value
The authors believe that this is the first time the equation is handled numerically using the fractional step method. Apart from the fact that the fractional step method substantially reduces computational time, it has the advantage of simplifying a complex process efficiently. This method permits the treatment of each segment of the original equation separately and piece them together, in a way that will be explained shortly, without destroying the properties of the equation.
Details
Keywords
C. Kassiotis, J.‐B. Colliat, A. Ibrahimbegovic and H.G. Matthies
The purpose of this paper is to study the partitioned solution procedure for thermomechanical coupling, where each sub‐problem is solved by a separate time integration scheme.
Abstract
Purpose
The purpose of this paper is to study the partitioned solution procedure for thermomechanical coupling, where each sub‐problem is solved by a separate time integration scheme.
Design/methodology/approach
In particular, the solution which guarantees that the coupling condition will preserve the stability of computations for the coupled problem is studied. The consideration is further generalized for the case where each sub‐problem will possess its particular time scale which requires different time step to be selected for each sub‐problem.
Findings
Several numerical simulations are presented to illustrate very satisfying performance of the proposed solution procedure and confirm the theoretical speed‐up of computations which follow from the adequate choice of the time step for each sub‐problem.
Originality/value
The paper confirms that one can make the most appropriate selection of the time step and carry out the separate computations for each sub‐problem, and then enforce the coupling which will preserve the stability of computations with such an operator split procedure.
Details
Keywords
Bartosz Górecki and Jacek Szumbarski
Both the importance of the natural convection in science and engineering and the shortage of publications in the field of numerical features of time-stepping schemes for the…
Abstract
Purpose
Both the importance of the natural convection in science and engineering and the shortage of publications in the field of numerical features of time-stepping schemes for the simulation of coupled heat and fluid flow problems motivate the present work. The paper aims to discuss these issues.
Design/methodology/approach
The paper presents the unconditionally stable time-stepping scheme for simulation of coupled problems of mass and heat transport. The paper is divided into two parts. The first part concerns the mathematical formulation of the scheme and discusses its implementation. The second part focuses on the numerical simulation and its results. A detailed investigation of the temporal order of the scheme with respect to the L2-norms of the errors of the pressure, velocity, temperature and divergence of velocity fields has also been given.
Findings
The work shows that it is possible to formulate a numerical scheme which is unconditionally stable with respect to the time step size. Moreover, application of the spectral element method for the spatial discretization results in a high order of approximation in space and very good overall accuracy. Furthermore, the investigation of the numerical features of the scheme showed that the formal temporal order of the scheme (formally second order) has been deferred very slightly and the order of 1.8-1.9 is achieved for all unknown fields.
Originality/value
The paper presents a new unconditionally stable scheme for simulation of unsteady flows with bidirectional coupling of heat transfer and the fluid flow. It also carefully investigates the numerical behaviour of the method.
Details
Keywords
Theodoros Zygiridis, Georgios Pyrialakos, Nikolaos Kantartzis and Theodoros Tsiboukis
The locally one-dimensional (LOD) finite-difference time-domain (FDTD) method features unconditional stability, yet its low accuracy in time can potentially become detrimental…
Abstract
Purpose
The locally one-dimensional (LOD) finite-difference time-domain (FDTD) method features unconditional stability, yet its low accuracy in time can potentially become detrimental. Regarding the improvement of the method’s reliability, existing solutions introduce high-order spatial operators, which nevertheless cannot deal with the augmented temporal errors. The purpose of the paper is to describe a systematic procedure that enables the efficient implementation of extended spatial stencils in the context of the LOD-FDTD scheme, capable of reducing the combined space-time flaws without additional computational cost.
Design/methodology/approach
To accomplish the goal, the authors introduce spatial derivative approximations in parametric form, and then construct error formulae from the update equations, once they are represented as a one-stage process. The unknown operators are determined with the aid of two error-minimization procedures, which equally suppress errors both in space and time. Furthermore, accelerated implementation of the scheme is accomplished via parallelization on a graphics-processing-unit (GPU), which greatly shortens the duration of implicit updates.
Findings
It is shown that the performance of the LOD-FDTD method can be improved significantly, if it is properly modified according to accuracy-preserving principles. In addition, the numerical results verify that a GPU implementation of the implicit solver can result in up to 100× acceleration. Overall, the formulation developed herein describes a fast, unconditionally stable technique that remains reliable, even at coarse temporal resolutions.
Originality/value
Dispersion-relation-preserving optimization is applied to an unconditionally stable FDTD technique. In addition, parallel cyclic reduction is adapted to hepta-diagonal systems, and it is proven that GPU parallelization can offer non-trivial benefits to implicit FDTD approaches as well.
Details
Keywords
Hooman Razmjoo and Masoud Movahhedi
In this paper, a modified meshless method, as one of the numerical techniques that has recently emerged in the area of computational electromagnetics, is extended to solving…
Abstract
Purpose
In this paper, a modified meshless method, as one of the numerical techniques that has recently emerged in the area of computational electromagnetics, is extended to solving time-domain wave equation. The paper aims to discuss these issues.
Design/methodology/approach
In space domain, the fields at the collocation points are expanded into a series of new Shepard's functions which have been suggested recently and are treated with a meshless method procedure. For time discretization of the second-order time-derivative, two finite-difference schemes, i.e. backward difference and Newmark-β techniques, are proposed.
Findings
Both schemes are implicit and always stable and have unconditional stability with different orders of accuracy and numerical dispersion. The unconditional stability of the proposed methods is analytically proven and numerically verified. Moreover, two numerical examples for electromagnetic field computation are also presented to investigate characteristics of the proposed methods.
Originality/value
The paper presents two unconditionally stable schemes for meshless methods in time-domain electromagnetic problems.
Details
Keywords
V.M Wheeler and K K Tamma
The purpose of this paper is to provide an overview and some recent advances in the models, analysis and simulation of thermal transport of phonons as related to the field of…
Abstract
Purpose
The purpose of this paper is to provide an overview and some recent advances in the models, analysis and simulation of thermal transport of phonons as related to the field of microscale/macroscale heat conduction in solids. The efforts focus upon a fairly comprehensive overview of the subject matter from a unified standpoint highlighting the various approximations inherent in the thermal models. Subsequently, the numerical formulations and illustrations using the current state-of-the-art are provided.
Design/methodology/approach
This paper is dedicated to the approximate solution to the relaxation time phonon Boltzmann equation (BE). While original contributions are pointed out and addressed appropriately, the efforts and contributions will be focussed on a relatively complete overview highlighting the field from one unified standpoint and clearly stating all assumptions that go into the approximations inherent to existing models. The contents will be divided as follows: In the first section the authors will give an overview of semi-classical phonon transport physics. Then the authors will discuss the equation of phonon radiative transport (EPRT) and its approximations—the ballistic-diffusive approximation (BDA) and the new heat equation (NHE). Next the authors derive and discuss the C-F model. A numerical discretization method valid for all models is then presented followed by results to numerical simulations and discussion.
Findings
From a unified treatment based on the introduction of an energy distribution function, the authors have derived the EPRT and its two well-known approximations: BDA and NHE. For completeness and to provide a vehicle for a general numerical discretization approach, the authors have also included analysis of the C-F model and the parabolic and hyperbolic descriptions of heat transfer along with it. The approximation of angular dependence of phonons in radiation-like descriptions of transport has been given special attention. The assumption of isotropy was found to be of paramount importance in the formulation of position space models for phononic thermal transport. For the thin film problem considered here, the NHE along with the proper boundary condition appears to be the best choice to approximate the phonon BE. Not only does it provide predictions that are in excellent agreement with EPRT, it does not require the discretization of phase space making it far more computationally efficient.
Originality/value
The authors hope this work will help dispel the idea that since Fourier’s law describes diffusion (under limiting assumptions) and it has shown to be ineffective in describing heat transfer for very thin films, that diffusion cannot describe heat transfer in thin films and one should look to a radiative description instead. If one considers diffusion in the sense of random motion, as invisaged by the original builders of the subject (Smoluchowski, Einstein, Ornstein et al.), instead of a temperature gradient, the idea that diffusion can govern thermal transport at this scale is not surprising. Indeed, the NHE is essentially a diffusion equation that describes the motion of particles up to the point of true randomness (isotropy) as well as thereafter.
Details
Keywords
Gerardo Abrugia and Umberto Amato
A mixed implicit semi Lagrangian finite difference‐finite volume method for numerical simulation of 2D air motion inside cylinders is derived and discussed. A conformal mapping…
Abstract
A mixed implicit semi Lagrangian finite difference‐finite volume method for numerical simulation of 2D air motion inside cylinders is derived and discussed. A conformal mapping from a physical (moving) domain to a computational (fixed) one is resorted in order to deal with a grid independent of time, making the numerical code very efficient. The numerical method is mass and energy conservative, unconditionally stable and at each timestep requires the solution of two well structured five‐band systems of linear equations. Its accuracy is first order in time and second one in space where the solution is smooth, while due to FCT space accuracy drops to the first order where the solution is steep. Stability of the method is proved both by a classical Von Newmann analysis and analysis of the matrices involved in the systems of linear equations. All these elements make the numerical method particularly fast. Numerical experiments are performed that show the influence of the maximum Courant number (with respect to the fluid speed) on the performance of the numerical method; moreover, comparison of simulations with a major existing code for engines is worked out.
Details
Keywords
M.J. Huntul and Mohammad Tamsir
The purpose of this paper is to reconstruct the potential numerically in the fourth-order Rayleigh–Love equation with boundary and nonclassical boundary conditions, from…
Abstract
Purpose
The purpose of this paper is to reconstruct the potential numerically in the fourth-order Rayleigh–Love equation with boundary and nonclassical boundary conditions, from additional measurement.
Design/methodology/approach
Although, the aforesaid inverse identification problem is ill-posed but has a unique solution. The authors use the cubic B-spline (CBS) collocation and Tikhonov regularization techniques to discretize the direct problem and to obtain stable as well as accurate solutions, respectively. The stability, for the discretized system of the direct problem, is also carried out by means of the von Neumann method.
Findings
The acquired results demonstrate that accurate as well as stable solutions for the a(t) are accessed for
Research limitations/implications
Since the noisy data are introduced, the investigation and analysis model real circumstances where the practical quantities are naturally infested with noise.
Practical implications
The acquired results demonstrate that accurate as well as stable solutions for the a(t) are accessed for
Originality/value
The potential term in the fourth-order Rayleigh–Love equation from additional measurement is reconstructed numerically, for the first time. The technique establishes that accurate, as well as stable solutions are obtained.
Details