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

Alfons G. Hoekstra

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

ResearcherUnclaimed source record
Computational Engineering, Finance, and Sciencephysics.comp-phComputationDistributed, Parallel, and Cluster Computingphysics.opticsApplicationsBiological Physicscond-mat.softcond-mat.stat-mechInformation TheoryMachine Learningmath.ITnlin.PSPerformancephysics.data-anphysics.med-ph

Catalog footprint

What is connected

9works
16topics
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

9 published item(s)

preprint2022arXiv

Convergence of the discrete dipole approximation. I. Theoretical analysis

We performed a rigorous theoretical convergence analysis of the discrete dipole approximation (DDA). We prove that errors in any measured quantity are bounded by a sum of a linear and quadratic term in the size of a dipole d, when the latter is in the range of DDA applicability. Moreover, the linear term is significantly smaller for cubically than for non-cubically shaped scatterers. Therefore, for small d errors for cubically shaped particles are much smaller than for non-cubically shaped. The relative importance of the linear term decreases with increasing size, hence convergence of DDA for large enough scatterers is quadratic in the common range of d. Extensive numerical simulations were carried out for a wide range of d. Finally we discuss a number of new developments in DDA and their consequences for convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The discrete dipole approximation: an overview and recent developments

We present a review of the discrete dipole approximation (DDA), which is a general method to simulate light scattering by arbitrarily shaped particles. We put the method in historical context and discuss recent developments, taking the viewpoint of a general framework based on the integral equations for the electric field. We review both the theory of the DDA and its numerical aspects, the latter being of critical importance for any practical application of the method. Finally, the position of the DDA among other methods of light scattering simulation is shown and possible future developments are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Uncertainty quantification of a three-dimensional in-stent restenosis model with surrogate modelling

In-Stent Restenosis is a recurrence of coronary artery narrowing due to vascular injury caused by balloon dilation and stent placement. It may lead to the relapse of angina symptoms or to an acute coronary syndrome. An uncertainty quantification of a model for In-Stent Restenosis with four uncertain parameters (endothelium regeneration time, the threshold strain for smooth muscle cells bond breaking, blood flow velocity and the percentage of fenestration in the internal elastic lamina) is presented. Two quantities of interest were studied, namely the average cross-sectional area and the maximum relative area loss in a vessel. Due to the computational intensity of the model and the number of evaluations required in the uncertainty quantification, a surrogate model, based on Gaussian process regression with proper orthogonal decomposition, was developed which subsequently replaced the original In-Stent Restenosis model in the uncertainty quantification. A detailed analysis of the uncertainty propagation and sensitivity analysis is presented. Around 11% and 16% of uncertainty are observed on the average cross-sectional area and maximum relative area loss respectively, and the uncertainty estimates show that a higher fenestration mainly determines uncertainty in the neointimal growth at the initial stage of the process. On the other hand, the uncertainty in blood flow velocity and endothelium regeneration time mainly determine the uncertainty in the quantities of interest at the later, clinically relevant stages of the restenosis process. The uncertainty in the threshold strain is relatively small compared to the other uncertain parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Semi-intrusive uncertainty propagation for multiscale models

A family of semi-intrusive uncertainty propagation (UP) methods for multiscale models is introduced. The methods are semi-intrusive in the sense that inspection of the model is limited up to the level of the single scale systems, and viewing these single scale components as black-boxes. The goal is to estimate uncertainty in the result of multiscale models at a reduced amount of time as compared to black-box Monte Carlo (MC). In the resulting semi-intrusive MC method, the required number of samples of an expensive single scale model is minimized in order to reduce the execution time for the overall UP. In the metamodeling approach, the expensive model component is replaced completely by a computationally much cheaper surrogate model. These semi-intrusive algorithms have been tested on two case studies based on reaction-diffusion dynamics. The results demonstrate that the proposed semi-intrusive methods can result in a significant reduction of the computational time for multiscale UP, while still computing accurately the estimates of uncertainties. The semi-intrusive methods can, therefore, be a valid alternative, when uncertainties of a multiscale model cannot be estimated by the black-box MC methods in a feasible amount of time.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Shear Thickening of Dense Suspensions: The Role of Friction

Shear thickening of particle suspensions is characterized by a transition between lubricated and frictional contacts between the particles. Using 3D numerical simulations, we study how the inter-particle friction coefficient influences the effective macroscopic friction coefficient and hence the microstructure and rheology of dense shear thickening suspensions. We propose expressions for effective friction coefficient in terms of distance to jamming for varying shear stresses and particle friction coefficient values. We find effective friction coefficient to be rather insensitive to interparticle friction, which is perhaps surprising but agrees with recent theory and experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Information geometric analysis of phase transitions in complex patterns: the case of the Gray-Scott reaction-diffusion model

The Fisher-Rao metric from Information Geometry is related to phase transition phenomena in classical statistical mechanics. Several studies propose to extend the use of Information Geometry to study more general phase transitions in complex systems. However, it is unclear whether the Fisher-Rao metric does indeed detect these more general transitions, especially in the absence of a statistical model. In this paper we study the transitions between patterns in the Gray-Scott reaction-diffusion model using Fisher information. We describe the system by a probability density function that represents the size distribution of blobs in the patterns and compute its Fisher information with respect to changing the two rate parameters of the underlying model. We estimate the distribution non-parametrically so that we do not assume any statistical model. The resulting Fisher map can be interpreted as a phase-map of the different patterns. Lines with high Fisher information can be considered as boundaries between regions of parameter space where patterns with similar characteristics appear. These lines of high Fisher information can be interpreted as phase transitions between complex patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-parametric estimation of Fisher information from real data

The Fisher Information matrix is a widely used measure for applications ranging from statistical inference, information geometry, experiment design, to the study of criticality in biological systems. Yet there is no commonly accepted non-parametric algorithm to estimate it from real data. In this rapid communication we show how to accurately estimate the Fisher information in a nonparametric way. We also develop a numerical procedure to minimize the errors by choosing the interval of the finite difference scheme necessary to compute the derivatives in the definition of the Fisher information. Our method uses the recently published "Density Estimation using Field Theory" algorithm to compute the probability density functions for continuous densities. We use the Fisher information of the normal distribution to validate our method and as an example we compute the temperature component of the Fisher Information Matrix in the two dimensional Ising model and show that it obeys the expected relation to the heat capacity and therefore peaks at the phase transition at the correct critical temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Distributed Multiscale Computing with MUSCLE 2, the Multiscale Coupling Library and Environment

We present the Multiscale Coupling Library and Environment: MUSCLE 2. This multiscale component-based execution environment has a simple to use Java, C++, C, Python and Fortran API, compatible with MPI, OpenMP and threading codes. We demonstrate its local and distributed computing capabilities and compare its performance to MUSCLE 1, file copy, MPI, MPWide, and GridFTP. The local throughput of MPI is about two times higher, so very tightly coupled code should use MPI as a single submodel of MUSCLE 2; the distributed performance of GridFTP is lower, especially for small messages. We test the performance of a canal system model with MUSCLE 2, where it introduces an overhead as small as 5% compared to MPI.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Flexible composition and execution of high performance, high fidelity multiscale biomedical simulations

Multiscale simulations are essential in the biomedical domain to accurately model human physiology. We present a modular approach for designing, constructing and executing multiscale simulations on a wide range of resources, from desktops to petascale supercomputers, including combinations of these. Our work features two multiscale applications, in-stent restenosis and cerebrovascular bloodflow, which combine multiple existing single-scale applications to create a multiscale simulation. These applications can be efficiently coupled, deployed and executed on computers up to the largest (peta) scale, incurring a coupling overhead of 1 to 10% of the total execution time.

0 citations0 reviewsNo code linkedOpen access