A DYNAMICALLY ADAPTIVE LATTICE BOLTZMANN METHOD FOR THERMAL CONVECTION PROBLEMS
|
|
- Katherine Eaton
- 5 years ago
- Views:
Transcription
1 Int. J. Appl. Math. Comput. Sci.,, Vol., No., DOI: A DYNAMICALLY ADAPTIVE LATTICE BOLTZMANN METHOD FOR THERMAL CONVECTION PROBLEMS KAI FELDHUSEN a,b, RALF DEITERDING c,, CLAUS WAGNER a,b a Institute of Aerodynamics and Flow Technology German Aerospace Center (DLR), 3773 Göttingen, Germany b Institute of Thermo- and Fluiddynamics Technische Universität Ilmenau, Ilmenau, Germany c Aerodynamics and Flight Mechanics Research Group University of Southampton, Highfield Campus, Southampton SO7 BJ, United Kingdom r.deiterding@soton.ac.uk Utilizing the Boussinesq approximation, a double-population incompressible thermal lattice Boltzmann method (LBM) for forced and natural convection in two and three space dimensions is developed and validated. A block-structured dynamic adaptive mesh refinement (AMR) procedure tailored for LBM is applied to enable computationally efficient simulations of moderate to high Rayleigh number flows which are characterized by a large scale disparity in boundary layers and free stream flow. As test cases, the analytically accessible problem of a two-dimensional (2D) forced convection flow through two porous plates and the non-cartesian configuration of a heated rotating cylinder are considered. The objective of the latter is to advance the boundary conditions for accurate treatment of curved boundaries and to demonstrate the effect on the solution. The effectiveness of the overall approach is demonstrated for the natural convection benchmark of a 2D cavity with differentially heated walls at Rayleigh numbers from 3 up to 8. To demonstrate the benefit of the used AMR procedure for three-dimensional (3D) problems, results from the natural convection in a cubic cavity at Rayleigh numbers from 3 up to 5 are compared with benchmark results. Keywords: Lattice Boltzmann method, adaptive mesh refinement, thermal convection, incompressible. Introduction In recent years, the lattice Boltzmann method (LBM) has emerged as a powerful alternative to traditional Navier-Stokes (NS) solvers (Chen and Doolen, 998) to predict thermal fluid flow (Guo et al., 22; Kuznik et al., 27; Peng et al., 23), turbulent fluid flow (Jonas et al., 26), multiphase fluid Corresponding author flow (Lee and Lin, 25; Yu and Fan, 29) and magnetohydrodynamics (Deller, 22). Instead of discretizing the NS equations directly, the LBM is based on solving a simplified version of the Boltzmann equation in a specifically chosen discrete phase space. Using a Chapman-Enskog expansion, it has been shown that the approach recovers the NS equations in the limit of a vanishing Knudsen number (Hähnel, 24). Originally proposed for
2 2 K. FELDHUSEN et al. the isothermal weakly compressible case, several method enhancements for incompressibility (He and Luo, 997; Qian et al., 992) as well as incorporation of a buoyancy-driven temperature field for thermal convection flows are available (He et al., 998; Qian, 993). In general, there are two different categories of thermal lattice Boltzmann models. For the multispeed approach, the number of discrete velocity directions will be increased and the equilibrium distribution function is supplemented by higher order velocity terms to solve the internal energy equation, cf. (McNamara and Alder, 993; Alexander et al., 993; Qian, 993). However, this model is reported to exhibit numerical instabilities, cf. (Chen and Teixeira, 2). Here, we have chosen to pursue the strictly incompressible double distribution function (DDF) approach proposed by Guo et al. (22) for 2D and the straightforward expansion to 3D by He et al. (24) and Azwadi Che Sidik and Syahrullail (29). While the original LBM is formulated on a uniform Cartesian grid, an increase of local resolution is particularly necessary in the thermal boundary layers close to heated objects and walls. Kuznik et al. (27) and Peng et al. (23) demonstrated the computational benefit of a non-uniform grid for a thermal DDF LBM method in two and three spatial dimensions for simulating thermal convection in Cartesian cavities. In both works, a static geometry transformation is applied to the discretization in order to stretch the Cartesian lattice in the cavity center and reduce the spacing continuously towards the walls. Solution adaptive meshing is not used and on-the-fly mesh adaptation seems to have been applied so far to DDF LBM methods only in the context of isothermal two-phase flows, cf. (Yu and Fan, 29). Our objective in this paper is to close this gap. We supplement a thermal DDF LBM method with solution adaptive, dynamic mesh refinement. While adaptive lattice Boltzmann methods in the past have used primarily isotropic refinement of individual cells, cf. (Chen et al., 26), we apply in here a block-based approach, which is more suitable for the regular transport step of the LBM and thereby computationally significantly more efficient. The underlying data structures including distributed memory parallelization are borrowed from the finite volume mesh refinement system AMROC (Deiterding, 2). In order to fit smoothly into AMROC, the DDF LBM is formulated on cell-centered data structures and not node-based as it is mainly used for LBM in order to simplify the implementation of physical boundary conditions. In addition, complex geometry boundary condition treatment for possibly moving structures is incorporated. The update of the non-uniform lattice and the dynamic refinement procedure are orchestrated with the recursive Berger-Collela algorithm (Berger and Colella, 988). While the efficiency of this algorithm is undisputed for time-explicit finite volume schemes, its application to LBM is a novelty. In summary, our adaptive method is uniquely designed for the efficient simulation of real-world thermal flow problems. In this paper, the underlying computational techniques are described and the required validation for well-understood thermal convection problems is provided. In Section 2, we discuss the details of the numerical method, including the advanced thermal lattice Boltzmann approach, the block-based AMR method and the treatment of geometrically complex boundaries in the originally Cartesian scheme. Section 3 presents the computational results, where the analytic solution of the 2D flow between two moving porous plates, the 2D flow around a rotating heated cylinder and the well-known benchmark case of a two-dimensional cavity with differentially heated walls are considered. The result section is closed presenting the solution of the flow in a 3D cubic cavity with differentially heated walls. The conclusions including a short outlook are given in Section Numerical method 2.. Thermal lattice Boltzmann scheme. The incompressible two-dimensional LBM constructed under Boussinesq approximation used in the present work has been proposed by Guo et al. (22). For the three-dimensional case the incompressible LBM operator by
3 A dynamically adaptive lattice Boltzmann method for thermal convection problems 3 e 6 e 3 e 7 e 2 Fig.. Numerical stencil of D2Q9 - Discrete velocity directions in a computational cell. He et al. (24) is applied. By using the Bhatnagar-Gross-Krook (BGK) collision model (Bhatnagar et al., 954), the lattice Boltzmann equation for the partial probability distribution function f i with force field term F i can be formulated as e4 e 5 e e 8 f i (x + ce i t, t + t) = f i (x, t) ( ) f i (x, t) f (eq) i (x, t) + tf i. () τ ν In the DDF approach, a set of corresponding lattice Boltzmann equations g i (x + ce i t, t + t) = g i (x, t) ( ) g i (x, t) g (eq) i (x, t) τ D (2) is introduced based on distribution functions g i that are used to convect the macroscopic scalar quantity, here temperature, with the flow field. In the latter, e i is the unit velocity vector in direction of the ith discrete velocity space direction, t and t denote the time and time step, x the position, x the spatial increment, and c = x/ t is the particle speed. The relaxations times are τ ν for the flow field and τ D for the temperature field. functions are denoted by f (eq) i The respective equilibrium distribution and g (eq) i. In the two-dimensional case, a model with nine discrete unit velocities is used to compute the flow field (D2Q9) and an operator with four discrete velocities for the temperature field (D2Q4). The orientation of the discrete unit length velocities e i used to compute the velocity fields are depicted in Fig.. In the three-dimensional case, an operator with nineteen unit velocities is used for the flow field (D3Q9) and a model with six discrete velocities for the temperature field (D3Q6). The extended version of the orientation of the discrete unit length velocities e i are given in (3). (, ) i=, e i = (±,,),(,±,),(,,±) i=,...,6, (±,±,),(±,,±),(,±,±) i=7,...,8 (3) The basic LBM algorithm is divided into the steps of transport (or streaming) and collision, which are applied basically identically to () and (2). The following transport step represents the advection of fluid particles along the corresponding discrete velocities and is T : fi (x + ce i t, t + t) = f i (x, t). (4) Relaxation of the distribution functions towards the local equilibrium is performed on the transported distribution functions in the collision step C : f i (, t + t) = f i (, t + t) ( ) (eq) fi (, t) f i (, t). (5) τ ν With the pressure p and the velocity vector u as independent variables, the specific equilibrium distribution function f (eq) i for the D2Q9 model is defined as (Guo et al., 22) 4σ p f (eq) c s 2 i (u), for i =, i = λ p c + s 2 i (u), for i =,..., 4, γ p c + s 2 i (u), for i = 5,..., 8, (6) where the parameters σ, λ, and γ satisfy λ + γ = σ and λ + 2γ = /2. The functions s i (u) depend on the macroscopic velocity vector u and the discrete velocity vector e i and obey [ ] s i (u)=ω i 3 e i u +4.5 (e i u) 2 c c 2.5 u 2 c 2, (7) where the coefficients are given by ω = 4/9, ω,...,4 = /9, and ω 5,...,8 = /36. Using (6) and (7), the macroscopic values for velocity and dynamic pressure are given as u = [ ] ce i f i, p = c2 f i + s (u). 4σ i> i> (8)
4 4 K. FELDHUSEN et al. For the D3Q9 model the parameters change to σ = /2, λ = /8, and γ = /36. Furthermore the weight coefficients are given by ω = /3, ω,...,6 = /8, and ω 7,...,8 = /36. For the D2Q4 model used to compute the temperature field, the equilibrium function g (eq) i is g (eq) i = T [ e i u c ], for i =,..., 4 (9) and the macroscopic temperature is T = 4 g i. i= Analogously, in the D3Q6 model of the temperature field, the equilibrium function reads g (eq) i = T [ e i u c ], for i =,..., 6 () and the macroscopic temperature T = 6 g i. i= Since the fluid is assumed to be incompressible, a linear dependency between temperature differences and gravitational forces is applied (Boussinesq approximation), cf. (Mohamad and Kuzmin, 2), which leads to the force term F i. The force in () acts only in the two direct vertical directions. For 2D, this can be expressed according to Fig. (Guo et al., 22) as with F i = 2 (δ i2 + δ i4 ) e i F () F = gβ (T T ref ), (2) where g and β are the acceleration vector of gravity and the coefficient of thermal expansion, respectively; T ref is the average temperature. The force term establishes the coupling between the lattice Boltzmann equations for the flow field () and the temperature field (2). Note that through a multiscale Chapman-Enskoq expansion, the incompressible Navier-Stokes equations can be derived from the discussed incompressible LBGK model. After neglecting the viscous heat dissipation and compression work carried out by the pressure, the temperature field obeys a passive scalar equation. In sum, the approximated incompressible equations in this work are, cf. (Guo et al., 22), u =, (3) u t + (uu) = p + ν 2 u + F, (4) T t + (ut ) = D 2 T. (5) The kinematic viscosity ν and the thermal diffusivity D are related to the dimensionless collision times by ν = 6 (2τ ν ) c x and D = 4 (2τ D ) c x. Introducing the physical speed of sound as c s = c/ 3 these expressions yield the relations τ ν = ν + c2 s t/2 c 2 s t, τ D = D c2 s t/2 3 2 c2 s t (6) which can be used to evaluate the dimensionless collision times in () and (2) for given macroscopic gas properties ν, D and time step t Adaptive mesh refinement. For local dynamic mesh adaptation we have adopted the block-structured AMR method proposed by Berger and Colella (988). This method was originally designed for time-explicit finite volume schemes for hyperbolic conservation laws, however, its recursive execution procedure and natural consideration of time step refinement make it equally applicable to lattice Boltzmann schemes, which is is not surprising as a hyperbolic constant velocity advection equation is the theoretical underpinning of the transport step (4). In order to fit smoothly into our existing, fully parallelized finite volume AMR software system AMROC (Deiterding, 2), we have implemented the LBM cell-based. In the block-based AMR approach, finite volume cells are clustered with a special algorithm into non-overlapping rectangular grids. The grids have a suitable layer of halo cells for synchronization and applying inter-level and physical boundary conditions. Refinement levels are integrated recursively starting from the coarsest level. With index l denoting the AMR level, the spatial mesh width x l and the time step t l are refined by the same factor r l, where we assume r l 2 for l > and,
5 A dynamically adaptive lattice Boltzmann method for thermal convection problems (a) f C,n i,in (b) f F,n i,in cf. Fig. 2(b). F,n 3. Complete transport f i whole fine mesh. 5 := T (f F,n i ) on i := Collision f F,n+/2 F,n C( f i ) is applied only in the interior cells (yellow in Fig. 2(b)). (c) T C,n ( f i,out ) f F,n+ i,out Fig. 2. Visualization of distributions involved in data exchange at a coarse (C) - fine (F ) boundary. The thick black lines indicate a physical boundary. (a) Coarse distributions going into fine grid; (b) incoming interpolated fine distributions in halos (top), outgoing distributions in halos after two fine-level transport steps (bottom); (c) averaged distributions replacing coarse values before update is repeated in cells next to boundary. r =. In the adaptive thermal LBM, it is of foremost importance that the dimensionless collision times of the DDF LBM are adjusted on a level basis according to (6) as the time step is recursively refined. In addition, the interface region requires a specialized treatment to ensure consistent transport of coarse-grid distributions into refined cells and of fine-grid distributions into the coarse cells adjacent to the boundaries of refined regions. Since the D2Q4 stencil is just a simplified version of the D2Q9 method, we restrict our description of the interface algorithm to the latter. Distinguishing between the transport and collision operators, T and C, respectively (cf. (4) and (5)), our method proceeds in the following steps if a refinement factor of 2 is considered:. Complete update on coarse grid: f C,n+ i := CT (f C,n i ) 2. Use coarse grid distributions f C,n i,in that propagate into the fine grid, cf. Fig. 2(a), to construct initial fine grid halo values f F,n i,in, F,n+/2 4. Repeat 3. to obtain f i := T (f F,n+/2 i ) and f F,n+ F,n+/2 i := C( f i ). 5. Average outgoing distributions from fine grid halos (Fig. 2(c)), that is in the inner halo layer and C,n to obtain f i,out. f F,n i,out f F,n+/2 i,out (outer halo layer) 6. Revert transport for averaged outgoing C,n distributions, f i,out := T C,n ( f i,out ), and overwrite those in the previous coarse grid time step. 7. Synchronization of f C,n i, level. f C,n i,out on entire 8. Repeat complete update on coarse grid cells next to coarse-fine boundary only: f C,n+ i := CT (f C,n C,n i, f i,out ) In this description and in Fig. 2, the time steps on the coarse level C are indexed by the superscript n, index F denotes the fine level and the subscripts in and out indicate distributions which are convected in- and outwards of the fine grid along the coarse-fine boundary. The overall algorithm is computationally equivalent to the method by Chen et al. (26) but explicitly tailored to the Berger-Collela recursion that updates coarse grids in their entirety before fine grids are computed. The complete update of the entire respective coarse mesh and subsequent correction is the basis of the computational efficiency of the Berger-Collela method; however, this approach has so far hardly been applied to lattice Boltzmann methods. Previous adaptive LBM, cf. (Chen et al., 26), update the fine grid before the respective coarse level and provide no apparent avenue for implementing time-interpolated fine level interface conditions. While not being used above, the benefit of interpolating in time the non-equilibrium portion of coarse-grid distributions crossing
6 6 K. FELDHUSEN et al. the coarse-fine interface in Step 4 has been demonstrated by Dupuis and Chopard (23) and will be considered in our implementation in the future Wall boundary treatment. The correct implementation of boundary condition is very important for numerical stability. For the considered test cases we need different implementations of boundary conditions for the velocity and temperature partial distribution functions. No-slip or adiabatic boundary conditions are realized via a bounce-back approach for the unknown partial distribution functions as described in (Succi, 2). To prescribe fixed macroscopic values on the wall in form of Dirichlet boundary conditions we use a second order extrapolation scheme from Guo et al. (22). The outflow boundary conditions are implemented via a linear propagation as prescribed in (Mohamad, 2). We use a set of halo cells around the computational domain to manipulate the unknown partial probability distribution functions in the transport step Curved boundary treatment. We represent non-cartesian boundaries implicitly on the adaptive Cartesian grid by utilizing a scalar level set function ϕ that stores the distance to the boundary surface. The boundary surface is located exactly at ϕ = and the boundary outer normal in every mesh point can be evaluated as n = ϕ/ ϕ, (Deiterding, 2). We treat a fluid cell as an embedded ghost cell if its midpoint satisfies ϕ <. In order to implement non-cartesian boundary conditions with the LBM, we have chosen to pursue for now a st order accurate ghost fluid approach. In our technique, the density distributions in embedded ghost cells are adjusted to model the boundary conditions of a non-cartesian reflective wall moving with velocity vector w before applying the unaltered LBM. The last step involves interpolation and mirroring of p, T, u, across the boundary to p, T and ū and modification of the macroscopic velocity vector in the immersed boundary cells to u = 2w ū, cf. (Deiterding, 2). From the newly constructed macroscopic values the distributions in the embedded ghost cells are simply set to (p, u ) and g eq i (T ). f eq i 3. Results For the setup of physical configurations it is useful to recall the definitions of the dimensionless Rayleigh and Prandtl number which is Ra = gβ T H3 νd, Pr = ν D. (7) The characteristic velocity U for thermal convection flows is generally set to the buoyancy velocity U = gβ T H, where H denotes a problem-dependent geometric height. A cell (j, k) is flagged for refinement if any of the scaled gradient relations φ j+,k φ j,k >ɛ φ, φ j,k+ φ j,k >ɛ φ, φ j+,k+ φ j,k >ɛ φ (8) is satisfied for a particular macroscopic component φ j,k and a prescribed limit ɛ φ. If not stated otherwise, ɛ T is set to % of maximum temperature and ɛ u, ɛ v, ɛ w are set to 5% of characteristic velocity. 3.. Porous Plate. In order to validate the basic numerical method, we selected the problem of forced thermal convection between two porous plates also employed by Guo et al. (22). This problem is set up as a Couette flow between two porous plates of which the upper is in motion. A constant flow is injected normal to the lower plate and leaves the domain through the top plate with the same rate. The bottom plate is cooled, while the upper plate is heated. The analytic solutions for the horizontal velocity and the temperature profile in steady state are ( e u Re y/h (y) = U e Re T (y) = T C + T ), (9) ( e RePr y/h e RePr ), (2) where U is the velocity of the upper plate. The Reynolds number Re is based on the injection
7 A dynamically adaptive lattice Boltzmann method for thermal convection problems y/h Ana: Re=5 LBM: Re=5 Ana: Re= LBM: Re= Ana: Re=2 LBM: Re= Velocity Error ave [%] Iteration (*) 7 Re=5 Re= Re=2 u/u y/h Ana: Re=5 LBM: Re=5 Ana: Re= LBM: Re= Ana: Re=2 LBM: Re= (T-T C )/(T H -T C ) Temperature Error ave [%] Re=5 Re= Re= Iteration (*) Fig. 3. Comparison of velocity and temperature distribution predicted for different Re in comparison with analytic solution. Fig. 4. Averaged L2-norm error for computed macroscopic velocity and temperature over iteration steps for different Re. Table. Spatial averaged error: porous plate problem. Re E ave (u) [%] E ave (T ) [%] velocity V and is given by Re = V H ν. We study three different configurations with varying Reynolds number. The Prandtl number is fixed and set to Pr =.7, which corresponds to air and the Rayleigh number is set to Ra =. The velocity of the upper plate is also fixed and set to U =.. Finally, the dimensionless relaxation time τ ν on the coarsest level is prescribed as τ ν = /.25. The simulations are performed for the Reynolds numbers Re = 5, and 2 using a base grid of cells. Successive embedded static refinement with four additional levels with refinement factors r,...,4 = 4 is realized in the complete computational domain [, 64] [, 32]. In detail, we have the finest resolution r 4 near the top and bottom boundaries [, 64] ([, 4] & [28, 32]), then r 2 in [, 64] ([4, 8] & [24, 28]) and r 3 in [, 64] ([8, 2] & [2, 24]). The coarsest refinement level r is in the center region [, 64] [2, 2]. The entire velocity field is initialized at rest as (, ) T and the temperature field to the constant value T C. We compare the numerical predictions of the velocity and temperature distributions with the analytic solution. Figure 3 plots the normalized numerical results vs. the analytic solutions. From the point of validation, the macroscopic values for the horizontal velocity and scalar temperature are being calculated in each cell midpoint along each vertical line. The macroscopic values in the cells are averaged along the horizontal lines. The L2-norm error of the averaged macroscopic quantities Φ are calculated with (2) and displayed for the last iteration step in Table. Φ ave (x i ) Φ (x i ) 2 i E ave (Φ) = (2) Φ (x i ) 2 i The agreement is obviously excellent and below 2% for all three cases. It is noteworthy that the error for the velocity is smaller than the one for the temperature. When increasing the
8 8 K. FELDHUSEN et al. u = U v = T = T C v =, u T y =, y = T H u =, v = ω v =, u T y =, y = u x = v x = T x = Fig. 5. Setup for the flow past the heated rotating cylinder. discrete velocity directions for the temperature distribution functions from 6 to 9, this error should decrease. Figure 4 plots the averaged error for the computed macroscopic velocity and temperature over the computational iteration steps. The convergence to a fixed value is obvious Fluid flow past a heated rotating cylinder. In order to test the dynamic adaptation capabilities and boundary conditions for embedded complex geometries, we study the setup of a two-dimensional fluid flow past a heated isothermal rotating cylinder. The origin of the coordinate system is located in the center of the cylinder. As shown in Figure 5, the left boundary is an inlet with constant temperature T C, zero vertical velocity and constant inflow velocity U. On the right hand side of the domain, an outlet is modeled by imposing zero horizontal gradient boundary conditions for velocity and temperature. Slip adiabatic wall boundary conditions are applied at the upper and lower boundary. The cylinder boundary is modeled as a no-slip wall, which is isothermally heated to the constant temperature T H and has the constant prescribed angular velocity Ω. In terms of the cylinder radius R = 5, the computational domain has the extensions [ 6R, 6R] [ 8R, 8R], which is sufficiently large to eliminate boundary influences on the solution (Yan and Zu, 28). A base grid of cells is used and three additional levels refined by the factors r = 2 for level and r 2,3 = 4 for the other levels are applied. The dynamic refinement is based on scaled gradients of the velocity components as well as the temperature. The entire velocity field is initialized as (U, ) T and the temperature field to the constant value T C. The Reynolds number is given by Re = 2U R/ν and is set to Re = 2, where U =. is used. The peripheral velocity V of the rotating cylinder is given by V = ΩR. With the parameter k = V/U =.5 prescribed, we can determine V and the angular velocity Ω. To allow the direct comparison to the experimental results by Coutanceau and Menard (985) the Prandtl number is set to Pr =.5 and all variables are normalized with the reference length R and U as velocity. Further, T T C T H T C defines the reference temperature and the time normalization factor follows as R/U. Figure 6 shows the dynamic adaption during the computation at four different time points by displaying streamlines and the domains of different mesh refinement levels. The onset of vortex shedding can be inferred. The finest refinement level (red) is located directly around the cylinder. Namely, where the boundary layers are located and detach from the cylinders surface. The unrefined regions, colored in blue, are in the outer regions of the domain. The refined levels move downstream with the shedding vortices and the cylinder wake increases over time. Figure 7 compares the temporal evolution of the velocity components along representative points on the x-axis obtained in the simulation and with data from the experiment, while Fig. 8 displays the time evolution of the scalar temperature versus numerical results reported by Lai and Yan (2). The latter adopted a finite volume method with non-orthogonal grids. Again, our simulation results are in good agreement with some differences in the u-velocity component at t = 8 when the vortex is shed (see Fig. 6). A possible explanation is our rather simple temperature operator with only four discrete unity directions and with the used boundary conditions for the curved boundary explained in Section 2.4. However, by using the bounce back scheme for curved moving boundaries from Bouzidi et al. (2) and Li et al. (23) with a global uniform mesh the differences are considerably reduced, cf. Fig. 9. Therefore, the next step is to implement the curved boundary treatment in the
9 A dynamically adaptive lattice Boltzmann method for thermal convection problems 9.5 u * t = t * =3 (LBM) t * =3 (Exp.) t * =8 (LBM) t * =8 (Exp.) x *.5 v * t = t * =3 (LBM) t * =3 (Exp.) t * =8 (LBM) t * =8 (Exp.) x * Fig. 7. Time evolution of the velocity components along the x-axis for Re = 2 and k =.5. t = 8.8 t * = 3 (LBM) t * = 3 (FVM) t * = 6 (LBM) t * = 6 (FVM) t * = 2 (LBM) t * = 2 (FVM).6 T *.4 t =2 Fig. 6. Evolution of the velocity field and the adaptive mesh refinement regions for Re = 2 and k = x * Fig. 8. Time evolution of the temperature along the x- axis for Re = 2, Pr =.5 and k =.5.
10 K. FELDHUSEN et al. T * t * =2 (AMR - OldBC) t * =2 (UNI - NewBC) t * = 2 (FVM) t * =6 (AMR - OldBC) t * =6 (UNI - NewBC) t * = 6 (FVM) x * Fig. 9. Comparison of simulation results with different used curved boundary conditions: Time evolution of the Temperature along the x-axis for Re = 2, Pr =.5 and k =.5. T = T H u = v = T/ y =, u =, v = y x H T/ y =, u =, v = g H T = T C u = v = Fig.. Configuration of the two dimensional cavity. AMR method Natural convection in a square 2D-cavity. In order to benchmark the overall method we employ a two-dimensional square cavity with differentially heated walls. At the vertical walls isothermal temperatures T H and T C are prescribed and adiabatic boundary conditions are applied at top and bottom. Further, at all four walls we prescribe no-slip boundary conditions for the velocity field. Figure depicts this setup. The flow is characterized by the Prandtl number Pr =.7 (air) and the Rayleigh numbers Ra = j with j = 3,..., 8 with accordingly increasing velocities U. The reference temperature is given by T ref = (T H + T C )/2. The simulations were terminated after reaching steady state. Two additional levels of refinement with r,2 = 2 are used and the base mesh has (H x ) 2 cells, whereby x = and H is given in the left column of Table 2. For simulations with Ra = 3,, 6 we use the defined refinement thresholds for horizontal and vertical velocity ɛ u, ɛ v with 2.5% of the characteristic velocity and % of the maximum temperature. The thresholds for Ra = 7 and 8 remain as previously stated. We compare our adaptive simulation results to published reference data by De Vahl Davis (983), who solved the NS equations on a uniform square mesh with a second order finite difference method, and by Guo et al. (22), who used the incompressible thermal LBGK approach presented above with a uniform mesh. Further results by Kuznik et al. (27), who used a D2Q9 DDF LBM approach with non-uniform mesh resolution, are listed in Table 2. Table 2 contains the obtained maximal horizontal velocity u max along the vertical center line at x = H/2 and the location y max of its occurrence and similarly for the horizontal center line at y = H/2, the maximal vertical velocity v max and its location x max. Furthermore, the average Nusselt number H T Nu ave = T x dy (22) x= is compared. Velocity values in Table 2 are normalized by the reference diffusion velocity D/H. As expected, u max, v max and Nu ave increase with increasing Rayleigh number Ra. Comparing the Nu numbers predicted by our adaptive method to the literature data, an agreement within 2 % is found for all Ra numbers. Figure shows the vertical velocity component in the horizontal mid-plane for all discussed Rayleigh numbers. The velocity profiles plotted in Fig. reveal the development of a boundary layer close to the heated/cooled walls with velocity maxima/minima whose values increase/decrease with increasing/decreasing Ra. This increase of the magnitude of the vertical velocity with increasing Ra is also reflected in Table 2. To give an impression of the flow solution, contours of the temperature fields and
11 A dynamically adaptive lattice Boltzmann method for thermal convection problems Table 2. Comparison of the simulation results: natural convection in the square cavity. Ref. u max y max v max x max Nu ave Ra = 3 a U =. b H = c d Ra = 4 a U =.2 b H = 5 c d Ra = 5 a U =.5 b H = 2 c d Ra = 6 a U =.5 b H = 2 c d Ra = 7 a U =.5 d H = 256 Ra = 8 a U =.5 d H = 256 v a = Present (LBM-AMROC), b = (De Vahl Davis, 983) (FDM - uniform), c = (Guo et al., 22) (LBM - uniform), d = (Kuznik et al., 27) (LBM - nonuniform) x/h Ra= 3 Ra= 4 Ra= 5 Ra= 6 Ra= 7 Ra= 8 Fig.. Vertical velocity in the horizontal mid-plane of the 2D cavity for different Rayleigh numbers. streamlines are presented in Fig. 2 for three considered Ra numbers. For all three Ra numbers the streamlines reflect that fluid rises at the heated wall and descends at the cooled wall. This generates a circulation around the center where the velocity is zero. For the lower Ra numbers the computed flow field are in good agreement with results reported in previous studies (De Vahl Davis, 983; Guo et al., 22; Azwadi Che Sidik and Irwan, 2; Kuznik et al., 27; Abdelhadi et al., 26). In the graph with the contours predicted for Ra = 7 the mesh refinement levels realized in the domain are additionally highlighted by colors. From the predominantly vertical isotherms obtained for the low Ra number case it can be concluded that the heat conduction dominates the heat transport between the heated walls. For larger Ra the isotherms are aligned more horizontally in the cavity s center due to the thinner boundary layers. The denser isotherms near the hot and cold wall further reflect the lower thermal boundary layer thickness for higher Rayleigh number. It is in this region where on-the-fly mesh resolution is particularly beneficial Natural convection in a cubic cavity. To benchmark the three-dimensional implementation of the method, we employ a 3D cubic cavity with differentially heated walls. As before, at the vertical walls the constant temperatures T H and T C are prescribed. At the bottom, top and front, back walls adiabatic boundary conditions are used for the temperature, while no-slip boundary conditions at all six walls are realized for the velocity fields. In summary, Fig. 3 represents this numerical setup. Again, the Prandtl number is Pr =.7 (air) and in the 3D simulations the Rayleigh numbers Ra = j is varied from j = 3,..., 5. Here, we focus on the flow for Ra 5, since for higher Ra the flow is expected to become unsteady and eventually turbulent. To benchmark our method for a turbulent flow is however beyond the scope of this paper. As discussed above, the buoyancy (reference) velocity U rises with increasing Ra and the reference temperature is given by T ref = (T H +T C )/2. Two additional levels of refinement with r = 2, r 2 = 4 are used and the base
12 2 K. FELDHUSEN et al. Ra = 3 Ra = 6 Ra = 7 Fig. 2. LBM results of natural convective flow in the square cavity for three Ra numbers. Left: contours of isotherms. Right: streamlines. H T H T C g y H H z x Fig. 3. Configuration of the three dimensional cavity. mesh has (H x ) 3 cells, whereby x = and H is given in the left column of Table 3. The adaptive mesh refinement obeys the scaled gradient criteria given above in (8). The used thresholds for all three velocity components are %, 2% and 5% of the reference velocity U for Ra = 3, 4 and 5, respectively. As before, % of T H is used as the temperature refinement threshold. The computed results are compared to published literature results after reaching steady state. Azwadi Che Sidik and Syahrullail (29) use a D3Q9 DDF LBM approach with D3Q6 operator for the temperature field and a uniform cubic mesh to get excellent numerical stability and accuracy. Peng et al. (23) use a three-dimensional incompressible LBM with DDF approach and two D3Q9 operators for the two fields and a non-uniform mesh resolution. Finally, Fusegi et al. (99) use a high-resolution, finite difference NS solver with a uniform mesh resolution result and obtain results which agree reasonably well with experimental measurements. Figure 4 visualizes the temperature isosurfaces in the cubic enclosure and the different mesh refinement levels in the symmetry plane for Ra = 4, 5. Near the heated walls, the isosurfaces are predominantly vertical. Notice, that the isosurfaces in the center of the cavity become more horizontally with increasing Ra. The reason is that the thermal boundary layer is becoming thinner. This observation is similar to that in the previous chapter. Note that the shaping of the mesh refinement levels for Ra = 4 is much more pronounced than for Ra = 5. As in the previous chapter, we compare the results in the symmetry plane z = H/2 in terms of maximal horizontal velocity u max along the vertical center line at x = H/2 and at the corresponding location y max of its occurrence and similarly for the horizontal center line at y = H/2, the maximal vertical velocity v max and its location x max. Furthermore, we use the average Nusselt number (22) for comparison. Our results are listed in Table 3. The velocity values in Table 3 are normalized with the reference velocity U. The Nusselt number increases with increasing Ra number, which means that the convective part of the heat transfer predominates the conduction. Comparing the Nu numbers predicted with our
13 A dynamically adaptive lattice Boltzmann method for thermal convection problems 3 Table 3. Comparison of the simulation results: natural convection in the cubic cavity. Ref. u max y max v max x max Nu ave Ra = 3 a U =. e H = 8 f g Ra = 4 a U =.2 e H = 8 f g Ra = 4 Ra = 5 a U =. e H = 9 f g a = Present (LBM-AMROC), e = Azwadi et al. (Azwadi Che Sidik and Syahrullail, 29) (LBM - uniform), f = Peng et al. (Peng et al., 23) (LBM - nonuniform), g = Fusegi et al. (Fusegi et al., 99) (NS - uniform) method to the literature, an agreement within 2% is found for all three Ra numbers, although the comparison of the horizontal velocity component shows larger differences. The reason for this might be a lack of dynamic mesh refinement near the upper and bottom walls. The mesh refinement is more pronounced near the heated and cooled walls, where the thinner thermal boundary layers are located. 4. Conclusions A novel two and three dimensional incompressible dynamically adaptive thermal lattice Boltzmann method on block-based hierarchical finite volume meshes with embedded complex geometric structures has been developed and validated. The agreement for a two-dimensional porous plate problem on a Cartesian grid is nearly perfect. Successful validation against analytic solutions of the Navier-Stokes equations, e.g., for a heated rotating cylinder for Pr =.5 has been achieved. While for this particular example the deviations in velocity and temperature were found to increase over time, a possible improvement could be the implementation of a bounce-back boundary condition for curved boundaries. For the benchmark of a two-dimensional heated cavity with Rayleigh numbers from Ra = 3 to 8, the predictions are in good agreement Ra = 5 Fig. 4. Simulation results of natural convective flow in the cubic cavity. Left: isosurfaces of temperature. Color indicates temperature from hot (red) to cold (blue). Right: mesh refinement levels. with published results. Our results in form of the computed Nusselt number reach an agreement within 2%. For higher Rayleigh numbers, the deviations in the considered quantities are greater in regions without refinement. The comparison for a three-dimensional heated cubic cavity with Rayleigh numbers from Ra = 3 to 5 against literature results delivers a good agreement as well. In terms of the Nusselt number, the agreement with literature results is again under 2%. A comprehensive analysis of CPU-time and memory savings by employing our unique block-based adaptive LBM will be conducted in the future. We will also take a closer look at how the results are influenced by the refinement criteria. Finally, extension and validation of the 3D approach to turbulent flows at higher Ra or Re numbers is planned. 5. Biographies The author biographies will be provided upon acceptance.
14 4 K. FELDHUSEN et al. References Abdelhadi, B., Hamza, G., Razik, B. and Raouache, E. (26). Natural convection and turbulent instability in cavity, WSEAS Transactions on Heat and Mass Transfer (2): Alexander, F. J., Chen, S. and Sterling, J. D. (993). Lattice boltzmann thermohydrodynamics, Physical Review E 47: Azwadi Che Sidik, N. and Irwan, M. (2). Simplified mesoscale lattice boltzmann numerical model for prediction of natural convection in a square enclosure filled with homogeneous porous media, WSEAS Transactions on Fluid Mechanics 3(5): Azwadi Che Sidik, N. and Syahrullail, S. (29). A three-dimension double-population thermal lattice bgk model for simulation of natural convection heat transfer in a cubic cavity, WSEAS Transactions on Mathematics 8(9): Berger, M. and Colella, P. (988). Local adaptive mesh refinement for shock hydrodynamics, Journal of Computational Physics 82: Bhatnagar, P., Gross, E. and Krook, M. (954). A model for collisional processes in gases I: small amplitude processes in charged and in neutral one-component systems, Physical Review. 94: Bouzidi, M., Firdaouss, M. and Lallemand, P. (2). Momentum transfer of a boltzmann-lattice fluid with boundaries, Physics of Fluids 3: Chen, H., Filippova, O., Hoch, J., Molvig, K., Shock, R., Teixeira, C. and Zhang, R. (26). Grid refinement in lattice Boltzmann methods based on volumetric formulation, Physica A 362: Chen, H. and Teixeira, C. (2). H-theorem and origins of instability in thermal lattice boltzmann models, Computer Physics Communications 29(): 2 3. Chen, S. and Doolen, G. (998). Lattice Boltzmann method for fluid flows, Annual Review of Fluid Mechanics 3: Coutanceau, M. and Menard, C. (985). Influence of rotation on the near-wake development behind an impulsively started circular cylinder, Journal of Fluid Mechanics 58: De Vahl Davis, G. (983). Natural convection of air in a square cavity a benchmark numerical solution, International Journal for Numerical Methods in Fluids 3: Deiterding, R. (2). Block-structured adaptive mesh refinement - theory, implementation and application, ESAIM Proceedings 34: Deller, P. (22). Lattice kinetic schemes for magnetohydrodynamics, Journal of Computational Physics 79(): Dupuis, A. and Chopard, B. (23). Theory and applications of an alternative lattice Boltzmann grid refinement algorithm, Physica E 67: Fusegi, T., Hyun, J., Kuwahara, K. and Farouk, B. (99). A numerical study of three-dimensional natural convection in a differentially heated cubical enclosure, International Journal of Heat and Mass Transfer 34(6): Guo, Z., Shi, B. and Zheng, C. (22). A coupled lattice BGK model for the Boussinesq equations, International Journal for Numerical Methods in Fluids 39: Hähnel, D. (24). Molekulare Gasdynamik, Springer. He, N.-Z., Wang, N.-C., Shi, B.-C. and Guo, Z.-L. (24). A unified incompressible lattice BGK model and its application to three-dimensional lid-driven cavity flow, Chinese Physics 3. He, X., Chen, S. and Doolen, G. (998). A novel thermal model for the lattice Boltzmann method in incompressible limit, Journal of Computational Physics 46: He, X. and Luo, L.-S. (997). Lattice Boltzmann model for the incompressible Navier-Stokes equation, Journal of Statistical Physics 88: Jonas, L., Chopard, B., Succi, S. and Toschi, F. (26). Numerical analysis of the average flow field in a turbulent lattice Boltzmann simulation, Physica A 32(): 6. Kuznik, F., Vareilles, J., Rusaouen, G. and Kraiss, G. (27). A double-population lattice Boltzmann method with non-uniform mesh for the simulation of natural convection in a square cavity, International Journal of Heat and Fluid Flow 28: Lai, H. and Yan, Y. (2). The effect of choosing dependent variables and cell-face velocities on convergence of the SIMPLE algorithm using non-orthohonal grids, International Journal of Numerical Methods for Heat & Fluid Flow : Lee, T. and Lin, C. (25). A stable discretization of the lattice Boltzmann equation for simulation of incompressible two-phase flows at high density ratio, Journal of Computational Physics 26(): 6 47.
15 A dynamically adaptive lattice Boltzmann method for thermal convection problems 5 Li, L., Mei, R. and Klausner, J. (23). Boundary conditions for thermal lattice Boltzmann equation method, Journal of Computational Physics 237: McNamara, G. and Alder, B. (993). Analysis of lattice Boltzmann treatment of hydrodynamics, Physica A 94(): Mohamad, A. (2). Lattice Boltzmann Method - Fundamentals and Engineering Applications with Computer Codes, Springer London. Mohamad, A. and Kuzmin, A. (2). A critical evaluation of force term in lattice Boltzmann method, natural convection problem, International Journal of Heat and Mass Transfer 53: Peng, Y., Shu, C. and Che, Y. (23). A 3d incompressible thermal lattice Boltzman model and its application to simulation natural convection in a cubic cavity, Journal of Computational Physics 93: Qian, Y. (993). Simulating thermohydrodynamics with lattice BGK models, Journal of Scientific Computing 8(3): Qian, Y., D Humires, D. and Lallemand, P. (992). Lattice BGK models for Navier-Stokes equation, Europhysics Letters 7: 479. Succi, S. (2). The Lattice Boltzmann Equation: For Fluid Dynamics and Beyond, Numerical Mathematics and Scientific Computation, Clarendon Press. Yan, Y. and Zu, Y. (28). Numerical simulation of heat transfer and fluid flow past a rotating isothermal cylinder - A LBM approach, International Journal of Heat and Mass Transfer 5: Yu, Z. and Fan, L.-S. (29). An interaction potential based lattice Boltzmann method with adaptive mesh refinement (AMR) for two-phase flow simulation, Journal of Computational Physics 7: Received:
Numerical 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
Lattice 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
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
Lattice 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
SIMULATION OF MIXED CONVECTIVE HEAT TRANSFER USING LATTICE BOLTZMANN METHOD
International Journal of Automotive and Mechanical Engineering (IJAME) ISSN: 2229-8649 (Print); ISSN: 2180-1606 (Online); Volume 2, pp. 130-143, July-December 2010 Universiti Malaysia Pahang DOI: http://dx.doi.org/10.15282/ijame.2.2010.3.0011
arxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
Simulation 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
An 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
A Compact and Efficient Lattice Boltzmann Scheme to Simulate Complex Thermal Fluid Flows
A Compact and Efficient Lattice Boltzmann Scheme to Simulate Complex Thermal Fluid Flows Tao Zhang 1[0000-0001-7216-0423] and Shuyu Sun 1[0000-0002-3078-864X] 1 Computational Transport Phenomena Laboratory
LATTICE BOLTZMANN SIMULATION OF FLUID FLOW IN A LID DRIVEN CAVITY
LATTICE BOLTZMANN SIMULATION OF FLUID FLOW IN A LID DRIVEN CAVITY M. Y. Gokhale, Ignatius Fernandes Maharashtra Institute of Technology, Pune 4 38, India University of Pune, India Email : mukundyg@yahoo.co.in,
68 Guo Wei-Bin et al Vol. 12 presented, and are thoroughly compared with other numerical data with respect to the Strouhal number, lift and drag coeff
Vol 12 No 1, January 2003 cfl 2003 Chin. Phys. Soc. 1009-1963/2003/12(01)/0067-08 Chinese Physics and IOP Publishing Ltd Lattice-BGK simulation of a two-dimensional channel flow around a square cylinder
Thermal 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
Simplified Mesoscale Lattice Boltzmann Numerical Model for Prediction of Natural Convection in a Square Enclosure filled with Homogeneous Porous Media
Simplified Mesoscale Lattice Boltzmann Numerical Model for Prediction of Natural Convection in a Square Enclosure filled with Homogeneous Porous Media C.S.NOR AZWADI M.A.M. IRWAN Faculty of Mechanical
Due Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
NATURAL CONVECTION FLOW IN A SQUARE CAVITY WITH INTERNAL HEAT GENERATION AND A FLUSH MOUNTED HEATER ON A SIDE WALL
Journal of Naval Architecture and Marine Engineering December, 2010 DOI: 10.3329/jname.v7i2.3292 http://www.banglajol.info NATURAL CONVECTION FLOW IN A SQUARE CAVITY WITH INTERNAL HEAT GENERATION AND A
On the influence of tube row number for mixed convection around micro tubes
Thessaloniki, Greece, 22-24 August 211 On the influence of tube row number for mixed convection around micro tubes Chuanshan DAI, Qiuxiang WANG, Biao LI * Corresponding author: Tel.: +86-22274183; Fax:
SELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE
Computational Thermal Sciences, 3 (1): 63 72 (2011) SELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE M. Zhao, 1, M. Yang, 1 M. Lu, 1 & Y. W. Zhang 2 1
Drag Force Simulations of Particle Agglomerates with the Lattice-Boltzmann Method
Drag Force Simulations of Particle Agglomerates with the Lattice-Boltzmann Method Christian Feichtinger, Nils Thuerey, Ulrich Ruede Christian Binder, Hans-Joachim Schmid, Wolfgang Peukert Friedrich-Alexander-Universität
APPRAISAL 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
Analysis of Turbulent Free Convection in a Rectangular Rayleigh-Bénard Cell
Proceedings of the 8 th International Symposium on Experimental and Computational Aerothermodynamics of Internal Flows Lyon, July 2007 Paper reference : ISAIF8-00130 Analysis of Turbulent Free Convection
FREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER
Free Convective Heat Transfer From an Object at Low Rayleigh Number FREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER Md. Golam Kader and Khandkar Aftab Hossain * Department of Mechanical
FINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT OF A HEATED SQUARE HOLLOW CYLINDER IN A LID-DRIVEN RECTANGULAR ENCLOSURE
Proceedings of the International Conference on Mechanical Engineering 2011 (ICME2011) 18-20 December 2011, Dhaka, Bangladesh ICME11-TH-014 FINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT
Predictor-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.
FINITE-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
Numerical 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
Using 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,
NUMERICAL STUDIES OF TRANSITION FROM STEADY TO UNSTEADY COUPLED THERMAL BOUNDARY LAYERS
International Journal of Computational Methods Vol. 11, Suppl. 1 (214) 13442 (15 pages) c World Scientific Publishing Company DOI: 1.1142/S2198762134427 NUMERICAL STUDIES OF TRANSITION FROM STEADY TO UNSTEADY
Lattice 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.:
A 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
Alternative and Explicit Derivation of the Lattice Boltzmann Equation for the Unsteady Incompressible Navier-Stokes Equation
ISSN (e): 2250 3005 Volume, 06 Issue, 12 December 2016 International Journal of Computational Engineering Research (IJCER) Alternative and Explicit Derivation of the Lattice Boltzmann Equation for the
Simulation 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
Analysis and boundary condition of the lattice Boltzmann BGK model with two velocity components
Analysis and boundary condition of the lattice Boltzmann BGK model with two velocity components Xiaoyi He and Qisu Zou arxiv:comp-gas/9507002v1 22 Jul 1995 Abstract In this paper, we study the two dimensional
Analysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle
Nonlinear Analysis: Modelling and Control, 2011, Vol. 16, No. 1, 89 99 89 Analysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle S. Parvin,
Study 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
Three-dimensional simulation of slip-flow and heat transfer in a microchannel using the lattice Boltzmann method
75 Three-dimensional simulation of slip-flow and heat transfer in a microchannel using the lattice Boltzmann method A. C. M. Sousa,2, M. Hadavand & A. Nabovati 3 Department of Mechanical Engineering, University
On pressure and velocity boundary conditions for the lattice Boltzmann BGK model
On pressure and velocity boundary conditions for the lattice Boltzmann BGK model Qisu Zou Theoretical Division, Los Alamos National Lab, Los Alamos, New Mexico 7545 and Department of Mathematics, Kansas
Coupling 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
Natural Convection in Parabolic Enclosure Heated from Below
www.ccsenet.org/mas Modern Applied Science Vol. 5, No. 3; June 011 Natural Convection in Parabolic Enclosure Heated from Below Dr. Ahmed W. Mustafa (Corresponding auther) University of Tikrit, College
Visualization of Natural Convection in Enclosure. Filled with Porous Medium by Sinusoidally. Temperature on the One Side
Applied Mathematical Sciences, Vol., 2012, no. 97, 801-812 Visualization of Natural Convection in Enclosure Filled with Porous Medium by Sinusoidally Temperature on the One Side Paweena Khansila Department
Validation 3. Laminar Flow Around a Circular Cylinder
Validation 3. Laminar Flow Around a Circular Cylinder 3.1 Introduction Steady and unsteady laminar flow behind a circular cylinder, representing flow around bluff bodies, has been subjected to numerous
Study of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
SIMULATION OF NATURAL CONVECTION HEAT TRANSFER USING NANOFLUID IN A CONCENTRIC ANNULUS
THERMAL SCIENCE: Year 17, Vol. 1, No., pp. 175-186 175 SIMULATION OF NATURAL CONVECTION HEAT TRANSFER USING NANOFLUID IN A CONCENTRIC ANNULUS by Keivan FALLAH a*, Atena GHADERI b, Nima SEDAGHATIZADEH c,
Effect of Buoyancy Force on the Flow Field in a Square Cavity with Heated from Below
International Journal of Discrete Mathematics 017; (): 43-47 http://www.sciencepublishinggroup.com/j/dmath doi: 10.11648/j.dmath.01700.13 Effect of Buoyancy Force on the Flow Field in a Square Cavity with
The Use of Lattice Boltzmann Numerical Scheme for Contaminant Removal from a Heated Cavity in Horizontal Channel
www.cfdl.issres.net Vol. 6 (3) September 2014 The Use of Lattice Boltzmann Numerical Scheme for Contaminant Removal from a Heated Cavity in Horizontal Channel Nor Azwadi Che Sidik C and Leila Jahanshaloo
NATURAL CONVECTION OF AIR IN TILTED SQUARE CAVITIES WITH DIFFERENTIALLY HEATED OPPOSITE WALLS
Proceedings of the International onference on Mechanical Engineering 0 (IME0 8-0 December 0, Dhaka, Bangladesh IME- NATURAL ONVETION OF AIR IN TILTED SQUARE AVITIES WIT DIFFERENTIALLY EATED OPPOSITE WALLS
VORTEX SHEDDING PATTERNS IN FLOW PAST INLINE OSCILLATING ELLIPTICAL CYLINDERS
THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1395-1399 1395 VORTEX SHEDDING PATTERNS IN FLOW PAST INLINE OSCILLATING ELLIPTICAL CYLINDERS by Li-Zhong HUANG a* and De-Ming NIE b a State Key Laboratory
SEMICLASSICAL 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
Lattice 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.
NON-DARCY POROUS MEDIA FLOW IN NO-SLIP AND SLIP REGIMES
THERMAL SCIENCE, Year 2012, Vol. 16, No. 1, pp. 167-176 167 NON-DARCY POROUS MEDIA FLOW IN NO-SLIP AND SLIP REGIMES by Antonio F. MIGUEL Geophysics Centre of Evora & Department of Physics, University of
Simulation 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
Lattice 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
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
The Study of Natural Convection Heat Transfer in a Partially Porous Cavity Based on LBM
Send Orders for Reprints to reprints@benthamscience.ae 88 The Open Fuels & Energy Science Journal, 014, 7, 88-93 Open Access The Study of Natural Convection Heat Transfer in a Partially Porous Cavity Based
Flow patterns and heat transfer in square cavities with perfectly conducting horizontal walls: the case of high Rayleigh numbers ( )
Advances in Fluid Mechanics VII 391 Flow patterns and heat transfer in square cavities with perfectly conducting horizontal walls: the case of high Rayleigh numbers (10 6 10 9 ) R. L. Frederick & S. Courtin
Lattice 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
Entropy 2011, 13, ; doi: /e OPEN ACCESS. Entropy Generation at Natural Convection in an Inclined Rectangular Cavity
Entropy 011, 13, 100-1033; doi:10.3390/e1305100 OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Article Entropy Generation at Natural Convection in an Inclined Rectangular Cavity Mounir
Research 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,
Coupled double-distribution-function lattice Boltzmann method for the compressible Navier-Stokes equations
PHYSICAL REVIEW E 76, 5675 7 Coupled double-distribution-function lattice Boltzmann method for the compressible Navier-Stokes equations Q. Li, Y. L. He,* Y. Wang, and W. Q. Tao State Key Laboratory of
UNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
MIXED CONVECTION IN A SQUARE CAVITY WITH A HEAT-CONDUCTING HORIZONTAL SQUARE CYLINDER
Suranaree J. Sci. Technol. Vol. 17 No. 2; April - June 2010 139 MIXED CONVECTION IN A SQUARE CAVITY WITH A HEAT-CONDUCTING HORIZONTAL SQUARE CYLINDER Md. Mustafizur Rahman 1 *, M. A. Alim 1 and Sumon Saha
EFFECT OF THE INLET OPENING ON MIXED CONVECTION INSIDE A 3-D VENTILATED CAVITY
THERMAL SCIENCE: Year 2018, Vol. 22, No. 6A, pp. 2413-2424 2413 EFFECT OF THE INLET OPENING ON MIXED CONVECTION INSIDE A 3-D VENTILATED CAVITY by Hicham DOGHMI *, Btissam ABOURIDA, Lahoucin BELARCHE, Mohamed
Effect of an adiabatic fin on natural convection heat transfer in a triangular enclosure
American Journal of Applied Mathematics 2013; 1(4): 78-83 Published online November 10, 2013 (http://www.sciencepublishinggroup.com/j/ajam) doi: 10.11648/j.ajam.20130104.16 Effect of an adiabatic fin on
NATURAL CONVECTION WITHIN TRAPEZOIDAL ENCLOSURE WITH TWO BAFFLES: EFFECT OF VARIOUS ANGLES OF INCLINATION
NATURAL CONVECTION WITHIN TRAPEZOIDAL ENCLOSURE WITH TWO BAFFLES: EFFECT OF VARIOUS ANGLES OF INCLINATION Éliton Fontana, eliton_fontana@hotmail.com Universidade Federal de Santa Catarina - UFSC Adriano
Natural Convection and Entropy Generation in a Porous Enclosure with Sinusoidal Temperature Variation on the Side Walls
Avestia Publishing Journal of Fluid Flow, Heat and Mass Transfer Volume 1, Year 14 Journal ISSN: 368-6111 DOI: 1.11159/jffhmt.14.4 Natural Convection and Entropy Generation in a Porous Enclosure with Sinusoidal
A Finite Element Analysis on MHD Free Convection Flow in Open Square Cavity Containing Heated Circular Cylinder
American Journal of Computational Mathematics, 2015, 5, 41-54 Published Online March 2015 in SciRes. http://www.scirp.org/journal/ajcm http://dx.doi.org/10.4236/ajcm.2015.51003 A Finite Element Analysis
LATTICE 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
Available online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (24 ) 55 556 th International Conference on Mechanical Engineering, ICME 23 Analysis of heat transfer and flow due to natural
NATURAL CONVECTION HEAT TRANSFER IN PARTIALLY OPEN ENCLOSURES CONTAINING AN INTERNAL LOCAL HEAT SOURCE
Brazilian Journal of Chemical Engineering ISSN 0104-6632 Printed in Brazil www.abeq.org.br/bjche Vol. 24, No. 03, pp. 375-388, July - September, 2007 NATURAL CONVECTION HEAT TRANSFER IN PARTIALLY OPEN
Calculations on a heated cylinder case
Calculations on a heated cylinder case J. C. Uribe and D. Laurence 1 Introduction In order to evaluate the wall functions in version 1.3 of Code Saturne, a heated cylinder case has been chosen. The case
DIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING
DIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING Yukinori Kametani Department of mechanical engineering Keio
Natural Convection in Vertical Channels with Porous Media and Adiabatic Extensions
Natural Convection in Vertical Channels with Porous Media and Adiabatic Extensions Assunta Andreozzi 1,a, Bernardo Buonomo 2,b, Oronzio Manca 2,c and Sergio Nardini 2,d 1 DETEC, Università degli Studi
MODELLING OF THE BOUNDARY CONDITION FOR MICRO CHANNELS WITH USING LATTICE BOLTZMANN METHOD (LBM)
Reports Awarded with "Best Paper" Crystal Prize 17 FRI-1.417-1-MEMBT-05 MODELLING OF THE BOUNDARY CONDITION FOR MICRO CHANNELS WITH USING LATTICE BOLTZMANN METHOD (LBM) Rsc. Asst. İlkay Çolpan, BSc Department
Contents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
Numerical Study of Free Convection Heat Transfer in a Square Cavity with a Fin Attached to Its Cold Wall
Heat Transfer Research, 2011, Vol. 42, No. 3 Numerical Study of Free Convection Heat Transfer in a Square Cavity with a Fin Attached to Its Cold Wall SAEID JANI, 1* MEYSAM AMINI, 2 and MOSTAFA MAHMOODI
Simulation of Lid-driven Cavity Flow by Parallel Implementation of Lattice Boltzmann Method on GPUs
Simulation of Lid-driven Cavity Flow by Parallel Implementation of Lattice Boltzmann Method on GPUs S. Berat Çelik 1, Cüneyt Sert 2, Barbaros ÇETN 3 1,2 METU, Mechanical Engineering, Ankara, TURKEY 3 METU-NCC,
arxiv:cond-mat/ v1 [cond-mat.stat-mech] 11 Dec 2002
arxiv:cond-mat/0212257v1 [cond-mat.stat-mech] 11 Dec 2002 International Journal of Modern Physics B c World Scientific Publishing Company VISCOUS FINGERING IN MISCIBLE, IMMISCIBLE AND REACTIVE FLUIDS PATRICK
Principles of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
Lattice 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,
Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection
Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection Ahmed Waheed Mustafa 1 Mays Munir Ismael 2 AL-Nahrain University College of Engineering Mechanical Engineering Department ahmedwah@eng.nahrainuniv.edu.iq
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
Meysam ATASHAFROOZ, Seyyed Abdolreza GANDJALIKHAN NASSAB, and Amir Babak ANSARI
THERMAL SCIENCE: Year 014, Vol. 18, No., pp. 479-49 479 NUMERICAL INVESTIGATION OF ENTROPY GENERATION IN LAMINAR FORCED CONVECTION FLOW OVER INCLINED BACKWARD AND FORWARD FACING STEPS IN A DUCT UNDER BLEEDING
Numerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus
Nonlinear Analysis: Modelling and Control, 2007, Vol. 12, No. 4, 541 552 Numerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus M. Sankar
Turbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
Combined Natural Convection and Thermal Radiation in an Inclined Cubical Cavity with a Rectangular Pins Attached to Its Active Wall
Periodicals of Engineering and Natural Sciences ISSN 2303-4521 Vol.5, No.3, November 2017, pp. 347~354 Available online at:http://pen.ius.edu.ba Combined Natural Convection and Thermal Radiation in an
Research Article Study of the Transient Natural Convection of a Newtonian Fluid inside an Enclosure Delimited by Portions of Cylinders
Research Journal of Applied Sciences, Engineering and Technology 7(5): 3069-3074, 04 DOI:0.906/rjaset.7.644 ISSN: 040-7459; e-issn: 040-7467 04 Maxwell Scientific Publication Corp. Submitted: August 7,
Sigma J Eng & Nat Sci 36 (1), 2018, Sigma Journal of Engineering and Natural Sciences Sigma Mühendislik ve Fen Bilimleri Dergisi
Sigma J Eng & Nat Sci 36 (1), 2018, 49-62 Sigma Journal of Engineering and Natural Sciences Sigma Mühendislik ve Fen Bilimleri Dergisi Research Article THE ENERGY EFFICIENT CONFIGURATIONS OF NATURAL CONVECTION
COMPARISON 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
6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions
Chapter 1 Earth Science Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions Project Representative Yozo Hamano Authors Ataru Sakuraba Yusuke Oishi
Project #1 Internal flow with thermal convection
Project #1 Internal flow with thermal convection MAE 494/598, Fall 2017, Project 1 (20 points) Hard copy of report is due at the start of class on the due date. The rules on collaboration will be released
Assessment of the Accuracy of the Multiple-Relaxation-Time Lattice Boltzmann Method for the Simulation of Circulating Flows
Mathematical Modelling and Applications 2017; 25): 47-51 http://www.sciencepublishinggroup.com/j/mma doi: 10.11648/j.mma.20170205.11 ISSN: 2575-1786 Print); ISSN: 2575-1794 Online) Assessment of the Accuracy
Computation of turbulent Prandtl number for mixed convection around a heated cylinder
Center for Turbulence Research Annual Research Briefs 2010 295 Computation of turbulent Prandtl number for mixed convection around a heated cylinder By S. Kang AND G. Iaccarino 1. Motivation and objectives
THERMAL PERFORMANCE EVALUATION OF AN INNOVATIVE DOUBLE GLAZING WINDOW
THERMAL PERFORMANCE EVALUATION OF AN INNOVATIVE DOUBLE GLAZING WINDOW Luigi De Giorgi, Carlo Cima, Emilio Cafaro Dipartimento di Energetica, Politecnico di Torino, Torino, Italy Volfango Bertola School
Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
Transport in Porous Media (2006) 64: 1 14 Springer 2006 DOI 10.1007/s11242-005-1126-6 Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
Chapter 7: Natural Convection
7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection
Natural Convection Heat Transfer from a Rectangular Block Embedded in a Vertical Enclosure
Chapter 5 Natural Convection Heat Transfer from a Rectangular Block Embedded in a Vertical Enclosure Xiaohui Zhang Additional information is available at the end of the chapter http://dx.doi.org/10.5772/52666
EXTENDED 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
Simulation 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
Grad 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,
CHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,