Example of a temperature profile (circles) and the corresponding particle number density profile (squares). The symbols report the simulated values and the lines correspond respectively to Eqs. (4) and (5).
(a) Temperature profile with h = 0.5 and walls with virtual particles. Open circles are simulation results, solid lines are the temperatures imposed at the walls, and dashed lines the actual boundary temperatures T(0) and T(L z ). Otherwise stated the employed parameters are h = 0.1, α = 120, ρ = 5, a cubic size L = 40, T c = 0.9, and T h = 1.1. (b) Dimensionless temperature jump T J in Eq. (8) as a function of α (circles, down-axis), and as a function of h (squares, up-axis). (c) Dependence of T J with the inverse of the system length L z in the temperature gradient direction, with α = 30 and L x = 20 = L y .
T J for walls with thermostats as a function of α, and h. Symbols and parameters similar to Fig. 2. The inset shows the detail of a temperature profile close to the wall for with α = 30, h = 0.1, ρ = 5, and L z = 40.
Illustration of the periodic simulation box in the presence of a temperature gradient.
(a) Temperature profile obtained from the velocity exchange algorithm with h = 0.1, α = 120, and ρ = 10. Symbols correspond to the measured temperatures, dashed-line is the estimated temperature profile from Eq. (10). (b) Velocity squared distribution in Eq. (9) for the temperatures T c = 0.9 and T h = 1.1, typically used for the cold and hot baths.
Thermal diffusivity for two values of ρ. (a) k T as a function of α for h = 0.1. (b) k T as a function of h for α = 120. The insets are a zoom-in for large values of α and small h, respectively. Lines correspond to the analytical approach in Eq. (11) and symbols to simulation results. Continuous lines correspond to ρ = 5 and dashes lines to ρ = 20. Triangles refer to simulations with walls and thermostats, circles to walls with virtual particles, and squares to the velocity exchange algorithm.
Relative deviation of the simulated thermal diffusivity k T, sim with respect to the analytical approach k T, an , Δk T /k T ≡ (k T, sim − k T, an )/k T, an a) as a function of α (b) as a function of h. Simulation values are those presented in Fig. 6 for walls and thermostats. Open symbols and dashed lines employ both analytical contributions in Eq. (11), while solid symbols and solid lines take into account the collisional contribution in Eq. (15).
Article metrics loading...
Full text loading...