Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lattice Boltzmann”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗

A Unified Theory of Non-Ideal Gas Lattice Boltzmann Models

A non-ideal gas lattice Boltzmann model is directly derived, in an a priori fashion, from the Enskog equation for dense gases. The model is rigorously obtained by a systematic procedure to discretize the Enskog equation (in the presence of an external force) in both phase space and time. The lattice Boltzmann model derived here is thermodynamically consistent and is free of the defects which exist in previous lattice Boltzmann models for non-ideal gases. The existing lattice Boltzmann models for non-ideal gases are analyzed and compared with the model derived here.

Luo, Li-Shi↗

Lattice Boltzmann model for compressible fluids

A lattice Boltzmann model is derived which simulates compressible fluids. By choosing the parameters of the equilibrium distribution appropriately, the sound speed (which may be set arbitrarily low), bulk viscosity, and kinematic viscosity can be selected. This model simulates compressible flows and can include shocks. With a proper rescaling and zero-sound speed, this model simulates Burgers's equation. The viscosity determined by a Chapman-Enskog expansion compares well with that measured form simulations. The exact solutions of Burgers's equation on the unit circle are compared to solutions of lattice Boltzmann model finding reasonable agreement.

Alexander, F. J.↗

Connection Between the Lattice Boltzmann Equation and the Beam Scheme

In this paper we analyze and compare the lattice Boltzmann equation with the beam scheme in details. We notice the similarity and differences between the lattice Boltzmann equation and the beam scheme. We show that the accuracy of the lattice Boltzmann equation is indeed second order in space. We discuss the advantages and limitations of lattice Boltzmann equation and the beam scheme. Based on our analysis, we propose an improved multi-dimensional beam scheme.

Xu, Kun↗

NASA/ONERA Collaboration on Small Hovering Rotor Broadband Noise Prediction using Lattice-Boltzmann and Structured Navier-Stokes Solvers

This work compares two lattice-Boltzmann method solvers, PowerFLOW and ProLB, and two structured Navier-Stokes solvers, OVERFLOW2 and FAST, used by NASA and ONERA, respectively, for the broadband noise prediction of an ideally twisted rotor as part of Implementing Arrangement number FR-0685-0, ‘Comparing Computational Fluid Dynamics Solvers for Broadband Noise Prediction.’ Predicted results are evaluated against measured data from both smooth and rough blade sets acquired in the Small Hover Anechoic Chamber at the NASA Langley Research Center. Aerodynamic thrust predictions are seen to agree more favorably with the rough-blade measurements, whereas torque predictions agree better with the smooth-blade measurements. A tonal noise comparison shows better agreement to the measured data with the two structured Navier-Stokes solvers than with the two lattice-Boltzmann solvers, which is thought to be caused by the different geometric discretization associated with each solver paradigm. Broadband noise comparisons show that both lattice-Boltzmann method solvers trend well with the smooth-blade measurements, with the exception of an overprediction by ProLB between 4 kHz and 15 kHz. OVERFLOW2 is seen to capture the measured nondeterministic tonal content between 3 kHz and 8 kHz on a narrowband spectral basis and FAST agrees well with the rough blades on a one-third octave band basis.

Christopher S. Thurman↗

A Lattice Boltzmann Method for Turbomachinery Simulations

Lattice Boltzmann (LB) Method is a relatively new method for flow simulations. The start point of LB method is statistic mechanics and Boltzmann equation. The LB method tries to set up its model at molecular scale and simulate the flow at macroscopic scale. LBM has been applied to mostly incompressible flows and simple geometry.

Hsu, A. T.↗

Theory of the Lattice Boltzmann Equation: Symmetry properties of Discrete Velocity Sets

In the lattice Boltzmann equation, continuous particle velocity space is replaced by a finite dimensional discrete set. The number of linearly independent velocity moments in a lattice Boltzmann model cannot exceed the number of discrete velocities. Thus, finite dimensionality introduces linear dependencies among the moments that do not exist in the exact continuous theory. Given a discrete velocity set, it is important to know to exactly what order moments are free of these dependencies. Elementary group theory is applied to the solution of this problem. It is found that by decomposing the velocity set into subsets that transform among themselves under an appropriate symmetry group, it becomes relatively straightforward to assess the behavior of moments in the theory. The construction of some standard two- and three-dimensional models is reviewed from this viewpoint, and procedures for constructing some new higher dimensional models are suggested.

Rubinstein, Robert↗

Multiple-Relaxation-Time Lattice Boltzmann Models in 3D

This article provides a concise exposition of the multiple-relaxation-time lattice Boltzmann equation, with examples of fifteen-velocity and nineteen-velocity models in three dimensions. Simulation of a diagonally lid-driven cavity flow in three dimensions at Re=500 and 2000 is performed. The results clearly demonstrate the superior numerical stability of the multiple-relaxation-time lattice Boltzmann equation over the popular lattice Bhatnagar-Gross-Krook equation.

dHumieres, Dominique↗

Simulation of gaseous diffusion in partially saturated porous media under variable gravity with lattice Boltzmann methods

Liquid distributions in unsaturated porous media under different gravitational accelerations and corresponding macroscopic gaseous diffusion coefficients were investigated to enhance understanding of plant growth conditions in microgravity. We used a single-component, multiphase lattice Boltzmann code to simulate liquid configurations in two-dimensional porous media at varying water contents for different gravity conditions and measured gas diffusion through the media using a multicomponent lattice Boltzmann code. The relative diffusion coefficients (D rel) for simulations with and without gravity as functions of air-filled porosity were in good agreement with measured data and established models. We found significant differences in liquid configuration in porous media, leading to reductions in D rel of up to 25% under zero gravity. The study highlights potential applications of the lattice Boltzmann method for rapid and cost-effective evaluation of alternative plant growth media designs under variable gravity.

NASA Center JSC↗

Wall-Modeled Lattice Boltzmann and Navier-Stokes Approaches for Separated Flows

Lattice Boltzmann (LB) and hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) methods within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are applied to NASA's Revolutionary Computational Aerosciences (RCA) standard test cases for separated flows. A detailed comparison between the performance and accuracy of the two emerging numerical methodologies for turbulence resolving simulations, i.e. the LB and hybrid RANS/LES methods will be presented. This contribution addresses the RCA technical challenge to identify and down-select critical turbulence, transition, and numerical method technologies for 40% reduction in predictive error for standard turbulence separated flow test cases. Results for the 2D NASA wall-mounted hump and the axisymmetric transonic bump including time-averaged pressure coefficient, skin friction, and velocity pro les, as well as resolved and modeled Reynolds stresses for both numerical approaches will be presented and differences between LB and hybrid RANS/LES will be discussed.

Lattice Boltzmann↗

Theory of the lattice Boltzmann Method: Dispersion, Dissipation, Isotropy, Galilean Invariance, and Stability

The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyper-viscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters which optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some 2D LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long wave-length limit (wave vector k = 0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.

Lallemand, Pierre↗

HLPW-4/GMGW-3: Wall-Modeled LES and Lattice-Boltzmann Technology Focus Group Workshop Summary

A summary of the nine submissions to the Wall-Modeled LES and Lattice-Boltzmann(WMLESLB) Technical Focus Group (TFG) at the 4th High lift Prediction Workshop is provided. The focus of this TFG was to assess the current capabilities of WMLES and Lattice-Boltzmann methods on a complex high-lift configuration across a wide range of angles of attack. Preliminary analysis of the submitted data suggests that>250𝑀spatial degrees of freedom are needed to accurately predict pitching moments at high angles of attack due to large pressure gradients present on the outboard slat and main element for𝛼 >17◦(corrected for free-air).While some scatter is reported in pitching moments at the low-angles of attack (𝛼 <11◦) - likely caused by differences in flap separation possibly due to low Reynolds number effects - excellent agreement is observed between the submissions near the𝐶𝐿,max state. Objective superiority over steady state RANS can be seen in terms of lack of excess outboard separation; a majority of the good quality WMLES and LB submissions predict wedge-shaped separation patterns consistent with the experimental oil flow. Differences in the onset of stall mechanism in the free-air configuration for𝛼 >20◦is reported with two distinct topologies observed. Topology A is characterized by the onset of corner-flow separation which progressively grows to produce a pitch break in free-air with an angle of attack increase from𝛼=20.55◦to𝛼=21.47◦+𝜖where𝜖is a small perturbation (varying between submissions and likely to be within±0.3◦). Topology Bis characterized by boundary layer weakness emanating from the inboard side of the wing-pylon juncture substantially larger than any weakness in the wing-body juncture (due to corner-flow);submissions in this category do not show any tendency for occurrence of a pitch break for the free-air configuration within the vicinity of𝛼=21.47◦. The in-tunnel simulations submitted by3 participants using different discretizations, grids and closure models show excellent agreement with the experiment in terms of a) integrated loads, b) surface flow-topology, and c) mechanism for the onset of inboard stall. Further evidence is provided to demonstrate both qualitative and quantitative superiority of all 3 WMLES submissions over a single steady state RANS submission to the workshop.

TTT↗

Towards Full Aircraft Airframe Noise Prediction: Lattice Boltzmann Simulations

Computational results for an 18%-scale, semi-span Gulfstream aircraft model are presented. Exa Corporation's lattice Boltzmann PowerFLOW(trademark) solver was used to perform time-dependent simulations of the flow field associated with this high-fidelity aircraft model. The simulations were obtained for free-air at a Mach number of 0.2 with the flap deflected at 39 deg (landing configuration). We focused on accurately predicting the prominent noise sources at the flap tips and main landing gear for the two baseline configurations, namely, landing flap setting without and with gear deployed. Capitalizing on the inherently transient nature of the lattice Boltzmann formulation, the complex time-dependent flow features associated with the flap were resolved very accurately and efficiently. To properly simulate the noise sources over a broad frequency range, the tailored grid was very dense near the flap inboard and outboard tips. Extensive comparison of the computed time-averaged and unsteady surface pressures with wind tunnel measurements showed excellent agreement for the global aerodynamic characteristics and the local flow field at the flap inboard and outboard tips and the main landing gear. In particular, the computed fluctuating surface pressure field for the flap agreed well with the measurements in both amplitude and frequency content, indicating that the prominent airframe noise sources at the tips were captured successfully. Gear-flap interaction effects were remarkably well predicted and were shown to affect only the inboard flap tip, altering the steady and unsteady pressure fields in that region. The simulated farfield noise spectra for both baseline configurations, obtained using a Ffowcs-Williams and Hawkings acoustic analogy approach, were shown to be in close agreement with measured values.

Khorrami, Mehdi R.↗

Lattice Boltzmann Method for 3-D Flows with Curved Boundary

In this work, we investigate two issues that are important to computational efficiency and reliability in fluid dynamics applications of the lattice, Boltzmann equation (LBE): (1) Computational stability and accuracy of different lattice Boltzmann models and (2) the treatment of the boundary conditions on curved solid boundaries and their 3-D implementations. Three athermal 3-D LBE models (D3QI5, D3Ql9, and D3Q27) are studied and compared in terms of efficiency, accuracy, and robustness. The boundary treatment recently developed by Filippova and Hanel and Met et al. in 2-D is extended to and implemented for 3-D. The convergence, stability, and computational efficiency of the 3-D LBE models with the boundary treatment for curved boundaries were tested in simulations of four 3-D flows: (1) Fully developed flows in a square duct, (2) flow in a 3-D lid-driven cavity, (3) fully developed flows in a circular pipe, and (4) a uniform flow over a sphere. We found that while the fifteen-velocity 3-D (D3Ql5) model is more prone to numerical instability and the D3Q27 is more computationally intensive, the 63Q19 model provides a balance between computational reliability and efficiency. Through numerical simulations, we demonstrated that the boundary treatment for 3-D arbitrary curved geometry has second-order accuracy and possesses satisfactory stability characteristics.

Mei, Renwei↗

Application of Lattice Boltzmann and Navier-Stokes Methods to NASA's Wall Mounted Hump

Lattice Boltzmann (LB) based Large Eddy Simulation (LES), Reynolds-averaged Navier-Stokes (RANS) as well as hybrid RANS/LES methods within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are applied to NASA's wall-mounted hump. Computational results are compared with experiments performed by Greenblatt et al. A detailed comparison between the accuracy and resolution requirements of the two approaches for turbulence resolving simulations, as well as the suitability of different grid paradigms (body-fitted curvilinear and block structured Cartesian) are presented. This test case is part of NASA's Revolutionary Computational Aerosciences (RCA) sub-project which addresses the technical challenge of predicting flow separation and reattachment accurately. Improvements in predictive accuracy by as much as 90% are demonstrated using LB as well as hybrid RANS/LES approaches compared to state-of-the-art steady state RANS simulations.

Lattice Boltzmann↗

Lattice Boltzmann and Navier-Stokes Cartesian CFD Approaches for Airframe Noise Predictions

Lattice Boltzmann (LB) and compressible Navier-Stokes (NS) equations based computational fluid dynamics (CFD) approaches are compared for simulating airframe noise. Both LB and NS CFD approaches are implemented within the Launch Ascent and Vehicle Aerodynamics (LAVA) framework. Both schemes utilize the same underlying Cartesian structured mesh paradigm with provision for local adaptive grid refinement and sub-cycling in time. We choose a prototypical massively separated, wake-dominated flow ideally suited for Cartesian-grid based approaches in this study - The partially-dressed, cavity-closed nose landing gear (PDCC-NLG) noise problem from AIAA's Benchmark problems for Airframe Noise Computations (BANC) series of workshops. The relative accuracy and computational efficiency of the two approaches are systematically compared. Detailed comments are made on the potential held by LB to significantly reduce time-to-solution for a desired level of accuracy within the context of modeling airframes noise from first principles.

Airframe Noise Prediction↗

Performance Enhancements for the Lattice-Boltzmann Solver in the LAVA Framework

Performance enhancements in NASA's recently developed Lattice Boltzmann solver within the Launch Ascent and Vehicle Aerodynamics (LAVA) framework are presented. Two key algorithmic developments are highlighted. A coarse-fine interface treatment that discretely conserves mass and momentum has been implemented and successfully verified and validated. Code optimizations targeting improved serial and parallel performance were presented. For a simple turbulent Taylor-Green Vortex problem, we were able to demonstrate a 2.3 times speedup over the baseline code for a single Skylake-SP node containing 40 physical cores, and a 2.14 times speedup for 64 nodes containing 2560 physical cores. In addition, we were able to show that the optimizations enabled us to scale the code almost perfectly to 20480 physical cores where, including ghost cells, the problem size was 10 billion cells.

LAVA↗