Grid refinement in LBM based on continuous distribution functions
|
|
- Monica Emerald Williamson
- 5 years ago
- Views:
Transcription
1 Grid refement LBM based on contuous distribution functions Denis Ricot 1 Simon Marié 1,2 Pierre Sagaut 2 1 Renault - Research, Material and Advanced Engeerg Department 2 d Alembert Institute, Université Paris VI 5th ICMMES june 2008
2 Grid refement problem For simulatg complex engeerg flows, mesh refement technique is necessary LB distribution functions g α are not conserved through a refement terface [Filippova & H -anel, 1998] Spatial and time terpolation methods can not be used directly on the distribution functions to evaluate the unknown data
3 Grid refement schemes Vertex centered approach with ([Yu et al., 2002][Dupuis & Chopard, 2003]) our without rescalg ([L & Lai, 2000]) rescalg and terpolations applied on the post-collision distribution functions ([Yu et al., 2002][Filippova & H -anel, 1998]) spatial terpolation performed before time terpolation [Yu et al., 2002][Dupuis & Chopard, 2003] Volumetric approach (cell centered) [Rohde, 2004][Chen et al., 2005 (PowerFLOW)] Explode (c f ) and coalesce (f c) distribution functions : mass conservative scheme No rescalg of distribution functions
4 Presentation outle Contuous distribution functions? Grid refement problem DVBE and LBE Numerical verification Rescalg procedure Overlappg nodes Algorithm description Grid parameters Pulse propagation uniform flow Convected vortex
5 Grid refement problem DVBE and LBE Numerical verification Grid refement problem After one standard LB timestep, some distribution functions are missg Distribution functions are known to be discontuous through an terface between meshes with different grid size : gα c (,, t) gα f (,, t)
6 Grid refement problem DVBE and LBE Numerical verification Base ideas of Filippova & H -anel (1998) Macroscopic properties of the fluid must be conserved (ν, c s) δx c /δt c = δx f /δt f τg c = 1 τ f 2 g + 1 (the relaxation time is not contuous) 4 Macroscopic variables of the flow must be conserved (ρ, u) g eq,c α = g eq,f α The shear stress must be contuous g neq,c α = m τ g c g neq,f τg f α (m = δxc = 2) δx f
7 Grid refement problem DVBE and LBE Numerical verification Why the distribution functions of LBM are not contuous? Boltzmann equation f t + c f i x i = f f eq ɛλ (BE) Quadrature formulae to calculate moments of distribution functions usg a discrete velocity set f α t + c α,i f α x i = 1 τ `fα fα eq (DVBE) By defition the distribution functions of DVBE f α are contuous space and time (and also ν = τθ 0, c s = θ 0, ρ = P f α, ρu = P c αf α)
8 Grid refement problem DVBE and LBE Numerical verification Integrate DVBE along the characteristic c α for a time terval δt The tegral of the BGK collision operator is approximated by the trapezium rule : f α (x + c αδt, t + δt) f α (x, t) = δt {fα (x + cαδt, t + δt) f eq α 2τ (x + c αδt, t + δt)+f α (x, t) fα eq Change of variable ([He et al., 1998][Dellar, 2001]) g α (x, t) =f α (x, t)+ δt 2τ `fα (x, t) fα eq (x, t) (x, t) + O δt 3 Fal Lattice Boltzmann Equation g α(x+c αδt, t +δt)= 1 δt «g α(x, t)+ δt fα eq (x, t) (2) τ g τ g (1) with τ g = τ + δt/2.
9 Grid refement problem DVBE and LBE Numerical verification Numerical verification of the lk between f α and g α LBE and DVBE computations have been performed on the same uniform grid, for the same flow problems and fluid parameters LB model : standard D2Q9, BGK collision operator, equilibrium function truncated at second order DVBE model : D2Q9 velocity model, BGK collision operator, equilibrium function truncated at second order 6th-order central fite difference scheme for space derivatives 5th-order Runge-Kutta scheme for time tegration Fluid is air, τ/δt = 0.93, same δt and δx for LBE and DVBE For two simple flows, comparison of f α, g α and d α = f α + δt eq (fα f 2τ α )
10 Grid refement problem DVBE and LBE Numerical verification Reference cases Acoustic pressure pulse an uniform flow 8 h >< ρ = 1 + a P exp `x1 ln 2 x 0 2 b P 1 + `x2 x 2 i 2 0 u 1 = U 0 >: u 2 = 0 Convected vortex 8 >< >: with a P = 0.03 and b P = 5 ρ = 1 h u 1 = U 0 + a T U 0 `x2 x 2 0 exp `x1 ln 2 b T x `x2 x 2 i 2 0 h u 2 = a T U 0 `x1 x 0 1 exp `x1 ln 2 b T x `x2 x 2 i 0 2 with a T = 0.5 etb T = 25
11 Grid refement problem DVBE and LBE Numerical verification Comparison of the macroscopic variables Density profile of the pressure pulse at time t = 50 along the le x 2 = x Velocity profile u 2 of the convected vortex 0.01 at time t = 100 along the le x 2 = x x/δx LBM; DVBE. x/δx The simulated results are exactly the same term of macroscopic variables gα eq = fα eq
12 Grid refement problem DVBE and LBE Numerical verification Distribution functions at a given pot as a function of time Pressure pulse, α = 6 Convected vortex, α = t/δt ( ):g α ( ):f α ( ):d α = f α + δt(f α f eq α )/2τ t/δt ( ):gα eq ( ):fα eq This numerical test confirms that g α = f α + δt eq (fα f 2τ α )
13 Rescalg procedure Overlappg nodes Algorithm description Conversion g α f α for a given x and t g α = f α + δt 2τ `fα fα eq <=> fα = 2τ 2τ + δt g α + δt «2τ f α eq By defition f α, τ, ρ = P f α, ρu = P c αf α ( f eq α ) are contuous 8 >< >: fα c = 2τ g c 2τ+δt c α + δtc f eq 2τ α fα f = 2τ g f 2τ+δt f α + δtf f eq 2τ α fα f = fα c ==> τ = τg f δtf = τ c 2 g δtc g c 2 α = 1 2 τ g f m = 2 gα f = 1 f ``2 τ 2 τ g g f gc +1 α + fα eq f ``2 τ g + 1 gα f fα eq Remark: these expressions immediately imply that g neq,c α = 2 τ g c g neq,f τ g f α
14 Rescalg procedure Overlappg nodes Algorithm description First approach (implemented with success our 3D LB code (L-BEAM)) Conversion gα c f α for needed terface pots Spatial and temporal terpolations on f α Conversion f α gα f Second approach (preferred this work) Conversion gα c gα f for needed terface pots Spatial and temporal terpolations on g α f In both cases Conversion step needs the equilibrium functions to be known
15 Rescalg procedure Overlappg nodes Algorithm description Fe grid pots must also be calculated as coarse grid pots Conversion g f α out (, t) g c α out (, t)
16 Rescalg procedure Overlappg nodes Algorithm description Conversion g c α `, t + 2δt f g f α `, t + 2δt f can be done because macroscopic variables can be now calculated at pots
17 Rescalg procedure Overlappg nodes Algorithm description Timestep t : all functions g c α and gf α are known everywhere
18 Rescalg procedure Overlappg nodes Algorithm description Timestep t : all functions g c α and gf α are known everywhere Stage 1 : convert g f α out (, t) g c α out (, t)
19 Rescalg procedure Overlappg nodes Algorithm description Timestep t : all functions g c α and gf α are known everywhere Stage 1 : convert g f α out (, t) g c α out (, t) Stage 2 : collide and propagate all pots the fe and coarse regions (cludg pots )
20 Rescalg procedure Overlappg nodes Algorithm description Timestep t : all functions g c α and gf α are known everywhere Stage 1 : convert g f α out (, t) g c α out (, t) Stage 2 : collide and propagate all pots the fe and coarse regions (cludg pots ) Stage 3 : convert g c α, t +2δt f g f α, t +2δt f
21 Rescalg procedure Overlappg nodes Algorithm description Injection scheme : gα f,t +δt f = g f α,t +2δt f Lear terpolation usg gα,t f +2δt f and g f α (,t) (that must be stored) Timestep t : all functions g c α and gf α are known everywhere Stage 1 : convert g f α out (, t) g c α out (, t) Stage 2 : collide and propagate all pots the fe and coarse regions (cludg pots ) Stage 3 : convert g c α, t +2δt f g f α, t +2δt f Stage 4 : time terpolation of g f α,t +δt f
22 Rescalg procedure Overlappg nodes Algorithm description Injection scheme : gα f,t +δt f = g f α,t +2δt f Lear terpolation usg gα,t f +2δt f and g f α (,t) (that must be stored) Stage 5 : spatial terpolation of gα f,t +δt f usg g f α 1,t +δt f and gα f + 1,t +δt f Timestep t : all functions g c α and gf α are known everywhere Stage 1 : convert g f α out (, t) g c α out (, t) Stage 2 : collide and propagate all pots the fe and coarse regions (cludg pots ) Stage 3 : convert g c α, t +2δt f g f α, t +2δt f Stage 4 : time terpolation of g f α,t +δt f
23 Rescalg procedure Overlappg nodes Algorithm description Injection scheme : gα f,t +δt f = g f α,t +2δt f Lear terpolation usg gα,t f +2δt f and g f α (,t) (that must be stored) Stage 5 : spatial terpolation of gα f,t +δt f usg g f α 1,t +δt f and gα f + 1,t +δt f Timestep t+δt f : all functions gα f are known the fe grid region Timestep t : all functions g c α and gf α are known everywhere Stage 1 : convert g f α out (, t) g c α out (, t) Stage 2 : collide and propagate all pots the fe and coarse regions (cludg pots ) Stage 3 : convert g c α, t +2δt f g f α, t +2δt f Stage 4 : time terpolation of g f α,t +δt f
24 Rescalg procedure Overlappg nodes Algorithm description Injection scheme : gα f,t +δt f = g f α,t +2δt f Lear terpolation usg gα,t f +2δt f and g f α (,t) (that must be stored) Stage 5 : spatial terpolation of gα f,t +δt f usg g f α 1,t +δt f and gα f + 1,t +δt f Timestep t : all functions g c α and gf α are known everywhere Stage 1 : convert g f α out (, t) g c α out (, t) Stage 2 : collide and propagate all pots the fe and coarse regions (cludg pots ) Stage 3 : convert g c α, t +2δt f g f α, t +2δt f Stage 4 : time terpolation of g f α,t +δt f Timestep t+δt f : all functions gα f are known the fe grid region Stage 6 : collide and propagate all pots the fe region
25 Rescalg procedure Overlappg nodes Algorithm description Injection scheme : gα f,t +δt f = g f α,t +2δt f Lear terpolation usg gα,t f +2δt f and g f α (,t) (that must be stored) Stage 5 : spatial terpolation of gα f,t +δt f usg g f α 1,t +δt f and gα f + 1,t +δt f Timestep t : all functions g c α and gf α are known everywhere Stage 1 : convert g f α out (, t) g c α out (, t) Stage 2 : collide and propagate all pots the fe and coarse regions (cludg pots ) Stage 3 : convert g c α, t +2δt f g f α, t +2δt f Stage 4 : time terpolation of g f α,t +δt f Timestep t+δt f : all functions gα f are known the fe grid region Stage 6 : collide and propagate all pots the fe region Stage 7 : recover g f α,t +2δt f from g c α,t +2δt f
26 Rescalg procedure Overlappg nodes Algorithm description Injection scheme : gα f,t +δt f = g f α,t +2δt f Lear terpolation usg gα,t f +2δt f and g f α (,t) (that must be stored) Stage 5 : spatial terpolation of gα f,t +δt f usg g f α 1,t +δt f and gα f + 1,t +δt f Timestep t : all functions g c α and gf α are known everywhere Stage 1 : convert g f α out (, t) g c α out (, t) Stage 2 : collide and propagate all pots the fe and coarse regions (cludg pots ) Stage 3 : convert g c α, t +2δt f g f α, t +2δt f Stage 4 : time terpolation of g f α,t +δt f Timestep t+δt f : all functions gα f are known the fe grid region Stage 6 : collide and propagate all pots the fe region Stage 7 : recover g f α,t +2δt f from g c α,t +2δt f Stage 8 : spatial terpolation of gα f,t +2δt f
27 Rescalg procedure Overlappg nodes Algorithm description Injection scheme : gα f,t +δt f = g f α,t +2δt f Lear terpolation usg gα,t f +2δt f and g f α (,t) (that must be stored) Stage 5 : spatial terpolation of gα f,t +δt f usg g f α 1,t +δt f and gα f + 1,t +δt f Timestep t : all functions g c α and gf α are known everywhere Stage 1 : convert g f α out (, t) g c α out (, t) Stage 2 : collide and propagate all pots the fe and coarse regions (cludg pots ) Stage 3 : convert g c α, t +2δt f g f α, t +2δt f Stage 4 : time terpolation of g f α,t +δt f Timestep t+δt f : all functions gα f are known the fe grid region Stage 6 : collide and propagate all pots the fe region Stage 7 : recover g f α,t +2δt f from g c α,t +2δt f Stage 8 : spatial terpolation of gα f,t +2δt f Timestep t +2δt f : all functions g c α and gf α are known everywhere
28 Grid parameters Pulse propagation uniform flow Convected vortex Coarse grid : 120x80, fe grid : Acoustic pulse or vortex is itialized outside or side the fe grid M = 0.2, ν = m 2 /s, c s = 340 m/s, δx f = 10 2 m All computations are also done with a uniform coarse grid (reference) 80 U 0 = y /δ x c x/ δ x c
29 Grid parameters Pulse propagation uniform flow Convected vortex Iso-contours of the density fluctuation 80 Uniform grid Injection (no time terpolation) Lear terpolation time y/ δ x c y /δ x c y /δ x c x/ δ x c x/ δ x c x/ δ x c Small pressure reflexion occurs at the grid terface. The error is reduced when lear time terpolation is used. There is an analytical solution for the propagation of the pulse uniform flow precise error quantification is possible
30 Grid parameters Pulse propagation uniform flow Convected vortex Simulation with the uniform grid 5 x ( ) reference analytical solution ( ) LB simulation ρ / ρ x/ δ x c
31 Grid parameters Pulse propagation uniform flow Convected vortex Two-grid simulation with jection (no time terpolation) 5 x ( ) reference analytical solution ( ) results on coarse pots ρ / ρ ( ) results on fe pots x/ δ x c
32 Grid parameters Pulse propagation uniform flow Convected vortex Two-grid simulation with lear terpolation time 5 x ( ) reference analytical solution ( ) results on coarse pots ρ / ρ ( ) results on fe pots x/ δ x c
33 Grid parameters Pulse propagation uniform flow Convected vortex Convergence rate Parameter b p gives the itial width of the pressure pulse, i.e. the spatial resolution 10 6 ( ) slope ( ) Uniform grid ( ) Two-grid simulation without time terpolation (jection) L ( ) Two-grid simulation with lear terpolation time /2 b p
34 Grid parameters Pulse propagation uniform flow Convected vortex Initial pulse the fe grid region 80 Uniform grid Injection (no time terpolation) Lear terpolation time y/ δ x c y /δ x c y /δ x c x/ δ x c x/ δ x c x/ δ x c
35 Grid parameters Pulse propagation uniform flow Convected vortex Convected vortex (lear terpolation time) Uniform grid Lear terpolation time Lear terpolation time y/ δ x c y/ δ x c y/ δ x c x/ δ x c x/ δ x c x/ δ x c Fe to coarse grid convection : spurious vorticity appears side the fe region due to non-physical reflexion at the grid terface
36 Conclusion New theoretical sight the lk between coarse and fe distribution functions fal rescalg expressions are fully equivalent to the previous approaches this new theoretical approach can be useful for other LB models (other collision operator such as MRT models) The proposed grid refement algorithm mimize the unknown functions that must be approximated with terpolations use of overlappg coarse nodes terpolations are done on the distribution functions time terpolation is done before spatial terpolation Accuracy of the grid refement scheme can be improved usg higher order terpolation schemes
External and Internal Incompressible Viscous Flows Computation using Taylor Series Expansion and Least Square based Lattice Boltzmann Method
Available online at http://ijim.srbiau.ac.ir/ Int. J. Industrial Mathematics (ISSN 2008-5621) Vol. 10, No. 2, 2018 Article ID IJIM-00726, 8 pages Research Article External and Internal Incompressible Viscous
More informationLATTICE BOLTZMANN MODELLING OF PULSATILE FLOW USING MOMENT BOUNDARY CONDITIONS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 5 June 28, Glasgow, UK LATTICE BOLTZMANN MODELLING OF PULSATILE FLOW USING MOMENT
More informationLattice Boltzmann Method for Moving Boundaries
Lattice Boltzmann Method for Moving Boundaries Hans Groot March 18, 2009 Outline 1 Introduction 2 Moving Boundary Conditions 3 Cylinder in Transient Couette Flow 4 Collision-Advection Process for Moving
More informationME615 Project Presentation Aeroacoustic Simulations using Lattice Boltzmann Method
ME615 Project Presentation Aeroacoustic Simulations using Lattice Boltzmann Method Kameswararao Anupindi Graduate Research Assistant School of Mechanical Engineering Purdue Universit December 11, 11 Outline...
More informationOn the lattice Boltzmann method for multiphase flows
Center for Turbulence Research Annual Research Briefs 2009 377 On the lattice Boltzmann method for multiphase flows By S. H. Kim AND H. Pitsch 1. Motivation and objectives The lattice Boltzmann (LB) method
More informationarxiv: v1 [physics.flu-dyn] 5 Oct 2017
Hermite regularization of the Lattice Boltzmann Method for open source computational aeroacoustics arxiv:1710.02065v1 [physics.flu-dyn] 5 Oct 2017 F. Brogi a),b) Department of Earth Science, University
More informationLattice Boltzmann Method for Fluid Simulations
1 / 16 Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas Simon Fraser University April 7, 2011 2 / 16 Ludwig Boltzmann and His Kinetic Theory of Gases The Boltzmann Transport
More informationComputers and Mathematics with Applications
Computers and Mathematics with Applications 59 (2010 2194 2199 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Toward
More informationDIRECT SIMULATION OF TRAILING-EDGE NOISE GENERATED BY A CONTROLLED DIFFUSION AIRFOIL USING A LATTICE-BOLTZMANN METHOD
DIRECT SIMULATION OF TRAILING-EDGE NOISE GENERATED BY A CONTROLLED DIFFUSION AIRFOIL USING A LATTICE-BOLTZMANN METHOD M. Sanjosé, S. Moreau Department of Mechanical Engineering Université de Sherbrooke
More informationPredictor-Corrector Finite-Difference Lattice Boltzmann Schemes
Applied Mathematical Sciences, Vol. 9, 2015, no. 84, 4191-4199 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.5138 Predictor-Corrector Finite-Difference Lattice Boltzmann Schemes G. V.
More informationEquivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme
Equivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme Pietro Asinari (1), Taku Ohwada (2) (1) Department of Energetics, Politecnico di Torino, Torino 10129,
More informationLattice Boltzmann Method
3 Lattice Boltzmann Method 3.1 Introduction The lattice Boltzmann method is a discrete computational method based upon the lattice gas automata - a simplified, fictitious molecular model. It consists of
More informationLattice Boltzmann model for the Elder problem
1549 Lattice Boltzmann model for the Elder problem D.T. Thorne a and M.C. Sukop a a Department of Earth Sciences, Florida International University, PC 344, University Park, 11200 SW 8th Street, Miami,
More informationarxiv: v1 [physics.flu-dyn] 26 May 2015
Increasing stability and accuracy of the lattice Boltzmann scheme: recursivity and regularization arxiv:1505.06900v1 [physics.flu-dyn] 26 May 2015 Orestis Malaspinas a,b a Institut Jean le Rond d Alembert,
More informationPOTENTIAL TURBULENCE MODEL PREDICTIONS OF FLOW PAST A TRIANGULAR CYLINDER USING AN UNSTRUCTURED STAGGERED MESH METHOD
POTENTIAL TURBULENCE MODEL PREDICTIONS OF FLOW PAST A TRIANGULAR CYLINDER USING AN UNSTRUCTURED STAGGERED MESH METHOD Xg Zhang Blair Perot Department of Mechanical and Industrial Engeerg, University of
More informationComputers & Fluids. Improving lattice Boltzmann simulation of moving particles in a viscous flow using local grid refinement
Computers & Fluids 00 (204) 2 Computers & Fluids Improving lattice Boltzmann simulation o moving particles in a viscous low using local grid reinement Songying Chen a,, Cheng Peng b, Yihua Teng c, Lian-Ping
More informationOn the stability of a relative velocity lattice Boltzmann scheme for compressible Navier-Stokes equations
On the stability of a relative velocity lattice Boltzmann scheme for compressible Navier-Stokes equations François Dubois, Tony Fevrier, Benjamin Graille To cite this version: François Dubois, Tony Fevrier,
More informationME615 Project Report Aeroacoustic Simulations using Lattice Boltzmann Method
ME615 Project Report Aeroacoustic Simulations using Lattice Boltzmann Method Kameswararao Anupindi School of Mechanical Engineering, Purdue Universit, West Lafaette, IN, 4796, USA Lattice Boltzmann method
More informationConnection Between the Lattice Boltzmann Equation and the Beam Scheme
NASA/CR-1999-209001 ICASE Report No. 99-10 Connection Between the Lattice Boltzmann Equation and the Beam Scheme Kun Xu Hong Kong University of Science and Technology, Kowloon, Hong Kong Li-Shi Luo ICASE,
More informationGeneralized Local Equilibrium in the Cascaded Lattice Boltzmann Method. Abstract
Accepted for publication on Physical Review E (R), code ERR1034 Generalized Local Equilibrium in the Cascaded Lattice Boltzmann Method Pietro Asinari Department of Energetics, Politecnico di Torino, Corso
More informationFree energy concept Free energy approach LBM implementation Parameters
BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH
More informationA hybrid method for hydrodynamic-kinetic flow - Part II - Coupling of hydrodynamic and kinetic models
A hybrid method for hydrodynamic-kinetic flow - Part II - Coupling of hydrodynamic and kinetic models Alessandro Alaia, Gabriella Puppo May 31, 2011 Abstract In this work we present a non stationary domain
More informationLattice Boltzmann approach to liquid - vapour separation
Lattice Boltzmann approach to liquid - vapour separation T.Biciușcă 1,2, A.Cristea 1, A.Lamura 3, G.Gonnella 4, V.Sofonea 1 1 Centre for Fundamental and Advanced Technical Research, Romanian Academy Bd.
More informationResearch of Micro-Rectangular-Channel Flow Based on Lattice Boltzmann Method
Research Journal of Applied Sciences, Engineering and Technology 6(14): 50-55, 013 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scientific Organization, 013 Submitted: November 08, 01 Accepted: December 8,
More informationCOMPARISON OF CPU AND GPU IMPLEMENTATIONS OF THE LATTICE BOLTZMANN METHOD
XVIII International Conference on Water Resources CMWR 2010 J. Carrera (Ed) c CIMNE, Barcelona, 2010 COMPARISON OF CPU AND GPU IMPLEMENTATIONS OF THE LATTICE BOLTZMANN METHOD James.E. McClure, Jan F. Prins
More informationLattice Boltzmann Modeling of Wave Propagation and Reflection in the Presence of Walls and Blocks
Lattice Boltzmann Modeling of Wave Propagation and Reflection in the Presence of Walls and Blocs Maysam Saidi, Hassan Basirat Tabrizi, Member, IAENG, and Reza Sepahi Samian Abstract Lattice Boltzmann method
More informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More informationApplication of the Lattice-Boltzmann method in flow acoustics
Application of the Lattice-Boltzmann method in flow acoustics Andreas Wilde Fraunhofer Institut für Integrierte Schaltungen, Außenstelle EAS Zeunerstr. 38, 01069 Dresden Abstract The Lattice-Boltzmann
More informationLattice Boltzmann Method Solver Documentation
Lattice Boltzmann Method Solver Documentation Release 0.0.1 Sayop Kim May 03, 2016 Contents 1 Code Instruction 3 1.1 Quick instruction for running the simulation.............................. 3 2 Results
More informationLevel Set-based Topology Optimization Method for Viscous Flow Using Lattice Boltzmann Method
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level Set-based Topology Optimization Method for Viscous Flow Using Lattice Boltzmann Method
More informationFINITE-DIFFERENCE IMPLEMENTATION OF LATTICE BOLTZMANN METHOD FOR USE WITH NON-UNIFORM GRIDS
7. ANKARA INTERNATIONAL AEROSPACE CONFERENCE AIAC-2013-143 11-13 September 2013 - METU, Ankara TURKEY FINITE-DIFFERENCE IMPLEMENTATION OF LATTICE BOLTZMANN METHOD FOR USE WITH NON-UNIFORM GRIDS Fatih ÇEVİK
More informationThe lattice Boltzmann method for contact line dynamics
The lattice Boltzmann method for contact line dynamics Sudhir Srivastava, J.H.M. ten Thije Boonkkamp, Federico Toschi April 13, 2011 Overview 1 Problem description 2 Huh and Scriven model 3 Lattice Boltzmann
More informationDiscontinuous shear thickening on dense non-brownian suspensions via lattice Boltzmann method
Discontinuous shear thickening on dense non-brownian suspensions via lattice Boltzmann method Pradipto and Hisao Hayakawa Yukawa Insitute for Theoretical Physics Kyoto University Rheology of disordered
More informationStudy on lattice Boltzmann method/ large eddy simulation and its application at high Reynolds number flow
Research Article Study on lattice Boltzmann method/ large eddy simulation and its application at high Reynolds number flow Advances in Mechanical Engineering 1 8 Ó The Author(s) 2015 DOI: 10.1177/1687814015573829
More informationAnalysis of the Lattice Boltzmann Method by the Truncated Moment System
by the Truncated Moment System Pietro Asinari, PhD Department of Energetics, Politecnico di Torino, Torino 10129, Italy, pietro.asinari@polito.it http://staff.polito.it/pietro.asinari Lab for Simulation,
More informationLattice Boltzmann Simulations of 2D Laminar Flows past Two Tandem Cylinders
Lattice Boltzmann Simulations of 2D Laminar Flows past Two Tandem Cylinders Alberto Mussa a, Pietro Asinari a,, and Li-Shi Luo b a Department of Energetics, Politecnico di Torino, Torino 10129, Italy Tel.:
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationTURBULENT VORTEX SHEDDING FROM TRIANGLE CYLINDER USING THE TURBULENT BODY FORCE POTENTIAL MODEL
Proceedgs of ASME FEDSM ASME 2 Fluids Engeerg Division Summer Meetg June 11-15, 2 Boston, Massachusetts FEDSM2-11172 TURBULENT VORTEX SHEDDING FROM TRIANGLE CYLINDER USING THE TURBULENT BODY FORCE POTENTIAL
More informationSemi-implicit-linearized Multiple-relaxation-time formulation of Lattice Boltzmann Schemes for Mixture Modeling. Abstract
Accepted for publication on Physical Review E Semi-implicit-linearized Multiple-relaxation-time formulation of Lattice Boltzmann Schemes for Mixture Modeling Pietro Asinari Department of Energetics, Politecnico
More informationAPPRAISAL OF FLOW SIMULATION BY THE LATTICE BOLTZMANN METHOD
APPRAISAL OF FLOW SIMULATION BY THE LATTICE BOLTZMANN METHOD Guillermo Izquierdo Bouldstridge Imperial College of London Department of Aeronautics Master s Thesis Supervisor: Dr. Joaquim Peiró September
More informationStructural stability of Lattice Boltzmann schemes
Structural stability of Lattice Boltzmann schemes Claire David, Pierre Sagaut To cite this version: Claire David, Pierre Sagaut. Structural stability of Lattice Boltzmann schemes. Physica A, Elsevier,
More informationProperties of the Lattice Boltzmann Method for Acoustics
15th AIAA/CEAS Aeroacoustics Conference (3th AIAA Aeroacoustics Conference) 11-13 May 29, Miami, Florida AIAA 29-3395 Properties of the Lattice Boltzmann Method for Acoustics Guillaume A. Brès, Franck
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationPhysical Modeling of Multiphase flow. Boltzmann method
with lattice Boltzmann method Exa Corp., Burlington, MA, USA Feburary, 2011 Scope Scope Re-examine the non-ideal gas model in [Shan & Chen, Phys. Rev. E, (1993)] from the perspective of kinetic theory
More informationThe lattice Boltzmann equation: background and boundary conditions
The lattice Boltzmann equation: background and boundary conditions Tim Reis Department of Mathematics Sciences University of Greenwich Leslie Comrie Seminar Series November 2016 Acknowledgements: P. Dellar
More informationThe lattice Boltzmann equation (LBE) has become an alternative method for solving various fluid dynamic
36th AIAA Fluid Dynamics Conference and Exhibit 5-8 June 2006, San Francisco, California AIAA 2006-3904 Direct and Large-Eddy Simulation of Decaying and Forced Isotropic Turbulence Using Lattice Boltzmann
More informationOn the Comparison Between Lattice Boltzmann Methods and Spectral Methods for DNS of Incompressible Turbulent Channel Flows on Small Domain Size
Advances in Applied Mathematics and Mechanics Adv. Appl. Math. Mech., Vol. 11, No. 3, pp. 598-67 DOI: 1.428/aamm.218.s4 June 219 On the Comparison Between Lattice Boltzmann Methods and Spectral Methods
More informationSEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS
The Seventh Nobeyama Workshop on CFD: To the Memory of Prof. Kuwahara Sanjo Kaikan, the University of Tokyo, Japan, September. 23-24, 2009 SEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS Jaw-Yen
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationSummary We develop an unconditionally stable explicit particle CFD scheme: Boltzmann Particle Hydrodynamics (BPH)
Summary Steps in BPH Space is divided into Cartesian cells Finite number of particles ~10 6-10 7 Particles fly freely between t n and t n+1 Mass, momentum and total energy are conserved Relax into a (LTE)
More informationSimulation of floating bodies with lattice Boltzmann
Simulation of floating bodies with lattice Boltzmann by Simon Bogner, 17.11.2011, Lehrstuhl für Systemsimulation, Friedrich-Alexander Universität Erlangen 1 Simulation of floating bodies with lattice Boltzmann
More informationProblem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 1600
Problem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 6 Overview This problem is aimed at testing the accuracy and the performance of high-order methods on the direct numerical simulation
More informationComputers and Mathematics with Applications. Investigation of the LES WALE turbulence model within the lattice Boltzmann framework
Computers and Mathematics with Applications 59 (2010) 2200 2214 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Investigation
More informationInteraction(s) fluide-structure & modélisation de la turbulence
Interaction(s) fluide-structure & modélisation de la turbulence Pierre Sagaut Institut Jean Le Rond d Alembert Université Pierre et Marie Curie- Paris 6, France http://www.ida.upmc.fr/~sagaut GDR Turbulence
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 19, 2013 http://acousticalsociety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Engineering Acoustics Session 1pEAa: Active and Passive Control of Fan
More informationNumerical Investigation of Fluid and Thermal Flow in a Differentially Heated Side Enclosure walls at Various Inclination Angles
Numerical Investigation of Fluid and Thermal Flow in a Differentially Heated Side Enclosure walls at Various Inclination Angles O. A. SHAHRUL Faculty of Mechanical & Manufacturing Engineering Universiti
More informationMulti species Lattice Boltzmann Models and Applications to Sustainable Energy Systems
Multi species Lattice Boltzmann Models and Applications to Sustainable Energy Systems Pietro Asinari, PhD Dipartimento di Energetica (DENER), Politecnico di Torino (POLITO), Torino 10129, Italy e-mail:
More informationGrad s approximation for missing data in lattice Boltzmann simulations
Europhysics Letters PREPRINT Grad s approximation for missing data in lattice Boltzmann simulations S. Ansumali 2, S. S. Chikatamarla 1 and I. V. Karlin 1 1 Aerothermochemistry and Combustion Systems Laboratory,
More informationDesigning Correct Fluid Hydrodynamics on A Rectangular Grid using MRT Lattice Boltzmann Approach
Computers & Fluids 00 04 3 Computers & Fluids Designing Correct Fluid Hydrodynamics on A Rectangular Grid using MRT Lattice Boltzmann Approach Yuan Zong a,b, Zhaoli Guo c, Lian-Ping Wang b,c, a State key
More informationOn the relation between lattice variables and physical quantities in lattice Boltzmann simulations
On the relation between lattice variables and physical quantities in lattice Boltzmann simulations Michael Junk Dirk Kehrwald th July 6 Kaiserslautern, Germany CONTENTS Contents Initial situation Problem
More informationNumerical Investigation on Shear Driven Cavity Flow by the Constrained Interpolated Profile Lattice Boltzmann Method
Numerical Investigation on Shear Driven Cavity Flow by the Constrained Interpolated Profile Lattice Boltzmann Method *C. S. NOR AZWADI, M. H. Al-MOLA and S. AGUS Department of Thermo-fluid Universiti Teknologi
More informationRound Cavity Noise Simulations using Lattice-Boltzmann Solver
Round Cavity Noise Simulations using Lattice-Boltzmann Solver Christophe Coreixas CERFACS, 42 Avenue G. Coriolis, 31507 Toulouse, France A new three-dimensional CFD solver, named LaBS, based on Lattice-Boltzmann
More informationLattice Boltzmann methods Summer term Cumulant-based LBM. 25. July 2017
Cumulant-based LBM Prof. B. Wohlmuth Markus Muhr 25. July 2017 1. Repetition: LBM and notation Discretization SRT-equation Content 2. MRT-model Basic idea Moments and moment-transformation 3. Cumulant-based
More informationNumerical Simulation of Flows about a Stationary and a Free-Falling Cylinder Using Immersed Boundary-Finite Difference Lattice Boltzmann Method
Numerical Simulation of Flows about a Stationary and a Free-Falling Cylinder Using Immersed Boundary-Finite Difference Lattice Boltzmann Method Roberto Rojas 1*, Kosuke Hayashi 1, Takeshi Seta and Akio
More informationA generic, mass conservative local grid renement technique for lattice-boltzmann schemes
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2006; 51:439 468 Published online 20 December 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/d.1140
More informationPREDICTION OF INTRINSIC PERMEABILITIES WITH LATTICE BOLTZMANN METHOD
PREDICTION OF INTRINSIC PERMEABILITIES WITH LATTICE BOLTZMANN METHOD Luís Orlando Emerich dos Santos emerich@lmpt.ufsc.br Carlos Enrique Pico Ortiz capico@lmpt.ufsc.br Henrique Cesar de Gaspari henrique@lmpt.ufsc.br
More informationarxiv:physics/ v1 [physics.comp-ph] 2 Nov 2006
Three-Dimensional Multi-Relaxation Time arxiv:physics/0611016v1 [physics.comp-ph] 2 Nov 2006 (MRT) Lattice-Boltzmann Models for Multiphase Flow Kannan N. Premnath, John Abraham M.J. Zucrow Labs., School
More informationZFS - Flow Solver AIA NVIDIA Workshop
ZFS - Flow Solver AIA NVIDIA Workshop Andreas Lintermann, Matthias Meinke, Wolfgang Schröder Institute of Aerodynamics and Chair of Fluid Mechanics RWTH Aachen University Outline Softwareproject ZFS (applications
More informationLattice Boltzmann method for adiabatic acoustics
369, 2371 2380 doi:10.1098/rsta.2011.0109 Lattice Boltzmann method for adiabatic acoustics BY YANBING LI AND XIAOWEN SHAN* Exa Corporation, Burlington, MA 01803, USA The lattice Boltzmann method (LBM)
More informationFIELD SCALE RESERVOIR SIMULATION THROUGH A LATTICE BOLTZMANN FRAMEWORK. A Dissertation ZACHARY ROSS BENAMRAM
FIELD SCALE RESERVOIR SIMULATION THROUGH A LATTICE BOLTZMANN FRAMEWORK A Dissertation by ZACHARY ROSS BENAMRAM Submitted to the Office of Graduate and Professional Studies of Texas A&M University in partial
More informationROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT
ROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT PETER C. LICHTNER 1, QINJUN KANG 1 1 Los Alamos National Laboratory,
More informationarxiv: v1 [physics.flu-dyn] 10 Aug 2015
Slip velocity of lattice Boltzmann simulation using bounce-back boundary scheme Jianping Meng, Xiao-Jun Gu, and David R Emerson Scientific Computing Department, STFC Daresbury laboratory, arxiv:1508.009v1
More informationSimulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang,
Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang, yfwang09@stanford.edu 1. Problem setting In this project, we present a benchmark simulation for segmented flows, which contain
More informationA Momentum Exchange-based Immersed Boundary-Lattice. Boltzmann Method for Fluid Structure Interaction
APCOM & ISCM -4 th December, 03, Singapore A Momentum Exchange-based Immersed Boundary-Lattice Boltzmann Method for Fluid Structure Interaction Jianfei Yang,,3, Zhengdao Wang,,3, and *Yuehong Qian,,3,4
More informationLattice Boltzmann simulations of a time-dependent natural convection problem
Lattice Boltzmann simulations of a time-dependent natural convection problem B. Trouette To cite this version: B. Trouette. Lattice Boltzmann simulations of a time-dependent natural convection problem.
More informationSimulation of 2D non-isothermal flows in slits using lattice Boltzmann method
Simulation of 2D non-isothermal flows in slits using lattice Boltzmann method João Chambel Leitão Department of Mechanical Engineering, Instituto Superior Técnico, Lisboa, Portugal Abstract: The present
More information[N175] Development of Combined CAA-CFD Algorithm for the Efficient Simulation of Aerodynamic Noise Generation and Propagation
The 32nd International Congress and Exposition on Noise Control Engineering Jeju International Convention Center, Seogwipo, Korea, August 25-28, 2003 [N175] Development of Combined CAA-CFD Algorithm for
More informationLattice Boltzmann simulations of a time-dependent natural convection problem
Lattice Boltzmann simulations of a time-dependent natural convection problem Benoît Trouette To cite this version: Benoît Trouette. Lattice Boltzmann simulations of a time-dependent natural convection
More informationAeroacoustic simulation of automotive ventilation outlets
Aeroacoustic simulation of automotive ventilation outlets J.-L. Adam a, D. Ricot a, F. Dubief a and C. Guy b a Renault SAS, 1 avenue du golf, 78288 Guyancourt, France b Ligeron, Les Algorithmes Bâtiment
More informationSolution of the Hirota Equation Using Lattice-Boltzmann and the Exponential Function Methods
Advanced Studies in Theoretical Physics Vol. 11, 2017, no. 7, 307-315 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/astp.2017.7418 Solution of the Hirota Equation Using Lattice-Boltzmann and the
More informationLattice Boltzmann Method
Chapter 5 Lattice Boltzmann Method 5.1 Pseudo-kinetic computational methods Understanding the flow properties of complex fluids in complex topologies is of importance to technological applications and
More informationAn adaptive lattice Boltzmann method for predicting wake fields behind wind turbines
An adaptive lattice Boltzmann method for predicting wake fields behind wind turbines Ralf Deiterding and Stephen L. Wood Abstract The crucial components of a dynamically adaptive, parallel lattice Boltzmann
More informationUsing LBM to Investigate the Effects of Solid-Porous Block in Channel
International Journal of Modern Physics and Applications Vol. 1, No. 3, 015, pp. 45-51 http://www.aiscience.org/journal/ijmpa Using LBM to Investigate the Effects of Solid-Porous Bloc in Channel Neda Janzadeh,
More informationA one-point boundary condition for the lattice Boltzmann method
manuscript A one-point boundary condition for the lattice Boltzmann method Michael Junk and Zhaoxia Yang FB Mathematik und Statistik, Universität Konstanz, Fach D194, 78457 Konstanz, Germany (Dated: September
More informationSimulation of Rarefied Gas Flow in Slip and Transitional Regimes by the Lattice Boltzmann Method
www.cfdl.issres.net Vol. 2 (2) June 2010 Simulation of Rarefied Gas Flow in Slip and Transitional Regimes by the Lattice Boltzmann Method N. Azwadi C. Sidik 1C, N.C. Horng 1, M.A. Mussa 2, and S. Abdullah
More informationLATTICE BOLTZMANN METHOD AND DIFFUSION IN MATERIALS WITH LARGE DIFFUSIVITY RATIOS
THERMAL SCIENCE: Year 27, Vol. 2, No. 3, pp. 73-82 73 LATTICE BOLTZMANN METHOD AND DIFFUSION IN MATERIALS WITH LARGE DIFFUSIVITY RATIOS by Edouard WALTHER *, Rachid BENNACER, and Caroline DE SA LMT, ENS
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationThermal lattice Bhatnagar-Gross-Krook model for flows with viscous heat dissipation in the incompressible limit
PHYSICAL REVIEW E 70, 066310 (2004) Thermal lattice Bhatnagar-Gross-Krook model for flows with viscous heat dissipation in the incompressible limit Yong Shi, T. S. Zhao,* Z. L. Guo Department of Mechanical
More informationIMPLEMENTING THE LATTICE-BOLTZMANN
IMPLEMENTING THE LATTICE-BOLTZMANN METHOD A RESEARCH ON BOUNDARY CONDITION TECHNIQUES by S.C. Wetstein in partial fulfillment of the requirements for the degree of Bachelor of Science in Applied Physics
More informationOn the Stability and Accuracy of the BGK, MRT and RLB Boltzmann Schemes for the Simulation of Turbulent Flows
Commun. Comput. Phys. doi: 10.4208/cicp.OA-2016-0229 Vol. 23, No. 3, pp. 846-876 March 2018 On the Stability and Accuracy of the BGK, MRT and RLB Boltzmann Schemes for the Simulation of Turbulent Flows
More informationA Framework for Hybrid Parallel Flow Simulations with a Trillion Cells in Complex Geometries
A Framework for Hybrid Parallel Flow Simulations with a Trillion Cells in Complex Geometries SC13, November 21 st 2013 Christian Godenschwager, Florian Schornbaum, Martin Bauer, Harald Köstler, Ulrich
More informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More informationCoupling of Lattice Boltzmann Equation and Finite Volume Method to Simulate Heat Transfer in a Square Cavity
Copyright 2009 Tech Science Press FDMP, vol.5, no.3, pp.283-295, 2009 Coupling of Lattice Boltzmann Equation and Finite Volume Method to Simulate Heat Transfer in a Square Cavity Ahmed Mezrhab 1 and Hassan
More informationPore-scale lattice Boltzmann simulation of laminar and turbulent flow through a sphere pack
Pore-scale lattice Boltzmann simulation of laminar and turbulent flow through a sphere pack Ehsan Fattahi a,b, Christian Waluga, Barbara Wohlmuth a, Ulrich Rüde b, Michael Manhart c, Rainer Helmig d a
More informationOutline. Unit Two Work and Paths. Unit Two Goals. Unit Two Goals. More Unit Two Goals. Analysis of Work. Unit 2 Work and paths September 2, 2010
Unit Work and paths Septemer, 00 Outle Unit Two Work and aths Mechanical Engeerg 70 Thermodynamics Larry Caretto Septemer, 00 Solution to quiz Assessment results Unit two work and paths Work, W = path
More informationEXTENDED FREE SURFACE FLOW MODEL BASED ON THE LATTICE BOLTZMANN APPROACH
METALLURGY AND FOUNDRY ENGINEERING Vol. 36, 2010, No. 2 Micha³ Szucki*, Józef S. Suchy***, Pawe³ ak*, Janusz Lelito**, Beata Gracz* EXTENDED FREE SURFACE FLOW MODEL BASED ON THE LATTICE BOLTZMANN APPROACH
More informationImproved treatment of the open boundary in the method of lattice Boltzmann equation
Improved treatment of the open boundary in the method of lattice Boltzmann equation Dazhi Yu, Renwei Mei and Wei Shyy Department of Mechanical and Aerospace Engineering, University of Florida, Gainesville,
More informationLATTICE BOLTZMANN AND FINITE VOLUME SIMULATIONS OF MULTIPHASE FLOW IN BGA ENCAPSULATION PROCESS
LATTICE BOLTZMANN AND FINITE VOLUME SIMULATIONS OF MULTIPHASE FLOW IN BGA ENCAPSULATION PROCESS Aizat Abas, MZ Abdullah, MHH Ishak, Nurfatin AS and Soon Fuat Khor Advanced Packaging and SMT, School of
More informationADVANCES IN RADIATION TRANSPORT MODELING USING LATTICE BOLTZMANN METHODS. Richard McCulloch. B.S., Utah State University, 2013 A THESIS
ADVANCES IN RADIATION TRANSPORT MODELING USING LATTICE BOLTZMANN METHODS by Richard McCulloch B.S., Utah State University, 2013 A THESIS submitted in partial fulfillment of the requirements for the degree
More informationSimplified Hyperbolic Moment Equations
Simplified Hyperbolic Moment Equations Julian Koellermeier and Manuel Torrilhon Abstract Hyperbolicity is a necessary property of model equations for the solution of the BGK equation to achieve stable
More information