Stroboscopic averaging of highly oscillatory nonlinear wave equations

In this paper, we are concerned with stroboscopic averaging for highly oscillatory evolution equations posed in a Banach space. Using Taylor expansion, we construct a non‐oscillatory high‐order system whose solution remains exponentially close to the exact one over a long time. We then apply this result to the nonlinear wave equation in one dimension. We present the stroboscopic averaging method, which is a numerical method introduced by Chartier, Murua and Sanz‐Serna, and apply it to our problem. Finally, we conclude by presenting the qualitative and quantitative efficiency of this numerical method for some nonlinear wave problem. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical Solution for the Helmholtz Equation with Mixed Boundary Condition

We consider the numerical solution for the Helmholtz equation in R~2 with mixed boundary conditions.The solvability of this mixed boundary value problem is estab- lished by the boundary integral equation method.Based on the Green formula,we express the solution in terms of the boundary data.The key to the numerical real- ization of this method is the computation of weakly singular integrals.Numerical performances show the validity and feasibility of our method.The numerical schemes proposed in this paper have been applied in the realization of probe method for inverse scattering problems.     

Power series solutions for the KPP equation

In this paper we solve the KPP equation by a non numerical method. To this end we find power series solutions where the coefficients are computed recursively. We also prove convergence of the series and illustrate the method by few examples.      

Phase-Field Modeling and Numerical Simulation for Ice Melting

In this paper, we propose a mathematical model and present numerical simulations for ice melting phenomena. The model is based on the phase-field modeling for the crystal growth. To model ice melting, we ignore anisotropy in the crystal growth model and introduce a new melting term. The numerical solution algorithm is a hybrid method which uses both the analytic and numerical solutions. We perform various computational experiments. The computational results confirm the accuracy and efficiency of the proposed method for ice melting.     

利用随机模拟方法求解第二类积分方程

给出一种求解第二类Fredholm和Volterra积分方程的数值算法,算法在数值积分技术的基础上使用Monte Carlo随机模拟方法求积分方程的近似解.通过数值例子证明了该算法是有效的.     

Well-posedness and finite element approximation of time dependent generalized bioconvective flow

We consider the finite element approximation of a time dependent generalized bioconvective flow. The underlying system of partial differential equations consists of incompressible Navier–Stokes type convection equations coupled with an equation describing the transport of micro-organisms. The viscosity of the fluid is assumed to be a function of the concentration of the micro-organisms. We show the existence and uniqueness of the weak solution of the system in two dimensions and construct numerical approximations based on the finite element method, for which we obtain error estimates. In addition, we conduct several numerical experiments to demonstrate the accuracy of the numerical method and perform simulations of the bioconvection pattern formations based on realistic model parameters to demonstrate the validity of the proposed numerical algorithm.     

Implicit numerical approximation scheme for the fractional Fokker-Planck equation

In this paper, we consider an anomalous subdiffusion process, governed by fractional Fokker-Planck equation. An effective numerical method for approximating Fokker-Planck equation in a bounded domain is presented. The stability and convergence of the numerical method are analyzed. Some numerical examples are presented to show the application of the present technique. The numerical results exhibit the good performance of our theoretical analysis.     

Helmholtz方程Cauchy问题的间接积分方程方法

 本文处理二维和三维Helmholtz方程的边界数据复原问题.通过利用位势理论近似问题的解,导出了解决Cauchy问题的一种非迭代积分方程方法.为了处理形成问题的不适定性,采用了Tikhonov正则化结合Morozov偏差原理的方法,并且给出了算法的收敛性和误差估计,最后给出了二维和三维的数值算例.通过数值算例我们检验了源点和边界之间距离的关系,算法关于噪声、源点数目的数值收敛性,稳定性.     

Explicit Hybrid Numerical Method for the Allen-Cahn Type Equations on Curved Surfaces

We present a simple and fast explicit hybrid numerical scheme for the motion by mean curvature on curved surfaces in three-dimensional (3D) space. We numerically solve the Allen-Cahn (AC) and conservative Allen-Cahn (CAC) equations on a triangular surface mesh. We use the operator splitting method and an explicit hybrid numerical method. For the AC equation, we solve the diffusion term using a discrete Laplace-Beltrami operator on the triangular surface mesh and solve the reaction term using the closed-form solution, which is obtained using the separation of variables. Next, for the CAC equation, we additionally solve the time-space dependent Lagrange multiplier using an explicit scheme. Our numerical scheme is computationally fast and efficient because we use an explicit hybrid numerical scheme. We perform various numerical experiments to demonstrate the robustness and efficiency of the proposed scheme.     

A meshless based numerical technique for traveling solitary wave solution of Boussinesq equation

In this paper, we employ the boundary-only meshfree method to find out numerical solution of the classical Boussinesq equation in one dimension. The proposed method in the current paper is a combination of boundary knot method and meshless analog equation method. The boundary knot technique is an integration free, boundary-only, meshless method which is used to avoid the known disadvantages of the method of fundamental solution. Also, we use the meshless analog equation method to replace the nonlinear governing equation with an equivalent nonhomogeneous linear equation. A predictor-corrector scheme is proposed to solve the resulted differential equation of the collocation. The numerical results and conclusions are obtained for both the ‘good’ and the ‘bad’ Boussinesq equations.     

Analysis and numerical solution of Benjamin-Bona-Mahony equation with moving boundary

In this work we present the existence, the uniqueness and numerical solutions for a mathematical model associated with equations of Benjamin-Bona-Mahony type in a domain with moving boundary. We apply the Galerkin method, multiplier techniques, energy estimates and compactness results to obtain the existence and uniqueness. For numerical solutions, we shall employ the finite element method together with the Crank-Nicolson method. Some numerical experiments are presented to show the moving boundary for the problem.     

Mathematical Modeling and Simulation of Antibubble Dynamics

In this study, we propose a mathematical model and perform numerical simulations for the antibubble dynamics. An antibubble is a droplet of liquid surrounded by a thin film of a lighter liquid, which is also in a heavier surrounding fluid. The model is based on a phase-field method using a conservative Allen-Cahn equation with a space-time dependent Lagrange multiplier and a modified Navier-Stokes equation. In this model, the inner fluid, middle fluid and outer fluid locate in specific diffusive layer regions according to specific phase filed (order parameter) values. If we represent the antibubble with conventional binary or ternary phase-field models, then it is difficult to have stable thin film. However, the proposed approach can prevent nonphysical breakup of fluid film during the simulation. Various numerical tests are performed to verify the efficiency of the proposed model.     

On the numerical solution of Korteweg-de Vries equation by the iterative splitting method

In this paper, we apply the method of iterative operator splitting on the Korteweg-de Vries (KdV) equation. The method is based on first, splitting the complex problem into simpler sub-problems. Then each sub-equation is combined with iterative schemes and solved with suitable integrators. Von Neumann analysis is performed to achieve stability criteria for the proposed method applied to the KdV equation. The numerical results obtained by iterative splitting method for various initial conditions are compared with the exact solutions. It is seen that they are in a good agreement with each other.     

A fast collocation method for acoustic scattering in shallow oceans

In this paper, we investigate the numerical solution of the integral equation of the second kind reduced by acoustic scattering in shallow oceans with Dirichlet condition. Based on analyzing the singularity of the truncating kernel with a sum of infinite series, using our trigonometric interpolatory wavelets and collocation method, we obtain the numerical solution which possesses a fast convergence rate like o(2^−j). Moreover, the entries of the stiffness matrix can be obtained by FFT, which lead the computational complexity to decrease obviously.     

Numerical studies on nonlinear Schrödinger equations by spectral collocation method with preconditioning

In this study, we use the spectral collocation method with preconditioning to solve various nonlinear Schrödinger equations. To reduce round-off error in spectral collocation method we use preconditioning. We study the numerical accuracy of the method. The numerical results obtained by this way have been compared with the exact solution to show the efficiency of the method.     

二维刚性边值问题的小波-吉尔解法

本文利用微分方程数值解的离散小波表示，讨论了此类方程在满足一定初始条件和边值条件下，在一个方向上利用小波伽辽金法，另一方向上利用吉尔方法进行求解，提出了一种解二维刚性初，边值问题的小波数值算法，计算结果表明，利用该方法所求得的数值解精度高，而且由小波特有的性质，它特别适用于求解带有奇异摄动的刚性问题。     

On approximation in the fictitious domain method for a bidimensional transport equation

In this article a numerical method for solving a two‐dimensional transport equation in the stationary case is presented. Using the techniques of the variational calculus, we find the approximate solution for a homogeneous boundary‐value problem that corresponds to a square domain D₂. Then, using the method of the fictitious domain, we extend our algorithm to a boundary value problem for a set D that has an arbitrary shape. In this approach, the initial computation domain D (called physical domain) is immersed in a square domain D₂. We prove that the solution obtained by this method is a good approximation of the exact solution. The theoretical results are verified with the help of a numerical example. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Local Discontinuous Galerkin Methods for Three Classes of Nonlinear Wave Equations

In this paper, we further develop the local discontinuous Galerkin method to solve three classes of nonlinear wave equations formulated by the general KdV-Burgers type equations, the general fifth-order KdV type equations and the fully nonlinear K(n, n, n) equations, and prove their stability for these general classes of nonlinear equations. The schemes we present extend the previous work of Yan and Shu [30, 31] and of Levy, Shu and Yan [24] on local discontinuous Galerkin method solving partial differential equations with higher spatial derivatives. Numerical examples for nonlinear problems are shown to illustrate the accuracy and capability of the methods. The numerical experiments include stationary solitons, soliton interactions and oscillatory solitary wave solutions.The numerical experiments also include the compacton solutions of a generalized fifthorder KdV equation in which the highest order derivative term is nonlinear and the fully nonlinear K (n, n, n) equations.     

ZK-BBM方程的多辛Preissmann格式

本文研究一类非线性ZK-BBM方程的初值问题.利用Hamilton系统的多辛Preissmann方法,获得ZK-BBM方程初值问题的数值结果,数值结果表明该多辛离散格式具有较好的长时间数值稳定性.     

ID-WAVELETS METHOD FOR HAMMERSTEIN INTEGRAL EQUATIONS

1.IntroductionInthispaperwewillconsiderthenumericalsolutionsofthenon--linearintegralequationsofHammersteintype:wheref,kandgaregivenfunctionandyistheunknown.TherehasbeenmuchinterestinthisproblemsinceHammersteinintegralequations,whichcamefromtheelectromagneticfluiddynamics,yieldsstrongphysicalbackground.Moreover,theFredholmintegralequationsofsecondkindarethespecialcaseoftheHammersteinintegralequations.In[6,p.700]thestandardcollocationmethodisappliedtoobtaintheapproximationsolutionofEq.(1).Int…