Search results for "ordinary differential equation"
showing 10 items of 98 documents
Performance modeling of epidemic routing
2006
In this paper, we develop a rigorous, unified framework based on ordinary differential equations (ODEs) to study epidemic routing and its variations. These ODEs can be derived as limits of Markovian models under a natural scaling as the number of nodes increases. While an analytical study of Markovian models is quite complex and numerical solution impractical for large networks, the corresponding ODE models yield closed-form expressions for several performance metrics of interest, and a numerical solution complexity that does not increase with the number of nodes. Using this ODE approach, we investigate how resources such as buffer space and the number of copies made for a packet can be tra…
Adaptive Gaussian particle method for the solution of the Fokker-Planck equation
2012
The Fokker-Planck equation describes the evolution of the probability density for a stochastic ordinary differential equation (SODE). A solution strategy for this partial differential equation (PDE) up to a relatively large number of dimensions is based on particle methods using Gaussians as basis functions. An initial probability density is decomposed into a sum of multivariate normal distributions and these are propagated according to the SODE. The decomposition as well as the propagation is subject to possibly large numeric errors due to the difficulty to control the spatial residual over the whole domain. In this paper a new particle method is derived, which allows a deterministic error…
On attracting sets in artificial networks: cross activation
2018
Mathematical models of artificial networks can be formulated in terms of dynamical systems describing the behaviour of a network over time. The interrelation between nodes (elements) of a network is encoded in the regulatory matrix. We consider a system of ordinary differential equations that describes in particular also genomic regulatory networks (GRN) and contains a sigmoidal function. The results are presented on attractors of such systems for a particular case of cross activation. The regulatory matrix is then of particular form consisting of unit entries everywhere except the main diagonal. We show that such a system can have not more than three critical points. At least n–1 eigenvalu…
The McShane, PU and Henstock integrals of Banach valued functions
2002
Some relationships between the vector valued Henstock and McShane integrals are investigated. An integral for vector valued functions, defined by means of partitions of the unity (the PU-integral) is studied. In particular it is shown that a vector valued function is McShane integrable if and only if it is both Pettis and PU-integrable. Convergence theorems for the Henstock variational and the PU integrals are stated. The families of multipliers for the Henstock and the Henstock variational integrals of vector valued functions are characterized.
Constructing transient response probability density of non-linear system through complex fractional moments
2014
Abstract The probability density function for transient response of non-linear stochastic system is investigated through the stochastic averaging and Mellin transform. The stochastic averaging based on the generalized harmonic functions is adopted to reduce the system dimension and derive the one-dimensional Ito stochastic differential equation with respect to amplitude response. To solve the Fokker–Plank–Kolmogorov equation governing the amplitude response probability density, the Mellin transform is first implemented to obtain the differential relation of complex fractional moments. Combining the expansion form of transient probability density with respect to complex fractional moments an…
Estimation of the mean crystal size and the moments of the crystal size distribution in batch crystallization processes
2016
International audience; A cascade high gain observer is designed to estimate the first four leading moments of the crystal size distribution (CSD) and the mean crystal size in batch crystallization processes. The proposed observer is based on a well-known transformation of the partial differential equation describing the CSD to a set of ordinary differential equations (the method of moments). Due to numerical difficulties resulting from the important differences in the magnitudes of the moments, a set of new variables is computed to allow a good estimation of the moments and thus the mean crystal size. In this work, only solute concentration and crystallizer temperature are used to estimate…
High-order Runge–Kutta–Nyström geometric methods with processing
2001
Abstract We present new families of sixth- and eighth-order Runge–Kutta–Nystrom geometric integrators with processing for ordinary differential equations. Both the processor and the kernel are composed of explicitly computable flows associated with non trivial elements belonging to the Lie algebra involved in the problem. Their efficiency is found to be superior to other previously known algorithms of equivalent order, in some case up to four orders of magnitude.
On the Stochastic Response of a Fractionally-damped Duffing Oscillator
2012
A numerical method is presented to compute the response of a viscoelastic Duffing oscillator with fractional derivative damping, subjected to a stochastic input. The key idea involves an appropriate discretization of the fractional derivative, based on a preliminary change of variable, that allows to approximate the original system by an equivalent system with additional degrees of freedom, the number of which depends on the discretization of the fractional derivative. Unlike the original system that, due to the presence of the fractional derivative, is governed by non-ordinary differential equations, the equivalent system is governed by ordinary differential equations that can be readily h…
The MAST-edge centred lumped scheme for the flow simulation in variably saturated heterogeneous porous media
2012
A novel methodology is proposed for the solution of the flow equation in a variably saturated heterogeneous porous medium. The computational domain is descretized using triangular meshes and the governing PDEs are discretized using a lumped in the edge centres numerical technique. The dependent unknown variable of the problem is the piezometric head. A fractional time step methodology is applied for the solution of the original system, solving consecutively a prediction and a correction problem. A scalar potential of the flow field exists and in the prediction step a MArching in Space and Time (MAST) formulation is applied for the sequential solution of the Ordinary Differential Equation of…