Access the full text.
Sign up today, get DeepDyve free for 14 days.
We simulate the two-dimensional ﬂuid ﬂow around National Advisory Committee for Aeronautics (NACA) 0012 airfoil using a hybrid lattice Boltzmann method (HLBM), which combines the standard lattice Boltz- mann method with an unstructured ﬁnite-volume formulation. The aim of the study is to assess the numerical performances and the robustness of the computational method. To this purpose, after providing a conver- gence study to estimate the overall accuracy of the method, we analyze the numerical solution for diﬀerent values of the angle of attack at a Reynolds number equal to 10 . Subsequently, ﬂow ﬁelds at Reynolds numbers up to 10 are computed for a zero angle of attack conﬁguration. A grid reﬁnement scheme is applied to the uniformly spaced component of the overlapping grid system to further enhance the numerical eﬃciency of the model. The results demonstrate the capability of the HLBM to achieve high accuracy near solid curved walls, thus providing a viable alternative in the realm of oﬀ-lattice Boltzmann methods based on body-ﬁtted mesh. 1 Introduction The lattice Boltzmann methods (LBM) [1, 2, 3, 4, 5, 6] are a relatively recent approach to computational ﬂuid dynamics (CFD), which has been proven to be successful in a broad range of applications, from turbulence [7, 8, 9], to multiphase and free-surface ﬂows [10, 11], as well as to non-Newtonian ﬂows , ﬂuid-structure interaction problems , porous media [14, 15, 16] and beyond. The original formulation of LBM is based on uniform cartesian grids. Such a feature makes the algorithm particularly easy to implement and the nu- merical scheme extremely eﬃcient; moreover, LBM algorithms are characterized by an high level of scalability on parallel processing systems. Despite these appealing peculiarities, the uniform space discretization aﬀects the capability of LBM to solve problems involving complex geometries and multi-scale problems. This may constitute an obstacle for many typical practical engineering applications. To overcome this issue and extend the applicability of LBM, many variants of the original LB scheme have been proposed in literature. First chronological appearances of ﬁnite-diﬀerence based schemes are those of Higuera and Succi  and Nannelli and Succi , where LBM is applied on irregular lattices. Numerical implementations on unstructured grids have been ﬁrst introduced by Peng et al. [19, 20] and Xi et al.  who proposed a cell-vertex ﬁnite-volume ap- proach to numerically solve the integral form of the lattice Boltzmann equation. Other ﬁnite-volume techniques are those presented by Patil and Lakshmisha [22, 23, 24], Chen and Schaefer , and some of the authors of the present paper [26, 27, 28, 29, 30]. A diﬀerent strategy to solve the lattice Boltzmann equation on non- uniform and unstructured meshes is based on Galerkin-type ﬁnite element methods and it has been introduced in [31, 32, 33, 34]. Other approaches are the general characteristics based oﬀ-lattice Boltzmann scheme proposed by Bardow et al. , the discrete uniﬁed gas kinetic scheme (DUGKS)  and the semi-lagrangian approach of Kramer et al. . In  a complete comparative study is performed to evaluate the performance of several explicit oﬀ-lattice Boltzmann methods. Although the mentioned schemes enhance signiﬁcantly the geometrical ﬂexibility of LBM, in general, they in- crease the complexity of the method itself, partially losing some key advantages. Moreover, their computational eﬃciency is considerably lower than that of traditional LBM. To overcome this limitation, we recently proposed an hybrid lattice Boltzmann approach (HLBM) based on the coupling between a uniform grid model and an Electronic address: firstname.lastname@example.org; Corresponding author arXiv:2006.10487v1 [physics.flu-dyn] 18 Jun 2020 unstructured body-ﬁtted grid model . The HLBM retains the outstanding advantages of traditional LBM while bringing the ﬂexibility of an unstructured discretization scheme. In this work, we present a numerical investigation on the unsteady behavior of the ﬂuid ﬂow over a stationary NACA 0012 airfoil in the low Reynolds numbers regime, with the aim of further analyze the capabilities of the HLBM. Results demonstrate that the new method is robust and eﬃcient and, for this reason, it can be regarded as a viable approach to ﬂuid dynamic problems involving curved geometries. 2 The hybrid strategy In this section, we brieﬂy recall the principles and the key points of the HLBM. The HLBM is constructed on the overlapping between a uniform cartesian grid, where the traditional LBM is applied, with an unstructured grid made of triangular elements, where a ﬁnite-volume lattice Boltzmann formulation is applied. Both the numerical methods are based on the single-time relaxation approach. The structured uniform grid covers the whole computational domain, while the unstructured one is deﬁned in a conﬁned region, that is, around the solid body. For a distribution function f (x, t), representing the probability of ﬁnding a ﬂuid particle at position x and time t that is moving along the i-th lattice direction with a discrete speed c , the lattice Boltzmann equation reads as follows: Δt s eq f (x + c Δt , t + Δt )− f (x, t) = − [f (x, t)− f (x, t)] , (1) i i s s i i Equation 1 is solved numerically on a uniform Cartesian grid, with τ and Δt being deﬁned as the relaxation s s time and the time step, respectively, related to the ’structured’ scheme. The equilibrium distribution function eq f (x, t) is calculated as follows: 2 2 c · u (x, t) [c · u (x, t)] [u (x, t)] eq i i f (x, t) = w ρ (x, t) 1 + + − (2) 2 4 2 c 2c 2c s s s where c is the lattice speed of sound, the parameters w are a set of weights normalized to unity, and s i ρ (x, t) and u (x, t) are the macroscopic variables, namely ﬂuid density and velocity, respectively. For an ideal incompressible ﬂuid ﬂow, equation 1 reproduces the Navier-Stokes equations, with pressure p = ρc and the 2 Δt kinematic viscosity being deﬁned as ν = c (τ − ). The ﬁnite-volume lattice Boltzmann method adopted on the unstructured grid is a cell-vertex type scheme. The following equation applies at each node P of the grid: 0 0 K K X X Δt u eq f (P, t + Δt ) = f (P, t) + Δt S f (P , t)− C [f (P , t)− f (P , t)]. (3) i u i u ik i k ik i k k k=0 k=0 Details for derivation of such equation are provided in [26, 39], therefore are not reported here. In equation 3, τ and Δt are the relaxation time and the time step, respectively, related to the ’unstructured’ scheme, k = 0 u u denotes the pivotal node P and the summations run over the nodes P connected to P ; the quantities S and k ik C represent the streaming and collisional matrices of the i-th population related to the k-th node, respectively; ik the equilibrium distribution function is deﬁned by equation (2). The theoretical kinematic viscosity associated with this scheme is ν = c τ . Within the hybrid framework, equations (1) and (3) are solved consistently in time, with the exchange of information taking place at predeﬁned interpolation nodes, which are placed over the overlapping grids region, in terms of both macroscopic variables and distribution functions. The following equations describe the post- collision for the structured and unstructured interpolation nodes, respectively: 2 (τ − Δt ) eq,u s s neq,u f = f + f (4) i i i 2τ − Δt s s 0 0 K K X X Δt Δt s u eq,s neq,s eq,? neq,? neq,? f = f + 1− f + Δt S f + f − C f (5) u ik ik i i i ik ik ik 2τ τ s u k=0 k=0 neq where superscripts s and u refer to the ’structured’ and the ’unstructured’ nodes, respectively, and f is eq,? neq,? the non-equilibrium distribution function. The quantities f and f in the summation terms of equation ik ik (5) are deﬁned as follows: Δt eq,? eq,s neq,? neq,s f = f , f = 1− f (6) ik ik ik ik 2τ 2 or eq,? eq,u neq,? neq,u f = f , f = f (7) ik ik ik ik depending on whether the k-th node is an interpolation node or not, respectively. The distribution functions eq,u neq,u eq,s neq,s f and f of the right-hand side of equation (4) and f and f of the right-hand side of equation i i i i (5) are evaluated by interpolation procedure. 3 The interpolation scheme For traditional multi-domain grid reﬁnement, the information communication on grid transitions signiﬁcantly aﬀects the accuracy of the numerical scheme . Similarly, in the HLBM, the exchange of information between structured and unstructured grid is of crucial importance from the accuracy point of view. Here we describe the interpolation scheme used for reconstructing both the macroscopic variables and the distribution functions from one grid to the other. Later, we will provide a convergence study, in order to assess the overall accuracy of the present HLBM. The hybrid grid presents interpolation nodes of two kinds: structured and unstructured. The quantities on the structured interpolation nodes are reconstructed by considering four triangular elements of the unstructured mesh, as depicted in Figure 1. φ(x, y) Figure 1: Interpolation scheme representation for the structured nodes. The generic variable φ(x, y) is approx- imated by a linear combination of the shape functions deﬁned on the six unstructured donor nodes. A set of 2nd order shape functions is deﬁned on the six unstructured nodes surrounding the structured interpolation node. Therefore, the polygon composed of the four selected triangles represents a quadratic element, and the generic variable φ(x, y) on the structured node is approximated by a linear combination of the shape functions ψ , as follows: φ(x, y) = φ ψ (x, y). (8) k k k=1 Similarly, in order to reconstruct the quantities on the unstructured interpolation nodes, we select a square region, which is composed of four structured elements, as depicted in Figure 2. 3 φ φ φ 7 8 9 φ φ φ 4 5 6 φ(x, y) φ φ φ 1 2 3 Figure 2: Interpolation scheme representation for the unstructured nodes. The generic variable φ(x, y) is approximated by a linear combination of the shape functions deﬁned on the nine structured donor nodes. On the nine structured nodes surrounding the unstructured interpolation node we deﬁne a set of 2nd order shape functions. Therefore, the generic variable φ(x, y) on the unstructured node is evaluated by means of a linear combination of these shape functions, as follows: φ(x, y) = φ ψ (x, y). (9) k k k=1 4 Problem setup In this study, we simulate a laminar ﬂow past a NACA 0012 airfoil, which is characterized by a symmetric proﬁle whit 12% thickness to chord length ratio. NACA airfoils have been extensively used for the validation of numerical schemes, thanks to the availability of experimental data for several diﬀerent proﬁles. Here, with the aim of analyzing the inﬂuence of the angle of attack on the development and the evolution of vortices from the airfoil surface, a parametric study is conducted at Re = 1000 for a wide range of angles of attack. Moreover, ﬂow ﬁelds are computed at higher values of Reynolds number, speciﬁcally Re = 2000, 5000 and 10000, at zero angle of attack. The geometry of the two-dimensional computational domain adopted for the present analysis is illustrated in Figure 3 where C indicates the chord length. The airfoil is placed symmetrically with respect to bottom and top boundaries. Figure 3: Domain conﬁguration for simulation of ﬂuid ﬂow around NACA 0012 airfoil. In order to further enhance the computational eﬃciency of the hybrid method, we implement the grid reﬁne- ment technique of Filippova and Hanel  within the structured sub-model scheme. In particular, we consider four structured reﬁnement levels. Each structured level is constituted by a rectangular box, which represents 4 a resolution region in which the structured mesh has uniform discretization. The grid reﬁnement is performed by dividing the lattice spacing by a factor of two, with the inner boxes having higher resolution over the outer box in which they are contained. The innermost structured region, which is characterized by a unitary lattice spacing, overlaps the unstructured mesh surrounding the airfoil. No-slip wall boundary condition is applied at the airfoil surface, uniform velocity and zero pressure gradient are set at the inlet, ﬁxed pressure value and zero velocity gradient at the outlet, and open-boundary conditions are applied to the top and bottom walls of the domain. The unstructured mesh surrounding the solid body is constituted by an ellipse. This choice is purely arbitrary, and it has been taken after testing other elementary shapes and assessing the non-dependency of the results from the unstructured mesh conﬁguration. Figure 4 shows a detail for a generic conﬁguration of the selected unstructured discretization around the airfoil. In the ﬁgure, the overlapping area between the ﬁnest structured grid and the unstructured grid is represented as the region enclosed by two red lines. The interpolation nodes, which are responsible of the exchange of information between the two grid structures, are placed along this region. Also the internal ﬁctitious curve deﬁning such an interface region is constituted by an ellipse. Figure 4: Detail of the hybrid mesh around NACA 0012 airfoil. The interface region is deﬁned as the region between the two red curves. 5 Convergence study In order to estimate the overall accuracy of the numerical method and to verify that the adopted interpolation scheme is compatible with the order of accuracy of the standard LBM we preliminary perform a convergence study. For this analysis we consider ﬁve hybrid meshes of diﬀerent resolution. In particular, the following chord lengths are chosen: C = 48, 64, 96, 128, 144. The average element-size at wall is the same for all the unstructured meshes. Considering the Reynolds number based on the free-stream velocity and the airfoil chord, the analysis is carried out at Re = 500. The angle of attack of the airfoil is set to α = 0 . All the simulations are performed with same time-step for the unstructured routine. The results are provided in terms of percentage error for the drag coeﬃcient, with reference to the value computed by XFOIL - Subsonic Airfoil Development System. XFOIL is an interactive program for the design and analysis of subsonic isolated airfoils, developed by Mark Drela at Massachusetts Institute of Technology . In Figure 5 we show the deviation of the computed mean drag coeﬃcient from the reference value. 5 100 numerical results ﬁtting 2nd-order slope 0.1 30 60 90 120 150 180 210 C [l.u.] Figure 5: Convergence of the HLBM. The numerical results are represented by open circles; the solid thick line is the ﬁtting curve. The second-order slope is also provided as a guide (dashed line). We notice that the order of convergence of the present scheme is about equal to 2. This result indicates that the interpolation scheme preserves the second order accuracy of the standard LBM. 6 Results and discussion 6.1 Eﬀect of the angle of attack ◦ ◦ Simulations are performed for a range of angles of attack, α, between 0 and 29 . All the cases are studied at the ﬁxed value of Re = 1000. For this analysis, the chord length of the airfoil is set to be equal to 512 lattice units. For each hybrid mesh, the total number of structured nodes is approximately equal to 3.3· 10 , while the number of unstructured nodes is signiﬁcantly lower, that is about 2.2· 10 . Figure 6 illustrates the instantaneous velocity and pressure ﬁelds for Re = 1000 at several angle of attack, ◦ ◦ ◦ ◦ ◦ ◦ namely: α= 0 , 8 , 10 , 20 and 28 . For angles of attack lower than 4 , no ﬂow separation is observed. However, beyond this particular value of α, a recirculation region arise at the rear of the airfoil. The unsteady ◦ ◦ vortex shedding phenomenon is observed at about α = 8 , as for α < 8 the numerical solution converges to steady state. This evidence is in agreement with results available in literature [43, 44, 45, 46]. A further increasing of the angle of attack gives rise to leading-edge vortexes and ﬂow instabilities at the rear of the airfoil become larger. Then, starting from an angle of attack of 24 degrees, the vortex shedding becomes non-periodic. err [%] D (a) (b) (c) (d) (e) (f) (g) (h) (i) (j) ◦ ◦ Figure 6: Instantaneous velocity and pressure ﬁelds at Re = 1000: (a),(b): α = 0 ; (c),(d): α = 8 ; (e),(f): ◦ ◦ ◦ α = 10 ; (g),(h): α = 20 ; (i),(j): α = 28 . Figure 7 shows the instantaneous streamline patterns for the ﬂuid ﬂow around the airfoil, at 4, 10 and 28 degrees, respectively. 7 (a) (b) (c) Figure 7: Instantaneous streamline patterns of the ﬂow around NACA0012 airfoil at Re = 1000: (a): α = 4 ; ◦ ◦ (b): α = 10 ; (c): α = 28 . We note that, for α = 4 , ﬂow separation occurs and a recirculation area is observable at the rear of the airfoil. From Figure 7, it is also evident the shedding of the trailing-edge vortexes for α = 10 . As the angle of attack is increased a large-scale ﬂow separation on the surface of the airfoil takes place and the ﬂow becomes highly unsteady. This behavior is already visible at α = 28 . To asses the aerodynamic performance of the airfoil, drag and lift coeﬃcients are calculated as follows: C = , (10) 1 2 ρ u C C = , (11) ρ u C where the quantities with subscript 0 refer to the far ﬁeld conditions, while F and F are the drag and D L lift force, respectively, namely the horizontal and the vertical components of the aerodynamic force acting on the solid body. Moreover, we quantify the main features of the ﬂow by computing the mean pressure coeﬃcient and the skin friction coeﬃcient as follows: 8 p− p C = , (12) ρ u C = , (13) ρ u 2 0 where τ denotes the mean wall shear stress. We emphasize that the HLBM allows a direct computation of the forces on the surface of the body. In fact, thanks to the presence of the unstructured body-ﬁtted mesh, both pressure and stress tensor are locally available at the computational nodes deﬁning the geometry of the body. This implies that, to compute the aerodynamics variables, interpolation routines are not needed to be implemented. In Figures 8 and 9 we show the mean pressure coeﬃcient and the mean skin friction coeﬃcient, respectively, at α = 8 , along the non-dimensional chord C. The results are compared with those available from literature. 1.5 present Kurtulus  1.0 0.5 p 0.0 −0.5 −1.0 −1.5 0 0.2 0.4 0.6 0.8 1 Figure 8: Mean pressure coeﬃcient distribution over the upper and lower surfaces of the NACA 0012 airfoil at α = 8 : comparison between present results and those reported in . 0.40 present 0.35 Kurtulus  0.30 0.25 0.20 0.15 0.10 0.05 0.00 −0.05 0 0.2 0.4 0.6 0.8 1 Figure 9: Mean skin friction coeﬃcient distribution over the upper and lower surfaces of the NACA 0012 airfoil at α = 8 : comparison between present results and those reported in . 9 Despite a slight overestimation of the maximum value of the mean skin friction coeﬃcient on both lower and upper surfaces of the airfoil, the overall agreement between our computed trends and those from literature is satisfactory. Also, the skin friction distribution provides a meaningful way to determine the separation point. In fact, this can be inferred from the location of zero skin friction with a negative gradient. The obtained result for C indicates that the location of the separation point on the upper surface of the airfoil is approximately equal to 0.41C, which is in line with the value provided by Kurtulus  (0.4C). In Figures 10 and 11 the time-averaged drag and lift coeﬃcients are reported as a function of the angle of attack, respectively. The ﬁgures shows also a comparison with literature data. 1.4 present study Liu et al.  1.2 Kurtulus  Khalid and Akhtar  0.8 0.6 0.4 0.2 0 5 10 15 20 25 30 α[ ] Figure 10: Lift coeﬃcient as a function of the angle of attack, at Re = 1000. Present results are compared with those from the numerical works of Liu et al.  (triangles), Kurtulus  (squares) and Khalid and Akhtar.  (stars). present study Kurtulus  0.8 0.6 0.4 0.2 0 5 10 15 20 25 30 α[ ] Figure 11: Drag coeﬃcient as a function of the angle of attack, at Re = 1000. Present results are compared with the numerical ﬁndings of Kurtulus . Overall, a good agreement between the solution obtained by HLBM and values from literature is achieved. In particular, the lift coeﬃcient calculated by Liu et al. seems to be slightly overestimated at higher angles of attacks while the values computed by Kurtulus and by Khalid and Akhtar oscillates after a critical values is reached. 10 The phenomenon of stall associated with airfoils is caused by a massive ﬂow separation, which leads to a sudden drop in lift when the angle of attack is slightly increased beyond a critical value. This behavior seems to be properly captured by our numerical simulations. In particular, in this analysis, the numerical results suggest that stall occurs for a critical value of α around 26 , corresponding to a sharp drop in lift forces acting on the airfoil. Moreover, Figure 11 shows that the drag coeﬃcient remains fairly constant until a particular value of angle of attack is reached. This value corresponds, approximatively, to α = 8 , when vortex shedding starts to appear. For angles of attack greater than α = 8 , the drag force increases with α, as expected. For the drag coeﬃcient, we observe a good agreement between our ﬁndings and those available in the literature in the entire range of angles of attack analyzed. However, we notice a slight discrepancy at α = 29 . Our analysis suggests that, for this particular case, the drag coeﬃcient is close to that at 26 , corresponding to the critical value, while Kurtulus obtains a signiﬁcant higher value. This diﬀerence is however acceptable, and can be attributed to the irregular behavior of the ﬂow in this regime. Figure 12 depicts the time history of the lift coeﬃcient. 2.5 ◦ ◦ ◦ ◦ α= 4 α= 10 α= 16 α= 28 1.5 0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 t[l.u.] ·10 ◦ ◦ ◦ ◦ Figure 12: Time-history of lift coeﬃcient for α = 4 , 10 , 16 and 28 . The numerical results demonstrate that, below an angle of attack of 8 degrees, the lift coeﬃcient converges to a constant value. This behavior corresponds to a steady state solution, since no ﬂuctuations in the wake of the airfoil occur. However, due to the occurring of a trailing-edge vortex shedding, at α = 8 the lift coeﬃcient begins to oscillate periodically. Above this particular value of α the ﬂow becomes unsteady and the amplitude of the lift coeﬃcient diﬀers from zero. As a matter of example, in Figure 12, for α equal to 10 and 16 degrees, the lift coeﬃcient is shown to exhibit a periodic trend. Moreover, it is possible to note that, as the angle of attack is increased, the amplitude of oscillation of the lift coeﬃcient becomes larger. The angle of attack for which the stall phenomenon occurs is known as stall angle. At this critical value, the developing of leading-edge vortices on the surface of airfoil and the massive ﬂow separation lead to irregular oscillations of lift coeﬃcient. Finally, the aerodynamic characteristic C /C against the angle of attack is plotted in Figure 13. L D 11 3 2.5 1.5 present study 0.5 Kurtulus  0 5 10 15 20 25 30 α[ ] Figure 13: Lift-to-drag ratio vs angle of attack α: comparison between present results and those reported in . From the analysis, it emerges that the maximum aerodynamic eﬃciency occurs between 10 and 12 degrees, since the ratio C /C reaches a pick. Present computational results perfectly ﬁt numerical literature data . L D 6.2 Eﬀect of the Reynolds number In this section we discuss the ﬂow ﬁelds computed for Re = 2000, 5000 and 10000 at an angle of attack of α = 0 . The same computational domain conﬁguration illustrated in Figure 3 is employed, while the chord length is taken to be equal to 1024 lattice units, referring to the ﬁnest structured reﬁnement level. In this case, the 6 4 hybrid mesh is composed of approximately 12.6· 10 structured nodes and 8.8· 10 unstructured nodes. Figures 14 and 15 show velocity and pressure ﬁelds for Re = 10000. Figure 14: Instantaneous velocity ﬁeld at Re = 10000. C /C L D Figure 15: Instantaneous pressure ﬁeld at Re = 10000. In order to validate the solution, a comparison between the results obtained by HLBM simulations and those calculated from XFOIL is made. Table 1 reports the values for the drag coeﬃcient as a function of Reynolds number. Re HLBM XFOIL 1000 0.119 0.119 2000 0.084 0.084 5000 0.052 0.054 10000 0.037 0.040 Table 1: Drag coeﬃcients obtained by HLBM in comparison with values calculated by XFOIL for NACA0012 airfoil at α = 0. We note that the two solutions agree with each other very well for the entire range of Reynolds number investigated. 6.3 Numerical performances In order to estimate the gain in eﬃciency of the HLBM with respect to the standard LBM, we consider an hybrid mesh which is composed of a uniformly spaced region and the same unstructured mesh employed in the analysis described in 6.1. Also, we consider the same domain size and set of numerical parameters. A practical way to quantify the diﬀerence of computational cost between HLBM and LBM is provided by the following deﬁnition of speedup (gain) indicator, introduced in : T 1 + χs n G = = . (14) T g eq In equation 14, T is the wall-clock time per algorithm iteration associated to the HLBM, while T represents the h eq wall-clock time per algorithm iteration with reference to an equivalent LBM implementation on uniform spaced mesh, that is a standard LBM with similar resolution around the solid body. The parameter χ in equation 14 indicates the ratio of the number of unstructured nodes to that of structured ones, g is the ratio of the number of structured nodes which would be required by an equivalent LBM to the number of structured nodes actually used in the HLBM, while s represents the ratio of the processing speed for the standard LBM to the speciﬁc processing speed of the unstructured method. The unstructured ﬁnite-volume scheme is computationally more demanding than the standard LBM one, with typical values of s ranging between 3 and 5, as conﬁrmed also in the present study. By considering that the resolution of the unstructured mesh in proximity of the solid −4 body is around 0.8 lattice units, we ﬁnd g ' 1.6. Also, for the present hybrid mesh, we ﬁnd χ ' 1.5 · 10 . Moreover, in all the performed simulations, we set the number of sub-iterations required by the unstructured method as n = 10. By setting such a set of constants in equation 14, we obtain G ' 0.63. This value indicates that, for the case under consideration, the hybrid method runs faster than an equivalent standard LBM since, by deﬁnition, a value of G less than 1 denotes the regime where the computational cost of HLBM is lower than that of standard LBM, when a similar resolution around the body is considered. 13 To conclude, the numerical eﬃciency of the HLBM is, in general, higher than the one corresponding to the standard LBM. This is due to the presence of the unstructured mesh that, despite being computationally more expensive to solve than the uniform spaced one which is inherent to standard LBM, represents a small portion of the overall discretized domain. Therefore, the unstructured method provides a very ﬂexible strategy for local grid reﬁnement with major potential for multiscale applications. 7 Conclusions The HLBM capability to properly predict the ﬂuid ﬂow behavior of a complex case study, such as that over an airfoil in static stall, has been assessed through this work. The analysis, performed for a NACA 0012 airfoil at relatively low Reynolds numbers and diﬀerent angles of attack, shows that the hybrid method is able to provide accurate results. Boundary layer separation, static stall, as well as the other physical phenomena involved, were captured by the numerical simulations. In particular, we found that, for a Reynolds number equal to 10 , the vortex shedding regime starts at an angle of attack of 8 , while the critical angle of attack, for which the stall phenomenon occurs, is around 26 . For angles of attack greater than this value, we observe a sharp drop in lift forces acting on the airfoil. Such results are in line with those available in literature. Overall, the involved ﬂuid ﬂow phenomena seems to be well predicted and the computation of drag and lift forces acting on the body, in the whole range of parameters analyzed, leads to satisfactory results. This achievement strengthens the robustness of HLBM thus projecting it as a potential tool for dealing with complex geometries. In fact, the presence of a body-conformed unstructured mesh allows high accuracy solution of the boundary layer regions, since a precise deﬁnition of the solid body proﬁle as well as a ﬂexible reﬁning of the discretized domain is achievable. Moreover, thanks to the feature of LBM of allowing a local computation of the stress tensor at nodes of the computational domain, the HLBM prevents from the need of implementing boundary-ﬁtting procedures for the computation of hydrodynamic forces on solid bodies. This property, which is partially lost in traditional lattice Boltzmann approaches, is here fully retained. Further, the method presents high computational eﬃciency, since the unstructured model is conﬁned to a relatively small portion of the computational domain. These features, exploited throughout this study, are instrumental to the implementation of schemes for solving multi-scale problems. The HLBM can be regarded also as an interesting alternative to traditional computational ﬂuid dynamic approaches as well as to other oﬀ-lattice Boltzmann methods to the simulation of turbulent ﬂows. In fact, to properly capture all the main features of a fully developed turbulent ﬂow, high resolution is required in those region of the ﬂuid domain which are characterized by the smallest scales. On the other hand, the numerical eﬃciency of the method plays a key role, therefore the geometric ﬂexibility of a non-uniform mesh is often desirable. To extend the applicability of the HLBM to turbulent ﬂows, the implementation of a three- dimensional scheme is required. Extensions of the unstructured method used in this work are already present in literature . Therefore, further studies will be aimed to develop suitable interpolation schemes for three- dimensioanl interfaces. Acknowledgments The numerical simulations were performed on Zeus HPC facility, at the University of Naples "Parthenope"; Zeus HPC has been realized through the Italian Government Grant PAC01_00119 MITO - Informazioni Multime- diali per Oggetti Territoriali, with Prof. Elio Jannelli as the Scientiﬁc Responsible. One of the authors (Sauro Succi) was partially supported by the European Research Council under the Euro- pean Union’s Horizon 2020 Framework Programme (No. FP 2014-2020) ERC Grant Agreement No. 739964 (COPMAT). 14 References  S. Succi, The Lattice Boltzmann Equation for Fluid Dynamics and Beyond, Clarendon Press, 2001.  S. Chen, G. D. Doolen, Lattice Boltzmann method for ﬂuid ﬂows, Annu. Rev. Fluid Mech. 30 (1998) 329–364. doi:10.1146/annurev.fluid.30.1.329.  S. Succi, Lattice Boltzmann 2038, EPL 109 (2015) 50001. doi:10.1209/0295-5075/109/50001.  R. Benzi, S. Succi, M. Vergassola, The Lattice Boltzmann Equation: Theory and Applications, Phys. Rep. 222 (1992) 145–197. doi:10.1016/0370-1573(92)90090-M.  F. Higuera, S. Succi, Simulating the ﬂow around a circular cylinder with a Lattice Boltzmann Equation, Europhys. Lett. 8 (1989) 517–521.  G. McNamara, G. Zanetti, Use of the Boltzmann Equation to simulate Lattice-Gas Automata, Phys. Rev. Lett. 61 (1988) 2332–2335.  I. Karlin, F. BÃ¶sch, S. Chikatamarla, Gibbs principle for the lattice-kinetic theory of ﬂuid dynamics, Phys. Rev. E 90 (2014) 031302(R). doi:10.1103/PhysRevE.90.031302.  F. BÃ¶sch, S. Chikatamarla, I. Karlin, Entropic multirelaxation lattice Boltzmann models for turbulent ﬂows, Phys. Rev. E 92 (2015) 043309. doi:10.1103/PhysRevE.92.043309.  B. Dorschner, S. Chikatamarla, I. Karlin, Transitional ﬂows with the entropic lattice Boltzmann method, J. Fluid Mech. 824 (2017) 388–412. doi:10.1017/jfm.2017.356.  G. Falcucci, S. Ubertini, D. Chiappini, S. Succi, Modern lattice Boltzmann methods for multiphase mi- croﬂows, IMA Journal of Applied Mathematics 76 (5) (2011) 712–725. doi:10.1093/imamat/hxr014.  S. D. Francesco, C. Biscarini, P. Manciola, Numerical simulation of water free-surface ﬂows through a front- tracking lattice Boltzmann approach, Journal of Hydroinformatics 17 (2015) 1–6. doi:10.2166/hydro. 2014.028.  G. Di Ilio, D. Chiappini, G. Bella, A comparison of numerical methods for non-Newtonian ﬂuid ﬂows in a sudden expansion, Int. J. Mod. Phys. C 27 (2016) 1650139. doi:10.1142/S0129183116501394.  A. De Rosis, G. Falcucci, S. Ubertini, F. Ubertini, Aeroelastic study of ﬂexible ﬂapping wings by a cou- pled lattice Boltzmann-ﬁnite element approach with immersed boundary method, Journal of Fluids and Structures 49 (2014) 516–533. doi:10.1016/j.jfluidstructs.2014.05.010.  D. Chiappini, Numerical simulation of natural convection in open-cells metal foams, International Journal of Heat and Mass Transfer 117 (2018) 527–537. doi:10.1016/j.ijheatmasstransfer.2017.10.022.  D. Chiappini, G. Bella, A. Festuccia, A. Simoncini, Direct Numerical Simulation of an Open-Cell Metallic Foam through Lattice Boltzmann Method, Communications in Computational Physics 18 (3) (2015) 707– 722. doi:10.4208/cicp.191114.270315a.  A. Zarghami, S. Di Francesco, C. Biscarini, Porous substrate eﬀects on thermal ﬂows through a rev-scale ﬁnite volume Lattice Boltzmann Model, International Journal of Modern Physics C 25 (2014) 1350086. doi:10.1142/S0129183113500861.  F. Higuera, S. Succi, R. Benzi, Lattice gas dynamics with enhanced collisions, Europhys. Lett. 9 (1989) 345–349.  F. Nannelli, S. Succi, The Lattice Boltzmann Equation on Irregular Lattices, J. Stat. Phys. 68 (1992) 401–407.  G. Peng, H. Xi, C. Duncan, S. Chou, Lattice Boltzmann method on irregular meshes, Phys. Rev. E 58 (1998) R4124–R4127.  G. Peng, H. Xi, C. Duncan, S. Chou, A ﬁnite volume scheme for the Lattice Boltzmann Method on unstructured meshes, Phys. Rev. E 59 (1999) 4675–4682.  H. Xi, G. Peng, S. Chou, Finite-volume lattice Boltzmann method, Phys. Rev. E 59 (1999) 6202–6205.  D. Patil, K. Lakshmisha, Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh, Journal of Computational Physics 228 (2009) 5262-5279. doi:10.1016/j.jcp.2009.04.008. 15  D. Patil, K. Lakshmisha, Two-dimensional ﬂow past circular cylinders using ﬁnite volume lattice Boltzmann formulation, Int. J. Numer. Meth. Fluids 69 (2012) 1149–1164. doi:10.1002/fld.2637.  D. Patil, Chapman-Enskog analysis for ﬁnite-volume formulation of lattice Boltzmann equation, Physica A 392 (2013) 2701-2712. doi:10.1016/j.physa.2013.02.016.  L. Chen, L. Schaefer, A uniﬁed and preserved Dirichlet boundary treatment for the cell-centered ﬁnite volume discrete Boltzmann method, Physics of Fluids 27 (2015) 027104. doi:10.1063/1.4907782.  S. Ubertini, G. Bella, S. Succi, Lattice Boltzmann method on unstructured grids: Further developments, Phys. Rev. E 68 (2003) 016701. doi:10.1103/PhysRevE.68.016701.  S. Ubertini, S. Succi, G. Bella, Lattice Boltzmann schemes without coordinates, Phil. Trans. R. Soc. Lond. A 362 (2004) 1763–1771.  S. Ubertini, G. Bella, S. Succi, Unstructured Lattice Boltzmann equation with memory, Math. Comp. Sim. 72 (2006) 237–241.  A. Zarghami, C. Biscarini, S. Succi, S. Ubertini, Hydrodynamics in porous media: a ﬁnite volume lattice Boltzmann study, Journal of Scientiﬁc Computing 59 (2014) 80–103. doi:doi.org/10.1007/ s10915-013-9754-4.  A. Zarghami, S. Ubertini, S. Succi, Finite-volume lattice Boltzmann modeling of thermal transport in nanoﬂuids, Computers and Fluids 77 (2013) 56–65. doi:10.1016/j.compfluid.2013.02.018.  T. Lee, C. Lin, A Characteristic Galerkin Method for Discrete Boltzmann Equation, J. Comp. Physics 171 (2001) 336–356. doi:10.1006/jcph.2001.6791.  T. Lee, C. Lin, An Eulerian description of the streaming process in the lattice Boltzmann equation, J. Comp. Physics 185 (2003) 445–471. doi:10.1016/S0021-9991(02)00065-7.  M. Min, T. Lee, A spectral-element discontinuous Galerkin lattice Boltzmann method for nearly incom- pressible ﬂows, J. Comp. Physics 230 (2011) 245–259. doi:10.1016/j.jcp.2010.09.024.  S. Patel, T. Lee, A new splitting scheme to the discrete Boltzmann equation for non-ideal gases on non- uniform meshes, J. Comp. Physics 327 (2016) 799–809. doi:10.1016/j.jcp.2016.09.060.  A. Bardow, I. Karlin, A. Gusev, General characteristic-based algorithm for oﬀ-lattice Boltzmann simula- tions, Europhys. Lett. 75 (2006) 434-440. doi:10.1209/epl/i2006-10138-1.  L. Zhu, P. Wang, Z. Guo, Performance evaluation of the general characteristics based oﬀ-lattice Boltzmann scheme and DUGKS for low speed continuum ﬂows, Journal of Computational Physics 333 (2017) 227-246. doi:10.1016/j.jcp.2016.11.051.  A. KrÃ€mer, K. KÃŒllmer, D. Reith, W. Joppich, H. Foysi, Semi-Lagrangian oﬀ-lattice Boltzmann method for weakly compressible ﬂows, Phys. Rev. E 95 (2017) 023305. doi:10.1103/PhysRevE.95.023305.  P. Rao, L. Schaefer, Numerical stability of explicit oﬀ-lattice Boltzmann schemes: A comparative study, Journal of Computational Physics 285 (2015) 251–264. doi:10.1016/j.jcp.2015.01.017.  G. Di Ilio, D. Chiappini, S. Ubertini, G. Bella, S. Succi, Hybrid lattice Boltzmann method on overlapping grids, Phys. Rev. E 95 (2017) 013309. doi:10.1103/PhysRevE.95.013309.  D. Lagrava, O. Malaspinas, J. Latt, B. Chopard, Advances in multi-domain lattice Boltzmann grid reﬁne- ment, Journal of Computational Physics 231 (2012) 4808–4822. doi:10.1016/j.jcp.2012.03.015.  O. Filippova, D. Hänel, Grid Reﬁnement for Lattice-BGK Models, J. Comp. Physics 147 (1998) 219–228.  M. Drela, XFOIL: An Analysis and Design System for Low Reynolds Number Airfoils, Proceeding of the Conference: Low Reynolds Number Airfoil Aerodynamics, Notre Dame, Indiana, USA 5-7 (1989) 1–12. doi:10.1007/978-3-642-84010-4_1.  D. Kurtulus, On the unsteady behavior of the ﬂow around NACA 0012 airfoil with steady external condi- tions at Re = 1000, Int. J. Micro Air Veh. 7 (2015) 301–326.  R. Huang, J. Wu, J. Jeng, R. Chen, Surface ﬂow and vortex shedding of an impulsively started wing, J. Fluid Mech. 441 (2001) 265–292. 16  K. Ohmi, M. Coutanceau, O. Daube, T. Loc, Further experiments on vortex formation around an oscillating and translating airfoil at large incidences, J. Fluid Mech. 225 (1991) 607–630.  Y. Liu, K. Li, J. Zhang, H. Wang, L. Liu, Numerical bifurcation analysis of static stall of airfoil and dynamic stall under unsteady perturbation, Commun. Nonlinear Sci. Numer. Simulat. 17 (2012) 3427–3434.  M. Khalid, I. Akhtar, Characteristics of ﬂow past a symmetric airfoil at low reynolds number: a nonlinear perspective, 2012 ASME International Mechanical Engineering Congress and Exposition 7 (2012) 167–175. doi:10.1115/IMECE2012-87389.  N. Rossi, S. Ubertini, G. Bella, S. Succi, Unstructured Lattice Boltzmann method in three dimensions, Int. J. Numer. Meth. Fluids 49 (2005) 619–633. doi:10.1002/fld.1018.
Physics – arXiv (Cornell University)
Published: Jun 18, 2020
Access the full text.
Sign up today, get DeepDyve free for 14 days.