Smoldering, characterized as a slow, low-temperature, and flameless reaction process, represents one of the most persistent types of combustion phenomena. While oxygen supply is one of the governing mechanisms controlling smoldering, the specific oxygen threshold or smothering limit remains inadequately investigated. Herein, we built a physics-based 1-D computational model integrating heat-and-mass transfer and 5-step heterogeneous chemistry to investigate the oxygen threshold or smothering limit of smoldering propagation in porous pine needle beds subjected to forced internal oxidizer flow. Simulation results revealed that, the required oxidizer flow velocity or oxygen supply rate increased as the oxygen concentration decreased, and the predicted limiting oxygen concentration (LOC) was about 3 %, agreeing well with the experimental observations and theoretical analysis. Moreover, the required airflow velocity was predicted to increase as the fuel density and environmental temperature decreased, or the moisture content increased, and the predicted maximum moisture content capable of supporting smoldering was about 110 %. At the smothering limit, the modeled minimum smoldering temperature and propagation rate were around 300 °C and 0.5 cm/h. This work helps deepen our understanding of the limiting conditions of smoldering combustion, thus improving the mitigation strategy of smoldering fire and the efficiency of applied smoldering systems. Novelty and significance statementsOxygen supply is one of the key mechanisms governing the smoldering propagation. In the literature, scattered studies have examined the limiting oxygen concentration (LOC) of smoldering under different conditions, leading to disparate results even for the same fuels. In our previous work, for the first time, we developed a tubular smoldering reactor capable of precisely controlling the flow of oxidizer, and successfully quantified the exact oxygen supply limit of smoldering combustion. However, no computational model was established specifically for the oxygen threshold of smoldering with forced internal oxygen supply.The novelty of this research is the establishment of the first-ever computational model to investigate the oxygen threshold or smothering limit of smoldering propagation in porous pine needle beds subjected to forced internal oxidizer flow. It is significant because (1) it contributes to a fundamental understanding of smoldering combustion; (2) it investigates the role of fuel properties and environmental conditions in smoldering dynamics and the oxygen supply limits; (3) it assists in optimizing the applied smoldering systems and enriching prevention strategies for smoldering fires.