本文旨在通过考虑太阳等离子体含有显著比例中性气体的事实，利用真实的三维（3D）太阳磁对流建模来加深对太阳色球层加热机制的理解。为此，我们在三个不同的分辨率下进行了相同物理体积的太阳模拟，即5.76×5.76×2.3 Mm³（其中1.4 Mm位于光学表面之上），分别为20×20×14、10×10×7和5×5×3.5 km³。在所有三种分辨率下，我们比较了包含/不包含双极扩散的时间序列模拟，双极扩散是由于中性气体导致的主要非理想加热机制。我们还比较了三种不同磁化条件下的模拟：（1）小尺度发电机情况；（2）初始植入的垂直磁场为50 G；（3）初始植入的垂直磁场为200 G，但并非所有情况在所有分辨率下都可用。我们发现随着分辨率的提高，模拟的平均磁化程度增加。高分辨率下的平均磁能通量也增加，这意味着在更高分辨率下模拟盒中的磁能量更多。双极扩散作用于相对较大的尺度，实际上可以通过最高分辨率模拟的网格尺度进行数值解析。我们采用两种方法评估双极扩散的数值解析位置：（i）评估模拟的数值扩散并与物理双极扩散进行比较；（ii）将由双极扩散给出的空间尺度与网格分辨率进行数量级比较。在解析位置处，我们比较了包含/不包含双极扩散的模拟平均温度，并得出结论，在包含双极扩散的情况下，经过1200秒的模拟时间后，等离子体平均温度约高600 K。

The aim of this paper is to improve our understanding of the heating mechanisms of the solar chromosphere via realistic three-dimensional (3D) modeling of solar magneto-convection, considering the fact that solar plasma contains a significant fraction of neutral gas. For that we performed simulations of the same physically volume of the Sun, namely 5.76x5.76x2.3 Mm^3 (with 1.4 Mm being above the optical surface), at three different resolutions: 20x20x14, 10x10x7 and 5x5x3.5 km^3. At all three resolutions we compare the time series of simulations with/without ambipolar diffusion, as the main non-ideal heating mechanism due to neutrals. We also compare simulations with three different magnetizations: (1) case of a small-scale dynamo; (2) an initially implanted vertical magnetic field of 50 G; (3) an initially implanted vertical field of 200 G, though not all of them are available at all resolutions. We obtain that the average magnetization of the simulations increases with improving resolution. So does the average magnetic Poynting flux, meaning that there is more magnetic energy in the simulation box at higher resolutions. Ambipolar diffusion operates at relatively large scales, which can be actually numerically resolved with the grid scale of the highest resolution simulations as the ones reported here. We consider two ways of evaluating where the ambipolar scales are numerically resolved: (i) a method to evaluate the numerical diffusion of the simulations and compare it to the physical ambipolar diffusion; (ii) an order of magnitude comparison of spatial scales given by the ambipolar diffusion to our grid resolution. At the resolved locations we compare the average temperature in the simulations with/without ambipolar diffusion, and conclude that the plasma is on average about 600 K hotter after 1200 sec of simulation time when the ambipolar diffusion is included.