Symmetries of compatibility conditions for systems of differential equations

When symmetries of differential equations are applied, various types of associated systems of equations appear. Compatibility conditions of the associated systems expressed in the form of differential equations inherit Lie symmetries of the initial equations. Invariant solutions to compatibility systems are known as orbits of partially invariant and generic solutions involved in the Lie group foliation of differential equations and so on. In some cases Bäcklund transformations and differential substitutions connecting quotient equations for compatibility conditions and initial systems naturally arise. Besides, Ovsiannikov's orbit method for finding partially invariant solutions is essentially based on such symmetries.     

Derivation of the equations of nonisothermal acoustics in elastic porous media

We consider the problem of the joint motion of a thermoelastic solid skeleton and a viscous thermofluid in pores, when the physical process lasts for a few dozens of seconds. These problems arise in describing the propagation of acoustic waves. We rigorously derive the homogenized equations (i.e., the equations not containing fast oscillatory coefficients) which are different types of nonclassical acoustic equations depending on relations between the physical parameters and the homogenized heat equation. The proofs are based on Nguetseng’s two-scale convergence method.     

On the weak discontinuity surfaces and tapes of equation systems in fluid mechanics

The weak discontinuity surfaces for a system of quasi-linear differential equations of higher order are developed. The classification of equation systems in fluid mechanics is based on the propagative weak discontinuity surfaces. Types of equations for different flow models are discussed. The conclusion is as follows:(a) For incompressible nonviscous flow, incompressible viscous flow and compressible viscous flow, the types of equations are all parabolic in the unsteady case and elliptic in the steady case.(b) For compressible nonviscous flow, the type of equations is hyperbolic in the unsteady case or steady supersonic case, and the type is elliptic in the steady subsonic case.     

A method for the numerical solution of the Painlevé equations

A numerical method for solving the Cauchy problem for all the six Painlevé equations is proposed. The difficulty of solving these equations is that the unknown functions can have movable (that is, dependent on the initial data) singular points of the pole type. Moreover, the Painlevé III–VI equations may have singularities at points where the solution takes certain finite values. The positions of all these singularities are not a priori known and are determined in the process of solving the equation. The proposed method is based on the transition to auxiliary systems of differential equations in neighborhoods of the indicated points. The equations in these systems and their solutions have no singularities at the corresponding point and its neighborhood. Such auxiliary equations are derived for all Painlevé equations and for all types of singularities. Efficient criteria for transition to auxiliary systems are formulated, and numerical results illustrating the potentials of the method are presented.     

A method for calculating the Painlevé transcendents

A numerical method for solving the Cauchy problem for all the six Painlevé equations is proposed. The difficulty of solving these equations is that the unknown functions can have movable (that is, dependent on the initial data) singular points of the pole type. Moreover, the Painlevé III–VI equations may have singularities at points where the solution takes certain finite values. The positions of all these singularities are not a priori known and are determined in the process of solving the equation. The proposed method is based on the transition to auxiliary systems of differential equations in neighborhoods of the indicated points. The equations in these systems and their solutions have no singularities at the corresponding point and its neighborhood. Such auxiliary equations are derived for all Painlevé equations and for all types of singularities. Efficient criteria for transition to auxiliary systems are formulated, and numerical results illustrating the potentials of the method are presented.     

Lie symmetry analysis, optimal systems and exact solutions to the fifth-order KdV types of equations

In this paper, the Lie symmetry analysis is performed on the fifth-order KdV types of equations which arise in modeling many physical phenomena. The similarity reductions and exact solutions are obtained based on the optimal system method. Then the exact analytic solutions are considered by using the power series method.     

（２＋１）维非线性色散长波方程的相似约化和解析解

首先借助于Ｍａｔｈｅｍａｔｉｃａ软件,将Ｃｌａｒｋｓｏｎ和Ｋｒｕｓｋａｌ引入的直接约化法推广并应用于（２＋１） 维偏微分方程组情形 （２＋１） 维非线性色散长波方程,获得了该方程的六种类型的相似约化和若干解析解,其中包括ＰａｉｎｌｅｖｅⅡ型方程和孤子解．然后基于文［５］的结论,通过引入新的级数变换,获得了该方程的有理分式解析解．这种方法也适合于其它的微分方程．     

WEB-Spline–Based Multigrid Methods for the Stationary Stokes Problem

We discuss a multigrid technique in solving a large system of linear algebraic equations arising in the approximation of Stokes equations by a new strategy based on weighted extended B-spline (WEB-spline) methods. Three types of WEB-spline–based Stokes elements satisfying the inf-sup condition are considered. First for a linear-constant type of Stokes element, we give the detailed multigrid algorithm and its convergence proof. The convergence proof of the multigrid algorithm for a bubble-stabilized WEB-spline–based Stokes element is dealt with separately. Multigrid method in the case of bubble-condensed variational form is simplified using the techniques from the bubble-stabilized case.     

Non-monotone traveling waves and entire solutions for a delayed nonlocal dispersal equation

This paper is concerned with the traveling waves and entire solutions for a delayed nonlocal dispersal equation with convolution- type crossing-monostable nonlinearity. We first establish the existence of non-monotone traveling waves. By Ikehara’s Tauberian theorem, we further prove the asymptotic behavior of traveling waves, including monotone and non-monotone ones. Then, based on the obtained asymptotic behavior, the uniqueness of the traveling waves is proved. Finally, the entire solutions are considered. By introducing two auxiliary monostable equations and establishing some comparison arguments for the three equations, some new types of entire solutions are constructed via the traveling wavefronts and spatially independent solutions of the auxiliary equations.     

A Novel Numerical Approach to Time-Fractional Parabolic Equations with Nonsmooth Solutions

This paper is concerned with numerical solutions of time-fractional parabolic equations. Due to the Caputo time derivative being involved, the solutions of equations are usually singular near the initial time $t = 0$ even for a smooth setting. Based on a simple change of variable $s = t^β$, an equivalent $s$-fractional differential equation is derived and analyzed. Two types of finite difference methods based on linear and quadratic approximations in the $s$-direction are presented, respectively, for solving the $s$-fractional differential equation. We show that the method based on the linear approximation provides the optimal accuracy$\mathcal{O}(N ^{−(2−α)})$ where $N$ is the number of grid points in temporal direction. Numerical examples for both linear and nonlinear fractional equations are presented in comparison with $L1$ methods on uniform meshes and graded meshes, respectively. Our numerical results show clearly the accuracy and efficiency of the proposed methods.     

The numerical solution of an evolution problem of second order in time on a closed smooth boundary

We consider an initial value problem for the second-order differential equation with a Dirichlet-to-Neumann operator coefficient. For the numerical solution we carry out semi-discretization by the Laguerre transformation with respect to the time variable. Then an infinite system of the stationary operator equations is obtained. By potential theory, the operator equations are reduced to boundary integral equations of the second kind with logarithmic or hypersingular kernels. The full discretization is realized by Nyström's method which is based on the trigonometric quadrature rules. Numerical tests confirm the ability of the method to solve these types of nonstationary problems.     

Predictor-corrector pseudospectral methods for stochastic partial differential equations with additive white noise

Commonly used finite-difference numerical schemes show some deficiencies in the integration of certain types of stochastic partial differential equations with additive white noise. In this paper efficient predictor-corrector spectral schemes to integrate these equations are discussed. They are all based on the discretization of the system in Fourier space. The nonlinear terms are treated using a pseudospectral approach so as to speed up the computations without a significant loss of accuracy. The proposed schemes are applied to solve, both in one and two spatial dimensions, two paradigmatic continuum models arising in the context of nonequilibrium dynamics of growing interfaces: the Kardar-Parisi-Zhang and Lai-Das Sarma-Villain equations. Numerical results about the Lai-Das Sarma-Villain equation in two spatial dimensions have not been previously reported in the literature.     

Diagram equations of the theory of fully developed turbulence

Skeleton diagram equations of turbulence theory — the Dyson equations and the equations for vertices of three types — are obtained nonperturbatively. Their derivation is based on the use of an equation in functional derivatives for the characteristic functional of a hydrodynamic system described by Navier-Stokes equations in the presence of an external random force. The iterative solution of these equations reproduces the perturbation series for second moments that is usually obtained in a more complicated way and also the series for the third moments.Institute of Problems in Mechanics, Russian Academy of Sciences. Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 101, No. 1, pp. 28–37, October, 1994.     

Reverse-time continuation of elastic waves

Methods based on downward continuation of an oscillation field have received wide recognition in the processing of multichannel seismic prospecting data. Physically, the idea behind this approach is clear: a wavefield observed on some surface is continued into the medium and backward in time. Mathematically, all the continuation algorithms used are based on scalar wave equation modeling, which takes proper account of the wave nature of oscillations of some types, but defies consideration of their vector nature. It is well known that a system of equations of dynamic elasticity theory (Lame equations) is a more adequate model for describing seismic oscillations. In this paper, a continuation of a field of elastic oscillations in an inhomogeneous isotropic medium is considered.     

Instability of parametrically second- and third-subharmonic resonances governed by nonlinear Shrödinger equations with complex coefficients

A theoretical analysis of the parametric harmonic response of two resonant modes is made based on a cubic nonlinear system. The analysis based on the method of multiple scales. Two types of the modified nonlinear Schrödinger equations with complex coefficients are derived to govern the resonance wave. One of these equations contains the first derivatives in space for a complex-conjugate type as well as a linear complex-conjugate term that is valid in the second-harmonic resonance cases. The second parametric equation contains a complex-conjugate type which is valid at the third-subharmonic resonance case. Estimates of nonlinear coefficients are made. The resulting equations have an interesting in many dynamical and physical cases. Temporal modulational method is confirmed to discuss the stability behavior at both parametric second- and third-harmonic resonance cases. Furthermore, the Benjamin–Feir instability is discussed for the sideband perturbation. The instability behavior at the sharp resonance is examined and the existence of the instability is found.     

Construction of splitting schemes based on transition operator approximation

The stability analysis of approximate solutions to unsteady problems for partial differential equations is usually based on the use of the canonical form of operator-difference schemes. Another possibility widely used in the analysis of methods for solving Cauchy problems for systems of ordinary differential equations is associated with the estimation of the norm of the transition operator from the current time level to a new one. The stability of operator-difference schemes for a first-order model operator-differential equation is discussed. Primary attention is given to the construction of additive schemes (splitting schemes) based on approximations of the transition operator. Specifically, classical factorized schemes, componentwise splitting schemes, and regularized operator-difference schemes are related to the use of a certain multiplicative transition operator. Additive averaged operator-difference schemes are based on an additive representation of the transition operator. The construction of second-order splitting schemes in time is discussed. Inhomogeneous additive operator-difference schemes are constructed in which various types of transition operators are used for individual splitting operators.     

Solutions of Semilinear Elliptic Equations in with Conical&Shaped Level Sets

This article deals with the questions of the existence, of the uniqueness and of the qulitative properties of solutions of semiliner elliptic equations. Three types of conical conditions at infinity are successively considered. this defines three frameworks: the weak framework, the strong are based on different kinds of sliding methods and, following the ideas of Berestycki, Nirenberg and Vega, on comparison principles in cones     

有限域上二次型序列与几何序列的互相关性

§1Introduction Letpbeaprimenumber,q=pm,andFqdenotethefinitefieldwithqelements.Fo anyn≥1,Trmnm(·)isthetracefunctionfromFqntoFq.LetαbeaprimitiveelementofFq and{α1,α2,...,αn}beabasisofFqnoverFq.Definition1.ForanonlinearfunctionffromFqtoFp,thesequenceS={Si}qn-1i=1withit termdefinedby Si=f(Trmnm(αi))(1iscalledageometricsequence.SuchageometricsequenceShasperioddividingqn-1.Geometricsequences includingm-sequence[1],GMWsequence[2,3],cascadedGMWsequence[4,5]andman others[6],mayhavelar…