Search results

It is of great significance to study the influence of subgrade filling on permafrost temperature field in permafrost area for the smooth construction and safe operation of railway.

The paper builds up the model for the hydrothermal coupling calculation of permafrost using finite element software COMSOL to study how permafrost temperature field changes in the short term after subgrade filling, on which basis it proposes the method of calculation for the concave distortion of freezing front in the subgrade-covered area.

The results show that the freezing front below the subgrade center sinks due to the thermal effect of subgrade filling, which will trigger hydrothermal erosion in case of sufficient moisture inflows, leading to the thawing settlement or the cracking of the subgrade, etc. The heat output of soil will be hindered the most in case of July filling, in which case the sinking and the distortion of the freezing front is found to be the most severe, which the recovery of the permafrost temperature field, the slowest, constituting the most unfavorable working condition. The concave distortion of the freezing front in the subgrade area increases with the increase in temperature difference between the filler and ground surface, the subgrade height, the subgrade width and the volumetric thermal capacity of filler, while decreases with the increase of the thermal conductivity of filler. Therefore, the filler chose for engineering project shall be of small volumetric thermal capacity, low initial temperature and high thermal conductivity whenever possible.

The concave distortion of the freezing front under different working conditions at different times after filling can be calculated using the method proposed.

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.