1
|
Ambruş VE, Bazzanini L, Gabbana A, Simeoni D, Succi S, Tripiccione R. Fast kinetic simulator for relativistic matter. NATURE COMPUTATIONAL SCIENCE 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] [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
|
Schullian O, Antila HS, Heazlewood BR. A variable time step self-consistent mean field DSMC model for three-dimensional environments. J Chem Phys 2022; 156:124309. [PMID: 35364882 DOI: 10.1063/5.0083033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A self-consistent mean field direct simulation Monte Carlo (SCMFD) algorithm was recently proposed for simulating collision environments for a range of one-dimensional model systems. This work extends the one-dimensional SCMFD approach to three dimensions and introduces a variable time step (3D-vt-SCMFD), enabling the modeling of a considerably wider range of different collision environments. We demonstrate the performance of the augmented method by modeling a varied set of test systems: ideal gas mixtures, Poiseuille flow of argon, and expansion of gas into high vacuum. For the gas mixtures, the 3D-vt-SCMFD method reproduces the properties (mean free path, mean free time, collision frequency, and temperature) in excellent agreement with theoretical predictions. From the Poiseuille flow simulations, we extract flow profiles that agree with the solution to the Navier-Stokes equations in the high-density limit and resemble free molecular flow at low densities, as expected. The measured viscosity from 3D-vt-SCMF is ∼15% lower than the theoretical prediction from Chapman-Enskog theory. The expansion of gas into vacuum is examined in the effusive regime and at the hydrodynamic limit. In both cases, 3D-vt-SCMDF simulations produce gas beam density, velocity, and temperature profiles in excellent agreement with analytical models. In summary, our tests show that 3D-vt-SCMFD is robust and computationally efficient, while also illustrating the diversity of systems the SCMFD model can be successfully applied to.
Collapse
Affiliation(s)
- O Schullian
- Department of Theory and Bio-Systems, Max Planck Institute of Colloids and Interfaces, Potsdam-Golm Science Park, 14476 Potsdam, Germany
| | - H S Antila
- Department of Theory and Bio-Systems, Max Planck Institute of Colloids and Interfaces, Potsdam-Golm Science Park, 14476 Potsdam, Germany
| | - B R Heazlewood
- Department of Physics, University of Liverpool, Liverpool L69 7ZE, United Kingdom
| |
Collapse
|
3
|
Kallikounis NG, Dorschner B, Karlin IV. Multiscale semi-Lagrangian lattice Boltzmann method. Phys Rev E 2021; 103:063305. [PMID: 34271620 DOI: 10.1103/physreve.103.063305] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2021] [Accepted: 05/20/2021] [Indexed: 11/07/2022]
Abstract
We present a multi-scale lattice Boltzmann scheme, which adaptively refines particles' velocity space. Different velocity sets of lower and higher order are consistently and efficiently coupled, allowing us to use the higher-order model only when and where needed. This includes regions of high Mach or high Knudsen numbers. The coupling procedure of discrete velocity sets consists of either a projection of the higher-order populations onto the lower-order lattice or lifting of the lower-order populations to the higher-order velocity space. Both lifting and projection are local operations, which enable a flexible adaptive velocity set. The proposed scheme is formulated for both a static and an optimal, co-moving reference frame, in the spirit of the recently introduced Particles on Demand method. The multi-scale scheme is validated with an advection of an athermal vortex and in a jet flow setup. The performance of the proposed scheme is further investigated in the shock structure problem and a high-Knudsen-number Couette flow, typical examples of highly non-equilibrium flows in which the order of the velocity set plays a decisive role. The results demonstrate that the proposed multi-scale scheme can operate accurately, with flexibility in terms of the underlying models and with reduced computational requirements.
Collapse
Affiliation(s)
- N G Kallikounis
- Department of Mechanical and Process Engineering, ETH Zurich, 8092 Zurich, Switzerland
| | - B Dorschner
- Department of Mechanical and Process Engineering, ETH Zurich, 8092 Zurich, Switzerland
| | - I V Karlin
- Department of Mechanical and Process Engineering, ETH Zurich, 8092 Zurich, Switzerland
| |
Collapse
|
4
|
Gan Y, Xu A, Zhang G, Zhang Y, Succi S. Discrete Boltzmann trans-scale modeling of high-speed compressible flows. Phys Rev E 2018; 97:053312. [PMID: 29906918 DOI: 10.1103/physreve.97.053312] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2018] [Indexed: 06/08/2023]
Abstract
We present a general framework for constructing trans-scale discrete Boltzmann models (DBMs) for high-speed compressible flows ranging from continuum to transition regime. This is achieved by designing a higher-order discrete equilibrium distribution function that satisfies additional nonhydrodynamic kinetic moments. To characterize the thermodynamic nonequilibrium (TNE) effects and estimate the condition under which the DBMs at various levels should be used, two measures are presented: (i) the relative TNE strength, describing the relative strength of the (N+1)th order TNE effects to the Nth order one; (ii) the TNE discrepancy between DBM simulation and relevant theoretical analysis. Whether or not the higher-order TNE effects should be taken into account in the modeling and which level of DBM should be adopted is best described by the relative TNE intensity and/or the discrepancy rather than by the value of the Knudsen number. As a model example, a two-dimensional DBM with 26 discrete velocities at Burnett level is formulated, verified, and validated.
Collapse
Affiliation(s)
- Yanbiao Gan
- North China Institute of Aerospace Engineering, Langfang 065000, China
- College of Mathematics and Informatics & FJKLMAA, Fujian Normal University, Fuzhou 350007, China
| | - Aiguo Xu
- National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, P.O. Box 8009-26, Beijing 100088, China
- Center for Applied Physics and Technology, MOE Key Center for High Energy Density Physics Simulations, College of Engineering, Peking University, Beijing 100871, China
| | - Guangcai Zhang
- National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, P.O. Box 8009-26, Beijing 100088, China
| | - Yudong Zhang
- National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, P.O. Box 8009-26, Beijing 100088, China
- Key Laboratory of Transient Physics, Nanjing University of Science and Technology, Nanjing 210094, China
| | - Sauro Succi
- Center for Life Nano Science at La Sapienza, Fondazione Istituto Italiano di Tecnologia, Viale Regina Margherita 295, 00161 Roma, Italy
- Physics Department and Institute for Applied Computational Science, John A. Paulson School of Applied Science and Engineering, Harvard University, Oxford Street 29, Cambridge, Massachusetts 02138, USA
| |
Collapse
|
5
|
Coveney PV, Boon JP, Succi S. Bridging the gaps at the physics-chemistry-biology interface. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 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] [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
|
6
|
Succi S. Chimaera simulation of complex states of flowing matter. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2016; 374:rsta.2016.0151. [PMID: 27698031 PMCID: PMC5052734 DOI: 10.1098/rsta.2016.0151] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 08/03/2016] [Indexed: 05/27/2023]
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
|