Search results

The purpose of this paper is to present a fast and bare bones implementation of a numerical method for quickly simulating turbulent thermal flows on GPUs. The work also validates earlier research showing that the lattice Boltzmann method (LBM) method is suitable for complex thermal flows.

A dual lattice hydrodynamic (D3Q27) thermal (D3Q7) multiple-relaxation time LBM model capable of thermal DNS calculations is implemented in CUDA.

The model has the same computational performance compared to earlier publications of similar LBM solvers. The solver is validated against three benchmark cases for turbulent thermal flow with available data and is shown to be in excellent agreement.

The combination of a D3Q27 and D3Q7 stencil for a multiple relaxation time -LBM has, to the authors’ knowledge, not been used for simulations of thermal flows. The code is made available in a public repository under a free license.

The purpose of this paper is the examination of fluid flow around NACA0012 airfoil, with the aim of the numerical validation between the experimental results in the wind tunnel and the Lattice Boltzmann method (LBM) analysis, for the medium Reynolds number (Re = 191,000). The LBM–large Eddy simulation (LES) method described in this paper opens up opportunities for faster computational fluid dynamics (CFD) analysis, because of the LBM scalability on high performance computing architectures, more specifically general purpose graphics processing units (GPGPUs), pertaining at the same time the high resolution LES approach.

Process starts with data collection in open-circuit wind tunnel experiment. Furthermore, the pressure coefficient, as a comparative variable, has been used with varying angle of attack (2°, 4°, 6° and 8°) for both experiment and LBM analysis. To numerically reproduce the experimental results, the LBM coupled with the LES turbulence model, the generalized wall function (GWF) and the cumulant collision operator with D3Q27 velocity set has been used. Also, a mesh independence study has been provided to ensure result congruence.

The proposed LBM methodology is capable of highly accurate predictions when compared with experimental data. Besides, the special significance of this work is the possibility of experimental and CFD comparison for the same domain dimensions.

Considering the quality of results, root-mean-square error (RMSE) shows good correlations both for airfoil’s upper and lower surface. More precisely, maximal RMSE for the upper surface is 0.105, whereas 0.089 for the lower surface, regarding all angles of attack.

This study aims at investigating the heat transfer characteristics of a nonsquare enclosure when hydrodynamic resistance is altered discontinuously along its inner surface. Particularly, it focuses on investigating how several essential factors collaboratively influence the natural convection, including the Rayleigh number (Ra), the aspect ratio (AR), the nanoparticle volume fraction (ϕ) and the locations of changing hydrodynamic resistance.

To achieve these objectives, an L-shaped enclosure of various AR is adopted, while zero local shear resistance is applied and modeled by stress-free (SF) patches of four distinct arrangements (corresponding to Cases 1–4). The nanofluid is modeled by Buongiorno’s two-phase model. The effects are explored using an in-house numerical framework based on a hybrid lattice Boltzmann-finite difference method with the total variation minimization scheme.

The results show that when Ra is sufficiently large, i.e. Ra = 10⁵, SF patches can generally enhance the heat transfer performance regardless of other factors. However, the ways of achieving those enhancements are different, which mainly depend on the arrangement of the SF patches and AR but are nearly independent of ϕ. The maximum improvement of heat transfer can be achieved in Case 3 with AR = 0.6, Ra = 10⁵ and ϕ = 0.04, where the averaged Nusselt number is enhanced by 8.89%.

This study presents a new scenario where the SF patches of various arrangements are applied to enhance the nanofluid natural convection of a nonsquared enclosure, and it reveals how the improvement is achieved and cooperatively affected by several important factors.

Present study aims to extend the lattice Boltzmann method (LBM) to simulate radiation in geometries with curved boundaries, as the first step to simulate radiation in complex porous media. In recent years, researchers have increasingly explored the use of porous media to improve the heat transfer processes. The lattice Boltzmann method (LBM) is one of the most effective techniques for simulating heat transfer in such media. However, the application of the LBM to study radiation in complex geometries that contain curved boundaries, as found in many porous media, has been limited.

The numerical evaluation of the effect of the radiation-conduction parameter and extinction coefficient on temperature and incident radiation distributions demonstrates that the proposed LBM algorithm provides highly accurate results across all cases, compared to those found in the literature or those obtained using the finite volume method (FVM) with the discrete ordinates method (DOM) for radiative information.

For the case with a conduction-radiation parameter equal to 0.01, the maximum relative error is 1.9% in predicting temperature along vertical central line. The accuracy improves with an increase in the conduction-radiation parameter. Furthermore, the comparison between computational performances of two approaches reveals that the LBM-LBM approach performs significantly faster than the FVM-DOM solver.

The difficulty of radiative modeling in combined problems involving irregular boundaries has led to alternative approaches that generally increase the computational expense to obtain necessary radiative details. To address the limitations of existing methods, this study presents a new approach involving a coupled lattice Boltzmann and first-order blocked-off technique to efficiently model conductive-radiative heat transfer in complex geometries with participating media. This algorithm has been developed using the parallel lattice Boltzmann solver.

The purpose of this paper is to develop a coupled lattice Boltzmann model for the simulation of the freezing process in unsaturated porous media.

In the developed model, the porous structure with complexity and disorder was generated by using a stochastic growth method, and then the Shan-Chen multiphase model and enthalpy-based phase change model were coupled by introducing a freezing interface force to describe the variation of phase interface. The pore size of porous media in freezing process was considered as an influential factor to phase transition temperature, and the variation of the interfacial force formed with phase change on the interface was described.

The larger porosity (0.2 and 0.8) will enlarge the unfrozen area from 42 mm to 70 mm, and the rest space of porous medium was occupied by the solid particles. The larger specific surface area (0.168 and 0.315) has a more fluctuated volume fraction distribution.

The concept of interfacial force was first introduced in the solid–liquid phase transition to describe the freezing process of frozen soil, enabling the formulation of a distribution equation based on enthalpy to depict the changes in the water film. The increased interfacial force serves to diminish ice formation and effectively absorb air during the freezing process. A greater surface area enhances the ability to counteract liquid migration.

The paper aims to expand the scope of application of the lattice Boltzmann method (LBM), especially in the field of aircraft engineering. The traditional LBM is usually applied to incompressible flows at a low Reynolds number, which is not sufficient to satisfy the needs of aircraft engineering. Devoted to tackling the defect, the paper proposes a developed LBM combining the subgrid model and the multiple relaxation time (MRT) approach. A multilayer adaptive Cartesian grid method to improve the computing efficiency of the traditional LBM is also employed.

The subgrid model and the multilayer adaptive Cartesian grid are introduced into MRT-LBM for simulations of incompressible flows at a high Reynolds number. Validated by several typical flow simulations, the numerical methods in this paper can efficiently study the flows under high Reynolds numbers.

Some numerical simulations for the lid-driven flow of cavity, flow around iced GLC305, LB606b and ONERA-M6 are completed. The paper presents the investigation results, indicating that the methods are accurate and effective for the separated flow after icing.

LBM is developed with the addition of the subgrid model and the MRT method. A numerical strategy is proposed using a multilayer adaptive Cartesian grid method and its treatment of boundary conditions. The paper refers to innovative algorithm developments and applications to the aircraft engineering, especially for iced wing simulations with flow separations.

This paper aims to establish a lattice Boltzmann (LB) method for solid-liquid phase transition (SLPT) from the pore scale to the representative elementary volume (REV) scale. By applying this method, detailed information about heat transfer and phase change processes within the pores can be obtained, while also enabling the calculation of larger-scale SLPT problems, such as shell-and-tube phase change heat storage systems.

Three-dimensional (3D) pore-scale enthalpy-based LB model is developed. The computational input parameters at the REV scale are derived from calculations at the pore scale, ensuring consistency between the two scales. The approaches to reconstruct the 3D porous structure and determine the REV of metal foam were discussed. The implementation of conjugate heat transfer between the solid matrix and the solid−liquid phase change material (SLPCM) for the proposed model is developed. A simple REV-scale LB model under the local thermal nonequilibrium condition is presented. The method of bridging the gap between the pore-scale and REV-scale enthalpy-based LB models by the REV is given.

This coupled method facilitates detailed simulations of flow, heat transfer and phase change within pores. The approach holds promise for multiscale calculations in latent heat storage devices with porous structures. The SLPT of the heat sinks for electronic device thermal control was simulated as a case, demonstrating the efficiency of the present models in designing and optimizing SLPT devices.

A coupled pore-scale and REV-scale LB method as a numerical tool for investigating phase change in porous materials was developed. This innovative approach allows for the capture of details within pores while addressing computations over a large domain. The LB method for simulating SLPT from the pore scale to the REV scale was given. The proposed method addresses the conjugate heat transfer between the SLPCM and the solid matrix in the enthalpy-based LB model.

The purpose of this study is the hydrothermal analysis of the natural convection phenomenon within the heat exchanger containing nanofluids using the lattice Boltzmann method (LBM).

The thermal conductivity as well as dynamic viscosity of the CuO–water nanofluid is estimated using the Koo-Kleinstreuer-Li model. The LBM has been used with unique modifications to make it flexible with the curved boundaries. The local as well as total entropy generation assessment, local Nusselt variation, as well as heatline visualization are used.

The solid volume percentage of the CuO–water nanofluid, a range of Rayleigh numbers (Ra) and thermal settings of internal operational fins and bodies are all factors that have been thoroughly researched to determine their effects on entropy production, heat transfer efficiency and nanofluid flow.

The originality of this work is using a novel numerical method (i.e. curved boundary LBM) as well as the local/volumetric second law analysis for the application of heat exchanger hydrothermal analysis.

Accurately evaluating fluid flow behaviors and determining permeability for deforming porous media is time-consuming and remains challenging. This paper aims to propose a novel machine-learning method for the rapid estimation of permeability of porous media at different deformation stages constrained by hydro-mechanical coupling analysis.

A convolutional neural network (CNN) is proposed in this paper, which is guided by the results of finite element coupling analysis of equilibrium equation for mechanical deformation and Boltzmann equation for fluid dynamics during the hydro-mechanical coupling process [denoted as Finite element lattice Boltzmann model (FELBM) in this paper]. The FELBM ensures the Lattice Boltzmann analysis of coupled fluid flow with an unstructured mesh, which varies with the corresponding nodal displacement resulting from mechanical deformation. It provides reliable label data for permeability estimation at different stages using CNN.

The proposed CNN can rapidly and accurately estimate the permeability of deformable porous media, significantly reducing processing time. The application studies demonstrate high accuracy in predicting the permeability of deformable porous media for both the test and validation sets. The corresponding correlation coefficients (R²) is 0.93 for the validation set, and the R² for the test set A and test set B are 0.93 and 0.94, respectively.

This study proposes an innovative approach with the CNN to rapidly estimate permeability in porous media under dynamic deformations, guided by FELBM coupling analysis. The fast and accurate performance of CNN underscores its promising potential for future applications.

The purpose of this paper is mainly to know: (1) the sound absorption coefficient of porous composite structures constituted by a new kind of lightweight ceramic foam and perforated plate; (2) the availability of an equivalent porous material model, recently proposed by the present author, to these composite structures in sound absorption.

A kind of lightweight ceramic foam with bulk density of 0.38–0.56 g·cm-3 was produced by means of molding, drying and sintering. The effect of stainless steel perforated plate on sound absorption performance of the ceramic foam was investigated by means of JTZB absorption tester.

The results indicate that the sound absorption performance could be obviously changed by adding the stainless steel perforated plate in front of the porous samples and the air gap in back of the porous samples. Adding the perforated plate to the porous sample with a relatively large pore size, the sound absorption performance could be evidently improved for the composite structure. When the air gap is added to the composite structure, the first absorption peak shifts to the lower frequency, and the sound absorption coefficient could increase in the low frequency range.

Based on the equivalent porous material model and the “perforated plate with air gap” model, the sound absorption performance of the composite structures can be simulated conveniently to a great extent by using Johnson-Champoux-Allard model.