BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Giovanni Samaey

Giovanni Samaey appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysisphysics.comp-phcond-mat.mes-hallmath.APmath.DSnlin.PSQuantitative MethodsTissues and Organs

Catalog footprint

What is connected

17works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2021arXiv

Projective Integration Schemes for Hyperbolic Moment Equations

In this paper, we apply projective integration methods to hyperbolic moment models of the Boltzmann equation and the BGK equation, and investigate the numerical properties of the resulting scheme. Projective integration is an explicit, asymptotic-preserving scheme that is tailored to problems with large spectral gaps between slow and (one or many) fast eigenvalue clusters of the model. The spectral analysis of a linearized moment model clearly shows spectral gaps and reveals the multi-scale nature of the model for which projective integration is a matching choice. The combination of the non-intrusive projective integration method with moment models allows for accurate, but efficient simulations with significant speedup, as demonstrated using several 1D and 2D test cases with different collision terms, collision frequencies and relaxation times.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A micro-macro Markov chain Monte Carlo method for molecular dynamics using reaction coordinate proposals I: direct reconstruction

We introduce a new micro-macro Markov chain Monte Carlo method (mM-MCMC) to sample invariant distributions of molecular dynamics systems that exhibit a time-scale separation between the microscopic (fast) dynamics, and the macroscopic (slow) dynamics of some low-dimensional set of reaction coordinates. The algorithm enhances exploration of the state space in the presence of metastability by allowing larger proposal moves at the macroscopic level, on which a conditional accept-reject procedure is applied. Only when the macroscopic proposal is accepted, the full microscopic state is reconstructed from the newly sampled reaction coordinate value and is subjected to a second accept/reject procedure. The computational gain stems from the fact that most proposals are rejected at the macroscopic level, at low computational cost, while microscopic states, once reconstructed, are almost always accepted. We analytically show convergence and discuss the rate of convergence of the proposed algorithm, and numerically illustrate its efficiency on a number of standard molecular test cases. We also investigate the effect of the choice of different numerical parameters on the efficiency of the resulting mM-MCMC method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Multilevel Monte Carlo Asymptotic-Preserving Particle Method for Kinetic Equations in the Diffusion Limit

We propose a multilevel Monte Carlo method for a particle-based asymptotic-preserving scheme for kinetic equations. Kinetic equations model transport and collision of particles in a position-velocity phase-space. With a diffusive scaling, the kinetic equation converges to an advection-diffusion equation in the limit of zero mean free path. Classical particle-based techniques suffer from a strict time-step restriction to maintain stability in this limit. Asymptotic-preserving schemes provide a solution to this time step restriction, but introduce a first-order error in the time step size. We demonstrate how the multilevel Monte Carlo method can be used as a bias reduction technique to perform accurate simulations in the diffusive regime, while leveraging the reduced simulation cost given by the asymptotic-preserving scheme. We describe how to achieve the necessary correlation between simulation paths at different levels and demonstrate the potential of the approach via numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A high-order relaxation method with projective integration for solving nonlinear systems of hyperbolic conservation laws

We present a general, high-order, fully explicit relaxation scheme which can be applied to any system of nonlinear hyperbolic conservation laws in multiple dimensions. The scheme consists of two steps. In a first (relaxation) step, the nonlinear hyperbolic conservation law is approximated by a kinetic equation with stiff BGK source term. Then, this kinetic equation is integrated in time using a projective integration method. After taking a few small (inner) steps with a simple, explicit method (such as direct forward Euler) to damp out the stiff components of the solution, the time derivative is estimated and used in an (outer) Runge-Kutta method of arbitrary order. We show that, with an appropriate choice of inner step size, the time step restriction on the outer time step is similar to the CFL condition for the hyperbolic conservation law. Moreover, the number of inner time steps is also independent of the stiffness of the BGK source term. We discuss stability and consistency, and illustrate with numerical results (linear advection, Burgers' equation and the shallow water and Euler equations) in one and two spatial dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

DDE-BIFTOOL Manual - Bifurcation analysis of delay differential equations

DDEBIFTOOL is a collection of Matlab routines for numerical bifurcation analysis of systems of delay differential equations with discrete constant and state-dependent delays. The package supports continuation and stability analysis of steady state solutions and periodic solutions. Further one can compute and continue several local and global bifurcations: fold and Hopf bifurcations of steady states; folds, period doublings and torus bifurcations of periodic orbits; and connecting orbits between equilibria. To analyse the stability of steady state solutions, approximations are computed to the rightmost, stability-determining roots of the characteristic equation which can subsequently be used as starting values in a Newton procedure. For periodic solutions, approximations to the Floquet multipliers are computed. The manual describes the structure of the package, its routines, and its data and method parameter structures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Telescopic projective integration for kinetic equations with multiple relaxation times

We study a general, high-order, fully explicit numerical method for simulating kinetic equations with a BGK-type collision model with multiple relaxation times. In that case, the problem is stiff and its spectrum consists of multiple separated eigenvalue clusters. Projective integration methods are explicit integration schemes that first take a few small (inner) steps with a simple, explicit method, after which the solution is extrapolated forward in time over a large (outer) time step. These are very efficient schemes, provided there are only two clusters of eigenvalues. Telescopic projective integration methods generalize the idea of projective integration methods by constructing a hierarchy of projective levels. Here, we show how telescopic projective integration methods can be used to efficiently integrate kinetic equations with multiple relaxation times. We show that the required number of projective levels depends on the number of clusters, which in turn depends on the stiffness of the BGK source term. The size of the outer level time step only depends on the slowest time scale present in the model and is independent of the stiffness of the problem. We discuss stability and illustrate the approach with simulations in one and two spatial dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Variance-reduced multiscale simulation of slow-fast stochastic differential equations

We study a variance reduction strategy based on control variables for simulating the averaged macroscopic behavior of a stochastic slow-fast system. We assume that this averaged behavior can be written in terms of a few slow degrees of freedom, and that the fast dynamics is ergodic for every fixed value of the slow variable. The time derivative for the averaged dynamics can then be approximated by a Markov chain Monte Carlo method. The variance-reduced scheme that is introduced here uses the previous time instant as a control variable. We analyze the variance and bias of the proposed estimator and illustrate its performance when applied to a linear and nonlinear model problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variance-reduced simulation of stochastic agent-based models for tumor growth

We investigate a hybrid PDE/Monte Carlo technique for the variance reduced simulation of an agent-based multiscale model for tumor growth. The variance reduction is achieved by combining a simulation of the stochastic agent-based model on the microscopic scale with a deterministic solution of a simplified (coarse) partial differential equation (PDE) on the macroscopic scale as a control variable. We show that this technique is able to significantly reduce the variance with only the (limited) additional computational cost associated with the deterministic solution of the coarse PDE. We illustrate the performance with numerical experiments in different regimes, both in the avascular and vascular stage of tumor growth.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A high-order asymptotic-preserving scheme for kinetic equations using projective integration

We investigate a high-order, fully explicit, asymptotic-preserving scheme for a kinetic equation with linear relaxation, both in the hydrodynamic and diffusive scalings in which a hyperbolic, resp. parabolic, limiting equation exists. The scheme first takes a few small (inner) steps with a simple, explicit method (such as direct forward Euler) to damp out the stiff components of the solution and estimate the time derivative of the slow components. These estimated time derivatives are then used in an (outer) Runge-Kutta method of arbitrary order. We show that, with an appropriate choice of inner step size, the time-step restriction on the outer time step is similar to the stability condition for the limiting macroscopic equation. Moreover, the number of inner time steps is also independent of the scaling parameter. We analyse stability and consistency, and illustrate with numerical results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Noise reduction in coarse bifurcation analysis of stochastic agent-based models: an example of consumer lock-in

We investigate coarse equilibrium states of a fine-scale, stochastic agent-based model of consumer lock-in in a duopolistic market. In the model, agents decide on their next purchase based on a combination of their personal preference and their neighbours' opinions. For agents with independent identically-distributed parameters and all-to-all coupling, we derive an analytic approximate coarse evolution-map for the expected average purchase. We then study the emergence of coarse fronts when spatial segregation is present in the relative perceived quality of products. We develop a novel Newton-Krylov method that is able to compute accurately and efficiently coarse fixed points when the underlying fine-scale dynamics is stochastic. The main novelty of the algorithm is in the elimination of the noise that is generated when estimating Jacobian-vector products using time-integration of perturbed initial conditions. We present numerical results that demonstrate the convergence properties of the numerical method, and use the method to show that macroscopic fronts in this model destabilise at a coarse symmetry-breaking bifurcation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the acceleration of spatially distributed agent-based computations: a patch dynamics scheme

In recent years, individual-based/agent-based modeling has been applied to study a wide range of applications, ranging from engineering problems to phenomena in sociology, economics and biology. Simulating such agent-based models over extended spatiotemporal domains can be prohibitively expensive due to stochasticity and the presence of multiple scales. Nevertheless, many agent-based problems exhibit smooth behavior in space and time on a macroscopic scale, suggesting that a useful coarse-grained continuum model could be obtained. For such problems, the equation-free framework [16-18] can significantly reduce the computational cost. Patch dynamics is an essential component of this framework. This scheme is designed to perform numerical simulations of an unavailable macroscopic equation on macroscopic time and length scales; it uses appropriately initialized simulations of the fine-scale agent-based model in a number of small "patches", which cover only a fraction of the spatiotemporal domain. In this work, we construct a finite-volume-inspired conservative patch dynamics scheme and apply it to a financial market agent-based model based on the work of Omurtag and Sirovich [22]. We first apply our patch dynamics scheme to a continuum approximation of the agent-based model, to study its performance and analyze its accuracy. We then apply the scheme to the agent-based model itself. Our computational experiments indicate that here, typically, the patch dynamics-based simulation requires only 20% of the full agent-based simulation in space, and need occur over only 10% of the temporal domain.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A micro-macro parareal algorithm: application to singularly perturbed ordinary differential equations

We introduce a micro-macro parareal algorithm for the time-parallel integration of multiscale-in-time systems. The algorithm first computes a cheap, but inaccurate, solution using a coarse propagator (simulating an approximate slow macroscopic model), which is iteratively corrected using a fine-scale propagator (accurately simulating the full microscopic dynamics). This correction is done in parallel over many subintervals, thereby reducing the wall-clock time needed to obtain the solution, compared to the integration of the full microscopic model. We provide a numerical analysis of the algorithm for a prototypical example of a micro-macro model, namely singularly perturbed ordinary differential equations. We show that the computed solution converges to the full microscopic solution (when the parareal iterations proceed) only if special care is taken during the coupling of the microscopic and macroscopic levels of description. The convergence rate depends on the modeling error of the approximate macroscopic model. We illustrate these results with numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Towards an efficient multiscale modeling of low-dimensional reactive systems: study of numerical closure procedures

In this paper, we present a study on how to develop an efficient multiscale simulation strategy for the dynamics of chemically active systems on low-dimensional supports. Such reactions are encountered in a wide variety of situations, ranging from heterogeneous catalysis to electrochemical or (membrane) biological processes, to cite a few. We analyzed in this context different techniques within the framework of an important multiscale approach known as the equation free method (EFM), which "bridges the multiscale gap" by building microscopic configurations using macroscopic-level information only. We hereby considered two simple reactive processes on a one-dimensional lattice, the simplicity of which allowed for an in-depth understanding of the parameters controlling the efficiency of this approach. We demonstrate in particular that it is not enough to base the EFM on the time evolution of the average concentrations of particles on the lattice, but that the time evolution of clusters of particles has to be included as well. We also show how important it is for the accuracy of this method to carefully choose the procedure with which microscopic states are constructed, starting from the measured macroscopic quantities. As we also demonstrate that some errors cannot be corrected by increasing the number of observed macroscopic variables, this work points towards which procedures should be used in order to generate efficient and reliable multiscale simulations of these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A micro/macro algorithm to accelerate Monte Carlo simulation of stochastic differential equations

We present and analyze a micro/macro acceleration technique for the Monte Carlo simulation of stochastic differential equations (SDEs) in which there is a separation between the (fast) time-scale on which individual trajectories of the SDE need to be simulated and the (slow) time-scale on which we want to observe the (macroscopic) function of interest. The method performs short bursts of microscopic simulation using an ensemble of SDE realizations, after which the ensemble is restricted to a number of macroscopic state variables. The resulting macroscopic state is then extrapolated forward in time and the ensemble is projected onto the extrapolated macroscopic state. We provide a first analysis of its convergence in terms of extrapolation time step and number of macroscopic state variables. The effects of the different approximations on the resulting error are illustrated via numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Individual-based models for bacterial chemotaxis in the diffusion asymptotics

We discuss velocity-jump models for chemotaxis of bacteria with an internal state that allows the velocity jump rate to depend on the memory of the chemoattractant concentration along their path of motion. Using probabilistic techniques, we provide a pathwise result that shows that the considered process converges to an advection-diffusion process in the (long-time) diffusion limit. We also (re-)prove using the same approach that the same limiting equation arises for a related, simpler process with direct sensing of the chemoattractant gradient. Additionally, we propose a time discretization technique that retains these diffusion limits exactly, i.e., without error that depends on the time discretization. In the companion paper \cite{variance}, these results are used to construct a coupling technique that allows numerical simulation of the process with internal state with asymptotic variance reduction, in the sense that the variance vanishes in the diffusion limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simulating individual-based models of bacterial chemotaxis with asymptotic variance reduction

We discuss variance reduced simulations for an individual-based model of chemotaxis of bacteria with internal dynamics. The variance reduction is achieved via a coupling of this model with a simpler process in which the internal dynamics has been replaced by a direct gradient sensing of the chemoattractants concentrations. In the companion paper \cite{limits}, we have rigorously shown, using a pathwise probabilistic technique, that both processes converge towards the same advection-diffusion process in the diffusive asymptotics. In this work, a direct coupling is achieved between paths of individual bacteria simulated by both models, by using the same sets of random numbers in both simulations. This coupling is used to construct a hybrid scheme with reduced variance. We first compute a deterministic solution of the kinetic density description of the direct gradient sensing model; the deviations due to the presence of internal dynamics are then evaluated via the coupled individual-based simulations. We show that the resulting variance reduction is \emph{asymptotic}, in the sense that, in the diffusive asymptotics, the difference between the two processes has a variance which vanishes according to the small parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Asymptotic-preserving projective integration schemes for kinetic equations in the diffusion limit

We investigate a projective integration scheme for a kinetic equation in the limit of vanishing mean free path, in which the kinetic description approaches a diffusion phenomenon. The scheme first takes a few small steps with a simple, explicit method, such as a spatial centered flux/forward Euler time integration, and subsequently projects the results forward in time over a large time step on the diffusion time scale. We show that, with an appropriate choice of the inner step size, the time-step restriction on the outer time step is similar to the stability condition for the diffusion equation, whereas the required number of inner steps does not depend on the mean free path. We also provide a consistency result. The presented method is asymptotic-preserving, in the sense that the method converges to a standard finite volume scheme for the diffusion equation in the limit of vanishing mean free path. The analysis is illustrated with numerical results, and we present an application to the Su-Olson test.

0 citations0 reviewsNo code linkedOpen access