Direct numerical simulation (DNS) and large eddy simulation (LES) were performed on the wall‐bounded flow at Re_τ = 180 using lattice Boltzmann method (LBM) and multiple GPUs (Graphic Processing Units). In the DNS, 8 K20M GPUs were adopted. The maximum number of meshes is 6.7 × 10⁷, which results in the nondimensional mesh size of Δ⁺ = 1.41 for the whole solution domain. It took 24 hours for GPU‐LBM solver to simulate 3 × 10⁶ LBM steps. The aspect ratio of resolution domain was tested to obtain accurate results for DNS. As a result, both the mean velocity and turbulent variables, such as Reynolds stress and velocity fluctuations, perfectly agree with the results of Kim et al. (1987) when the aspect ratios in streamwise and spanwise directions are 8 and 2, respectively. As for the LES, the local grid refinement technique was tested and then used. Using 1.76 × 10⁶ grids and Smagorinsky constant (C_s) = 0.13, good results were obtained. The ability and validity of LBM on simulating turbulent flow were verified.

The present study aims to reveal the near-wall vortices and the effects of flow states in coupled-domain transpiration cooling using a recursive regularized thermal lattice Boltzmann method (RR-TLBM). Large-eddy simulations of turbulent flow and heat transfer have been conducted on high-resolution computational grids using a desktop-level computer with CUDA 11.6. Results indicate that the near-wall flow structures present spatial characteristics along the streamwise direction. The vortex evolution promotes the downstream heat dissipation, even though turbulence impairs the effective cooling area. The spanwise evolvement of vortices strengthens the mixing of coolant and hot gas, and small-scale structures are beneficial for turbulent heat transfer. Moreover, the transition onset occurs earlier at higher Reynolds numbers, and it weakens the downstream cooling. The cooling performance of the derived coolant film is improved as the Reynolds number varies from 5 × 10³ to 3 × 10⁴ with a blowing ratio of F = 10%, whereas the local cooling is impaired at the high Reynolds numbers exceeding 5 × 10⁴. The variation in flow states has little influence on the cooling performance at the Reynolds numbers larger than 3 × 10⁶. On the other hand, our in-house RR-TLBM solver is highly stable and efficient for the simulation of flow and heat transfer with high Reynolds numbers. Simultaneously, a high computational performance of 1127 million lattices updated per second is achieved for our simulation of a coupled-domain turbulent flow and heat transfer, using the desktop-level computer with three Tesla V100 graphics processing units.

The structural anisotropy of open-cell foam leads to the anisotropy of effective thermal conductivity (ETC). To quantitatively analyze the effect of structural anisotropy on the anisotropy of ETC, a new predicting model for the ETC of anisotropic open-cell foam was proposed based on an anisotropy tetrakaidecahedron cell (ATC). Feret diameters in three orthogonal directions obtained by morphological analysis of real foam structures were used to characterize the anisotropy of ATC. To validate our proposed anisotropic model, the ETCs of real foam structures in three orthogonal directions predicted by it were compared with the numerical results, for which the structures of numerical models are reconstructed by X-ray computed tomography (X-CT). Using the present anisotropic model, the influences of the thermal conductivity ratio (TCR) and porosity of the foams on the anisotropic ratios of ETCs are also investigated. Results show that there is good consistency between the ETCs obtained by the anisotropic model and the numerical method. The maximum relative errors between them are 2.84% and 13.57% when TCRs are 10 and 100, respectively. The present anisotropic model can not only predict the ETCs in different orthogonal directions but also quantitatively predict the anisotropy of ETC. The anisotropies of the ETCs decrease with porosity because the proportion of the foam skeleton decreases. However, the anisotropies of ETCs increase with TCR, and there exist asymptotic values in anisotropic ratios of ETCs as TCR approaches infinity and they are equal to the relative Feret diameters in different orthogonal directions.

Stress development is one of the significant factors leading to the failure of thermal barrier coating (TBC) systems. In this work, stress development in the two phase mixed zone named phase transition layer (PTL), which grows between the thermally grown oxide (TGO) and the bond coat (BC), is investigated by using two different homogenization models. A constitutive equation of the PTL based on the Reuss model is proposed to study the stresses in the PTL. The stresses computed with the proposed constitutive equation are compared with those obtained with Voigt model-based equation in detail. The stresses based on the Voigt model are slightly higher than those based on the Reuss model. Finally, a further study is carried out to explore the influence of phase transition proportions on the stress difference caused by homogenization models. Results show that the stress difference becomes more evident with the increase of the PTL thickness ratio in the TGO.

The turbulent flow field induced by a round jet in crossflow is calculated by parallel direct numerical simulation (DNS) on multi-GPU clusters. The DNS is based on the lattice- Boltzmannmethod.With...