A bivariate Chebyshev spectral collocation quasilinearization method for nonlinear evolution parabolic equations.
ABSTRACT: This paper presents a new method for solving higher order nonlinear evolution partial differential equations (NPDEs). The method combines quasilinearisation, the Chebyshev spectral collocation method, and bivariate Lagrange interpolation. In this paper, we use the method to solve several nonlinear evolution equations, such as the modified KdV-Burgers equation, highly nonlinear modified KdV equation, Fisher's equation, Burgers-Fisher equation, Burgers-Huxley equation, and the Fitzhugh-Nagumo equation. The results are compared with known exact analytical solutions from literature to confirm accuracy, convergence, and effectiveness of the method. There is congruence between the numerical results and the exact solutions to a high order of accuracy. Tables were generated to present the order of accuracy of the method; convergence graphs to verify convergence of the method and error graphs are presented to show the excellent agreement between the results from this study and the known results from literature.
Project description:A numerical solution of the modified Burgers' equation (MBE) is obtained by using quartic B-spline subdomain finite element method (SFEM) over which the nonlinear term is locally linearized and using quartic B-spline differential quadrature (QBDQM) method. The accuracy and efficiency of the methods are discussed by computing L 2 and L ∞ error norms. Comparisons are made with those of some earlier papers. The obtained numerical results show that the methods are effective numerical schemes to solve the MBE. A linear stability analysis, based on the von Neumann scheme, shows the SFEM is unconditionally stable. A rate of convergence analysis is also given for the DQM.
Project description:The rarefactive KdV solitary waves in a dusty plasma have been extensively studied analytically and found experimentally in the previous works. Though the envelope solitary wave described by a nonlinear Schrödinger equation (NLSE) has been proposed by using the reductive perturbation method, it is first verified by using the particle-in-cell (PIC) numerical method in this paper. Surprisingly, there is no phase shift after the head on collision between two envelope solitary waves, while it is sure that there are phase shifts of two colliding KdV solitary waves after head on collision.
Project description:A mechanism for dispersion to automatically arise from the dispersionless Whitham Modulation equations (WMEs) is presented, relying on the use of a moving frame. The speed of this is chosen to be one of the characteristics which emerge from the linearization of the Whitham system, and assuming these are real (and thus the WMEs are hyperbolic) morphs the WMEs into the Korteweg-de Vries (KdV) equation in the boosted coordinate. Strikingly, the coefficients of the KdV equation are universal, in the sense that they are determined by abstract properties of the original Lagrangian density. Two illustrative examples of the theory are given to illustrate how the KdV may be constructed in practice. The first being a revisitation of the derivation of the KdV equation from shallow water flows, to highlight how the theory of this paper fits into the existing literature. The second is a complex Klein-Gordon system, providing a case where the KdV equation may only arise with the use of a moving frame.
Project description:This work presents the novel Leal-polynomials (LP) for the approximation of nonlinear differential equations of different kind. The main characteristic of LPs is that they satisfy multiple expansion points and its derivatives as a mechanism to replicate behaviour of the nonlinear problem, giving more accuracy within the region of interest. Therefore, the main contribution of this work is that LP satisfies the successive derivatives in some specific points, resulting more accurate polynomials than Taylor expansion does for the same degree of their respective polynomials. Such characteristic makes of LPs a handy and powerful tool to approximate different kind of differential equations including: singular problems, initial condition and boundary-valued problems, equations with discontinuities, coupled differential equations, high-order equations, among others. Additionally, we show how the process to obtain the polynomials is straightforward and simple to implement; generating a compact, and easy to compute, expression. Even more, we present the process to approximate Gelfand's equation, an equation of an isothermal reaction, a model for chronic myelogenous leukemia, Thomas-Fermi equation, and a high order nonlinear differential equations with discontinuities getting, as result, accurate, fast and compact approximate solutions. In addition, we present the computational convergence and error studies for LPs resulting convergent polynomials and error tendency to zero as the order of LPs increases for all study cases. Finally, a study of CPU time shows that LPs require a few nano-seconds to be evaluated, which makes them suitable for intensive computing applications.
Project description:Analytic approximate solutions using Optimal Homotopy Perturbation Method (OHPM) are given for steady boundary layer flow over a nonlinearly stretching wall in presence of partial slip at the boundary. The governing equations are reduced to nonlinear ordinary differential equation by means of similarity transformations. Some examples are considered and the effects of different parameters are shown. OHPM is a very efficient procedure, ensuring a very rapid convergence of the solutions after only two iterations.
Project description:The paper presents the optimal homotopy perturbation method, which is a new method to find approximate analytical solutions for nonlinear partial differential equations. Based on the well-known homotopy perturbation method, the optimal homotopy perturbation method presents an accelerated convergence compared to the regular homotopy perturbation method. The applications presented emphasize the high accuracy of the method by means of a comparison with previous results.
Project description:An analytic study was conducted on coupled partial differential equations. We formally derived new solitary wave solutions of generalized coupled system of Zakharov-Kuznetsov (ZK) and KdV equations by using modified extended tanh method. The traveling wave solutions for each generalized coupled system of ZK and KdV equations are shown in form of periodic, dark, and bright solitary wave solutions. The structures of the obtained solutions are distinct and stable.
Project description:Continuum solvation modeling based upon the Poisson-Boltzmann equation (PBE) is widely used in structural and functional analysis of biomolecules. In this work, we propose a charge-central interpretation of the full nonlinear PBE electrostatic interactions. The validity of the charge-central view or simply charge view, as formulated as a vacuum Poisson equation with effective charges, was first demonstrated by reproducing both electrostatic potentials and energies from the original solvated full nonlinear PBE. There are at least two benefits when the charge-central framework is applied. First the convergence analyses show that the use of polarization charges allows a much faster converging numerical procedure for electrostatic energy and forces calculation for the full nonlinear PBE. Second, the formulation of the solvated electrostatic interactions as effective charges in vacuum allows scalable algorithms to be deployed for large biomolecular systems. Here, we exploited the charge-view interpretation and developed a particle-particle particle-mesh (P3M) strategy for the full nonlinear PBE systems. We also studied the accuracy and convergence of solvation forces with the charge-view and the P3M methods. It is interesting to note that the convergence of both the charge-view and the P3M methods is more rapid than the original full nonlinear PBE method. Given the developments and validations documented here, we are working to adapt the P3M treatment of the full nonlinear PBE model to molecular dynamics simulations.
Project description:The bidomain equations have been widely used to model the electrical activity of cardiac tissue. While it is well-known that implicit methods have much better stability than explicit methods, implicit methods usually require the solution of a very large nonlinear system of equations at each timestep which is computationally prohibitive. In this work, we present two fully implicit time integration methods for the bidomain equations: the backward Euler method and a second-order one-step two-stage composite backward differentiation formula (CBDF2) which is an L-stable time integration method. Using the backward Euler method as fundamental building blocks, the CBDF2 scheme is easily implementable. After solving the nonlinear system resulting from application of the above two fully implicit schemes by a nonlinear elimination method, the obtained nonlinear global system has a much smaller size, whose Jacobian is symmetric and possibly positive definite. Thus, the residual equation of the approximate Newton approach for the global system can be efficiently solved by standard optimal solvers. As an alternative, we point out that the above two implicit methods combined with operator splittings can also efficiently solve the bidomain equations. Numerical results show that the CBDF2 scheme is an efficient time integration method while achieving high stability and accuracy.
Project description:In this research article, we investigated a comprehensive analysis of time-dependent free convection electrically and thermally conducted water-based nanofluid flow containing Copper and Titanium oxide (Cu and TiO 2 ) past a moving porous vertical plate. A uniform transverse magnetic field is imposed perpendicular to the flow direction. Thermal radiation and heat sink terms are included in the energy equation. The governing equations of this flow consist of partial differential equations along with some initial and boundary conditions. The solution method of these flow interpreting equations comprised of two parts. Firstly, principal equations of flow are symmetrically transformed to a set of nonlinear coupled dimensionless partial differential equations using convenient dimensionless parameters. Secondly, the Laplace transformation technique is applied to those non-dimensional equations to get the close form exact solutions. The control of momentum and heat profile with respect to different associated parameters is analyzed thoroughly with the help of graphs. Fluid accelerates with increasing Grashof number (Gr) and porosity parameter (K), while increasing values of heat sink parameter (Q) and Prandtl number (Pr) drop the thermal profile. Moreover, velocity and thermal profile comparison for Cu and TiO 2 -based nanofluids is graphed.