Spectral Ewald Acceleration of Stokesian Dynamics for polydisperse suspensions

In this work we develop the Spectral Ewald Accelerated Stokesian Dynamics (SEASD), a novel computational method for dynamic simulations of polydisperse colloidal suspensions with full hydrodynamic interactions. SEASD is based on the framework of Stokesian Dynamics (SD) with extension to compressible solvents, and uses the Spectral Ewald (SE) method [Lindbo & Tornberg, J. Comput. Phys. 229 (2010) 8994] for the wave-space mobility computation. To meet the performance requirement of dynamic simulations, we use Graphic Processing Units (GPU) to evaluate the suspension mobility, and achieve an order of magnitude speedup compared to a CPU implementation. For further speedup, we develop a novel far-field block-diagonal preconditioner to reduce the far-field evaluations in the iterative solver, and SEASD-nf, a polydisperse extension of the mean-field Brownian approximation of Banchio & Brady [J. Chem. Phys. 118 (2003) 10323]. We extensively discuss implementation and parameter selection strategies in SEASD, and demonstrate the spectral accuracy in the mobility evaluation and the overall $\mathcal{O}(N\log N)$ computation scaling. We present three computational examples to further validate SEASD and SEASD-nf in monodisperse and bidisperse suspensions: the short-time transport properties, the equilibrium osmotic pressure and viscoelastic moduli, and the steady shear Brownian rheology. Our validation results show that the agreement between SEASD and SEASD-nf is satisfactory over a wide range of parameters, and also provide significant insight into the dynamics of polydisperse colloidal suspensions.Comment: 39 pages, 21 figure

Short-time diffusion in concentrated bidisperse hard-sphere suspensions

Diffusion in bidisperse Brownian hard-sphere suspensions is studied by Stokesian Dynamics (SD) computer simulations and a semi-analytical theoretical scheme for colloidal short-time dynamics, based on Beenakker and Mazur's method [Physica 120A, 388 (1983) & 126A, 349 (1984)]. Two species of hard spheres are suspended in an overdamped viscous solvent that mediates the salient hydrodynamic interactions among all particles. In a comprehensive parameter scan that covers various packing fractions and suspension compositions, we employ numerically accurate SD simulations to compute the initial diffusive relaxation of density modulations at the Brownian time scale, quantified by the partial hydrodynamic functions. A revised version of Beenakker and Mazur's $\delta\gamma$-scheme for monodisperse suspensions is found to exhibit surprisingly good accuracy, when simple rescaling laws are invoked in its application to mixtures. The so-modified $\delta\gamma$ scheme predicts hydrodynamic functions in very good agreement with our SD simulation results, for all densities from the very dilute limit up to packing fractions as high as $40\%$.Comment: 12 pages, 6 figure

Deductive Verification of Unmodified Linux Kernel Library Functions

This paper presents results from the development and evaluation of a deductive verification benchmark consisting of 26 unmodified Linux kernel library functions implementing conventional memory and string operations. The formal contract of the functions was extracted from their source code and was represented in the form of preconditions and postconditions. The correctness of 23 functions was completely proved using AstraVer toolset, although success for 11 functions was achieved using 2 new specification language constructs. Another 2 functions were proved after a minor modification of their source code, while the final one cannot be completely proved using the existing memory model. The benchmark can be used for the testing and evaluation of deductive verification tools and as a starting point for verifying other parts of the Linux kernel.Comment: 18 pages, 2 tables, 6 listings. Accepted to ISoLA 2018 conference. Evaluating Tools for Software Verification trac

Diquark in Nona-quark States

We study the nonaquark states $S^0(3115)$ and $S^+(3140)$ which are reported by KEK-PS (Phys.Lett. B597 (2004) 236; nucl-ex/0310018) by means of the quark model with diquark correlation. The nonaquark states form $\bf{1},\bf{8},\bf{10},\bar{\bf{10}},\bf{27},\bar{\bf{35}}$ SU(3) multiplets. The flavor wave functions of all the nonaquark states are constructed through the standard tensor technique. The mass spectrum is studied by using Gell-Mann-Okubo mass formula. Some nonaquark mass sum rules are obtained. We further investigate the decay of $S^0(3115)$ and $S^+(3140)$ under the assumption of "fall-apart" mechanism. It has been found that the main decay mode is $\Sigma NN$ rather than $\Lambda NN$ which is consistent with experiment. Also we have uniquely determine the flavor wave function of $S^0(3115)$ which belong to $\bf{27}$-plet with the quantum number $Y=2,I=1,I_z=-1$. Whereas the exotic states $S^+(3140)$ can belong to either $\bf{27}$-plet or $\bar{\bf{35}}$-plet. In the exact $SU(3)^{flavor}\times SU(3)^{color}\times SU(2)^{spin}$ limit, both $S^0(3115)$ and $S^+(3140)$ belong to ${\bf 27}$-plet with negative parity. We predict that its flavor structure can be determined by measuring the branch fractions of its decay channels. The experiments to check this prediction are expected.Comment: 23 pages, 8 figures, to appear in Phys. Rev.

Challenges of Primary Frequency Control and Benefits of Primary Frequency Response Support from Electric Vehicles

As the integration of wind generation displaces conventional plants, system inertia provided by rotating mass declines, causing concerns over system frequency stability. This paper implements an advanced stochastic scheduling model with inertia-dependent fast frequency response requirements to investigate the challenges on the primary frequency control in the future Great Britain electricity system. The results suggest that the required volume and the associated cost of primary frequency response increase significantly along with the increased capacity of wind plants. Alternative measures (e.g. electric vehicles) have been proposed to alleviate these concerns. Therefore, this paper also analyses the benefits of primary frequency response support from electric vehicles in reducing system operation cost, wind curtailment and carbon emissions

Relating the maximum radial stress on pile shaft to pile base resistance

An approximate analytic relationship is developed between the maximum radial stress on the shaft of a displacement pile in sand and the base resistance of the pile. Using the cavity expansion analogy, together with a confined failure mechanism, the ratio between the two quantities, defined as a factor S t, is established as a function of the friction angle, shear stiffness, compressibility and mean effective stress of the sand near the pile tip. It is shown that, given otherwise identical input parameters, the value of S t will decrease with increasing friction angle, and with decreasing mean stress level. It also tends to decrease with an increase in relative density. It is predicted that S t has typical values between 0.03 and 0.05, in broad agreement with the range of empirically derived values in the literature. The relationship also predicts that S t may take much higher values (~0.1) for piles installed in dense sand or in highly compressible sand. Because of the analytical nature, the established relationship provides useful insights into the mechanisms involved and important implications for design practice.published_or_final_versio

Oxidized (non)-regenerated cellulose affects fundamental cellular processes of wound healing.

In this study we investigated how hemostats such as oxidized regenerated cellulose (ORC, TABOTAMP) and oxidized non-regenerated cellulose (ONRC, RESORBA CELL) influence local cellular behavior and contraction of the extracellular matrix (ECM). Human stromal fibroblasts were inoculated in vitro with ORC and ONRC. Cell proliferation was assayed over time, and migration was evaluated by Live Cell imaging microscopy. Fibroblasts grown in collagen-gels were treated with ORC or ONRC, and ECM contraction was measured utilizing a contraction assay. An absolute pH decline was observed with both ORC and ONRC after 1 hour. Mean daily cell proliferation, migration and matrix contraction were more strongly inhibited by ONRC when compared with ORC (p < 0.05). When control media was pH-lowered to match the lower pH values typically seen with ORC and ONRC, significant differences in cell proliferation and migration were still observed between ONRC and ORC (p < 0.05). However, in these pH conditions, inhibition of matrix contraction was only significant for ONRC (p < 0.05). We find that ORC and ONRC inhibit fibroblast proliferation, migration and matrix contraction, and stronger inhibition of these essential cellular processes of wound healing were observed for ONRC when compared with ORC. These results will require further validation in future in vivo experiments to clarify the clinical implications for hemostat use in post-surgical wound healing

Constant Stress and Pressure Rheology of Colloidal Suspensions

We study the constant stress and pressure rheology of dense hard-sphere colloidal suspensions using Brownian dynamics simulation. Expressing the flow behavior in terms of the friction coefficient—the ratio of shear to normal stress—reveals a shear arrest point from the collapse of the rheological data in the non-Brownian limit. The flow curves agree quantitatively (when scaled) with the experiments of Boyer et al. [Phys. Rev. Lett. 107, 188301 (2011)]. Near suspension arrest, both the shear and the incremental normal viscosities display a universal power law divergence. This work shows the important role of jamming on the arrest of colloidal suspensions and illustrates the care needed when conducting and analyzing experiments and simulations near the flow-arrest transition

Short-time transport properties of bidisperse suspensions and porous media: a Stokesian Dynamics study

We present a comprehensive computational study of the short-time transport properties of bidisperse neutral colloidal suspensions and the corresponding porous media. Our study covers bidisperse particle size ratios up to $4$, and total volume fractions up to and beyond the monodisperse hard-sphere close packing limit. The many-body hydrodynamic interactions are computed using conventional Stokesian Dynamics (SD) via a Monte-Carlo approach. We address suspension properties including the short-time translational and rotational self-diffusivities, the instantaneous sedimentation velocity, the wavenumber-dependent partial hydrodynamic functions, and the high-frequency shear and bulk viscosities; and porous media properties including the permeability and the translational and rotational hindered diffusivities. We carefully compare the SD computations with existing theoretical and numerical results. For suspensions, we also explore the range of validity of various approximation schemes, notably the Pairwise Additive (PA) approximations with the Percus-Yevick structural input. We critically assess the strengths and weaknesses of the SD algorithm for various transport properties. For very dense systems, we discuss in detail the interplay between the hydrodynamic interactions and the structures due to the presence of a second species of a different size.Comment: 26 pages, 26 figure