Algebraic dynamics solutions and algebraic dynamics algorithm for nonlinear ordinary differential equations

The problem of preserving fidelity in numerical computation of nonlinear ordinary differential equations is studied in terms of preserving local differential structure and approximating global integration structure of the dynamical system. The ordinary differential equations are lifted to the corresponding partial differential equations in the framework of algebraic dynamics, and a new algorithm—algebraic dynamics algorithm is proposed based on the exact analytical solutions of the ordinary differential equations by the algebraic dynamics method. In the new algorithm, the time evolution of the ordinary differential system is described locally by the time translation operator and globally by the time evolution operator. The exact analytical piece-like solution of the ordinary differential equations is expressed in terms of Taylor series with a local convergent radius, and its finite order truncation leads to the new numerical algorithm with a controllable precision better than Runge Kutta Algorithm and Symplectic Geometric Algorithm.     

A bubble-stabilized least-squares finite element method for steady MHD duct flow problems at high Hartmann numbers

In this paper we devise a stabilized least-squares finite element method using the residual-free bubbles for solving the governing equations of steady magnetohydrodynamic duct flow. We convert the original system of second-order partial differential equations into a first-order system formulation by introducing two additional variables. Then the least-squares finite element method using C⁰$C^0$ linear elements enriched with the residual-free bubble functions for all unknowns is applied to obtain approximations to the first-order system. The most advantageous features of this approach are that the resulting linear system is symmetric and positive definite, and it is capable of resolving high gradients near the layer regions without refining the mesh. Thus, this approach is possible to obtain approximations consistent with the physical configuration of the problem even for high values of the Hartmann number. Before incoorperating the bubble functions into the global problem, we apply the Galerkin least-squares method to approximate the bubble functions that are exact solutions of the corresponding local problems on elements. Therefore, we indeed introduce a two-level finite element method consisting of a mesh for discretization and a submesh for approximating the computations of the residual-free bubble functions. Numerical results confirming theoretical findings are presented for several examples including the Shercliff problem.     

Prescribing a multistage analytical method to a prey-predator dynamical system

This article discusses the effectiveness of a fresh analytical method in solving a prey-predator problem, which is described as a system of two nonlinear ordinary differential equations. The method of interest is the multistage variational iteration method (MVIM), which provides a slight modification of the classical variational iteration method (VIM). We shall compare solutions of the classical VIM along with MVIM and match them against the conventional numerical method, Runge-Kutta (RK4) (fourth-order).     

Rational Chebyshev pseudospectral approach for solving Thomas-Fermi equation

In this Letter we propose a pseudospectral method for solving Thomas-Fermi equation which is a nonlinear ordinary differential equation on semi-infinite interval. This approach is based on rational Chebyshev pseudospectral method. This method reduces the solution of this problem to the solution of a system of algebraic equations. Comparison with some numerical solutions shows that the present solution is highly accurate.     

Application of the Generalized Differential Quadrature Method in Solving Burgers'Equations

The aim of this paper is to obtain numerical solutions of the one-dimensional, two-dimensional and coupled Burgers' equations through the generalized differential quadrature method (GDQM). The polynomial-based differential quadrature (PDQ) method is employed and the obtained system of ordinary differential equations is solved via the total variation diminishing Runge-Kutta (TVD-RK) method. The numerical solutions are satisfactorily coincident with
the exact solutions. The method can compete against the methods applied in the literature.     

Rational Legendre pseudospectral approach for solving nonlinear differential equations of Lane–Emden type

Lane–Emden equation is a nonlinear singular equation in the astrophysics that corresponds to the polytropic models. In this paper, a pseudospectral technique is proposed to solve the Lane–Emden type equations on a semi-infinite domain. The method is based on rational Legendre functions and Gauss–Radau integration. The method reduces solving the nonlinear ordinary differential equation to solve a system of nonlinear algebraic equations. The comparison of the results with the other numerical methods shows the efficiency and accuracy of this method.     

Numerical methods for solving one-dimensional cochlear models in the time domain

In this article, a robust numerical solution method for one-dimensional (1-D) cochlear models in the time domain is presented. The method has been designed particularly for models with a cochlear partition having nonlinear and active mechanical properties. The model equations are discretized with respect to the spatial variable by means of the principle of Galerkin to yield a system of ordinary differential equations in the time variable. To solve this system, several numerical integration methods concerning stability and computational performance are compared. The selected algorithm is based on a variable step size fourth-order Runge-Kutta scheme; it is shown to be both more stable and much more efficient than previously published numerical solution techniques.     

A shifted Jacobi collocation algorithm for wave type equations with non-local conservation conditions

In this paper, we propose an efficient spectral collocation algorithm to solve numerically wave type equations subject to initial, boundary and non-local conservation conditions. The shifted Jacobi pseudospectral approximation is investigated for the discretization of the spatial variable of such equations. It possesses spectral accuracy in the spatial variable. The shifted Jacobi-Gauss-Lobatto (SJ-GL) quadrature rule is established for treating the non-local conservation conditions, and then the problem with its initial and non-local boundary conditions are reduced to a system of second-order ordinary differential equations in temporal variable. This system is solved by two-stage forth-order A-stable implicit RK scheme. Five numerical examples with comparisons are given. The computational results demonstrate that the proposed algorithm is more accurate than finite difference method, method of lines and spline collocation approach     

Birkhoffian Symplectic Scheme for a Quantum System

In this paper, a classical system of ordinary differential equations is built to describe a kind of n-dimensional quantum systems. The absorption spectrum and the density of the states for the system are defined from the points of quantum view and classical view. From the Birkhoffian form of the equations, a Birkhoffian symplectic scheme is derived for solving n-dimensional equations by using the generating function method. Besides the Birkhoffian structure- preserving, the new scheme is proven to preserve the discrete local energy conservation law of the system with zero vector f . Some numerical experiments for a 3-dimensional example show that the new scheme can simulate the general Birkhoffian system better than the implicit midpoint scheme, which is well known to be symplectic scheme for Hamiltonian system.     

两点边值问题解的一种有效算法

本文讨论非线性微分方程边值问题的数值解。对非线性项在局部予以线性化后,再应用打靶方法求解,可加快收敛过程;同时对线性化的函数值采取插值逼近,进一步减少了计算量。本文算法格式简便、编程容易。若辅助内、外存交换技术,利用本文算法,可在微机上完成较大规模复杂问题的分析。算例表明,本文算法大大快于用牛顿法求解一些差分格式方程的收敛速度。     

Numerical analysis of dispersion characteristics of flat dielectric waveguides

A method is proposed for calculation of the dispersion characteristics of flat gradient dielectric waveguides. The method is based on numerical integration of a system of ordinary differential equations obtained by differentiation of solutions of the initial boundary-value problem with respect to a parameter.Samara State University. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 36, No. 2, pp. 143–147, February, 1993.     

Comparative study of travelling-wave and numerical solutions for the coupled short pulse （CSP） equation

The Lie symmetry analysis is performed for the coupled short plus (CSP) equation. We derive the infinitesimals that admit the classical symmetry group. Five types arise depending on the nature of the Lie symmetry generator. In all types, we find reductions in terms of system of ordinary differential equations, and exact solutions of the CSP equation are derived, which are compared with numerical solutions using the classical fourth-order Runge-Kutta scheme.     

New exact travelling wave solutions of the generalized zakharov equations

A direct algebraic method is introduced for constructing exact travelling wave solutions of nonlinear partial differential equations with complex phases. The scheme is implemented for obtaining multiple soliton solutions of the generalized Zakharov equations, and then new exact travelling wave solutions with complex phases are obtained. In addition, by using new exact solutions of an auxiliary ordinary differential equation, new exact travelling wave solutions for the generalized Zakharov equations are obtained.     

A high-order cell-centered Lagrangian scheme for compressible fluid flows in two-dimensional cylindrical geometry

The goal of this paper is to present high-order cell-centered schemes for solving the equations of Lagrangian gas dynamics written in cylindrical geometry. A node-based discretization of the numerical fluxes is obtained through the computation of the time rate of change of the cell volume. It allows to derive finite volume numerical schemes that are compatible with the geometric conservation law (GCL). Two discretizations of the momentum equations are proposed depending on the form of the discrete gradient operator. The first one corresponds to the control volume scheme while the second one corresponds to the so-called area-weighted scheme. Both formulations share the same discretization for the total energy equation. In both schemes, fluxes are computed using the same nodal solver which can be viewed as a two-dimensional extension of an approximate Riemann solver. The control volume scheme is conservative for momentum, total energy and satisfies a local entropy inequality in its first-order semi-discrete form. However, it does not preserve spherical symmetry. On the other hand, the area-weighted scheme is conservative for total energy and preserves spherical symmetry for one-dimensional spherical flow on equi-angular polar grid. The two-dimensional high-order extensions of these two schemes are constructed employing the generalized Riemann problem (GRP) in the acoustic approximation. Many numerical tests are presented in order to assess these new schemes. The results obtained for various representative configurations of one and two-dimensional compressible fluid flows show the robustness and the accuracy of our new schemes.     

反应扩散模型在图灵斑图中的应用及数值模拟

反应扩散方程模型常被用于描述生物学中斑图的形成.从反应扩散模型出发,理论推导得到GiererMeinhardt模型的斑图形成机理,解释了非线性常微分方程系统的稳定常数平衡态在加入扩散项后会发生失稳并产生图灵斑图的过程.通过计算该模型,得到图灵斑图产生的参数条件.数值方法中采用一类有效的高精度数值格式,即在空间离散条件下采用Chebyshev谱配置方法,在时间离散条件下采用紧致隐积分因子方法.该方法结合了谱方法和紧致隐积分因子方法的优点,具有精度高、稳定性好、存储量小等优点.数值模拟表明,在其他条件一定的情况下,系统控制参数κ取不同值对于斑图的产生具有重要的影响,数值结果验证了理论结果.     

Solving coupled nonlinear Schrödinger equations via a direct discontinuous Galerkin method

In this work, we present the direct discontinuous Galerkin (DDG) method for the one-dimensional coupled nonlinear Schrödinger (CNLS) equation. We prove that the new discontinuous Galerkin method preserves the discrete mass conservations corresponding to the properties of the CNLS system. The ordinary differential equations obtained by the DDG space discretization is solved via a third-order stabilized Runge-Kutta method. Numerical experiments show that the new DDG scheme gives stable and less diffusive results and has excellent long-time numerical behaviors for the CNLS equations.     

High-Order Numerical Method for Two-Point Boundary Value Problems

In this paper the two-point boundary value problem is transformed into general first-order ordinary differential equation system through introduction of conditions of an integral character to supplement the simultaneous set of first-order equations. A new discrete approximation of a high-order compact difference scheme is presented for the first-order system. It is a block-bidiagonal profile and removes the limits of other high-order discrete schemes at the interval ends. The numerical tests of a seventh-order compact difference scheme show that the proposed scheme is very convenient and efficient for linear and nonlinear two-point boundary value problems.     

On bifurcations in a chaotically stirred excitable medium

We study a one-dimensional filamental model of a chaotically stirred excitable medium. In a numerical simulation we systematically explore its rich bifurcation scenarios involving saddle-nodes, Hopf bifurcations and hysteresis loops. The bifurcations are described in terms of two parameters signifying the excitability of the reacting medium and the strength of the chaotic stirring, respectively. The solution behaviour, in particular at the bifurcation points, is analytically described by means of a nonperturbative variational method. Using this method we reduce the partial differential equations to either algebraic equations for stationary solutions and bifurcations, or to ordinary differential equations in the case of non-stationary solutions and bifurcations. We present numerical simulations corroborating our analytical results.     

Domain decomposition solution of nonlinear two-dimensional parabolic problems by random trees

A domain decomposition method is developed for the numerical solution of nonlinear parabolic partial differential equations in any space dimension, based on the probabilistic representation of solutions as an average of suitable multiplicative functionals. Such a direct probabilistic representation requires generating a number of random trees, whose role is that of the realizations of stochastic processes used in the linear problems. First, only few values of the sought solution inside the space-time domain are computed (by a Monte Carlo method on the trees). An interpolation is then carried out, in order to approximate interfacial values of the solution inside the domain. Thus, a fully decoupled set of sub-problems is obtained. The algorithm is suited to massively parallel implementation, enjoying arbitrary scalability and fault tolerance properties. Pruning the trees is shown to increase appreciably the efficiency of the algorithm. Numerical examples conducted in 2D, including some for the KPP equation, are given.     

Application of He's homotopy perturbation method to boundary layer flow and convection heat transfer over a flat plate

In this Letter, the problem of forced convection over a horizontal flat plate is presented and the homotopy perturbation method (HPM) is employed to compute an approximation to the solution of the system of nonlinear differential equations governing on the problem. It has been attempted to show the capabilities and wide-range applications of the homotopy perturbation method in comparison with the previous ones in solving heat transfer problems. The obtained solutions, in comparison with the exact solutions admit a remarkable accuracy. A clear conclusion can be drawn from the numerical results that the HPM provides highly accurate numerical solutions for nonlinear differential equations.