Particle-Based Monte-Carlo Simulations of Steady-State Mass Transport at Intermediate Péclet Numbers

Abstract

Conventional approaches for simulating steady-state distributions of particles under diffusive and advective transport at high Péclet numbers involve solving the diffusion and advection equations in at least two dimensions. Here, we present an alternative computational strategy by combining a particle-based rather than a field-based approach with the initialisation of particles in proportion to their flux. This method allows accurate prediction of the steady state and is applicable even at high Péclet numbers where traditional particle-based Monte-Carlo methods starting from randomly initialised particle distributions fail. We demonstrate that generating a flux of particles according to a predetermined density and velocity distribution at a single fixed time and initial location allows for accurate simulation of mass transport under flow. Specifically, upon initialisation in proportion to their flux, these particles are propagated individually and detected by summing up their Monte-Carlo trajectories in predefined detection regions. We demonstrate quantitative agreement of the predicted concentration profiles with the results of experiments performed with fluorescent particles in microfluidic channels under continuous flow. This approach is computationally advantageous and readily allows non-trivial initial distributions to be considered. In particular, this method is highly suitable for simulating advective and diffusive transport in microfluidic devices.

Secondary nucleation and elongation occur at different sites on Alzheimer’s amyloid-β aggregates

“…The fabrication and the operation of the microfluidic diffusion device used in the present studies have been described in previous papers ( 33 , 48 ). Briefly, the microfluidic chips were fabricated by using standard soft lithography.…”

β-Synuclein suppresses both the initiation and amplification steps of α-synuclein aggregation via competitive binding to surfaces

“…Microfluidic diffusion devices were fabricated according to previously published protocols3746. Devices were cast using polydimethylsiloxane (PDMS) from a silicon wafer master imprinted with 25  μ m high channels.…”

Cooperative Assembly of Hsp70 Subdomain Clusters

“… 39 The cysteine residues of SBD641 were conjugated to an AlexaFluor488 dye for visualization of the sample diffusion profiles, and images were taken at 12 different positions under steady-state flow conditions using the experimental procedure shown in Figure 1 b. The diffusion profiles acquired in this way were then fitted to a linear combination of simulated basis functions calculated for ideal monodisperse solutions 44 ( Figure 1 c), yielding the average diffusion coefficients of the particles in the solution. 37 , 39 , 49 …”

“… 39 The basis functions used in this analysis were simulated in the range of radii relevant to the samples, namely, 0.25–10 nm in 0.25 nm increments. 44 …”