Refine
Year of publication
- 2011 (38) (remove)
Document Type
- Article (29)
- Monograph/Edited Volume (4)
- Doctoral Thesis (3)
- Preprint (2)
Keywords
Institute
- Institut für Mathematik (38) (remove)
We consider the problem of discrete time filtering (intermittent data assimilation) for differential equation models and discuss methods for its numerical approximation. The focus is on methods based on ensemble/particle techniques and on the ensemble Kalman filter technique in particular. We summarize as well as extend recent work on continuous ensemble Kalman filter formulations, which provide a concise dynamical systems formulation of the combined dynamics-assimilation problem. Possible extensions to fully nonlinear ensemble/particle based filters are also outlined using the framework of optimal transportation theory.
We are interested in modeling some two-level population dynamics, resulting from the interplay of ecological interactions and phenotypic variation of individuals (or hosts) and the evolution of cells (or parasites) of two types living in these individuals. The ecological parameters of the individual dynamics depend on the number of cells of each type contained by the individual and the cell dynamics depends on the trait of the invaded individual. Our models are rooted in the microscopic description of a random (discrete) population of individuals characterized by one or several adaptive traits and cells characterized by their type. The population is modeled as a stochastic point process whose generator captures the probabilistic dynamics over continuous time of birth, mutation and death for individuals and birth and death for cells. The interaction between individuals (resp. between cells) is described by a competition between individual traits (resp. between cell types). We look for tractable large population approximations. By combining various scalings on population size, birth and death rates and mutation step, the single microscopic model is shown to lead to contrasting nonlinear macroscopic limits of different nature: deterministic approximations, in the form of ordinary, integro- or partial differential equations, or probabilistic ones, like stochastic partial differential equations or superprocesses. The study of the long time behavior of these processes seems very hard and we only develop some simple cases enlightening the difficulties involved.
We develop a multigrid, multiple time stepping scheme to reduce computational efforts for calculating complex stress interactions in a strike-slip 2D planar fault for the simulation of seismicity. The key elements of the multilevel solver are separation of length scale, grid-coarsening, and hierarchy. In this study the complex stress interactions are split into two parts: the first with a small contribution is computed on a coarse level, and the rest for strong interactions is on a fine level. This partition leads to a significant reduction of the number of computations. The reduction of complexity is even enhanced by combining the multigrid with multiple time stepping. Computational efficiency is enhanced by a factor of 10 while retaining a reasonable accuracy, compared to the original full matrix-vortex multiplication. The accuracy of solution and computational efficiency depend on a given cut-off radius that splits multiplications into the two parts. The multigrid scheme is constructed in such a way that it conserves stress in the entire half-space.
We analyze a general class of difference operators H(epsilon) = T(epsilon) + V(epsilon) on l(2)((epsilon Z)(d)), where V(epsilon) is a one-well potential and epsilon is a small parameter. We construct formal asymptotic expansions of WKB-type for eigenfunctions associated with the low lying eigenvalues of H(epsilon). These are obtained from eigenfunctions or quasimodes for the operator H(epsilon), acting on L(2)(R(d)), via restriction to the lattice (epsilon Z)(d).
Asymptotic first exit times of the chafee-infante equation with small heavy-tailed levy noise
(2011)
This article studies the behavior of stochastic reaction-diffusion equations driven by additive regularly varying pure jump Levy noise in the limit of small noise intensity. It is shown that the law of the suitably renormalized first exit times from the domain of attraction of a stable state converges to an exponential law of parameter 1 in a strong sense of Laplace transforms, including exponential moments. As a consequence, the expected exit times increase polynomially in the inverse intensity, in contrast to Gaussian perturbations, where this growth is known to be of exponential rate.
In this paper, we propose a derivative-free method for recovering symmetric and non-symmetric potential functions of inverse Sturm-Liouville problems from the knowledge of eigenvalues. A class of boundary value methods obtained as an extension of Numerov's method is the major tool for approximating the eigenvalues in each Broyden iteration step. Numerical examples demonstrate that the method is able to reduce the number of iteration steps, in particular for non-symmetric potentials, without accuracy loss.
Processes with independent increments are characterized via a duality formula, including Malliavin derivative and difference operators. This result is based on a characterization of infinitely divisible random vectors by a functional equation. A construction of the difference operator by a variational method is introduced and compared to approaches used by other authors for L´evy processes involving the chaos decomposition. Finally we extend our method to characterize infinitely divisible random measures.
The problem of an ensemble Kalman filter when only partial observations are available is considered. In particular, the situation is investigated where the observational space consists of variables that are directly observable with known observational error, and of variables of which only their climatic variance and mean are given. To limit the variance of the latter poorly resolved variables a variance-limiting Kalman filter (VLKF) is derived in a variational setting. The VLKF for a simple linear toy model is analyzed and its range of optimal performance is determined. The VLKF is explored in an ensemble transform setting for the Lorenz-96 system, and it is shown that incorporating the information of the variance of some unobservable variables can improve the skill and also increase the stability of the data assimilation procedure.