Travelling Waves for a Density Dependent Diffusion Nagumo Equation over the Real Line

We consider the density dependent diffusion Nagumo equation, where the diffusion coefficient is a simple power function. This equation is used in modelling electrical pulse propagation in nerve axons and in population genetics (amongst other areas). In the present paper, the δ-expansion method is applied to a travelling wave reduction of the problem, so that we may obtain globally valid perturbation solutions (in the sense that the perturbation solutions are valid over the entire infinite domain, not just locally; hence the results are a generalization of the local solutions considered recently in the literature). The resulting boundary value problem is solved on the real line subject to conditions at z → ±∞. Whenever a perturbative method is applied, it is important to discuss the accuracy and convergence properties of the resulting perturbation expansions. We compare our results with those of two different numerical methods (designed for initial and boundary value problems, respectively) and deduce that the perturbation expansions agree with the numerical results after a reasonable number of iterations. Finally, we are able to discuss the influence of the wave speed c and the asymptotic concentration value α on the obtained solutions. Upon recasting the density dependent diffusion Nagumo equation as a two-dimensional dynamical system, we are also able to discuss the influence of the nonlinear density dependence (which is governed by a power-law parameter m) on oscillations of the travelling wave solutions.     

Two-component Fermi systems: II. Superfluid coupled cluster theory

In this second paper of a series the coupled cluster method (CCM) or exp(S) formalism is applied to two-component Fermi superfluids using a Bardeen-Cooper-Schrieffer (BCS) ground state as a zeroth-order approximation. We concentrate on developing the formalism necessary for carrying out eventual numerical calculations on realistic superconducting systems. We do this by generalising the one-component formalism in an appropriate manner and by using the results in the first paper of this series, where we studied two-component Fermi fluids. We stress the previous successes of the CCM, both from the point of view of analytic and numerical results, and we further indicate its potential for studying superconductivity. We restrict ourselves here to a so-called ring plus single particle energy (RING+SPE) approximation for general potentials and show how it can be formulated as a set of four coupled, bilinear integral equations for the cluster-integrated amplitudes. These latter amplitudes are themselves derived from the four-point functions of the system which provide a measure of the two-particle/two-hole component in the true ground-state wavefunction with respect to the BCS model state. We indicate how to obtain possible analytic solutions.     

Numerical and analytical studies of nonequilibrium fluctuation-induced transport processes

We present a numerical simulation algorithm that is well suited for the study of noise-induced transport processes. The algorithm has two advantages over standard techniques: (1) it preserves the property of detailed balance for systems in equilibrium and (2) it provides an efficient method for calculating nonequilibrium currents. Numerical results are compared with exact solutions from two different types of correlation ratchets, and are used to verify the results of perturbation calculations done on a three-state ratchet.     

Direct Perturbation Method for Derivative Nonlinear Schrödinger Equation

We extend Lou's direct perturbation method for solving the nonlinear Schrödinger equation to the case of the derivative nonlinear Schrödinger equation (DNLSE). By applying this method, different types of perturbation solutions are obtained. Based on these approximate solutions, the analytical forms of soliton parameters, such as the velocity, the width and the initial position, are carried out and the effects of perturbation on solitons are analyzed at the same time. A numerical simulation of perturbed DNLSE finally verifies the results of the perturbation method.     

Approximate Bound States Solution of the Hellmann Potential

The Hellmann potential, which is a superposition of an attractive Coulomb potential -a/r and a Yukawa potential b e^-δr/r, is often used to compute bound-state normalizations and energy levels of neutral atoms. By using the generalized parametric Nikiforov-Uvarov (NU) method, we have obtained the approximate analytical solutions of the radial Schrödinger equation (SE) for the Hellmann potential. The energy eigenvalues and corresponding eigenfunctions are calculated in closed forms. Some numerical results are presented, which show good agreement with a numerical amplitude phase method and also those previously obtained by other methods. As a particular case, we find the energy levels of the pure Coulomb potential.     

Effect of magnetic field on the stimulated photon echo in ytterbium vapors

Polarization of the stimulated photon echo (SPE) at the 0 ? 1 transition in ytterbium vapors in the presence of a longitudinal magnetic field of 0–40 G is experimentally and theoretically studied. The SPE is generated using three light pulses with identical linear polarizations, so that the SPE polarization is the same at zero magnetic field. In the presence of a weak magnetic field, the SPE polarization vector rotates around the magnetic field vector and the depolarization of the SPE signal takes place. Each of the SPE polarization components exhibits biharmonic oscillations depending on the magnetic field. In the presence of a strong magnetic field, these oscillations vanish and the SPE becomes depolarized. The experimental data are in qualitative agreement with the results of the numerical calculations performed with the method of the evolution operator for the finite-duration excitation pulses. The application of the results for the processing of optical data is discussed.     

Padé approximations of quantized-vortex solutions of the Gross–Pitaevskii equation

Quantized vortices are important topological excitations in Bose–Einstein condensates. The Gross–Pitaevskii equation is a widely accepted theoretical tool. High accuracy quantized-vortex solutions are desirable in many numerical and analytical studies. We successfully derive the Padéapproximate solutions for quantized vortices with winding numbers ω = 1, 2, 3, 4, 5, 6 in the context of the Gross–Pitaevskii equation for a uniform condensate. Compared with the numerical solutions, we find that(1) they approximate the entire solutions quite well from the core to infinity;(2) higher-order Padé approximate solutions have higher accuracy;(3) Padé approximate solutions for larger winding numbers have lower accuracy. The healing lengths of the quantized vortices are calculated and found to increase almost linearly with the winding number. Based on experiments performed with ~(87)Rb cold atoms, the healing lengths of quantized vortices and the number of particles within the healing lengths are calculated, and they may be checked by experiment. Our results show that the Gross–Pitaevskii equation is capable of describing the structure of quantized vortices and physics at length scales smaller than the healing length.     

Numerical Solutions of Variable Coefficient Higher-Order Partial Differential Equations Arising in Beam Models

In this work, an efficient and robust numerical scheme is proposed to solve the variable coefficients’ fourth-order partial differential equations (FOPDEs) that arise in Euler–Bernoulli beam models. When partial differential equations (PDEs) are of higher order and invoke variable coefficients, then the numerical solution is quite a tedious and challenging problem, which is our main concern in this paper. The current scheme is hybrid in nature in which the second-order finite difference is used for temporal discretization, while spatial derivatives and solutions are approximated via the Haar wavelet. Next, the integration and Haar matrices are used to convert partial differential equations (PDEs) to the system of linear equations, which can be handled easily. Besides this, we derive the theoretical result for stability via the Lax–Richtmyer criterion and verify it computationally. Moreover, we address the computational convergence rate, which is near order two. Several test problems are given to measure the accuracy of the suggested scheme. Computations validate that the present scheme works well for such problems. The calculated results are also compared with the earlier work and the exact solutions. The comparison shows that the outcomes are in good agreement with both the exact solutions and the available results in the literature.     

Numerical soliton-like solutions of the potential Kadomtsev–Petviashvili equation by the decomposition method

In this Letter we present an Adomian's decomposition method (shortly ADM) for obtaining the numerical soliton-like solutions of the potential Kadomtsev–Petviashvili (shortly PKP) equation. We will prove the convergence of the ADM. We obtain the exact and numerical solitary-wave solutions of the PKP equation for certain initial conditions. Then ADM yields the analytic approximate solution with fast convergence rate and high accuracy through previous works. The numerical solutions are compared with the known analytical solutions.     

A new approach for modelling the damped Helmholtz oscillator:applications to plasma physics and electronic circuits

In this paper,a new approach is devoted to find novel analytical and approximate solutions to the damped quadratic nonlinear Helmholtz equation(HE)in terms of the Weiersrtrass elliptic function.The exact solution for undamped HE(integrable case)and approximate/semi-analytical solution to the damped HE(non-integrable case)are given for any arbitrary initial conditions.As a special case,the necessary and sufficient condition for the integrability of the damped HE using an elementary approach is reported.In general,a new ansatz is suggested to find a semi-analytical solution to the non-integrable case in the form of Weierstrass elliptic function.In addition,the relation between the Weierstrass and Jacobian elliptic functions solutions to the integrable case will be derived in details.Also,we will make a comparison between the semi-analytical solution and the approximate numerical solutions via using Runge-Kutta fourth-order method,finite difference method,and homotopy perturbation method for the first-two approximations.Furthermore,the maximum distance errors between the approximate/semi-analytical solution and the approximate numerical solutions will be estimated.As real applications,the obtained solutions will be devoted to describe the characteristics behavior of the oscillations in RLC series circuits and in various plasma models such as electronegative complex plasma model.     

A Hybrid Data-Differencing and Compression Algorithm for the Automotive Industry

We propose an innovative delta-differencing algorithm that combines software-updating methods with LZ77 data compression. This software-updating method relates to server-side software that creates binary delta files and to client-side software that performs software-update installations. The proposed algorithm creates binary-differencing streams already compressed from an initial phase. We present a software-updating method suitable for OTA software updates and the method’s basic strategies to achieve a better performance in terms of speed, compression ratio or a combination of both. A comparison with publicly available solutions is provided. Our test results show our method, Keops, can outperform an LZMA (Lempel–Ziv–Markov chain-algorithm) based binary differencing solution in terms of compression ratio in two cases by more than 3% while being two to five times faster in decompression. We also prove experimentally that the difference between Keops and other competing delta-creator software increases when larger history buffers are used. In one case, we achieve a three times better performance for a delta rate compared to other competing delta rates.     

Bound states resulting from interaction of the non-relativistic particles with the multiparameter potential

In this study, we present the analytical solutions of bound states for the Schrodinger equation with the multiparameter potential containing the different types of physical potentials via the asymptotic iteration method by applying the Pekeristype approximation to the centrifugal potential. For any n and l(states) quantum numbers, we derive the relation that gives the energy eigenvalues for the bound states numerically and the corresponding normalized eigenfunctions. We also plot some graphics in order to investigate effects of the multiparameter potential parameters on the energy eigenvalues.Furthermore, we compare our results with the ones obtained in previous works and it is seen that our numerical results are in good agreement with the literature.     

Numerical simulation of the regularized long wave equation by He's homotopy perturbation method

In this Letter, we present the homotopy perturbation method (shortly HPM) for obtaining the numerical solution of the RLW equation. We obtain the exact and numerical solutions of the Regularized Long Wave (RLW) equation for certain initial condition. The initial approximation can be freely chosen with possible unknown constants which can be determined by imposing the boundary and initial conditions. Comparison of the results with those of other methods have led us to significant consequences. The numerical solutions are compared with the known analytical solutions.     

Propagation of shape-preserving optical pulses in inhomogeneously broadened multi-level systems

We study stable propagation of multiple shape-preserving optical pulses in an inhomogeneously broadened multi-level atomic medium. By analytically solving the Maxwell-Schr?dinger equations governing the evolution of N coupled optical fields and atomic amplitudes we show that N pulsed optical waves coupling to (N+1)-levels can be automatically matched with the same soliton waveform and identical yet very slow propagation velocity. Several sets of coupled soliton solutions for two different (N+1)-level models are given and their stability is studied by using a numerical simulation.     

Application of homotopy perturbation method to the RLW and generalized modified Boussinesq equations

In this Letter, He's homotopy perturbation method (HPM) is implemented for finding the solitary-wave solutions of the regularized long-wave (RLW) and generalized modified Boussinesq (GMB) equations. We obtain numerical solutions of these equations for the initial conditions. We will show that the convergence of the HPM is faster than those obtained by the Adomian decomposition method (ADM). The obtained solutions, in comparison with the exact solutions admit a remarkable accuracy. A clear conclusion can be drawn from the numerical results that the HPM provides highly accurate numerical solutions for nonlinear differential equations.     

Effect of Bacterial Memory Dependent Growth by Using Fractional Derivatives Reaction-Diffusion Chemotactic Model

In this paper, numerical solutions of a reaction-diffusion chemotactic model of fractional orders for bacterial growth will be present. A new solution is constructed in power series. The fractional derivatives are described in the Caputo sense. We compare the experimental result obtained with those obtained by simulation of the chemotactic model without fractional derivatives. The results show that the solution continuously depends on the time-fractional derivative. The resulting solutions spread faster than the classical solutions and may exhibit asymmetry, depending on the fractional derivative used. We present results of numerical simulations to illustrate the method, and investigate properties of numerical solutions. The Adomian’s decomposition method (ADM) is used to find the approximate solution of fractional ‘reaction-diffusion chemotactic model. Numerical results show that the approach is easy to implement and accurate when applied to partial differential equations of fractional order.     

Explicit multi-symplectic extended leap-frog methods for Hamiltonian wave equations

In this paper, we study the integration of Hamiltonian wave equations whose solutions have oscillatory behaviors in time and/or space. We are mainly concerned with the research for multi-symplectic extended Runge–Kutta–Nyström (ERKN) discretizations and the corresponding discrete conservation laws. We first show that the discretizations to the Hamiltonian wave equations using two symplectic ERKN methods in space and time respectively lead to an explicit multi-symplectic integrator (Eleap-frogI). Then we derive another multi-symplectic discretization using a symplectic ERKN method in time and a symplectic partitioned Runge–Kutta method, which is equivalent to the well-known Störmer–Verlet method in space (Eleap-frogII). These two new multi-symplectic schemes are extensions of the leap-frog method. The numerical stability and dispersive properties of the new schemes are analyzed. Numerical experiments with comparisons are presented, where the two new explicit multi-symplectic methods and the leap-frog method are applied to the linear wave equation and the Sine–Gordon equation. The numerical results confirm the superior performance and some significant advantages of our new integrators in the sense of structure preservation.     

串级结构Z扫描特性的分析

对非线性介质串级结构的Z扫描特性进行了分析，分别给出了解析解和数值解. 解析分析主要是基于高斯分解方法（GDM）和分布透镜方法的组合；数值分析采用了Crank-Nicholson有限差分方法和快速汉克尔变换方法(FHT). 同时，也给出了这两种数值分析的使用条件和方法. 将解析结果、数值结果和实验结果进行了比较，结果显示它们具有很好的一致性. 关键词： 串级结构 Z扫描 有限差分方法 快速汉克尔变换     

水平变化波导中的全局矩阵耦合简正波方法

提出了一种新的水平变化波导中声场的耦合简正波求解方法,该方法能够处理二维点源和线源问题,提供声场的双向解。该方法利用全局矩阵(DGM)一次性求解耦合模式的系数,消除了传播矩阵递推求解中存在的误差累积问题;此外,改善了现有模型中对距离函数的归一化方法,从而避免了泄露模式指数增长导致的数值溢出问题。本文还给出了绝对软海底理想波导中耦合矩阵的闭合表达式,并分析了单个阶梯下简正波耦合现象。此外,本文还计算了理想楔形波导中的声传播问题(ASA标准问题),并与解析解及COUPLE07计算结果进行了比较,结果表明该方法是一种稳定、精确的水平变化波导中的声场计算方法。      

The shallow water equations in Lagrangian coordinates

Recent advances in the collection of Lagrangian data from the ocean and results about the well-posedness of the primitive equations have led to a renewed interest in solving flow equations in Lagrangian coordinates. We do not take the view that solving in Lagrangian coordinates equates to solving on a moving grid that can become twisted or distorted. Rather, the grid in Lagrangian coordinates represents the initial position of particles, and it does not change with time. We apply numerical methods traditionally used to solve differential equations in Eulerian coordinates, to solve the shallow water equations in Lagrangian coordinates. The difficulty with solving in Lagrangian coordinates is that the transformation from Eulerian coordinates results in solving a highly nonlinear partial differential equation. The non-linearity is mainly due to the Jacobian of the coordinate transformation, which is a precise record of how the particles are rotated and stretched. The inverse Jacobian must be calculated, thus Lagrangian coordinates cannot be used in instances where the Jacobian vanishes. For linear (spatial) flows we give an explicit formula for the Jacobian and describe the two situations where the Lagrangian shallow water equations cannot be used because either the Jacobian vanishes or the shallow water assumption is violated. We also prove that linear (in space) steady state solutions of the Lagrangian shallow water equations have Jacobian equal to one. In the situations where the shallow water equations can be solved in Lagrangian coordinates, accurate numerical solutions are found with finite differences, the Chebyshev pseudospectral method, and the fourth order Runge–Kutta method. The numerical results shown here emphasize the need for high order temporal approximations for long time integrations.