An implicit method for turbulent boundary layers simulation

Here, we describe a numerical method which combines a mixed finite volume/finite element formulation on unstructured triangulations and a linearized implicit method for time advancing. A particular attention has been paid to the source term implicitation. This method is applied to three typical low-Reynolds number $k-\varepsilon$ models, namely those of Lam-Bremhorst, Nagano-Tagawa and Speziale-Abid-Anderson. The selected test cases consist of a subsonic flow in a duct and a supersonic flow over a flat plate