In recent years, the use of Reynolds-Averged Navier-Stokes (RANS) based CFD within industry has led to a sea change in the design and development approaches used across all aspects of engineering, allowing more efficient design, reducing development costs, and bringing products to market more rapidly and economically than was previously possible. Because of computer run time constraints, RANS approaches will remain the principal method for representing the effects of turbulence for engineering application for the foreseeable future. However, the turbulence models embodied within such approaches are still dependent on the flow configuration being examined, requiring the use of experimental data in their specification, and case-by-case adjustment to accurately predict particular types of flow. Turbulent shear flows with free boundaries display an intermittent character, such that the fluctuations of velocity rapidly change from rotational to irrotational, and vice versa. Intermittency is important in many practical flows, for example, in boundary layer transition which remains a major challenge because of its relevance to the design of modern compressors and turbines, and airfoils. Intermittency is also important in the safe and efficient operation of many combustion devices, which are dependent on ignition occurring in flows with significant intermittency, as well as in hazard and risk assessments concerning the potential ignition of flammable releases on chemical and process plant. In all these applications, the enhanced performance of these devices, and improved safety, relies on a detailed understanding of turbulent flow and ignition processes, and hence intermittency. The majority of engineering turbulence models in use today were derived for fully developed flows, and hence cannot be expected to accurately predict free shear flows where the outer regions are contaminated with irrotational flow. The development of turbulence models that accommodate intermittency effects is therefore also of fundamental importance. Additionally, the inclusion of such effects helps to generalize existing turbulence models, and provide more accurate predictions of the velocity and scalar fields of interest to the applications noted above. The necessity for such work has also been identified by both the EROFTAC Special Interest Group on Transition Modelling, and the Isaac Newton Institute for Mathematical Sciences Programme on Turbulence, particularly in relation to the need for greater universality in engineering models of turbulence, and the improved modelling of the influence of intermittency, particularly on scalar fields. This proposal concerns the use of an hierarchical approach to the development of more accurate engineering models of intermittent turbulent flows through the co-ordinated application of direct numerical simulation and large eddy simulation to inform the formulation and validation of second-moment turbulence closures that are more fundamentally based that existing approaches. The work proposed will also improve our overall understanding of such flows, the universality of engineering models of turbulence, and will validate the models developed for application to a wide range of practical flows.