Search

Highlights • A finite-buffer vacation queue with batch Markovian arrival process has been analyzed. • The server is subjected to serve under gated-limited service discipline. • Proposed analysis is based on the successive substitution and the supplementary variable method. • The results have been matched with the corresponding infinite-buffer queue. • Numerical results are presented for different service- and vacation-time distributions. Abstract This paper analyzes the finite-buffer single server queue with vacation(s). It is assumed that the arrivals follow a batch Markovian arrival process (BMAP) and the server serves customers according to a non-exhaustive type gated-limited service discipline. It has been also considered that the service and vacation distributions possess rational Laplace-Stieltjes transformation (LST) as these types of distributions may approximate many other distributions appeared in queueing literature. Among several batch acceptance/rejection strategies, the partial batch acceptance strategy is discussed in this paper. The service limit L (1 ≤ L ≤ N) is considered to be fixed, where N is the buffer-capacity excluding the one in service. It is assumed that in each busy period the server continues to serve until either L customers out of those that were waiting at the start of the busy period are served or the queue empties, whichever occurs first. The queue-length distribution at vacation termination/service completion epochs is determined by solving a set of linear simultaneous equations. The successive substitution method is used in the steady-state equations embedded at vacation termination/service completion epochs. The distribution of the queue-length at an arbitrary epoch has been obtained using the supplementary variable technique. The queue-length distributions at pre-arrival and post-departure epoch are also obtained. The results of the corresponding infinite-buffer queueing model have been analyzed briefly and matched with the previous model. Net profit function per unit of time is derived and an optimal service limit and buffer-capacity are obtained from a maximal expected profit. Some numerical results are presented in tabular and graphical forms. [ABSTRACT FROM AUTHOR]

In this paper, we significantly extend the applicability of state-of-the-art ELDP (equations for linearizing discrete product terms) method by providing a new linearization to handle more complicated non-linear terms involving both of discrete and bounded continuous variables. A general class of “representable programming problems” is formally proposed for a much wider range of engineering applications. Moreover, by exploiting the logarithmic feature embedded in the discrete structure, we present an enhanced linear reformulation model which requires half an order fewer equations than the original ELDP. Computational experiments on various engineering design problems support the superior computational efficiency of the proposed linearization reformulation in solving engineering optimization problems with discrete and bounded continuous variables. [ABSTRACT FROM AUTHOR]

In this paper, the dynamics of a viscoelastic plate resting on a viscoelastic Winkler foundation and traversed by a moving mass is studied. The Laplace transform is employed to derive the governing equation of the problem. Thereafter, an analytical-numerical method is proposed in order to determine the dynamic response of the plate. The method is based on transforming the governing partial differential equation into a new solvable system of linear ordinary differential equations. To that extent, the proposed solution proves to be applicable to plates made of any viscoelastic material and with various boundary conditions. Moreover, the moving mass may travel at any arbitrary trajectory with no restriction. The capability of the method is demonstrated through several illustrative examples, wherein a comprehensive parametric study is performed to investigate the effect of various parameters on the response of the system. A comparison between the results obtained by the complete moving mass approach with those coming from the approximate solutions is also carried out. The comparative study suggests that neglecting the inertia of the moving mass in the problems where a relatively heavy mass is traveling at high speeds on a plate with a low level of viscosity leads, in general, to an unsafe design. Additionally, ignoring the effect of convective accelerations in the moving mass formulation results in an overly uneconomical design in most cases. However, for large values of viscosity, the accuracy of the approximate methods further improves, and those methods can be used to reasonably estimate the response. [ABSTRACT FROM AUTHOR]

Abstract: Many physical phenomena are modeled by nonclassical hyperbolic boundary value problems with nonlocal boundary conditions. In this paper, the problem of solving the one-dimensional wave equation subject to given initial and non-local boundary conditions is considered. These non-local conditions arise mainly when the data on the boundary cannot be measured directly. Several finite difference methods with low order have been proposed in other papers for the numerical solution of this one dimensional non-classic boundary value problem. Here, we derive a new family of efficient three-level algorithms with higher order to solve the wave equation and also use a Simpson formula with higher order to approximate the integral conditions. Additionally, the fourth-order formula is also adapted to nonlinear equations, in particular to the well-known nonlinear Klein–Gordon equations which many physical problems can be modeled with. Numerical results are presented and are compared with some existing methods showing the efficiency of the new algorithms. [Copyright &y& Elsevier]

Abstract: In this paper, we employ the boundary-only meshfree method to find out numerical solution of the classical Boussinesq equation in one dimension. The proposed method in the current paper is a combination of boundary knot method and meshless analog equation method. The boundary knot technique is an integration free, boundary-only, meshless method which is used to avoid the known disadvantages of the method of fundamental solution. Also, we use the meshless analog equation method to replace the nonlinear governing equation with an equivalent nonhomogeneous linear equation. A predictor–corrector scheme is proposed to solve the resulted differential equation of the collocation. The numerical results and conclusions are obtained for both the ‘good’ and the ‘bad’ Boussinesq equations. [Copyright &y& Elsevier]

Abstract: In this paper, we consider a backward heat problem that appears in many applications. This problem is ill-posed. The solution of the problem as the solution exhibits unstable dependence on the given data functions. Using a new regularization method, we regularize the problem and get some new error estimates. Some numerical tests illustrate that the proposed method is feasible and effective. This work is a generalization of many recent papers, including the earlier paper [A new regularized method for two dimensional nonhomogeneous backward heat problem, Appl. Math. Comput. 215(3) (2009) 873–880] and some other authors such as Chu-Li Fu et al. , Campbell et al. . [Copyright &y& Elsevier]

Abstract: This paper presents a generic high dimensional model representation (HDMR) method for approximating the system response in terms of functions of lower dimensions. The proposed approach, which has been previously applied for problems dealing only with random variables, is extended in this paper for problems in which physical properties exhibit spatial random variation and may be modelled as random fields. The formulation of the extended HDMR is similar to the spectral stochastic finite element method in the sense that both of them utilize Karhunen–Loève expansion to represent the input, and lower-order expansion to represent the output. The method involves lower dimensional HDMR approximation of the system response, response surface generation of HDMR component functions, and Monte Carlo simulation. Each of the low order terms in HDMR is sub-dimensional, but they are not necessarily translating to low degree polynomials. It is an efficient formulation of the system response, if higher-order variable correlations are weak, allowing the physical model to be captured by the first few lower-order terms. Once the approximate form of the system response is defined, the failure probability can be obtained by statistical simulation. The proposed approach decouples the finite element computations and stochastic computations, and consecutively the finite element code can be treated as a black box, as in the case of a commercial software. Numerical examples are used to illustrate the features of the extended HDMR and to compare its performance with full scale simulation. [Copyright &y& Elsevier]

Abstract: In this paper, we proposed a modified model of Koiter’s type for nonlinearly elastic shell. The change of metric tensor and the change of curvature tensor play an important role in constructing the linearly and nonlinearly elastic shell model of Koiter’s type. The approximate expressions of them once were proposed by Ciarlet. In this paper, the exact full expressions of the change of metric tensor and the change of curvature tensor are provided by tensor analysis. The former coincides with Ciarlet’s expression. And the latter is more exact than Ciarlet’s expression. Thus the modified model is better than Ciarlet’s model. At the same time, a numerical experiment of special hemispherical shell is provided to validate the modified model of Koiter’s type. [Copyright &y& Elsevier]

The paper presents the approximation strategy for solutions and derivatives of solutions of boundary value problems on the example of two-dimensional solids. The effectiveness of the proposed strategy lies in the fact that it gives possibility to calculate solutions and their derivatives continuously at all points of the boundary and the domain, irrespective of the method used to solve the boundary problem and regardless of the type of the problem. The strategy has been developed in order to: (1) improve the accuracy of solutions (displacement) and their derivatives (strains, stresses) in the vicinity of the boundary, where they are affected by errors, (2) make the ability to directly obtain the derivatives of solutions (strains, stresses) on the boundary, (3) avoid computing strongly singular integrals present in the integral identity used for obtaining solutions or their derivatives (e.g. stresses in plasticity problems). The different variants of the proposed strategy have been developed and their accuracy has been verified considering examples with analytical solutions. [ABSTRACT FROM AUTHOR]

This paper presents a multi-domain Legendre–Gauss pseudospectral method for approximate solutions of the fractional Volterra’s model for population growth of a species in a closed system. Fractional Volterra’s model is a fractional Volterra integro-differential equation. The fractional derivative is considered in the Caputo sense. In this method, the fractional integro-differential equation is first replaced with a singular Volterra integro-differential equation (SVIDE). Then, by choosing a step-size, the replaced problem is converted into a sequence of SVIDEs in subintervals. The difficulty in SVIDEs due to singularity is overcome here by utilizing integration by parts. The obtained problems in subintervals are then step by step reduced to systems of algebraic equations using collocation. We give some numerical applications to show validity and high accuracy of the proposed technique. [ABSTRACT FROM AUTHOR]

Abstract: In this paper, a high-order compact stencil for solving the convection–diffusion equation in two dimensions is proposed. The convection and diffusion terms are both approximated by means of radial basis functions (RBFs) that are constructed over rectangular stencils. Salient features here are that (i) integration is employed to construct local RBF approximations; and (ii) through the constants of integration, values of the convection–diffusion equation at several selected nodes on the stencil are also enforced. Numerical results indicate that (i) the inclusion of the governing equation into the stencil leads to a significant improvement in accuracy; (ii) when the convection dominates, accurate solutions are obtained at a regime of the RBF width which makes the RBFs peaked; and (iii) high levels of accuracy are achieved using relatively coarse grids. [Copyright &y& Elsevier]

Abstract: In this paper, we provide a new modification of the variational iteration method (MVIM) for solving van der Pol equations. The modification couples the classical variational iteration method with He’s polynomials, where the He’s polynomials are applied to the approximate solution and the initial condition to eliminate secular terms. For the large ɛ, the numerical results demonstrate that the modification method get an accurate approximate period than the other presented methods. [Copyright &y& Elsevier]

Abstract: In this paper, we propose a new high accuracy numerical method of O(k 2 + k 2 h 2 + h 4) based on off-step discretization for the solution of 3-space dimensional non-linear wave equation of the form u tt = A(x,y,z,t)u xx + B(x,y,z,t)u yy + C(x,y,z,t)u zz + g(x,y,z,t,u,u x ,u y ,u z ,u t ), 0< x,y,z <1,t >0 subject to given appropriate initial and Dirichlet boundary conditions, where k >0 and h >0 are mesh sizes in time and space directions respectively. We use only seven evaluations of the function g as compared to nine evaluations of the same function discussed in [3,4]. We describe the derivation procedure in details of the algorithm. The proposed numerical algorithm is directly applicable to wave equation in polar coordinates and we do not require any fictitious points to discretize the differential equation. The proposed method when applied to a telegraphic equation is also shown to be unconditionally stable. Comparative numerical results are provided to justify the usefulness of the proposed method. [Copyright &y& Elsevier]

Abstract: In this paper, we consider the weak solutions of hyperbolic problems subject to inhomogeneous Dirichlet and Neumann boundary conditions. Using Fourier–Galerkin method, we obtain approximate solutions of the problems and test the obtained results on numerical examples by MAPLE®. [Copyright &y& Elsevier]

Abstract: A stabilized finite element method for the time-dependent Stokes equations based on Crank–Nicolson scheme is considered in this paper. The method combines the Crank–Nicolson scheme with a stabilized finite element method which uses the lowest equal-order element pair, i.e., the stabilized finite element method is applied for the spatial approximation and the time discretization is based on the Crank–Nicolson scheme. Moreover, we present optimal error estimates and prove that the scheme is unconditionally stable and convergent. Finally, numerical tests confirm the theoretical results of the presented method. [Copyright &y& Elsevier]

Abstract: In this paper, the solution of fuzzy linear system (FLS) is investigated based on a -level expansion. To this end, -level of FLS is solved for calculating the core of fuzzy solution and then its spreads are obtained by solving an optimization problem with a special objective function. Using our proposed method, if the computed solution satisfies the FLS, it can be as exact or approximated solution. Finally, the existence of solution of FLS is proved in details and some numerical examples are solved to illustrate the accuracy and capability of the method. [Copyright &y& Elsevier]

Abstract: In this paper, synchronization problems for a general complex networks are investigated by Takagi–Sugeno (T–S) fuzzy theory. A novel concept named linear approximation method is firstly proposed to solve the synchronization problems for T–S fuzzy complex networks. This novel method can linearize the system into some time-delay subsystems, which can effectively simplify the complicated system and then be easy to acquire the synchronization results. Since the expression based on linear matrix inequality (LMI) is used, the synchronization criteria can be easily checked in practice. Numerical simulation examples are provided to verify the effectiveness and the applicability of the proposed method. [Copyright &y& Elsevier]

Abstract: The present paper is devoted to solve the problem of identifying an unknown heat source depending simultaneously on both space and time variables. This problem is transformed into an optimization problem and the uniqueness of minimum element is proved rigorously. Then a variational formulation for solving the optimization problem is given. A conjugate gradient method and a finite difference method are used to solve the variational problem. Some numerical examples are also provided to show the efficiency of the proposed method. [Copyright &y& Elsevier]

Abstract: This paper aims to present a general framework of the cubic parametric spline functions to develop a numerical method for obtaining an approximate solution for the time fractional Burgers’ equation. The truncation error of the method is theoretically analyzed. Using Von Neumann method, the proposed method is also shown to be conditionally stable. Two numerical examples are then included to illustrate the practical implementation of the proposed method. The obtained results reveal that the proposed technique is very effective, convenient and quite accurate to such considered problems. [Copyright &y& Elsevier]

Abstract: In this paper, a system of M (⩾2) singularly perturbed semilinear reaction–diffusion equations is considered. To obtain a high order approximation to the solution of this system, we propose a hybrid numerical method that employs a generalized Shishkin mesh with the Numerov discretization in the boundary layer regions and either a non-equidistant generalization of the Numerov discretization or classical central differences in the outer region. It is proved that the method is almost fourth order convergent in the maximum norm uniformly with respect to the perturbation parameter. Numerical experiments support the theoretical results and demonstrate the effectiveness of the method. [Copyright &y& Elsevier]

Abstract: In this paper, we study the initial-boundary value problem of the usual Rosenau-RLW equation by finite difference method. We design a conservative numerical scheme which preserves the original conservative properties for the equation. The scheme is three-level and linear-implicit. The unique solvability of numerical solutions has been shown. Priori estimate and second order convergence of the finite difference approximate solutions are discussed by discrete energy method. Numerical results demonstrate that the scheme is efficient and accurate. [Copyright &y& Elsevier]

Abstract: In this paper, we successfully use 1D model to approximate the 3D problems. Firstly the PDEs (Partial Differential Equations) forms of Koiter’s Model for 2D linear elastic shell is proposed on special curvature coordinate system, i.e., spherical-coordinate system. Then the ODEs (Ordinary Differential Equations) forms of Koiter’s Model for 2D linear elastic shell is proposed under the assumption that the shell is axis-symmetric. Finally, we do numerical experiments to verify validity and accuracy of 1D models. [Copyright &y& Elsevier]

Abstract: The Klein–Gordon–Schrödinger equations describe a classical model of the interaction between conservative complex neutron field and neutral meson Yukawa in quantum field theory. In this paper, we study the long-time behavior of solutions for the Klein–Gordon–Schrödinger equations. We propose the Chebyshev pseudospectral collocation method for the approximation in the spatial variable and the explicit Runge–Kutta method in time discretization. In comparison with the single domain, the domain decomposition methods have good spatial localization and generate a sparse space differentiation matrix with high accuracy. In this study, we choose an overlapping multidomain scheme. The obtained numerical results show the Pseudospectral multidomain method has excellent long-time numerical behavior and illustrate the effectiveness of the numerical scheme in controlling two particles. Some comparisons with single domain pseudospectral and finite difference methods will be also investigated to confirm the efficiency of the new procedure. [Copyright &y& Elsevier]

Abstract: This paper considers bi-criteria scheduling on a single parallel-batch machine to minimizing two regular scheduling criteria that are non-decreasing in the job completion times. For minimizing the total weighted completion time and the makespan simultaneously, we prove that the problem is NP-hard in the ordinary sense and possesses a fully polynomial-time approximation scheme. For minimizing the weighted total number of tardy jobs and the makespan simultaneously, we provide a pseudo-polynomial-time algorithm and a fully polynomial-time approximation scheme. Furthermore, we identify that all Pareto-optimal solutions for and can be found in pseudo-polynomial time, respectively. [Copyright &y& Elsevier]

Abstract: In this paper, we consider numerical approximation of generalised Burgers–Fisher equation using the pseudo-spectral method. For the time discretization we apply Crank–Nicolson /leapfrog scheme. The space discretization is based on Legendre Galerkin formulation while the Chebyshev–Gauss–Lobatto (CGL) nodes are used in practical computation, which is called “Chebyshev–Legendre” method. The stability and convergence are rigorously set up. Numerical experiments are presented to demonstrate the effectiveness of the method and to confirm the theoretical results. [Copyright &y& Elsevier]

Abstract: In this paper a numerical method for solving the nonlinear age-structured population models is presented which is based on Bernstein polynomials approximation. Operational matrices of integration, differentiation, dual and product are introduced and are utilized to reduce the age-structured population problem to the solution of algebraic equations. The method in general is easy to implement, and yields good results. Illustrative examples are included to demonstrate the validity and applicability of the new technique. [Copyright &y& Elsevier]

Abstract: In this paper, we study the eventual periodicity of the initial boundary value problem (IBVP) for Korteweg-de Vries equation posed on a bounded domain. We show that if the boundary forcing is periodic of period τ, then the solution u of the IBVP at each spatial point becomes eventually time-periodic of period τ. In order to exhibit eventual periodicity, we approximate the solution of the IBVP using the Adomian decomposition method. We compare our work with the approximate solution of IBVP obtained by the homotopy perturbation method and present numerical experiments using Mathematica. [Copyright &y& Elsevier]

Abstract: In this paper the combined integral method is applied to a simple one-dimensional ablation problem. One of the drawbacks of heat balance integral methods is how to choose the approximating function. It is common to use a polynomial form but even then it is not clear what the power of the highest order term should be. Previous studies have determined exponents either from exact solutions or from expansions valid over short time scales; neither approach is satisfactory nor very accurate for larger times. We combine the heat balance and refined integral methods to determine this exponent as part of the solution process, and conclude that it is in fact time-dependent in the ablation stage. From comparing the approximate solutions with numerical and exact analytical solutions whenever possible, we show that this new method greatly improves the accuracy on standard methods, without overcomplicating the method. [Copyright &y& Elsevier]

Abstract: As an alternative to classic boundary conditions, new mean boundary conditions for image restoration problem have been recently introduced by Shi and Chang [Y. Shi, Q. Chang, Acceleration methods for image restoration problem with different boundary conditions, Applied Numerical Mathematics 58 (2008) 602–614]. In this paper, we propose an efficient scheme for computing the Kronecker product approximations of blurring matrices with the new mean boundary conditions. Our scheme does not require the symmetry condition of point spread functions. Detailed experiments in image restoration are given to demonstrate the efficiency of our scheme. [Copyright &y& Elsevier]

Abstract: In this paper, we state and prove a new formula expressing explicitly the derivatives of shifted Chebyshev polynomials of any degree and for any fractional-order in terms of shifted Chebyshev polynomials themselves. We develop also a direct solution technique for solving the linear multi-order fractional differential equations (FDEs) with constant coefficients using a spectral tau method. The spatial approximation with its fractional-order derivatives (described in the Caputo sense) are based on shifted Chebyshev polynomials T L,n (x) with x ∈(0, L), L >0 and n is the polynomial degree. We presented a shifted Chebyshev collocation method with shifted Chebyshev–Gauss points used as collocation nodes for solving nonlinear multi-order fractional initial value problems. Several numerical examples are considered aiming to demonstrate the validity and applicability of the proposed techniques and to compare with the existing results. [Copyright &y& Elsevier]

Abstract: Recently, fuzzy linear regression is considered by Mosleh et al. . In this paper, a novel hybrid method based on fuzzy neural network for approximate fuzzy coefficients (parameters) of fuzzy polynomial regression models with fuzzy output and crisp inputs, is presented. Here a neural network is considered as a part of a large field called neural computing or soft computing. Moreover, in order to find the approximate parameters, a simple algorithm from the cost function of the fuzzy neural network is proposed. Finally, we illustrate our approach by some numerical examples. [Copyright &y& Elsevier]

Abstract: In this paper, a novel hybrid method based on fuzzy neural network for approximate solution of fuzzy linear systems of the form Ax = Bx + d, where A and B are two square matrices of fuzzy coefficients, x and d are two fuzzy number vectors, is presented. Here a neural network is considered as a part of a large field called neural computing or soft computing. Moreover, in order to find the approximate solution, a simple and fast algorithm from the cost function of the fuzzy neural network is proposed. Finally, we illustrate our approach by some numerical examples. [Copyright &y& Elsevier]

Abstract: This paper deals with the problem of estimating the parameters for fractional Ornstein–Uhlenbeck processes from discrete observations when the Hurst parameter H is known. Both the drift and the diffusion coefficient estimators of discrete form are obtained based on approximating integrals via Riemann sums with Hurst parameter H ∈(1/2,3/4). By adapting the stochastic integral representation to the fractional Brownian motion, these two estimators can be efficiently computed by the use of computer software. Numerical examples are presented to examine the performance of our method. An application to real data is also presented to show how to apply this method in practice. [Copyright &y& Elsevier]

Abstract: In this paper, quadratic B-spline functions are used for solution of 2-D elastic problems. Because B-spline functions are directly used as basis function, there is no need to use meshes and nodes in function approximation. In order to improve the computational efficiency, different scales are used for sub-domains of entire problem domain in function approximation. The modified variational form and Lagrange multipliers method are used for coupling of different scale in function approximation. Compared with meshless methods and other wavelet based methods, this multi-scale B-spline-based method is simple and easy to work with for numerical analysis. Furthermore, the computational efficiency of the multi-scale method is much higher than that of single scale approach. The numerical examples of 2-D elastic problems indicate that the present method is effective and stable for solving complicated problems. [Copyright &y& Elsevier]

Abstract: In this paper we propose a numerical method to solve the Cauchy problem based on the viscous shallow water equations in an horizontally moving domain. More precisely, we are interested in a flooding and drying model, used to modelize the overflow of a river or the intrusion of a tsunami on ground. We use a nonconservative form of the two-dimensional shallow water equations, in eight velocity formulation and we build a numerical approximation, based on the Arbitrary lagrangian eulerian formulation, in order to compute the solution in the moving domain. [Copyright &y& Elsevier]

Abstract: In this paper, we investigate a Cauchy problem associated with Helmholtz-type equation in an infinite “strip”. This problem is well known to be severely ill-posed. The optimal error bound for the problem with only nonhomogeneous Neumann data is deduced, which is independent of the selected regularization methods. A framework of a modified Tikhonov regularization in conjunction with the Morozov’s discrepancy principle is proposed, it may be useful to the other linear ill-posed problems and helpful for the other regularization methods. Some sharp error estimates between the exact solutions and their regularization approximation are given. Numerical tests are also provided to show that the modified Tikhonov method works well. [Copyright &y& Elsevier]

Abstract: This paper presents the implementation of element free Galerkin method for the stress analysis of structures having cracks at the interface of two dissimilar materials. The material discontinuity at the interface has been modeled using a jump function with a jump parameter that governs its strength. The jump function enriches the approximation by the addition of special shape function that contains discontinuities in the derivative. The trial and test functions of the weak form are constructed using moving least-square interpolants in each material domain. An intrinsic enrichment criterion with enriched basis has been used to model the crack tip stress fields. The mixed mode (complex) stress intensity factors for bi-material interface cracks are numerically evaluated using the modified domain form of interaction integral. The numerical results are obtained for edge and center cracks lying at the bi-material interface, and are found to be in good agreement with the reference solutions for the interfacial crack problems. [Copyright &y& Elsevier]

Abstract: In this paper, we propose a new multilevel univariate multiquadric (MQ) quasi-interpolation approach with higher approximation order compared with the existing MQ quasi-interpolations. The proposed approach includes two schemes, which are based on inverse multiquadric radial basis function (IMQ-RBF) interpolation, and Wu & Schaback’s MQ quasi-interpolation operator . Error analysis shows that our operators produce higher accuracy. Numerical examples demonstrate that the proposed MQ quasi-interpolation schemes are valid. [ABSTRACT FROM AUTHOR]

Abstract: Lotfi et al. [Solving a full fuzzy linear programming using lexicography method and fuzzy approximate solution, Appl. Math. Modell. 33 (2009) 3151–3156] pointed out that there is no method in literature for finding the fuzzy optimal solution of fully fuzzy linear programming (FFLP) problems and proposed a new method to find the fuzzy optimal solution of FFLP problems with equality constraints. In this paper, a new method is proposed to find the fuzzy optimal solution of same type of fuzzy linear programming problems. It is easy to apply the proposed method compare to the existing method for solving the FFLP problems with equality constraints occurring in real life situations. To illustrate the proposed method numerical examples are solved and the obtained results are discussed. [ABSTRACT FROM AUTHOR]

Abstract: This paper presents a meshless method, which replaces the inhomogeneous biharmonic equation by two Poisson equations in terms of an intermediate function. The solution of the Poisson equation with the intermediate function as the right-hand term may be written as a sum of a particular solution and a homogeneous solution of a Laplace equation. The intermediate function is approximated by a series of radial basis functions. Then the particular solution is obtained via employing Kansa’s method, while the homogeneous solution is approximated by using the boundary radial point interpolation method by means of boundary integral equations. Besides, the proposed meshless method, in conjunction with the analog equation method, is further developed for solving generalized biharmonic-type problems. Some numerical tests illustrate the efficiency of the method proposed. [ABSTRACT FROM AUTHOR]

Abstract: In this paper, we discuss an -Galerkin mixed finite element method for a class of heat transport equations in the one-dimensional case. Optimal error estimates in -norm and -norm are obtained for both semi-discrete and fully discrete approximation scheme. Moreover, a numerical example is presented to illustrate the effectiveness of the method. [Copyright &y& Elsevier]

Abstract: This paper presents a simple model for the study of the dissolution of calcite and precipitation of gypsum in acidic solutions. The numerical predictions of calcium and hydroxyl ion concentration have been used to develop approximate analytical expressions for the growth rate of calcium in solution, and the hydroxyl ions concentration. These expressions have been used in turn to determine the evolution of the kinetic rate constants, and to make predictions that show good agreement with experimental data. [Copyright &y& Elsevier]

Abstract: This paper presents an efficient numerical algorithm for approximate solutions of a fractional population growth model in a closed system. The time-fractional derivative is considered in the Caputo sense. The algorithm is based on Adomian’s decomposition approach and the solutions are calculated in the form of a convergent series with easily computable components. Then the Padé approximants are effectively used in the analysis to capture the essential behavior of the population u(t) of identical individuals. [Copyright &y& Elsevier]