BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Lennaert van Veen

Lennaert van Veen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.DSphysics.flu-dynBiological Physicsmath.NANeurons and CognitionDistributed, Parallel, and Cluster Computingnlin.PSphysics.comp-phphysics.med-phTissues and Organs

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2020arXiv

A brief history of simple invariant solutions in turbulence

When studying fluid mechanics in terms of instability, bifurcation and invariant solutions one quickly finds out how little can be done by pen and paper. For flows on sufficiently simple domains and under sufficiently simple boundary conditions, one may be able to predict the parameter values at which the base flow becomes unstable and the basic properties of the secondary flow. On more complicated domains and under more realistic boundary conditions, such questions can usually only be addressed by numerical means. Moreover, for a wide class of elementary parallel shear flows the base flow remains stable in the presence of sustained turbulent motion. In such flows, secondary solutions often appear with finite amplitude and completely unconnected to the base flow. Only using techniques from computational dynamical systems can such behaviour be explained. Many of these techniques, such as for the detection and classification of bifurcations and for the continuation in parameters of equilibria and time-periodic solutions, were developed in the late 1970s for dynamical systems with few degrees of freedom. The application to fluid dynamics or, to be more precise, to spatially discretized Navier-Stokes flow, is far from straightforward. In this historical review chapter, we follow the development of this field of research from the valiant naivety of the early 1980s to the open challenges of today.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

I'm a doctor, not a mathematician! Homeostasis as a proportional-integral control system

The distinction between "healthy" and "unhealthy" patients is commonly based on single, discrete values taken at an isolated point in time (e.g., blood pressure or core temperature). Perhaps a more robust and insightful diagnosis can be obtained by studying the functional interdependence of such indicators and the homeostasis that controls them. This requires quasi-continuous measurements and a procedure to map the data onto a parsimonious control model with a degree of universality. The current research illustrates this approach using glucose homeostasis as a target. Data were obtained from 41 healthy subjects wearing over-the-counter glucose monitors, and projected onto a simple proportional-integral (PI) controller, widely used in engineering applications. The indicators quantifying the control function are clustered for the great majority of subjects, while a few outliers exhibit less efficient homeostasis. Practical implications for healthcare and education are further discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Green's function-based time stepping for the Kuramoto-Sivashinsky initial-boundary value problem

Both theoretical and numerical studies of the Kuramoto-Sivashinsky equation have mostly considered periodic boundary conditions. In this setting, the Fourier decomposition of the solution is central to theoretical ideas, such as renormalization group arguments, as well as to numerical solution, allowing for the construction of accurate and efficient time-steppers using standard pseudo-spectral methods. In contrast, fixed boundary conditions induce boundary layers and necessitate the use of non-uniform grids, usually generated by orthogonal polynomials. On such bases, numerical differentiation is ill-conditioned and can potentially lead to a catastrophic blow-up of round-off error. In this paper, we use ideas recently explored by Viswanath (J. Comput. Phys. 251(2013), pp. 414-431) to completely eliminate numerical differentiation and linear solving from the time-stepping algorithm. We use the Green's function-based method to investigate elements of the Kuramoto-Sivashinsky dynamics over a range of five decades of the viscosity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Periodic solutions to a mean-field model for electrocortical activity

We consider a continuum model of electrical signals in the human cortex, which takes the form of a system of semilinear, hyperbolic partial differential equations for the inhibitory and excitatory membrane potentials and the synaptic inputs. The coupling of these components is represented by sigmoidal and quadratic nonlinearities. We consider these equations on a square domain with periodic boundary conditions, in the vicinity of the primary transition from a stable equilibrium to time-periodic motion through an equivariant Hopf bifurcation. We compute part of a family of standing wave solutions, emanating from this point.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

PAMPAC: A Parallel Adaptive Method for Pseudo-Arclength Continuation

Pseudo-arclength continuation is a well-established method for generating a numerical curve approximating the solution of an underdetermined system of nonlinear equations. It is an inherently sequential predictor-corrector method in which new approximate solutions are extrapolated from previously converged results and then iteratively refined. Convergence of the iterative corrections is guaranteed only for sufficiently small prediction steps. In high-dimensional systems, corrector steps are extremely costly to compute and the prediction step-length must be adapted carefully to avoid failed steps or unnecessarily slow progress. We describe a parallel method for adapting the step-length employing several predictor-corrector sequences of different step lengths computed concurrently. In addition, the algorithm permits intermediate results of unconverged correction sequences to seed new predictions. This strategy results in an aggressive optimization of the step length at the cost of redundancy in the concurrent computation. We present two examples of convoluted solution curves of high-dimensional systems showing that speed-up by a factor of two can be attained on a multi-core CPU while a factor of three is attainable on a small cluster.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A homoclinic tangle on the edge of shear turbulence

Experiments and simulations lend mounting evidence for the edge state hypothesis on subcritical transition to turbulence, which asserts that simple states of fluid motion mediate between laminar and turbulent shear flow as their stable manifolds separate the two in state space. In this Letter we describe a flow homoclinic to a time-periodic edge state. Its existence explains turbulent bursting through the classical Smale-Birkhoff theorem. During a burst, vortical structures and the associated energy dissipation are highly localized near the wall, in contrast to the familiar regeneration cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Significance of Simple Invariant Solutions in Turbulent Flows

Recent remarkable progress in computing power and numerical analysis is enabling us to fill a gap in the dynamical systems approach to turbulence. One of the significant advances in this respect has been the numerical discovery of simple invariant sets, such as nonlinear equilibria and periodic solutions, in well-resolved Navier--Stokes flows. This review describes some fundamental and practical aspects of dynamical systems theory for the investigation of turbulence, focusing on recently found invariant solutions and their significance for the dynamical and statistical characterization of low-Reynolds-number turbulent flows. It is shown that the near-wall regeneration cycle of coherent structures can be reproduced by such solutions. The typical similarity laws of turbulence, i.e. the Prandtl wall law and the Kolmogorov law for the viscous range, as well as the pattern and intensity of turbulence-driven secondary flow in a square duct can also be represented by these simple invariant solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Metabifurcation analysis of a mean field model of the cortex

Mean field models (MFMs) of cortical tissue incorporate salient features of neural masses to model activity at the population level. One of the common aspects of MFM descriptions is the presence of a high dimensional parameter space capturing neurobiological attributes relevant to brain dynamics. We study the physiological parameter space of a MFM of electrocortical activity and discover robust correlations between physiological attributes of the model cortex and its dynamical features. These correlations are revealed by the study of bifurcation plots, which show that the model responses to changes in inhibition belong to two families. After investigating and characterizing these, we discuss their essential differences in terms of four important aspects: power responses with respect to the modeled action of anesthetics, reaction to exogenous stimuli, distribution of model parameters and oscillatory repertoires when inhibition is enhanced. Furthermore, while the complexity of sustained periodic orbits differs significantly between families, we are able to show how metamorphoses between the families can be brought about by exogenous stimuli. We unveil links between measurable physiological attributes of the brain and dynamical patterns that are not accessible by linear methods. They emerge when the parameter space is partitioned according to bifurcation responses. This partitioning cannot be achieved by the investigation of only a small number of parameter sets, but is the result of an automated bifurcation analysis of a representative sample of 73,454 physiologically admissible sets. Our approach generalizes straightforwardly and is well suited to probing the dynamics of other models with large and complex parameter spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On matrix-free computation of 2D unstable manifolds

Recently, a flexible and stable algorithm was introduced for the computation of 2D unstable manifolds of periodic solutions to systems of ordinary differential equations. The main idea of this approach is to represent orbits in this manifold as the solutions of an appropriate boundary value problem. The boundary value problem is under determined and a one parameter family of solutions can be found by means of arclength continuation. This family of orbits covers a piece of the manifold. The quality of this covering depends on the way the boundary value problem is discretised, as do the tractability and accuracy of the computation. In this paper, we describe an implementation of the orbit continuation algorithm which relies on multiple shooting and Newton-Krylov continuation. We show that the number of time integrations necessary for each continuation step scales only with the number of shooting intervals but not with the number of degrees of freedom of the dynamical system. The number of shooting intervals is chosen based on linear stability analysis to keep the conditioning of the boundary value problem in check. We demonstrate our algorithm with two test systems: a low-order model of shear flow and a well-resolved simulation of turbulent plane Couette flow.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The quasi-periodic doubling cascade in the transition to weak turbulence

The quasi-periodic doubling cascade is shown to occur in the transition from regular to weakly turbulent behaviour in simulations of incompressible Navier-Stokes flow on a three-periodic domain. Special symmetries are imposed on the flow field in order to reduce the computational effort. Thus we can apply tools from dynamical systems theory such as continuation of periodic orbits and computation of Lyapunov exponents. We propose a model ODE for the quasi-period doubling cascade which, in a limit of a perturbation parameter to zero, avoids resonance related problems. The cascade we observe in the simulations is then compared to the perturbed case, in which resonances complicate the bifurcation scenario. In particular, we compare the frequency spectrum and the Lyapunov exponents. The perturbed model ODE is shown to be in good agreement with the simulations of weak turbulence. The scaling of the observed cascade is shown to resemble the unperturbed case, which is directly related to the well known doubling cascade of periodic orbits.

0 citations0 reviewsNo code linkedOpen access