Refine
Document Type
- Doctoral thesis (10)
Has Fulltext
- yes (10)
Is part of the Bibliography
- no (10)
Keywords
- Turbulente Strömung (10) (remove)
Numerical investigation and extension of quadrature-based moment methods for population balances
(2023)
Particulate systems can be described by a number density function (NDF) with respect to a vector of internal coordinates. The evolution of the NDF is governed by the typically high-dimensional population balance equation (PBE). A common approach to reduce the dimensionality of the problem is to solve only for a set of moments instead of the NDF. The derived system of moment equations, however, includes unclosed integral terms that still contain the unknown NDF. One way to close the system of moment equations is to approximate the unclosed integral terms using a Gaussian quadrature computed from the moments. The procedure of taking a set of moments to compute a Gaussian quadrature, which is, in turn, used to close the moment equations, is known as the quadrature method of moments (QMOM). It gave rise to an entire family of methods, the quadrature-based moment methods (QBMMs), which are the primary focus of this work. The presented research can be divided into three major parts. The first part involves the formulation of a common Lagrangian droplet breakup model for QBMMs and the numerical investigation with the QMOM as well as the more sophisticated extended QMOM (EQMOM). The results indicate that the approximations are reasonably accurate when at least six moment equations are solved, with the EQMOM providing no advantages for the investigated configurations. In the second part, a quadrature-based moment model for the effects of fluid turbulence on particle velocities is formulated. The resulting moment equations contain non-smooth integrands that are the source of large errors when using common QBMMs. As an alternative, the Gauss/anti-Gauss QMOM (GaG-QMOM) is proposed that uses the average of a Gaussian and an anti-Gaussian quadrature. Numerical studies show that the GaG-QMOM is able to significantly reduce the previously observed large errors. Another novelty proposed in this context is the modification of the second-order strong-stability preserving Runge-Kutta method to guarantee the preservation of moment realizability in the presence of phase-space diffusion. The third part is concerned with the numerical exploration of the core algorithm of most QBMMs in terms of performance and accuracy. The algorithm consists of, first, computing the recurrence coefficients of the orthogonal polynomials associated with a set of moments, second, solving a symmetric tridiagonal eigenvalue problem to obtain the quadrature nodes and weights, and third, evaluating the integral terms in the moment equations. The results indicate that the contribution of the first step to compute the recurrence coefficients from moments to the overall computational costs is negligible. Instead, the primary focus should be on the fast solution of the eigenvalue problem and, possibly, on the efficient implementation of the moment source term evaluation, which becomes important when second-order processes are considered.
The Taylor-Couette (TC) flow, the flow confined between two concentric independently rotating cylinders, is used as a perfect model to investigate shear flow over concave surfaces and one of the paradigmatic systems of the physics of fluids. In this thesis, an experimental investigation of the turbulent TC flow in a very wide gap geometry with a radius ratio 𝜂 = 0.1 is performed. The physical and dynamic behavior of the flow is studied in a geometry that has rarely been investigated before the current study, which makes this study unique. The study aims to understand the effect of curvature on the TC flow, particularly in cases where the circumferential length of the inner cylinder is smaller than the gap width. The flow is studied in the different rotation regimes: counter-rotating, co-rotating, and purely inner cylinder rotating regimes up to shear Reynolds numbers Re_s≤ 150000. The flow field has been qualitatively studied using visualization techniques. By probing the different flow parameters, familiar coherent TC flow patterns appear, in addition to newly observed patterns we assume only exist for very wide gap TC flows. For a more detailed quantitative study, a time-resolved velocity field measurement has been conducted using the High-speed Particle Image Velocimetry technique through the system end plate. The radial and azimuthal velocity components in the 2D horizontal plane are measured at different axial positions, in order to scan the axial variance of the flow. The recorded flow field is used to compute the angular momentum transport in terms of the quasi-Nusselt number (Nu_ω). The results show a maximum of Nu_ω for low counter-rotating rates of −0.011 ≤ μ_max ≤ −0.0077, which is associated with large-scale structures that span the entire gap. Moreover, the Nu_ω decreases for counter-rotation rates higher than μ_max until it reaches a minimum value and then tends to increase again for higher counter-rotation cases. The space-time behavior of the turbulent flow field for the high counter-rotating cases shows the existence of newly observed patterns next to the outer cylinder wall that propagates inward, enhancing the angular momentum transport and resulting in a second maximum in transport for higher counter-rotating rates. For the pure rotating inner cylinder, the momentum transport scaling with the shear rate Nu_ω ∼ Re_s^α has been studied, and it shows a transition in scaling to 𝛼 = −0.76 for all flows with Re_(s )≥ 25000. This new scaling reveals the transition of the flow from the classical turbulent regime to the ultimate one, where this transition is accompanied by a clear change in the flow behavior. Moreover, the flow in the co-rotating regime and particularly in the centrifugal stable regime (𝜇 > +0.01) is investigated. The Velocity measurements show the presence of disturbed flow near the inner cylinder, where the measured velocity profiles showed a clear deviation from those predicted by laminar flow for flows up to 𝜇 = +0.04.
Precession driven flows are believed to play a relevant role in planetary dynamics, such as in atmospheric phenomena, and as a complementary energy source for homogeneous dynamo action, i.e. the self-generation of planetary magnetic fields. Precessional motion occurs when a body rotates around an axis, which itself is rotating around another axis. The main influence of this forcing mechanism is a gyroscopic effect on the fluid flow which gives rise to a wavy dynamics even in the laminar regime. If the forcing magnitude is strong enough the flow goes through a series of phenomena such as instabilities, resonant interactions between waves, and transition to turbulence whose occurrence depends on the container shape and the angle between the two axis. Although many phenomena have a satisfactory explanation, others still remain elusive and merit further investigations.
The interest in moderate to large forcing is particularly motivated by the need of theoretical supports for the upcoming DRESDYN (DREsden Sodium facility for DYNnamo and thermohydraulic studies) precession experiment, whose main purpose is to test the capability of a precessing fluid system to achieve a dynamo effect. Here, the possibility to generate a magnetic field is connected to the emergence of three large scale structures in the bulk flow: a directly forced standing wave, poloidal vortices, and a geostrophic axisymmetric flow.
In this thesis we use numerical simulations to study and understand the flow behavior in a fluid-filled precessing cylinder. We use two types of approaches: a global study to investigate large scale phenomena and the resulting magnetohydrodynamics behavior, and a local model to analyze and unveil the properties of turbulence forced by precession. The bulk flow behavior present different responses with respect to the sense of motion: while prograde precession shows a steep transition to turbulence when increasing the forcing magnitude with a marked breakdown of the directly forced mode, retrograde precession presents a much smoother change. A related distinction has been found also for the dynamo action, which is more likely to occur for perpendicular and retrograde precession. The precession driven turbulence is a complex scenario determined by the coexistence of geostrophic vortices (called also condensates), a typical feature of rotating turbulence prone to an inverse cascade of energy, and small scale 3D waves characterized by a direct energy cascade. We observe the interaction of these two structures as being governed by a clear hierarchy.
The incompressible temporally developing turbulent boundary layer (TBL) and spatially developing turbulent boundary layer (SBL) with and without blowing is analyzed using the map-based stochastic one-dimensional turbulence (ODT) model. An understanding of these idealized flows is of fundamental relevance for boundary layer-type problems, which are frequently encountered in several applications, from atmospheric sciences to engineering. In the ODT model, the flow variables are resolved on all scales along a wall-normal, 1-D domain. These variables are evolved by a deterministic process representing molecular diffusion and by a stochastic process modeling the effect of turbulent advection and pressure fluctuations. Due to the reduction in dimensionality, the model is particularly appropriate for high Reynolds number flow. It is shown that the ODT model is able to capture salient features of the turbulent boundary layer-type flows by comparing the results with various available reference direct numerical simulation (DNS), large eddy simulation (LES) and experimental results. The comparison is presented for the mean velocity profiles, turbulent velocity fluctuation profiles (up to fourth order), the skin friction coefficient and shape factor for different bulk (Reb) and momentum Reynolds numbers (Reθ) using fixed model parameters. The influence of the model parameters is also discussed for various momentum Reynolds numbers for each investigated flow configuration. The results discussed in this thesis suggest that the ODT model is an economical and reasonably accurate approach for the simulation of turbulent boundary layer flows.
Stratified vortices can be found from small to large scales in geophysical and astrophysical flows. On the one hand, tornadoes and hurricanes can lead to devastation and even a large
number of casualties. On the other hand, vortices can distribute heat and momentum in the atmosphere which is important for a habitable environment on Earth. In the astrophysical context, accretion disks (from which solar systems are formed) can be seen as stratified vortices. In such systems, understanding the mechanisms that can result in an outward transport of angular momentum is a central problem. For a planet or star to be formed in a disk, angular momentum has to be carried away from its center to allow matter aggregation by gravity; otherwise, its rotation speed would be far too large, avoiding this matter aggregation (and the consequent star formation) to happen. In such gas systems, turbulence is the most likely mechanism to achieve such a large angular momentum transport. However, it was shown that the flow profile of accretion disks is stable with respect to purely shear instabilities, and the question arises about how the turbulence can be generated. Among other candidates, the strato-rotational instability (SRI) has attracted attention in recent years. The SRI is a purely hydrodynamic instability that can be modeled by a classical Taylor-Couette (TC) system with stable density stratification due to axial salinity or temperature gradients.
In this thesis, a combined experimental and high-performance computing study of new specific behaviors of the strato-Rotational Instability (SRI) is performed. The density stratification causes a change in the marginal instability transition when compared to classical non-stratified TC systems, making the flow unstable in regions where – without stratification – it would be stable. This characteristic makes the SRI a relevant phenomenon in planetary and astrophysical applications, particularly in accretion disk theory.
Despite many advances in the understanding of strato-rotational flows, the confrontation of experimental data with non-linear numerical simulations remains relevant, since it involves
linear aspects and non-linear interactions of SRI modes which still need to be better understood. These comparisons also reveal new non-linear phenomena and patterns not yet observed in the SRI, that can contribute to our understanding of geophysical flows.
This thesis is an attempt to assess some of the effects that electroquasistatic body forces exert on turbulent internally forced convective flows. In order to do that, a stochastic turbulence model is employed, namely, the One-Dimensional Turbulence (ODT) model.
The reduced dimensionality of ODT demands a reduction of the Navier-Stokes equations (and in this case, also the Maxwell equations), into a 1-D system. This is done by performing an asymptotic analysis in terms of the nondimensional numbers of the flow. Also, a validation step due to the relative novel character of the cylindrical ODT formulation is done for an incompressible and constant properties flow regime and a variable density flow regime. The validation is presented for both the temporal (T-ODT) and a novel spatial (S-ODT) formulation in both planar and cylindrical geometries. Results in the constant property case show that wall normal (and radial) profiles, in both the T-ODT and S-ODT formulations, show good agreement with each other and to the data of Direct Numerical Simulations (DNSs). For the evaluated variable density heated pipe flow, gradients at the wall can be better reproduced with S-ODT.
After validating the model, ODT is applied first into a planar configuration which emulates the flow in a wire-plate Electrostatic Precipitator (ESP). For this flow, the additional input energy due to the electroquasistatic body force has an effect on the modification of the bulk velocity, and subsequently, the skin friction coefficient. Some qualitative DNS trends are confirmed with ODT, such as the localized increase of the Reynolds stress, as a consequence of increased eddy activity close to the discharge electrodes. Next, the results of ODT simulations in a cylindrical wire-tube ESP are presented. Here, ODT results are compared to experimental results. ODT results for global integral quantities such as the streamwise pressure gradient and the Nusselt number enhancement ratio are able to match in a reasonable way the experimental results. The competing relevance between the EHD contribution to turbulence by momentum, and by affecting the temperature and density due to the Joule heating effect is also analyzed, showing the leading order relevance of the former one. Specifically for the Nusselt number results, the sensitivity of the EHD flow to transition effects is shown to be very significant.
This thesis may open the door to a vast new field of phenomena which can not only serve for the further validation of the ODT model against DNSs or experiments, but also for the real use of ODT in applications which are so far inaccessible for traditional DNSs.
Application of ODT to turbulent combustion problems in incompressible and compressible regimes
(2016)
The one-dimensional turbulence (ODT) model is applied to a reactant - to - product counterflow configuration as well as to a shock tube configuration in non-reactive flow and in deflagration and detonation regimes. The model employed herein solves conservation equations for momentum, energy, and species on a one dimensional (1D) domain corresponding to the line spanning the domain between nozzle orifice centers in the counterflow configuration and corresponding to the tube length in the shock tube configuration. The effects of turbulent mixing are modeled via a stochastic process, while the Kolmogorov and reactive length and time scales are explicitly resolved.
In the counterflow configuration, comparisons between model and DNS results for spatial mean and root-mean-square (RMS) velocity, temperature, and major and minor species profiles are shown. The ODT approach shows qualitatively and quantitatively reasonable agreement with the DNS data. Scatter plots and statistics conditioned on temperature are also compared for heat release rate and all species. ODT is able to capture the range of results depicted by DNS. However, conditional statistics show signs of underignition.
To carry out the shock tube simulations, the ODT methodology is extended to include an efficient compressible implementation and a model for capturing shock-induced turbulence is presented. The necessary algorithmic changes to include compressibility effects are highlighted and the model for capturing shock-turbulence interaction is presented. To validate the compressible solver, results for Sod’s shock tube problem are compared against a finite volume Riemann solver. To validate the model for shock-turbulence interaction, comparisons for a non-reactive and a reactive case are presented. First, results of a shock traveling from light (air) to heavy (SF6) with reshock have been simulated to match mixing width growth data of experiments and turbulent kinetic energy results from LES. Then, for one-step chemistry calibrated to represent an acetylene/air mixture, the interaction of a shock wave with an expanding flame front is simulated, and results with 2D simulation (2D-sim) data for flame brush formation and ensuing deflagration-to-detonation transitions (DDT) are compared. Results for the Sod shock tube comparison show that the shock speed and profile are captured accurately. Results for the nonreactive shock-reshock problem show that interface growth at all simulated Mach numbers is captured accurately and that the turbulent kinetic energy agrees in order of magnitude with LES data. The reactive shock tube results show that the flame brush thickness compares well to 2D-sim data and that the approximate location and timing of the DDT can be captured. The known sensitivity of DDT characteristics to details of individual Wow realizations, seen also in ODT, implies that model agreement can be quantified only by comparing Wow ensembles, which are presently unavailable other than in an ODT run-to-run sensitivity study that is reported herein.
Investigations of high Reynolds number pipe flow is up to now a great challenge due to the complex mechanisms which appear in pipe flow turbulence. Hence, suitable experimental facilities are necessary to resolve turbulent dynamics and therewith to provide the knowledge for the understanding of such a simple shear flow. For this reason the recent thesis deals with conceptual design and setup of a new high Reynolds number pipe test facility further on named CoLaPipe - Cottbus Large Pipe. It also comprises first investigations on pipe flow obtained from the new CoLaPipe, which can be classified into 1.)calibration measurements to put the facility into service and 2.)continuative measurements to provide experimental results helping to understand pipe flow.
The first results within the CoLaPipe show that this new experimental facility is suitable to investigate turbulence at high Reynolds numbers, where this conclusion can be drawn from intensive investigations on the development length of the flow either for natural and artificial transition. From further experiments on the evaluation of the wall friction velocity using different estimation methods great difficulties and variations in the calculated values are obtained. These deviations are directly related to the scaling behavior of the mean and fluctuating velocity, which is also shown within this thesis and intensively discussed.
Among the discussion of the setup of the new CoLaPipe and the first experimental results this thesis contains a broad literature review with the focus on high and very high Reynolds numbers. Nevertheless, pipe flow at low and moderate Reynolds numbers is described as well.
The continual optimization process for more efficiency of industrial flows has raised the need for providing deeper understanding of turbulence. These details can be provided by direct numerical simulation (DNS), which is impossible for most flows with current computers. Therefore, progress in optimizing Reynolds averaged Navier-Stokes (RANS) and large eddy simulation (LES) modeling strategies will need to continue. Another ansatz is the reduction to 2D or 1D models to reduce the numerical cost.
One dimensional turbulence (ODT) as presented by A. R. Kerstein is a new modeling strategy that reduces the 3D simulation to a 1D line of sight through the flow region. Due to the higher resolution afforded by the 1D model, it is possible to simulate even the smallest scales and to provide insight into turbulence statistics.
To assess the advantages and disadvantages of the model, ODT has to be validated against several flows. Within this thesis, ODT is validated against the channel flow, the passive scalar transport and the channel flow with a fluctuating pressure gradient. These flows are simplified test cases for the phenomena present in single-phase industrial flows. ODT produces meaningful results for friction Reynolds numbers up to Reτ = 6·10⁵ and for Prandtl numbers from Pr = 0.025 to 50. Statistics of the wall shear stress are presented and the influence of pressure fluctuations is discussed.
Based on these channel results, the non-breaking and breaking jet are simulated. While the former is a simplified case of a free-surface flow, the latter is of primary interest for spray formation and fuel injection. Detailed statistics of the TKE budgets and the breakup are presented. As the last case, the cloud top of a stratocumulus-topped boundary layer (STBL) was simulated. The case combines the interaction of an active and a passive scalar. It further combines the simulation of a stable and an unstable stratified region that suppresses and enhances turbulence respectively. The simulations reproduce the entrainment velocity and generate comparable mean and flux profiles compared to DNSs.
Die Interaktion einer turbulenten Zuströmung mit der Vorderkante eines aerodynamischen Profils ist ein dominierender Schallentstehungsmechanismus, der beispielsweise an Ventilatoren in Klimaanlagen, Windkrafträdern und Rotor-Stator-Konfigurationen in Turbomaschinen auftritt. In der vorliegenden Arbeit wird die breitbandige Schallentstehung infolge dieses Mechanismus experimentell untersucht. Der Versuchsaufbau im aeroakustischen Windkanal besteht aus einem Profil und einem Zylinder mit kreisrundem Querschnitt, der zur Generierung der turbulenten Profilzuströmung stromauf des Profils installiert ist (Zylinder-Profil-Konfiguration). Die Profilzuströmung wird nacheinander durch Nutzung von fünf Zylindern unterschiedlicher Durchmesser und durch das Einstellen fünf verschiedener Abstände zwischen Zylinder und Profilvorderkante variiert. Insgesamt wird die aeroakustische Schallentstehung an der Profilvorderkante für 50 Zylinder-Profil-Konfigurationen untersucht, wobei die Zylinder-Zuströmgeschwindigkeit für jede Konfiguration zusätzlich in 26 Schritten zwischen 26 m/s und 72 m/s vorgegeben wird. Der von den Experimenten abgedeckte Machzahlbereich liegt zwischen 0,075 und 0,21 und der Bereich der mit dem Zylinderdurchmesser gebildeten Reynoldszahl liegt zwischen 8,3·10³ und 7,3·10⁴. Der generierte Schall wurde mit einem Mikrofonarray aufgezeichnet. Der infolge der turbulenten Zuströmung an der Profilvorderkante generierte Schall wird während der Analyse der Mikrofonarraymessdaten mit erweiterten Beamforming-Algorithmen (Entfaltungsverfahren) von den anderen Schallquellen getrennt. Somit ist sichergestellt, dass ausschließlich die Schallentstehung an der Vorderkante des Profils und hiermit die Schallentstehung durch turbulente Zuströmung untersucht wird. Insgesamt wurden über 2000 einzelne akustische Messungen durchgeführt, die zusammengenommen ein Datenvolumen von ca. 1,7 TByte aufweisen. Neben den akustischen Messungen wurden für ausgewählte Zylinder-Profil-Konfigurationen Hitzdrahtmessungen durchgeführt, um die Profilzuströmung im Bereich der Profilvorderkante zu charakterisieren. Auf Basis der Analyse der umfangreichen Messdaten wird mit Hilfe einer Dimensionsanalyse ein empirisches Modell zur Abschätzung der Schallentstehung bei inhomogener turbulenter Zuströmung entwickelt. Zudem werden einige für homogene Zuströmung vorhandene Modelle für inhomogene Zuströmbedingungen adaptiert. Die Eignung der getesteten Modelle wird anhand des Vergleichs von Modell- und Messergebnissen diskutiert.