A spectral collocation algorithm for two-point boundary value problem in fiber Raman amplifier equations

A novel algorithm implementing Chebyshev spectral collocation (pseudospectral) method in combination with Newton’s method is proposed for the nonlinear two-point boundary value problem (BVP) arising in solving propagation equations in fiber Raman amplifier. Moreover, an algorithm to train the known linear solution for use as a starting solution for the Newton iteration is proposed and successfully implemented. The exponential accuracy obtained by the proposed Chebyshev pseudospectral method is demonstrated on a case of the Raman propagation equations with strong nonlinearities. This is in contrast to algebraic accuracy obtained by typical solvers used in the literature. The resolving power and the efficiency of the underlying Chebyshev grid are demonstrated in comparison to a known BVP solver.     

Variational iteration method for solving integro-differential equations

In this Letter, the variational iteration method is applied to solve integro-differential equations. Some examples are given to illustrate the effectiveness of the method, the results show that the method provides a straightforward and powerful mathematical tool for solving various integro-differential equations.     

Approximate analytical solutions of generalized linear and nonlinear reaction-diffusion equations in an infinite domain

This Letter deals with the solution of unified fractional reaction-diffusion systems in an infinite domain. The results are obtained in compact and elegant forms in terms of generalized Mittag-Leffler functions, which are suitable for numerical computation.     

A linear nonconforming finite element method for Maxwell’s equations in two dimensions. Part I: Frequency domain

We suggest a linear nonconforming triangular element for Maxwell’s equations and test it in the context of the vector Helmholtz equation. The element uses discontinuous normal fields and tangential fields with continuity at the midpoint of the element sides, an approximation related to the Crouzeix–Raviart element for Stokes. The element is stabilized using the jump of the tangential fields, giving us a free parameter to decide. We give dispersion relations for different stability parameters and give some numerical examples, where the results converge quadratically with the mesh size for problems with smooth boundaries. The proposed element is free from spurious solutions and, for cavity eigenvalue problems, the eigenfrequencies that correspond to well-resolved eigenmodes are reproduced with the correct multiplicity.     

A hybrid multilevel method for high-order discretization of the Euler equations on unstructured meshes

Higher order discretization has not been widely successful in industrial applications to compressible flow simulation. Among several reasons for this, one may identify the lack of tailor-suited, best-practice relaxation techniques that compare favorably to highly tuned lower order methods, such as finite-volume schemes. In this paper we investigate solution algorithms in conjunction with high-order Spectral Difference discretization for the Euler equations, using such techniques as multigrid and matrix-free implicit relaxation methods. In particular we present a novel hybrid multilevel relaxation method that combines (optionally matrix-free) implicit relaxation techniques with explicit multistage smoothing using geometric multigrid. Furthermore, we discuss efficient implementation of these concepts using such tools as automatic differentiation.     

求解非线性差分方程孤立波解的直接代数法

推广了求解非线性差分方程孤立波解的直接代数法.用此方法研究了Hybrid晶格方程，借助于符号计算Maple,得到它的新孤波解.这种方法也可用于求解其他的差分方程. 关键词： 微分-差分方程 Hybrid晶格方程 行波解 孤     

Artificial boundary conditions for the numerical solution of the Euler equations by the discontinuous galerkin method

We present artificial boundary conditions for the numerical simulation of compressible flows using high-order accurate discretizations with the discontinuous Galerkin (DG) finite element method. The construction of the proposed boundary conditions is based on characteristic analysis and applied for boundaries with arbitrary shape and orientation. Numerical experiments demonstrate that the proposed boundary treatment enables to convect out of the computational domain complex flow features with little distortion. In addition, it is shown that small-amplitude acoustic disturbances could be convected out of the computational domain, with no significant deterioration of the overall accuracy of the method. Furthermore, it was found that application of the proposed boundary treatment for viscous flow over a cylinder yields superior performance compared to simple extrapolation methods.     

An adaptive SPH method for strong shocks

We propose an alternative SPH scheme to usual SPH Godunov-type methods for simulating supersonic compressible flows with sharp discontinuities. The method relies on an adaptive density kernel estimation (ADKE) algorithm, which allows the width of the kernel interpolant to vary locally in space and time so that the minimum necessary smoothing is applied in regions of low density. We have performed a von Neumann stability analysis of the SPH equations for an ideal gas and derived the corresponding dispersion relation in terms of the local width of the kernel. Solution of the dispersion relation in the short wavelength limit shows that stability is achieved for a wide range of the ADKE parameters. Application of the method to high Mach number shocks confirms the predictions of the linear analysis. Examples of the resolving power of the method are given for a set of difficult problems, involving the collision of two strong shocks, the strong shock-tube test, and the interaction of two blast waves.     

求一类非线性偏微分方程解析解的一种简洁方法

通过引入一个变换和选准试探函数，进行求偏导数运算,将非线性偏微分方程化为代数方程，然后用待定系数法确定相应的常数，最后得到其解析 解.不难看出，这种方法特别简洁. 关键词： 非线性偏微分方程 试探函数 待定系数法 解析解     

Hybridizable discontinuous Galerkin methods for the time-harmonic Maxwell’s equations

We present two hybridizable discontinuous Galerkin (HDG) methods for the numerical solution of the time-harmonic Maxwell’s equations. The first HDG method explicitly enforces the divergence-free condition and thus necessitates the introduction of a Lagrange multiplier. It produces a linear system for the degrees of freedom of the approximate traces of both the tangential component of the vector field and the Lagrange multiplier. The second HDG method does not explicitly enforce the divergence-free condition and thus results in a linear system for the degrees of freedom of the approximate trace of the tangential component of the vector field only. For both HDG methods, the approximate vector field converges with the optimal order of k + 1 in the L²-norm, when polynomials of degree k are used to represent all the approximate variables. We propose elementwise postprocessing to obtain a new H^curl-conforming approximate vector field which converges with order k + 1 in the H^curl-norm. We present extensive numerical examples to demonstrate and compare the performance of the HDG methods.     

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.     

Generalization of the FDTD algorithm for simulations of hydrodynamic nonlinear Drude model

In this paper we present a numerical method for solving a three-dimensional cold-plasma system that describes electron gas dynamics driven by an external electromagnetic wave excitation. The nonlinear Drude dispersion model is derived from the cold-plasma fluid equations and is coupled to the Maxwell’s field equations. The Finite-Difference Time-Domain (FDTD) method is applied for solving the Maxwell’s equations in conjunction with the time-split semi-implicit numerical method for the nonlinear dispersion and a physics based treatment of the discontinuity of the electric field component normal to the dielectric-metal interface. The application of the proposed algorithm is illustrated by modeling light pulse propagation and second-harmonic generation (SHG) in metallic metamaterials (MMs), showing good agreement between computed and published experimental results.     

Interior penalty DG methods for Maxwell’s equations in dispersive media

In this paper, we develop a fully-discrete interior penalty discontinuous Galerkin method for solving the time-dependent Maxwell’s equations in dispersive media. The model is described by a vector integral–differential equation. Our scheme is proved to be unconditionally stable and achieve optimal error estimates in both L² norm and energy norm. The scheme is implemented and numerical results supporting our analysis are presented.     

一类非线性微分方程的近似解析解

从理论上研究了一类广义扩散方程的求解问题. 利用相似变换和解析拆分技巧给出了求解该类非线性微分方程近似解的一种有效方法, 方程的解可以表示为一个收敛的幂级数. 近似解结果和数值结果非常符合，证明了所提出的方法的准确性和可靠性, 该方法可以用于解决其他科学和工程技术问题. 关键词： 广义扩散方程 非线性边界值问题 解析拆分 近似解析解     

Non-conformal domain decomposition method with second-order transmission conditions for time-harmonic electromagnetics

Non-overlapping domain decomposition (DD) methods provide efficient algorithms for solving time-harmonic Maxwell equations. It has been shown that the convergence of DD algorithms can be improved significantly by using high order transmission conditions. In this paper, we extend a newly developed second-order transmission condition (SOTC), which involves two second-order transverse derivatives, to facilitate fast convergence in the non-conformal DD algorithms. However, the non-conformal nature of the DD methods introduces an additional technical difficulty, which results in poor convergence in many real-life applications. To mitigate the difficulty, a corner-edge penalty method is proposed and implemented in conjunction with the SOTC to obtain truly robust solver performance. Numerical results verify the analysis and demonstrate the effectiveness of the proposed methods on a few model problems. Finally, drastically improved convergence, compared to the conventional Robin transmission condition, was observed for an electrically large problem of practical interest.     

用三角函数法获得非线性Boussinesq方程的广义孤子解

找到一个合适的代换——三角函数法，将非线性Boussinesq微分方程转换为非线性代数方程组.用吴消元法求解该非线性代数方程组，从而获得一般形式Boussinesq微分方程的广义孤子解. 关键词： Boussinesq方程 吴消元法 非线性代数方程组 孤子解     

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.     

Unified compact ADI methods for solving nonlinear viscous and nonviscous wave equations

In this paper, two unified alternating direction implicit (ADI) methods, based on the combination of fourth-order compact difference for the approximations of the second spatial derivatives with approximation factorization of difference operators, are presented for solving a two-dimensional (2D) and three-dimensional (3D) nonlinear viscous and nonviscous wave equations, respectively. By the discrete energy method, it is shown that their solutions converge to exact solutions with an order of two in time and four in space in L²- and H¹-norms. Finally, numerical findings testify the computational efficiency of the algorithms.     

A numerical method for one-dimensional nonlinear sine-Gordon equation using multiquadric quasi-interpolation

In this paper, we use a univariate multiquadric quasi-interpolation scheme to solve the one-dimensional nonlinear sine-Gordon equation that is related to many physical phenomena. We obtain a numerical scheme by using the derivative of the quasi-interpolation to approximate the spatial derivative and a difference scheme to approximate the temporal derivative. The advantage of the obtained scheme is that the algorithm is very simple so that it is very easy to implement. The results of numerical experiments are presented and compared with analytical solutions to confirm the good accuracy of the presented scheme.     

Application of He's variational iteration method to nonlinear heat transfer equations

Instead of finding a small parameter for solving nonlinear problems through perturbation method, a new analytical method called He's variational iteration method (VIM) is introduced to be applied to solve nonlinear heat transfer equations in this Letter. In this research, variational iteration method is used to solve an unsteady nonlinear convective-radiative equation and a nonlinear convective-radiative-conduction equation containing two small parameters of ε₁ and ε₂ and evaluate the efficiency of straight fins. VIM can apply to the nonlinear equations with boundary or initial conditions defined in different points just with developing the correction functional using the extra parameters such as Cn, as used in this Letter.