A Numerical Method for Computing Time‐Periodic Solutions in Dissipative Wave Systems

A numerical method is proposed for computing time‐periodic and relative time‐periodic solutions in dissipative wave systems. In such solutions, the temporal period, and possibly other additional internal parameters such as the propagation constant, are unknown priori and need to be determined along with the solution itself. The main idea of the method is to first express those unknown parameters in terms of the solution through quasi‐Rayleigh quotients, so that the resulting integrodifferential equation is for the time‐periodic solution only. Then this equation is computed in the combined spatiotemporal domain as a boundary value problem by Newton‐conjugate‐gradient iterations. The proposed method applies to both stable and unstable time‐periodic solutions; its numerical accuracy is spectral; it is fast‐converging; its memory use is minimal; and its coding is short and simple. As numerical examples, this method is applied to the Kuramoto–Sivashinsky equation and the cubic‐quintic Ginzburg–Landau equation, whose time‐periodic or relative time‐periodic solutions with spatially periodic or spatially localized profiles are computed. This method also applies to systems of ordinary differential equations, as is illustrated by its simple computation of periodic orbits in the Lorenz equations. MATLAB codes for all numerical examples are provided in the Appendices to illustrate the simple implementation of the proposed method.     

Weighted pseudo‐almost periodic functions on time scales with applications to cellular neural networks with discrete delays

In this paper, we first propose a concept of weighted pseudo‐almost periodic functions on time scales and study some basic properties of weighted pseudo‐almost periodic functions on time scales. Then, we establish some results about the existence of weighted pseudo‐almost periodic solutions to linear dynamic equations on time scales. Finally, as an application of our results, we study the existence and global exponential stability of weighted pseudo‐almost periodic solutions for a class of cellular neural networks with discrete delays on time scales. The results of this paper are completely new. Copyright © 2016 John Wiley & Sons, Ltd.     

Conservative schemes for the symmetric Regularized Long Wave equations

In this paper, we study the Symmetric Regularized Long Wave (SRLW) equations by finite difference method. We design some numerical schemes which preserve the original conservative properties for the equations. The first scheme is two-level and nonlinear-implicit. Existence of its difference solutions are proved by Brouwer fixed point theorem. It is proved by the discrete energy method that the scheme is uniquely solvable, unconditionally stable and second-order convergent for U in L_∞ norm, and for N in L₂ norm on the basis of the priori estimates. The second scheme is three-level and linear-implicit. Its stability and second-order convergence are proved. Both of the two schemes are conservative so can be used for long time computation. However, they are coupled in computing so need more CPU time. Thus we propose another three-level linear scheme which is not only conservative but also uncoupled in computation, and give the numerical analysis on it. Numerical experiments demonstrate that the schemes are accurate and efficient.     

Integrable Discrete Model for One‐Dimensional Soil Water Infiltration

We propose an integrable discrete model of one‐dimensional soil water infiltration. This model is based on the continuum model by Broadbridge and White, which takes the form of nonlinear convection–diffusion equation with a nonlinear flux boundary condition at the surface. It is transformed to the Burgers equation with a time‐dependent flux term by the hodograph transformation. We construct a discrete model preserving the underlying integrability, which is formulated as the self‐adaptive moving mesh scheme. The discretization is based on linearizability of the Burgers equation to the linear diffusion equation, but the naïve discretization based on the Euler scheme which is often used in the theory of discrete integrable systems does not necessarily give a good numerical scheme. Taking desirable properties of a numerical scheme into account, we propose an alternative discrete model that produces solutions with similar accuracy to direct computation on the original nonlinear equation, but with clear benefits regarding computational cost.     

Analysis of algebraic systems arising from fourth‐order compact discretizations of convection‐diffusion equations

We study the properties of coefficient matrices arising from high‐order compact discretizations of convection‐diffusion problems. Asymptotic convergence factors of the convex hull of the spectrum and the field of values of the coefficient matrix for a one‐dimensional problem are derived, and the convergence factor of the convex hull of the spectrum is shown to be inadequate for predicting the convergence rate of GMRES. For a two‐dimensional constant‐coefficient problem, we derive the eigenvalues of the nine‐point matrix, and we show that the matrix is positive definite for all values of the cell‐Reynolds number. Using a recent technique for deriving analytic expressions for discrete solutions produced by the fourth‐order scheme, we show by analyzing the terms in the discrete solutions that they are oscillation‐free for all values of the cell Reynolds number. Our theoretical results support observations made through numerical experiments by other researchers on the non‐oscillatory nature of the discrete solution produced by fourth‐order compact approximations to the convection‐diffusion equation. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 155–178, 2002; DOI 10.1002/num.1041     

Bifurcating periodic solution for a class of first-order nonlinear delay differential equation after Hopf bifurcation

In this paper, we predict the accurate bifurcating periodic solution for a general class of first-order nonlinear delay differential equation with reflectional symmetry by constructing an approximate technique, named residue harmonic balance. This technique combines the features of the homotopy concept with harmonic balance which leads to easy computation and gives accurate prediction on the periodic solution to the desired accuracy. The zeroth-order solution using just one Fourier term is applied by solving a set of nonlinear algebraic equations containing the delay term. The unbalanced residues due to Fourier truncation are considered iteratively by solving linear equations to improve the accuracy and increase the number of Fourier terms of the solutions successively. It is shown that the solutions are valid for a wide range of variation of the parameters by two examples. The second-order approximations of the periodic solutions are found to be in excellent agreement with those obtained by direct numerical integration. Moreover, the residue harmonic balance method works not only in determining the amplitude but also the frequency of the bifurcating periodic solution. The method can be easily extended to other delay differential equations.     

On the numerical analysis of finite element and dirichlet-to-neumann methods for nonlinear exterior transmission problems

We provide the numerical analysis of the combination of finite elements and Dirichlet-to-Neumann mappings (based on boundary integral operators) for a class of nonlinear exterior transmission problems whose weak formulations reduce to Lipschitz-continuous and strongly monotone operator equations. As a model we consider a nonlinear second order elliptic equation in divergence form in a bounded inner region of the plane, coupled with the Laplace equation in the corresponding unbounded exterior part. A discrete Galerkin scheme is presented by using linear finite elements on a triangulation of the domain, and then applying numerical quadrature and analytical formulae to evaluate all the linear, bilinear and semilinear forms involved. We prove the unique solvability of the discrete equations, and show the strong convergence of the approximate solutions. Furthermore, assuming additional regularity on the solution of the continuous operator equation, the asymptotic rate of convergence O(h) is also derived. Finally, numerical experiments are presented, which confirm the convergence results.     

Computation of periodic solutions in maximal monotone dynamical systems with guaranteed consistency

In this paper, we study a class of set-valued dynamical systems that satisfy maximal monotonicity properties. This class includes linear relay systems, linear complementarity systems, and linear mechanical systems with dry friction under some conditions. We discuss two numerical schemes based on time-stepping methods for the computation of the periodic solutions when these systems are periodically excited. We provide formal mathematical justifications for the numerical schemes in the sense of consistency, which means that the continuous-time interpolations of the numerical solutions converge to the continuous-time periodic solution when the discretization step vanishes. The two time-stepping methods are applied for the computation of the periodic solution exhibited by a power electronic converter and the corresponding methods are compared in terms of approximation accuracy and computation time.     

Conservative compact finite difference scheme for the N‐coupled nonlinear Klein–Gordon equations

In this article, a compact finite difference method is developed for the periodic initial value problem of the N‐coupled nonlinear Klein–Gordon equations. The present scheme is proved to preserve the total energy in the discrete sense. Due to the difficulty in obtaining the priori estimate from the discrete energy conservation law, the cut‐off function technique is employed to prove the convergence, which shows the new scheme possesses second order accuracy in time and fourth order accuracy in space, respectively. Additionally, several numerical results are reported to confirm our theoretical analysis. Lastly, we apply the reliable method to simulate and study the collisions of solitary waves numerically.     

Numerical solution of the free‐surface viscous flow on a horizontal rotating elliptical cylinder

The numerical solution of the free‐surface fluid flow on a rotating elliptical cylinder is presented. Up to the present, research has concentrated on the circular cylinder for which steady solutions are the main interest. However, for noncircular cylinders, such as the ellipse, steady solutions are no longer possible, but there will be periodic solutions in which the solution is repeated after one full revolution of the cylinder. It is this new aspect that makes the investigation of noncircular cylinders novel. Here we consider both the time‐dependent and periodic solutions for zero Reynolds number fluid flow. The numerical solution is expedited by first mapping the fluid film domain onto a rectangle such that the position of the free‐surface is determined as part of the solution. For the time‐dependent case a simple time‐marching method of lines approach is adopted. For the periodic solution the discretised nonlinear equations have to be solved simultaneously over a time period. The resulting large system of equations is solved using Newton's method in which the form of the Jacobian enables a straightforward decomposition to be implemented, which makes matrix inversion manageable. In the periodic case all derivatives have been approximated pseudospectrally with the time derivative approximated by a differentiation matrix which has been specially derived so that the weight of fluid is algebraically conserved. Of interest is the solution for which the weight of fluid is at its maximum possible value, and this has been obtained by increasing the weight until a consistency break‐down occurs. Time‐dependent solutions do not produce the periodic solution after a long time‐scale but have protuberances which are constantly appearing and disappearing. Periodic solutions exhibit spectral accuracy solutions and maximum supportable weight solutions have been obtained for ranges of eccentricity and angular velocity. The maximum weights are less than and approximately proportional to those obtained for the circular case. The shapes of maximum weight solutions is distinctly different from sub‐maximum weight solutions. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Entropy dissipative one‐leg multistep time approximations of nonlinear diffusive equations

New one‐leg multistep time discretizations of nonlinear evolution equations are investigated. The main features of the scheme are the preservation of the non‐negativity and the entropy dissipation structure of the diffusive equations. The key ideas are to combine Dahlquist's G‐stability theory with entropy dissipation methods and to introduce a nonlinear transformation of variables, which provides a quadratic structure in the equations. It is shown that G‐stability of the one‐leg scheme is sufficient to derive discrete entropy dissipation estimates. The general result is applied to a cross‐diffusion system from population dynamics and a nonlinear fourth‐order quantum diffusion model, for which the existence of semidiscrete weak solutions is proved. Under some assumptions on the operator of the evolution equation, the second‐order convergence of solutions is shown. Moreover, some numerical experiments for the population model are presented, which underline the theoretical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1119–1149, 2015     

粘弹性方程全离散化有限体积元格式及数值模拟

本文利用有限体积元方法研究二维粘弹性方程, 给出一种时间二阶精度的全离散化有限体积元格式, 并给出这种全离散化有限体积元解的误差估计, 最后用数值例子验证数值结果与理论结果是相吻合的. 通过与有限元方法和有限差分方法相比较, 进一步说明了全离散化有限体积元格式是求解二维粘弹性方程数值解的最有效方法之一.     

Convergence of an unconditionally stable semi-implicit scheme for linear incompressible MHD equations

Semi-implicit methods have been introduced by Harned et al. to solve magneto-hydrodynamic equations (MHD) with numerical schemes which are unconditionally stable with respect to fast and shear Alfven modes. They prove the stability of their scheme for linear ideal MHD equations with periodic boundary conditions, and with some technical assumptions. In this paper, we prove convergence of the numerical approximation (time discretization), under the same hypothesis, but looking for solutions on any regular bounded open set of R³ with appropriate boundary conditions, and introducing finite resistivity and viscosity.     

Convergence analysis for parallel‐in‐time solution of hyperbolic systems

Parallel‐in‐time algorithms have been successfully employed for reducing time‐to‐solution of a variety of partial differential equations, especially for diffusive (parabolic‐type) equations. A major failing of parallel‐in‐time approaches to date, however, is that most methods show instabilities or poor convergence for hyperbolic problems. This paper focuses on the analysis of the convergence behavior of multigrid methods for the parallel‐in‐time solution of hyperbolic problems. Three analysis tools are considered that differ, in particular, in the treatment of the time dimension: (a) space–time local Fourier analysis, using a Fourier ansatz in space and time; (b) semi‐algebraic mode analysis, coupling standard local Fourier analysis approaches in space with algebraic computation in time; and (c) a two‐level reduction analysis, considering error propagation only on the coarse time grid. In this paper, we show how insights from reduction analysis can be used to improve feasibility of the semi‐algebraic mode analysis, resulting in a tool that offers the best features of both analysis techniques. Following validating numerical results, we investigate what insights the combined analysis framework can offer for two model hyperbolic problems, the linear advection equation in one space dimension and linear elasticity in two space dimensions.     

Two new approaches to computing Hopf bifurcation problems

We consider the computation of Hopf bifurcation for ordinary differential equations. Two new extended systems are given for the calculation of Hopf bifurcation problems: the first is composed of differential-algebraic equations with index 1, the other consists of differential equations by using a symmetry inherited from the autonomous system of ordinary differential equations. Both methods are especially suitable for calculating bifurcating periodic solutions since they transform the Hopf bifurcation problem into regular nonlinear boundary value problems which are very easy to implement. The bifurcation solutions become isolated solutions of the extended system so that our methods work both in the subcritical and supercritical case. The extended systems are based on an additional parameter ε; practical experience shows that one gets convergence for ε sufficiently large so that a substantial part of the bifurcating branch can be computed. The two methods are illustrated by numerical examples and compared with other procedures.     

Stability and convergence of optimum spectral non‐linear Galerkin methods

Our objective in this article is to present some numerical schemes for the approximation of the 2‐D Navier–Stokes equations with periodic boundary conditions, and to study the stability and convergence of the schemes. Spatial discretization can be performed by either the spectral Galerkin method or the optimum spectral non‐linear Galerkin method; time discretization is done by the Euler scheme and a two‐step scheme. Our results show that under the same convergence rate the optimum spectral non‐linear Galerkin method is superior to the usual Galerkin methods. Finally, numerical example is provided and supports our results. Copyright © 2001 John Wiley & Sons, Ltd.     

Conjugate gradient techniques for the optimal control evolution dam problem

This paper considers the numerical simulation of optimal control evolution dam problem by using conjugate gradient method.The paper considers the free boundary value problem related to time dependent fluid flow in a homogeneous earth rectangular dam.The dam is taken to be sufficiently long that the flow is considered to be two dimensional.On the left and right walls of the dam there is a reservoir of fluid at a level dependent on time.This problem can be transformed into a variational inequality on a fixed domain.The numerical techniques we use are based on a linear finite element method to approximate the state equations and a conjugate gradient algorithm to solve the discrete optimal control problem.This algorithm is based on Armijo's rule in the unconstrained optimization theory.The convergence of the discrete optimal solutions to the continuous optimal solutions,and the convergence of the conjugate gradient algorithm are proved.A numerical example is given to determine the location of the minimum surface     

Finding all Solutions of Systems of Nonlinear Equations Using the Dual Simplex Method

An efficient algorithm is proposed for finding all solutions of nonlinear equations using linear programming (LP). This algorithm is based on a simple test (termed the LP test) for nonexistence of a solution to a system of nonlinear equations in a given region. In the conventional LP test, the system of nonlinear equations is transformed into an LP problem, to which the simplex method is applied. However, although the LP test is very powerful, it requires many pivotings for each region. In this paper, we use the dual simplex method in the LP test, which makes the average number of pivotings per region much smaller (less than one, for example) and makes the algorithm very efficient. By numerical examples, it is shown that the proposed algorithm can find all solutions of systems of 200 nonlinear equations in practical computation time.     

Explicit periodic travelling waves for a discrete lambda–omega reaction–diffusion system

In 1973, Kopell and Howard introduced a λ–ω reaction–diffusion system and found an explicit family of periodic travelling wave solutions lying on circles with radius less than 1. Since λ–ω systems represent universal models for studying chemical processes, and onset of turbulent behaviour, etc., explicit solutions of λ–ω systems with delays or discrete λ–ω systems can be of further help when the only method for obtaining other solutions is through numerical computation. There are now much investigations of various λ–ω systems. However, it is of interest to note that none attempts to find explicit travelling wave solutions. In this paper, we investigate the existence of explicit solutions for the simplest Euler scheme of a λ–ω system with delays or advancements which is described as a coupled pair of partial difference equations. We are able to provide necessary as well as sufficient conditions for the existence of numerical periodic travelling wave solutions. Additionally, we also provide some examples to show that our explicit solutions are qualitatively different from those found by Kopell and Howard and hence they may be of interests for specialists in the area of reaction–diffusion systems.     

A Jacobi-collocation method for solving second kind Fredholm integral equations with weakly singular kernels

In this work, we propose a Jacobi-collocation method to solve the second kind linear Fredholm integral equations with weakly singular kernels. Particularly, we consider the case when the underlying solutions are sufficiently smooth. In this case, the proposed method leads to a fully discrete linear system. We show that the fully discrete integral operator is stable in both infinite and weighted square norms. Furthermore, we establish that the approximate solution arrives at an optimal convergence order under the two norms. Finally, we give some numerical examples, which confirm the theoretical prediction of the exponential rate of convergence.