1
|
Ambruş VE, Bazzanini L, Gabbana A, Simeoni D, Succi S, Tripiccione R. Fast kinetic simulator for relativistic matter. Nat Comput Sci 2022; 2:641-654. [PMID: 38177272 DOI: 10.1038/s43588-022-00333-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/14/2022] [Accepted: 09/08/2022] [Indexed: 01/06/2024]
Abstract
Relativistic kinetic theory is ubiquitous to several fields of modern physics, finding application at large scales in systems in astrophysical contexts, all of the way down to subnuclear scales and into the realm of quark-gluon plasmas. This motivates the quest for powerful and efficient computational methods that are able to accurately study fluid dynamics in the relativistic regime as well as the transition to beyond hydrodynamics-in principle all of the way down to ballistic regimes. We present a family of relativistic lattice kinetic schemes for the efficient simulation of relativistic flows in both strongly (fluid) and weakly (rarefied gas) interacting regimes. The method can deal with both massless and massive particles, thereby encompassing ultra- and mildly relativistic regimes alike. The computational performance of the method for the simulation of relativistic flows across the aforementioned regimes is discussed in detail, along with prospects of future applications.
Collapse
Affiliation(s)
- V E Ambruş
- Institut für Theoretische Physik, Johann Wolfgang Goethe-Universität, Frankfurt, Germany
- Department of Physics, West University of Timişoara, Timisoara, Romania
| | - L Bazzanini
- Università di Ferrara and INFN-Ferrara, Ferrara, Italy
| | - A Gabbana
- Eindhoven University of Technology, Eindhoven, Netherlands
| | - D Simeoni
- Università di Ferrara and INFN-Ferrara, Ferrara, Italy.
- Bergische Universität Wuppertal, Wuppertal, Germany.
- University of Cyprus, Physics department, Nicosia, Cyprus.
| | - S Succi
- Center for Life, Nano & Neuro Science @ La Sapienza, Italian Institute of Technology, Roma, Italy
- Department of Physics, Harvard University, Cambridge, MA, USA
| | - R Tripiccione
- Università di Ferrara and INFN-Ferrara, Ferrara, Italy
| |
Collapse
|
2
|
Tiribocchi A, Montessori A, Durve M, Bonaccorso F, Lauricella M, Succi S. Dynamics of polydisperse multiple emulsions in microfluidic channels. Phys Rev E 2021; 104:065112. [PMID: 35030928 DOI: 10.1103/physreve.104.065112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2021] [Accepted: 12/09/2021] [Indexed: 06/14/2023]
Abstract
Multiple emulsions are a class of soft fluid in which small drops are immersed within a larger one and stabilized over long periods of time by a surfactant. We recently showed that, if a monodisperse multiple emulsion is subject to a pressure-driven flow, a wide variety of nonequilibrium steady states emerges at late times, whose dynamics relies on a complex interplay between hydrodynamic interactions and multibody collisions among internal drops. In this work, we use lattice Boltzmann simulations to study the dynamics of polydisperse double emulsions driven by a Poiseuille flow within a microfluidic channel. Our results show that their behavior is critically affected by multiple factors, such as initial position, polydispersity index, and area fraction occupied within the emulsion. While at low area fraction inner drops may exhibit either a periodic rotational motion (at low polydispersity) or arrange into nonmotile configurations (at high polydispersity) located far from each other, at larger values of area fraction they remain in tight contact and move unidirectionally. This decisively conditions their close-range dynamics, quantitatively assessed through a time-efficiency-like factor. Simulations also unveil the key role played by the capsule, whose shape changes can favor the formation of a selected number of nonequilibrium states in which both motile and nonmotile configurations are found.
Collapse
Affiliation(s)
- A Tiribocchi
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185 Rome, Italy
| | - A Montessori
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185 Rome, Italy
| | - M Durve
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy
| | - F Bonaccorso
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185 Rome, Italy
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy
- Department of Physics and INFN, University of Rome "Tor Vergata," Via della Ricerca Scientifica, 00133 Rome, Italy
| | - M Lauricella
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185 Rome, Italy
| | - S Succi
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185 Rome, Italy
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy
- Institute for Applied Computational Science, John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA
| |
Collapse
|
3
|
Tiribocchi A, Montessori A, Lauricella M, Bonaccorso F, Brown KA, Succi S. Microscale modelling of dielectrophoresis assembly processes. Philos Trans A Math Phys Eng Sci 2021; 379:20200407. [PMID: 34455845 DOI: 10.1098/rsta.2020.0407] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/04/2021] [Indexed: 06/13/2023]
Abstract
This work presents a microscale approach for simulating the dielectrophoresis assembly of polarizable particles under an external electric field. The model is shown to capture interesting dynamical and topological features, such as the formation of chains of particles and their incipient aggregation into hierarchical structures. A quantitative characterization in terms of the number and size of these structures is also discussed. This computational model could represent a viable numerical tool to study the mechanical properties of particle-based hierarchical materials and suggest new strategies for enhancing their design and manufacture. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.
Collapse
Affiliation(s)
- A Tiribocchi
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini, 19, 00185 Rome, Italy
| | - A Montessori
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini, 19, 00185 Rome, Italy
| | - M Lauricella
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini, 19, 00185 Rome, Italy
| | - F Bonaccorso
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini, 19, 00185 Rome, Italy
- Department of Physics and INFN, University of Rome 'Tor Vergata', Via della Ricerca Scientifica, 1 00133 Rome, Italy
| | - K A Brown
- Department of Mechanical Engineering, Boston University, Boston, MA 02215, USA
| | - S Succi
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini, 19, 00185 Rome, Italy
- Institute for Applied Computational Science, John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA
| |
Collapse
|
4
|
Tiribocchi A, Montessori A, Lauricella M, Bonaccorso F, Succi S, Aime S, Milani M, Weitz DA. The vortex-driven dynamics of droplets within droplets. Nat Commun 2021; 12:82. [PMID: 33398018 PMCID: PMC7782531 DOI: 10.1038/s41467-020-20364-0] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2020] [Accepted: 11/20/2020] [Indexed: 01/29/2023] Open
Abstract
Understanding the fluid-structure interaction is crucial for an optimal design and manufacturing of soft mesoscale materials. Multi-core emulsions are a class of soft fluids assembled from cluster configurations of deformable oil-water double droplets (cores), often employed as building-blocks for the realisation of devices of interest in bio-technology, such as drug-delivery, tissue engineering and regenerative medicine. Here, we study the physics of multi-core emulsions flowing in microfluidic channels and report numerical evidence of a surprisingly rich variety of driven non-equilibrium states (NES), whose formation is caused by a dipolar fluid vortex triggered by the sheared structure of the flow carrier within the microchannel. The observed dynamic regimes range from long-lived NES at low core-area fraction, characterised by a planetary-like motion of the internal drops, to short-lived ones at high core-area fraction, in which a pre-chaotic motion results from multi-body collisions of inner drops, as combined with self-consistent hydrodynamic interactions. The onset of pre-chaotic behavior is marked by transitions of the cores from one vortex to another, a process that we interpret as manifestations of the system to maximize its entropy by filling voids, as they arise dynamically within the capsule.
Collapse
Affiliation(s)
- A. Tiribocchi
- grid.25786.3e0000 0004 1764 2907Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, Roma, 00161 Italy ,grid.5326.20000 0001 1940 4177Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, 00185 Italy
| | - A. Montessori
- grid.5326.20000 0001 1940 4177Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, 00185 Italy
| | - M. Lauricella
- grid.5326.20000 0001 1940 4177Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, 00185 Italy
| | - F. Bonaccorso
- grid.25786.3e0000 0004 1764 2907Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, Roma, 00161 Italy ,grid.5326.20000 0001 1940 4177Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, 00185 Italy
| | - S. Succi
- grid.25786.3e0000 0004 1764 2907Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, Roma, 00161 Italy ,grid.5326.20000 0001 1940 4177Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, 00185 Italy ,grid.38142.3c000000041936754XInstitute for Applied Computational Science, John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138 USA
| | - S. Aime
- grid.38142.3c000000041936754XInstitute for Applied Computational Science, John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138 USA ,grid.15736.360000 0001 1882 0021Matiére Molle et Chimie, Ecole Supérieure de Physique et Chimie Industrielles, Paris, 75005 France
| | - M. Milani
- grid.4708.b0000 0004 1757 2822Universitá degli Studi di Milano, via Celoria 16, Milano, 20133 Italy
| | - D. A. Weitz
- grid.38142.3c000000041936754XInstitute for Applied Computational Science, John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138 USA ,grid.38142.3c000000041936754XDepartment of Physics, Harvard University, Cambridge, MA 02138 USA
| |
Collapse
|
5
|
Abstract
We present a mechanistic model of drug release from a multiple emulsion into an external surrounding fluid. We consider a single multilayer droplet where the drug kinetics are described by a pure diffusive process through different liquid shells. The multilayer problem is described by a system of diffusion equations coupled via interlayer conditions imposing continuity of drug concentration and flux. Mass resistance is imposed at the outer boundary through the application of a surfactant at the external surface of the droplet. The two-dimensional problem is solved numerically by finite volume discretization. Concentration profiles and drug release curves are presented for three typical round-shaped (circle, ellipse, and bullet) droplets and the dependency of the solution on the mass transfer coefficient at the surface analyzed. The main result shows a reduced release time for an increased elongation of the droplets.
Collapse
Affiliation(s)
- G Pontrelli
- Istituto per le Applicazioni del Calcolo, CNR, Via dei Taurini 19, 00185 Rome, Italy
| | - E J Carr
- School of Mathematical Sciences, Queensland University of Technology (QUT), Brisbane, Australia
| | - A Tiribocchi
- Istituto per le Applicazioni del Calcolo, CNR, Via dei Taurini 19, 00185 Rome, Italy
- Italian Institute of Technology, CNLS@Sapienza, Rome, Italy
| | - S Succi
- Istituto per le Applicazioni del Calcolo, CNR, Via dei Taurini 19, 00185 Rome, Italy
- Italian Institute of Technology, CNLS@Sapienza, Rome, Italy
| |
Collapse
|
6
|
Gabbana A, Simeoni D, Succi S, Tripiccione R. Probing bulk viscosity in relativistic flows. Philos Trans A Math Phys Eng Sci 2020; 378:20190409. [PMID: 32564720 PMCID: PMC7333951 DOI: 10.1098/rsta.2019.0409] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 02/26/2020] [Indexed: 06/11/2023]
Abstract
We derive an analytical connection between kinetic relaxation rate and bulk viscosity of a relativistic fluid in d spatial dimensions, all the way from the ultra-relativistic down to the near non-relativistic regime. Our derivation is based on both Chapman-Enskog asymptotic expansion and Grad's method of moments. We validate our theoretical results against a benchmark flow, providing further evidence of the correctness of the Chapman-Enskog approach; we define the range of validity of this approach and provide evidence of mounting departures at increasing Knudsen number. Finally, we present numerical simulations of transport processes in quark-gluon plasmas, with special focus on the effects of bulk viscosity which might prove amenable to future experimental verification. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.
Collapse
Affiliation(s)
- A. Gabbana
- Università di Ferrara and INFN-Ferrara, 44122 Ferrara, Italy
- Bergische Universität Wuppertal, 42119 Wuppertal, Germany
| | - D. Simeoni
- Università di Ferrara and INFN-Ferrara, 44122 Ferrara, Italy
- Bergische Universität Wuppertal, 42119 Wuppertal, Germany
- University of Cyprus, 1678 Nicosia, Cyprus
| | - S. Succi
- Center for Life Nano Science @ La Sapienza, Italian Institute of Technology, Viale Regina Elena 295, 00161 Roma, Italy
- Istituto Applicazioni del Calcolo, National Research Council of Italy, Via dei Taurini 19, 00185 Roma, Italy
| | - R. Tripiccione
- Università di Ferrara and INFN-Ferrara, 44122 Ferrara, Italy
| |
Collapse
|
7
|
Montessori A, Tiribocchi A, Bonaccorso F, Lauricella M, Succi S. Lattice Boltzmann simulations capture the multiscale physics of soft flowing crystals. Philos Trans A Math Phys Eng Sci 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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
8
|
Tiribocchi A, Montessori A, Miliani S, Lauricella M, La Rocca M, Succi S. Microvorticity fluctuations affect the structure of thin fluid films. Phys Rev E 2019; 100:042606. [PMID: 31770937 DOI: 10.1103/physreve.100.042606] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2019] [Indexed: 11/07/2022]
Abstract
The dynamic interaction of complex fluid interfaces is highly sensitive to near-contact interactions occurring at the scale of ten of nanometers. Such interactions are difficult to analyze because they couple self-consistently to the dynamic morphology of the evolving interface, as well as to the hydrodynamics of the interstitial fluid film. In this work, we show that, above a given magnitude threshold, near-contact interactions trigger nontrivial microvorticity patterns, which in turn affect the effective near-contact interactions, giving rise to persistent fluctuating ripples at the fluid interface. In such a regime, near-contact interactions may significantly affect the macroscopic arrangement of emulsion configurations, such as those arising in soft-flowing microfluidic crystals.
Collapse
Affiliation(s)
- A Tiribocchi
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy.,Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, Italy
| | - A Montessori
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, Italy
| | - S Miliani
- Department of Engineering, Roma Tre University, Via Vito Volterra 62, Rome 00146, Italy
| | - M Lauricella
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, Italy
| | - M La Rocca
- Department of Engineering, Roma Tre University, Via Vito Volterra 62, Rome 00146, Italy
| | - S Succi
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy.,Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, Rome, Italy.,Institute for Applied Computational Science, John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 01451, USA
| |
Collapse
|
9
|
Gabbana A, Simeoni D, Succi S, Tripiccione R. Relativistic dissipation obeys Chapman-Enskog asymptotics: Analytical and numerical evidence as a basis for accurate kinetic simulations. Phys Rev E 2019; 99:052126. [PMID: 31212459 DOI: 10.1103/physreve.99.052126] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2019] [Indexed: 11/07/2022]
Abstract
We present an analytical derivation of the transport coefficients of a relativistic gas in (2+1) dimensions for both Chapman-Enskog (CE) asymptotics and Grad's expansion methods. We further develop a systematic calibration method, connecting the relaxation time of relativistic kinetic theory to the transport parameters of the associated dissipative hydrodynamic equations. Comparison of our analytical results and numerical simulations shows that the CE method correctly captures dissipative effects, while Grad's method does not, in agreement with previous analyses performed in the (3+1)-dimensional case. These results provide a solid basis for accurately calibrated computational studies of relativistic dissipative flows.
Collapse
Affiliation(s)
- A Gabbana
- Department of Physics and Earth Sciences, Università di Ferrara and INFN-Ferrara, I-44122 Ferrara, Italy.,Chair of Applied Mathematics and Numerical Analysis, Bergische Universität Wuppertal, D-42119 Wuppertal, Germany
| | - D Simeoni
- Department of Physics and Earth Sciences, Università di Ferrara and INFN-Ferrara, I-44122 Ferrara, Italy.,Chair of Applied Mathematics and Numerical Analysis, Bergische Universität Wuppertal, D-42119 Wuppertal, Germany.,Department of Physics, University of Cyprus, CY-1678 Nicosia, Cyprus
| | - S Succi
- Center for Life Nano Science @ La Sapienza, Italian Institute of Technology, Viale Regina Elena 295, I-00161 Roma, Italy.,Istituto Applicazioni del Calcolo, National Research Council of Italy, Via dei Taurini 19, I-00185 Roma, Italy
| | - R Tripiccione
- Department of Physics and Earth Sciences, Università di Ferrara and INFN-Ferrara, I-44122 Ferrara, Italy
| |
Collapse
|
10
|
Abstract
We outline the main ideas behind the numerical modelling of soft flowing crystals, paying special attention to their application to microfluidic devices for the design of novel mesoscale porous materials. This article is part of the theme issue 'Multiscale modelling, simulation and computing: from the desktop to the exascale'.
Collapse
Affiliation(s)
- A. Montessori
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185 Rome, Italy
- Department of Engineering, University of Rome ‘Roma Tre’, Via della Vasca Navale 79, 00141 Rome, Italy
| | - M. Lauricella
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185 Rome, Italy
| | - S. Succi
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185 Rome, Italy
- Center for Life Nanoscience at la Sapienza, Istituto Italiano di Tecnologia, viale Regina Elena 295, 00161 Rome, Italy
- Institute for Applied Computational Science, Harvard John A. Paulson School of Engineering and Applied Sciences, Cambridge, MA 02138, USA
| |
Collapse
|
11
|
Tiribocchi A, Bonaccorso F, Lauricella M, Melchionna S, Montessori A, Succi S. Curvature dynamics and long-range effects on fluid-fluid interfaces with colloids. Soft Matter 2019; 15:2848-2862. [PMID: 30816901 DOI: 10.1039/c8sm02396d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We investigate the dynamics of a phase-separating binary fluid, containing colloidal dumbbells anchored to the fluid-fluid interface. Extensive lattice Boltzmann-immersed boundary method simulations reveal that the presence of soft dumbbells can significantly affect the curvature dynamics of the interface between phase-separating fluids, even though the coarsening dynamics is left nearly unchanged. In addition, our results show that the curvature dynamics exhibits distinct non-local effects, which might be exploited for the design of new soft mesoscale materials. We point out that the inspection of the statistical dynamics of the curvature can disclose new insights into local inhomogeneities of the binary fluid configuration, as a function of the volume fraction and aspect ratio of the dumbbells.
Collapse
Affiliation(s)
- A Tiribocchi
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy. and Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185, Rome, Italy.
| | - F Bonaccorso
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy.
| | - M Lauricella
- Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185, Rome, Italy.
| | - S Melchionna
- ISC-CNR, Istituto Sistemi Complessi, Università Sapienza, P.le A. Moro 2, 00185 Rome, Italy.
| | - A Montessori
- Department of Engineering, University of Rome, "Roma Tre" Via Vito Volterra 62, 00146 Rome, Italy.
| | - S Succi
- Center for Life Nano Science@La Sapienza, Istituto Italiano di Tecnologia, 00161 Roma, Italy. and Istituto per le Applicazioni del Calcolo CNR, via dei Taurini 19, 00185, Rome, Italy. and Institute for Applied Computational Science, John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA.
| |
Collapse
|
12
|
Gabbana A, Polini M, Succi S, Tripiccione R, Pellegrino FMD. Prospects for the Detection of Electronic Preturbulence in Graphene. Phys Rev Lett 2018; 121:236602. [PMID: 30576199 DOI: 10.1103/physrevlett.121.236602] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/07/2018] [Indexed: 06/09/2023]
Abstract
Based on extensive numerical simulations, accounting for electrostatic interactions and dissipative electron-phonon scattering, we propose experimentally realizable geometries capable of sustaining electronic preturbulence in graphene samples. In particular, preturbulence is predicted to occur at experimentally attainable values of the Reynolds number between 10 and 50, over a broad spectrum of frequencies between 10 and 100 GHz.
Collapse
Affiliation(s)
- A Gabbana
- Università di Ferrara and INFN-Ferrara, I-44122 Ferrara, Italy
- Bergische Universität Wuppertal, D-42119 Wuppertal, Germany
| | - M Polini
- Istituto Italiano di Tecnologia, Graphene Labs, Via Morego 30, I-16163 Genova, Italy
| | - S Succi
- Center for Life Nano Science at La Sapienza, Italian Institute of Technology, Viale Regina Elena 295, I-00161 Roma, Italy
- Istituto Applicazioni del Calcolo, National Research Council of Italy, Via dei Taurini 19, I-00185 Roma, Italy
| | - R Tripiccione
- Università di Ferrara and INFN-Ferrara, I-44122 Ferrara, Italy
| | - F M D Pellegrino
- Dipartimento di Fisica e Astronomia, Università di Catania, Via S. Sofia 64, I-95123 Catania, Italy
- INFN, Sezione Catania, I-95123 Catania, Italy
| |
Collapse
|
13
|
Abstract
We present a numerical scheme to solve the Wigner equation, based on a lattice discretization of momentum space. The moments of the Wigner function are recovered exactly, up to the desired order given by the number of discrete momenta retained in the discretization, which also determines the accuracy of the method. The Wigner equation is equipped with an additional collision operator, designed in such a way as to ensure numerical stability without affecting the evolution of the relevant moments of the Wigner function. The lattice Wigner scheme is validated for the case of quantum harmonic and anharmonic potentials, showing good agreement with theoretical results. It is further applied to the study of the transport properties of one- and two-dimensional open quantum systems with potential barriers. Finally, the computational viability of the scheme for the case of three-dimensional open systems is also illustrated.
Collapse
Affiliation(s)
- S Solórzano
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Wolfgang-Pauli-Strasse 27, HIT, CH-8093 Zürich, Switzerland
| | - M Mendoza
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Wolfgang-Pauli-Strasse 27, HIT, CH-8093 Zürich, Switzerland
| | - S Succi
- Istituto per le Applicazioni del Calcolo C.N.R., Via dei Taurini, 19 00185 Rome, Italy and Institute for Advanced Computational Science, Harvard University, 29 Oxford Street, Cambridge, Massachusetts 02138, USA
| | - H J Herrmann
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Wolfgang-Pauli-Strasse 27, HIT, CH-8093 Zürich, Switzerland
| |
Collapse
|
14
|
Abstract
Despite a long record of intense effort, the basic mechanisms by which dissipation emerges from the microscopic dynamics of a relativistic fluid still elude complete understanding. In particular, several details must still be finalized in the pathway from kinetic theory to hydrodynamics mainly in the derivation of the values of the transport coefficients. In this paper, we approach the problem by matching data from lattice-kinetic simulations with analytical predictions. Our numerical results provide neat evidence in favor of the Chapman-Enskog [The Mathematical Theory of Non-Uniform Gases, 3rd ed. (Cambridge University Press, Cambridge, U.K., 1970)] procedure as suggested by recent theoretical analyses along with qualitative hints at the basic reasons why the Chapman-Enskog expansion might be better suited than Grad's method [Commun. Pure Appl. Math. 2, 331 (1949)0010-364010.1002/cpa.3160020403] to capture the emergence of dissipative effects in relativistic fluids.
Collapse
Affiliation(s)
- A Gabbana
- INFN-Ferrara, Università di Ferrara, Via Saragat 1, I-44122 Ferrara, Italy
| | - M Mendoza
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Schafmattstraße 6, HIF, CH-8093 Zürich, Switzerland
| | - S Succi
- Istituto per le Applicazioni del Calcolo C.N.R., Via dei Taurini, 19, 00185 Rome, Italy and Institute for Applied Computational Science, John Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA
| | - R Tripiccione
- INFN-Ferrara, Università di Ferrara, Via Saragat 1, I-44122 Ferrara, Italy
| |
Collapse
|
15
|
Abstract
We present a systematic derivation of relativistic lattice kinetic equations for finite-mass particles, reaching close to the zero-mass ultrarelativistic regime treated in the previous literature. Starting from an expansion of the Maxwell-Jüttner distribution on orthogonal polynomials, we perform a Gauss-type quadrature procedure and discretize the relativistic Boltzmann equation on space-filling Cartesian lattices. The model is validated through numerical comparison with standard tests and solvers in relativistic fluid dynamics such as Boltzmann approach multiparton scattering and previous relativistic lattice Boltzmann models. This work provides a significant step towards the formulation of a unified relativistic lattice kinetic scheme, covering both massive and near-massless particles regimes.
Collapse
Affiliation(s)
- A Gabbana
- Università di Ferrara and INFN-Ferrara, Via Saragat 1, I-44122 Ferrara, Italy
| | - M Mendoza
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Schafmattstraße 6, HIF, CH-8093 Zürich, Switzerland
| | - S Succi
- Istituto per le Applicazioni del Calcolo C.N.R., Via dei Taurini, 19 00185 Rome, Italy
| | - R Tripiccione
- Università di Ferrara and INFN-Ferrara, Via Saragat 1, I-44122 Ferrara, Italy
| |
Collapse
|
16
|
Abstract
In this work, a hybrid lattice Boltzmann method (HLBM) is proposed, where the standard lattice Boltzmann implementation based on the Bhatnagar-Gross-Krook (LBGK) approximation is combined together with an unstructured finite-volume lattice Boltzmann model. The method is constructed on an overlapping grid system, which allows the coexistence of a uniform lattice nodes spacing and a coordinate-free lattice structure. The natural adaptivity of the hybrid grid system makes the method particularly suitable to handle problems involving complex geometries. Moreover, the provided scheme ensures a high-accuracy solution near walls, given the capability of the unstructured submodel of achieving the desired level of refinement in a very flexible way. For these reasons, the HLBM represents a prospective tool for solving multiscale problems. The proposed method is here applied to the benchmark problem of a two-dimensional flow past a circular cylinder for a wide range of Reynolds numbers and its numerical performances are measured and compared with the standard LBGK ones.
Collapse
Affiliation(s)
- G Di Ilio
- University of Naples "Parthenope," Centro Direzionale Isola C4, 80133 Naples, Italy
| | - D Chiappini
- University of Rome "Niccolò Cusano," Via don Carlo Gnocchi 3, 00166 Rome, Italy
| | - S Ubertini
- University of Tuscia, Largo dell'Università snc, 01100, Viterbo, Italy
| | - G Bella
- University of Rome "Tor Vergata," Via del Politecnico 1, 00133 Rome, Italy
| | - S Succi
- Istituto Applicazioni Calcolo, CNR, Via dei Taurini 19, 00185 Rome, Italy
| |
Collapse
|
17
|
Di Staso G, Clercx HJH, Succi S, Toschi F. Lattice Boltzmann accelerated direct simulation Monte Carlo for dilute gas flow simulations. Philos Trans A Math Phys Eng Sci 2016; 374:rsta.2016.0226. [PMID: 27698045 DOI: 10.1098/rsta.2016.0226] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 06/27/2016] [Indexed: 05/08/2023]
Abstract
Hybrid particle-continuum computational frameworks permit the simulation of gas flows by locally adjusting the resolution to the degree of non-equilibrium displayed by the flow in different regions of space and time. In this work, we present a new scheme that couples the direct simulation Monte Carlo (DSMC) with the lattice Boltzmann (LB) method in the limit of isothermal flows. The former handles strong non-equilibrium effects, as they typically occur in the vicinity of solid boundaries, whereas the latter is in charge of the bulk flow, where non-equilibrium can be dealt with perturbatively, i.e. according to Navier-Stokes hydrodynamics. The proposed concurrent multiscale method is applied to the dilute gas Couette flow, showing major computational gains when compared with the full DSMC scenarios. In addition, it is shown that the coupling with LB in the bulk flow can speed up the DSMC treatment of the Knudsen layer with respect to the full DSMC case. In other words, LB acts as a DSMC accelerator.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.
Collapse
Affiliation(s)
- G Di Staso
- Department of Applied Physics, and J.M. Burgers Centre for Fluid Dynamics, Eindhoven University of Technology, Den Dolech 2, 5600 MB Eindhoven, The Netherlands
| | - H J H Clercx
- Department of Applied Physics, and J.M. Burgers Centre for Fluid Dynamics, Eindhoven University of Technology, Den Dolech 2, 5600 MB Eindhoven, The Netherlands
| | - S Succi
- Istituto per le Applicazioni del Calcolo, Consiglio Nazionale delle Ricerche, Via dei Taurini 19, 00185 Rome, Italy
| | - F Toschi
- Department of Applied Physics, and J.M. Burgers Centre for Fluid Dynamics, Eindhoven University of Technology, Den Dolech 2, 5600 MB Eindhoven, The Netherlands Istituto per le Applicazioni del Calcolo, Consiglio Nazionale delle Ricerche, Via dei Taurini 19, 00185 Rome, Italy Department of Mathematics and Computer Science, Eindhoven University of Technology, Den Dolech 2, 5600 MB Eindhoven, The Netherlands
| |
Collapse
|
18
|
Coveney PV, Boon JP, Succi S. Bridging the gaps at the physics-chemistry-biology interface. Philos Trans A Math Phys Eng Sci 2016; 374:rsta.2016.0335. [PMID: 27698047 PMCID: PMC5052737 DOI: 10.1098/rsta.2016.0335] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 08/17/2016] [Indexed: 05/13/2023]
Affiliation(s)
- P V Coveney
- Centre for Computational Science, University College London, Gordon Street, London WC1H 0AJ, UK
| | - J P Boon
- Physics Department, Université Libre de Bruxelles, Campus Plaine, CP 231, Avenue F.D. Roosevelt 50, 1050 Bruxelles, Belgium
| | - S Succi
- Istituto Applicazioni del Calcolo-CNR, Viale del Policlinico 19, 00185 Roma, Italy Institute for Applied Computational Science, Harvard J. Paulson School of Engineering and Applied Sciences, Harvard University, 29 Oxford Street, Cambridge, MA 02138, USA
| |
Collapse
|
19
|
Abstract
We discuss a unified mesoscale framework (chimaera) for the simulation of complex states of flowing matter across scales of motion. The chimaera framework can deal with each of the three macro-meso-micro levels through suitable 'mutations' of the basic mesoscale formulation. The idea is illustrated through selected simulations of complex micro- and nanoscale flows.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.
Collapse
Affiliation(s)
- S Succi
- Istituto Applicazioni del Calcolo-CNR, Viale del Policlinico 19, 00185 Roma, Italy Institute of Applied Computational Science, J. Paulson School of Applied Science and Engineering, Harvard University, 29 Oxford Street, Cambridge MA 02138, USA
| |
Collapse
|
20
|
Abstract
We investigate Poiseuille channel flow through intrinsically curved media, equipped with localized metric perturbations. To this end, we study the flux of a fluid driven through the curved channel in dependence of the spatial deformation, characterized by the parameters of the metric perturbations (amplitude, range, and density). We find that the flux depends only on a specific combination of parameters, which we identify as the average metric perturbation, and derive a universal flux law for the Poiseuille flow. For the purpose of this study, we have improved and validated our recently developed lattice Boltzmann model in curved space by considerably reducing discrete lattice effects.
Collapse
Affiliation(s)
- J-D Debus
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Wolfgang-Pauli-Strasse 27, HIT, CH-8093 Zürich, Switzerland
| | - M Mendoza
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Wolfgang-Pauli-Strasse 27, HIT, CH-8093 Zürich, Switzerland
| | - S Succi
- Instituto per le Applicazioni del Calcolo C.N.R., Via dei Taurini, 19 00185, Rome, Italy
| | - H J Herrmann
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Wolfgang-Pauli-Strasse 27, HIT, CH-8093 Zürich, Switzerland
| |
Collapse
|
21
|
Benzi R, Sbragaglia M, Bernaschi M, Succi S, Toschi F. Cooperativity flows and shear-bandings: a statistical field theory approach. Soft Matter 2016; 12:514-530. [PMID: 26486875 DOI: 10.1039/c5sm01862e] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Cooperativity effects have been proposed to explain the non-local rheology in the dynamics of soft jammed systems. Based on the analysis of the free-energy model proposed by L. Bocquet, A. Colin and A. Ajdari, Phys. Rev. Lett., 2009, 103, 036001, we show that cooperativity effects resulting from the non-local nature of the fluidity (inverse viscosity) are intimately related to the emergence of shear-banding configurations. This connection materializes through the onset of inhomogeneous compact solutions (compactons), wherein the fluidity is confined to finite-support subregions of the flow and strictly zero elsewhere. The compacton coexistence with regions of zero fluidity ("non-flowing vacuum") is shown to be stabilized by the presence of mechanical noise, which ultimately shapes up the equilibrium distribution of the fluidity field, the latter acting as an order parameter for the flow-noflow transitions occurring in the material.
Collapse
Affiliation(s)
- R Benzi
- Department of Physics and INFN, University of "Tor Vergata", Via della Ricerca Scientifica 1, 00133 Rome, Italy.
| | - M Sbragaglia
- Department of Physics and INFN, University of "Tor Vergata", Via della Ricerca Scientifica 1, 00133 Rome, Italy.
| | - M Bernaschi
- Istituto per le Applicazioni del Calcolo CNR, Via dei Taurini 19, 00185 Rome, Italy
| | - S Succi
- Istituto per le Applicazioni del Calcolo CNR, Via dei Taurini 19, 00185 Rome, Italy
| | - F Toschi
- Istituto per le Applicazioni del Calcolo CNR, Via dei Taurini 19, 00185 Rome, Italy and Department of Physics and Department of Mathematics and Computer Science and J. M. Burgerscentrum, Eindhoven University of Technology, 5600 MB, Eindhoven, Netherlands
| |
Collapse
|
22
|
Montessori A, Prestininzi P, La Rocca M, Succi S. Publisher's Note: Lattice Boltzmann approach for complex nonequilibrium flows [Phys. Rev. E 92, 043308 (2015)]. Phys Rev E Stat Nonlin Soft Matter Phys 2015; 92:069901. [PMID: 26764860 DOI: 10.1103/physreve.92.069901] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/19/2015] [Indexed: 06/05/2023]
Abstract
This corrects the article DOI: 10.1103/PhysRevE.92.043308.
Collapse
|
23
|
Sega M, Sbragaglia M, Biferale L, Succi S. The importance of chemical potential in the determination of water slip in nanochannels. Eur Phys J E Soft Matter 2015; 38:127. [PMID: 26614498 DOI: 10.1140/epje/i2015-15127-y] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/20/2015] [Accepted: 10/30/2015] [Indexed: 06/05/2023]
Abstract
We investigate the slip properties of water confined in graphite-like nanochannels by non-equilibrium molecular dynamics simulations, with the aim of identifying and analyze separately the influence of different physical quantities on the slip length. In a system under confinement but connected to a reservoir of fluid, the chemical potential is the natural control parameter: we show that two nanochannels characterized by the same macroscopic contact angle--but a different microscopic surface potential--do not exhibit the same slip length unless the chemical potential of water in the two channels is matched. Some methodological issues related to the preparation of samples for the comparative analysis in confined geometries are also discussed.
Collapse
Affiliation(s)
- M Sega
- Institute of Computational Physics, University of Vienna, Sensengasse 8/9, 1090, Vienna, Austria.
| | - M Sbragaglia
- Department of Physics and INFN, University of "Tor Vergata", Via della Ricerca Scientifica 1, 00133, Rome, Italy
| | - L Biferale
- Department of Physics and INFN, University of "Tor Vergata", Via della Ricerca Scientifica 1, 00133, Rome, Italy
| | - S Succi
- Istituto per le Applicazioni del Calcolo CNR, Via dei Taurini 19, 00185, Rome, Italy
| |
Collapse
|
24
|
Montessori A, Prestininzi P, La Rocca M, Succi S. Lattice Boltzmann approach for complex nonequilibrium flows. Phys Rev E Stat Nonlin Soft Matter Phys 2015; 92:043308. [PMID: 26565365 DOI: 10.1103/physreve.92.043308] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/11/2015] [Indexed: 05/08/2023]
Abstract
We present a lattice Boltzmann realization of Grad's extended hydrodynamic approach to nonequilibrium flows. This is achieved by using higher-order isotropic lattices coupled with a higher-order regularization procedure. The method is assessed for flow across parallel plates and three-dimensional flows in porous media, showing excellent agreement of the mass flow with analytical and numerical solutions of the Boltzmann equation across the full range of Knudsen numbers, from the hydrodynamic regime to ballistic motion.
Collapse
Affiliation(s)
- A Montessori
- Department of Engineering, University of Rome, "Roma Tre" Via Vito Volterra 62, 00146 Rome, Italy
| | - P Prestininzi
- Department of Engineering, University of Rome, "Roma Tre" Via Vito Volterra 62, 00146 Rome, Italy
| | - M La Rocca
- Department of Engineering, University of Rome, "Roma Tre" Via Vito Volterra 62, 00146 Rome, Italy
| | - S Succi
- Istituto per le Applicazioni del Calcolo, CNR Via dei Taurini 19, 00185 Rome, Italy
| |
Collapse
|
25
|
Mezzacapo A, Sanz M, Lamata L, Egusquiza IL, Succi S, Solano E. Quantum Simulator for Transport Phenomena in Fluid Flows. Sci Rep 2015; 5:13153. [PMID: 26278968 PMCID: PMC4538376 DOI: 10.1038/srep13153] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2015] [Accepted: 07/21/2015] [Indexed: 12/01/2022] Open
Abstract
Transport phenomena still stand as one of the most challenging problems in computational physics. By exploiting the analogies between Dirac and lattice Boltzmann equations, we develop a quantum simulator based on pseudospin-boson quantum systems, which is suitable for encoding fluid dynamics transport phenomena within a lattice kinetic formalism. It is shown that both the streaming and collision processes of lattice Boltzmann dynamics can be implemented with controlled quantum operations, using a heralded quantum protocol to encode non-unitary scattering processes. The proposed simulator is amenable to realization in controlled quantum platforms, such as ion-trap quantum computers or circuit quantum electrodynamics processors.
Collapse
Affiliation(s)
- A Mezzacapo
- Department of Physical Chemistry, University of the Basque Country UPV/EHU, Apartado 644, E-48080 Bilbao, Spain
| | - M Sanz
- Department of Physical Chemistry, University of the Basque Country UPV/EHU, Apartado 644, E-48080 Bilbao, Spain
| | - L Lamata
- Department of Physical Chemistry, University of the Basque Country UPV/EHU, Apartado 644, E-48080 Bilbao, Spain
| | - I L Egusquiza
- Department of Theoretical Physics and History of Science, University of the Basque Country UPV/EHU, Apartado 644, E-48080 Bilbao, Spain
| | - S Succi
- Istituto per le Applicazioni del Calcolo "M. Picone" CNR, I-00185 Rome, Italy.,Institute for Applied Computational Science, Harvard University, Oxford Street, 33, 02138 Cambridge, USA
| | - E Solano
- Department of Physical Chemistry, University of the Basque Country UPV/EHU, Apartado 644, E-48080 Bilbao, Spain.,IKERBASQUE, Basque Foundation for Science, Maria Diaz de Haro 3, 48013 Bilbao, Spain
| |
Collapse
|
26
|
Mohseni F, Mendoza M, Succi S, Herrmann HJ. Lattice Boltzmann model for resistive relativistic magnetohydrodynamics. Phys Rev E Stat Nonlin Soft Matter Phys 2015; 92:023309. [PMID: 26382548 DOI: 10.1103/physreve.92.023309] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/30/2015] [Indexed: 06/05/2023]
Abstract
In this paper, we develop a lattice Boltzmann model for relativistic magnetohydrodynamics (MHD). Even though the model is derived for resistive MHD, it is shown that it is numerically robust even in the high conductivity (ideal MHD) limit. In order to validate the numerical method, test simulations are carried out for both ideal and resistive limits, namely the propagation of Alfvén waves in the ideal MHD and the evolution of current sheets in the resistive regime, where very good agreement is observed comparing to the analytical results. Additionally, two-dimensional magnetic reconnection driven by Kelvin-Helmholtz instability is studied and the effects of different parameters on the reconnection rate are investigated. It is shown that the density ratio has a negligible effect on the magnetic reconnection rate, while an increase in shear velocity decreases the reconnection rate. Additionally, it is found that the reconnection rate is proportional to σ-1/2, σ being the conductivity, which is in agreement with the scaling law of the Sweet-Parker model. Finally, the numerical model is used to study the magnetic reconnection in a stellar flare. Three-dimensional simulation suggests that the reconnection between the background and flux rope magnetic lines in a stellar flare can take place as a result of a shear velocity in the photosphere.
Collapse
Affiliation(s)
- F Mohseni
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Wolfgang-Pauli-Strasse 27, HIT, CH-8093 Zürich, Switzerland
| | - M Mendoza
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Wolfgang-Pauli-Strasse 27, HIT, CH-8093 Zürich, Switzerland
| | - S Succi
- Istituto per le Applicazioni del Calcolo C.N.R., Via dei Taurini, 19 00185 Rome, Italy and Freiburg Institute for Advanced Studies, Albertstrasse 19, D-79104 Freiburg, Germany
| | - H J Herrmann
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Wolfgang-Pauli-Strasse 27, HIT, CH-8093 Zürich, Switzerland
- Departamento de Física, Universidade Federal do Ceará, Campus do Pici, 60455-760 Fortaleza, Ceará, Brazil
| |
Collapse
|
27
|
Benzi R, Sbragaglia M, Scagliarini A, Perlekar P, Bernaschi M, Succi S, Toschi F. Internal dynamics and activated processes in soft-glassy materials. Soft Matter 2015; 11:1271-1280. [PMID: 25560202 DOI: 10.1039/c4sm02341b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Plastic rearrangements play a crucial role in the characterization of soft-glassy materials, such as emulsions and foams. Based on numerical simulations of soft-glassy systems, we study the dynamics of plastic rearrangements at the hydrodynamic scales where thermal fluctuations can be neglected. Plastic rearrangements require an energy input, which can be either provided by external sources, or made available through time evolution in the coarsening dynamics, in which the total interfacial area decreases as a consequence of the slow evolution of the dispersed phase from smaller to large droplets/bubbles. We first demonstrate that our hydrodynamic model can quantitatively reproduce such coarsening dynamics. Then, considering periodically oscillating strains, we characterize the number of plastic rearrangements as a function of the external energy-supply, and show that they can be regarded as activated processes induced by a suitable "noise" effect. Here we use the word noise in a broad sense, referring to the internal non-equilibrium dynamics triggered by spatial random heterogeneities and coarsening. Finally, by exploring the interplay between the internal characteristic time-scale of the coarsening dynamics and the external time-scale associated with the imposed oscillating strain, we show that the system exhibits the phenomenon of stochastic resonance, thereby providing further credit to the mechanical activation scenario.
Collapse
Affiliation(s)
- R Benzi
- Department of Physics and INFN, University of "Tor Vergata", Via della Ricerca Scientifica 1, 00133 Rome, Italy.
| | | | | | | | | | | | | |
Collapse
|
28
|
Mendoza M, Succi S, Herrmann HJ. Kinetic formulation of the Kohn-Sham Equations for ab initio electronic structure calculations. Phys Rev Lett 2014; 113:096402. [PMID: 25215997 DOI: 10.1103/physrevlett.113.096402] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/10/2013] [Indexed: 06/03/2023]
Abstract
We introduce a new connection between density functional theory and kinetic theory. In particular, we show that the Kohn-Sham equations can be reformulated as a macroscopic limit of the steady-state solution of a suitable single-particle kinetic equation. We derive a Boltzmann-like equation for a gas of quasiparticles, where the potential plays the role of an external source that generates and destroys particles, so as to drive the system towards its ground state. The ions are treated as classical particles by using either the Born-Oppenheimer dynamics or by imposing concurrent evolution with the electronic orbitals. In order to provide quantitative support to our approach, we implement a discrete (lattice) kinetic model and compute the exchange and correlation energies of simple atoms and the geometrical configuration of the methane molecule. Moreover, we also compute the first vibrational mode of the hydrogen molecule, with both Born-Oppenheimer and concurrent dynamics. Excellent agreement with values in the literature is found in all cases.
Collapse
Affiliation(s)
- M Mendoza
- Computational Physics for Engineering Materials, Institute for Building Materials, ETH Zürich, Schafmattstrasse 6, HIF, CH-8093 Zürich, Switzerland
| | - S Succi
- Istituto per le Applicazioni del Calcolo C.N.R., Via dei Taurini, 19 00185 Rome, Italy and Institute for Advanced Computational Science, Harvard University, 29 Oxford Street, Cambridge, Massachusetts 02138, USA
| | - H J Herrmann
- Computational Physics for Engineering Materials, Institute for Building Materials, ETH Zürich, Schafmattstrasse 6, HIF, CH-8093 Zürich, Switzerland and Departamento de Física, Universidade Federal do Ceará, Campus do Pici, 60455-760 Fortaleza, Ceará, Brazil
| |
Collapse
|
29
|
Benzi R, Sbragaglia M, Perlekar P, Bernaschi M, Succi S, Toschi F. Direct evidence of plastic events and dynamic heterogeneities in soft-glasses. Soft Matter 2014; 10:4615-4624. [PMID: 24827455 DOI: 10.1039/c4sm00348a] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
By using fluid-kinetic simulations of confined and concentrated emulsion droplets, we investigate the nature of space non-homogeneity in soft-glassy dynamics and provide quantitative measurements of the statistical features of plastic events in the proximity of the yield-stress threshold. Above the yield stress, our results show the existence of a finite stress correlation scale, which can be mapped directly onto the cooperativity scale, recently introduced in the literature to capture non-local effects in the soft-glassy dynamics. In this regime, the emergence of a separate boundary (wall) rheology with higher fluidity than the bulk is highlighted in terms of near-wall spontaneous segregation of plastic events. Near the yield stress, where the cooperativity scale cannot be estimated with sufficient accuracy, the system shows a clear increase of the stress correlation scale, whereas plastic events exhibit intermittent clustering in time, with no preferential spatial location. A quantitative measurement of the space-time correlation associated with the motion of the interface of the droplets is key to spot the elastic rigidity of the system.
Collapse
Affiliation(s)
- R Benzi
- Department of Physics and INFN, University of "Tor Vergata", Via della Ricerca Scientifica 1, 00133 Rome, Italy.
| | | | | | | | | | | |
Collapse
|
30
|
Amiri Delouei A, Nazari M, Kayhani MH, Succi S. Non-Newtonian unconfined flow and heat transfer over a heated cylinder using the direct-forcing immersed boundary-thermal lattice Boltzmann method. Phys Rev E Stat Nonlin Soft Matter Phys 2014; 89:053312. [PMID: 25353919 DOI: 10.1103/physreve.89.053312] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/09/2013] [Indexed: 06/04/2023]
Abstract
In this study, the immersed boundary-thermal lattice Boltzmann method has been used to simulate non-Newtonian fluid flow over a heated circular cylinder. The direct-forcing algorithm has been employed to couple the off-lattice obstacles and on-lattice fluid nodes. To investigate the effect of boundary sharpness, two different diffuse interface schemes are considered to interpolate the velocity and temperature between the boundary and computational grid points. The lattice Boltzmann equation with split-forcing term is applied to consider the effects of the discrete lattice and the body force to the momentum flux, simultaneously. A method for calculating the Nusselt number based on diffuse interface schemes is developed. The rheological and thermal properties of non-Newtonian fluids are investigated under the different power-law indices and Reynolds numbers. The effect of numerical parameters on the accuracy of the proposed method has been investigated in detail. Results show that the rheological and thermal properties of non-Newtonian fluids in the presence of a heated immersed body can be suitably captured using the immersed boundary thermal lattice Boltzmann method.
Collapse
Affiliation(s)
- A Amiri Delouei
- Mechanical Engineering Department, Shahrood University of Technology, Shahrood, Iran
| | - M Nazari
- Mechanical Engineering Department, Shahrood University of Technology, Shahrood, Iran
| | - M H Kayhani
- Mechanical Engineering Department, Shahrood University of Technology, Shahrood, Iran
| | - S Succi
- IAC-CNR, Rome, Via dei Taurini 19, 00185, Roma, Italy
| |
Collapse
|
31
|
Montessori A, Falcucci G, Prestininzi P, La Rocca M, Succi S. Regularized lattice Bhatnagar-Gross-Krook model for two- and three-dimensional cavity flow simulations. Phys Rev E Stat Nonlin Soft Matter Phys 2014; 89:053317. [PMID: 25353924 DOI: 10.1103/physreve.89.053317] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/11/2013] [Indexed: 06/04/2023]
Abstract
We investigate the accuracy and performance of the regularized version of the single-relaxation-time lattice Boltzmann equation for the case of two- and three-dimensional lid-driven cavities. The regularized version is shown to provide a significant gain in stability over the standard single-relaxation time, at a moderate computational overhead.
Collapse
Affiliation(s)
- A Montessori
- Department of Engineering, University of Rome "Roma Tre," Via della Vasca Navale 79, 00141 Rome, Italy
| | - G Falcucci
- Department of Engineering, University of Naples "Parthenope," Centro Direzionale, Isola C4, 80149 Naples, Italy
| | - P Prestininzi
- Department of Engineering, University of Rome "Roma Tre," Via della Vasca Navale 79, 00141 Rome, Italy
| | - M La Rocca
- Department of Engineering, University of Rome "Roma Tre," Via della Vasca Navale 79, 00141 Rome, Italy
| | - S Succi
- Istituto per le Applicazioni del Calcolo, CNR, Via dei Taurini 19, 00185 Rome, Italy
| |
Collapse
|
32
|
Abstract
We present a computational study of the transport properties of campylotic (intrinsically curved) media. It is found that the relation between the flow through a campylotic media, consisting of randomly located curvature perturbations, and the average Ricci scalar of the system, exhibits two distinct functional expressions, depending on whether the typical spatial extent of the curvature perturbation lies above or below the critical value maximizing the overall scalar of curvature. Furthermore, the flow through such systems as a function of the number of curvature perturbations is found to present a sublinear behavior for large concentrations, due to the interference between curvature perturbations leading to an overall less curved space. We have also characterized the flux through such media as a function of the local Reynolds number and the scale of interaction between impurities. For the purpose of this study, we have also developed and validated a new lattice Boltzmann model.
Collapse
Affiliation(s)
- M Mendoza
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Wolfgang-Pauli-Strasse 27, HIT, CH-8093 Zürich (Switzerland)
| | | | | |
Collapse
|
33
|
Fillion-Gourdeau F, Herrmann HJ, Mendoza M, Palpacelli S, Succi S. Formal analogy between the Dirac equation in its Majorana form and the discrete-velocity version of the Boltzmann kinetic equation. Phys Rev Lett 2013; 111:160602. [PMID: 24182245 DOI: 10.1103/physrevlett.111.160602] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/01/2013] [Indexed: 06/02/2023]
Abstract
We point out a formal analogy between the Dirac equation in Majorana form and the discrete-velocity version of the Boltzmann kinetic equation. By a systematic analysis based on the theory of operator splitting, this analogy is shown to turn into a concrete and efficient computational method, providing a unified treatment of relativistic and nonrelativistic quantum mechanics. This might have potentially far-reaching implications for both classical and quantum computing, because it shows that, by splitting time along the three spatial directions, quantum information (Dirac-Majorana wave function) propagates in space-time as a classical statistical process (Boltzmann distribution).
Collapse
Affiliation(s)
- F Fillion-Gourdeau
- Centre de Recherches Mathématiques, Université de Montréal, Montréal, Canada H3T 1J4
| | | | | | | | | |
Collapse
|
34
|
Mendoza M, Araújo NAM, Succi S, Herrmann HJ. Transition in the equilibrium distribution function of relativistic particles. Sci Rep 2012; 2:611. [PMID: 22937220 PMCID: PMC3430878 DOI: 10.1038/srep00611] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2012] [Accepted: 08/01/2012] [Indexed: 11/23/2022] Open
Abstract
We analyze a transition from single peaked to bimodal velocity distribution in a relativistic fluid under increasing temperature, in contrast with a non-relativistic gas, where only a monotonic broadening of the bell-shaped distribution is observed. Such transition results from the interplay between the raise in thermal energy and the constraint of maximum velocity imposed by the speed of light. We study the Bose-Einstein, the Fermi-Dirac, and the Maxwell-Jüttner distributions, and show that they all exhibit the same qualitative behavior. We characterize the nature of the transition in the framework of critical phenomena and show that it is either continuous or discontinuous, depending on the group velocity. We analyze the transition in one, two, and three dimensions, with special emphasis on twodimensions, for which a possible experiment in graphene, based on the measurement of the Johnson-Nyquist noise, is proposed.
Collapse
Affiliation(s)
- M Mendoza
- Computational Physics for Engineering Materials, IfB, ETH Zürich, Schafmattstrasse 6, CH-8093 Zürich, Switzerland.
| | | | | | | |
Collapse
|
35
|
|
36
|
Karlin IV, Succi S, Chikatamarla SS. Comment on "Numerics of the lattice Boltzmann method: effects of collision models on the lattice Boltzmann simulations". Phys Rev E Stat Nonlin Soft Matter Phys 2011; 84:068701. [PMID: 22304222 DOI: 10.1103/physreve.84.068701] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/14/2011] [Revised: 09/06/2011] [Indexed: 05/31/2023]
Abstract
Critical comments on the entropic lattice Boltzmann equation (ELBE), by Li-Shi Luo, Wei Liao, Xingwang Chen, Yan Peng, and Wei Zhang in Ref. , are based on simulations, which make use of a model that, despite being referred to as the ELBE by the authors, is in fact equivalent to the standard lattice Bhatnagar-Gross-Krook equation for low Mach number simulations. In this Comment, a concise review of the ELBE is provided and illustrated by means of a three-dimensional turbulent flow simulation, which highlights the subgrid features of the ELBE.
Collapse
Affiliation(s)
- I V Karlin
- Energy Technology Research Group, School of Engineering Sciences, University of Southampton, Southampton GB-SO17 1BJ, United Kingdom.
| | | | | |
Collapse
|
37
|
Benzi R, Bernaschi M, Sbragaglia M, Succi S. Heterogeneous diffuse interfaces: a new mechanism for arrested coarsening in binary mixtures. Heterogeneous diffuse interfaces. Eur Phys J E Soft Matter 2011; 34:93. [PMID: 21947891 DOI: 10.1140/epje/i2011-11093-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/18/2011] [Revised: 06/15/2011] [Accepted: 07/15/2011] [Indexed: 05/31/2023]
Abstract
We discuss the dynamics of binary fluid mixtures in which surface tension density is allowed to become locally negative within the interface, while still preserving positivity of the overall surface tension (heterogeneous diffuse interface). Numerical simulations of two-dimensional Ginzburg-Landau phase field equations implementing such mechanism and including hydrodynamic motion, show evidence of dynamically arrested domain coarsening. Under specific conditions on the functional form of the surface tension density, dynamical arrest can be interpreted in terms of the collective dynamics of metastable, non-linear excitations of the density field, named compactons, as they are localized to finite-size regions of configuration space and strictly zero elsewhere. Aside from compactons, the heterogeneous diffuse interface scenario appears to provide a robust mechanism for the interpretation of many aspects of soft-glassy behaviour in binary fluid mixtures.
Collapse
Affiliation(s)
- R Benzi
- Physics Department, University of Roma, Tor Vergata, via della Ricerca Scientifica, 1, 00133, Roma, Italy
| | | | | | | |
Collapse
|
38
|
Benzi R, Sbragaglia M, Bernaschi M, Succi S. Phase-field model of long-time glasslike relaxation in binary fluid mixtures. Phys Rev Lett 2011; 106:164501. [PMID: 21599369 DOI: 10.1103/physrevlett.106.164501] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/02/2010] [Indexed: 05/30/2023]
Abstract
We present a new phase-field model for binary fluids, exhibiting typical signatures of soft-glassy behavior, such as long-time relaxation, aging, and long-term dynamical arrest. The present model allows the cost of building an interface to vanish locally within the interface, while preserving positivity of the overall surface tension. A crucial consequence of this property, which we prove analytically, is the emergence of free-energy minimizing density configurations, hereafter named "compactons," to denote their property of being localized to a finite-size region of space and strictly zero elsewhere (no tails). Thanks to compactness, any arbitrary superposition of compactons still is a free-energy minimizer, which provides a direct link between the complexity of the free-energy landscape and the morphological complexity of configurational space.
Collapse
Affiliation(s)
- R Benzi
- Department of Physics and INFN, University of Tor Vergata, Rome, Italy
| | | | | | | |
Collapse
|
39
|
Abstract
We provide numerical evidence that electronic preturbulent phenomena in graphene could be observed, under current experimental conditions, through current fluctuations, echoing the detachment of vortices past localized micron-sized impurities. Vortex generation, due to micron-sized constriction, is also explored with special focus on the effects of relativistic corrections to the normal Navier-Stokes equations. These corrections are found to cause a delay in the stability breakout of the fluid as well as a small shift in the vortex shedding frequency.
Collapse
Affiliation(s)
- M Mendoza
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Schafmattstrasse 6, HIF, CH-8093 Zürich, Switzerland.
| | | | | |
Collapse
|
40
|
Dellar PJ, Lapitski D, Palpacelli S, Succi S. Isotropy of three-dimensional quantum lattice Boltzmann schemes. Phys Rev E Stat Nonlin Soft Matter Phys 2011; 83:046706. [PMID: 21599333 DOI: 10.1103/physreve.83.046706] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/08/2010] [Indexed: 05/30/2023]
Abstract
Numerical simulations with previous formulations of the quantum lattice Boltzmann (QLB) scheme in three spatial dimensions showed significant lack of isotropy. In two or more spatial dimensions the QLB approach relies upon operator splitting to decompose the time evolution into a sequence of applications of the one-dimensional QLB scheme along coordinate axes. Each application must be accompanied by a rotation of the wave function into a basis of chiral eigenstates aligned along the relevant axis. The previously observed lack of isotropy was due to an inconsistency in the application of these rotations. Once this inconsistency is removed, the QLB scheme is shown to exhibit isotropic behavior to within a numerical error that scales approximately linearly with the lattice spacing. This establishes the viability of the QLB approach in two and three spatial dimensions.
Collapse
Affiliation(s)
- P J Dellar
- OCIAM, Mathematical Institute, 24-29 St Giles', Oxford OX1 3LB, United Kingdom.
| | | | | | | |
Collapse
|
41
|
Gizzi A, Bernaschi M, Bini D, Cherubini C, Filippi S, Melchionna S, Succi S. Three-band decomposition analysis of wall shear stress in pulsatile flows. Phys Rev E Stat Nonlin Soft Matter Phys 2011; 83:031902. [PMID: 21517520 DOI: 10.1103/physreve.83.031902] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/23/2010] [Revised: 11/23/2010] [Indexed: 05/30/2023]
Abstract
Space-time patterns of wall shear stress (WSS) resulting from the numerical simulation of pulsating hemodynamic flows in semicoronal domains are analyzed, in the case of both regular semicoronal domains and semicoronal domains with bumpy insertions, mimicking aneurysm-like geometries. A new family of cardiovascular risk indicators, which we name three-band diagrams (TBDs), are introduced, as a sensible generalization of the two standard indicators, i.e., the time-averaged WSS and the oscillatory shear index. TBDs provide a handy access to additional information contained in the dynamic structure of the WSS signal as a function of the physiological risk threshold, thereby allowing a quick visual assessment of the risk sensitivity to individual fluctuations of the physiological risk thresholds. Due to its generality, TBD analysis is expected to prove useful for a wide host of applications in science, engineering, and medicine, where risk assessment plays a central role.
Collapse
Affiliation(s)
- A Gizzi
- Nonlinear Physics and Mathematical Modeling Lab, University Campus Bio-Medico, I-00128 Rome, Italy
| | | | | | | | | | | | | |
Collapse
|
42
|
Mendoza M, Boghosian BM, Herrmann HJ, Succi S. Fast lattice Boltzmann solver for relativistic hydrodynamics. Phys Rev Lett 2010; 105:014502. [PMID: 20867451 DOI: 10.1103/physrevlett.105.014502] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2009] [Indexed: 05/29/2023]
Abstract
A lattice Boltzmann formulation for relativistic fluids is presented and numerically validated through quantitative comparison with recent hydrodynamic simulations of relativistic fluids. In order to illustrate its capability to handle complex geometries, the scheme is also applied to the case of a three-dimensional relativistic shock wave, generated by a supernova explosion, impacting on a massive interstellar cloud. This formulation opens up the possibility of exporting the proven advantages of lattice Boltzmann methods, namely, computational efficiency and easy handling of complex geometries, to the context of (mildly) relativistic fluid dynamics at large, from quark-gluon plasmas up to supernovae with relativistic outflows.
Collapse
Affiliation(s)
- M Mendoza
- ETH Zürich, Computational Physics for Engineering Materials, Institute for Building Materials, Schafmattstrasse 6, HIF, CH-8093 Zürich, Switzerland.
| | | | | | | |
Collapse
|
43
|
Ubertini S, Asinari P, Succi S. Three ways to lattice Boltzmann: a unified time-marching picture. Phys Rev E Stat Nonlin Soft Matter Phys 2010; 81:016311. [PMID: 20365464 DOI: 10.1103/physreve.81.016311] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/25/2009] [Revised: 08/31/2009] [Indexed: 05/29/2023]
Abstract
It is shown that the lattice Boltzmann equation (LBE) corresponds to an explicit Verlet time-marching scheme for a continuum generalized Boltzmann equation with a memory delay equal to a half time step. This proves second-order accuracy of LBE with respect to this generalized equation, with no need of resorting to any implicit time-marching procedure (Crank-Nicholson) and associated nonlinear variable transformations. It is also shown, and numerically demonstrated, that this equivalence is not only formal, but it also translates into a complete equivalence of the corresponding computational schemes with respect to the hydrodynamic equations. Second-order accuracy with respect to the continuum kinetic equation is also numerically demonstrated for the case of the Taylor-Green vortex. It is pointed out that the equivalence is however broken for the case in which mass and/or momentum are not conserved, such as for chemically reactive flows and mixtures. For such flows, the time-centered implicit formulation may indeed offer a better numerical accuracy.
Collapse
Affiliation(s)
- S Ubertini
- Dipartimento per le Tecnologie (DiT), Centro Direzionale, Università di Napoli Parthenope, Napoli, Italy
| | | | | |
Collapse
|
44
|
Chibbaro S, Costa E, Dimitrov DI, Diotallevi F, Milchev A, Palmieri D, Pontrelli G, Succi S. Capillary filling in microchannels with wall corrugations: a comparative study of the Concus-Finn criterion by continuum, kinetic, and atomistic approaches. Langmuir 2009; 25:12653-12660. [PMID: 19817349 DOI: 10.1021/la901993r] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
We study the impact of wall corrugations in microchannels on the process of capillary filling by means of three broadly used methods: computational fluid dynamics (CFD), lattice Boltzmann equations (LBE), and molecular dynamics (MD). The numerical results of these approaches are compared and tested against the Concus-Finn (CF) criterion, which predicts pinning of the contact line at rectangular ridges perpendicular to flow for contact angles of theta > 45 degrees . Whereas for theta = 30, 40 (no flow), and 60 degrees (flow) all methods are found to produce data consistent with the CF criterion, at theta = 50 degrees the numerical experiments provide different results. Whereas the pinning of the liquid front is observed both in the LB and CFD simulations, MD simulations show that molecular fluctuations allow front propagation even above the critical value predicted by the deterministic CF criterion, thereby introducing a sensitivity to the obstacle height.
Collapse
Affiliation(s)
- S Chibbaro
- Dip. Fisica, Universitá di Roma TRE, Roma, Italy.
| | | | | | | | | | | | | | | |
Collapse
|
45
|
Falcucci G, Chiatti G, Succi S, Mohamad AA, Kuzmin A. Rupture of a ferrofluid droplet in external magnetic fields using a single-component lattice Boltzmann model for nonideal fluids. Phys Rev E Stat Nonlin Soft Matter Phys 2009; 79:056706. [PMID: 19518592 DOI: 10.1103/physreve.79.056706] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2009] [Revised: 03/23/2009] [Indexed: 05/27/2023]
Abstract
A nonisotropic tensorial extension of the single-component Shan-Chen pseudopotential Lattice Boltzmann method for nonideal fluids is presented. Direct comparison with experimental data shows that this extension is able to capture relevant features of ferrofluid behavior, such as the deformation and subsequent rupture of a liquid droplet as a function of an externally applied magnetic field. The present model offers an economic lattice-kinetic pathway to the simulation of complex ferrofluid hydrodynamics.
Collapse
Affiliation(s)
- G Falcucci
- Department of Mechanical and Industrial Engineering, University of Roma Tre, Via della Vasca Navale 79, 00146 Rome, Italy
| | | | | | | | | |
Collapse
|
46
|
Benzi R, Sbragaglia M, Succi S, Bernaschi M, Chibbaro S. Mesoscopic lattice Boltzmann modeling of soft-glassy systems: Theory and simulations. J Chem Phys 2009. [DOI: 10.1063/1.3216105] [Citation(s) in RCA: 78] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
47
|
Adhikari R, Succi S. Duality in matrix lattice Boltzmann models. Phys Rev E Stat Nonlin Soft Matter Phys 2008; 78:066701. [PMID: 19256975 DOI: 10.1103/physreve.78.066701] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/01/2007] [Revised: 09/04/2008] [Indexed: 05/27/2023]
Abstract
The notion of duality between the hydrodynamic and kinetic (ghost) variables of lattice kinetic formulations of the Boltzmann equation is introduced. It is suggested that this notion can serve as a guideline in the design of matrix versions of the lattice Boltzmann equation in a physically transparent and computationally efficient way.
Collapse
Affiliation(s)
- R Adhikari
- The Institute of Mathematical Sciences, CIT Campus, Tharamani, Chennai 600113, India
| | | |
Collapse
|
48
|
Diotallevi F, Biferale L, Chibbaro S, Puglisi A, Succi S. Front pinning in capillary filling of chemically coated channels. Phys Rev E Stat Nonlin Soft Matter Phys 2008; 78:036305. [PMID: 18851141 DOI: 10.1103/physreve.78.036305] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/08/2008] [Indexed: 05/26/2023]
Abstract
The dynamics of capillary filling in the presence of chemically coated heterogeneous boundaries is investigated both theoretically and numerically. In particular, by mapping the equations of front motion onto the dynamics of a dissipative driven oscillator, an analytical criterion for front pinning is derived under the condition of diluteness of the coating spots. The criterion is tested against two-dimensional lattice Boltzmann simulations and found to provide satisfactory agreement as long as the width of the front interface remains much thinner than the typical heterogeneity scale of the chemical coating.
Collapse
Affiliation(s)
- F Diotallevi
- Istituto per le Applicazioni del Calcolo CNR V. Policlinico 137, 00161 Roma, Italy.
| | | | | | | | | |
Collapse
|
49
|
Palpacelli S, Succi S. Quantum lattice Boltzmann simulation of expanding Bose-Einstein condensates in random potentials. Phys Rev E Stat Nonlin Soft Matter Phys 2008; 77:066708. [PMID: 18643398 DOI: 10.1103/physreve.77.066708] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/18/2008] [Indexed: 05/26/2023]
Abstract
The phenomenon of Anderson localization in expanding one-dimensional Bose-Einstein condensates is investigated by numerically solving the Gross-Pitaevskii equation with a random speckle potential. To this purpose, a quantum lattice Boltzmann (QLB) method is used, and compared with a standard Crank-Nicolson scheme. The QLB simulations show evidence of Anderson localization even for relatively low-energy condensates, with a healing length as large as one-tenth of the Thomas-Fermi length. Moreover, very long-time simulations, lasting up to 15 000 optical confinement periods, indicate that the Anderson localization degrades in time, although at a very slow pace. In particular, the inverse localization length is found to decay according to a t;{-1/3} law. This lends support to the idea that localized wave functions, although not strictly ground states, represent extremely long-lived metastable states of the expanding condensate.
Collapse
Affiliation(s)
- S Palpacelli
- Dipartimento di Matematica, Università Roma Tre, Largo San Leonardo Murialdo 1, 00146 Roma, Italy
| | | |
Collapse
|
50
|
Succi S. Applied Lattice Boltzmann Method for Transport Phenomena, Momentum, Heat and Mass Transfer. A. A. Mohamad Sure Printing, Calgary, AB April 2007. CAN J CHEM ENG 2008. [DOI: 10.1002/cjce.5450850617] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|