06 Fakultät Luft- und Raumfahrttechnik und Geodäsie

Permanent URI for this collectionhttps://elib.uni-stuttgart.de/handle/11682/7

Browse

Search Results

Now showing 1 - 10 of 27
  • Thumbnail Image
    ItemOpen Access
    Linear stability investigations of three-dimensional disturbances in the boundary layer over anisotropic compliant walls
    (2023) Zengl, Marcus; Rist, Ulrich (apl. Prof. Dr.-Ing.)
    In dieser Arbeit werden dreidimensionale Störungen in der Grenzschicht über anisotropen nachgiebigen Wänden mit linearer Stabilitätstheorie untersucht. Ein oberflächenbasiertes Modell wird verwendet, um die nachgiebige Wand abzubilden. Hierbei wird das anisotrope Wandmodell von Carpenter erweitert, um einen zusätzlichen Schiebewinkel der Wand bezüglich der Strömungsrichtung einzubringen. Basierend auf diesem Wandmodell wird eine Randbedingung für die Lineare Stabilitätstheorie hergeleitet. Aufgrund der Tatsache, dass diese Randbedingung die Orr-Sommerfeld- und Squire-Gleichung koppelt, wurden zwei neuartige Lösungsverfahren, ein Schießverfahren und ein Matrixlöser, für diesen besonderen Umstand entwickelt. Der Schießlöser transformiert das zugrunde gelegte Eigenwertproblem in ein Randwertproblem und verwendet ein klassisches Schießverfahren zur Lösung des Problems. Um das numerisch steife Problem mit seinem parasitärem Fehlerwachstum zu berücksichtigen beinhaltet das Lösungsverfahren eine Gram-Schmid Orthonormierungsroutine. Durch eine neuartige Skalierung der Phase des zu minimierenden Residuums wird das zeitliche und räumliche Modell robust und performant für gegebene Eigenmoden gelöst. Das durch die gekoppelte Orr-Sommerfeld- and Squire-Gleichung entstehende Eigenwertproblem wird auch mit einer Matrix-basierenden Methode gelöst. Das durch die nachgiebige Wand entstehende zeitliche quadratische Eigenwertproblem wird dabei berücksichtigt. Hierbei wird eine pseudospektrale Diskretisierung mit Chebyshev-Kollokation verwendet. Besonders betrachtet wird die Formulierung des diskretisierten Problems auf seine numerischen Fehler. Die numerische Genauigkeit der Lösungsverfahren wird genau überprüft, um die Gitterunabhängigkeit der Ergebnisse sicherzustellen. Um das Potenzial der nachgiebigen Wände zur Verzögerung des laminar-turbulenten Umschlags zu untersuchen, wurde die Vorgehensweise von Carpenter [15] übernommen. Carpenter optimierte die Parameter der nachgiebigen Wand so, dass Tollmien-Schlichting (TS) Moden so weit wie möglich abgeschwächt werden, während Fluid-Struktur (FISI) Moden grenzwertig stabil bleiben. Dieses Vorgehen wurde ausgewählt, weil Fluid-Struktur Moden absolut instabil sein können, was zu sofortigem Strömungsumschlag führen kann. Stabilitätsrechnungen wurden ausgeführt für zwei Sätze von Wandparametern, die Carpenter mit seinem zweidimensionalen Rahmenwerk optimiert hat. Hierbei wurden nicht nur dreidimensionale Störungen betrachtet, sondern es wurde auch der Einfluss des neu eingebrachten Schiebewinkels der nachgiebigen Wand untersucht. Die Ergebnisse wurden bezüglich der zeitlichen Anfachung der TS- und FISI-Moden, und bezüglich des mit N-Faktoren vorhergesagten Umschlagspunkts beurteilt. Es wird gezeigt, dass dreidimensionale Störungen bestimmte N-Faktoren vor ihren zweidimensionalen Pendants erreichen. Die vorhergesagte laminare Länge ist etwas kürzer als mit zweidimensionalen Verfahren vorhergesagt. Es scheint als ob der eingebrachte Schiebewinkel für die untersuchten Parametersätze keinen Vorteil bezüglich Laminarhaltung bringt. Schließlich wurden optimale Störungen berechnet, um das Transiente Energiewachstum für die anisotrope nachgiebige Wand zu untersuchen. Hierbei wurden die Anfangsverteilungen von Eigenmoden so optimiert, dass deren Überlagerung ein maximales Energiewachstum für eine vorgegebene Zeit erfährt. Die Einhüllende dieser optimalen Störungen wird dann für variierende Wellenzahlen in Strömungs- und Spannweitenrichtung, und variierende Wachstumszeit berechnet. Die Ergebnisse zeigen kein durch die nachgiebige Wand hervorgerufenes relevantes transientes Wachstum. Es wird gezeigt, dass der klassische Mechanismus für transientes Wachstum, der bei der steifen Wand dominiert, nicht verändert wird.
  • Thumbnail Image
    ItemOpen Access
    Laminar-to-turbulent transition in airfoil boundary layer flows at oscillating inflow conditions
    (2023) Ohno, Duncan; Rist, Ulrich (apl. Prof. Dr.-Ing.)
    This thesis numerically investigates the laminar-to-turbulent transition in boundary-layer flows on natural laminar flow airfoils under oscillating inflow conditions. Large-scale fluctuations in the form of periodic vertical gusts generate an oscillating pressure gradient, resulting in a complex transient behavior of the boundary layer. Under these conditions, two scenarios are investigated: an attached flow with natural Tollmien-Schlichting (TS) wave transition and a boundary-layer flow featuring a laminar separation bubble (LSB). The study aims to provide a deeper understanding of the transient mechanisms involved as well as the basis for new transition prediction methods for unsteady conditions. Direct numerical simulations (DNS) are performed where the gust disturbance is imposed on the fully-resolved transitional boundary layers via unsteady boundary conditions. In this novel approach, transient base flows are generated in advance with unsteady Reynolds-averaged Navier-Stokes (URANS) simulations of entire unsteady airfoil flows in conjunction with the disturbance velocity approach (DVA) to introduce sinusoidal gusts. The spatio-temporal evolution of the modal disturbances is analyzed using the continuous wavelet transform (CWT), which is then compared with linear stability theory (LST) by employing a trajectory-following method for transient flows. Several physical effects responsible for the transient characteristics of the studied flows are identified and observations from previous experimental studies are classified. The results demonstrate that the quasi-steady linear theory adequately predicts the transient behavior of the convective modal disturbances for a wide range of gust perturbations. The so-called convective-transition mode with a subsequent calmed region is found for cases with a high degree of unsteadiness. This study is the first to provide a physical explanation for the occurrence of this mechanism for natural transition.
  • Thumbnail Image
    ItemOpen Access
    A time-accurate inflow coupling for zonal LES
    (2023) Blind, Marcel P.; Kleinert, Johannes; Lutz, Thorsten; Beck, Andrea
    Generating turbulent inflow data is a challenging task in zonal large eddy simulation (zLES) and often relies on predefined DNS data to generate synthetic turbulence with the correct statistics. The more accurate, but more involved alternative is to use instantaneous data from a precursor simulation. Using instantaneous data as an inflow condition allows to conduct high fidelity simulations of subdomains of, e.g. an aircraft including all non-stationary or rare events. In this paper, we introduce a toolchain that is capable of interchanging highly resolved spatial and temporal data between flow solvers with different discretization schemes. To accomplish this, we use interpolation algorithms suitable for scattered data in order to interpolate spatially. In time, we use one-dimensional interpolation schemes for each degree of freedom. The results show that we can get stable simulations that map all flow features from the source data into a new target domain. Thus, the coupling is capable of mapping arbitrary data distributions and formats into a new domain while also recovering and conserving turbulent structures and scales. The necessary time and space resolution requirements can be defined knowing the resolution requirements of the used numerical scheme in the target domain.
  • Thumbnail Image
    ItemOpen Access
    Gust alleviation by spanwise load control applied on a forward and backward swept wing
    (2023) Klug, Lorenz; Ullah, Junaid; Lutz, Thorsten; Streit, Thomas; Heinrich, Ralf; Radespiel, Rolf
    The present paper investigates the feasibility of gust load alleviation at transonic speeds on a backward swept and a forward swept transport aircraft configuration. Spanwise-distributed control surfaces at the leading and trailing edges are employed to control gust-induced wing bending as well as wing torsion moments. The deflection amplitude and temporal flap actuation are determined by a novel scheme that builds on the aerodynamic strip theory. The aerodynamic effectiveness of the actuators is taken from a data base, computed from either 2D infinite swept wing simulations, or from yawed computations that take the effects of boundary-layer cross flow and the local sweep angle of the control surface into account. The present numerical flow simulations reveal that careful application of control laws at the trailing edge alleviates wing bending moments caused by strong vertical gusts by 85-90%, for both aircraft configurations. The application of leading-edge flaps introduces significant nonlinear aerodynamic interactions, that make the control of torsional moments comparably challenging. Here, the present results indicate that about 60% of wing torsion loads due to strong gusts can be removed.
  • Thumbnail Image
    ItemOpen Access
    A reinforcement learning based slope limiter for second‐order finite volume schemes
    (2023) Schwarz, Anna; Keim, Jens; Chiocchetti, Simone; Beck, Andrea
    Hyperbolic equations admit discontinuities in the solution and thus adequate and physically sound numerical schemes are necessary for their discretization. Second‐order finite volume schemes are a popular choice for the discretization of hyperbolic problems due to their simplicity. Despite the numerous advantages of higher‐order schemes in smooth regions, they fail at strong discontinuities. Crucial for the accurate and stable simulation of flow problems with discontinuities is the adequate and reliable limiting of the reconstructed slopes. Numerous limiters have been developed to handle this task. However, they are too dissipative in smooth regions or require empirical parameters which are globally defined and test case specific. Therefore, this paper aims to develop a new slope limiter based on deep learning and reinforcement learning techniques. For this, the proposed limiter is based on several admissibility constraints: positivity of the solution and a relaxed discrete maximum principle. This approach enables a slope limiter which is independent of a manually specified global parameter while providing an optimal slope with respect to the defined admissibility constraints. The new limiter is applied to several well‐known shock tube problems, which illustrates its broad applicability and the potential of reinforcement learning in numerics.
  • Thumbnail Image
    ItemOpen Access
    Investigation of a realistic flap modeling using a combination of Chimera method and grid deformation on a wing fuselage configuration
    (2023) Hillebrand, Marco; Müller, Jens; Ullah, Junaid; Lutz, Thorsten
    Flap deflections of an aircraft wing for active load alleviation within CFD simulations are realized using pure grid deformation due to time saving and low modeling complexity. In this case, spanwise gaps are neglected, which are present in reality during a flap deflection. Another possibility to realize the deflections is the combination of pure grid deformation and Chimera method, which allows the modeling of the gap between flap and wing or consecutive flaps. The overall aim of this work is the analysis of the aerodynamic effects caused by the different modeling approaches realizing leading and trailing edge flap deflections. The comparison of the modeling methods is investigated on the DLR LEISA configuration, which is a generic wing‐fuselage configuration. For active gust load alleviation, the leading edge flaps are deflected downward and the trailing edge flaps are deflected upward. Due to the downward deflection of the leading edge flaps, vortices are formed using the combined Chimera method as a result of the gap consideration. These vortices lead to a local drag increase resulting in a difference between both modeling methods in the spanwise as well as global drag coefficient. With the pure grid deformation these vortices do not occur. Due to the upward trailing edge deflection, the combined Chimera method leads to a pressure compensation via the effective gap enlargement, which is not present in the pure grid deformation. Overall, the combined Chimera method offers a good possibility to model the induced drag as well as the pressure compensation at a large flap deflection.
  • Thumbnail Image
    ItemOpen Access
    Experiments on the laminar to turbulent transition under unsteady inflow conditions
    (2023) Romblad, Jonas; Krämer, Ewald (Prof. Dr.-Ing.)
    Natural laminar flow (NLF) airfoils are key to the performance of sailplanes and wind turbines. They provide a significant reduction of friction drag by delaying the transition from laminar to turbulent boundary layer. However, the most common method for transition prediction in NLF airfoil design, the e^n method (Mack 1977), has limited capabilities for taking inflow turbulence into account. The current work employs wind tunnel experiments to study how the transition on an NLF airfoil is affected by free-stream turbulence. The effect of small- and large-scale turbulence is studied separately, as well as in combination. In the wind tunnel, turbulence grids generate small-scale turbulence, and a gust generator induces inflow angle oscillations corresponding to large-scale turbulence. The study includes a detailed characterization of the turbulence generated by grids placed in the settling chamber of the wind tunnel. The Reynolds number Re = 3400000 and the airfoil pressure distribution is matched to cruise or dash flight of general aviation aircraft. The results are compared with direct numerical simulation, linear stability theory (LST) and flight measurements. The results show that small-scale turbulence does have an influence on the transition location in the investigated range of turbulence level, 0.01% < Tu < 0.11%. The modified e^n method (Mack 1977) captures the general trend, but the sensitivity to Tu is airfoil dependent. The effects of 2D, single-mode inflow angle oscillations are investigated in the range of reduced frequency 0.06 < kappa < 1.7. In this range, the transition process changes from quasi-steady to clearly unsteady, but a fully convective transition mode is not formed. This is an intermediate range of unsteady flow in which trajectory-following LST is able to capture the main features of the unsteady transition process. No significant interaction between the effects of small- and large-scale turbulence are observed in the investigated range of Tu and kappa, indicating that the effects can be superposed.
  • Thumbnail Image
    ItemOpen Access
    The near-wake development of a wind turbine operating in stalled conditions : part 1: assessment of numerical models
    (2024) Weihing, Pascal; Cormier, Marion; Lutz, Thorsten; Krämer, Ewald
    This study comprehensively investigates the near-wake development of a model wind turbine operating at a low tip-speed ratio in stalled conditions. In the present paper, part 1, different ways of representing the turbine, which include a full geometrical representation and modeling by means of the actuator line method, and different approaches for the modeling of turbulence are assessed. The simulation results are compared with particle image velocimetry (PIV) measurements from the MEXICO and New MEXICO experiments. A highly resolved numerical setup was created and a higher-order numerical scheme was applied to target an optimal resolution of the tip vortex development and the wakes of the blades. Besides the classical unsteady Reynolds-averaged methodology, a recently developed variant of the detached-eddy simulation (DES) was employed, which features robust shielding capabilities of the boundary layers and enhanced transition to a fully developed large-eddy simulation (LES) state. Two actuator line simulations were performed in which the aerodynamic forces were either evaluated by means of tabulated data or imposed from the averaged blade loads of the simulation with full blade geometry. The purpose is to distinguish between the effects of the force projection and the force calculation in the underlying blade-element method on the blade wake development. With the hybrid Reynolds-averaged Navier-Stokes (RANS)-LES approach and the geometrically fully resolved rotor blade, the details of the flow of the detached blade wake could be resolved. The prediction of the wake deficit also agreed very well with the experimental data. Furthermore, the strength and size of the blade tip vortices were correctly predicted. With the linear unsteady Reynolds-averaged Navier-Stokes (URANS) model, the wake deficit could also be described correctly, yet the size of the tip vortices was massively overestimated. The actuator line method, when fed with forces from the fully resolved simulation, provides very similar results in terms of wake deficit and tip vortices to its fully resolved parent simulation. However, using uncorrected two-dimensional polars shows significant deviations in the wake topology of the inner blade region. This shows that the application in such flow conditions requires models for rotational augmentation. In part 2 of the study, to be published in another paper, the development and the dynamics of the early tip vortex formation are detailed.
  • Thumbnail Image
    ItemOpen Access
    Assessment of low‐frequency aeroacoustic emissions of a wind turbine under rapidly changing wind conditions based on an aero‐servo‐elastic CFD simulation
    (2023) Wenz, Florian; Maas, Oliver; Arnold, Matthias; Lutz, Thorsten; Krämer, Ewald
    A meteorologically challenging situation that represents a demanding control task (rotational speed, pitch and yaw) for a wind turbine is presented and its implementation in a simulation is described. A high-fidelity numerical process chain, consisting of the computational fluid dynamics (CFD) solver FLOWer, the multi-body system (MBS) software SIMPACK and the Ffowcs Williams-Hawkings code ACCO, is used. With it, the aerodynamic, servoelastic and aeroacoustic (<20 Hz) behaviour of a generic wind turbine during a meteorological event with strong and rapid changes in wind speed and direction is investigated. A precursor simulation with the meteorological model system PALM is deployed to generate realistic inflow data. The simulated strong controller response of the wind turbine and the resulting aeroelastic behaviour are analysed. Finally, the low-frequency sound emissions are evaluated and the influence of the different operating and flow parameters during the variable inflow is assessed. It is observed that the wind speed and, linked to it, the rotational speed as well as the turbulence intensity are the main influencing factors for the emitted low-frequency sound power of the wind turbine. Yawed inflow, on the other hand, has little effect unless it changes the operational mode to load reduction, resulting in a swap of the main emitter from the blades to the tower.