Using operational matrix for solving nonlinear class of mixed Volterra–Fredholm integral equations

This paper presents a computational technique for the solution of the nonlinear mixed Volterra–Fredholm integral equations of the second kind. Using the properties of three‐dimensional modification of hat functions, these are types of equations to a nonlinear system of algebraic equations. Also, convergence results and error analysis are discussed. The efficiency and accuracy of the proposed method is illustrated by numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.     

Mixed two‐grid finite difference methods for solving one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations

The aim of this paper is to propose mixed two‐grid finite difference methods to obtain the numerical solution of the one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations. The finite difference equations at all interior grid points form a large‐sparse linear system, which needs to be solved efficiently. The solution cost of this sparse linear system usually dominates the total cost of solving the discretized partial differential equation. The proposed method is based on applying a family of finite difference methods for discretizing the spatial and time derivatives. The obtained system has been solved by two‐grid method, where the two‐grid method is used for solving the large‐sparse linear systems. Also, in the proposed method, the spectral radius with local Fourier analysis is calculated for different values of h and Δt. The numerical examples show the efficiency of this algorithm for solving the one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical solution of nonlinear stochastic Itô‐Volterra integral equations driven by fractional Brownian motion

In this paper, an efficient numerical technique is applied to provide the approximate solution of nonlinear stochastic Itô‐Volterra integral equations driven by fractional Brownian motion with Hurst parameter . The proposed method is based on the operational matrices of modification of hat functions (MHFs) and the collocation method. In this approach, by approximating functions that appear in the integral equation by MHFs and using Newton's‐Cotes points, nonlinear integral equation is transformed to nonlinear system of algebraic equations. This nonlinear system is solved by using Newton's numerical method, and the approximate solution of integral equation is achieved. Some theorems related to error estimate and convergence analysis of the suggested scheme are also established. Finally, 2 illustrative examples are included to confirm applicability, efficiency, and accuracy of the proposed method. It should be noted that this scheme can be used to solve other appropriate problems, but some modifications are required.     

Chebyshev‐Legendre spectral method for solving the two‐dimensional vorticity equations with homogeneous Dirichlet conditions

The Chebyshev‐Legendre spectral method for the two‐dimensional vorticity equations is considered. The Legendre Galerkin Chebyshev collocation method is used with the Chebyshev‐Gauss collocation points. The numerical analysis results under the L²‐norm for the Chebyshev‐Legendre method of one‐dimensional case are generalized into that of the two‐dimensional case. The stability and optimal order convergence of the method are proved. Numerical results are given. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Numerical solution of nonlinear partial quadratic integro‐differential equations of fractional order via hybrid of block‐pulse and parabolic functions

In this paper, an effective numerical approach based on a new two‐dimensional hybrid of parabolic and block‐pulse functions (2D‐PBPFs) is presented for solving nonlinear partial quadratic integro‐differential equations of fractional order. Our approach is based on 2D‐PBPFs operational matrix method together with the fractional integral operator, described in the Riemann–Liouville sense. The main characteristic behind this approach is to reduce such problems to those of solving systems of algebraic equations, which greatly simplifies the problem. By using Newton's iterative method, this system is solved, and the solution of fractional nonlinear partial quadratic integro‐differential equations is achieved. Convergence analysis and an error estimate associated with the proposed method is obtained, and it is proved that the numerical convergence order of the suggested numerical method is O(h³) . The validity and applicability of the method are demonstrated by solving three numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the exact solutions much easier.     

Mathematical and numerical modellization of large‐scale oceanic waves

This paper is devoted to the theoretical and numerical study of a method which computes the variability of current and density in an oceanic domain. The equations are of Navier–Stokes type for the velocity and of transport‐diffusion type for the density. They are linearized around a given mean circulation and modified by physical assumptions including hydrostatic approximation. The existence and uniqueness of a solution are proved for two sets of equations: first the three‐dimensional problem and then the two‐dimensional cyclic problem derived by assuming a sinusoïdal x‐dependence for the perturbation of the mean flow. The latter corresponds to a modellization of tropical instability waves which are illustrated by the ‘El Nino’ phenomenon. These two problems differ from classical ones because of hydrostatic approximation, boundary conditions imposed by the oceanic domain and complex‐valued functions for the cyclic case. A numerical model is developed for the two‐dimensional cyclic equations. Time discretization is performed by the characteristics method; space discretization uses Q₁ finite elements. Numerical results are presented in a realistic case corresponding to the tropical Pacific Ocean. Copyright © 1999 John Wiley & Sons. Ltd.     

The use of cubic B‐spline scaling functions for solving the one‐dimensional hyperbolic equation with a nonlocal conservation condition

Problems for parabolic partial differential equations with nonlocal boundary conditions have been studied in many articles, but boundary value problems for hyperbolic partial differential equations have so far remained nearly uninvestigated. In this article a numerical technique is presented for the solution of a nonclassical problem for the one‐dimensional wave equation. This method uses the cubic B‐spline scaling functions. Some numerical results are reported to support our study. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Efficient numerical schemes for solving the self‐consistent field equations of flexible–semiflexible diblock copolymers

We present two efficient iterative schemes for solving the self‐consistent field equations of flexible–semiflexible diblock copolymers. One is a semi‐implicit scheme developed by employing asymptotic expansion, and the other is a hybrid scheme combining the robustness of the steepest descent method with the efficiency of the conjugate gradient method. In our position‐one‐dimensional and position‐two‐dimensional numerical experiments, we demonstrate that these schemes are much more efficient than the steepest descent method. Copyright © 2013 John Wiley & Sons, Ltd.     

On the stabilizing effect of convection in three‐dimensional incompressible flows

We investigate the stabilizing effect of convection in three‐dimensional incompressible Euler and Navier‐Stokes equations. The convection term is the main source of nonlinearity for these equations. It is often considered destabilizing although it conserves energy due to the incompressibility condition. In this paper, we show that the convection term together with the incompressibility condition actually has a surprising stabilizing effect. We demonstrate this by constructing a new three‐dimensional model that is derived for axisymmetric flows with swirl using a set of new variables. This model preserves almost all the properties of the full three‐dimensional Euler or Navier‐Stokes equations except for the convection term, which is neglected in our model. If we added the convection term back to our model, we would recover the full Navier‐Stokes equations. We will present numerical evidence that seems to support that the three‐dimensional model may develop a potential finite time singularity. We will also analyze the mechanism that leads to these singular events in the new three‐dimensional model and how the convection term in the full Euler and Navier‐Stokes equations destroys such a mechanism, thus preventing the singularity from forming in a finite time. © 2008 Wiley Periodicals, Inc.     

The Runge‐Kutta DG finite element method and the KFVS scheme for compressible flow simulations

This article presents a new type of second‐order scheme for solving the system of Euler equations, which combines the Runge‐Kutta discontinuous Galerkin (DG) finite element method and the kinetic flux vector splitting (KFVS) scheme. We first discretize the Euler equations in space with the DG method and then the resulting system from the method‐of‐lines approach will be discretized using a Runge‐Kutta method. Finally, a second‐order KFVS method is used to construct the numerical flux. The proposed scheme preserves the main advantages of the DG finite element method including its flexibility in handling irregular solution domains and in parallelization. The efficiency and effectiveness of the proposed method are illustrated by several numerical examples in one‐ and two‐dimensional spaces. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

The use of Chebyshev cardinal functions for solution of the second‐order one‐dimensional telegraph equation

A numerical technique is presented for the solution of the second order one‐dimensional linear hyperbolic equation. This method uses the Chebyshev cardinal functions. The method consists of expanding the required approximate solution as the elements of Chebyshev cardinal functions. Using the operational matrix of derivative, the problem is reduced to a set of algebraic equations. Some numerical examples are included to demonstrate the validity and applicability of the technique. The method is easy to implement and produces very accurate results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A meshless method based on the dual reciprocity method for one‐dimensional stochastic partial differential equations

This article describes a new meshless method based on the dual reciprocity method (DRM) for the numerical solution of one‐dimensional stochastic heat and advection–diffusion equations. First, the time derivative is approximated by the time–stepping method to transforming the original stochastic partial differential equations (SPDEs) into elliptic SPDEs. The resulting elliptic SPDEs have been approximated with the new method, which is a combination of radial basis functions (RBFs) method and the DRM method. We have used inverse multiquadrics (IMQ) and generalized IMQ (GIMQ) RBFs, to approximate functions in the presented method. The noise term has been approximated at the source points, at each time step. The developed formulation is verified in two test problems with investigating the convergence and accuracy of numerical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 292–306, 2016     

Numerical scheme for solving system of fractional partial differential equations with Volterra‐type integral term through two‐dimensional block‐pulse functions

In the current study, an approximate scheme is established for solving the fractional partial differential equations (FPDEs) with Volterra integral terms via two‐dimensional block‐pulse functions (2D‐BPFs). According to the definitions and properties of 2D‐BPFs, the original problem is transformed into a system of linear algebra equations. By dispersing the unknown variables for these algebraic equations, the numerical solutions can be obtained. Besides, the proof of the convergence of this system is given. Finally, several numerical experiments are presented to test the feasibility and effectiveness of the proposed method.     

Two‐dimensional shifted Legendre polynomial collocation method for electromagnetic waves in dielectric media via almost operational matrices

In this paper, a numerical solution of fractional partial differential equations (FPDEs) for electromagnetic waves in dielectric media will be discussed. For the solution of FPDEs, we developed a numerical collocation method using an algorithm based on two‐dimensional shifted Legendre polynomials approximation, which is proposed for electromagnetic waves in dielectric media. By implementing the partial Riemann–Liouville fractional derivative operators, two‐dimensional shifted Legendre polynomials approximation and its operational matrix along with collocation method are used to convert FPDEs first into weakly singular fractional partial integro‐differential equations and then converted weakly singular fractional partial integro‐differential equations into system of algebraic equation. Some results concerning the convergence analysis and error analysis are obtained. Illustrative examples are included to demonstrate the validity and applicability of the technique. Copyright © 2017 John Wiley & Sons, Ltd.     

Semi‐implicit schemes for transient Navier–Stokes equations and eddy viscosity models

This study presents two computational schemes for the numerical approximation of solutions to eddy viscosity models as well as transient Navier–Stokes equations. The eddy viscosity model is one example of a class of Large Eddy Simulation models, which are used to simulate turbulent flow. The first approximation scheme is a first order single step method that treats the nonlinear term using a semi‐implicit discretization. The second scheme employs a two step approach that applies a Crank–Nicolson method for the nonlinear term while also retaining the semi‐implicit treatment used in the first scheme. A finite element approximation is used in the spatial discretization of the partial differential equations. The convergence analysis for both schemes is discussed in detail, and numerical results are given for two test problems one of which is the two dimensional flow around a cylinder. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Composite spectral method for solution of the diffusion equation with specification of energy

Many physical subjects are modeled by nonclassical parabolic boundary value problems with nonlocal boundary conditions replacing the classic boundary conditions. In this article, we introduce a new numerical method for solving the one‐dimensional parabolic equation with nonlocal boundary conditions. The approximate proposed method is based upon the composite spectral functions. The properties of composite spectral functions consisting of terms of orthogonal functions are presented and are utilized to reduce the problem to some algebraic equations. The method is easy to implement and yields very accurate result. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Analytic solutions of a (2 + 1)‐dimensional variable‐coefficient Broer–Kaup system

Based on a Riccati equation and one of its new generalized solitary solutions constructed by the Exp‐function method, new analytic solutions with free parameters and arbitrary functions of a (2 + 1)‐dimensional variable‐coefficient Broer–Kaup system are obtained. These free parameters and arbitrary functions reveal that the (2 + 1)‐dimensional variable‐coefficient Broer–Kaup system has rich spatial structures. As an illustrative example, two new spatial structures are shown by setting the arbitrary functions as different Jacobi elliptic functions. Compared with tanh‐function method and its extensions, the method proposed in this paper is more powerful and it can be applied to other nonlinear evolution equations. Copyright © 2010 John Wiley & Sons, Ltd.     

Efficient numerical methods for pricing American options under stochastic volatility

Five numerical methods for pricing American put options under Heston's stochastic volatility model are described and compared. The option prices are obtained as the solution of a two‐dimensional parabolic partial differential inequality. A finite difference discretization on nonuniform grids leading to linear complementarity problems with M‐matrices is proposed. The projected SOR, a projected multigrid method, an operator splitting method, a penalty method, and a componentwise splitting method are considered. The last one is a direct method while all other methods are iterative. The resulting systems of linear equations in the operator splitting method and in the penalty method are solved using a multigrid method. The projected multigrid method and the componentwise splitting method lead to a sequence of linear complementarity problems with one‐dimensional differential operators that are solved using the Brennan and Schwartz algorithm. The numerical experiments compare the accuracy and speed of the considered methods. The accuracies of all methods appear to be similar. Thus, the additional approximations made in the operator splitting method, in the penalty method, and in the componentwise splitting method do not increase the error essentially. The componentwise splitting method is the fastest one. All multigrid‐based methods have similar rapid grid independent convergence rates. They are about two or three times slower that the componentwise splitting method. On the coarsest grid the speed of the projected SOR is comparable with the multigrid methods while on finer grids it is several times slower. ©John Wiley & Sons, Inc. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Convergence analysis of a vertex‐centered finite volume scheme for a copper heap leaching model

In this paper a two‐dimensional solute transport model is considered to simulate the leaching of copper ore tailing using sulfuric acid as the leaching agent. The mathematical model consists in a system of differential equations: two diffusion–convection‐reaction equations with Neumann boundary conditions, and one ordinary differential equation. The numerical scheme consists in a combination of finite volume and finite element methods. A Godunov scheme is used for the convection term and an P1‐FEM for the diffusion term. The convergence analysis is based on standard compactness results in L². Some numerical examples illustrate the effectiveness of the scheme. Copyright © 2009 John Wiley & Sons, Ltd.     

The finite volume method based on the Crouzeix–Raviart element for a fracture model

In this paper, we present and analyze a finite volume method based on the Crouzeix–Raviart element for the coupled fracture model, where the fluid flow is governed by Darcy's law in the one‐dimensional fracture and two‐dimensional surrounding matrix. In the numerical scheme, the pressure in the matrix and fracture is respectively approximated by the Crouzeix–Raviart elements and piecewise constant functions, and then the velocity is calculated by piecewise constant functions element by element. The existence and uniqueness of the numerical solution are discussed, and optimal order error estimates for both the pressure p and the velocity u are proved on general triangulations. We finally carry out numerical experiments, and results confirm our theoretical analysis.