Application of an Acceleration Scheme for an Age-Structured Diffusion Model

In this paper we propose an optimized algorithm, which is faster compared to previously described finite difference acceleration scheme, namely the Modified Super-Time-Stepping (Modified STS) scheme for age- structured population models with diffusion.

Increasing the Calculation Speed of an Acceleration Scheme for an Age-Structured Diffusion Model

Дойчин Бояджиев, Галена Пеловска - В статията се предлага оптимизиран алгоритъм, който е по-бърз в сравнение с по- рано описаната ускорена (модифицирана STS) диференчна схема за възрастово структуриран популационен модел с дифузия. Запазвайки апроксимацията на модифицирания STS алгоритъм, изчислителното времето се намаля почти два пъти. Това прави оптимизирания метод по-предпочитан за задачи с нелинейност или с по-висока размерност.

Asynchronous growth and competition in a two-sex age-structured population model

Asynchronous exponential growth has been extensively studied in population dynamics. In this paper we find out the asymptotic behaviour in a non-linear age-dependent model which takes into account sexual reproduction interactions. The main feature of our model is that the non-linear process converges to a linear one as the solution becomes large, so that the population undergoes asynchronous growth. The steady states analysis and the corresponding stability analysis are completely made and are summarized in a bifurcation diagram according to the parameter R0. Furthermore the effect of intraspecific competition is taken into account, leading to complex dynamics around steady states.

Tuberculosis in Cape Town: An age-structured transmission model.

BACKGROUND Tuberculosis (TB) is the leading cause of death in South Africa. The burden of disease varies by age, with peaks in TB notification rates in the HIV-negative population at ages 0-5, 20-24, and 45-49 years. There is little variation between age groups in the rates in the HIV-positive population. The drivers of this age pattern remain unknown. METHODS We developed an age-structured simulation model of Mycobacterium tuberculosis (Mtb) transmission in Cape Town, South Africa. We considered five states of TB progression: susceptible, infected (latent TB), active TB, treated TB, and treatment default. Latently infected individuals could be re-infected; a previous Mtb infection slowed progression to active disease. We further considered three states of HIV progression: HIV negative, HIV positive, on antiretroviral therapy. To parameterize the model, we analysed treatment outcomes from the Cape Town electronic TB register, social mixing patterns from a Cape Town community and used literature estimates for other parameters. To investigate the main drivers behind the age patterns, we conducted sensitivity analyses on all parameters related to the age structure. RESULTS The model replicated the age patterns in HIV-negative TB notification rates of Cape Town in 2009. Simulated TB notification rate in HIV-negative patients was 1000/100,000 person-years (pyrs) in children aged <5 years and decreased to 51/100,000 in children 5-15 years. The peak in early adulthood occurred at 25-29 years (463/100,000 pyrs). After a subsequent decline, simulated TB notification rates gradually increased from the age of 30 years. Sensitivity analyses showed that the dip after the early adult peak was due to the protective effect of latent TB and that retreatment TB was mainly responsible for the rise in TB notification rates from the age of 30 years. CONCLUSION The protective effect of a first latent infection on subsequent infections and the faster progression in previously treated patients are the key determinants of the age-structure of TB notification rates in Cape Town.

Harvest control rule vs optimal harvesting of an age-structured population: The case of the Ibero-Atlantic Sardine Fishery

Research Masters

A stochastic spatially structured epidemic model with diffusive processes

We developed a stochastic lattice model to describe the vector-borne disease (like yellow fever or dengue). The model is spatially structured and its dynamical rules take into account the diffusion of vectors. We consider a bipartite lattice, forming a sub-lattice of human and another occupied by mosquitoes. At each site of lattice we associate a stochastic variable that describes the occupation and the health state of a single individual (mosquito or human). The process of disease transmission in the human population follows a similar dynamic of the Susceptible-Infected-Recovered model (SIR), while the disease transmission in the mosquito population has an analogous dynamic of the Susceptible-Infected-Susceptible model (SIS) with mosquitos diffusion. The occurrence of an epidemic is directly related to the conditional probability of occurrence of infected mosquitoes (human) in the presence of susceptible human (mosquitoes) on neighborhood. The probability of diffusion of mosquitoes can facilitate the formation of pairs Susceptible-Infected enabling an increase in the size of the epidemic. Using an asynchronous dynamic update, we study the disease transmission in a population initially formed by susceptible individuals due to the introduction of a single mosquito (human) infected. We find that this model exhibits a continuous phase transition related to the existence or non-existence of an epidemic. By means of mean field approximations and Monte Carlo simulations we investigate the epidemic threshold and the phase diagram in terms of the diffusion probability and the infection probability.

Why do scientists migrate? A diffusion model

The article is intended to improve our understanding of the reasons underlying the intellectual migration of scientists from well-known cognitive domains to nascent scientific fields. To that purpose we present, first, a number of findings from the sociology of science that give different insights about this phenomenon. We then attempt to bring some of these insights together under the conceptual roof of an actor-based approach linking expected utility and diffusion theory. Intellectual migration is regarded as the rational choice of scientists who decide under uncertainty and on the base of a number of decision-making variables, which define probabilities, costs, and benefits of the migration.

A Hull and White formula for a general stochastic volatility jump-diffusion model with applications to the study of the short-time behavior of the implied volatility

In this paper, generalizing results in Alòs, León and Vives (2007b), we see that the dependence of jumps in the volatility under a jump-diffusion stochastic volatility model, has no effect on the short-time behaviour of the at-the-money implied volatility skew, although the corresponding Hull and White formula depends on the jumps. Towards this end, we use Malliavin calculus techniques for Lévy processes based on Løkka (2004), Petrou (2006), and Solé, Utzet and Vives (2007).

Analysis of a finite volume method for a cross-diffusion model in population dynamics

The main goal of this paper is to propose a convergent finite volume method for a reactionâeuro"diffusion system with cross-diffusion. First, we sketch an existence proof for a class of cross-diffusion systems. Then the standard two-point finite volume fluxes are used in combination with a nonlinear positivity-preserving approximation of the cross-diffusion coefficients. Existence and uniqueness of the approximate solution are addressed, and it is also shown that the scheme converges to the corresponding weak solution for the studied model. Furthermore, we provide a stability analysis to study pattern-formation phenomena, and we perform two-dimensional numerical examples which exhibit formation of nonuniform spatial patterns. From the simulations it is also found that experimental rates of convergence are slightly below second order. The convergence proof uses two ingredients of interest for various applications, namely the discrete Sobolev embedding inequalities with general boundary conditions and a space-time $L^1$ compactness argument that mimics the compactness lemma due to Kruzhkov. The proofs of these results are given in the Appendix.

Diffusion Model Applied to Postfeeding Larval Dispersal in Blowflies (Diptera: Calliphoridae)

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)