The Loewner driving function of trajectory arcs of quadratic differentials

We obtain a first order differential equation for the driving function of the chordal Loewner differential equation in the case where the domain is slit by a curve which is a trajectory arc of a certain type of quadratic differential. In particular this includes the case when the curve is a path on the square, triangle or hexagonal lattice in the upper half-plane or, indeed, in any domain with boundary on the lattice. We also demonstrate how we use this to calculate the driving function numerically. Equivalent results for other variants of the Loewner differential equation are also obtained.     

The Power Matrix, coadjoint action and quadratic differentials

The coefficients of a quadratic differential which is changing under the Loewner flow satisfy a well-known differential system studied by Schiffer, Schaeffer and Spencer, and others. By work of Roth, this differential system can be interpreted as Hamilton's equations. We apply the power matrix to interpret this differential system in terms of the coadjoint action of the matrix group on the dual of its Lie algebra. As an application, we derive a set of integral invariants of Hamilton's equations which is in a certain sense complete. In function theoretic terms, these are expressions in the coefficients of the quadratic differential and Loewner map which are independent of the parameter in the Loewner flow.     

A parareal approach of semi‐linear parabolic equations based on general waveform relaxation

We present a parareal approach of semi‐linear parabolic equations based on general waveform relaxation (WR) at the partial differential equation (PDE) level. An algorithm for initial‐boundary value problem and two algorithms for time‐periodic boundary value problem are constructed. The convergence analysis of three algorithms are provided. The results show that the algorithm for initial‐boundary value problem is superlinearly convergent while both algorithms for the time‐periodic boundary value problem linearly converge to the exact solutions at most. Numerical experiments show that the parareal algorithms based on general WR at the PDE level, compared with the parareal algorithm based on the classical WR at the ordinary differential equations (ODEs) level (the PDEs is discretized into ODEs), require much fewer number of iterations to converge.     

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.     

Structured populations with distributed recruitment: from PDE to delay formulation

In this work, first, we consider a physiologically structured population model with a distributed recruitment process. That is, our model allows newly recruited individuals to enter the population at all possible individual states, in principle. The model can be naturally formulated as a first‐order partial integro‐differential equation, and it has been studied extensively. In particular, it is well posed on the biologically relevant state space of Lebesgue integrable functions. We also formulate a delayed integral equation (renewal equation) for the distributed birth rate of the population. We aim to illustrate the connection between the partial integro‐differential and the delayed integral equation formulation of the model utilising a recent spectral theoretic result. In particular, we consider the equivalence of the steady state problems in the two different formulations, which then lead us to characterise irreducibility of the semigroup governing the linear partial integro‐differential equation. Furthermore, using the method of characteristics, we investigate the connection between the time‐dependent problems. In particular, we prove that any (non‐negative) solution of the delayed integral equation determines a (non‐negative) solution of the partial differential equation and vice versa. The results obtained for the particular distributed states at birth model then lead us to present some very general results, which establish the equivalence between a general class of partial differential and delay equation, modelling physiologically structured populations. Copyright © 2016 John Wiley & Sons, Ltd.     

On the asymptotic positivity of solutions for the extended Fisher–Kolmogorov equation with nonlinear diffusion

The objective of this paper aims to prove positivity of solutions for a semilinear dissipative partial differential equation with non‐linear diffusion. The equation is a generalized model of the well‐known Fisher–Kolmogorov equation and represents a class of dissipative partial differential equations containing differential operators of higher order than the Laplacian. It arises in a variety of meaningful physical situations including gas flows, diffusion of an electron–ion plasma and the dynamics of biological populations whose mobility is density dependent. In all these situations, the solutions of the equation must be positive functions. Copyright © 2002 John Wiley & Sons, Ltd.     

Bernstein Ritz‐Galerkin method for solving an initial‐boundary value problem that combines Neumann and integral condition for the wave equation

In this article, the Ritz‐Galerkin method in Bernstein polynomial basis is implemented to give an approximate solution of a hyperbolic partial differential equation with an integral condition. We will deal here with a type of nonlocal boundary value problem, that is, the solution of a hyperbolic partial differential equation with a nonlocal boundary specification. The nonlocal conditions arise mainly when the data on the boundary cannot be measured directly. The properties of Bernstein polynomial and Ritz‐Galerkin method are first presented, then Ritz‐Galerkin method is used to reduce the given hyperbolic partial differential equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique presented in this article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A multigrid compact finite difference method for solving the one‐dimensional nonlinear sine‐Gordon equation

The aim of this paper is to propose a multigrid method to obtain the numerical solution of the one‐dimensional nonlinear sine‐Gordon equation. 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 compact finite difference scheme of fourth‐order for discretizing the spatial derivative and the standard second‐order central finite difference method for the time derivative. The proposed method uses the Richardson extrapolation method in time variable. The obtained system has been solved by V‐cycle multigrid (VMG) method, where the VMG method is used for solving the large sparse linear systems. The numerical examples show the efficiency of this algorithm for solving the one‐dimensional sine‐Gordon equation. Copyright © 2014 John Wiley & Sons, Ltd.     

A FINITE ELEMENT COLLOCATION METHOD FOR TWO-PHASE INCOMPRESSIBLE IMMISCIBLE PROBLEMS

Two-phase,incompressible,immiscible flow in porous media is governed by a coupled system of nonlinear partial differential equations.The pressure equation is elliptic, whereas the concentration equation is parabolic,and both are treated by the collocation scheme.Existence and uniqueness of solutions of the algorithm are proved.A optimal convergence analysis is given for the method.     

An optimal‐order estimate for MMOC‐MFEM approximations to porous medium flow

Mathematical models used to describe porous medium flow lead to coupled systems of time‐dependent partial differential equations. Standard methods tend to generate numerical solutions with nonphysical oscillations or numerical dispersion along with spurious grid‐orientation effect. The MMOC‐MFEM time‐stepping procedure, in which the modified method of characteristics (MMOC) is used to solve the transport equation and a mixed finite element method (MFEM) is used for the pressure equation, simulates porous medium flow accurately even if large spatial grids and time steps are used. In this article we prove an optimal‐order error estimate for a family of MMOC‐MFEM approximations. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

李群方法在渗流力学的应用

试图用李群方法来分析流体及渗流的运动规律.对于流形上流体、渗流力学方程的研究,物理空间的流动中的拓扑结构只要具有李群的性质,便可以此来进行流动分析.这是将李群理论直接、直地应用于渗流力学的一种方法.李群方法将众多求解特定类型的渗流微分方程方法统一到共同的概念之下.李群无穷小变换方法为寻找微分方程的闭合形式的解提供的广泛的应用,补充了求解渗流力学方程的数学物理技巧.     

New approach to the Lax‐Wendroff modified differential equation for linear and nonlinear advection

The paper presents an enhanced analysis of the Lax‐Wendroff difference scheme—up to the eighth‐order with respect to time and space derivatives—of the modified‐partial differential equation (MDE) of the constant‐wind‐speed advection equation. The modified equation has been so far derived mainly as a fourth‐order equation. The Π ‐form of the first differential approximation (differential approximation or equivalent equation) derived by expressing the time derivatives in terms of the space derivatives is used for presenting the MDE. The obtained coefficients at higher order derivatives are analyzed for indications of the character of the dissipative and dispersive errors. The authors included a part of the stencil applied for determining the modified differential equation up to the eighth‐order of the analyzed modified differential equation for the second‐order Lax‐Wendroff scheme. Neither the derived coefficients at the space derivatives of order p ∈ (7 – 8) in the modified differential equation for the Lax‐Wendroff difference scheme nor the results of analyses on the basis of these coefficients of the group velocity, phase shift errors, or dispersive and dissipative features of the scheme have been published. The MDEs for 2 two‐step variants of the Lax‐Wendroff type difference schemes and the MacCormack predictor–corrector scheme (see MacCormack's study) constructed for the scalar hyperbolic conservation laws are also presented in this paper. The analysis of the inviscid Burgers equation solution with the initial condition in a form of a shock wave has been discussed on their basis. The inviscid Burgers equation with the source is also presented. The theory of MDE started to develop after the paper of C. W. Hirt was published in 1968.     

An algorithmic method for showing existence of nontrivial non-classical symmetries of partial differential equations without solving determining equations

In this paper, based on differential characteristic set theory and the associated algorithm (also called Wu?s method), an algorithmic method is presented to decide on the existence of a nontrivial non-classical symmetry of a given partial differential equation without solving the corresponding nonlinear determining system. The theory and algorithm give a partial answer for the open problem posed by P.A. Clarkson and E.L. Mansfield in [21] on non-classical symmetries of partial differential equations. As applications of our algorithm, non-classical symmetries and corresponding invariant solutions are found for several evolution equations.     

扩充偏微分方程(组)守恒律和对称的辅助方程方法及微分形式吴方法的应用

首先,我们给出了引入伴随方程(组)扩充原方程(组)的策略使给定偏微分方程(组)的扩充方程组具有对应泛瓯即,成为Lagrange系统的方法,以此为基础提出了作为偏微分方程(组)传统守恒律和对称概念的一种推广-偏微分方程(组)扩充守恒律和扩充对称的概念;其次,以得到的Lagrange系统为基础给定了确定原方程(组)扩充守恒律和扩充对称的方法,从而达到扩充给定偏微分方程(组)的首恒律和对称的目的;第三,提出了适用于一般形式微分方程(组)的计算固有守恒律的方法;第四,实现以上算法过程中,我们先把计算(扩充)守恒律和对称问题均归结为求解超定线性齐次偏微分方程组(确定方程组)的问题.然后,对此关键问题我们提出了用微分形式吴方法处理的有效算法;最后,作为方法的应用我们计算确定了非线性电报方程组在内的五个发展方程(组)的新守恒律和对称,同时也说明了方法的有效性.     

Exact solutions of nonlinear time fractional partial differential equations by sub‐equation method

In this article, the sub‐equation method is presented for finding the exact solutions of a nonlinear fractional partial differential equations. For this, the fractional complex transformation method has been used to convert fractional‐order partial differential equation to ordinary differential equation. The fractional derivatives are described in Jumarie's the modified Riemann–Liouville sense. We apply to this method for the nonlinear time fractional differential equations. With the aid of symbolic computation, a variety of exact solutions for them are obtained. Copyright © 2014 John Wiley & Sons, Ltd.     

广义Boussinesq方程的多辛方法

广义Boussinesq方程作为一类重要的非线性方程有着许多有趣的性质,基于Hamilton空间体系的多辛理论研究了广义Boussinesq方程的数值解法,构造了一种等价于多辛Box格式的新隐式多辛格式,该格式满足多辛守恒律、局部能量守恒律和局部动量守恒律．对广义Boussinesq方程孤子解的数值模拟结果表明,该多辛离散格式具有较好的长时间数值稳定性．     

Fourth‐order compact scheme with local mesh refinement for option pricing in jump‐diffusion model

The value of a contingent claim under a jump‐diffusion process satisfies a partial integro‐differential equation. A fourth‐order compact finite difference scheme is applied to discretize the spatial variable of this equation. It is discretized in time by an implicit‐explicit method. Meanwhile, a local mesh refinement strategy is used for handling the nonsmooth payoff condition. Moreover, the numerical quadrature method is exploited to evaluate the jump integral term. It guarantees a Toeplitz‐like structure of the integral operator such that a fast algorithm is feasible. Numerical results show that this approach gives fourth‐order accuracy in space. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011