In this work, we use a boundary layer resolved model to improve a Lagrangian point particle model to simulate the combustion of single iron particles. By resolving the full boundary layer, mass and heat transfer are accurately modeled, including Stefan flow. Therefore, the model is suitable to improve point particle models. This work focuses on the first stage of iron combustion, which lasts up to the maximum temperature. Temperature- and composition-dependent properties are used and phase transitions from solid to liquid and liquid to gas are taken into account. The Nusselt and Sherwood correlations are investigated in conditions typical for iron particle combustion. It is found that the $1/2$-film temperature is the best film rule to use to model heat- and mass transfer for iron particle combustion. The boundary layer resolved model is used to validate the point particle models. Then, the model is systematically elaborated by including a temperature-dependent particle density, slip velocity and Stefan flow. The individual and combined effect of these phenomena on the burn duration are investigated. Including all these effects decreases the time to maximum temperature by around $25\mathrm{\%}$. Furthermore, it is shown that if one neglects physical phenomena like slip and Stefan flow, but uses the $1/3$-film rule instead of the $1/2$-film rule, errors cancel and still reasonable agreement is obtained with experiments.

As iron powder nowadays attracts research attention as a carbon-free, circular energy carrier, molecular dynamics (MD) simulations can be used to better understand the mechanisms of liquid iron oxidation at elevated temperatures. However, prudence must be practiced in the selection of a reactive force field. This work investigates the influence of currently available reactive force fields (ReaxFFs) on a number of properties of the liquid iron-oxygen (Fe-O) system derived (or resulting) from MD simulations. Liquid Fe-O systems are considered over a range of oxidation degrees Z_O, which represents the molar ratio of O/(O + Fe), with 0 < Z_O < 0.6 and at a constant temperature of 2000 K, which is representative of the combustion temperature of micrometric iron particles burning in air. The investigated properties include the minimum energy path, system structure, (im)miscibility, transport properties, and the mass and thermal accommodation coefficients. The properties are compared to experimental values and thermodynamic calculation results if available. Results show that there are significant differences in the properties obtained with MD using the various ReaxFF parameter sets. Based on the available experimental data and equilibrium calculation results, an improved ReaxFF is required to better capture the properties of a liquid Fe-O system.