a b s t r a c tIn this work the simulation of velocity and temperature distributions inside a refrigerated vehicle is evaluated. For this purpose a 3D double distribution lattice Boltzmann method (LBM) with the Bhatnagar-Gross-Krook (BGK) collision operator is coupled by the buoyancy force calculated with the Boussinesq approximation. This LBM is extended by a Smagorinsky subgrid method, which numerically stabilizes the BGK scheme for low resolutions and high Reynolds and Rayleigh numbers. Besides validation against the two benchmark cases porous plate and natural convection in a square cavity evaluated at resolutions of y + ≈ 2 for Ra numbers between 10 3 and 10 10 , the method and its implementation are tested via comparison with experimental data for a refrigerated vehicle at Re ≈ 53 000.The aim of the investigation is to provide a deeper understanding of the refrigerated vehicle's insulation processes including its thermal performance under turbulent flow conditions. Therefore, we extend this method by the half lattice division scheme for conjugate heat transfer to simulate in the geometry of a refrigerated vehicle including its insulation walls. This newly developed method combination enables us to accurately predict velocity and temperature distributions inside the cooled loading area, while spatially resolving the heat flux through the insulation walls. We simulate the time dependent heating process of the open door test and validate against measurements at four characteristic velocity and 13 temperature positions in the truck. (M. Gaedtke). are complex material requirements which, in addition to the insulation, include static stability, load securing as well as resistance against vibrational stress.According to Smale [2] most simulations and models for the representation of a velocity and temperature distribution in the field of applied cooling of geometries between 1974 and 2006 have shown low accuracy and large deviations from experimental data. Tabsoba et al. [3] studied the influence of the k-ϵ turbulence model against the Reynolds stress model, where they indicated the k-ϵ model to fail at the prediction of certain flow characteristics in ventilated enclosures. Ambaw [4] summarized studies on the cooling of harvested food, indicating some progress being made between 2006 and 2013 with the accuracy of the models. With simulations mainly based on Reynolds averaging turbulence models as the k-ϵ or shear stress transport (SST) model [5][6][7][8], he also came to the conclusion, that a clear increase in the accuracy of turbulent 3D simulations within complex cooled geometries compared with experimental data has still to come.According to James [9] first simulations of the open door tests have been carried out by Tso et al. [10]. Tso's simulated temperatures deviate from experimental recordings by up to 4 K at a temperature difference of about 26 K. Son [11] showed an approximation to the validation data up to 1.13 K in the simulation of the velocity and temperature distribution in the interior of a filled...