The numerical modelling of the VKI-Longshot facility remains a challeng-ing task as it requires multi-physical numerical methods in order to simulate all the components. In the current dissertation, numerical tools were developed in order to study each component of the facility separately and a deep investigations of each stage of the shot were performed. This helped to better understand the different processes involved in the flow development inside this hypersonic wind tunnel. However the numerical computation of different regions of the facility treated as independent from each others remains an approximation at best.The accuracy of the rebuilding code for determining the free stream conditions and the total enthalpy in the VKI-Longshot facility was investigated by using a series of unsteady numerical computations of axisymmetric hypersonic flow over a heat flux probe. Good agreement was obtained between the numerical results and the measured data for both the stagnation pressure and the heat flux dur- ing the useful test time.The driver-driven part of the Longshot facility was modelled using the quasi one-dimensional Lagrangian solver L1d2. The three main conditions used for the experiments —low, medium and high Reynolds number —were considered.The chambrage effect due to the junction between the driver and the driven tubes in the VKI-Longshot facility was investigated. The computation showed great ben- efit of the chambrage in increasing the speed of the piston and thus the final compression ratio of the test gas.Two dimensional simulations of the flow in the driver and the driven tube were performed using Arbitrary Lagrangian Eulerian (ALE) solver in COOLFLuiD. A parallel multi-domain strategy was developed in order to integrate the moving piston within the computational domain.The computed pressure in the reservoir is compared to the one provided by the experiment and good agreement was obtained for both con- editions.Finally, an attempt was made to compute the starting process of the flow in the contoured nozzle. The transient computation of the flow showed how the primary shock initiates the flow in the nozzle before reaching the exit plan at time of 1.5 [ms] after the diaphragm rupture. The complex interactions of the reflected shocks in the throat raise the temperature above 9500 [K] which was not expected. Chemical dissociation of Nitrogen was not taken into account during this transient investigation which may play a key role considering the range of temperature reached near the throat. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Le travail réalisé dans cette thèse introduit le concept de l'utilisation de la méthode Level-Set pour simuler l'interface Glace/Air au cours du temps lors du processus de givrage en vol des aéronefs. Pour cela, un code de givrage tri-dimensionnel multi-blocs et parallélisé a été implémenté au sein du solveur NSMB (Navier-Stokes-Multi-Blocks). Il comprend notamment un module de calcul des trajectoires des gouttelettes par une approche Eulérienne compatible avec l'utilisation de grilles chimères et un module thermodynamique pour le calcul des masses de glace incluant deux modèles différents : un modèle algébrique itératif et un modèle à dérivées partielles. Une attention particulière a été portée sur la vérification du code de givrage implémenté en comparant systématiquement, si possible, les résultats obtenus avec les données expérimentales et numériques existantes dans la littérature. Pour cette raison, le module de déformation de maillage existant dans NSMB a été intégré au code implémenté afin de pouvoir simuler le givrage par une méthode traditionnelle. Enfin, un nouveau principe pour le suivi de l'interface glace/air est introduit via l'utilisation d'une méthode Level-Set. Puisque dans ce travail de thèse nous nous intéressons particulièrement au concept, la méthode Level-Set développée est d'ordre un et est résolue implicitement. On montrera cependant que des résultats valides sont obtenus avec une telle approximation. / This thesis introduces the concept of the Level-Set method for simulating the evolution through time of the ice/air interface during the process of in-flight aircraft icing. For that purpose, a three-dimensionnal multi-block and parallelized icing code have been implemented in the NSMB flow solver (Navier-Stokes-Multi-Blocks). It includes a module for calculating the droplet trajectories by an Eulerian approach compatible with the use of chimera grids and a thermodynamic module to calculate the ice masses including two different models : an iterative algebraic model and a PDE model. Particular attention was paid to the validation of the icing code irnplemented by comparing results with existing experimental and numerical data in the literature. For this reason, the existing mesh deformation algorithm in NSMB was integrated into the code to simulate icing by a traditional method. Finally a new principle to track the ice/air interface is introduced using the Level-Set method. Since we are particularly interested in the concept, the Level-Set method developped is first order and solved implicitly. However it will be shown that valid results are obtained with such an approximation.

Au sein des opérations de retraitement du combustible usé, la colonne pulsée à garnissage est l'appareil d'extraction liquide-liquide principalement utilisé. Dans un contexte de compétitivité économique et de raréfaction des ressources, l'efficacité de ces appareils est devenue un enjeu pour l'industriel. Afin d'améliorer leur rendement à travers un meilleur dimensionnement, la connaissance de la vitesse de glissement entre les phases de l'émulsion est nécessaire.Les travaux menés et présentés dans ce manuscrit s'articulent autour de la modélisation physique et numérique de l'hydrodynamique de l'émulsion ainsi que de sa caractérisation expérimentale.Dans ce travail, une modélisation d'approche eulérienne, inspirée des travaux de D. Lhuillier, permet l'obtention d'un modèle bi-fluide couplé à une équation d'évolution de la surface d'échange (aire interfaciale). La résolution du modèle s'effectue par éléments finis sous le logiciel CAST3M. A l'issue des calculs, le modèle montre sa capacité à restituer le comportement de l'émulsion et permet l'obtention des vitesses de glissement. Dans une optique de validation expérimentale du modèle, des expérimentations sur deux installations sont menées. Celles-ci font notamment intervenir un couplage entre les méthodes de vélocimétrie par image de particules et de fluorescence induite par laser afin d'obtenir les vitesses de chaque phases et le taux de rétention de la phase dispersée. Un algorithme de détection et de suivi de gouttes est développé afin d'obtenir la vitesse de la phase dispersée et sa fraction volumique. La confrontation de ces aux résultats numériques permet une première qualification encourageante du modèle. / In the core of spent fuel reprocessing operations, the pulsed columns with packing are the liquid-liquid extraction apparati mainly used. The context of economical competiveness and scarce resources, industrials are driven to improve the efficiency of these processes. Pulsed column efficiency is bound to the amount of available exchange surface, which depends on geometrical parameters of the column and the operating conditions. A better design would improve the efficiency. In this aim the knowledge of the interphase slip velocity is necessary. The work presented in this thesis revolves around physical and numerical modelling of the hydrodynamics of the emulsion and its experimental characterization.In this work, a eulerian approach, based on the work of D.Lhuillier, allows to obtain a two-fluid model coupled with an evolution equation of the exchange surface (interfacial area). We use finite elements method to solve this model along with CAST3M software. Numerical simulations have shown the model abilities to correctly reproduce the emulsion behaviour and to obtain the slip velocity.In order to experimentally validate the model, we carried out two types of experimentation. Particles images velocimetry coupled to laser induced fluorescence are involved to obtain velocities of each phases and the dispersed phase volume fraction. We developed a tracking algorithm to obtain the dispersed phase velocity and the hold up. These results, such as velocities and strain rate tensor, have been used in a first validation of the model.

Les écoulements diphasiques turbulents sont présents dans de nombreux systèmes industriels (moteur à piston, turbines à gaz, moteurs fusée...). La compréhension fine de telles configurations s'avèrent de nos jours nécessaire pour limiter notamment les émissions de polluants et de gaz à effet de serre, et la consommation des énergies fossiles. Nous nous intéressons ici à la simulation aux grandes échelles des écoulements diphasiques turbulents, permettant de capturer une large partie du spectre de la turbulence, et ainsi être capable de prédire des phénomènes instables ou transitoires. La phase dispersée est ici modélisée par une approche eulérienne, en raison de ses avantages dans le contexte du calcul haute performance. Le travail de cette thèse a consisté à étendre le formalisme eulérien existant dans le code AVBP à la simulation de sprays polydisperses dans des écoulements turbulents. Pour cela, le Formalisme Eulérien Mésoscopique (FEM) a été couplé à une approche Multi-fluide. Cette nouvelle approche, intitulée Formalisme Eulérien Mésoscopique Multi-fluide (FEMM), a été évaluée sur des cas simples canoniques, permettant de bien caractériser le comportement autant en terme de dynamique turbulente que d'effets polydisperses. Les stratégies numériques disponibles dans le code de calcul AVBP sont aussi analysées, afin d'en cerner les limites pour la simulation eulérienne d'une phase liquide. Ce nouveau formalisme est finalement appliqué à la configuration aéronautique MERCATO, pour laquelle on dispose de résultats numériques obtenus avec d'autres approches (FEM et approche lagrangienne), et de résultats expérimentaux. Un accord satisfaisant avec l'expérience est montré pour toutes les approches, même si le FEM, monodisperse, obtient de moins bon résultats en terme de fluctuations. D'autres résultats expérimentaux s'avèrent nécessaires pour évaluer les approches et déterminer quelle est la plus prédictive pour cette configuration, notamment concernant la fraction massique de kerosene, autant en phase liquide qu'en phase gazeuse. / Turbulent two-phase flows are encountered in several industrial devices (piston engine, gas turbine, rocket engine...). A fine understanding of such configurations is mandatory to face problems of pollutant emissions, greenhouse gas, and fossil fuel rarefaction. The Large Eddy Simulation seems to be a good candidate. This kind of simulation captures a wide part of turbulence spectrum, and thus allows to predict instabilities and transient phenomena. The dispersed phase is simulated using an Eulerian approach, which seems to be more suitable than lagrangian methods for High Performance Computing. The present work consists in the extension to polydisperse flows of the existing eulerian formalism in the AVBP code. The Mesoscopic Eulerian Formalism (MEF) is coupled with the Multifluid approach. This new formalism, called Multifluid Mesoscopic Eulerian Formalism, is evaluated on simple test cases, showing the ability of such approach to capture turbulent and polydisperse effects. Numerical strategies available in AVBP are also evaluated, in order to emphasize on their limiting aspects for the eulerian simulation of a dispersed phase. The new formalism is finally applied to the simulation of the aeronautical configuration called MERCATO. Several experimental results are available, as well as numerical results using FEM and lagrangian approach. Results show a good agreement between experiments and numerical results, even if FEM results are worse concerning the fluctuations. New experimental results are necessary to determine which is the best approach, especially in terms of liquid and gas kerosene mass fraction.

This thesis presents the development of a numerical water quality model using a general paradigm of reaction-based approaches. In a reaction-based approach, all conceptualized biogeochemical processes are transformed into a reaction network. Through the decomposition of species governing equations via Gauss-Jordan column reduction of the reaction network, (1) redundant fast reactions and irrelevant kinetic reactions are removed from the system, which alleviates the problem of unnecessary and erroneous formulation and parameterization of these reactions, and (2) fast reactions and slow reactions are decoupled, which enables robust numerical integrations. The system of species transport equations is transformed to reaction-extent transport equations, which is then approximated with two subsets: algebraic equations and kinetic-variables transport equations. As a result, the model alleviates the needs of using simple partitions for fast reactions. With the diagonalization strategy, it makes the inclusion of arbitrary number of fast and kinetic reactions relatively easy, and, more importantly, it enables the formulation and parameterization of kinetic reactions one by one. To demonstrate the general paradigm, QAUL2E was recasted in the mode of a reaction network. The model then was applied to the Loxahatchee estuary to study its response to a hypothetical biogeochemical loading from its surrounding drainage. Preliminary results indicated that the model can simulate four interacting biogeochemical processes: algae kinetics, nitrogen cycle, phosphorus cycle, and dissolved oxygen balance.

<p dir="ltr">Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, volume 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the pre-chamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within 10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (p_r^o) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above p_r^o= 2.2. Case with p_r^o= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.</p>

No description available.

This dissertation explores an analogue of the prime number theorem for polynomials over finite fields as well as its connection to the necklace factorization algorithm T-transform and the string complexity measure T-complexity. Specifically, a precise asymptotic expansion for the prime polynomial counting function is derived. The approximation given is more accurate than previous results in the literature while requiring very little computational effort. In this context asymptotic series expansions for Lerch transcendent, Eulerian polynomials, truncated polylogarithm, and polylogarithms of negative integer order are also provided. The expansion formulas developed are general and have applications in numerous areas other than the enumeration of prime polynomials. A bijection between the equivalence classes of aperiodic necklaces and monic prime polynomials is utilized to derive an asymptotic bound on the maximal T-complexity value of a string. Furthermore, the statistical behaviour of uniform random sequences that are factored via the T-transform are investigated, and an accurate probabilistic model for short necklace factors is presented. Finally, a T-complexity based conditional string complexity measure is proposed and used to define the normalized T-complexity distance that measures similarity between strings. The T-complexity distance is proven to not be a metric. However, the measure can be computed in linear time and space making it a suitable choice for large data sets. / Graduate / 0544 0984 0405 / nrebenich@gmail.com

In power plants using large utility coal-fired boilers for generation of electricity the coal is pulverised in coal mills and then it has to be pneumatically transported and distributed to a larger number of burners (e.g. 30-40) circumferentially arranged in several rows around the burning chamber of the boiler. Besides the large pipework flow splitting devices are necessary for distribution of an equal amount of pulverised fuel (PF) to each of the burners. So called trifurcators (without inner fittings or guiding vanes) and ''riffle'' type bifurcators are commonly used to split the gas-coal particle flow into two or three pipes/channels with an equal amount of PF mass flow rate in each outflow cross section of the flow splitting device. These PF flow splitting devices are subject of a number of problems. First of all an uneven distribution of PF over the burners of a large utility boiler leads to operational and maintenance problems, increased level of unburned carbon and higher rates of NOX emissions. Maldistribution of fuel between burners caused by non uniform concentration of the PF (particle roping) in pipe and channel bends prior to flow splitting devices leads to uncontrolled differences in the fuel to air ratio between burners. This results in localised regions in the furnace which are fuel rich, where insufficient air causes incomplete combustion of the fuel. Other regions in the furnace become fuel lean, forming high local concentrations of NOX due to the high local concentrations of O2. Otherwise PF maldistribution can impact on power plant maintenance in terms of uneven wear on PF pipework, flow splitters as well as the effects on boiler panels (PF deposition, corrosion, slagging). In order to address these problems in establishing uniform PF distribution over the outlet cross sections of flow splitting devices in the pipework of coal-fired power plants the present paper deals with numerical prediction and analysis of the complex gas and coal particle (PF) flow through trifurcators and ''riffle'' type bifurcators. The numerical investigation is based on a 3-dimensional Eulerian- Lagrangian approach (MISTRAL/PartFlow-3D) developed by Frank et al. The numerical method is capable to predict isothermal, incompressible, steady gas- particle flows in 3-dimensional, geometrically complex flow geometries using boundary fitted, block-structured, numerical grids. Due to the very high numerical effort of the investigated gas-particle flows the numerical approach has been developed with special emphasis on efficient parallel computing on clusters of workstations or other high performance computing architectures. Besides the aerodynamically interaction between the carrier fluid phase and the PF particles the gas-particle flow is mainly influenced by particle-wall interactions with the outer wall boundaries and the inner fittings and guiding vanes of the investigated flow splitting devices. In order to allow accurate quantitative prediction of the motion of the disperse phase the numerical model requires detailed information about the particle-wall collision process. In commonly used physical models of the particle-wall interaction this is the knowledge or experimental prediction of the restitution coefficients (dynamic friction coefficient, coefficient of restitution) for the used combination of particle and wall material, e.g. PF particles on steel. In the present investigation these parameters of the particle-wall interaction model have been obtained from special experiments in two test facilities. Basic experiments to clarify the details of the particle-wall interaction process were made in a test facility with a spherical disk accelerator. This test facility furthermore provides the opportunity to investigate the bouncing process under normal pressure as well as under vacuum conditions, thus excluding aerodynamically influences on the motion of small particles in the near vicinity of solid wall surfaces (especially under small angles of attack). In this experiments spherical glass beads were used as particle material. In a second test facility we have investigated the real impact of non-spherical pulverised fuel particles on a steel/ceramic target. In this experiments PF particles were accelerated by an injector using inert gas like e.g. CO2 or N2 as the carrier phase in order to avoid dust explosion hazards. The obtained data for the particle-wall collision models were compared to those obtained for glass spheres, where bouncing models are proofed to be valid. Furthermore the second test facility was used to obtain particle erosion rates for PF particles on steel targets as a function of impact angles and velocities. The results of experimental investigations has been incorporated into the numerical model. Hereafter the numerical approach MISTRAL/PartFlow-3D has been applied to the PF flow through a ''riffle'' type bifurcator. Using ICEM/CFD-Hexa as grid generator a numerical mesh with approximately 4 million grid cells has been designed for approximation of the complex geometry of the flow splitting device with all its interior fittings and guiding vanes. Based on a predicted gas flow field a large number of PF particles are tracked throughout the flow geometry of the flow-splitter. Besides mean quantities of the particle flow field like e.g. local particle concentrations, mean particle velocities, distribution of mean particle diameter, etc. it is now possible to obtain information about particle erosion on riffle plates and guiding vanes of the flow splitting device. Furthermore the influence of different roping patterns in front of the flow splitter on the uniformness of PF mass flow rate splitting after the bifurcator has been investigated numerically. Results show the efficient operation of the investigated bifurcator in absence of particle roping, this means under conditions of an uniform PF particle concentration distribution in the inflow cross section of the bifurcator. If particle roping occurs and particle concentration differs over the pipe cross section in front of the bifurcator the equal PF particle mass flow rate splitting can be strongly deteriorated in dependence on the location and intensity of the particle rope or particle concentration irregularities. The presented results show the importance of further development of efficient rope splitting devices for applications in coal-fired power plants. Numerical analysis can be used as an efficient tool for their investigation and further optimisation under various operating and flow conditions.

The transport or the separation of solid particles or droplets suspended in a fluid flow is a common task in mechanical and process engineering. To improve machinery and physical processes (e.g. for coal combustion, reduction of NO_x and soot) an optimization of complex phenomena by simulation applying the fundamental conservation equations is required. Fluid-particle flows are characterized by the ratio of density of the two phases gamma=rho_P/rho_F, by the Stokes number St=tau_P/tau_F and by the loading in terms of void and mass fraction. Those numbers (Stokes number, gamma) define the flow regime and which relevant forces are acting on the particle. Dependent on the geometrical configuration the particle-wall interaction might have a heavy impact on the mean flow structure. The occurrence of particle-particle collisions becomes also more and more important with the increase of the local void fraction of the particulate phase. With increase of the particle loading the interaction with the fluid phase can not been neglected and 2-way or even 4-way coupling between the continous and disperse phases has to be taken into account. For dilute to moderate dense particle flows the Euler-Lagrange method is capable to resolve the main flow mechanism. An accurate computation needs unfortunately a high number of numerical particles (1,...,10^7) to get the reliable statistics for the underlying modelling correlations. Due to the fact that a Lagrangian algorithm cannot be vectorized for complex meshes the only way to finish those simulations in a reasonable time is the parallization applying the message passing paradigma. Frank et al. describes the basic ideas for a parallel Eulererian-Lagrangian solver, which uses multigrid for acceleration of the flow equations. The performance figures are quite good, though only steady problems are tackled. The presented paper is aimed to the numerical prediction of time-dependend fluid-particle flows using the simultanous particle tracking approach based on the Eulerian-Lagrangian and the particle-source-in-cell (PSI-Cell) approach. It is shown in the paper that for the unsteady flow prediction efficiency and load balancing of the parallel numerical simulation is an even more pronounced problem in comparison with the steady flow calculations, because the time steps for the time integration along one particle trajectory are very small per one time step of fluid flow integration and so the floating point workload on a single processor node is usualy rather low. Much time is spent for communication and waiting time of the processors, because for cold flow particle convection not very extensive calculations are necessary. One remedy might be a highspeed switch like Myrinet or Dolphin PCI/SCI (500 MByte/s), which could balance the relative high floating point performance of INTEL PIII processors and the weak capacity of the Fast-Ethernet communication network (100 Mbit/s) of the Chemnitz Linux Cluster (CLIC) used for the presented calculations. Corresponding to the discussed examples calculation times and parallel performance will be presented. Another point is the communication of many small packages, which should be summed up to bigger messages, because each message requires a startup time independently of its size. Summarising the potential of such a parallel algorithm, it will be shown that a Beowulf-type cluster computer is a highly competitve alternative to the classical main frame computer for the investigated Eulerian-Lagrangian simultanous particle tracking approach.