In order to meet future emission norms for small carbureted SI engines, such as those used on motorcycles in India, there is a need to study mixture preparation, specifically the two-phase flow exiting the carburetor and entering the inlet manifold. A fully functional, modular experimental rig is designed and erected for performing both qualitative and quantitative flow visualization. The vibrations of the engine are minimized to reduce their effect on the flow. A special, optically accessible tube of square cross-section is added between the carburetor and the inlet manifold, to enable the visualization of flow at the exit of the carburetor. An electronic circuit to obtain a signal for the engine crank angle and convert it to a standard TTL pulse, for use on standard imaging systems to capture cycle resolved-images is also designed. The flow in the optical section is qualitatively visualized using high and low speed cameras. The resulting images and movies show two modes of fuel transport within the inlet manifold, one of which is in the form of a dense cloud of fine fuel droplets during some part of the intake stroke. The second mode is in the form of a film at all times in the cycle, along the lower surface of the inlet manifold during idling and along vertical walls under loaded conditions. Recirculation is seen on the vertical walls of the manifold during idling and under load. Finally, the thickness of the fuel film in the optical section at the exit of the carburetor is measured, using PLIF. This part of the study also reveals that there is a film on upper surface of the optical section, at all loads and speeds. This film is lesser than the resolution of measurement for low loads, and increases to 0.5 mm in the case of highest load and speed attained at full throttle. In contrast to the loaded conditions, during idling, the film occurs on the lower surface of the manifold and its thickness is highest (1 mm.). The film is also present throughout the cycle during idling and all load-speed conditions, suggesting that the mixture that goes into the engine has a significant part of fuel in liquid form.

With the age of hypersonic flight imminent just beyond the horizon, researchers are working hard at designing work-arounds for all the major problems as well as the minor quirks associated with it. One such issue, seemingly innocuous but one that could be potentially deadly, is the problem of interference heating due to surface protuberances. Although an ideal design of the external surfaces of a high-speed aircraft dictates complete smoothness to reduce drag, this is not always possible in reality. Control surfaces, sheet joints, cable protection pads etc. generate surface discontinuities of varying geometries, in the form of both protrusions as well as cavities. These discontinuities are most often small in dimension, comparable to the local boundary layer thickness at that location. Such protuberances always experience high rates of heat transfer, and therefore should be appropriately shielded. However, thermal shielding of the protrusions alone is not a full solution to the problem at hand. The interference caused to the boundary layer by the flow causes the generation of local hot spots in the vicinity of the protuberances, which should be properly mapped and adequately addressed. The work presented in this thesis aims at locating and measuring the heat flux values at these hot spots near the protrusions, and possibly formulating empirical correlations to predict the hot spot heat flux for a given set of flow conditions and protrusion geometry. Experimental investigations were conducted on a flat plate model and a cone model, with interchangeable sharp and blunt nose tips, with attached 3D protuberances. Platinum thin-film sensors were placed around the protrusion so that the heat fluxes could be measured in its vicinity and the hottest spot located. These experiments were carried out at five different hypersonic free stream flow conditions generated using two shock tunnels, one of the conventional type, and the other of the free-piston driven type. The geometry of the protrusions, i.e., the height and the deflection angle, was also parametrically varied to study its effect on the hot spot heat flux. The results thus obtained for the flat plate case were compared to existing correlations in the open literature from a similar previous study at a much higher Reynolds number range. Since a mismatch was observed between the results of the current experiments and the existing correlations, a new empirical correlation has been developed to predict the hot spot heat flux, that is valid within the range of flow conditions studied here. A similar attempt was made for the case of the cone model, for which no previous correlations exist in the open literature. However, a global correlation covering the entire range of flow conditions used here could not be formed. A correlation that is valid for just one out of the five flow conditions used here is presented for the cones with sharp and blunt nose tips separately. Schlieren flow visualization was carried out to obtain a better understanding of the shock structures near the protuberances on both models. For most cases, where the protrusion height and deflection angle were large enough to cause flow separation immediately upstream of the protuberance, a separation shock was manifested which deflected some part of the boundary layer above the protuberance, while the rest of the fluid in the boundary layer entered a recirculating region in the separated zone before escaping to the side. Some preliminary computational analysis was conducted which confirmed this qualitatively. However, the quantitative match of surface heat flux between the simulations and experiments were not encouraging. Schlieren visualization revealed that for the flat plate case, the foot of the separation shock was located at a distance of 10.5 to 12 times the protrusion height ahead of it, whereas in the case of the sharp cone, it was at a distance of 9 to 10.5 times the protrusion height. The unsteady nature of the separation shock was also captured and addressed. Some preliminary experiments on boundary layer tripping were also conducted, the results of which have been presented here. From this analysis, it has become evident that a single global correlation cannot be formed which could be used for a wide range of flow conditions to predict the hot spot heat flux in interference interactions. The entire range of conditions that may be encountered during hypersonic flight has to be broken down into sections, and the interference heating pattern should be studied in each of these sections individually. By doing so, a series of different correlations can be formed at the varying flow conditions which will then be available for high-speed aircraft designers.

No description available.

In this thesis I investigate characteristics of turbulent flow in a channel where one of the walls has riblets, superhydrophobic microribs, or a hybrid surface with traditional riblets built on a superhydrophobic microrib surface. PIV measurements are used to find the velocity profile, the turbulent statistics, and shear stress profile in the rectangular channel with one wall having a structured test surface. Both riblets and superhydrophobic surfaces can each provide a reduction in the wall shear stress in a turbulent channel flow. Characterizing the features of the flow using particle image velocimetry (PIV) is the focus of this research. Superhydrophobicity results from the combination of a hydrophobic coating applied to a surface with microrib structures, resulting in a very low surface energy, such that the fluid does not penetrate in between the structures. The micro-rib structures are aligned in the streamwise flow direction. The riblets are larger than the micro-rib structure by an order of magnitude and protrude into the flow. All the test surfaces were produced on silicon wafers using photolithographic techniques. Pressure in the channel is maintained below the Laplace pressure for all testing, creating sustainable air pockets between the microribs. Velocity profiles, turbulent statistics, shear stress profiles, and friction factors are presented. Measurements were acquired for Reynolds numbers ranging from 4.5x10^3 to 2.0x10^4. Modest drag reductions were observed for the riblet surfaces. Substantial drag increase occurred over the superhydrophobic surfaces. The hybrid surfaces showed the greatest drag reduction. Turbulence production was strongly reduced during riblet and hybrid tests.

No description available.

La thèse présente la modélisation, l'étude expérimentale et le contrôle actif des instabilités aéroacoustiques rencontrées en cavité sous écoulement turbulent à faible nombre de Mach. On propose une formulation problème de stabilité de l'interface fluide séparant deux écoulements uniformes de vitesse différente en integrant les effets acoustiques. Les modes d'instabilité de l'interface sont alors étudiés en fonction du nombre de Mach et de la configuration géométrique. Une maquette comportant une cavité de dimensions réglables est ensuite étudiée en soufflerie à l'aide de mesures de pression. Ces données valident en partie l'approche analytique adoptée. On conçoit alors un dispositif de contrôle des modes d'instabilité, appliqué en particulier dans le cas de leur couplage avec l'acoustique de la veine d'essais. Enfin, un système de vélocimétrie par images de particules synchronisé sur les modes d'oscillation permet de valider l'étude théorique et la stratégie de contrôle.

Aerodynamic cavities are common features on hypersonic vehicles which are caused in both large and small scale features like surface defects, pitting, gap in joints etc. In the hypersonic regime, the presence of such cavities alters the flow phenomenon considerably and heating rates adjacent to the discontinuities can be greatly enhanced due to the diversion of flow. Since the 1960s, a great deal of theoretical and experimental research has been carried out on cavity flow physics and heating. However, most of the studies have been done to characterize the effect downstream and within the cavity. In the present study, a series of were carried out in the shock tunnel to investigate the heating characteristics, upstream and on the lateral side of the cavity. Heat flux measurement has been done using indigenously developed high resistance platinum thin film gauges. High resistance gauges, as contrary to the conventionally used low resistance gauges were showing good response to the extremely low heat flux values on a flat plate with sharp leading edge. The experimental measurements of heat done on a flat plate with sharp leading edge using these gauges show good match with theoretical relation by Crabtree et al. Flow visualization using high speed camera with the cavity model and shock structures visualized were similar to reported in supersonic cavity flow. This also goes to state that in spite of the fluctuating shear layer-the main feature of hypersonic flow over a cavity ,reasonable studies can be done within the short test time of shock tunnel. Numerical Simulations by solving the Navier-Stokes equation, using the commercially available CFD package FLUENT 13.0.0 has been done to complement the experimental studies.

This thesis constitutes detailed computational investigation on ow through the passages of a centrifugal compressor used for turbocharging applications. Given the dynamic nature of operation of the turbocharger, it becomes necessary to under- stand the ow that occurs within the blade passages and its e ect on performance. CFD is an established computational technique wherein the ow is dissected to fun- damental levels and a detailed picture is presented, application of this technique with limited and diverse sense towards understanding of ows through a turbocharger compressor has been successfully carried out by many before. This work presented attempts to address many of the lacuna reported and carries forward the work of several researchers to ll in the gaps. The complexity of the geometry of the blade shape poses many challenges in model- ing within the virtual space, an e ective way to overcome the obstacles is presented as a part of this work. Grid generation of the impeller and casing are discussed and adaptive approach is followed with generation of hexahedral grids for the impeller whereas tetrahedral for the casing. Since the grids of the impeller and its casing are di erent, ways of interfacing between the two domains in a CFD environment is discussed. An industry standard implicit 3D RANS solver was used to carry out the simula- tions. The importance of use of boundary conditions for the domain at unsteady operating points is presented in detail. On the choice made for turbulence model that governs the validity of the solution obtained, an extensive literature survey of the relevant topic as applicable for centrifugal compressors is presented and logic of the choice made for the present work is discussed. Menter's two equation SST-k! model emerges as the clear choice to be used even though the di erence in perfor- mance predictions by other turbulence models are insigni cant. Dynamics of ow at optimum design point, surge and choke of the compressor are presented in detail. With the geometry modeled with a tip clearance and the casing included within the simulation environment, it can be seen that the performance predicted is closer to actual at all operating points. A study of behavior of the compressor at extreme o design points is carried out and it can be seen that it depicts the trends that are seen in experimental works available in open literature. The distortion of pressure within the vaneless di user and the inviscid nature of the ow within the volute space are e ectively captured and an in depth analysis is carried out to uncover new patterns. A parametric study involving important geometric features such as the tip clearance and wrap angles are conducted leading to discovery of anomalies. The work summarizes to point out that the investigation carried out with the CFD simulations comprehensively leads to uncovering of ow dynamics within a complex system such as the centrifugal compressor within the limits of numerical analysis.

The context of this thesis is the numerical simulation of turbulent flows at moderate Reynolds numbers and the improvement of the capabilities of an in-house 3D unsteady and incompressible flow solver called SFELES to simulate such flows.<p>In addition to this abstract, this thesis includes five other chapters.<p><p>The second chapter of this thesis presents the numerical methods implemented in the two CFD solvers used as part of this work, namely SFELES and PHASTA.<p><p>The third chapter concentrates on the implementation of a new library called FlexMG. This library allows the use of various types of iterative solvers preconditioned by algebraic multigrid methods, which require much less memory to solve linear systems than a direct sparse LU solver available in SFELES. Multigrid is an iterative procedure that relies on a series of increasingly coarser approximations of the original 'fine' problem. The underlying concept is the following: low wavenumber errors on fine grids become high wavenumber errors on coarser levels, which can be effectively removed by applying fixed-point methods on coarser levels.<p>Two families of algebraic multigrid preconditioners have been implemented in FlexMG, namely smooth aggregation-type and non-nested finite element-type. Unlike pure gridless multigrid, both of these families use the information contained in the initial fine mesh. A hierarchy of coarse meshes is also needed for the non-nested finite element-type multigrid so that our approaches can be considered as hybrid. Our aggregation-type multigrid is smoothed with either a constant or a linear least square fitting function, whereas the non-nested finite element-type multigrid is already smooth by construction. All these multigrid preconditioners are tested as stand-alone solvers or coupled with a GMRES (Generalized Minimal RESidual) method. After analyzing the accuracy of the solutions obtained with our solvers on a typical test case in fluid mechanics (unsteady flow past a circular cylinder at low Reynolds number), their performance in terms of convergence rate, computational speed and memory consumption is compared with the performance of a direct sparse LU solver as a reference. Finally, the importance of using smooth interpolation operators is also underlined in this work.<p><p>The fourth chapter is devoted to the study of subgrid scale models for the large eddy simulation (LES) of turbulent flows.<p>It is well known that turbulence features a cascade process by which kinetic energy is transferred from the large turbulent scales to the smaller ones. Below a certain size, the smallest structures are dissipated into heat because of the effect of the viscous term in the Navier-Stokes equations.<p>In the classical formulation of LES models, all the resolved scales are used to model the contribution of the unresolved scales. However, most of the energy exchanges between scales are local, which means that the energy of the unresolved scales derives mainly from the energy of the small resolved scales.<p>In this fourth chapter, constant-coefficient-based Smagorinsky and WALE models are considered under different formulations. This includes a classical version of both the Smagorinsky and WALE models and several scale-separation formulations, where the resolved velocity field is filtered in order to separate the small turbulent scales from the large ones. From this separation of turbulent scales, the strain rate tensor and/or the eddy viscosity of the subgrid scale model is computed from the small resolved scales only. One important advantage of these scale-separation models is that the dissipation they introduce through their subgrid scale stress tensor is better controlled compared to their classical version, where all the scales are taken into account without any filtering. More precisely, the filtering operator (based on a top hat filter in this work) allows the decomposition u' = u - ubar, where u is the resolved velocity field (large and small resolved scales), ubar is the filtered velocity field (large resolved scales) and u' is the small resolved scales field. <p>At last, two variational multiscale (VMS) methods are also considered.<p>The philosophy of the variational multiscale methods differs significantly from the philosophy of the scale-separation models. Concretely, the discrete Navier-Stokes equations have to be projected into two disjoint spaces so that a set of equations characterizes the evolution of the large resolved scales of the flow, whereas another set governs the small resolved scales. <p>Once the Navier-Stokes equations have been projected into these two spaces associated with the large and small scales respectively, the variational multiscale method consists in adding an eddy viscosity model to the small scales equations only, leaving the large scales equations unchanged. This projection is obvious in the case of a full spectral discretization of the Navier-Stokes equations, where the evolution of the large and small scales is governed by the equations associated with the low and high wavenumber modes respectively. This projection is more complex to achieve in the context of a finite element discretization. <p>For that purpose, two variational multiscale concepts are examined in this work.<p>The first projector is based on the construction of aggregates, whereas the second projector relies on the implementation of hierarchical linear basis functions.<p>In order to gain some experience in the field of LES modeling, some of the above-mentioned models were implemented first in another code called PHASTA and presented along with SFELES in the second chapter.<p>Finally, the relevance of our models is assessed with the large eddy simulation of a fully developed turbulent channel flow at a low Reynolds number under statistical equilibrium. In addition to the analysis of the mean eddy viscosity computed for all our LES models, comparisons in terms of shear stress, root mean square velocity fluctuation and mean velocity are performed with a fully resolved direct numerical simulation as a reference.<p><p>The fifth chapter of the thesis focuses on the numerical simulation of the 3D turbulent flow over a re-entry Apollo-type capsule at low speed with SFELES. The Reynolds number based on the heat shield is set to Re=10^4 and the angle of attack is set to 180º, that is the heat shield facing the free stream. Only the final stage of the flight is considered in this work, before the splashdown or the landing, so that the incompressibility hypothesis in SFELES is still valid.<p>Two LES models are considered in this chapter, namely a classical and a scale-separation version of the WALE model. Although the capsule geometry is axisymmetric, the flow field in its wake is not and induces unsteady forces and moments acting on the capsule. The characterization of the phenomena occurring in the wake of the capsule and the determination of their main frequencies are essential to ensure the static and dynamic stability during the final stage of the flight. <p>Visualizations by means of 3D isosurfaces and 2D slices of the Q-criterion and the vorticity field confirm the presence of a large meandering recirculation zone characterized by a low Strouhal number, that is St≈0.15.<p>Due to the detachment of the flow at the shoulder of the capsule, a resulting annular shear layer appears. This shear layer is then affected by some Kelvin-Helmholtz instabilities and ends up rolling up, leading to the formation of vortex rings characterized by a high frequency. This vortex shedding depends on the Reynolds number so that a Strouhal number St≈3 is detected at Re=10^4.<p>Finally, the analysis of the force and moment coefficients reveals the existence of a lateral force perpendicular to the streamwise direction in the case of the scale-separation WALE model, which suggests that the wake of the capsule may have some <p>preferential orientations during the vortex shedding. In the case of the classical version of the WALE model, no lateral force has been observed so far so that the mean flow is thought to be still axisymmetric after 100 units of non-dimensional physical time.<p><p>Finally, the last chapter of this work recalls the main conclusions drawn from the previous chapters. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

❚✐t❧❡✿ ◗✉❛♥t✉♠ t✉r❜✉❧❡♥❝❡ ✐♥ s✉♣❡r✢✉✐❞ ❤❡❧✐✉♠ st✉❞✐❡❞ ❜② ♣❛rt✐❝❧❡ tr❛❝❦✐♥❣ ✈❡❧♦❝✐♠❡tr② ✈✐s✉❛❧✐③❛t✐♦♥ t❡❝❤♥✐q✉❡ ❆✉t❤♦r✿ ❘◆❉r✳ ❉❛♥✐❡❧ ❉✉❞❛ ❉❡♣❛rt♠❡♥t✿ ❉❡♣❛rt♠❡♥t ♦❢ ▲♦✇ ❚❡♠♣❡r❛t✉r❡ P❤②s✐❝s ❙✉♣❡r✈✐s♦r✿ ♣r♦❢✳ ❘◆❉r✳ ▲❛❞✐s❧❛✈ ❙❦r❜❡❦✱ ❉r❙❝✳ ❆❜str❛❝t✿ ❚❤❡ P❛rt✐❝❧❡ ❚r❛❝❦✐♥❣ ❱❡❧♦❝✐♠❡tr② ✈✐s✉❛❧✐③❛t✐♦♥ t❡❝❤♥✐q✉❡ ✉s✐♥❣ ♠✐❝✲ r♦♠❡t❡r s✐③❡ s♦❧✐❞ ❞❡✉t❡r✐✉♠ ♣❛rt✐❝❧❡s ❛s tr❛❝❡rs ❤❛s ❜❡❡♥ ❛♣♣❧✐❡❞ t♦ st✉❞② ♦s❝✐❧✲ ❧❛t♦r② ✢♦✇s ♦❢ ❍❡ ■■✱ ✇❤✐❝❤ ✐s ❛ q✉❛♥t✉♠ ✢✉✐❞ ✇✐t❤ q✉❛♥t✐③❡❞ ✈♦rt✐❝✐t②✱ ❛s ✇❡❧❧ ❛s ✢♦✇s ♦❢ ❍❡ ■✱ ✇❤✐❝❤ ✐s ❛ ❝❧❛ss✐❝❛❧ ✈✐s❝♦✉s ❧✐q✉✐❞✱ ❢♦❝✉s✐♥❣ ♦♥ t❤❡ s✐♠✐❧❛r✐t✐❡s ❛♥❞ ❞✐✛❡r❡♥❝❡s ❜❡t✇❡❡♥ t❤❡ q✉❛♥t✉♠ ❛♥❞ ❝❧❛ss✐❝❛❧ ✢♦✇s✳ ❚❤r❡❡ ❡①♣❡r✐♠❡♥ts ❛r❡ ❞❡s❝r✐❜❡❞✿ t❤❡ ✢♦✇ ♣❛st ❛ ❧❛r❣❡✲❛♠♣❧✐t✉❞❡ ❧♦✇✲❢r❡q✉❡♥❝② ♦s❝✐❧❧❛t✐♥❣ ♦❜st❛❝❧❡ ✐♥ t❤❡ ❢♦r♠ ♦❢ ❛ ♣r✐s♠❀ t❤❡ st❡❛❞② str❡❛♠✐♥❣ ✢♦✇ ❞✉❡ t♦ ❛ s♠❛❧❧✲❛♠♣❧✐t✉❞❡ ❧❛r❣❡✲ ❢r❡q✉❡♥❝② ♦s❝✐❧❧❛t✐♥❣ q✉❛rt③ t✉♥✐♥❣ ❢♦r❦ ✲ ❛ ✇✐❞❡❧② ✉s❡❞ t♦♦❧ t♦ st✉❞② q✉❛♥t✉♠ t✉r❜✉❧❡♥❝❡❀ ❛♥❞ t❤❡ ♣r♦❞✉❝t✐♦♥ ♦❢ ❝❛✈✐t❛t✐♦♥ ✐♥ t❤❡ ✈✐❝✐♥✐t② ♦❢ ❛ ❢❛st✲♦s❝✐❧❧❛t✐♥❣ t✉♥✐♥❣ ❢♦r❦✳ ❚❤❡ ♠❛✐♥ ♦✉t❝♦♠❡ ✐s t❤❡ ♦❜s❡r✈❛t✐♦♥ t❤❛t t❤❡s❡ ✢♦✇s ❛r❡ s✐♠✐❧❛r ✐♥ ❍❡ ■ ❛♥❞ ✐♥ ❍❡ ■■ ❛t ❧❛r❣❡ ❧❡♥❣t❤✲s❝❛❧❡s✱ ✇❤❡r❡❛s ❛t s♠❛❧❧ s❝❛❧❡s✱ t❤❡② ❡①❤✐❜✐t t♦t❛❧❧② ❞✐✛❡r❡♥t st❛t✐st✐❝❛❧ ♣r♦♣❡rt✐❡s✳ ▼♦r❡♦✈❡r✱ ✐♥ ❍❡ ■■✱ t❤❡s❡ s♠❛❧❧ s❝❛❧❡ st❛✲ t✐st✐❝❛❧ ♣r♦♣❡rt✐❡s ❛r❡ ✉♥✐✈❡rs❛❧ ✐♥ t❤❛t t❤❡② ❞♦ ♥♦t ❞❡♣❡♥❞ ♦♥ t❤❡ t②♣❡ ♦❢ t❤❡ ✐♠♣♦s❡❞ ♠❡❛♥ ✢♦✇ ♦❢ t❤❡ s✉♣❡r✢✉✐❞ ❛♥❞ ♥♦r♠❛❧...