L’allumage est un point essentiel dans le dimensionnement des moteurs-fusées, et il nécessite de prendre en compte plusieurs phénomènes physiques très distincts qui sont autant de challenges numériques. Le premier point abordé pendant cette thèse est la modélisation et la simulation par Simulation aux Grandes Échelles d’un allumeur de type VINCI. Des gaz chauds, riches en oxygène, sont délivrés de façon supersonique dans une chambre remplie d’hydrogène faisant apparaître un jet fortement sous-détendu et de multiples interactions choc/choc ou choc/flamme. Les premiers instants du processus d’allumage sont ici détaillés. Le second point abordé est la modélisation et la simulation numérique de la combustion H2/O2 à haute pression. En particulier, les effets d’une diffusion non-idéale sont étudiés dans le cas de flammes de prémélange 1D et sur la configuration 2D de type ‘splitter plate’. Un impact de la modélisation sur les espèces produites et le champ de température est ici mis en lumière. / Ignition is a key point in the design of liquid rocket engine (LRE), and it requires to take into account several distinct physical phenomena that constitute numerical challenges. The first point addressed during this thesis is the modeling and simulation using Large Eddy Simulation of a LRE igniter in a configuration close to VINCI rocket engine. The hot gases from the igniter, rich in oxygen, are delivered at supersonic speeds in a chamber filled with hydrogen. Such configuration creates under-expanded jets with multiple shock/shock or shock/flame interactions. A focus is done on the ignition process. The second point addressed is the modeling and simulation of high pressure H2/O2 combustion which occurs. In particular, the effects of non-ideal diffusion are studied through a 1D premixed flames and a 2D splitter plate configuration. An impact of modeling on the species produced and the temperature field is highlighted.

The sound produced by military aircraft is dominated by noise generated by the turbulent mixing of the jetted exhaust with the ambient air. This jet noise has the potential to annoy the community and pose a hearing loss risk for military personnel. The goal of this dissertation is to characterize spatiospectral features in the field produced by full-scale military aircraft that are not traditionally seen at the laboratory scale and identify potential noise mechanisms for these features. Measurements of two military aircraft jet noise fields are found to be best described as a superposition of spatiospectral lobes, whose relative amplitudes dictate the overall directivity at each engine power. Near-field acoustical holography techniques are applied to one of the military aircraft measurements to characterize the behavior of the lobes as a function of engine power. The simulated jet noise of a highly heated laboratory-scale jet is then analyzed to compare with the military aircraft measurement and is found to only partially contain the spatiospectral lobe phenomenon. Application of near to far field coherence tracing and near-field acoustical holography to the simulations provides validation of the methods used on the military aircraft and illuminate potential source mechanisms that may explain the presence of the spatiospectral lobes.

Réaliser un système de combustion à petite échelle reste aujourd’hui un défi. L’augmentation du rapport surface/volume favorise les pertes thermiques, contribue à la diminution du temps de séjour et limite la turbulence. Le premier objectif de cette thèse est de comprendre les phénomènes physiques intervenant dans un brûleur centimétrique tourbillonnaire de 8 x10 x 8 mm3 (millimètre cube) et mettre au point des outils numériques adaptés. L’écoulement réactif méthane/air est étudié au moyen de simulations numériques LES. La combustion ne consomme pas l’intégralité du carburant, entraînant un rendement de combustion de l’ordre de 50% et d’importantes émissions de polluants. Le deuxième objectif est d’adapter les performances de ce brûleur. L’enrichissement en hydrogène a montré une amélioration sensible du rendement et une réduction des émissions polluantes. Plusieurs configurations géométriques de la chambre ont aussi été étudiées, ce qui a permis de dégager des axes d’améliorations. / Designing a meso-scale combustion system remains a challenging scientific and technological issue. Increasing the surface-to-volume ratio promotes wall heat losses, reduces the residence time and turbulence intensity. The main objective of this thesis is to understand the physical phenomena involved in the centimetre-sized asymmetric whirl cubic burner of 8 x 10 x 8 mm3 (millimètre cube) and develop specific adapted numerical tools. The methane/air reactive flow is studied using detailed LES. While fuel and air are injected separately, combustion takes place in the premixed regime. However combustion is far from being complete, causing low combustion efficiency and significant emissions of pollutants. The second objective is to adapt in the best possible way the performances of this burner. Hydrogen enrichment of the fuel mixture showed significant efficiency enhancement and reduced pollutant emissions. Several other combustor geometries are also studied, paving the way for future improvement.

The main objective of this study was to evaluate the performance of three Large Eddy Simulation (LES) subgrid scale (SGS) models on a strongly heated, low Mach number upward gas flow in a vertical pipe with forced convection. The models chosen for this study were the Smagorinsky-Lilly Dynamic model (SLD), the Kinetic Energy Transport model (KET), and the Wall-Adaptive Local-Eddy viscosity model (WALE). The used heating rate was sufficiently large to cause properties to vary significantly in both the radial and streamwise directions. All simulations were carried out using the commercial software FLUENT. The effect of inlet turbulence generation techniques was considered as well. Three inlet turbulence generation techniques were compared, namely, the Spectral Synthesizer Method (SSM), the Vortex Method (VM), and the Generator (GEN) technique. A user-defined function (UDF) was written to implement the GEN technique into the solver; the SSM and VM techniques were already build-in. All simulation and solver settings were validated by performing computational simulations of isothermal fully developed pipe flow and results were compared to available experimental and Direct Numerical Simulation (DNS) data. For isothermal boundary conditions, among the three inlet turbulence generation techniques, the GEN technique produced results which best matched the experimental and DNS results. All three LES SGS models performed equally well when coupled with the GEN technique for the study of isothermal pipe flow. However, all models incorrectly predicted the behavior of radial and circumferential velocity fluctuations near the wall and the GEN technique proved to be the most computationally expensive. For simulations with longer computational domain, the effect of the inlet turbulence generation technique diminishes. However, results suggest that both the SLD and KET models need shorter computational domains to recover proper LES behavior when coupled with the VM technique in comparison to the WALE SGS model with the same turbulence inlet generation technique. For high heat flux simulations all SGS models were coupled with the VM technique to decrease the computational effort to obtain statistically steady-state solution. For comparative purposes, one simulation was carried out using the WALE and GEN techniques. All simulations equally significantly underpredicted the streamwise temperature distribution along the pipe wall as well as in the radial directions at various streamwise locations. These effects are attributed to the overpredicted streamwise velocity components and incorrect behavior of both the radial and circumferential velocity components in the near wall region for all subgrid scale models.

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第19685号 / 工博第4140号 / 新制||工||1639(附属図書館) / 32721 / 京都大学大学院工学研究科機械理工学専攻 / (主査)教授 小森 悟, 教授 中部 主敬, 教授 稲室 隆二 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

An atmospheric pressure swirl-stabilized methane-air burner has been developed as a test platform for nanosecond repetitively pulsed (NRP) discharge plasma-assisted combustion research. Qualitative flame and plasma discharge characterizations were conducted with high-speed video and low-light ICCD imagery, along with a modal acoustic analysis of the entire assembly. A large eddy simulation (LES) of the burner was created using the commercial solver Ansys Fluent to investigate the plasma effects on swirl-stabilized turbulent combustion. A modified version of the solver's premixed combustion mechanism is presented along with a phenomenological plasma discharge model to simulate plasma-assisted combustion. Cold flow particle image velocimetry (PIV) data were collected to validate the non-reacting flow field and assess non-reacting NRP discharge effects. Optical emission spectroscopy (OES) measurements of the second positive system (SPS) of nitrogen mapped temperature characteristics of NRP discharge bursts for comparison to time-resolved simulation data. Finally, time-averaged CH* chemiluminescence data were collected to qualitatively assess the effects of plasma on the experimental burner and simulated flame structure. Overall, the phenomenologically-based combustion mechanism proposed in this work shows good agreement with several experimental observations and provides a promising framework for future plasma-assisted combustion modeling.

No description available.

The implementation of a hybrid spectral/finite-element discretization on the unsteady, incompressible, Navier-Stokes equations with a semi-implicit time-stepping method, an explicit treatment of the advective terms, and an implicit treatment of the pressure and viscous terms leads to an algorithm capable of calculating 3D flows over complex 2D geometries. This also results in multiple Fourier mode linear systems which must be solved at every timestep, which naturally leads to two parallelization approaches: Fourier space partitioning, where each processor individually and simultaneously solves a linear system, and physical space partitioning, where all processors collectively solve each linear system, sequentially advancing through Fourier modes. These two parallelization approaches are compared based upon computational cost using multiple solvers: direct sparse LU, smoothed aggregation AMG, and single-level ILUT preconditioned GMRES; and on two supercomputers of different memory architecture(distributed and shared memory). This study revealed Fourier space partitioning outperforms physical space partitioning in all problems analyzed, and scales more efficiently as well. These differences were more dramatic on the distributed memory platform than the shared memory platform. Another study compares the previously mentioned solvers along with one additional solver, pointwise AMG, in Fourier space partitioning without parallelization to better understand computational scaling for problems with large meshes. It was found that the direct sparse LU solver performed well in terms of computational time, scaled linearly, but had very high memory usage which scaled in a super-linear manner. The single-level ILUT preconditioned GMRES solver required the least amount of memory, which also scaled linearly, but only had acceptable performance in terms of computational time for coarse meshes. Both AMG methods scaled linearly in both memory usage and time, and were comparable to the direct sparse LU solver in terms of computational time. The results of these studies are particularly useful for implementation of this algorithm on challenging and complex flows, especially direct numerical and large-eddy simulations. Reducing computational cost allows the analysis and understanding of more flows of practical interest.

The aim of the present study is to develop subgrid-scale models that are relevant for complex flows and combustion. A stochastic model based on a stochastic Smagorinsky constant with adjustable variance and time scale is proposed. The stochastic model is shown to provide for backscatter of both kinetic energy and scalar variance without causing numerical instabilities. A new subgrid-scale scalar flux model is developed using the same kind of methodology that leads to the explicit algebraic scalar flux model, EASFM, for RANS. The new model predicts the anisotropy of the subgrid-scales in a more realistic way than the eddy diffusion model. Both new models were tested in rotating homogeneous shear flow with a passive scalar. Rogallo’s method of moving the frame with the mean flow to enable periodic boundary conditions was used to simulate homogeneous shear flow. / QC 20101119

This thesis describes the numerical predictions of turbine film cooling interactions using Large Eddy Simulations. In most engineering industrial applications, the Reynolds-Averaged Navier-Stokes equations, usually paired with two-equation models such as k-Greek lowercase letter epsilon] or k-Greek lowercase letter omega], are utilized as an inexpensive method for modeling complex turbulent flows. By resolving the larger, more influential scale of turbulent eddies, the Large Eddy Simulation has been shown to yield a significant increase in accuracy over traditional two-equation RANS models for many engineering flows. In addition, Large Eddy Simulations provide insight into the unsteady characteristics and coherent vortex structures of turbulent flows. Discrete hole film cooling is a jet-in-cross-flow phenomenon, which is known to produce complex turbulent interactions and vortex structures. For this reason, the present study investigates the influence of these jet-crossflow interactions in a time-resolved unsteady simulation. Because of the broad spectrum of length scales present in moderate and high Reynolds number flows, such as the present topic, the high computational cost of Direct Numerical Simulation was excluded from possibility.