The limitation of nitrogen oxides emissions is nowadays a challenge in several engineering fields. Recent European regulations have reduced the maximum NOx emissions and therefore forced the glass production sector to develop emission reduction strategies. Two different systems have been developed within the framework of the European LIFE project and are currently applied to glass regenerative furnaces: the Waste Gas Recirculation (WGR) and the Hybrid Air Staging (HyAS). The above systems are primary NOx reduction strategies because they both operate to control the combustion evolution. Both WGR and HyAS systems have been conceived with the extensive use of Computational Fluid Dynamics (CFD) models: design strategies for both systems have been developed based on the use of CFD and are currently under use by glass furnace designers. In the present work, the CFD procedures routinely used for the design of the above systems are described. The systems effectiveness, due to the harsh conditions in the industrial installation, can be tested with oxygen concentration measurements inside the regenerators. The oxygen concentration is correlated to the flame evolution and therefore to the nitrogen oxides formation. For the above reason, the models have been validated with experimental data from pilot furnaces using measured values of O2 mole fraction. The CFD procedures are described in the paper together with their application to different configurations.