Search

This work proposes stochastic partial differential equations (SPDEs) as a practical tool to replicate clustering effects of more detailed particle-based dynamics. Inspired by membrane-mediated receptor dynamics on cell surfaces, we formulate a stochastic particle-based model for diffusion and pairwise interaction of particles, leading to intriguing clustering phenomena. Employing numerical simulation and cluster detection methods, we explore the approximation of the particle-based clustering dynamics through mean-field approaches. We find that SPDEs successfully reproduce spatiotemporal clustering dynamics, not only in the initial cluster formation period, but also on longer time scales where the successive merging of clusters cannot be tracked by deterministic mean-field models. The computational efficiency of the SPDE approach allows us to generate extensive statistical data for parameter estimation in a simpler model that uses a Markov jump process to capture the temporal evolution of the cluster number.

Bitumen residues in the oil sand tailings can be a threat to the environment that separating them from tailings before disposal is crucial. However, low bitumen concentration in the tailings slurry and the complex transport characteristics of the four-phase mixture make the process difficult. This study establishes an Eulerian-Eulerian CFD model for an industrial-scale oil sand tailings pipeline. A comprehensive sensitivity analysis was conducted on the selection of carrier-solid and solid-bitumen drag models. The combination of small and large particle sizes (i.e., 75 & 700 um) and bitumen droplet size (i.e., 400 um) provided good agreement with field data in velocity profiles and pressure drop. The validated model was subsequently extended to investigate the influence of the secondary phase (i.e., bitumen droplets and bubbles) on flow characteristics in a tailing pipeline. The investigation covered a range of bitumen droplet size (100-400 um), bitumen fraction (0.0025-0.1), bubble size (5-1000 um), and bubble fraction (0.0025-0.3) and their influences on the velocity, solids, and bitumen distribution are revealed. For an optimum bubble size of 500 um, a maximum recovery of 59% from the top 50 % and 83 % from the top 75 % of the pipe cross-section was obtained. The present study demonstrates the preferential distribution of bitumen and provides valuable insight on bitumen recovery from an industrial-scale tailings pipeline., Comment: 51, 14 figures

The dynamics of molecules are governed by rare event transitions between long-lived (metastable) states. To explore these transitions efficiently, many enhanced sampling protocols have been introduced that involve using simulations with biases or changed temperatures. Two established statistically optimal estimators for obtaining unbiased equilibrium properties from such simulations are the multistate Bennett Acceptance Ratio (MBAR) and the transition-based reweighting analysis method (TRAM). Both MBAR and TRAM are solved iteratively and can suffer from long convergence times. Here we introduce stochastic approximators (SA) for both estimators, resulting in SAMBAR and SATRAM, which are shown to converge faster than their deterministic counterparts, without significant accuracy loss. Both methods are demonstrated on different molecular systems.

Understanding the flow behavior of complex concentrated slurries is of tremendous importance for industrial waste management. In this study, the transport of three-phase oil sands tailings in a horizontal pipeline is simulated via the mixture multiphase model coupled with the kinetic theory of granular flow. The solid particles and bitumen droplets are conveyed via a non-Newtonian carrier fluid in a turbulent regime inside an industrial-scale pipeline. The simulation results showed exceptional agreement with the field data, with errors of <3.5% for velocity distribution and <15% for the pressure drop. A systematic parametric investigation was performed for a wide range of flow conditions, showing that the majority of bitumen droplets reside at the top region of the pipe. Our findings may help design an effective process for the separation of bitumen residues during pipeline transport., Comment: 49 pages, 20 figures

Asphaltene precipitation is a crucial phase separation phenomenon in the oil industry, especially in paraffinic froth treatment to extract bitumen from oil sands ores. This work reveals the formation of particles at 0.2 to 0.4 um in radius, defined as a primary sub-micron particle (PSMP), which is ubiquitous from diffusive mixing between asphaltene solution and any of 23 types of precipitants examined in our experiments. The yield and quantity of asphaltene particles are affected not only by the Hildebrand solubility parameter of the precipitants but also by the diffusion coefficients of the asphaltene solution and the precipitant. The Population Balance Model (PBM) with the Hildebrand solubility parameter has been used to model particle size distribution. Good agreement has been achieved between numerical predictions and the experimental data. It indicates that the colloid theory can describe the size distribution of PSMP and larger aggregates. Therefore, this study provides new insight into the mechanism for the dependence of yield and size distribution of the precipitated asphaltene particles on the composition of precipitants and adding inhibitors. Controlling of asphaltene yield and size distribution may be applied to the process of hydrocarbon separation or asphaltene precipitation prevention., Comment: Preprint submitted to Fuel

The great challenge with biological membrane systems is the wide range of scales involved, from nanometers and picoseconds for individual lipids, to the micrometers and beyond millisecond for cellular signalling processes. While solvent-free coarse-grained membrane models are convenient for large-scale simulations, and promising to provide insight into slow processes involving membranes, these models usually have unrealistic kinetics. One major obstacle is the lack of an equally convenient way of introducing hydrodynamic coupling without significantly increasing the computational cost of the model. To address this, we introduce a framework based on anisotropic Langevin dynamics, for which major in-plane and out-of-plane hydrodynamic effects are modeled via friction and diffusion tensors from analytical or semi-analytical solutions to Stokes hydrodynamic equations. Using this framework, we obtain accurate dispersion relations for planar membrane patches, both free-standing and in the vicinity of a wall. We also briefly discuss how non-equilibrium dynamics is affected by hydrodynamic interactions., Comment: Submitted to The Journal of Chemical Physics

We present a simple and computationally efficient coarse-grained and solvent-free model for simulating lipid bilayer membranes. In order to be used in concert with particle-based reaction-diffusion simulations, the model is purely based on interacting and reacting particles, each representing a coarse patch of a lipid monolayer. Particle interactions include nearest-neighbor bond-stretching and angle-bending, and are parameterized so as to reproduce the local membrane mechanics given by the Helfrich energy density over a range of relevant curvatures. In-plane fluidity is implemented with Monte Carlo bond-flipping moves. The physical accuracy of the model is verified by five tests: (i) Power spectrum analysis of equilibrium thermal undulations is used to verify that the particle-based representation correctly captures the dynamics predicted by the continuum model of fluid membranes. (ii) It is verified that the input bending stiffness, against which the potential parameters are optimized, is accurately recovered. (iii) Isothermal area compressibility modulus of the membrane is calculated and is shown to be tunable to reproduce available values for different lipid bilayers, independent of the bending rigidity. (iv) Simulation of two-dimensional shear flow under a gravity force is employed to measure the effective in-plane viscosity of the membrane model, and show the possibility of modeling membranes with specified viscosities. (v) Interaction of the bilayer membrane with a spherical nanoparticle is modeled as a test case for large membrane deformations and budding involved in cellular processes such as endocytosis..., Comment: 13 pages, 8 figures, 3 tables, manuscript accepted for publication at The Journal of Chemical Physics