Do you want to publish a course? Click here

Multiscale method based on coupled lattice-Boltzmann and Langevin-dynamics for direct simulation of nanoscale particle/polymer suspensions in complex flows

81   0   0.0 ( 0 )
 Added by Zixiang Liu
 Publication date 2018
  fields Physics
and research's language is English




Ask ChatGPT about the research

A hybrid computational method coupling the lattice-Boltzmann (LB) method and a Langevin-dynamics (LD) method is developed to simulate nanoscale particle and polymer (NPP) suspensions in the presence of both thermal fluctuation and long-range many-body hydrodynamic interactions (HI). Brownian motion of the NPP is explicitly captured by a stochastic forcing term in the LD method. The LD method is two-way coupled to the non-fluctuating LB fluid through a discrete LB forcing source distribution to capture the long-range HI. To ensure intrinsically linear scalability with respect to the number of particles, an Eulerian-host algorithm for short-distance particle neighbor search and interaction is developed and embedded to LB-LD framework. The validity and accuracy of the LB-LD approach are demonstrated through several sample problems. The simulation results show good agreements with theory and experiment. The LB-LD approach can be favorably incorporated into complex multiscale computational frameworks for efficiently simulating multiscale, multicomponent particulate suspension systems such as complex blood suspensions.



rate research

Read More

Most biological fluids are viscoelastic, meaning that they have elastic properties in addition to the dissipative properties found in Newtonian fluids. Computational models can help us understand viscoelastic flow, but are often limited in how they deal with complex flow geometries and suspended particles. Here, we present a lattice Boltzmann solver for Oldroyd-B fluids that can handle arbitrarily-shaped fixed and moving boundary conditions, which makes it ideally suited for the simulation of confined colloidal suspensions. We validate our method using several standard rheological setups, and additionally study a single sedimenting colloid, also finding good agreement with literature. Our approach can readily be extended to constitutive equations other than Oldroyd-B. This flexibility and the handling of complex boundaries holds promise for the study of microswimmers in viscoelastic fluids.
A cascaded lattice Boltzmann (LB) approach based on central moments and multiple relaxation times to simulate thermal convective flows, which are driven by buoyancy forces and/or swirling effects, in the cylindrical coordinate system with axial symmetry is presented. In this regard, the dynamics of the axial and radial momentum components along with the pressure are represented by means of the 2D Navier-Stokes equations with geometric mass and momentum source terms in the pseudo Cartesian form, while the evolutions of the azimuthal momentum and the temperature field are each modeled by an advection-diffusion type equation with appropriate local source terms. Based on these, cascaded LB schemes involving three distribution functions are formulated to solve for the fluid motion in the meridian plane using a D2Q9 lattice, and to solve for the azimuthal momentum and the temperature field each using a D2Q5 lattice. The geometric mass and momentum source terms for the flow fields and the energy source term for the temperature field are included using a new symmetric operator splitting technique, via pre-collision and post-collision source steps around the cascaded collision step for each distribution function. These result in a particularly simple and compact formulation to directly represent the effect of various geometric source terms consistently in terms of changes in the appropriate zeroth and first order moments. Simulations of several complex buoyancy-driven thermal flows and including rotational effects in cylindrical geometries using the new axisymmetric cascaded LB schemes show good agreement with prior benchmark results for the structures of the velocity and thermal fields as well as the heat transfer rates given in terms of the Nusselt numbers.
Janus phoretic colloids (JPs) self-propel as a result of self-generated chemical gradients and exhibit spontaneous nontrivial dynamics within phoretic suspensions, on length scales much larger than the microscopic swimmer size. Such collective dynamics arise from the competition of (i) the self-propulsion velocity of the particles, (ii) the attractive/repulsive chemically-mediated interactions between particles and (iii) the flow disturbance they introduce in the surrounding medium. These three ingredients are directly determined by the shape and physico-chemical properties of the colloids surface. Owing to such link, we adapt a recent and popular kinetic model for dilute suspensions of chemically-active JPs where the particles far-field hydrodynamic and chemical signatures are intrinsically linked and explicitly determined by the design properties. Using linear stability analysis, we show that self-propulsion can induce a wave-selective mechanism for certain particles configurations consistent with experimental observations. Numerical simulations of the complete kinetic model are further performed to analyze the relative importance of chemical and hydrodynamic interactions in the nonlinear dynamics. Our results show that regular patterns in the particle density are promoted by chemical signaling but prevented by the strong fluid flows generated collectively by the polarized particles, regardless of their chemotactic or antichemotactic nature (i.e. for both puller and pusher swimmers).
A series of benchmarks based on the physical situation of phase inversion between two immiscible liquids is presented. These benchmarks aim at progressing toward the direct numerical simulation of two-phase flows. Several CFD codes developed in French laboratories and using either Volume-of-Fluid or Level-Set interface tracking methods are used to provide physical solutions of the benchmarks, convergence studies and code comparisons. Two typical configurations are retained, with integral scale Reynolds numbers of 13.700 and 433.000, respectively. The physics of the problem are probed through macroscopic quantities such as potential and kinetic energies, or enstrophy. In addition, scaling laws for the temporal decay of the kinetic energy are derived to check the physical relevance of the simulations. Finally the droplet size distribution is probed. Additional test problems are also reported to estimate the influence of viscous effects in the vicinity of the interface.
The squirmer is a simple yet instructive model for microswimmers, which employs an effective slip velocity on the surface of a spherical swimmer to describe its self-propulsion. We solve the hydrodynamic flow problem with the lattice Boltzmann (LB) method, which is well-suited for time-dependent problems involving complex boundary conditions. Incorporating the squirmer into LB is relatively straight-forward, but requires an unexpectedly fine grid resolution to capture the physical flow fields and behaviors accurately. We demonstrate this using four basic hydrodynamic tests: Two for the far-field flow---accuracy of the hydrodynamic moments and squirmer-squirmer interactions---and two that require the near field to be accurately resolved---a squirmer confined to a tube and one scattering off a spherical obstacle---which LB is capable of doing down to the grid resolution. We find good agreement with (numerical) results obtained using other hydrodynamic solvers in the same geometries and identify a minimum required resolution to achieve this reproduction. We discuss our algorithm in the context of other hydrodynamic solvers and present an outlook on its application to multi-squirmer problems.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا