Search results

This paper aims to develop a new high accuracy numerical method based on off-step non-polynomial spline in tension approximations for the solution of Burgers-Fisher and coupled nonlinear Burgers’ equations on a graded mesh. The spline method reported here is third order accurate in space and second order accurate in time. The proposed spline method involves only two off-step points and a central point on a graded mesh. The method is two-level implicit in nature and directly derived from the continuity condition of the first order space derivative of the non-polynomial tension spline function. The linear stability analysis of the proposed method has been examined and it is shown that the proposed two-level method is unconditionally stable for a linear model problem. The method is directly applicable to problems in polar systems. To demonstrate the strength and utility of the proposed method, the authors have solved the generalized Burgers-Huxley equation, generalized Burgers-Fisher equation, coupled Burgers-equations and parabolic equation in polar coordinates. The authors show that the proposed method enables us to obtain the high accurate solution for high Reynolds number.

In this method, the authors use only two-level in time-direction, and at each time-level, the authors use three grid points for the unknown function u(x,t) and two off-step points for the known variable x in spatial direction. The methodology followed in this paper is the construction of a non-polynomial spline function and using its continuity properties to obtain consistency condition, which is third order accurate on a graded mesh and fourth order accurate on a uniform mesh. From this consistency condition, the authors derive the proposed numerical method. The proposed method, when applied to a linear equation is shown to be unconditionally stable. To assess the validity and accuracy, the method is applied to solve several benchmark problems, and numerical results are provided to demonstrate the usefulness of the proposed method.

The paper provides a third order numerical scheme on a graded mesh and fourth order spline method on a uniform mesh obtained directly from the consistency condition. In earlier methods, consistency conditions were only second order accurate. This brings an edge over other past methods. Also, the method is directly applicable to physical problems involving singular coefficients. So no modification in the method is required at singular points. This saves CPU time and computational costs.

There are no limitations. Obtaining a high accuracy spline method directly from the consistency condition is a new work. Also being an implicit method, this method is unconditionally stable.

Physical problems with singular and non-singular coefficients are directly solved by this method.

The paper develops a new method based on non-polynomial spline approximations of order two in time and three (four) in space, which is original and has lot of value because many benchmark problems of physical significance are solved in this method.

Mathematical and numerical models are developed to study the melting of a Phase Change Material (PCM) inside a 2D cavity. The bottom of the cell is heated at constant and uniform temperature or heat flux, assuming that the rest of the cavity is completely adiabatic. The paper used suitable numerical methods to follow the interface temporal evolution with a good accuracy. The purpose of this paper is to show how the evolution of the latent energy absorbed to melt the PCM depends on the temperature imposed on the lower wall of the cavity.

The problem is written with non-homogeneous boundary conditions. Momentum and energy equations are numerically solved in space by a spectral collocation method especially oriented to this situation. A Crank-Nicolson scheme permits the resolution in time.

The results clearly show the evolution of multicellular regime during the process of fusion and the kinetics of phase change depends on the boundary condition imposed on the bottom cell wall. Thus the charge and discharge processes in energy storage cells can be controlled by varying the temperature in the cell PCM. Substantial modifications of the thermal convective heat and mass transfer are highlighted during the transient regime. This model is particularly suitable to follow with a good accuracy the evolution of the solid/liquid interface in the process of storage/release energy.

The time-dependent physical properties that induce non-linear coupled unsteady terms in Navier-Stokes and energy equations are not taken into account in the present model. The present model is actually extended to these coupled situations. This problem requires smoother geometries. One can try to palliate this disadvantage by constructing smoother approximations of non-smooth geometries. The augmentation of polynomials developments orders increases strongly the computing time. When the external heat flux or temperature imposed at the PCM is much greater than the temperature of the PCM fusion, one must choose carefully some data to assume the algorithms convergence.

Among the areas where this work can be used, are: buildings where the PCM are used in insulation and passive cooling; thermal energy storage, the PCM stores energy by changing phase, solid to liquid (fusion); cooling and transport of foodstuffs or pharmaceutical or medical sensitive products, the PCM is used in the food industry, pharmaceutical and medical, to minimize temperature variations of food, drug or sensitive materials; and the textile industry, PCM materials in the textile industry are used in microcapsules placed inside textile fibres. The PCM intervene to regulate heat transfer between the body and the outside.

The paper's originality is reflected in the precision of its results, due to the use of a high-accuracy numerical approximation based on collocation spectral methods, and the choice of Chebyshev polynomials basis in both axial and radial directions.

The purpose of this paper is to develop a combined method for three-dimensional incompressible flow and heat transfer by the spectral collocation method (SCM) and the artificial compressibility method (ACM), and further to study the performance of the combined method SCM-ACM for three-dimensional incompressible flow and heat transfer.

The partial differentials in space are discretized by the SCM with Chebyshev polynomial and Chebyshev–Gauss–Lobbatto collocation points. The unsteady artificial compressibility equations are solved to obtain the steady results by the ACM. Three-dimensional exact solutions with trigonometric function form and exponential function form are constructed to test the accuracy of the combined method.

The SCM-ACM is developed successfully for three-dimensional incompressible flow and heat transfer with high accuracy that the minimum value of variance can reach. The accuracy increases exponentially along with time marching steps. The accuracy is also improved exponentially with the increasing of nodes before stable accuracy is achieved, while it keeps stably with the increasing of the time step. The central processing unit time increases exponentially with the increasing of nodes and decreasing of the time step.

It is difficult for the implementation of the implicit scheme by the developed SCM-ACM. The SCM-ACM can be used for solving unsteady impressible fluid flow and heat transfer.

The SCM-ACM is applied for two classic cases of lid-driven cavity flow and natural convection in cubic cavities. The present results show good agreement with the published results with much fewer nodes.

The combined method SCM-ACM is developed, firstly, for solving three-dimensional incompressible fluid flow and heat transfer by the SCM and ACM. The performance of SCM-ACM is investigated. This combined method provides a new choice for solving three-dimensional fluid flow and heat transfer with high accuracy.

The main purpose of this paper is to implement the piecewise spectral-variational iteration method (PSVIM) to solve the nonlinear Lane-Emden equations arising in mathematical physics and astrophysics.

This method is based on a combination of Chebyshev interpolation and standard variational iteration method (VIM) and matching it to a sequence of subintervals. Unlike the spectral method and the VIM, the proposed PSVIM does not require the solution of any linear or nonlinear system of equations and analytical integration.

Some well-known classes of Lane-Emden type equations are solved as examples to demonstrate the accuracy and easy implementation of this technique.

In this paper, a new and efficient technique is proposed to solve the nonlinear Lane-Emden equations. The proposed method overcomes the difficulties arising in calculating complicated and time-consuming integrals and terms that are not needed in the standard VIM.

This paper aims to describe a numerical procedure for approximating the potential distribution for a harmonic current point source, which is either buried in horizontally stratified multilayer earth, or positioned in the air. The procedure is very efficient and general. The total number of layers and the source position in relation to the medium model layers are completely arbitrary.

The efficiency of the computation procedure is based on the successful application of the numerical approximation of two kernel functions of the integral expression for the potential distribution within an arbitrarily chosen layer of the medium model. Each kernel function of the observed layer is approximated using a linear combination of 15 real exponential functions. Using these approximations and the analytical integration based on the Weber integral, a simple expression for numerical approximation of potential distribution within boundaries of the observed medium layer is given. Potential retardation is taken into account approximately.

The numerical procedure developed for the approximation of potential distribution for a harmonic current point source, which is positioned arbitrarily in air or in horizontally stratified multilayer earth, is efficient, numerically stable and generally applicable.

Numerical model developed for the harmonic current point source is the basis of a wider numerical models for computation of the harmonic and transient fields of earthing system, which consists of earthing grids buried in horizontally stratified multilayer earth and metallic structures in the air.

This is efficient and numerically stable frequency dependent harmonic current point source model. Potential retardation, which has been neglected at the first step of the approximation, is subsequently added to the potential expression in such a way that the Helmholtz differential equation has been approximately solved without introducing the Sommerfeld integrals.

In acceptance sampling, the hypergeometric operating characteristic (OC) function (so called type-A OC) is used to be approximated by the binomial or Poisson OC function, which actually reduce computational effort, but do not provide suffcient approximation results. The purpose of this paper is to examine binomial- and Poisson-type approximations to the hypergeometric distribution, in order to find a simple but accurate approximation that can be successfully applied in acceptance sampling.

The authors present a new binomial-type approximation for the type-A OC function, and derive its properties. Further, the authors compare this approximation via an extensive numerical study with other common approximations in terms of variation distance and relative efficiency under various conditions on the parameters including limiting cases.

The introduced approximation generates best numerical results over a wide range of parameter values, and ensures arithmetic simplicity of the binomial distribution and high accuracy to meet requirements regarding acceptance sampling problems. Additionally, it can considerably reduce the computational effort in relation to the type-A OC function and therefore is strongly recommended for calculating sampling plans.

The newly presented approximation provides a remarkably close fit to the type-A OC function, is discrete and needs no correction for continuity, and is skewed in the same direction by roughly the same amount as the exact OC. Due to less factorials, this OC in general involves lower powers than the type-A OC function. Moreover, the binomial-type approximation is easy to fit to the conventional statistical computing packages.

The purpose of this paper is to extend the generalized finite‐difference calculus of flexible local approximation methods (FLAME) to problems where local analytical solutions are unavailable.

FLAME uses accurate local approximations of the solution to generate difference schemes with small consistency errors. When local analytical approximations are too complicated, semi‐analytical or numerical ones can be used instead. In the paper, this strategy is applied to electrostatic multi‐particle simulations and to electromagnetic wave propagation and scattering. The FLAME basis is constructed by solving small local finite‐element problems or, alternatively, by a local multipole‐multicenter expansion. As yet another alternative, adaptive FLAME is applied to problems of wave propagation in electromagnetic (photonic) crystals.

Numerical examples demonstrate the high rate of convergence of new five‐ and nine‐point schemes in 2D and seven‐ and 19‐point schemes in 3D. The accuracy of FLAME is much higher than that of the standard FD scheme. This paves the way for solving problems with a large number of particles on relatively coarse grids. FLAME with numerical bases has particular advantages for the multi‐particle model of a random or quasi‐random medium.

Irregular stencils produced by local refinement may adversely affect the accuracy. This drawback could be rectified by least squares FLAME, where the number of stencil nodes can be much greater than the number of basis functions, making the method more robust and less sensitive to the irregularities of the stencils.

Previous applications of FLAME were limited to purely analytical basis functions. The present paper shows that numerical bases can be successfully used in FLAME when analytical ones are not available.

The purpose of this paper is to implement variational iteration method (VIM) and homotopy perturbation method (HPM) to solve modified Camassa‐Holm (mCH) and modified Degasperis‐Procesi (mDP) equations.

Perturbation method is a traditional method depending on a small parameter which is difficult to be found for real‐life nonlinear problems. To overcome the difficulties and limitations of the above method, two new ones have recently been introduced by He, i.e. VIM and HPM. In this paper, mCH and mDP equations are solved through these methods.

To assess the accuracy of the solutions, the comparison of the obtained results with the exact solutions reveals that both methods are tremendously effective.

The paper shows that VIM and HPM can be implemented to solve mCH and mDP equations.

A simple and efficient re‐analysis procedure is presented for large‐scale structural systems. The procedure is based on using a mixed formulation with the fundamental unknowns consisting of both stress and displacement parameters. The other key elements of the procedure are: (a) lumping of the large number of design variables into a single tracing parameter; (b) operator splitting or additive decomposition of the different arrays in the finite element equations of the modified structure into the corresponding arrays of the original structure plus correction terms; and (c) application of a reduction method through the successive use of the finite element method and the classical Bubnov‐Galerkin technique. The finite element method is first used to generate a few approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Bubnov—Galerkin technique. The re‐analysis procedure is applied to the linear static and free vibration problems of plate and shell structures. Changes in both the sizing and shape (configuration) design variables are considered. The high accuracy of the proposed technique, for sizable changes in the design variables, is demonstrated by means of numerical examples of composite plates and shells.

The purpose of this paper is to propose a new three-node triangular element in the framework of the numerical manifold method (NMM), which is designated by Trig3-MLScns.

The formulation uses the improved parametric shape functions of classical triangular elements (Trig3-0) to construct the partition of unity (PU) and the moving least square (MLS) interpolation method to construct the local approximation function.

Compared with the classical three-node element (Trig3-0), the Trig3-MLScns element has a higher order of approximations, much better accuracy and continuous nodal stress. Moreover, the linear dependence problem associated with many PU-based methods with high-order approximations is eliminated in the present element. A number of numerical examples indicate the high accuracy and robustness of the Trig3-MLScns element.

The proposed element inherits the individual merits of the NMM and the MLS.