A new domain decomposition algorithm for generalized Burger’s–Huxley equation based on Chebyshev polynomials and preconditioning

In this study, we use the spectral collocation method using Chebyshev polynomials for spatial derivatives and fourth order Runge–Kutta method for time integration to solve the generalized Burger’s–Huxley equation (GBHE). To reduce round-off error in spectral collocation (pseudospectral) method we use preconditioning. Firstly, theory of application of Chebyshev spectral collocation method with preconditioning (CSCMP) and domain decomposition on the generalized Burger’s–Huxley equation presented. This method yields a system of ordinary differential algebric equations (DAEs). Secondly, we use fourth order Runge–Kutta formula for the numerical integration of the system of DAEs. The numerical results obtained by this way have been compared with the exact solution to show the efficiency of the method.     

An adaptive domain decomposition method for the Hamilton–Jacobi–Bellman equation

In this paper, we propose an efficient algorithm for a Hamilton–Jacobi–Bellman equation governing a class of optimal feedback control and stochastic control problems. This algorithm is based on a non-overlapping domain decomposition method and an adaptive least-squares collocation radial basis function discretization with a novel matrix inversion technique. To demonstrate the efficiency of this method, numerical experiments on test problems with up to three states and two control variables have been performed. The numerical results show that the proposed algorithm is highly parallelizable and its computational cost decreases exponentially as the number of sub-domains increases.     

A computational meshless method for the generalized Burger’s–Huxley equation

A numerical solution of the generalized Burger’s–Huxley equation, based on collocation method using Radial basis functions (RBFs), called Kansa’s approach is presented. The numerical results are compared with the exact solution, Adomian decomposition method (ADM) and Variational iteration method (VIM). Highly accurate and efficient results are obtained by RBFs method. Excellent agreement with the exact solution is observed while better (or same) accuracy is obtained than other numerical schemes cited in this work.     

The application of the modified variational iteration method on the generalized Fisher’s equation

In a recent paper, Abassy et al. (J. Comput. Appl. Math. 207:137–147, 2007) proposed a modified variational iteration method (MVIM) for a special kind of nonlinear differential equations. In this paper, we consider variational iteration method (VIM) and MVIM (proposed in Abassy et al., J. Comput. Appl. Math. 207:137–147, 2007) to obtain an approximate series solution to the generalized Fisher’s equation which converges to the exact solution in the region of convergence. It is also shown that the application of VIM to the generalized Fisher’s equation leads to calculation of unneeded terms for series solution. Therefore, we use MVIM to overcome this disadvantage. Comparison of error between VIM and MVIM is made. The results show that the MVIM is more effective than the VIM.     

A continuation–domain decomposition algorithm for bifurcation problems

We study numerical solution branches of certain parameter-dependent problems defined on compact domains with various boundary conditions. The finite differences combined with the domain decomposition method are exploited to discretize the partial differential equations. We propose efficient numerical algorithms for solving the associated linear systems and for the detection of bifurcation points. Sample numerical results are reported. This revised version was published online in August 2006 with corrections to the Cover Date.     

On generalized Rubel’s equation

We solve generalized the generalized Rubel equation on the space of analytic functions in domains.     

A non-overlapping domain decomposition dual reciprocity method for solving the forward–backward heat equation in two-dimension

We apply a boundary element dual reciprocity method (DRBEM) to the numerical solution of the forward–backward heat equation in a two-dimensional case. The method is employed for the spatial variable via the fundamental solution of the Laplace equation and the Crank–Nicolson finite difference scheme is utilized to treat the time variable. The physical domain is divided into two non-overlapping subdomains resulting in two standard forward and backward parabolic equations. The subproblems are then treated by the underlying method assuming a virtual boundary in the interface and starting with an initial approximate solution on this boundary followed by updating the solution by an iterative procedure. In addition, we show that the time discrete scheme is unconditionally stable and convergent using the energy method. Furthermore, some computational aspects will be suggested to efficiently deal with the formulation of the proposed method. Finally, two forward–backward problems, for which the exact solution is available, will be numerically solved for two different domains to demonstrate the efficiency of the proposed approach.     

A criterion for the unique solvability of the dirichlet spectral problem in a cylindrical domain for a multidimensional Lavrent’ev-Bitsadze equation

We obtain a criterion for the unique solvability of the spectral problem in a cylindrical domain for a multidimensional Lavrent’ev-Bitsadze equation.     

Exact and numerical solutions for non-linear Burger’s equation by VIM

In this article variational iteration method (VIM), established by He in (1999), is considered to solve nonlinear Bergur’s equation. This method is a powerful tool for solving a large number of problems. Using variational iteration method, it is possible to find the exact solution or a closed approximate solution of a problem. Comparing the results with those of Adomian’s decomposition and finite difference methods reveals significant points. To illustrate the ability and reliability of the method, some examples are provided.     

Global solutions for the generalized Camassa–Holm equation

In this paper, we study the Cauchy problem of the generalized Camassa–Holm equation. Firstly, we prove the existence of the global strong solutions provide the initial data satisfying a certain sign condition. Then, we obtain the existence and the uniqueness of the global weak solutions under the same sign condition of the initial data.     

The Stefan problem for the Fisher–KPP equation

We study the Fisher–KPP equation with a free boundary governed by a one-phase Stefan condition. Such a problem arises in the modeling of the propagation of a new or invasive species, with the free boundary representing the propagation front. In one space dimension this problem was investigated in Du and Lin (2010) [11], and the radially symmetric case in higher space dimensions was studied in Du and Guo (2011) [10]. In both cases a spreading-vanishing dichotomy was established, namely the species either successfully spreads to all the new environment and stabilizes at a positive equilibrium state, or fails to establish and dies out in the long run; moreover, in the case of spreading, the asymptotic spreading speed was determined. In this paper, we consider the non-radially symmetric case. In such a situation, similar to the classical Stefan problem, smooth solutions need not exist even if the initial data are smooth. We thus introduce and study the “weak solution” for a class of free boundary problems that include the Fisher–KPP as a special case. We establish the existence and uniqueness of the weak solution, and through suitable comparison arguments, we extend some of the results obtained earlier in Du and Lin (2010) [11] and Du and Guo (2011) [10] to this general case. We also show that the classical Aronson–Weinberger result on the spreading speed obtained through the traveling wave solution approach is a limiting case of our free boundary problem here.     

New patterns of travelling waves in the generalized Fisher–Kolmogorov equation

We prove the existence and uniqueness of a family of travelling waves in a degenerate (or singular) quasilinear parabolic problem that may be regarded as a generalization of the semilinear Fisher–Kolmogorov–Petrovski–Piscounov equation for the advance of advantageous genes in biology. Depending on the relation between the nonlinear diffusion and the nonsmooth reaction function, which we quantify precisely, we investigate the shape and asymptotic properties of travelling waves. Our method is based on comparison results for semilinear ODEs.     

A Fourier spectral method for fractional-in-space Cahn–Hilliard equation

In this paper, a fractional extension of the Cahn–Hilliard (CH) phase field model is proposed, i.e. the fractional-in-space CH equation. The fractional order controls the thickness and the lifetime of the interface, which is typically diffusive in integer order case. An unconditionally energy stable Fourier spectral scheme is developed to solve the fractional equation with periodic or Neumann boundary conditions. This method is of spectral accuracy in space and of second-order accuracy in time. The main advantages of this method are that it yields high precision and high efficiency. Moreover, an extra stabilizing term is added to obey the energy decay property while maintaining accuracy and simplicity. Numerical experiments are presented to confirm the accuracy and effectiveness of the proposed method.     

A generalized Duffing equation with the Coulomb’s friction law and Signorini–type contact conditions

This work provides mathematical and numerical analyses for a spring–mass system, in which Signorini–type contact conditions and Coulomb’s friction law with thermal effects are taken into consideration. The motion of a mass attached to a viscoelastic (Kelvin–Voigt type) nonlinear spring is described by a generalized Duffing equation. Signorini contact conditions are understood as extended complementarity conditions (CCs), where convolution is incorporated, allowing to consider thermal aspects of an obstacle. We prove the existence of global weak solutions for the highly nonlinear differential equation system with all the conditions, based on the regularized differential equation and the normal compliance condition with the standard mollifier. In addition, we investigate what side effects produce higher singularities of contact forces in dynamic contact problems, which is also supported by numerical evidences. Numerical schemes are proposed and then several groups of data are selected for the display of our numerical simulations.