The turbulent flow characteristics in a coaxial injector were investigated by the nonlinear $k-{varepsilon}-f_{mu}$ model of Park et al.[1] and large eddy simulation (LES). In order to analyze the geometric effects on the scalar mixing for nonreacting variable-density flows, several recessed lengths and momentum flux ratios are selected at a constant Reynolds number. The nonlinear $k-{varepsilon}-f_{mu}$?? model proposed the meaningful characteristics for various momentum flux ratios and recess lengths. The LES results showed the changes of small-scale structures by the recess. When the inner jet was recessed, the development of turbulent kinetic energy became faster than that of non-recessed case. Also, the mixing characteristics were mainly influenced by the variation of shear rates, but the local mixing was changed by the adoption of recess.