No description available.

No description available.

This thesis deals with the Navier-Stokes equations for real, compressible fluid with first and second viscosity. The method of expansion into a series of eigenmodes of vibration is chosen to solve the Navier-Stokes equations. The general relations of the expansion coefficients and the eigenmode shapes are given for 1D-, 2D- and 3D-flows. The specific formulae of eigenmodes of vibration are determined only for 1D-flow. The final form of the pressure function is analysed using Matlab software.

No description available.

Dans ce manuscrit, nous parlons des matériaux composites à matrice céramiques (CMCs) qui sont envisagés pour intégrer les chambres de combustion de futurs moteurs aéronautiques civils. Pour faire face des conditions extrêmes, ces matériaux possèdent la particularité de s’auto-protéger vis-à-vis de l’oxydation par la formation d’un oxyde passivant qui limite la diffusion des espèces oxydantes au sein des fissures matricielles. Nous modélisons l’écoulement d’un oxyde dans une fissure par l’équation de Navier-Stokes, puis les mettons sous forme non dimensionnelles, et les dérivations de deux types de modèles sont intéressantes : les modèles de Saint-Venant et les modèles de lubrification. Ensuit nous nous engageons à chercher l’existence de solution faible de l’approximation de lubrification d’ordre 4 obtenue précédent dans le cas uni-dimensionnel. Enfin nous précisons la limite entre les équations de Saint-Venant et l’équation de lubrification. / In this work, we are interested in the ceramic matrix composite materials(CMCs) who will be used to integrate the combustion chambers of future civil aeronautical engines. To face extreme conditions, these materials possess the peculiarity to auto-protect itself towards the oxidation by the formation of an oxide passivate which limits the distribution of the oxidizing species within the matrix cracks. We model the flow of an oxide in a crack by the Navier-Stokes equation, then put them under an asymptotic analysis in order to get two types of asymptotic models : models of Saint-Venant (Shallow water model) and lubrication models. Next we are interested in looking for the existence of weak solution to the one-dimensional approximated lubrication equation of order 4 obtained before. Finally we talk about the limit between the Saint-Venant equations and the lubrication equation.

Clouds are ubiquitous and ever-changing feature of the outdoors. They are an integral factor in Earth's weather systems. Component of water circulation in a nature and a strong indicator of weather patterns and changes. Clouds are important component of the visual simulation of any outdoor scene, but the complexity of cloud formation, dynamics, and light interaction makes cloud simulation and rendering difficult in real time. In an interactive flight simulation, users would like to fly in and around realistic, volumetric clouds. I Present, two main ways of representation cloud in computer graphic, where one way use for modelling Perlin noise and second one use for modelling fluid dynamic system.

Nesta tese é apresentado um esquema novo de alta resolução upwind (denominado TDPUS-C3) para reconstrução de fluxos numéricos para leis de conservação não lineares e problemas relacionados em DFC. O esquema é baseado nos critérios de estabilidade CBC e TVD e desenvolvido utilizando condições de diferenciabilidade \'C POT. 3\'. Além disso, é realiozada a implementação da associação do esquema TDPLUS-C3 com a modelagem de turbulência RNG \'\\kappa - \\epsilon\'. O propósito é obter soluções numéricas de sistemas hiperbólicos de leis de conservação para dinâmica dos gases e equações de Navier-Stokes para escoamento incompreensível de fluidos newtonianos e não newtonianos (viscoelásticos). Fazendo o uso do esquema TDPUS-C3, a precisão global dos métodos numéricos é verificada acessando o erro em problemas teste (benchmark) 1D e 2D. Um estudo comparativo entre os resultados do esquema TDPUS-C3 e os esquemas upwind convencionais para leis de conservação hiperbólicas complexas é também realizado. A Associação das modelagens numéricas (upwinding mais RNG \'\\kappa - \\epsilon\') é , então, examinada na simulação de escoamentos turbulentos de fluidos newtonianos envolvendo superfícies livres móveis, usando a metodologia URANS. No geral, em termos do comportamento global, concordância satisfatória é observada / In this thesis, a new high-resolution upwind scheme (named TDPUS-C3) for reconstruction of numerical fluxes for nonlinear conservation laws and related CFD problems in presented. The scheme is based on CBC and TVD stability criteria and developed by employing differentiability condictions (\'C POT. 3\'). In additon, the implementation of an association of the TDPUS-C3 scheme with the RNG \'\\kappa - \\epsilon\' turbulence modelling is also performed. The purpose is to obtain numerical solutions of systems of hyperbolic conservation laws for gas dynamics and Navier-Stokes equations for incompressible flow of Newtonian and non-Newtonian (viscoelstic) fluids. By using the TDPUS-C3 scheme, the global accuracy of the numerical methods is verified by assessing the error on 1D and 2D benchmark test cases. A comparative study between the TDPUS-C3 scheme and convectional upwind schemes to solve standard and complex hyperbolic conservation laws is also accomplished. The association of the numerical modelling (upwinding plus RNG \'\\kappa - epsilon\') is then examined in the simulation of turbulent Newtonian fluid flows involving moving free surfaces, by using URANS methodology. Overall, satisfactory agreement is found in terms of the overall behaviour

The fluid flow and tracer transport in a single rock fracture during shear processes has been an important issue in rock mechanics and is investigated in this thesis using Finite Element Method (FEM) and streamline particle tracking method, considering evolutions of aperture and transmissivity with shear displacement histories under different normal stresses, based on laboratory tests. The distributions of fracture aperture and its evolution during shear were calculated from the initial aperture fields, based on the laser-scanned surface roughness features of replicas of rock fracture specimens, and shear dilations measured during the coupled shear-flow-tracer tests in laboratory performed using a newly developed testing apparatus in Nagasaki University, Nagasaki, Japan. Three rock fractures of granite with different roughness characteristics were used as parent samples from which nine plaster replicas were made and coupled shear-flow tests was performed under three normal loading conditions (two levels of constant normal loading (CNL) and one constant normal stiffness (CNS) conditions). In order to visualize the tracer transport, transparent acrylic upper parts and plaster lower parts of the fracture specimens were manufactured from an artificially created tensile fracture of sandstone and the coupled shear-flow tests with fluid visualization was performed using a dye tracer injected from upstream and a CCD camera to record the dye movement. A special algorithm for treating the contact areas as zero-aperture elements was used to produce more accurate flow field simulations by using FEM, which is important for continued simulations of particle transport, but was often not properly treated in literature. The simulation results agreed well with the flow rate data obtained from the laboratory tests, showing that complex histories of fracture aperture and tortuous flow channels with changing normal stresses and increasing shear displacements, which were also captured by the coupled shear-flow tests of fracture specimens with visualization of the fluid flow. From the obtained flow velocity fields, the particle transport was predicted by the streamline particle tracking method with calculated flow velocity fields (vectors) from the flow simulations, obtaining results such as flow velocity profiles, total flow rates, particle travel time, breakthrough curves and the Péclet number, Pe, respectively. The fluid flow in the vertical 2-D cross-sections of a rock fracture was also simulated by solving both Navier-Stokes (NS) and Reynolds equations, and the particle transport was predicted by streamline particle tracking method. The results obtained using NS and Reynolds equations were compared to illustrate the degree of the validity of the Reynolds equation for general applications in practice since the later is mush more computationally efficient for large scale problems. The flow simulation results show that the total flow rate and the flow velocity predicted by NS equations are quite different from that as predicted by the Reynolds equation. The results show that a roughly 5-10 % overestimation on the flow rate is produced when the Reynolds equation is used, and the ideal parabolic velocity profiles defined by the local cubic law, when Reynolds equation is used, is no longer valid, especially when the roughness feature of the fracture surfaces changes with shear. These deviations of flow rate and flow velocity profiles across the fracture aperture have a significant impact on the particle transport behavior and the associated properties, such as the travel time and Péclet number. The deviations increase with increasing flow velocity and become more significant when fracture aperture geometry changes with shear. The scientific findings from these studies provided new insights to the physical behavior of fluid flow and mass transport in rock fractures which is the scientific basis for many rock mechanics problems at the fundamental level, and with special importance to rock engineering problems such as geothermal energy extraction (where flow rate in fractures dominates the productivity of a geothermal energy reservoir) and nuclear waste repositories (where radioactive nuclides transport through fractures dominates the final safety evaluations) in fractured rocks. / Vätskeflödet och spårämnestransporten i en enskild bergsspricka under skjuvningsprocesser har varit ett viktigt ämne inom bergmekanik. I denna avhandling undersöks ämnet med hjälp av finita element metoden (FEM) och en strömlinjebaserad partikelspårningsmetod. Hänsyn tas till utveckling av öppningar och transmissivitet med skjuvningens förflyttningshistoria under olika normala belastningar baserat på laboratorietester. Fördelningen av spricköppningar och deras utveckling under skjuvning beräknades från de initiala öppningsfälten baserat på det laserscannade provets ytas grovhetskännetecken sam tskjuvningsöppningar uppmätta under de kopplade skjuvning-flöde-spårämneslaboratorietesterna som utförts med nyutvecklad testapparatur i Nagasaki Universitet i Nagasaki, Japan. Tre bergssprickor i granit med olika grovhetskarakteristika användes som utgångsprover från vilka nio gipskopior gjordes. Kopplade skjuvning-flödes tester utfördes sedan under tre normala belastningstillstånd (två nivåer med konstant normal last (KNL) och en konstant normal styvhetstillstånd (KNS). För att visualisera spårämnestransporten tillverkades en transparent övre del av sprickproverna av akryl och en nedre del av gipsbaserat på en kostgjord spänningsspricka i sandsten och de kopplade skjuvning-flödes testerna med vätskevisualisering utfördes med färgspårämne injekterat uppströms och en CCD kamera monterad ovanför för att registrera färgens rörelse. En särskild algoritm användes för att behandla kontaktytorna som nollöppningsämnen användes för att åstadkomma mer exakta flödesfältssimuleringar med FEM. Detta är viktigt för kontinuerliga simuleringar av partikelflöden men uppmärksammas oftast inte tillräckligt i litteraturen. Simuleringsresultaten överensstämde väl med de flödesnivådata som erhölls från laboratorietesterna vilket visade att komplexa historier av spricköppningar och invecklade flöden överensstämde med ändrade normala belastningar och ökande skjuvningsförflyttningar, vilket även fångades av de kopplade skjuvning-flödestesterna av sprickproverna genom visualisering av vätskeflödet. Från de erhållna flödesfälten förutsågs partikeltransporten genom en strömlinjebaserad partikelspårningsmetod med kalkylerade flödeshastighetsfält (vektorer) från flödessimuleringarna genom vilka resultat som flödeshastighetsprofiler, totala flödesnivåer,partikeltransporttid, genombrottskurvor samt Pécletnumret, Pe, erhölls. Vätskeflödet i det vertikala tvådimensionella tvärsnittet av en bergsspricka simulerades även genom att både Navier-Stokes (NS) och Reynoldsekvationerna löstes och partikeltransporten förutsågs genom den strömlinjebaserade partikelspårningsmetoden. Resultaten som erhöllsmed NS och Reynoldsekvationerna jämfördes för att illustrera graden av tillförlitlighet för Reynoldsekvationen för allmänna tillämpningar i praktiken då den senare är betydligt mer beräkningseffektiv för storskaliga problem. Resultaten från flödessimuleringarna visar att den totala flödesnivån och den totala flödeshastigheten förutsedda med NS ekvationer är helt annorlunda motsvarande värden som förutsågs med Reynoldsekvationen. Resultaten visar att en ca 5-10 % för hög uppskattning av flödesnivån erhålls då Reynoldsekvationen används och de ideala parabola hastighetsprofilerna, som definieras av den lokala kubiklagen när Reynoldsekvationen används, inte längre är giltiga särskilt när sprickytornas grovhetskarakteristika ändras med skjuvning. De här avvikelserna i flödesnivå och flödeshastighetsprofiler längs med spricköppningen har en betydande påverkan på partikeltransportuppträdande och de tillhörande egenskaperna såsom rörelsetid och Pécletnummer. Avvikelserna ökar med ökande flödeshastighet och blir mer signifikanta när spricköppningarnas geometri ändras med skjuvning. Forskningsresultaten från dessa studier gav nya insikter i de fysiska uppträdandet av vätskeflöde och masstransporter i bergssprickor vilket är den vetenskapliga basen för många bergmekanikproblem på grundläggande nivå och som har särskild vikt för bergstekniksproblem såsom geotermisk energiutvinning (där flödesnivå i sprickor dominerar produktiviteten för en geotermisk energikälla) och kärnavfallsförvaringsplatser (där transporten av radioaktiva nuklider genom sprickor dominerar den slutgiltigasäkerhetsutvärderingen) i sprickigt berg. / QC 20100803

A finite-volume procedure, comprising a gradient-reconstruction technique and a multidimensional limiter, has been proposed for upwind algorithms on unstructured grids. The high-resolution strategy, with its inherent dependence on a wide computational stencil, does not suffer from a catastrophic loss of accuracy on a grid with poor connectivity as reported recently is the case with many unstructured-grid limiting procedures. The continuously-differentiable limiter is shown to be effective for strong discontinuities, even on a grid which is composed of highly-distorted triangles, without adversely affecting convergence to steady state. Numerical experiments involving transient computations of two-dimensional scalar convection to steady-state solutions of Euler and Navier-Stokes equations demonstrate the capabilities of the new procedure.

This thesis is concerned with the study of parallel mixing of two dissimilar gases under compressible conditions in the confined environment. A number of numerical studies are reported in the literature for the compressible mixing of two streams of gases where (1) both the streams are of similar gases at the same temperatures, (2) both the streams are at different temperatures with similar gases, and (3) dissimilar gases are with nearly equal temperatures. The combination of dissimilar gases at large temperature difference, mixing under compressible conditions, as in the case of scramjet propulsion, has not been adequately addressed numerically. Also many of the earlier studies have used two dimensional numerical simulation and showed good match with the experimental results on mixing layers that are inherently three dimensional in nature. In the present study, both two-dimensional (2-d) and three dimensional (3-d) studies are reported and in particular the effect of side wall on the three dimensionality of the flow field is analyzed, and the reasons of the good match of two dimensional simulations with experimental results have been discussed. Both two dimensional and three dimensional model free simulations have been conducted for a flow configuration on which experimental results are available. In this flow configuration, the mixing duct has a rectangular cross section with height to width ratio of 0.5. In the upper part of the duct hydrogen gas at a temperature of 103 K is injected through a single manifold of two Ludweig tubes and in the lower part of the duct nitrogen gas at a temperature of 2436 K is supplied through an expansion tube, both the gases are at Mach numbers of 3.1 and 4.0 respectively. Measurements in the experiment are limited to wall pressures and heat flux. The choice of this experimental condition gives an opportunity to study the effect of large temperature difference on the mixing of two dissimilar gases with large molecular weights under compressible conditions. Both two dimensional and three dimensional model free simulations are carried out using higher order numerical scheme (4th order spatial and 2nd order temporal) to understand the structure and evolution of supersonic confined mixing layer of similar and dissimilar gases. Two dimensional simulations are carried out by both SPARK (finite difference method) and OpenFOAM (finite volume method based open source software that was specially picked out and put together), while 3D model free simulations are carried out by OpenFOAM. A fine grid structure with higher grid resolution near the walls and shear layer is chosen. The effect of forcing of fluctuations on the inlet velocity shows no appreciable change in the fully developed turbulent region of the flow. The flow variables are averaged after the attainment of statistical steady state established through monitoring the concentration of inert species introduced in the initial guess. The effect of side wall on the flow structure on the mixing layer is studied by comparing the simulation results with and without side wall. Two dimensional simulations show a good match for the growth rate of shear layer and experimental wall pressures. Three dimensional simulations without side wall shows 14% higher growth rate of shear layer than that of two dimensional simulations. The wall pressures predicted by these three dimensional simulations are also lower than that predicted using two dimensional simulations (6%) and experimental (9%) results in the downstream direction of the mixing duct. Three dimensionality of the flow is thought of as a cause for these differences. Simulations with the presence of side wall show that there is no remarkable difference of three dimensionality of the flow in terms of the variables and turbulence statistics compared to the case without side walls. However, the growth rate of shear layer and wall surface pressures matches well with that predicted using two dimensional simulations. It has been argued that this good match in shear layer growth rate occurs due to formation of oblique disturbances in presence of side walls that are considered responsible for the decrease in growth rate in 3-d mixing layers. The wall pressure match is argued to be good because of hindrance from side wall in the distribution of momentum in third direction results in higher wall pressure. The effect of dissimilar gases at large temperature difference on the growth rate reduction in compressible conditions is studied. Taking experimental conditions as baseline case, simulations are carried out for a range of convective Mach numbers. Simulations are also carried out for the same range of convective Mach numbers considering the mixing of similar gases at the same temperature. The normalized growth rates with incompressible counterpart for both the cases show that the dissimilar gas combination with large temperature difference shows higher growth rate. This result confirms earlier stability analysis that predicts increased growth rate for such cases. The growth rate reduction of a compressible mixing layer is argued to occur due to reduced pressure strain term in the Reynolds stress equation. This reduction also requires the pressure and density fluctuation correlation to be very near to unity. This holds good for a mixing layer formed between two similar gases at same temperature. For dissimilar gases at different temperatures this assumption does not hold well, and pressure-density correlation coefficient shows departure from unity. Further analysis of temperature density correlation factor, and temperature fluctuations shows that the changes in density occur predominantly due to temperature effects, than due to pressure effects. The mechanism of density variations is found to be different for similar and dissimilar gases, while for similar gases the density variations are due to pressure variations. For dissimilar gases density variation is also affected by temperature variations in addition to pressure variations. It has been observed that the traditional k-ε turbulence model within the RANS (Reynolds Averaged Navier Stokes) framework fails to capture the growth rate reduction for compressible shear layers. The performance of k-ε turbulence model is tested for the mixing of dissimilar gases at large temperature difference. For the experimental test case the shear layer growth rate and wall pressures show good match with other model free simulations. Simulations are further carried out for a range of convective Mach numbers keeping the mixing gases and their temperatures same. It has been observed that a drop in the growth rate is well predicted by RANS simulations. Further, the compressibility option has been removed and it has been observed that for the density and temperature difference, even for incompressible case, the drop in growth rate exists. This behaviour shows that the decrease in growth rate is mainly due to the interaction of temperature and species mass fraction on density. Also it can be inferred that RANS with k-ε turbulence model is able to capture the compressible shear layer growth rate for dissimilar gases at high temperature difference. The mixing of heat and species is governed by the values of turbulent Prandtl and Schmidt numbers respectively. These numbers have been observed to vary for different flow conditions, while affecting the flow field considerable in the form of temperature and species distribution. Model free simulations are carried out on an incompressible convective Mach number mixing layer, and the results are compared with that of a compressible mixing layer to study the effect of compressibility on the values of turbulent Prandtl / Schmidt numbers. It has been observed that both turbulent Prandtl and Schmidt numbers show an almost constant value in the mixing layer region for incompressible case. While, for a compressible case, both turbulent Prandtl and Schmidt numbers show a continuous variation within the mixing layer. However, the turbulent Lewis number is observed to be near unity for both incompressible and compressible cases. The thesis is composed of 8 chapters. An introduction of the subject with critical and relevant literature survey is presented in chapter 1. Chapter 2 describes the mathematical formulation and assumptions along with solution methodology needed for the simulations. Chapter 3 deals with the two and three dimensional model free simulations of the non reacting mixing layer. The effect of the presence of side wall is studied in chapter 4. Chapter 5 deals with the effect of compressibility on the mixing of two dissimilar gases at largely different temperatures. The performance of k-ε turbulence model is checked for dissimilar gases in Chapter 6. Chapter 7 is concerned with the effect of compressibility on turbulent Prandtl and Schmidt numbers. Finally concluding remarks are presented in chapter 8. The main aim of this thesis is the exploration of parallel mixing of dissimilar gases under compressible conditions for both two and three dimensional cases. The outcome of the thesis is (a) a finding that the presence of sidewall in a mixing duct does not make flow field two dimensional, instead it causes the formation of oblique disturbances and the shear layer growth rate is reduced, (b) that it has been shown that the growth rates of dissimilar gases are affected far more by large temperature difference than by compressibility as in case of similar gases, (c) that the growth rates of compressible shear layers formed between dissimilar gases are better predicted using k-εturbulence model and (d) that for compressible mixing conditions the turbulent Prandtl and Schmidt numbers vary continuously in the mixing layer region necessitating the use of some kind of model instead of assuming constant values.