A method for solving nonlinear equations

A direct method for solving constant coefficient nonlinear evolutional problems of mathematical physics extending the Hirota method to the case of degeneracy of k-parameter is described. An application of this method to the problem of finding the solutions describing coupled states in φ⁴-theory with damping is considered.     

Transformation methods for solving nonlinear field equations

The collection of extended canonical transformations of first-order contact manifolds is studied. This collection is shown to form a group under target-source composition and to contain the group of all first prolongations of point transformation of the underlying graph space and all isogroups of completely integrable horizontal ideals. Extended canonical transformations are compared and contrasted with Bäcklund transformations. These results are used to construct an extended Hamilton-Jacobi method for systems of nonlinear PDE. The collection of all extended canonical transformations is also shown to contain infinitely many one-parameter families of transformations, but there is no Lie group structure that contains these one-parameter families, in general. Conditions are obtained under which a one-parameter family of extended canonical transformations will map a solution of the fundamental ideal that characterizes a given system of PDE into a one-parameter family of solutions. These results are applied to the -Gordon equation _x1 = () and to the Navier-Stokes equations.     

An efficient method for solving highly anisotropic elliptic equations

Solving elliptic PDEs in more than one dimension can be a computationally expensive task. For some applications characterized by a high degree of anisotropy in the coefficients of the elliptic operator, such that the term with the highest derivative in one direction is much larger than the terms in the remaining directions, the discretized elliptic operator often has a very large condition number – taking the solution even further out of reach using traditional methods. This paper will demonstrate a solution method for such ill-behaved problems. The high condition number of the D-dimensional discretized elliptic operator will be exploited to split the problem into a series of well-behaved one and (D − 1)-dimensional elliptic problems. This solution technique can be used alone on sufficiently coarse grids, or in conjunction with standard iterative methods, such as Conjugate Gradient, to substantially reduce the number of iterations needed to solve the problem to a specified accuracy. The solution is formulated analytically for a generic anisotropic problem using arbitrary coordinates, hopefully bringing this method into the scope of a wide variety of applications.     

A deep learning method for solving third-order nonlinear evolution equations

It has still been difficult to solve nonlinear evolution equations analytically. In this paper, we present a deep learning method for recovering the intrinsic nonlinear dynamics from spatiotemporal data directly. Specifically, the model uses a deep neural network constrained with given governing equations to try to learn all optimal parameters. In particular, numerical experiments on several third-order nonlinear evolution equations, including the Korteweg–de Vries (KdV) equation, modified KdV equation, KdV–Burgers equation and Sharma–Tasso–Olver equation, demonstrate that the presented method is able to uncover the solitons and their interaction behaviors fairly well.     

A deep learning method for solving high-order nonlinear soliton equations

We propose an effective scheme of the deep learning method for high-order nonlinear soliton equations and explore the influence of activation functions on the calculation results for higher-order nonlinear soliton equations. The physics-informed neural networks approximate the solution of the equation under the conditions of differential operator, initial condition and boundary condition. We apply this method to high-order nonlinear soliton equations, and verify its efficiency by solving the fourth-order Boussinesq equation and the fifth-order Korteweg–de Vries equation. The results show that the deep learning method can be used to solve high-order nonlinear soliton equations and reveal the interaction between solitons.     

Boundary regularity for some nonlinear elliptic degenerate equations

We consider the nonlinear elliptic degenerate equation (1) $$ - x^2 \left( {\frac{{\partial ^2 u}}{{\partial x^2 }} + \frac{{\partial ^2 u}}{{\partial y^2 }}} \right) + 2u = f(u)in\Omega _a ,$$ where $$\Omega _a = \left\{ {(x,y) \in \mathbb{R}^2 ,0< x< a,\left| y \right|< a} \right\}$$ for some constanta>0 andf is aC ^∞ functions on ? such thatf(0)=f′(0)=0. Our main result asserts that: ifu∈C \((\bar \Omega _a )\) satisfies (2) $$u(0,y) = 0for\left| y \right|< a,$$ thenx ^?2 u(x,y)∈C ^∞ \(\left( {\bar \Omega _{a/2} } \right)\) and in particularu∈C ^∞ \(\left( {\bar \Omega _{a/2} } \right)\) .     

On spectral methods for solving nonlinear acoustics equations

Two very efficient methods for obtaining approximate solutions to nonlinear acoustics equations are discussed. I proposed these methods earlier, but they are still little known. The first method is based on expanding an unknown function into a Taylor series with respect to the coordinate (evolution variable) and on approximate summation of the terms of this series in all orders up to the infinite order. This series can be summed completely only in particular cases, e.g., for a simple wave. It has been noted that the partial summation technique is implemented more easily if all the terms of the series are represented as corresponding topological diagrams. The second method is based on introducing a “nonlinear” phase delay (proportional to the wave amplitude) for the temporal variable in linear solutions of the problem. The application technique of these methods is illustrated by obtaining approximate solutions of the Burgers equation.     

非线性离散微分方程的双曲函数法求解

本文推广了双曲函数方法用于求解非线性离散系统。求解离散的(2+1)维Toda系统和离散的mKdV系统，成功地得到了离散钟型孤立子、离散冲击波型孤立子及一些新的精确行波解。     

Improved algorithm for solving nonlinear parabolized stability equations

Due to its high computational efficiency and ability to consider nonparallel and nonlinear effects, nonlinear parabolized stability equations(NPSE) approach has been widely used to study the stability and transition mechanisms. However,it often diverges in hypersonic boundary layers when the amplitude of disturbance reaches a certain level. In this study, an improved algorithm for solving NPSE is developed. In this algorithm, the mean flow distortion is included into the linear operator instead of into the nonlinear forcing terms in NPSE. An under-relaxation factor for computing the nonlinear terms is introduced during the iteration process to guarantee the robustness of the algorithm. Two case studies, the nonlinear development of stationary crossflow vortices and the fundamental resonance of the second mode disturbance in hypersonic boundary layers, are presented to validate the proposed algorithm for NPSE. Results from direct numerical simulation(DNS) are regarded as the baseline for comparison. Good agreement can be found between the proposed algorithm and DNS, which indicates the great potential of the proposed method on studying the crossflow and streamwise instability in hypersonic boundary layers.     

A nonconforming combination for solving elliptic problems with interfaces

Nonconforming combinations are provided for solving interface problems of elliptic equations. In these approaches, the Ritz-Galerkin method with particular solutions is used for the part of a solution domain where there are interface singular points; and the conventional finite element method is used for the rest of the solution domain. In addition, admissible functions chosen are constrained to be continuous only at the element nodes on the common boundary of the subdomains. Error bounds are derived in the Sobolev norms, and numerical experiments are given for solving a model interface problem of the equation, −Δu + U = 0. Moreover, a significant coupling relation, L + 1 = O(|ln h|), is found for interface problems by using the nonconforming combinations, where (L + 1) is the total number of particular solutions used in the Ritz-Galerkin method, and h is the maximal boundary length of triangular elements in the finite element method.     

谐振管内非线性驻波的有限体积数值算法

为了扩展谐振管内非线性驻波在工程中的应用, 以及克服现有数值计算方法仅局限于求解直圆柱形和指数形谐振管内非线性驻波的问题. 根据变截面的非稳态可压缩热黏性流体Navier-Stokes方程和空间守恒方程, 并基于求解压力速度耦合方程的半隐式算法和交错网格技术, 构建一种能够计算任意形状轴对称谐振管受活塞驱动时内部非线性驻波的有限体积算法. 分别对圆柱形、指数形和圆锥形谐振管内的非线性驻波进行仿真计算. 通过与现有试验结果以及数值仿真结果的对比, 验证了该方法的正确性.并获得除驻波声压之外的另外一些新的物理结果, 包括速度、密度、温度的瞬时变化.在直圆柱形谐振管内产生冲击声压波, 速度波形中出现钉状结构.而在指数形和圆锥形谐振管内产生高声压幅值的驻波, 没有出现冲击波, 速度波形中均未发现钉状结构. 计算结果表明谐振管内非线性驻波的物理属性与谐振管形状之间有密切关系.     

Covolume-upwind finite volume approximations for linear elliptic partial differential equations

In this paper, we study covolume-upwind finite volume methods on rectangular meshes for solving linear elliptic partial differential equations with mixed boundary conditions. To avoid non-physical numerical oscillations for convection-dominated problems, nonstandard control volumes (covolumes) are generated based on local Peclet’s numbers and the upwind principle for finite volume approximations. Two types of discretization schemes with mass lumping are developed with use of bilinear or biquadratic basis functions as the trial space respectively. Some stability analyses of the schemes are presented for the model problem with constant coefficients. Various examples are also carried out to numerically demonstrate stability and optimal convergence of the proposed methods.     

Spectral splitting method for nonlinear Schrödinger equations with singular potential

We consider the time-dependent one-dimensional nonlinear Schrödinger equation with pointwise singular potential. By means of spectral splitting methods we prove that the evolution operator is approximated by the Lie evolution operator, where the kernel of the Lie evolution operator is explicitly written. This result yields a numerical procedure which is much less computationally expensive than multi-grid methods previously used. Furthermore, we apply the Lie approximation in order to make some numerical experiments concerning the splitting of a soliton, interaction among solitons and blow-up phenomenon.     

Chebyshev pseudospectral solution of second-order elliptic equations with finite element preconditioning

A new Chebyshev pseudospectral algorithm for second-order elliptic equations using finite element preconditioning is proposed and tested on various problems. Bilinear and biquadratic Lagrange elements are considered as well as bicubic Hermite elements. The numerical results show that bilinear elements produce spectral accuracy with the minimum computational work. L-shaped regions are treated by a subdomain approach.     

Employment of Jacobian elliptic functions for solving problems in nonlinear dynamics of microtubules

We show how Jacobian elliptic functions （JEFs） can be used to solve ordinary differential equations （ODEs） describing the nonlinear dynamics of microtubules （MTs）. We demonstrate that only one of the JEFs can be used while the remaining two do not represent the solutions of the crucial differential equation. We show that a kinkbtype soliton moves along MTs. Besides this solution, we also discuss a few more solutions that may or may not have physical meanings. Finally, we show what kind of ODE can be solved by using JEFs.     

非线性发展方程的丰富的Jacobi椭圆函数解

通过把十二个Jacobi椭圆函数分类成四组，提出了新的广泛的Jacobi椭圆函数展开法，利用这一方法求得了非线性发展方程的丰富的Jacobi椭圆函数双周期解.当模数m→0或1时，这些解退化为相应的三角函数解或孤立波解和冲击波解. 关键词： 非线性发展方程 Jacobi椭圆函数 双周期解 行波解     

非线性波动方程的Jacobi椭圆函数包络周期解

应用Jacobi椭圆函数展开法求得了一类非线性波方程的包络周期解,而且用这种方法得到的周期解在一定条件下可以退化为包络冲击波解或包络孤立波解 关键词： Jacobi椭圆函数 非线性方程 包络周期解 包络孤立波解     

The modified simple equation method for solving some fractional-order nonlinear equations

Nonlinear fractional differential equations are encountered in various fields of mathematics, physics, chemistry, biology, engineering and in numerous other applications. Exact solutions of these equations play a crucial role in the proper understanding of the qualitative features of many phenomena and processes in various areas of natural science. Thus, many effective and powerful methods have been established and improved. In this study, we establish exact solutions of the time fractional biological population model equation and nonlinear fractional Klein–Gordon equation by using the modified simple equation method.