-
Flow force calculation in Lattice Boltzmann Method
Authors:
Shaurya Kaushal,
Sauro Succi,
Santosh Ansumali
Abstract:
We revisit force evaluation methodologies on rigid solid particles suspended in a viscous fluid and simulated via lattice Boltzmann method (LBM). We point out the non-commutativity of streaming and collision operators in the force evaluation procedure and provide a theoretical explanation for this observation. Based on this analysis, we propose a discrete force calculation scheme with enhanced acc…
▽ More
We revisit force evaluation methodologies on rigid solid particles suspended in a viscous fluid and simulated via lattice Boltzmann method (LBM). We point out the non-commutativity of streaming and collision operators in the force evaluation procedure and provide a theoretical explanation for this observation. Based on this analysis, we propose a discrete force calculation scheme with enhanced accuracy. The proposed scheme is essentially a lattice version of the Reynolds Transport Theorem (RTT) in the context of the lattice Boltzmann formulation. Besides maintaining satisfactory levels of reliability and accuracy, the method also handles force evaluation on complex geometries in a simple and transparent way. We run simulations for NACA0012 airfoil for a range of Reynolds numbers ranging from 100 to 0.5x10^6 and show that the current approach significantly reduces the grid size requirement for accurate force evaluation.
△ Less
Submitted 19 October, 2022;
originally announced October 2022.
-
LES/DNS of flow past T106 LPT cascade using a higher-order LB model
Authors:
Naliganahalli Hanumantharayappa Maruthi,
Chakradhar Thantanapally,
Manjusha Namburi,
Viswanathan Kumaran,
Santosh Ansumali
Abstract:
The main objective of the present work is to assess higher-order Entropic Lattice Boltzmann Method (ELBM) for separated and transitional flows without the use of any explicit turbulence model. For this, we chose to simulate two cases of T106 Low-Pressure Turbine (LPT) cascade -- T106A and T106C -- representing incompressible and compressible flow regimes respectively. These results are obtained us…
▽ More
The main objective of the present work is to assess higher-order Entropic Lattice Boltzmann Method (ELBM) for separated and transitional flows without the use of any explicit turbulence model. For this, we chose to simulate two cases of T106 Low-Pressure Turbine (LPT) cascade -- T106A and T106C -- representing incompressible and compressible flow regimes respectively. These results are obtained using our company's in-house higher-order ELBM transonic solver. We have carried out two sets of simulations for both the test cases. One with a clean inlet and the other with an inlet disturbance given by white Gaussian noise superimposed on the inlet velocity. For the clean inlet case, the flow remains laminar on the entire blade surface for both the test cases. It undergoes transition on the suction side for the inlet disturbance case. The pressure coefficient for the T106A and the isentropic Mach number on the blade surface for the T106C matches well with the experimental results. Also, the qualitative comparison of flow features in terms of the Q-criterion is in good agreement with the earlier computational results reported in literature.
△ Less
Submitted 27 March, 2022;
originally announced April 2022.
-
Two-fluid kinetic theory for dilute polymer solutions
Authors:
Shiwani Singh,
Ganesh Subramanian,
Santosh Ansumali
Abstract:
We provide a Boltzmann-type kinetic description for dilute polymer solutions based on two-fluid theory. This Boltzmann-type description uses a quasi-equilibrium based relaxation mechanism to model collisions between a polymer dumbbell and a solvent molecule. The model reproduces the desired macroscopic equations for the polymer-solvent mixture. The proposed kinetic scheme leads to a numerical algo…
▽ More
We provide a Boltzmann-type kinetic description for dilute polymer solutions based on two-fluid theory. This Boltzmann-type description uses a quasi-equilibrium based relaxation mechanism to model collisions between a polymer dumbbell and a solvent molecule. The model reproduces the desired macroscopic equations for the polymer-solvent mixture. The proposed kinetic scheme leads to a numerical algorithm which is along the lines of the lattice Boltzmann method. Finally, the algorithm is applied to describe the evolution of a perturbed Kolmogorov flow profile, whereby we recover the major elastic effect exhibited by a polymer solution, specifically, the suppression of the original inertial instability.
△ Less
Submitted 25 January, 2022;
originally announced January 2022.
-
Reduced kinetic model of polyatomic gases
Authors:
Praveen Kumar Kolluru,
Mohammad Atif,
Santosh Ansumali
Abstract:
Kinetic models of polyatomic gas typically account for the internal degrees of freedom at the level of the two-particle distribution function. However, close to the hydrodynamic limit, the internal (rotational) degrees of freedom tend to be well represented just by rotational kinetic energy density. We account for the rotational energy by augmenting the Ellipsoidal-statistical BGK (ES-BGK) model,…
▽ More
Kinetic models of polyatomic gas typically account for the internal degrees of freedom at the level of the two-particle distribution function. However, close to the hydrodynamic limit, the internal (rotational) degrees of freedom tend to be well represented just by rotational kinetic energy density. We account for the rotational energy by augmenting the Ellipsoidal-statistical BGK (ES-BGK) model, an extension of the Bhatnagar-Gross- Krook (BGK) model, at the level of the single-particle distribution function with an advection-diffusion-relaxation equation for the rotational energy. This reduced model respects the H theorem and recovers the compressible hydrodynamics for polyatomic gases as its macroscopic limit. As required for a polyatomic gas model, this extension of the ES-BGK model has not only correct specific heat ratio but also allows for three independent tunable transport coefficients: thermal conductivity, shear viscosity, and bulk viscosity. We illustrate the effectiveness of the model via a lattice Boltzmann method implementation.
△ Less
Submitted 13 January, 2022;
originally announced January 2022.
-
Estimating Hidden Asymptomatics, Herd Immunity Threshold and Lockdown Effects using a COVID-19 Specific Model
Authors:
Shaurya Kaushal,
Abhineet Singh Rajput,
Soumyadeep Bhattacharya,
M. Vidyasagar,
Aloke Kumar,
Meher K. Prakash,
Santosh Ansumali
Abstract:
A quantitative COVID-19 model that incorporates hidden asymptomatic patients is developed, and an analytic solution in parametric form is given. The model incorporates the impact of lockdown and resulting spatial migration of population due to announcement of lockdown. A method is presented for estimating the model parameters from real-world data. It is shown that increase of infections slows down…
▽ More
A quantitative COVID-19 model that incorporates hidden asymptomatic patients is developed, and an analytic solution in parametric form is given. The model incorporates the impact of lockdown and resulting spatial migration of population due to announcement of lockdown. A method is presented for estimating the model parameters from real-world data. It is shown that increase of infections slows down and herd immunity is achieved when symptomatic patients are 4-6\% of the population for the European countries we studied, when the total infected fraction is between 50-56 \%. Finally, a method for estimating the number of asymptomatic patients, who have been the key hidden link in the spread of the infections, is presented.
△ Less
Submitted 29 May, 2020;
originally announced June 2020.
-
Fokker-Planck model for binary mixtures
Authors:
Samarth Agrawal,
S. K. Singh,
Santosh Ansumali
Abstract:
In dilute gas kinetic theory, model collision dynamics such as Bhatnagar-Gross-Krook (BGK) model is often used to get a better insight and numerical modelling. BGK model and its variants assume that highly nonlinear collision term can be replaced by a simple relaxation dynamics towards Maxwell-Boltzmann distribution. Lebowitz et al., proposed an alternative framework for the collision model, known…
▽ More
In dilute gas kinetic theory, model collision dynamics such as Bhatnagar-Gross-Krook (BGK) model is often used to get a better insight and numerical modelling. BGK model and its variants assume that highly nonlinear collision term can be replaced by a simple relaxation dynamics towards Maxwell-Boltzmann distribution. Lebowitz et al., proposed an alternative framework for the collision model, known as Fokker-Planck model, where the relaxation of an arbitrary distribution towards Maxwell-Boltzmann distribution is modelled as a drift-diffusion process in velocity space. In the present manuscript, we extend the single component Fokker-Planck model to a binary gas mixture model over a large range of Schmidt numbers. We prove that these mixture models satisfy the necessary conservation laws and the H-theorem.
△ Less
Submitted 23 October, 2019;
originally announced October 2019.
-
Lattice Boltzmann model for weakly compressible flows
Authors:
Praveen Kumar Kolluru,
Mohammad Atif,
Manjusha Namburi,
Santosh Ansumali
Abstract:
We present an energy conserving lattice Boltzmann model based on a crystallographic lattice for simulation of weakly compressible flows. The theoretical requirements and the methodology to construct such a model are discussed. We demonstrate that the model recovers the isentropic sound speed in addition to the effects of viscous heating and heat flux dynamics. Several test cases for acoustics, the…
▽ More
We present an energy conserving lattice Boltzmann model based on a crystallographic lattice for simulation of weakly compressible flows. The theoretical requirements and the methodology to construct such a model are discussed. We demonstrate that the model recovers the isentropic sound speed in addition to the effects of viscous heating and heat flux dynamics. Several test cases for acoustics, thermal and thermoacoustic flows are simulated to show the accuracy of the proposed model.
△ Less
Submitted 15 September, 2019;
originally announced September 2019.
-
Lattice operators from discrete hydrodynamics
Authors:
Rashmi Ramadugu,
Sumesh P. Thampi,
Ronojoy Adhikari,
Sauro Succi,
Santosh Ansumali
Abstract:
We present a general scheme to derive lattice differential operators from the discrete velocities and associated Maxwell-Boltzmann distributions used in lattice hydrodynamics. Such discretizations offer built-in isotropy and recursive techniques to increase the convergence order. This provides a simple and elegant procedure to derive isotropic and accurate discretizations of differential operators…
▽ More
We present a general scheme to derive lattice differential operators from the discrete velocities and associated Maxwell-Boltzmann distributions used in lattice hydrodynamics. Such discretizations offer built-in isotropy and recursive techniques to increase the convergence order. This provides a simple and elegant procedure to derive isotropic and accurate discretizations of differential operators, which are expected to apply across a broad range of problems in computational physics.
△ Less
Submitted 5 August, 2012;
originally announced August 2012.
-
Isotropic discrete Laplacian operators from lattice hydrodynamics
Authors:
Sumesh P. Thampi,
Santosh Ansumali,
Ronojoy Adhikari,
Sauro Succi
Abstract:
We show that discrete schemes developed for lattice hydrodynamics provide an elegant and physically transparent way of deriving Laplacians with isotropic discretisation error. Isotropy is guaranteed whenever the Laplacian weights follow from the discrete Maxwell-Boltzmann equilibrium since these are, by construction, isotropic on the lattice. We also point out that stencils using as few as 15 poin…
▽ More
We show that discrete schemes developed for lattice hydrodynamics provide an elegant and physically transparent way of deriving Laplacians with isotropic discretisation error. Isotropy is guaranteed whenever the Laplacian weights follow from the discrete Maxwell-Boltzmann equilibrium since these are, by construction, isotropic on the lattice. We also point out that stencils using as few as 15 points in three dimensions, generate isotropic Laplacians. These computationally efficient Laplacians can be used in cell-dynamical and hybrid lattice Boltzmann simulations, in favor of popular anisotropic Laplacians, which make use of larger stencils. The method can be extended to provide discretisations of higher order and for other differential operators, such the gradient, divergence and curl.
△ Less
Submitted 15 February, 2012;
originally announced February 2012.
-
Hydrodynamics beyond Navier-Stokes: Exact solution to the lattice Boltzmann hierarchy
Authors:
S. Ansumali,
I. V. Karlin,
S. Arcidiacono,
A. Abbas,
N. I. Prasianakis
Abstract:
Exact solution to the hierarchy of nonlinear lattice Boltzmann (LB) kinetic equations in the stationary planar Couette flow is found at non-vanishing Knudsen numbers. A new method of solving LB kinetic equations which combines the method of moments with boundary conditions for populations enables to derive closed-form solutions for all higher-order moments. Convergence of results suggests that t…
▽ More
Exact solution to the hierarchy of nonlinear lattice Boltzmann (LB) kinetic equations in the stationary planar Couette flow is found at non-vanishing Knudsen numbers. A new method of solving LB kinetic equations which combines the method of moments with boundary conditions for populations enables to derive closed-form solutions for all higher-order moments. Convergence of results suggests that the LB hierarchy with larger velocity sets is the novel way to approximate kinetic theory.
△ Less
Submitted 16 March, 2007;
originally announced March 2007.
-
Renormalization of the Lattice Boltzmann Hierarchy
Authors:
I. V. Karlin,
S. Ansumali
Abstract:
Is it possible to solve Boltzmann-type kinetic equations using only a small number of particles velocities? We introduce a novel techniques of solving kinetic equations with (arbitrarily) large number of particle velocities using only a lattice Boltzmann method on standard, low-symmetry lattices. The renormalized kinetic equation is validated with the exact solution of the planar Couette flow at…
▽ More
Is it possible to solve Boltzmann-type kinetic equations using only a small number of particles velocities? We introduce a novel techniques of solving kinetic equations with (arbitrarily) large number of particle velocities using only a lattice Boltzmann method on standard, low-symmetry lattices. The renormalized kinetic equation is validated with the exact solution of the planar Couette flow at moderate Knudsen numbers.
△ Less
Submitted 16 March, 2007;
originally announced March 2007.
-
Thermodynamic Theory of Incompressible Hydrodynamics
Authors:
Santosh Ansumali,
Iliya V. Karlin,
Hans Christian Öttinger
Abstract:
The grand potential for open systems describes thermodynamics of fluid flows at low Mach numbers. A new system of reduced equations for the grand potential and the fluid momentum is derived from the compressible Navier-Stokes equations. The incompressible Navier-Stokes equations are the quasi-stationary solution to the new system. It is argued that the grand canonical ensemble is the unifying co…
▽ More
The grand potential for open systems describes thermodynamics of fluid flows at low Mach numbers. A new system of reduced equations for the grand potential and the fluid momentum is derived from the compressible Navier-Stokes equations. The incompressible Navier-Stokes equations are the quasi-stationary solution to the new system. It is argued that the grand canonical ensemble is the unifying concept for the derivation of models and numerical methods for incompressible fluids, illustrated here with a simulation of a minimal Boltzmann model in a microflow setup.
△ Less
Submitted 2 February, 2005; v1 submitted 24 May, 2004;
originally announced May 2004.
-
Kinetic Theory of Turbulence Modeling: Smallness Parameter, Scaling and Microscopic Derivation of Smagorinsky Model
Authors:
Santosh Ansumali,
Iliya V. Karlin,
Sauro Succi
Abstract:
A mean-field approach (filtering out subgrid scales) is applied to the Boltzmann equation in order to derive a subgrid turbulence model based on kinetic theory. It is demonstrated that the only Smagorinsky type model which survives in the hydrodynamic limit on the viscosity time scale is the so-called tensor-diffusivity model. Scaling of the filter-width with Reynolds number and Knudsen number i…
▽ More
A mean-field approach (filtering out subgrid scales) is applied to the Boltzmann equation in order to derive a subgrid turbulence model based on kinetic theory. It is demonstrated that the only Smagorinsky type model which survives in the hydrodynamic limit on the viscosity time scale is the so-called tensor-diffusivity model. Scaling of the filter-width with Reynolds number and Knudsen number is established. This sets the first rigorous step in deriving turbulence models from kinetic theory.
△ Less
Submitted 6 February, 2004; v1 submitted 27 October, 2003;
originally announced October 2003.