Modeling of the nitric oxide formation in the turbulent nonpremixed methane and syngas (CO/H-2/N-2) diffusion flame is studied using an implicit time stepping and multigrid technique. The chemical kinetic model for both methane-air in a Sudden-expansion combustor and syngas-air combustion in a laboratory combustor as well as in a gas turbine combustor is assumed to have 49 species and 229 finite-rate, reversible reaction steps. The standard k-epsilon turbulence model and the algebraic correlation closure model are applied to close the time-averaged Navier-Stokes and species equations (Liao at al., 1995) respectively. The computation requires about 250 time steps to reduce the residual by 3 orders of magnitude for the 3-D turbulent methane-air diffusion dame case on a 34 x 18 x 18 grid, which shows convergence rate is much faster than conventional iterative methods. Computational results with detailed chemistry are exhibited and some of them are compared with experimental data. Qualitative agreement between the computational results and experimental results is observed. A three-dimensional calculation for a gas turbine combustor shows the potential of this method for modeling the pollution formation in practical flows.