[en] Natural convection within simulated core debris beds has been experimentally and theoretically investigated. The effect of heating method on bed behavior has been found to be important. For directly-heated beds, variations of the downward and upward power fraction and Nusselt numbers with bed loading, power density, particle size, overlying fluid layer height and top surface boundary condition have been determined. Generalized correlations for the upward and downward Nusselt numbers as functions of the internal Rayleigh number have been obtained. Particle tracing techniques have been used to visualize the flow patterns within the bed and overlying fluid layer. The temperature distributions within the bed and overlying fluid layer have also been measured. The experimental data have been compared with COMMIX-lA predictions. Poor agreement has been obtained for both the integral quantities, i.e. downward and upward power fractions and Nusselt numbers, as well as the steady state velocity and temperature distributions. The code does not correctly predict either the magnitude or even the trend of the data