1
|
Sgouros AP, Theodorou DN. Development of a Meshless Kernel-Based Scheme for Particle-Field Brownian Dynamics Simulations. J Phys Chem B 2024; 128:6907-6921. [PMID: 38984836 PMCID: PMC11264276 DOI: 10.1021/acs.jpcb.4c01441] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Revised: 05/29/2024] [Accepted: 06/21/2024] [Indexed: 07/11/2024]
Abstract
We develop a meshless discretization scheme for particle-field Brownian dynamics simulations. The density is assigned on the particle level using a weighting kernel with finite support. The system's free energy density is derived from an equation of state (EoS) and includes a square gradient term. The numerical stability of the scheme is evaluated in terms of reproducing the thermodynamics (equilibrium density and compressibility) and dynamics (diffusion coefficient) of homogeneous samples. Using a reduced description to simplify our analysis, we find that numerical stability depends strictly on reduced reference compressibility, kernel range, time step in relation to the friction factor, and reduced external pressure, the latter being relevant under isobaric conditions. Appropriate parametrization yields precise thermodynamics, further improved through a simple renormalization protocol. The dynamics can be restored exactly through a trivial manipulation of the time step and friction coefficient. A semiempirical formula for the upper bound on the time step is derived, which takes into account variations in compressibility, friction factor, and kernel range. We test the scheme on realistic mesoscopic models of fluids, involving both simple (Helfand) and more sophisticated (Sanchez-Lacombe) equations of state.
Collapse
Affiliation(s)
- Aristotelis P. Sgouros
- School of Chemical Engineering, National Technical University of Athens (NTUA), GR-15780 Athens, Greece
| | - Doros N. Theodorou
- School of Chemical Engineering, National Technical University of Athens (NTUA), GR-15780 Athens, Greece
| |
Collapse
|
2
|
Silva DPF, Coelho RCV, Pagonabarraga I, Succi S, Telo da Gama MM, Araújo NAM. Lattice Boltzmann simulation of deformable fluid-filled bodies: progress and perspectives. SOFT MATTER 2024; 20:2419-2441. [PMID: 38420837 PMCID: PMC10933750 DOI: 10.1039/d3sm01648j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/06/2023] [Accepted: 02/20/2024] [Indexed: 03/02/2024]
Abstract
With the rapid development of studies involving droplet microfluidics, drug delivery, cell detection, and microparticle synthesis, among others, many scientists have invested significant efforts to model the flow of these fluid-filled bodies. Motivated by the intricate coupling between hydrodynamics and the interactions of fluid-filled bodies, several methods have been developed. The objective of this review is to present a compact foundation of the methods used in the literature in the context of lattice Boltzmann methods. For hydrodynamics, we focus on the lattice Boltzmann method due to its specific ability to treat time- and spatial-dependent boundary conditions and to incorporate new physical models in a computationally efficient way. We split the existing methods into two groups with regard to the interfacial boundary: fluid-structure and fluid-fluid methods. The fluid-structure methods are characterised by the coupling between fluid dynamics and mechanics of the flowing body, often used in applications involving membranes and similar flexible solid boundaries. We further divide fluid-structure-based methods into two subcategories, those which treat the fluid-structure boundary as a continuum medium and those that treat it as a discrete collection of individual springs and particles. Next, we discuss the fluid-fluid methods, particularly useful for the simulations of fluid-fluid interfaces. We focus on models for immiscible droplets and their interaction in a suspending fluid and describe benchmark tests to validate the models for fluid-filled bodies.
Collapse
Affiliation(s)
- Danilo P F Silva
- Centro de Física Teórica e Computacional, Faculdade de Ciências, Universidade de Lisboa, P-1749-016 Lisboa, Portugal.
- Departamento de Física, Faculdade de Ciências, Universidade de Lisboa, P-1749-016 Lisboa, Portugal
| | - Rodrigo C V Coelho
- Centro de Física Teórica e Computacional, Faculdade de Ciências, Universidade de Lisboa, P-1749-016 Lisboa, Portugal.
- Departamento de Física, Faculdade de Ciências, Universidade de Lisboa, P-1749-016 Lisboa, Portugal
| | - Ignacio Pagonabarraga
- Departament de Física de la Matèria Condensada, Universitat de Barcelona, Carrer de Martí Franqués 1, 08028 Barcelona, Spain
- Universitat de Barcelona Institute of Complex Systems (UBICS), Universitat de Barcelona, 08028 Barcelona, Spain
| | - Sauro Succi
- Center for Life Nano Science at La Sapienza, Istituto Italiano di Tecnologia, 295 Viale Regina Elena, I/00161 Roma, Italy
- Harvard Institute for Applied Computational Science, Cambridge, MA 02138, USA
| | - Margarida M Telo da Gama
- Centro de Física Teórica e Computacional, Faculdade de Ciências, Universidade de Lisboa, P-1749-016 Lisboa, Portugal.
- Departamento de Física, Faculdade de Ciências, Universidade de Lisboa, P-1749-016 Lisboa, Portugal
| | - Nuno A M Araújo
- Centro de Física Teórica e Computacional, Faculdade de Ciências, Universidade de Lisboa, P-1749-016 Lisboa, Portugal.
- Departamento de Física, Faculdade de Ciências, Universidade de Lisboa, P-1749-016 Lisboa, Portugal
| |
Collapse
|
3
|
Jäger T, Mokos A, Prasianakis NI, Leyer S. Validating the Transition Criteria from the Cassie-Baxter to the Wenzel State for Periodically Pillared Surfaces with Lattice Boltzmann Simulations. ACS OMEGA 2024; 9:10592-10601. [PMID: 38463292 PMCID: PMC10918652 DOI: 10.1021/acsomega.3c08862] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Revised: 02/08/2024] [Accepted: 02/12/2024] [Indexed: 03/12/2024]
Abstract
Microfabrication techniques allow the development and production of artificial superhydrophobic surfaces that possess a precisely controlled roughness at the micrometer level, typically achieved through the arrangement of micropillar structures in periodic patterns. In this work, we analyze the stability and energy barrier of droplets in the Cassie-Baxter (CB) state on such periodic patterns. In addition, we further develop a transition criterion using the CB equation and derive an improved version which allows predicting for which pillar geometries, equilibrium contact angles, and droplet volumes the CB state switches from a metastable to an unstable state. This enables a comparison with existing experiments and three-dimensional multiphase Lattice Boltzmann simulations for different pillar distances, two contact angles, and two droplet volumes, where a good agreement has been found.
Collapse
Affiliation(s)
- Tobias Jäger
- Department
of Engineering, Faculty of Science, Technology and Medicine, University of Luxembourg, Luxembourg L-1359, Luxembourg
| | - Athanasios Mokos
- Transport
Mechanisms Group, Laboratory for Waste Management, Paul Scherrer Institute, PSI, Villigen 5232, Switzerland
| | - Nikolaos I. Prasianakis
- Transport
Mechanisms Group, Laboratory for Waste Management, Paul Scherrer Institute, PSI, Villigen 5232, Switzerland
| | - Stephan Leyer
- Department
of Engineering, Faculty of Science, Technology and Medicine, University of Luxembourg, Luxembourg L-1359, Luxembourg
| |
Collapse
|
4
|
Gogoi D, Chauhan A, Puri S, Singh A. Segregation of fluids with polymer additives at domain interfaces: a dissipative particle dynamics study. SOFT MATTER 2023; 19:6433-6445. [PMID: 37403605 DOI: 10.1039/d3sm00504f] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/06/2023]
Abstract
This paper investigates the phase separation kinetics of ternary fluid mixtures composed of a polymeric component (C) and two simple fluids (A and B) using dissipative particle dynamics simulations with a system dimensionality of d = 3. We model the affinities between the components to enable the settling of the polymeric component at the interface of fluids A and B. Thus, the system evolves to form polymer coated morphologies, enabling alteration of the fluids' interfacial properties. This manipulation can be utilized across various disciplines, such as the stabilization of emulsions and foams, rheological control, biomimetic design, and surface modification. We probe the effects of various parameters, such as the polymeric concentration, chain stiffness, and length, on the phase separation kinetics of the system. The simulation results show that changes in the concentration of flexible polymers exhibit perfect dynamic scaling for coated morphologies. The growth rate decreases as the polymeric composition is increased due to reduced surface tension and restricted connectivity between A- and B-rich clusters. Variations in the polymer chain rigidity at fixed composition ratios and degrees of polymerization slow the evolution kinetics of AB fluids marginally, although the effect is more pronounced for perfectly rigid chains. Whereas flexible polymer chain lengths at fixed composition ratios slow down the segregation kinetics of AB fluids slightly, varying the chain lengths of perfectly rigid polymers leads to a significant deviation in the length scale and dynamic scaling for the evolved coated morphologies. The characteristic length scale follows a power-law growth with a growth exponent ϕ that shows a crossover from the viscous to the inertial hydrodynamic regime, where the values of ϕ depend on the constraints imposed on the system.
Collapse
Affiliation(s)
- Dorothy Gogoi
- School of Physical Sciences, Jawaharlal Nehru University, New Delhi-110067, India
| | - Avinash Chauhan
- Department of Physics, Indian Institute of Technology (BHU), Varanasi, Uttar Pradesh-221005, India.
| | - Sanjay Puri
- School of Physical Sciences, Jawaharlal Nehru University, New Delhi-110067, India
| | - Awaneesh Singh
- Department of Physics, Indian Institute of Technology (BHU), Varanasi, Uttar Pradesh-221005, India.
| |
Collapse
|
5
|
Qian JX, Wang J, Lu YQ. Motional consensus of self-propelled particles. Sci Rep 2023; 13:8169. [PMID: 37210399 DOI: 10.1038/s41598-023-35238-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2022] [Accepted: 05/15/2023] [Indexed: 05/22/2023] Open
Abstract
The motional consensus of self-propelled particles is studied in both noise-free cases and cases with noise by the standard Vicsek model. In the absence of noise, we propose a simple method, using grid-based technique and defining the normalized variance of the ratio of the number of particles locally to globally, to quantitatively study the movement pattern of the system by the spatial distribution of the particles and the degree of aggregation of particles. It is found that the weaker correlation of velocity leads to larger degree of aggregation of the particles. In the cases with noise, we quantify the competition between velocity alignment and noise by considering the difference of the variety of order parameter result from the velocity alignment and noise. The variation of the effect of noise on motional consensus is non-monotonic for the change of the probability distribution of noise from uniform to non-uniform. Our results may be useful and encourage further efforts in exploring the basic principles of collective motion.
Collapse
Affiliation(s)
- Jia-Xin Qian
- National Laboratory of Solid State Microstructures, Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing, 210093, China
- College of Engineering and Applied Sciences, Nanjing University, Nanjing, 210093, China
| | - Jun Wang
- National Laboratory of Solid State Microstructures, Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing, 210093, China.
- School of Physics, Nanjing University, Nanjing, 210093, China.
| | - Yan-Qing Lu
- National Laboratory of Solid State Microstructures, Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing, 210093, China.
- College of Engineering and Applied Sciences, Nanjing University, Nanjing, 210093, China.
| |
Collapse
|
6
|
Akalın AA, Dedekargınoğlu B, Choi SR, Han B, Ozcelikkale A. Predictive Design and Analysis of Drug Transport by Multiscale Computational Models Under Uncertainty. Pharm Res 2023; 40:501-523. [PMID: 35650448 PMCID: PMC9712595 DOI: 10.1007/s11095-022-03298-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2021] [Accepted: 05/17/2022] [Indexed: 01/18/2023]
Abstract
Computational modeling of drug delivery is becoming an indispensable tool for advancing drug development pipeline, particularly in nanomedicine where a rational design strategy is ultimately sought. While numerous in silico models have been developed that can accurately describe nanoparticle interactions with the bioenvironment within prescribed length and time scales, predictive design of these drug carriers, dosages and treatment schemes will require advanced models that can simulate transport processes across multiple length and time scales from genomic to population levels. In order to address this problem, multiscale modeling efforts that integrate existing discrete and continuum modeling strategies have recently emerged. These multiscale approaches provide a promising direction for bottom-up in silico pipelines of drug design for delivery. However, there are remaining challenges in terms of model parametrization and validation in the presence of variability, introduced by multiple levels of heterogeneities in disease state. Parametrization based on physiologically relevant in vitro data from microphysiological systems as well as widespread adoption of uncertainty quantification and sensitivity analysis will help address these challenges.
Collapse
Affiliation(s)
- Ali Aykut Akalın
- Department of Mechanical Engineering, Middle East Technical University, 06531, Ankara, Turkey
| | - Barış Dedekargınoğlu
- Department of Mechanical Engineering, Middle East Technical University, 06531, Ankara, Turkey
| | - Sae Rome Choi
- School of Mechanical Engineering, Purdue University, 585 Purdue Mall, West Lafayette, Indiana, 47907, USA
| | - Bumsoo Han
- School of Mechanical Engineering, Purdue University, 585 Purdue Mall, West Lafayette, Indiana, 47907, USA.
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, Indiana, USA.
- Center for Cancer Research, Purdue University, 585 Purdue Mall, West Lafayette, Indiana, 47907, USA.
| | - Altug Ozcelikkale
- Department of Mechanical Engineering, Middle East Technical University, 06531, Ankara, Turkey.
| |
Collapse
|
7
|
Hülsberg M, Klapp SHL. Depinning dynamics of confined colloidal dispersions under oscillatory shear. Phys Rev E 2023; 107:014603. [PMID: 36797876 DOI: 10.1103/physreve.107.014603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Accepted: 12/06/2022] [Indexed: 06/18/2023]
Abstract
Strongly confined colloidal dispersions under shear can exhibit a variety of dynamical phenomena, including depinning transitions and complex structural changes. Here, we investigate the behavior of such systems under pure oscillatory shearing with shear rate γ[over ̇](t)=γ[over ̇]_{0}cos(ωt), as it is a common scenario in rheological experiments. The colloids' depinning behavior is assessed from a particle level based on trajectories, obtained from overdamped Brownian dynamics simulations. The numerical approach is complemented by an analytic one based on an effective single-particle model in the limits of weak and strong driving. Investigating a broad spectrum of shear rate amplitudes γ[over ̇]_{0} and frequencies ω, we observe complete pinning as well as temporary depinning behavior. We discover that temporary depinning occurs for shear rate amplitudes above a frequency-dependent critical amplitude γ[over ̇]_{0}^{crit}(ω), for which we attain an approximate functional expression. For a range of frequencies, approaching γ[over ̇]_{0}^{crit}(ω) is accompanied by a strongly increasing settling time. Above γ[over ̇]_{0}^{crit}(ω), we further observe a variety of dynamical structures, whose stability exhibits an intriguing (γ[over ̇]_{0},ω) dependence. This might enable new perspectives for potential control schemes.
Collapse
Affiliation(s)
- Marcel Hülsberg
- Institut für Theoretische Physik, Hardenbergstrasse 36, Technische Universität Berlin, D-10623 Berlin, Germany
| | - Sabine H L Klapp
- Institut für Theoretische Physik, Hardenbergstrasse 36, Technische Universität Berlin, D-10623 Berlin, Germany
| |
Collapse
|
8
|
Ding M, Li L. Flow-Induced Translocation and Conformational Transition of Polymer Chains through Nanochannels: Recent Advances and Future Perspectives. Macromolecules 2021. [DOI: 10.1021/acs.macromol.1c00909] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Mingming Ding
- School of Chemical Engineering and Light Industry, Guangdong University of Technology, Guangzhou 510006, China
| | - Lianwei Li
- Food Science and Processing Research Center, College of Chemistry and Environmental Engineering, Shenzhen University, Shenzhen 518060, China
| |
Collapse
|
9
|
Smirnov MA, Tolmachev DA, Glova AD, Sokolova MP, Geydt PV, Lukasheva NV, Lyulin SV. Combined Use of Atomic Force Microscopy and Molecular Dynamics in the Study of Biopolymer Systems. POLYMER SCIENCE SERIES C 2021. [DOI: 10.1134/s1811238221020089] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
|
10
|
Floyd C, Chandresekaran A, Ni H, Ni Q, Papoian GA. Segmental Lennard-Jones interactions for semi-flexible polymer networks. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1910358] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Affiliation(s)
- Carlos Floyd
- Biophysics Program, University of Maryland, College Park, MD, USA
| | | | - Haoran Ni
- Biophysics Program, University of Maryland, College Park, MD, USA
| | - Qin Ni
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD, USA
| | - Garegin A. Papoian
- Biophysics Program, University of Maryland, College Park, MD, USA
- Department of Chemistry and Biochemistry, University of Maryland, College Park, MD, USA
- Institute for Physical Science and Technology, University of Maryland, College Park, MD, USA
| |
Collapse
|
11
|
Driven Transport of Dilute Polymer Solutions through Porous Media Comprising Interconnected Cavities. COLLOIDS AND INTERFACES 2021. [DOI: 10.3390/colloids5020022] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Driven transport of dilute polymer solutions through porous media has been simulated using a recently proposed novel dissipative particle dynamics method satisfying the no-penetration and no-slip boundary conditions. The porous media is an array of overlapping spherical cavities arranged in a simple cubic lattice. Simulations were performed for linear, ring, and star polymers with 12 arms for two cases with the external force acting on (I) both polymer and solvent beads to model a pressure-driven flow; (II) polymer beads only, similar to electrophoresis. When the external force is in the direction of a principal axis, the extent of change in the polymers’ conformation and their alignment with the driving force is more significant for case I. These effects are most pronounced for linear chains, followed by rings and stars at the same molecular weight. Moreover, the polymer mean velocity is affected by its molecular weight and architecture as well as the direction and strength of the imposed force.
Collapse
|
12
|
Liu Y, Widmer-Cooper A. A dissipative particle dynamics model for studying dynamic phenomena in colloidal rod suspensions. J Chem Phys 2021; 154:104120. [PMID: 33722052 DOI: 10.1063/5.0041285] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
A dissipative particle dynamics (DPD) model is developed and demonstrated for studying dynamics in colloidal rod suspensions. The solvent is modeled as conventional DPD particles, while individual rods are represented by a rigid linear chain consisting of overlapping solid spheres, which interact with solvent particles through a hard repulsive potential. The boundary condition on the rod surface is controlled using a surface friction between the solid spheres and the solvent particles. In this work, this model is employed to study the diffusion of a single colloid in the DPD solvent and compared with theoretical predictions. Both the translational and rotational diffusion coefficients obtained at a proper surface friction show good agreement with calculations based on the rod size defined by the hard repulsive potential. In addition, the system-size dependence of the diffusion coefficients shows that the Navier-Stokes hydrodynamic interactions are correctly included in this DPD model. Comparing our results with experimental measurements of the diffusion coefficients of gold nanorods, we discuss the ability of the model to correctly describe dynamics in real nanorod suspensions. Our results provide a clear reference point from which the model could be extended to enable the study of colloid dynamics in more complex situations or for other types of particles.
Collapse
Affiliation(s)
- Yawei Liu
- ARC Centre of Excellence in Exciton Science, School of Chemistry, University of Sydney, Sydney, New South Wales 2006, Australia
| | - Asaph Widmer-Cooper
- ARC Centre of Excellence in Exciton Science, School of Chemistry, University of Sydney, Sydney, New South Wales 2006, Australia
| |
Collapse
|
13
|
Montessori A, Tiribocchi A, Bonaccorso F, Lauricella M, Succi S. Lattice Boltzmann simulations capture the multiscale physics of soft flowing crystals. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2020; 378:20190406. [PMID: 32564719 PMCID: PMC7333952 DOI: 10.1098/rsta.2019.0406] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 02/26/2020] [Indexed: 06/11/2023]
Abstract
The study of the underlying physics of soft flowing materials depends heavily on numerical simulations, due to the complex structure of the governing equations reflecting the competition of concurrent mechanisms acting at widely disparate scales in space and time. A full-scale computational modelling remains a formidable challenge since it amounts to simultaneously handling six or more spatial decades in space and twice as many in time. Coarse-grained methods often provide a viable strategy to significantly mitigate this issue, through the implementation of mesoscale supramolecular forces designed to capture the essential physics at a fraction of the computational cost of a full-detail description. Here, we review some recent advances in the design of a lattice Boltzmann mesoscale approach for soft flowing materials, inclusive of near-contact interactions (NCIs) between dynamic interfaces, as they occur in high packing-fraction soft flowing crystals. The method proves capable of capturing several aspects of the rheology of soft flowing crystals, namely, (i) a 3/2 power-law dependence of the dispersed phase flow rate on the applied pressure gradient, (ii) the structural transition between an ex-two and ex-one (bamboo) configurations with the associated drop of the flow rate, (iii) the onset of interfacial waves once NCI is sufficiently intense. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.
Collapse
Affiliation(s)
- A. Montessori
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, Italy
| | - A. Tiribocchi
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, Italy
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Rome, Italy
| | - F. Bonaccorso
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Rome, Italy
| | - M. Lauricella
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, Italy
| | - S. Succi
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, Italy
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Rome, Italy
- Institute for Applied Computational Science, John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA
| |
Collapse
|
14
|
From CS, Sauret E, Galindo-Torres SA, Gu YT. Application of high-order lattice Boltzmann pseudopotential models. Phys Rev E 2020; 101:033303. [PMID: 32290007 DOI: 10.1103/physreve.101.033303] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2019] [Accepted: 02/10/2020] [Indexed: 11/07/2022]
Abstract
Higher-order lattice Boltzmann (LB) pseudopotential models have great potential for solving complex fluid dynamics in various areas of modern science. The discreteness of the lattice discretization makes these models an attractive choice due to their flexibility, capacity to capture hydrodynamic details, and inherent adaptability to parallel computations. Despite those advantages, the discreteness makes high-order LB models difficult to apply due to the larger lattice structure, for which basic fundamental properties, namely diffusion coefficient and contact angle, remain unknown. This work addresses this by providing general continuum solutions for those two basic properties and demonstrating these solutions to compare favorably against known theory. Various high-order LB models are shown to reproduce the sinusoidal decay of a binary miscible mixture accurately and consistently. Furthermore, these models are shown to reproduce neutral, hydrophobic, and hydrophilic contact angles. Discrete differences are shown to exist, which are captured at the discrete level and confirmed through droplet shape analysis. This work provides practical tools that allow for high-order LB pseudopotential models to be used to simulate multicomponent flows.
Collapse
Affiliation(s)
- C S From
- School of Mechanical, Medical and Process Engineering, Science and Engineering Faculty, Queensland University of Technology, QLD 4001, Australia
| | - E Sauret
- School of Mechanical, Medical and Process Engineering, Science and Engineering Faculty, Queensland University of Technology, QLD 4001, Australia
| | - S A Galindo-Torres
- School of Engineering, Westlake University, Hangzhou Zhejiang Province 310024, China.,Institute of Advanced Technology, Westlake Institute for Advanced Study, Hangzhou Zhejiang Province 310024, China
| | - Y T Gu
- School of Mechanical, Medical and Process Engineering, Science and Engineering Faculty, Queensland University of Technology, QLD 4001, Australia
| |
Collapse
|
15
|
Soleymani FA, Ripoll M, Gompper G, Fedosov DA. Dissipative particle dynamics with energy conservation: Isoenergetic integration and transport properties. J Chem Phys 2020; 152:064112. [PMID: 32061230 DOI: 10.1063/1.5119778] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Simulations of nano- to micro-meter scale fluidic systems under thermal gradients require consistent mesoscopic methods accounting for both hydrodynamic interactions and proper transport of energy. One such method is dissipative particle dynamics with energy conservation (DPDE), which has been used for various fluid systems with non-uniform temperature distributions. We propose an easily parallelizable modification of the velocity-Verlet algorithm based on local energy redistribution for each DPDE particle such that the total energy in a simulated system is conserved up to machine precision. Furthermore, transport properties of a DPDE fluid are analyzed in detail. In particular, an analytical approximation for the thermal conductivity coefficient is derived, which allows its a priori estimation for a given parameter set. Finally, we provide approximate expressions for the dimensionless Prandtl and Schmidt numbers, which characterize fluid transport properties and can be adjusted independently by a proper selection of model parameters. In conclusion, our results strengthen the DPDE method as a very robust approach for the investigation of mesoscopic systems with temperature inhomogeneities.
Collapse
Affiliation(s)
- Fatemeh A Soleymani
- Theoretical Soft Matter and Biophysics, Institute of Complex Systems and Institute for Advanced Simulation, Forschungszentrum Jülich, 52425 Jülich, Germany
| | - Marisol Ripoll
- Theoretical Soft Matter and Biophysics, Institute of Complex Systems and Institute for Advanced Simulation, Forschungszentrum Jülich, 52425 Jülich, Germany
| | - Gerhard Gompper
- Theoretical Soft Matter and Biophysics, Institute of Complex Systems and Institute for Advanced Simulation, Forschungszentrum Jülich, 52425 Jülich, Germany
| | - Dmitry A Fedosov
- Theoretical Soft Matter and Biophysics, Institute of Complex Systems and Institute for Advanced Simulation, Forschungszentrum Jülich, 52425 Jülich, Germany
| |
Collapse
|
16
|
|
17
|
Unraveling the Vascular Fate of Deformable Circulating Tumor Cells Via a Hierarchical Computational Model. Cell Mol Bioeng 2019. [DOI: 10.1007/s12195-019-00587-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
|
18
|
Hafner AE, Krausser J, Šarić A. Minimal coarse-grained models for molecular self-organisation in biology. Curr Opin Struct Biol 2019; 58:43-52. [PMID: 31226513 DOI: 10.1016/j.sbi.2019.05.018] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2019] [Revised: 05/13/2019] [Accepted: 05/19/2019] [Indexed: 01/19/2023]
Abstract
The molecular machinery of life is largely created via self-organisation of individual molecules into functional assemblies. Minimal coarse-grained models, in which a whole macromolecule is represented by a small number of particles, can be of great value in identifying the main driving forces behind self-organisation in cell biology. Such models can incorporate data from both molecular and continuum scales, and their results can be directly compared to experiments. Here we review the state of the art of models for studying the formation and biological function of macromolecular assemblies in living organisms. We outline the key ingredients of each model and their main findings. We illustrate the contribution of this class of simulations to identifying the physical mechanisms behind life and diseases, and discuss their future developments.
Collapse
Affiliation(s)
- Anne E Hafner
- Department of Physics and Astronomy, Institute for the Physics of Living Systems, University College London, London WC1E 6BT, UK
| | - Johannes Krausser
- Department of Physics and Astronomy, Institute for the Physics of Living Systems, University College London, London WC1E 6BT, UK
| | - Anđela Šarić
- Department of Physics and Astronomy, Institute for the Physics of Living Systems, University College London, London WC1E 6BT, UK.
| |
Collapse
|
19
|
Belardinelli D, Sbragaglia M, Benzi R, Ciliberto S. Lattice Boltzmann simulations of nonequilibrium fluctuations in a nonideal binary mixture. Phys Rev E 2019; 99:063302. [PMID: 31330737 DOI: 10.1103/physreve.99.063302] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2018] [Indexed: 06/10/2023]
Abstract
In recent years the lattice Boltzmann (LB) methodology has been fruitfully extended to include the effects of thermal fluctuations. So far, all studied cases pertain to equilibrium fluctuations, i.e., fluctuations with respect to an equilibrium background state. In this paper we take a step further and present results of fluctuating LB simulations of a binary mixture confined between two parallel walls in the presence of a constant concentration gradient in the wall-to-wall direction. This is a paradigmatic setup for the study of nonequilibrium (NE) fluctuations, i.e., fluctuations with respect to a nonequilibrium state. We analyze the dependence of the structure factors for the hydrodynamical fields on the wave vector q in both the directions parallel and perpendicular to the walls, highlighting the long-range (∼|q|^{-4}) nature of correlations in the NE framework. Results at the small scales (high wave numbers) quantitatively agree with the predictions of fluctuating hydrodynamics without fitting parameters. At larger scales (low wave numbers), however, results show finite-size effects induced by confinement and call for further studies aimed at controlling boundary conditions in the fluctuating LB framework as well as compressibility effects. Moreover, in the presence of a nonideal equation of state of the mixture, we also observe that the (spatially homogeneous) average pressure changes, due to a genuinely new contribution triggered by NE fluctuations. These NE pressure effects are studied at changing the system size and the concentration gradient. Taken all together, we argue that the results of this article are useful and instrumental to boost the applicability of the fluctuating LB methodology in the framework of NE fluctuations, possibly in conjunction with experiments.
Collapse
Affiliation(s)
- Daniele Belardinelli
- Department of Physics & INFN, University of Rome "Tor Vergata", Via della Ricerca Scientifica 1, 00133, Rome, Italy
| | - Mauro Sbragaglia
- Department of Physics & INFN, University of Rome "Tor Vergata", Via della Ricerca Scientifica 1, 00133, Rome, Italy
| | - Roberto Benzi
- Department of Physics & INFN, University of Rome "Tor Vergata", Via della Ricerca Scientifica 1, 00133, Rome, Italy
| | - Sergio Ciliberto
- Laboratoire de Physique de Ecole Normale Superieure de Lyon (CNRS UMR5672), 46 Allée d'Italie, 69364, Lyon, France
| |
Collapse
|
20
|
Campeggio J, Polimeno A, Zerbetto M. DiTe2: Calculating the diffusion tensor for flexible molecules. J Comput Chem 2018; 40:697-705. [DOI: 10.1002/jcc.25742] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2018] [Revised: 10/06/2018] [Accepted: 10/09/2018] [Indexed: 01/29/2023]
Affiliation(s)
| | - Antonino Polimeno
- Department of Chemical SciencesUniversity of Padua Padova 35131 Italy
| | - Mirco Zerbetto
- Department of Chemical SciencesUniversity of Padua Padova 35131 Italy
| |
Collapse
|
21
|
The sickle cell trait affects contact dynamics and endothelial cell activation in Plasmodium falciparum-infected erythrocytes. Commun Biol 2018; 1:211. [PMID: 30534603 PMCID: PMC6269544 DOI: 10.1038/s42003-018-0223-3] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2018] [Accepted: 11/06/2018] [Indexed: 11/08/2022] Open
Abstract
Sickle cell trait, a common hereditary blood disorder, protects carriers from severe disease in infections with the human malaria parasite Plasmodium falciparum. Protection is associated with a reduced capacity of parasitized erythrocytes to cytoadhere to the microvascular endothelium and cause vaso-occlusive events. However, the underpinning cellular and biomechanical processes are only partly understood and the impact on endothelial cell activation is unclear. Here, we show, by combining quantitative flow chamber experiments with multiscale computer simulations of deformable cells in hydrodynamic flow, that parasitized erythrocytes containing the sickle cell haemoglobin displayed altered adhesion dynamics, resulting in restricted contact footprints on the endothelium. Main determinants were cell shape, knob density and membrane bending. As a consequence, the extent of endothelial cell activation was decreased. Our findings provide a quantitative understanding of how the sickle cell trait affects the dynamic cytoadhesion behavior of parasitized erythrocytes and, in turn, endothelial cell activation.
Collapse
|