Search results

The purpose of this paper is to investigate the effects of geometrical parameters, eccentricity and perforated fins on natural convection heat transfer in a finned horizontal annulus using three-dimensional lattice Boltzmann flux solver.

Three-dimensional lattice Boltzmann flux solver is used in the present study for simulating conjugate heat transfer within an annulus. D3Q15 and D3Q7 models are used to solve the fluid flow and temperature field, respectively. The finite volume method is used to discretize mass, momentum and energy equations. The Chapman–Enskog expansion analysis is used to establish the connection between the lattice Boltzmann equation local solution and macroscopic fluxes. To improve the accuracy of the lattice Boltzmann method for curved boundaries, lattice Boltzmann equation local solution at each cell interface is considered to be independent of each other.

It is found that the maximum heat transfer rate occurs at low fin spacing especially by increasing the fin height and decreasing the internal-cylindrical distance. The effect of inner cylinder eccentricity is not much considerable (up to 5.2% enhancement) while the impact of fin eccentricity is more remarkable. Negative fin eccentricity further enhances the heat transfer rate compared to a positive fin eccentricity and the maximum heat transfer enhancement of 91.7% is obtained. The influence of using perforated fins is more considerable at low fin spacing although some heat transfer enhancements are observed at higher fin spacing.

The originality of this paper is to study three-dimensional natural convection in a finned-horizontal annulus using three-dimensional lattice Boltzmann flux solver, as well as to apply symmetry and periodic boundary conditions and to analyze the effect of eccentric annular fins (for the first time for air) and perforated annular fins (for the first time so far) on the heat transfer rate.

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.

Lattice Boltzmann method (LBM) has made great success in computational fluid dynamics, and this paper aims to establish an efficient simulation model for the polymer injection molding process using the LBM. The study aims to validate the capacity of the model for accurately predicting the injection molding process, to demonstrate the superior numerical efficiency in comparison with the current model based on the finite volume method (FVM).

The study adopts the stable multi-relaxation-time scheme of LBM to model the non-Newtonian polymer flow during the filling process. The volume of fluid method is naturally integrated to track the movement of the melt front. Additionally, a novel fractional-step thermal LBM is used to solve the convection-diffusion equation of the temperature field evolution, which is of high Peclet number. Through various simulation cases, the accuracy and stability of the present model are validated, and the higher numerical efficiency verified in comparison with the current FVM-based model.

The paper provides an efficient alternative to the current models in the simulation of polymer injection molding. Through the test cases, the model presented in this paper accurately predicts the filling process and successfully reproduces several characteristic phenomena of injection molding. Moreover, compared with the popular FVM-based models, the present model shows superior numerical efficiency, more fit for the future trend of parallel computing.

Limited by the authors’ hardware resources, the programs of the present model and the FVM-based model are run on parallel up to 12 threads, which is adequate for most simulations of polymer injection molding. Through the tests, the present model has demonstrated the better numerical efficiency, and it is recommended for the researcher to investigate the parallel performance on even larger-scale parallel computing, with more threads.

To the authors’ knowledge, it is for the first time that the lattice Boltzmann method is applied in the simulation of injection molding, and the proposed model does obviously better in numerical efficiency than the current popular FVM-based models.

The purpose of the current paper is to develop a numerical methodology, based on the immersed boundary-lattice Boltzmann computational framework, for the Neumann and Dirichlet boundary conditions in problems involving natural and forced convection heat transfer.

The direct forcing immersed boundary method is extended to study the heat transfer by incompressible flow within the thermal lattice Boltzmann method (LBM) computational framework. The direct forcing and heating immersed boundary-LBM introduces a heat source term to the thermal LBM to account for the heat transfer occurring at the immersed boundary. New numerical treatments for the Neumann type of boundary condition and for the calculation of the local Nusselt number are developed. The developed methodologies have been applied to flows around immersed bodies with natural and forced convection, including steady as well as unsteady flows.

Numerical experiments involving immersed bodies in natural and forced convection have been performed in order to assess the validity of the direct heating IB-LBM. The flow cases studied also include steady and transient flow phenomena. Flow velocity field and isotherms have been used for qualitative comparisons with existing, published results. The surface averaged Nusselt number, Strouhal number, and lift coefficient (for the unsteady flow cases) have been used for quantitative comparison with published results. The results show that there are satisfactory agreements, qualitatively and quantitatively, between the results obtained by using the present method and those previously published.

Limited application of immersed boundary to thermal flows within the LBM has been studied by researchers; the few past studies were limited to Dirichlet boundary conditions and/or using of feedback forcing and heating approaches. In the current paper, the direct forcing and heating approach was used which helps to eliminate the arbitrary constants used in the feedback approaches. The developed new numerical treatments for the Neumann type of boundary condition and for the calculation of the local Nusselt number eliminate the need to determine surface normal and temperature gradient in the normal direction for heat transfer calculation, which is particularly beneficial in cases with deforming or changing boundaries.

The Lattice Boltzmann (LB) method offers an alternative to conventional computational fluid dynamics (CFD) methods. However, its practical use for complex turbulent flows of engineering interest is still at an early stage. This paper aims to outline an LB wall-modeled large-eddy simulation (WMLES) solver.

The solver is dedicated to complex high-Reynolds flows in the context of WMLES. It relies on an improved LB scheme and can handle complex geometries on multi-resolution block structured grids.

Dynamic and acoustic characteristics of a turbulent airflow past a rod-airfoil tandem are examined to test the capabilities of this solver. Detailed direct comparisons are made with both experimental and numerical reference data.

This study allows assessing the potential of an LB approach for industrial CFD applications.

The paper targets on providing new experimental data for validation of the well-established mathematical models within the framework of the lattice Boltzmann method (LBM), which are applied to problems of casting processes in complex mould cavities.

An experimental campaign aiming at the free-surface flow within a system of narrow channels is designed and executed under well-controlled laboratory conditions. An in-house lattice Boltzmann solver is implemented. Its algorithm is described in detail and its performance is tested thoroughly using both the newly recorded experimental data and well-known analytical benchmark tests.

The benchmark tests prove the ability of the implemented algorithm to provide a reliable solution when the surface tension effects become dominant. The convergence of the implemented method is assessed. The two new experimentally studied problems are resolved well by simulations using a coarse computational grid.

A detailed set of original experimental data for validation of computational schemes for simulations of free-surface gravity-driven flow within a system of narrow channels is presented.

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 main purpose of this paper is to develop a novel thermal lattice Boltzmann method (LBM) based on finite volume (FV) formulation. Validation of the suggested formulation is performed by simulating plane Poiseuille, backward-facing step and flow over circular cylinder.

For this purpose, a cell-centered scheme is used to discretize the convection operator and the double distribution function model is applied to describe the temperature field. To enhance stability, weighting factors are defined as flux correctors on a D2Q9 lattice.

The introduction of pressure-temperature-dependent flux-control coefficients in the streaming operator, in conjunction with suitable boundary conditions, is shown to result in enhanced numerical stability of the scheme. In all cases, excellent agreement with the existing literature is found and shows that the presented method is a promising scheme in simulating thermo-hydrodynamic phenomena.

A stable and accurate FV formulation of the thermal DDF-LBM is presented.

This aim of this work is to investigate different modelling approaches for air-cooled data centres. The study employs three computational methods, which are based on finite element, finite volume and lattice Boltzmann methods and which are respectively implemented via commercial Multiphysics software, open-source computational fluid dynamics code and graphical processing unit-based code developed by the authors. The results focus on comparison of the three methods, all of which include models for turbulence, when applied to two rows of datacom racks with cool air supplied via an underfloor plenum.

This paper studies thermal airflows in a data centre by applying different numerical simulation techniques that are able to analyse the thermal airflow distribution for a simplified layout of datacom racks in the presence of a computer room air conditioner.

Good quantitative agreement between the three methods is seen in terms of the inlet temperatures to the datacom equipment. The computational methods are contrasted in terms of application to thermal management of data centres.

The work demonstrates how the different simulation techniques applied to thermal management of airflow in a data centre can provide valuable design and operational understanding. Basing the analysis on three very different computational approaches is new and would offer an informed understanding of their potential for a class of problems.

This study aims to feature the application of the artificial compressibility method (ACM) for the numerical prediction of two-dimensional (2D) axisymmetric swirling flows.

The respective academic numerical solver, named IGal2D, is based on the axisymmetric Reynolds-averaged Navier–Stokes (RANS) equations, arranged in a pseudo-Cartesian form, enhanced by the addition of the circumferential momentum equation. Discretization of spatial derivative terms within the governing equations is performed via unstructured 2D grid layouts, with a node-centered finite-volume scheme. For the evaluation of inviscid fluxes, the upwind Roe’s approximate Riemann solver is applied, coupled with a higher-order accurate spatial reconstruction, whereas an element-based approach is used for the calculation of gradients required for the viscous ones. Time integration is succeeded through a second-order accurate four-stage Runge-Kutta method, adopting additionally a local time-stepping technique. Further acceleration, in terms of computational time, is achieved by using an agglomeration multigrid scheme, incorporating the full approximation scheme in a V-cycle process, within an efficient edge-based data structure.

A detailed validation of the proposed numerical methodology is performed by encountering both inviscid and viscous (laminar and turbulent) swirling flows with axial symmetry. IGal2D is compared against the commercial software ANSYS fluent – by using appropriate metrics and characteristic flow quantities – but also against experimental measurements, confirming the proposed methodology’s potential to predict such flows in terms of accuracy.

This study provides a robust methodology for the accurate prediction of swirling flows by combining the axisymmetric RANS equations with ACM. In addition, a detailed description of the convective flux Jacobian is provided, filling a respective gap in research literature.