Thermische Energiespeicher tragen u. a. zur Erhöhung der Versorgungsicherheit in der Fernwärmeversorgung und zur Effizienzsteigerung des Fernwärmesystems (z. B. Flexibilisierung der Erzeuger, Speicherung überschüssiger Wärme, besserer hydraulischer Betrieb) bei. Dafür eignen sich Druckbehälter, sog. schlanke Heißwasserspeicher (Speichertyp b1). Die oben genannten Vorteile setzen einen effizienten Speicherbetrieb (niedrige interne und externe Speicherverluste) voraus. Dieser Beitrag beschäftigt sich mit der Minimierung der internen Verluste durch die Verbesserung des thermischen Schichtungsverhaltens. Eine thermische Schichtung mit einem möglichst schmalen Übergangsbereich zwischen heißer und kalter Zone ist ein Indikator für geringe Mischvorgänge während der Beladung. Die Minimierung dieser Mischungsvorgänge bei der Beladung nimmt eine Schlüsselrolle bei der Minimierung der internen Speicherverluste ein. Lohse und Brähmer untersuchten die Beladung mit herkömmlichen radialen Diffusor in schlanken Heißwasserspeicher mit numerischer Strömungssimulation. Die Arbeiten identifizieren aufgrund der schlanken Speicherform nachteilige Strömungseffekte wie z. B. einen ausgeprägten Wandstrahl. Dieser Wandstrahl regt Mischvorgänge an und damit steigen die internen Speicherverluste. Zur Überwindung dieser Strömungsproblematik schlägt die Beladung mit Drall vor. Die Untersuchungen von Oestreich zeigten das Strömungsverhalten im Diffusor und im Speicher, die Auswirkungen auf die thermische Schichtung sowie die Vorteilhaftigkeit. Dieser Beitrag soll eine detailliertere Beschreibung der Strömungsvorgänge liefern. Dieses Wissen ist unbedingt notwendig, um die Ursachen und Wirkungen bei der Beladung mit Drall und beim Aufbau der thermischen Schichtung besser zu verstehen. Die Modellierung und Simulation des Diffusors bzw. des Speichers erfolgen mit ANSYS CFX. Zur Auflösung turbulenter Strukturen findet die Large Eddy Simulation Anwendung. Dieser Artikel präsentiert erstmalig die Wirbelstrukturen im Diffusor mit Leitelementen zur Drallerzeugung. Die Speicherströmung weist ein ähnliches Verhalten zu bekannten Dichteströmungen (z. B. Lappen-Kluft-Struktur, Instabilitäten in den freien Scherschichten) auf, was bisher nicht bekannt war. Hohe Peclet-Zahlen (hohe Advektionsströme) im Speichermodell führen zu numerischer Instabilität der Simulation und erfordern deshalb erhöhten Diskretisierungsaufwand. / Thermal energy storage systems contribute, among other things, to increasing the security of supply in the district heating system and to improving the efficiency of the district heating system (e.g., making the generators more flexible, storing waste heat, better hydraulic operation). Pressure vessels, so-called slim hot water storage tanks (storage type b1) are suitable for this purpose. The above mentioned advantages require efficient storage operation (low internal and external storage losses). This paper deals with the minimization of internal losses by improving the thermal stratification behavior. Thermal stratification with a thermocline between hot and cold zone as narrow as possible is an indicator of low mixing processes during loading. Minimizing these mixing processes during loading takes a key role in minimizing internal storage losses. Lohse and Brähmer investigated loading with conventional radial diffuser in slim hot water storage tanks with numerical flow simulation. The work identifies adverse flow effects due to the slim tank shape, such as a wall jet. This wall jet stimulates mixing processes and thus increases the internal storage losses. To overcome this flow problem, Findeisen et al. proposes swirl loading. The investigations of Oestreich et al. showed the flow behavior in the diffuser and in the storage, the effects on the thermal stratification as well as the advantageousness. This paper aims to provide a more detailed description of the flow processes. This knowledge is essential to better understanding the causes and effects of swirl loading and the structure of thermal stratification. Modeling and simulation of the diffuser and storage, respectively, are performed using Ansys CFX. Large eddy simulation (LES) is applied to resolve turbulent structures. This paper presents for the first time the vortex structures in the diffuser with internal elements for swirl generation. The storage flow exhibits similar behavior to known density flows (e.g., head and nose formation, instabilities in the free shear layers), which was previously unknown. High Peclet numbers (high advection currents) in the storage model lead to numerical instability of the simulation and therefore require increased discretization efforts.

No description available.

This thesis investigates the impact of upstream duct convolution on the plume development for high speed jets. In particular, investigations are carried out into an unconventional aero-engine exhaust systems comprised of a modified convergent-divergent rectangular nozzle where the converging section of the nozzle includes an S-bend in the duct. The motivation for this work comes from both the military and civilian sectors of the aerospace industry. The growing interest into highly efficient engines in the civilian sector and increasing complexities involved in stealth technologies for military applications has led to new design constraints on aero-engine exhaust systems that require further research into flows through more complex duct geometries. Due to a lack of experimental data into this area in the open literature validation studies are undertaken into flows through an S-bend duct and exhaust plume development from a rectangular convergent-divergent nozzle. The validation work is simulated using RANS CFD with common industrial turbulence models as well as LES with artificial inlet conditions. Subsequently, a CFD investigation into three unconventional aero-engine exhaust systems, with over-expanded conditions, with differing angles of curvature across the converging S-bend is undertaken using both RANS and LES methodologies governed by the validation work. As the curvature of the S-bend was increased it was found that the thrust and effective NPR both decrease. Whilst these changes were within acceptable levels (with some optimisation) for a circumferential extent of up to 53.1 the losses became prohibitive large at extents. For the ducts with a greater circumferential extents separation was seen to occur at the throat of the nozzle; this changes the design parameters of the nozzle leading to a higher Mach number and could potentially be harnessed to improve performance of the engine creating a `variable throat' nozzle. The impact of using different numerical solvers to simulate the flow through an unconventional aero-engine exhaust system has also been considered. The use of LES has shown that the octagonal, hexahedral and trapezoidal shapes initially observed in the development of the plumes of the RANS cases are likely to be an artifact caused by the RANS solver, as would the transverse total pressure gradients observed in the RANS cases at the nozzle exit as they are both absent from all of the LES results. Likewise the implementation of realistic inlet conditions has a significant impact on the development of the plume, particularly in the length of the potential core and the number of shock cells.

Large Eddy Simulation (LES) is carried out using the Rolls-Royce Hydra CFD code in order to investigate and give further insight into highly turbulent, unsteady flow structures for automotive applications. LES resolves time dependent eddies that are modelled in the steady-state by Reynolds-Averaged Navier-Stokes (RANS) turbulence models. A standard Smagorinsky subgrid scale model is used to model the energy transfer between large and subgrid scales. Since Hydra is an unstructured algorithm, a variety of unstructured hexahedral, tetrahedral and hybrid grids are used for the different cases investigated. Due to the computational requirements of LES, the cases in this study replicate and analyse generic flow problems through simplified geometry, rather than modelling accurate race car geometry which would lead to infeasible calculations. The first case investigates the flow around a diffuser-equipped bluff body at an experimental Reynolds number of 1.01 times 10 to the power 6 based on model height and inlet velocity. LES is carried out on unstructured hexahedral grids of 10 million and 20 million nodes, with the latter showing improved surface pressure when compared to the experiments. Comparisons of velocity and vorticity between the LES and experiments at the diffuser exit plane show a good level of agreement. Flow visualisation of the vortices in the diffuser region and behind the model from the mean and instantaneous flow attempts to explain the relation or otherwise between the two. The main weakness of the simulation was the late laminar to turbulent transition in the underbody region. The size of the domain and high experimental Reynolds number make this case very challenging. After the challenges faced by the diffuser-equipped bluff body, the underbody region is isolated so that increased grid refinement can be achieved in this region and the calculation is run at a Reynolds number of 220, 000, reducing the computational requirement from the previous case. A vortex generator mounted onto a flat underbody at an onset angle to the flow is modelled to generate vortices that extend along the length of the underbody and its interaction with the ground is analysed. Since the vortex generator resembles a slender wing with an incidence to the flow, a delta wing study is presented as a preliminary step since literature on automotive vortex generators in ground effect is scarce. Results from the delta wing study which is run at an experimental Reynolds number of 1.56 times 10 to the power 6 are in very good agreement with previous experiments and Detached Eddy Simulation (DES) studies, giving improved detail and understanding. Axial velocity and vorticity contours at several chordwise stations show that the leading edge vortices are predicted very well by a 20 million node tetrahedral grid. Sub-structures that originate from the leading edge of the wing and form around the core of the leading edge vortex are also captured. Large Eddy Simulation for the flow around an underbody vortex generator over a smooth ground and a rough ground is presented. A hexahedral grid of 40 million nodes is used for the smooth ground case, whilst a 48 million node hybrid grid was generated for the rough ground case so that the detailed geometry near the ground could be captured by tetrahedral cells. The geometry for the rough surface is modelled by scanning a tarmac surface to capture the cavities and protrusions in the ground. This is the first time that a rough surface representing a tarmac road has been computed in a CFD simulation, so that its effect on vortex decay can be studied. Flow visualisation of the instantaneous flow has shown strong interaction with the ground and the results from this study have given an initial understanding in this area.

In the last 10-15 years, large-eddy simulation (LES) has become well established for non-reacting flows, and several successful models have been developed for the transfer of momentum and kinetic energy to the subfilter-scales (SFS). However, for reacting flows, LES is still undergoing significant development. In particular, for many premixed combustion applications, the chemical reactions are confined to propagating surfaces that are significantly thinner than the computational grids used in practical LES. In these situations, the chemical kinetics and its interaction with the turbulence are not resolved and must be entirely modelled. There is, therefore, a need for accurate and robust physical modelling of combustion at the subfilter-scales. In this thesis, modelled transport equations for progress variable and flame surface density (FSD) were implemented and coupled to the Favre-filtered Navier-Stokes equations for a compressible reactive thermally perfect mixture. In order to reduce the computational costs and increase the resolution of simulating combusting flows, a parallel adaptive mesh (AMR) refinement finite-volume algorithm was extended and used for the prediction of turbulent premixed flames. The proposed LES methodology was applied to the numerical solution of freely propagating flames in decaying isotropic turbulent flow and Bunsen-type flames. Results for both stoichiometric and lean flames are presented. Comparisons are made between turbulent flame structure predictions for methane, propane, hydrogen fuels, and other available numerical results and experimental data. Details of subfilter-scale modelling, numerical solution scheme, computational results, and capabilities of the methodology for predicting premixed combustion processes are included in the discussions. The current study represents the first application of a full transport equation model for the FSD to LES of a laboratory-scale turbulent premixed flame. The comparisons of the LES results of this thesis to the experimental data provide strong support for the validity of the modelled transport equation for the FSD. While the LES predictions of turbulent burning rate are seemingly correct for flames lying within the wrinkled and corrugated flamelet regimes and for lower turbulence intensities, the findings cast doubt on the validity of the flamelet approximation for flames within the thin reaction zones regime.

This study investigates the complex flow structure generated by riparian emergent vegetation along the edge of floodplain. Detailed velocity and boundary shear stress measurements were carried out for various arrangements of emergent rigid cylindric rods of 3 mm, 6 mm and 9 mm diameters and for three different rod densities. In addition, the impact of foliage on the flow field was assessed during a series of experiments where brushes were used instead of smooth rods. The results of these new experiments are first presented. In addition to the laboratory data, field data was obtained through Acoustic Doppler Current Profiler measurements for two flood events in a stretch of the river Rhône that can be approximated to a straight compound channel with vegetated banks. The analysis of the flow structure highlights the presence of strong secondary circulation and increased vorticity on the river banks. The rods on the edge of the floodplain increase significantly flow resistance, reducing velocity and decreasing boundary shear stress. Flow rate was seen to decrease with increasing vegetative density for all cases except when foliage was added. This suggests that an optimum threshold density, for which a smaller density would lead to an increased flow rate might exist. Wakes trailing downstream of the vegetation stem, planform coherent structures advected between the main channel and the floodplain, and eddying motion in the flow due to enhanced turbulence anisotropy are among the defining patterns observed in the studied compound channel flows with one line of emergent vegetation along the edge of the floodplain. The Shiono and Knight Method (SKM) was modified in order to account for the increased turbulence activity due to the rods. The drag force term was introduced in the same way as in the work of Rameshwaran and Shiono (2007). However, a new term was added to the transverse shear stress term in the form of an Elder formulation, incorporating a friction drag coefficient which can be derived from the experimental data. In this proposed version, the advection term was set to zero. Another version of the SKM, similar to Rameshwaran and Shiono (2007), was also tested with the addition of a local drag friction only applied in the rod region. The proposed SKM version without the advection term was favored as it can be more closely related to the experimental data and to physical processes. Finally, the capabilities of Telemac-2D were tested against the experimental data for various turbulence models. The Large Eddy Simulation turbulence model highlighted some unsteady flow patterns that were observed during experiments, while satisfactorily predicting the lateral velocity and boundary shear stress distributions.

A large eddy simulation (LES) model has been developed and validated for turbulent non-premixed and partially premixed combustion systems. LES based combustion modelling strategy has the ability to capture the detailed structure of turbulent flames and account for the effects of radiation heat loss. Effects of radiation heat loss is modelled by employing an enthalpy-defect based non-adiabatic flamelet model (NAFM) in conjunction with a steady non-adiabatic flamelet approach. The steady laminar flamelet model (SLFM) is used with multiple flamelet solutions through the development of pre-integrated look up tables. The performance of the non-adiabatic model is assessed against experimental measurements of turbulent CH4/H2 bluff-body stabilized and swirl stabilized jet flames carried out by the University of Sydney combustion group. Significant enhancements in the predictions of mean thermal structure have been observed with both bluff body and swirl stabilized flames by the consideration of radiation heat loss through the non-adiabatic flamelet model. In particular, mass fractions of product species like CO2 and H2O have been improved with the consideration of radiation heat loss. From the Sydney University data the HM3e flame was also investigated with SLFM using multiple flamelet strategy and reasonably fair amount of success has been achieved. In this work, unsteady flamelet/progress variable (UFPV) approach based combustion model which has the potential to describe both non-premixed and partially premixed combustion, has been developed and incorporated in an in-house LES code. The probability density function (PDF) for reaction progress variable and scalar dissipation rate is assumed to follow a delta distribution while mixture fraction takes the shape of a beta PDF. The performance of the developed model in predicting the thermal structure of a partially premixed lifted turbulent jet flame in vitiated co-flow has been evaluated. The UFPV model has been found to successfully predict the flame lift-off, in contrast SLFM results in a false attached flame. The mean lift-off height is however over-predicted by UFPV-δ function model by ~20% for methane based flame and under-predicted by ~50% for hydrogen based flame. The form of the PDF for the reaction progress variable and inclusion of a scalar dissipation rate thus seems to have a strong influence on the predictions of gross characteristics of the flame. Inclusion of scalar dissipation rate in the calculations appears to be successful in predicting the flame extinction and re-ignition phenomena. The beta PDF distribution for the reaction progress variable would be a true prospect for extending the current simulation to predict the flame characteristics to a higher degree.

The periods of time-varying turbulence in the atmospheric boundary layer, i.e.\ the morning and evening transitions, are often overlooked or highly idealised by dispersion models. These transitions make up a significant portion of the diurnal cycle and are known to affect the spread of pollution due to the properties of turbulence in the residual and stable layers, resulting in phenomena such as lofting, trapping, and fumigation.\\ Two main simulation techniques are presented for the purpose of modelling the dispersion of passive tracers in both convective and evening transitional boundary layers: Lagrangian stochastic (LS) modelling for 1D, inhomogeneous, non-stationary turbulence; and large-eddy simulation (LES) with a particle model tracing pollutant paths using a combination of the resolved flow velocities and a random displacement model to represent sub-grid scale motions.\\ In the convective boundary layer, LS models more accurately representing the state of turbulence, and including the effect of skewness, are shown to produce dispersion results in closer agreement with LES. By considering individual particle trajectories, a reflective top boundary in LS models is shown to produce un-physical, sharp changes in velocity and position. By applying a correction to the vertical velocity variance based on representing the stable potential temperature gradient above the boundary layer, particles are contained within the boundary layer in a physically accurate way. \\ An LS model for predicting dispersion in time-varying, skewed turbulence is developed and tested for various particle releases in transitional boundary layers with different rates of decay, showing an improvement in accuracy compared with previous LS models. Further improvement is made by applying a correction to the vertical velocity variance to represent the effect of a positive potential temperature gradient developing over the course of the transition. Finally, a developing stable boundary layer is shown to have a significant trapping effect on particles released near the surface. \\

La recherche en bruit de combustion est de nos jours majoritairement consacrée au développement d'outils de calcul du bruit rayonné par les flammes. Les méthodes actuelles de CFD telles que la LES ou la DNS sont capables de fournir le champ acoustique rayonné par des sources de bruit, mais elles sont cependant limitées à des domaines de faible taille, ceci dû à leur fort coût de calcul. Pour surmonter cette limitation, on a vu l'émergence de méthodes hybrides. Dans cette approche, les sources de bruit sont découplées du son rayonné. Les sources sont alors calculées par les méthodes de DNS et de LES tandis que l'analogie acoustique permet de calculer le son rayonné par des codes acoustiques, moins coûteux en temps de calcul.L'objet de cette étude est de développer un outil numérique sur la base de l'analogie acoustique de Phillips pour de faibles nombres de Mach. Il permet de prendre en compte l'impact des conditions limites sur le champ acoustique résultant. La LES et le code de calcul acoustique développé ont été utilisés pour évaluer le bruit produit par une flamme turbulente confinée. Les deux techniques donnent des résultats en accord tant que les bonnes quantités sont comparées: il a été observé que le signal de pression obtenu directement à partir de la LES contient une quantité non négligeable d'hydrodynamique, laquelle doit être négligée si on cherche à comparer seulement les champs acoustiques issus des deux codes.L'hypothèse d'un nombre de Mach faible est totalement réaliste si l'on considère l'écoulement présent dans une chambre de combustion. Elle conduit à des simplifications significatives lorsque les analogies acoustiques sont considérées. Cependant, cette hypothèse ne peut pas être utilisée pour l'écoulement en amont (entrée d'air, compresseur) ni en aval (turbine, tuyère) des chambres de combustion aéronautiques. Un outil numérique a été développé pour pallier ce problème. Il est basé sur les équations d'Euler Quasi-1D, qui prennent en compte des écoulements convectifs, non isentropiques et non isenthalpiques. Cet outil permet d'estimer les conditions limites acoustiques qui doivent être imposées sur les entrées/sorties d'une chambre de combustion pour prendre en compte la présence d'un écoulement de nombre de Mach non négligeable, alors que les calculs acoustiques sont eux effectués sous cette hypothèse fortement restrictive. / Today, much of the current effort in combustion noise is the development of efficient numerical tools to calculate the noise radiated by flames. Although unsteady CFD methods such as LES or DNS can directly provide the acoustic field radiated by noise sources, this evaluation is limited to small domains due to high computational costs. Hybrid methods have been developed to overcome this limitation. In these schemes, the noise sources are decoupled from the radiated sound. The sources are still calculated by DNS or LES codes whereas the radiated sound is evaluated by acoustic codes using an acoustic analogy.In the present study, a numerical tool based on the Phllips' analogy for low Mach numbers flows has been developed. This tool accounts for the role of the boundary conditions in the resulting acoustic field. Both LES and the acoustic code developed here are used to assess the noise produced by a turbulent confined flame of a turbulent swirled-stabilized staged combustor. Good agreements are obtained between both techniques as long as the good quantities are compared: the pressure signal obtained directly from LES contains a non negligible amount of hydrodynamics that must be removed when a suitable acoustics-acoustics comparison is sought. The low Mach number assumption is completely realistic when considering the flow within a combustion chamber; it also conducts to considerable simplifications when leading with acoustic analogies. However, it cannot be used for the upstream (air-intake, compressors) and downstream (turbines, nozzle) of an aeronautical combustion chamber. A numerical tool is developed based on the quasi-1D Linearized Euler Equations in order to account for convective, non-isentropic and non-isenthalpic flows. By means of this tool, it is possible to estimate the acoustic boundary conditions that should be imposed at the inlet/oultlet of a given combustion chamber when performing low-Mach number acoustic computations.

Numerical simulation allows investigation into the influence of separation distance and rotation on the performance of two vertical-axis hydrokinetic turbines. Compu- tational fluid dynamics is applied to calulate the lift and drag coefficients acting upon interacting NACA 0021 turbine blades for a Reynolds number of Red = 10, 000. To understand the effect of separation distance, large-eddy simulation of the flow around side-by-side and staggered cylinders, ReD = 3,000, and airfoils, Rec = 3,000, are also performed. Based upon the simulations, a drag reduction of 11.3% and 19.8% is determined for the downstream cylinder and airfoil, respectively. A reduction in Reynolds stresses is also observed for the staggered configuration compared to the side-by-side configuration. Due to computational resources of large-eddy simulation, the Reynolds averaged Navier-Stokes method is also applied to investigate the influence of separation distance and rotation on two vertical axis hydrokinetic turbines. The numerical simulations show that a drag reduction of 15.5% occurs when the non-dimensional spanwise and streamwise separation distances, based on turbine diameter, reach 1 and 2, respectively. / October 2016