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.

Part 1: Theoretic background information - Review of Hooke’s law for linear elastic materials - The strain energy density of linear elastic materials - Hyperelastic material - Material laws for hyperelastic materials - About selecting the material model and performing tests - Implementation of hyperelastic material laws in Mechanica - Defining hyperelastic material parameters in Mechanica - Test set-ups and specimen shapes of the supported material tests - The uniaxial compression test - Stress and strain definitions in the Mechanica LDA analysis Part 2: Application examples - A test specimen subjected to uniaxial loading - A volumetric compression test - A planar test - Influence of the material law Appendix - PTC Simulation Services Introduction - Dictionary Technical English-German / Teil 1: Theoretische Hintergrundinformation - Das Hookesche Gesetz für linear-elastische Werkstoffe - Die Dehnungsenergiedichte für linear-elastische Materialien - Hyperelastisches Material - Materialgesetze für Hyperelastizität - Auswählen des Materialgesetzes und Testdurchführung - Implementierung der hyperelastischen Materialgesetze in Mechanica - Definieren der hyperelastischen Materialparameter in Mechanica - Testaufbauten und Prüfkörper der unterstützten Materialtests - Der einachsige Druckversuch - Spannungs- und Dehnungsdefinition in der Mechanica-Analyse mit großen Verformungen Teil 2: Anwendungsbeispiele - Ein einachsig beanspruchter Prüfkörper - Ein volumetrischer Drucktest - Ein planarer Test - Einfluss des Materialgesetzes Anhang: - Kurzvorstellung der PTC Simulationsdienstleistungen - Wörterbuch technisches Englisch-Deutsch

Part 1: Theoretic background information - Review of Hooke’s law for linear elastic materials - The strain energy density of linear elastic materials - Hyperelastic material - Material laws for hyperelastic materials - About selecting the material model and performing tests - Implementation of hyperelastic material laws in Mechanica - Defining hyperelastic material parameters in Mechanica - Test set-ups and specimen shapes of the supported material tests - The uniaxial compression test - Stress and strain definitions in the Mechanica LDA analysis Part 2: Application examples - A test specimen subjected to uniaxial loading - A volumetric compression test - A planar test - Influence of the material law Appendix - PTC Simulation Services Introduction - Dictionary Technical English-German / Teil 1: Theoretische Hintergrundinformation - Das Hookesche Gesetz für linear-elastische Werkstoffe - Die Dehnungsenergiedichte für linear-elastische Materialien - Hyperelastisches Material - Materialgesetze für Hyperelastizität - Auswählen des Materialgesetzes und Testdurchführung - Implementierung der hyperelastischen Materialgesetze in Mechanica - Definieren der hyperelastischen Materialparameter in Mechanica - Testaufbauten und Prüfkörper der unterstützten Materialtests - Der einachsige Druckversuch - Spannungs- und Dehnungsdefinition in der Mechanica-Analyse mit großen Verformungen Teil 2: Anwendungsbeispiele - Ein einachsig beanspruchter Prüfkörper - Ein volumetrischer Drucktest - Ein planarer Test - Einfluss des Materialgesetzes Anhang: - Kurzvorstellung der PTC Simulationsdienstleistungen - Wörterbuch technisches Englisch-Deutsch

We present a variational method for unfolding of the cortex based on a user-chosen point of view as an alternative to more traditional global flattening methods, which incur more distortion around the region of interest. Our approach involves three novel contributions. The first is an energy function and its corresponding gradient flow to measure the average visibility of a region of interest of a surface from a given viewpoint. The second is an additional energy function and flow designed to preserve the 3D topology of the evolving surface. This latter contribution receives significant focus in this thesis as it is crucial to obtain the desired unfolding effect derived from the first energy functional and flow. Without it, the resulting topology changes render the unconstrained evolution uninteresting for the purpose of cortical visualization, exploration, and inspection. The third is a method that dramatically improves the computational speed of the 3D topology-preservation approach by creating a tree structure of the triangulated surface and using a recursion technique.