Adomian Decomposition Method and Pade Approximants for Nonlinear Differential-Difference Equations

Combining Adomian decomposition method （ADM） with Pade approximants, we solve two differentiaidifference equations （DDEs）： the relativistic Toda lattice equation and the modified Volterra lattice equation. With the help of symbolic computation Maple, the results obtained by ADM-Pade technique are compared with those obtained by using ADM alone. The numerical results demonstrate that ADM-Pade technique give the approximate solution with faster convergence rate and higher accuracy and relative in larger domain of convergence than using ADM.     

Identifying the temperature distribution in a parabolice quation with overspecified data using a multiquadric quasi-interpolation method

In this paper, we use a kind of univariate multiquadric quasi-interpolation to solve a parabolic equation with overspecified data, which has arisen in many physical phenomena. We obtain the numerical scheme by using the derivative of the quasi-interpolation to approximate the spatial derivative of the dependent variable and a simple forward difference to approximate the temporal derivative of the dependent variable. The advantage of the presented scheme is that the algorithm is very simple so it is very easy to implement. The results of the numerical experiment are presented and are compared with the exact solution to confirm the good accuracy of the presented scheme.     

Estimation of spatially distributed processes using mobile sensor networks with missing measurements

This paper investigates the estimation problem for a spatially distributed process described by a partial differential equation with missing measurements.The randomly missing measurements are introduced in order to better reflect the reality in the sensor network.To improve the estimation performance for the spatially distributed process,a network of sensors which are allowed to move within the spatial domain is used.We aim to design an estimator which is used to approximate the distributed process and the mobile trajectories for sensors such that,for all possible missing measurements,the estimation error system is globally asymptotically stable in the mean square sense.By constructing Lyapunov functionals and using inequality analysis,the guidance scheme of every sensor and the convergence of the estimation error system are obtained.Finally,a numerical example is given to verify the effectiveness of the proposed estimator utilizing the proposed guidance scheme for sensors.     

Accurate studies of dissociation energies and vibrational energies on alkali metals

This paper studies full vibrational spectra {Ev} and molecular dissociation energies De by using conventional least-squares （LS） fitting and an algebraic method （AM） proposed recently for 10 diatomic electronic states of ^7Li2, Na2, NaK and NaLi molecules based on some known experimental vibrational energies in a subset [Ev^expt] respectively. Studies show that： （1） although both the full AM spectrum {Ev^AM} and the LS spectrum {Ev^LS} can reproduce the known experimental energies in [Ev^expt], the {EAM} is superior to the {Ev^LS} in that the high-lying AM vibrational energies which may not be available experimentally have better or much better accuracy than those LS counterparts in {Ev^LS}, and so is the AM dissociation energy De^AM; （2） the main source of the errors in the data obtained by using the LS fitting is that the fitting which is just a pure mathematical process does not use any physical criteria that must be satisfied by the full vibrational spectrum, while the AM method does. This study suggests that when fitting or solving a physical equation using a set of source data, it is important not only to apply a proper mathematical tool, but also to use correct physical criteria which measure the physical properties of the data, kick out those data having bigger errors, and impose conditional convergence on the numerical process.     

Non—equilibrium extrapolation method for velocity and pressure boundary conditions in the lattice Boltzmann method

In this paper, we propose a new approach to implementing boundary conditions in the lattice Boltzmann method (LBM). The basic idea is to decompose the distribution function at the boundary node into its equilibrium and non-equilibrium parts, and then to approximate the non-equilibrium part with a first-order extrapolation of the non-equilibrium part of the distribution at the neighbouring fluid node. Schemes for velocity and pressure boundary conditions are constructed based on this method. The resulting schemes are of second-order accuracy. Numerical tests show that the numerical solutions of the LBM together with the present boundary schemes are in excellent agreement with the analytical solutions. Second-order convergence is also verified from the results. It is also found that the numerical stability of the present schemes is much better than that of the original extrapolation schemes proposed by Chen et al. (1996 Phys. Fluids 8 2527).     

Magnetic behaviors of cerium oxide-based thin films deposited using electrochemical method

Zn and Co multi-doped CeO2 thin films have been prepared using an anodic electrochemical method. The structures and magnetic behaviors are characterized by several techniques, in which the oxygen states in the lattice and the absorptive oxygen bonds at the surface are carefully examined. The absorptive oxygen bond is about 50% of the total oxygen bond by using a semi-quantitative method. The value of actual stoichiometry δ is close to 2. The experimental results indicate that the thin films are of a cerium oxide-based solid solution with few oxygen vacancies in the lattice and many absorptive oxygen bonds at the surface. Week ferromagnetic behaviors were evidenced by observed M–H hysteresis loops at room temperature. Furthermore, an evidence of relative ferromagnetic contributions was revealed by the temperature dependence of magnetization. It is believed that the ferromagnetic contributions exhibited in the M–H loops originate from the absorptive oxygen on the surface rather than the oxygen vacancies in the lattice.     

An improvement of the lattice theory of dislocation for a two-dimensional triangular crystal

The structure of dislocation in a two-dimensional triangular crystal has been studied theoretically on the basis of atomic interaction and lattice statics. The theory presented in this paper is an improvement to that published previously.Within a reasonable interaction approximation, a new dislocation equation is obtained, which remedies a fault existing in the lattice theory of dislocation. A better simplification of non-diagonal terms of the kernel is given. The solution of the new dislocation equation asymptotically becomes the same as that obtained in the elastic theory, and agrees with experimental data. It is found that the solution is formally identical with that proposed phenomenologically by Foreman et al, where the parameter can be chosen freely, but cannot uniquely determined from theory. Indeed, if the parameter in the expression of the solution is selected suitably, the expression can be well applied to describe the fine structure of the dislocation.     

A design strategy of achievable linear optics for a complex storage ring lattice

A design strategy is discussed in this paper,and it provides much convenience for effectively exploring achievable linear optics and globally investigating the flexibility of a complex lattice with super-periodicity.A matching method of fractional steps,which means separately finding the standard cell setting and the matching cell setting,is adopted to simplify the complexity of the linear beam optics design in the complex lattice.The multi-objective genetic algorithm is used to find most of all the stable linear optics,and reach a target solution after multi-generational propagation,both in the standard cell and the matching cell.A fitting algorithm with gradient information is used to restore the periodicity and symmetry of the lattice,and finely adjust the linear optics for further optimization.This design strategy is applied in the Shanghai Synchrotron Radiation Facility（SSRF） storage ring,and the results are presented.     

Folded novel accurate analytical and semi-analytical solutions of a generalized Calogero–Bogoyavlenskii–Schiff equation

This paper studies the analytical and semi-analytic solutions of the generalized Calogero–Bogoyavlenskii–Schiff(CBS) equation. This model describes the(2 + 1)–dimensional interaction between Riemann-wave propagation along the y-axis and the x-axis wave. The extended simplest equation(ESE) method is applied to the model, and a variety of novel solitarywave solutions is given. These solitary-wave solutions prove the dynamic behavior of soliton waves in plasma. The accuracy of the obtained solution is verified using a variational iteration(VI) semi-analytical scheme. The analysis and the match between the constructed analytical solution and the semi-analytical solution are sketched using various diagrams to show the accuracy of the solution we obtained. The adopted scheme's performance shows the effectiveness of the method and its ability to be applied to various nonlinear evolution equations.     

STOCHASTIC CONVERGENCE PROPERTIES OF THE TWO-PARAMETER ADAPTIVE LATTICE FILTERS

In this paper,the perturbation theory of the matrix eigenvalues,the operator spectral norm and a fixed-point theorem are used to research on the stochastic convergence properties of the two-parameter adaptivelattice filters.The relative difference value between the two parameters is a important factor.We obtainsome results:1.The singular values of the mean of elementary operator vary in a circle which has the radiusproportionated to ,they vary not only in magnitude,but also in direction due to; 2.Theeigenvalues of the mean square of elementary operator vary in a circle,its radius is proportional to also;3.The limits of stepsize are stricker than the ones of the one-parameter lattice filters;4.Zero-misadjustment can not be obtained.The misadjustment varys in the range where the center is equal to themisadjustment of the one-parameter lattice filters and the length is propertional to 5.The relationsbetween the misadjustment and order Nare neither linear nor exponential,but are in some conditions betweenthe     

Nonsensitive Homotopy-Pade Approach： Anharmonic Oscillators

In order to investigate a complicated physical system, it is convenient to consider a simple, easy to solve model, which is chosen to reflect as much physics as possible of the original system, as an ideal approximation. Motivated by this fundamental idea, we propose a novel asymptotic method, the nonsensitive homotopy-Pade approach. In this method, homotopy relations are constructed to link the original system with an ideal, solvable model. An artificial homotopy parameter is introduced to the homotopy relations as the normal perturbation parameter to generate the perturbation series, and is used to implement the Padd approximation. Meanwhile, some other auxiliary nonperturbative parameters, which are used to control the convergence of the perturbation series, are inserted to the approximants, and are fixed via the principle of minimal sensitivity. The method is used to study the eigenvalue problem of the quantum anharmonic oscillators. Highly accurate numerical results show its validity. Possible further studies on this method are also briefly discussed.     

Quantum phase transitions of fermionic atoms in an anisotropic triangular optical lattice

The effect of anisotropy caused by a confining potential on the properties of fermionic cold atoms in a triangular optical lattice is systematically investigated by using the dynamical cluster approximation combined with the continuous time quantum Monte–Carlo algorithm.The quantum phase diagrams which reflect the temperature–interaction relation and the competition between the anisotropic parameter and the interaction are presented with full consideration of the anisotropy of the system.Our results show that the system undergoes a transition from Fermi liquid to Mott insulator when the repulsive interaction reaches a critical value.The Kondo effect also can be observed in this system and the pseudogap is suppressed at low temperatures due to the Kondo effect.A feasible experiment protocol to observe these phenomena in an anisotropic triangular optical lattice with cold atoms is proposed,in which the hopping terms are closely related to the lattice confining potential and the atomic interaction can be adjusted via the Feshbach resonance.     

Microwave Atomic Clock in the Optical Lattice with Specific Frequency

A scheme for a microwave atomic clock is proposed for Cs or Rb atoms trapped in a blue detuned optical lattice. The ac Stark shift of the clock transition due to a trapping laser is calculated. We analyze it at some specific laser wavelength. Compared with the case of the fountain clock, the cavity related shifts, the collision shift and the Doppler effect are eliminated or suppressed dramatically in an atomic lattice clock. By analyzing various sources of clock uncertainty, a microwave atomic lattice clock with a high accuracy and small volume is feasible.     

Nonlinear optimization of the modern synchrotron radiation storage ring based on frequency map analysis

In this paper, we present a rule to improve the nonlinear solution with frequency map analysis （FMA）, and without frequently revisiting the optimization algorithm. Two aspects of FMA are emphasized. The first one is the tune shift with amplitude, which can be used to improve the solution of harmonic sextupoles, and thus obtain a large dynamic aperture. The second one is the tune diffusion rate, which can be used to select a quiet tune. Application of these ideas is carried out in the storage ring of the Shanghai Synchrotron Radiation Facility （SSRF）, and the detailed processes, as well as better solutions, are presented in this paper. Discussions about the nonlinear behaviors of off-momentum particles are also presented.     

On the rate of convergence of the Legendre spectral collocation method for multi-dimensional nonlinear Volterra-Fredholm integral equations

While the approximate solutions of one-dimensional nonlinear Volterra-Fredholm integral equations with smooth kermels are now well understood,no systematic studies of the numerical solutions of their multi-dimensional counterparts exist.In this paper,we provide an efficient numerical approach for the multi-dimensional nonlinear Volterra-Fredholm integral equations based on the multi-variate Legendre-collocation approach.Spectral collocation methods for multi-dimensional nonlinear integral equations are known to cause major difficulties from a convergence analysis point of view.Consequently,rigorous error estimates are provided in the weighted Sobolev space showing the exponential decay of the numerical errors.The existence and uniqueness of the numerical solution are established.Numerical experiments are provided to support the theoretical convergence analysis.The results indicate that our spectral collocation method is more flexible with better accuracy than the existing ones.     

Theory of passively mode-locked lasers with self-frequency shift and saturable absorber response

A general theoretical model for passively mode-locked lasers is presented, in which both the self-frequency shift and either a fast or a slow saturable absorber response are taken into account. An exact soliton-like solution and condition for its existence are obtained under a definite compatible condition. The stability of the solution is analyzed by using a variational method, and a parameter region, in which the solution is linearly stable, is acquired theoretically. To verify the theoretical predictions, a typical example is given for stable pulse propagation over a long distance. The numerical results show that the soliton-like solution is stable under some perturbations within the linearly stable region and an arbitrary Gaussian pulse converges to the exact soliton-like solution after evolution in a distance.     

Symmetry groups and spiral wave solution of a wave propagation equation

We study a third-order nonlinear evolution equation, which can be transformed to the modified KdV equation, using the Lie symmetry method. The Lie point symmetries and the one-dimensional optimal system of the symmetry algebras are determined. Those symmetries are some types of nonlocal symmetries or hidden symmetries of the modified KdV equation. The group-invariant solutions, particularly the travelling wave and spiral wave solutions, are discussed in detail, and a type of spiral wave solution which is smooth in the origin is obtained.     

Application of radial basis functions to evolution equations arising in image segmentation

In this paper, radial basis functions are used to obtain the solution of evolution equations which appear in variational level set method based image segmentation. In this method, radial basis functions are used to interpolate the implicit level set function of the evolution equation with a high level of accuracy and smoothness. Then, the original initial value problem is discretized into an interpolation problem. Accordingly, the evolution equation is converted into a set of coupled ordinary differential equations, and a smooth evolution can be retained. Compared with finite difference scheme based level set approaches, the complex and costly re-initialization procedure is unnecessary. Numerical examples are also given to show the efficiency of the method.     

Effects of wave propagation anisotropy on the wave focusing by negative refractive sonic crystal flat lenses

<正>In this study,wave propagation anisotropy in a triangular lattice crystal structure and its associated waveform shaping in a crystal structure are investigated theoretically.A directional variation in wave velocity inside a crystal structure is shown to cause bending wave envelcpes.The authors report that a triangular lattice sonic crystal possesses six numbers of a high symmetry direction,which leads to a wave convergence caused by wave velocity anisotropy inside the crystal.However,two of them are utilized mostly in wave focusing by an acoustic flat lens.Based on wave velocity anisotropy,the pseudo ideal imaging effect obtained in the second band of the flat lens is discussed.     

A mortar contact formulation using scaled boundary isogeometric analysis

正Contact analysis is recognized as being the most challenging problem in computational mechanics,because the functional system of contact problems is nonlinear and non-smooth and the convergence and accuracy of contact algorithms are difficult to guarantee.In the traditional finite element method(FEM)-based contact analysis[1,2],the contact body is spatially discretized,and the contact boundary is described using a low-order Lagrange interpolation polynomial.While