Discretized Multisplitting AOR Waveform Relaxation Algorithms for Initial Value Problem of Systems of ODEs

The multisplitting algorithm for solving large systems of ordinary differential equations on parallel computers was introduced by Jeltsch and Pohl in [1]. On fixed time intervals conver gence results could be derived if the subsystems are solving exactly.Firstly,in theis paper,we deal with an extension of the waveform relaxation algorithm by us ing multisplittin AOR method based on an overlapping block decomposition. We restricted our selves to equidistant timepoints and dealed with the case that an implicit integration method was used to solve the subsystems numerically in parallel. Then we have proved convergence of multi splitting AOR waveform relaxation algorithm on a fixed window containing a finite number of timepoints.

Computer Based Modeling of Vegetation Propagation across a Disturbed Space Linked to GIS

The paper is focused on computer simulation of natural vegetation propagation across two selected disturbed sites. Two sites located in the different environments, the abandoned sedimentation basin of a former pyrite ore mine and the ash deposits of a power station, were selected to illustrate the proposed spatio-temporal model. Aerial images assisted in identifying and monitoring the progress in the propagation of vegetation. Analysis of the aerial images was based on varying vegetation coverage explored by classification algorithms. A new approach is proposed entailing coupling of a local dynamic model and a spatial model for vegetation propagation. The local dynamic model describes vegetation growth using a logistic growth approach based on delayed variables. Vegetation propagation is described by rules related to seed and its dispersal phenomena on a local scale and on the scale of outlying spreading. The disturbed sites are divided into a grid of microsites. Each microsite is represented by a 5 m x 5 m square. A state variable in each microsite indicates the relative vegetation density on a scale from 0 （no vegetation） to 1 （long-term maximum of vegetation density）. Growth, local vegetation propagation and the effects of outlying vegetation propagation in each cell are described by an ordinary differential equation with delayed state variables. The grid of cells forms a set of ordinary differential equations. The abandoned sedimentation basin and the ash deposits are represented by grids of 185 x 345 and 212 x 266 cells, respectively. A few case-oriented studies are provided to show various predictions of vegetation propagation across two selected disturbed sites. The first case study simulates vegetation growing without spatial propagations and delayed variables in the spatio-temporal model. The second and the third case studies extend the previous study by including local and outlying vegetation propagation, respectively. The fourth case study explores delayed impacts in the logistic growth term and the delayed outcome by vegetation propagation across the disturbed space. The performed case-oriented studies confirm the applicability of the proposed spatio-temporal model to predict vegetation propagation in short-term successions and to estimate approximate vegetation changes in long-term development. As a result, it can be concluded that remotely sensed data are a valuable source of information for estimates of model parameters and provide an effective method for monitoring the progress of vegetation propagation across the selected sites, spaces disturbed by human activities.

A Novel Approach for Numerical Computation of Fokker-Planck Equation

This paper present an implementation of＂modified cubic B-spline differential quadrature method （MCB-DQM）＂ proposed by Arora ＆ Singh （Applied Mathematics and Computation Vol. 224（1） （2013） 161-177） for numerical computation of Fokker-Planck equations. The modified cubic B-splines are used as set of basis functions in the differential quadrature to compute the weighting coefficients for the spatial derivatives, which reduces Fokker-Planck equation into system of first-order ordinary differential equations （ODEs）, in time. The well known SSP-RK43 scheme is then applied to solve the resulting system of ODEs. The efficiency of proposed method has been confirmed by three examples having their exact solutions. This shows that MCB-DQM results are capable of achieving high accuracy. Advantage of the scheme is that it can be applied very smoothly to solve the linear or nonlinear physical problems, and a very less storage space is required which causes less accumulation of numerical errors.

基于目标流线的曲线方向寻优

本文对于无约束问题提出了沿曲线寻优的思想,推导了确定寻优曲线的常微分方程组。对该常微分方程组研制了近似解析和数值的实用算法.借助对偶规划,这一方法由无约束问题推广到约束问题.

Solution procedure of residue harmonic balance method and its applications

This paper presents a simple and rigorous solution procedure of residue harmonic balance for predicting the accurate approximation of certain autonomous ordinary differential systems.In this solution procedure,no small parameter is assumed.The harmonic residue of balance equation is separated in two parts at each step.The first part has the same number of Fourier terms as the present order of approximation and the remaining part is used in the subsequent improvement.The corrections are governed by linear ordinary differential equation so that they can be solved easily by means of harmonic balance method again.Three kinds of different differential equations involving general,fractional and delay ordinary differential systems are given as numerical examples respectively.Highly accurate limited cycle frequency and amplitude are captured.The results match well with the exact solutions or numerical solutions for a wide range of control parameters.Comparison with those available shows that the residue harmonic balance solution procedure is very effective for these autonomous differential systems.Moreover,the present method works not only in predicting the amplitude but also the frequency of bifurcated period solution for delay ordinary differential equation.

Construction of Dulac functions for mathematical models in population biology

In this paper, we present a method for constructing a Dulac function for mathematical models in population biology, in the form of systems of ordinary differential equations in the plane.

Nonlinear Exact Solutions of the 2-Dimensional Rotational Euler Equations for the Incompressible Fluid

In this paper, the Clarkson-Kruskal direct approach is employed to investigate the exact solutions of the 2-dimensionai rotationai Euler equations for the incompressible fluid. The application of the method leads to a system of completely solvable ordinary differential equations. Several special cases are discussed and novel nonlinear exact solutions with respect to variables x and y are obtained. It is＇of interest to notice that the pressure p is obtained by the second kind of curvilinear integral and the coefficients of the nonlinear solutions are solitary wave type functions like tanh（ kt /2 ） and sech （kt/2） due to the rotational parameter k ≠ O. Such phenomenon never appear in the classical Euler equations wherein the Coriolis force arising from the gravity and Earth＇s rotation is ignored. Finally, illustrative numerical figures are attached to show the behaviors that the exact solutions may exhibit.

Separation Transformation and New Exact Solutions of the (N+1)-dimensional Dispersive Double sine-Gordon Equation

In this paper, the separation transformation approach is extended to the （N ＋ 1）-dimensional dispersive double sine-Gordon equation arising in many physical systems such as the spin dynamics in the B phase of SHe superfluid. This equation is first reduced to a set of partial differential equations and a nonlinear ordinary differential equation. Then the general solutions of the set of partial differential equations are obta/ned and the nonlinear ordinary differential equation is solved by F-expansion method. Finally, many new exact solutions of the （N ＋ 1）-dimensional dispersive double sine-Gordon equation are constructed explicitly via the separation transformation. For the case of N 〉 2, there is an arbitrary function in the exact solutions, which may reveal more novel nonlinear structures in the high-dimensional dispersive double sine-Gordon equation.

THE DYNAMICAL MODELS OF ACTIVATED SLUDGE SYSTEM： STOCHASTIC CELLULAR AUTOMATON AND DIFFERENTIAL EQUATIONS

A stochastic celhflar automaton （CA） model for activated sludge system （ASS） is for- mulated by a series of transition functions upon realistic treatment processes, and it is tested by comparing with ordinary differential equations （ODEs） of ASS. CA system performed by empirical parameters can reflect the characteristics of fluctuation, com- plexity and strong non-linearity of ASS. The results show that the predictions of CA are approximately similar to the dynamical behaviors of ODEs. Based on the extreme experimental system with complete cell recycle in model validation, the dynamics of biomass and substrate are predicted accurately by CA, but the large errors exist in ODEs except for integrating more spatially complicated factors. This is due to that the strong mechanical stress from spatial crowding effect is ignored in ODEs, while CA system as a spatially explicit model takes account of local interactions. Despite its extremely simple structure, CA still can capture the essence of ASS better than ODEs, thus it would be very useful in predicting long-term dynamics in other similar systems.