AbstractA finite element formulation of enclosed turbulent diffusion flames is presented. A primitive variables approach is preferred in the analysis. A mixed interpolation is employed for the velocity and pressure. In the solution of the Navier‐Stokes equations, a segregated formulation is adopted, where the pressure discretization equation is obtained directly from the discretized continuity equation, considering the velocity‐pressure relationships in the discretized momentum equations. The state of turbulence is defined by a κ–ϵ model. Near solid boundaries, a wall function approach is employed. The combustion rates are estimated using the eddy dissipation concept. The expensive direct treatment of the integrodifferential equations of radiation is avoided by employing the moment method, which allows the derivation of an approximate local field equation for the radiation intensity. The proposed finite element model is verified by investigating a technical turbulent diffusion flame of semi‐industrial size, and comparing the results with experiments and finite difference pr