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.

Measuring the liquid velocity and turbulence parameters in multiphase flows is a challenging task. In general, measurements based on optical methods are hindered by the presence of the gas phase. In the present work, it is shown that this leads to a sampling bias. Here, particle image velocimetry (PIV) is used to measure the liquid velocity and turbulence in a bubble column for different gas volume flow rates. As a result, passing bubbles lead to a significant sampling bias, which is evaluated by the mean liquid velocity and Reynolds stress tensor components. To overcome the sampling bias a window averaging procedure that waits a time depending on the locally distributed velocity information (hold processor) is derived. The procedure is demonstrated for an analytical test function. The PIV results obtained with the hold processor are reasonable for all values. By using the new procedure, reliable liquid velocity measurements in bubbly flows, which are vitally needed for CFD validation and modeling, are possible. In addition, the findings are general and can be applied to other flow situations and measuring techniques.

Reservoir testing and analysis are fundamental tools in understanding reservoir hydraulics and hence forecasting reservoir responses. The quality of the analysis is very dependent on the conceptual model used in investigating the responses under different flowing conditions. The use of reservoir testing in the characterization and derivation of reservoir parameters is widely established, especially in conventional oil and gas reservoirs. However, with depleting conventional reserves, the quest for unconventional reservoirs to secure the increasing demand for energy is increasing; which has triggered intensive research in the fields of reservoir characterization. Gas hydrate reservoirs, being one of the unconventional gas reservoirs with huge energy potential, is still in the juvenile stage with reservoir testing as compared to the other unconventional reservoirs. The endothermic dissociation hydrates to gas and water requires addressing multiphase flow and heat energy balance, which has made efforts to develop reservoir testing models in this field difficult. As of now, analytically quantifying the effect on hydrate dissociation on rate and pressure transient responses are till date a huge challenge. During depressurization, the heat energy stored in the reservoir is used up and due to the endothermic nature of the dissociation; heat flux begins from the confining layers. For Class 3 gas hydrates, just heat conduction would be responsible for the heat influx and further hydrate dissociation; however, the moving boundary problem could also be an issue to address in this reservoir, depending on the equilibrium pressure. To address heat flux problem, a proper definition of the inner boundary condition for temperature propagation using a Clausius-Clapeyron type hydrate equilibrium model is required. In Class 1 and 2, crossflow problems would occur and depending on the layer of production, convective heat influx from the free fluid layer and heat conduction from the cap rock of the hydrate layer would be further issues to address. All these phenomena make the derivation of a suitable reservoir testing model very complex. However, with a strong combination of heat energy and mass balance techniques, a representative diffusivity equation can be derived. Reservoir testing models have been developed and responses investigated for different boundary conditions in normally pressured Class 3 gas hydrates, over-pressured Class 3 gas hydrates (moving boundary problem) and Class 1 and 2 gas hydrates (crossflow problem). The effects of heat flux on the reservoir responses have been addressed in detail.

The legislation has forced the vehicle industry to reduce tail-end emissions. The air pollutant nitrogen oxide (NOX) has been shown to have a negative impact on human health and the environment. One of the key technologies to reduce the levels of NOX emitted from a vehicle is by implementing an after-treatment system. The after-treatment system includes catalysts, a particle filter and an evaporation system. In the evaporation system a liquid jet containing a urea-water solution known as AdBlue is injected into the hot exhaust gases to evaporate into gaseous ammonia NH3 and water H2O. Then NH3 enters the Selective Catalytic Reduction (SCR) catalyst where it chemically reacts with NOX to form N2 and H2O. Problems can arise if an excessive amount of AdBlue is injected and a fluid film is formed on evaporation surfaces. At certain operating conditions the fluid film can crystallise and form solid deposits. The solid deposits can cause high back-pressure, material deterioration and ammonia slip. This project is done in collaboration with Scania CV AB. Scania is a world-leading manufacturer of heavy-duty vehicles, busses and engines. Scania works continuously to develop new simulation methods to capture the complex phenomena of AdBlue spray, wall film dynamics and risk of solid deposits, to use in the development process of new components. The aim of this project is to implement and evaluate a new method to predict the risk of crystallisation of urea (AdBlue) using the software Star-CCM+. Two different geometries are studied, a test rig and a Scania silencer. Different operating conditions, parameter settings and a speed-up method are analysed. During the project a base-line model has been created and the results have been compared with measurement results and the software AVL Fire. The results on the test rig show the effect of altering the mesh and important model parameters. Injected particles are grouped into parcels with the same properties. The number of parcels is a crucial factor for the wall film formation and should be sufficiently high to get a statistical representation of the droplet size distribution. The results from the real silencer show strong evaporation and thin wall film formation with the suggested method. The method is shown to be stable and the software is user-friendly. A speed-up method was investigated to decrease the computational time. The computational time was reduced by a factor 20. The outcome of this project is a guide for set-up of AdBlue spray and wall film simulations. Recommendations to future work includes further validation of the settings, investigation of the evaporation rate and droplet size distribution and the application to other cases. The next step is also to tune the critical thresholds for deposit risk assessment. / Lagstiftning har tvingat fordonstillverkare att minska avgasutsläppen. Luftföroreningen kväveoxid (NOX) har visat sig ha en negativ inverkan på människors hälsa och på miljön. En viktig teknik för att minska utsläppen av NOX ¨ar att implementera ett efterbehandlingssystem. Efterbehandlingssystemet tar hand om avgaserna genom substrat, filter och ett förångningssystem. I förångningssystemet sprutas en urea-vattenlösning, som kallas AdBlue, in i de heta avgaserna där den förångas till ammoniak NH3 och vatten H2O. Ammoniakgasen leds därefter in till SCR katalysatorn där den kemiskt reagerar med NOX och bildar kvävgas N2 och vattenånga. Problem kan uppstå om fel mängd AdBlue sprutas in, då kan vätska byggas upp på förångsningsytor, kristallisera och bilda avlagringar. Avlagringarna kan bygga upp en solid klump som kan orsaka ett högt mottryck, nedbrytning av material och ammoniakslip. Detta arbete är ett samarbete med Scania CV AB som är en världsledande producent av lastbilar, bussar och motorer. Scania arbetar kontinuerligt med att utveckla nya simuleringsvertyg för att beskriva uppkomsten av Urea avlagringar för att använda i utvecklingen av nya komponenter. Syftet med detta arbete är att implementera och utvärdera en ny metod för att prediktera klump mha simuleringsverktyget Star-CCM+. Två olika geometrier är använd i arbetet: en testrigg och en av Scanias ljuddämpare. Olika driftspunkter, parametrar och en uppsnabbad metod är studerade. Under projektets gång har en modell byggts upp och jämförts med mätningar och simuleringar från programvaran AVL Fire. Resultatet från simuleringarna på testriggen visar effekten av att variera olika parametrar. Partiklarna som sprutas in i systemet är grupperade i paket med liknande egenskaper. Antalet paket påverkar uppbyggnaden av väggfilm och det rekommenderas att denna parameter hålls hög för att statistiskt beskriva droppfördelningen av partiklar. Resultaten på ljuddämparen visar en stark förångning och en tunn väggfilm för samtliga driftspunkter. Den implementerade metoden har visat sig vara stabil och användarvänlig. En uppsnabbad metod har utvärderats för att minska beräkningstiden. Beräkningstiden kunde minskas med en faktor 20. Resultatet av arbetet är en guide för hur metoden implementeras och bör användas. Rekommendationer till framtida arbete är en fortsatt undersökning av parametrar, utvärdering av förångningsmodellen, validering av droppstorleksfördelningen och tillämpningen på andra geometrier. Nästa steg i utvecklingen skulle vara att kalibrera tröskelvärden för prediktering av klump.

This dissertation is a contribution to the development of a unified model of continuum mechanics, describing both fluids and elastic solids as a general continua, with a simple material parameter choice being the distinction between inviscid or viscous fluid, or elastic solids or visco-elasto-plastic media. Additional physical effects such as surface tension, rate-dependent material failure and fatigue can be, and have been, included in the same formalism. The model extends a hyperelastic formulation of solid mechanics in Eulerian coordinates to fluid flows by means of stiff algebraic relaxation source terms. The governing equations are then solved by means of high order ADER Discontinuous Galerkin and Finite Volume schemes on fixed Cartesian meshes and on moving unstructured polygonal meshes with adaptive connectivity, the latter constructed and moved by means of a in- house Fortran library for the generation of high quality Delaunay and Voronoi meshes. Further, the thesis introduces a new family of exponential-type and semi- analytical time-integration methods for the stiff source terms governing friction and pressure relaxation in Baer-Nunziato compressible multiphase flows, as well as for relaxation in the unified model of continuum mechanics, associated with viscosity and plasticity, and heat conduction effects. Theoretical consideration about the model are also given, from the solution of weak hyperbolicity issues affecting some special cases of the governing equations, to the computation of accurate eigenvalue estimates, to the discussion of the geometrical structure of the equations and involution constraints of curl type, then enforced both via a GLM curl cleaning method, and by means of special involution-preserving discrete differential operators, implemented in a semi-implicit framework. Concerning applications to real-world problems, this thesis includes simulation ranging from low-Mach viscous two-phase flow, to shockwaves in compressible viscous flow on unstructured moving grids, to diffuse interface crack formation in solids.