Page 1 Next

Displaying 1 – 20 of 44

Showing per page

Analysis of a semi-Lagrangian method for the spherically symmetric Vlasov-Einstein system

Philippe Bechouche, Nicolas Besse (2010)

ESAIM: Mathematical Modelling and Numerical Analysis

We consider the spherically symmetric Vlasov-Einstein system in the case of asymptotically flat spacetimes. From the physical point of view this system of equations can model the formation of a spherical black hole by gravitational collapse or describe the evolution of galaxies and globular clusters. We present high-order numerical schemes based on semi-Lagrangian techniques. The convergence of the solution of the discretized problem to the exact solution is proven and high-order error estimates...

Analysis of the accuracy and convergence of equation-free projection to a slow manifold

Antonios Zagaris, C. William Gear, Tasso J. Kaper, Yannis G. Kevrekidis (2009)

ESAIM: Mathematical Modelling and Numerical Analysis

In [C.W. Gear, T.J. Kaper, I.G. Kevrekidis and A. Zagaris, SIAM J. Appl. Dyn. Syst. 4 (2005) 711–732], we developed a class of iterative algorithms within the context of equation-free methods to approximate low-dimensional, attracting, slow manifolds in systems of differential equations with multiple time scales. For user-specified values of a finite number of the observables, the mth member of the class of algorithms ( m = 0 , 1 , ... ) finds iteratively an approximation of the appropriate zero of the (m+1)st...

Combining stochastic and deterministic approaches within high efficiency molecular simulations

Bruno Escribano, Elena Akhmatskaya, Jon Mujika (2013)

Open Mathematics

Generalized Shadow Hybrid Monte Carlo (GSHMC) is a method for molecular simulations that rigorously alternates Monte Carlo sampling from a canonical ensemble with integration of trajectories using Molecular Dynamics (MD). While conventional hybrid Monte Carlo methods completely re-sample particle’s velocities between MD trajectories, our method suggests a partial velocity update procedure which keeps a part of the dynamic information throughout the simulation. We use shadow (modified) Hamiltonians,...

Computing the differential of an unfolded contact diffeomorphism

Klaus Böhmer, Drahoslava Janovská, Vladimír Janovský (2003)

Applications of Mathematics

Consider a bifurcation problem, namely, its bifurcation equation. There is a diffeomorphism Φ linking the actual solution set with an unfolded normal form of the bifurcation equation. The differential D Φ ( 0 ) of this diffeomorphism is a valuable information for a numerical analysis of the imperfect bifurcation. The aim of this paper is to construct algorithms for a computation of D Φ ( 0 ) . Singularity classes containing bifurcation points with c o d i m 3 , c o r a n k = 1 are considered.

Continuation of invariant subspaces via the Recursive Projection Method

Vladimír Janovský, O. Liberda (2003)

Applications of Mathematics

The Recursive Projection Method is a technique for continuation of both the steady states and the dominant invariant subspaces. In this paper a modified version of the RPM called projected RPM is proposed. The modification underlines the stabilization effect. In order to improve the poor update of the unstable invariant subspace we have applied subspace iterations preconditioned by Cayley transform. A statement concerning the local convergence of the resulting method is proved. Results of numerical...

Discrete anisotropic curvature flow of graphs

Klaus Deckelnick, Gerhard Dziuk (2010)

ESAIM: Mathematical Modelling and Numerical Analysis

The evolution of n–dimensional graphs under a weighted curvature flow is approximated by linear finite elements. We obtain optimal error bounds for the normals and the normal velocities of the surfaces in natural norms. Furthermore we prove a global existence result for the continuous problem and present some examples of computed surfaces.

Edge of chaos in reaction diffusion CNN model

Angela Slavova, Ronald Tetzlaff (2017)

Open Mathematics

In this paper, we study the dynamics of a reaction-diffusion Cellular Nonlinear Network (RD-CNN) nodel in which the reaction term is represented by Brusselator cell. We investigate the RD-CNN dynamics by means of describing function method. Comparison with classical results for Brusselator equation is provided. Then we introduce a new RD-CNN model with memristor coupling, for which the edge of chaos regime in the parameter space is determined. Numerical simulations are presented for obtaining dynamic...

EKF-based dual synchronization of chaotic colpitts circuit and Chua’s circuit

Shaohua Hong, Zhiguo Shi, Kangsheng Chen (2008)

Kybernetika

In this paper, dual synchronization of a hybrid system containing a chaotic Colpitts circuit and a Chua’s circuit, connected by an additive white Gaussian noise (AWGN) channel, is studied via numeric simulations. The extended Kalman filter (EKF) is employed as the response system to achieve the dual synchronization. Two methods are proposed and investigated. The first method treats the combination of a Colpitts circuit and a Chua’s circuit as a higher- dimensional system, while the second method...

Energy-preserving Runge-Kutta methods

Elena Celledoni, Robert I. McLachlan, David I. McLaren, Brynjulf Owren, G. Reinout W. Quispel, William M. Wright (2009)

ESAIM: Mathematical Modelling and Numerical Analysis

We show that while Runge-Kutta methods cannot preserve polynomial invariants in general, they can preserve polynomials that are the energy invariant of canonical Hamiltonian systems.

Epidemiology of Dengue Fever: A Model with Temporary Cross-Immunity and Possible Secondary Infection Shows Bifurcations and Chaotic Behaviour in Wide Parameter Regions

Maíra Aguiar, Bob Kooi, Nico Stollenwerk (2008)

Mathematical Modelling of Natural Phenomena

Basic models suitable to explain the epidemiology of dengue fever have previously shown the possibility of deterministically chaotic attractors, which might explain the observed fluctuations found in empiric outbreak data. However, the region of bifurcations and chaos require strong enhanced infectivity on secondary infection, motivated by experimental findings of antibody-dependent-enhancement. Including temporary cross-immunity in such models, which is common knowledge among field researchers...

Currently displaying 1 – 20 of 44

Page 1 Next