06 Fakultät Luft- und Raumfahrttechnik und Geodäsie
Permanent URI for this collectionhttps://elib.uni-stuttgart.de/handle/11682/7
Browse
41 results
Search Results
Item Open Access Numerical study on the aerodynamic characteristics of the NACA 0018 airfoil at low Reynolds number for Darrieus wind turbines using the Transition SST model(2021) Rogowski, Krzysztof; Królak, Grzegorz; Bangga, GalihA symmetrical NACA 0018 airfoil is often used in such applications as small-to-medium scale vertical-axis wind turbines and aerial vehicles. A review of the literature indicates a large gap in experimental studies of this airfoil at low and moderate Reynolds numbers in the previous century. This gap has limited the potential development of classical turbulence models, which in this range of Reynolds numbers predict the lift coefficients with insufficiently accurate results in comparison to contemporary experimental studies. Therefore, this paper validates the aerodynamic performance of the NACA 0018 airfoil and the characteristics of the laminar separation bubble formed on its suction side using the standard uncalibrated four-equation Transition SST turbulence model and the unsteady Reynolds-averaged Navier-Stokes (URANS) equations. A numerical study was conducted for the chord Reynolds number of 160,000, angles of attack between 0 and 11 degrees, as well as for the free-stream turbulence intensity of 0.05%. The calculated lift and drag coefficients, aerodynamic derivatives, as well as the location and length of the laminar bubble quite well agree with the results of experimental measurements taken from the literature for validation. A sensitivity study of the numerical model was performed in this paper to examine the effects of the time-step size, geometrical parameters and mesh distribution around the airfoil on the simulation results. The airfoil data sets obtained in this work using the Transition SST and the k-ω SST turbulence models were used in the improved double multiple streamtube (IDMS) to calculate aerodynamic blade loads of a vertical-axis wind turbine. The characteristics of the normal component of the aerodynamic blade load obtained by the Transition SST approach are much better suited to the experimental data compared to the k-ω SST turbulence model.Item Open Access The effects of airfoil thickness on dynamic stall characteristics of high‐solidity vertical axis wind turbines(2021) Bangga, Galih; Hutani, Surya; Heramarwan, HenidyaThe flow physics of high solidity vertical axis wind turbines (VAWTs) is influenced by the dynamic stall effects. The present study is aimed at investigating the effects of airfoil thickness on the unsteady characteristics of high solidity VAWTs. Seven different national advisory committee for aeronautics (NACA) airfoils (0008, 0012, 0018, 0021, 0025, 0030, 0040) are investigated. A high fidelity computational fluid dynamics (CFD) approach is used to examine the load and flow characteristics in detail. Before the study is undertaken, the CFD simulation is validated with experimental data as well as large eddy simulation results with sound agreement. The investigation demonstrates that increasing the airfoil thickness is actually beneficial not only for suppressing the dynamic stall effects but also to improve the performance of high solidity turbines. Interestingly this is accompanied by a slight reduction in thrust component. The strength and radius of the dynamic stall vortex decrease with increasing airfoil thickness. The airfoil thickness strongly influences the pressure distributions during dynamic stall process, which is driven by the suction peak near the leading edge. The knowledge gained might be used by blade engineers for designing future turbines and for improving the accuracy of engineering models.Item Open Access A time-accurate inflow coupling for zonal LES(2023) Blind, Marcel P.; Kleinert, Johannes; Lutz, Thorsten; Beck, AndreaGenerating 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.Item Open Access Piloted simulation of the rotorcraft wind turbine wake interaction during hover and transit flights(2022) Štrbac, Alexander; Greiwe, Daniel Heinrich; Hoffmann, Frauke; Cormier, Marion; Lutz, ThorstenHelicopters are used for offshore wind farms for maintenance and support flights. The number of helicopter operations is increasing with the expansion of offshore wind energy, which stresses the point that the current German regulations have not yet been validated through scientific analysis. A collaborative research project between DLR, the Technical University of Munich, the University of Stuttgart and the University of Tübingen has been conducted to examine the sizes of the flight corridors on offshore wind farms and the lateral safety clearance for helicopter hoist operations at offshore wind turbines. This paper details the results of piloted helicopter simulations in a realistic offshore wind farm scenario. The far-wake of rotating wind turbines and the near-wake of non-rotating wind turbines have been simulated with high-fidelity computational fluid dynamics under realistic turbulent inflow conditions. The resulting flow fields have been processed by superposition during piloted simulations in the research flight simulator AVES to examine the flight corridors in transit flights and the lateral safety clearance in hovering flights. The results suggest a sufficient size for the flight corridor and sufficient lateral safety clearance at the offshore wind turbines in the considered scenarios.Item Open Access DNS of a turbulent boundary layer using inflow conditions derived from 4D-PTV data(2021) Appelbaum, Jason; Ohno, Duncan; Rist, Ulrich; Wenzel, ChristophUnsteady, 3D particle tracking velocimetry (PTV) data are applied as an inlet boundary condition in a direct numerical simulation (DNS). The considered flow case is a zero pressure gradient (ZPG) turbulent boundary layer (TBL) flow over a flat plate. The study investigates the agreement between the experimentally measured flow field and its simulated counterpart with a hybrid 3D inlet region. The DNS field inherits a diminishing contribution from the experimental field within the 3D inlet region, after which it is free to spatially evolve. Since the measurement does not necessarily provide a spectrally complete description of the turbulent field, the spectral recovery of the flow field is analyzed as the TBL evolves. The study summarizes the pre-processing methodology used to bring the experimental data into a form usable by the DNS as well as the numerical method used for simulation. Spectral and mean flow analysis of the DNS results show that turbulent structures with a characteristic length on the order of one average tracer particle nearest neighbor radius r¯NN or greater are well reproduced and stay correlated to the experimental field downstream of the hybrid inlet. For turbulent scales smaller than r¯NN, where experimental data are sparse, a relatively quick redevelopment of previously unresolved turbulent energy is seen. The results of the study indicate applicability of the approach to future DNS studies in which specific upstream or far field boundary conditions (BCs) are required and may provide the utility of decreasing high initialization costs associated with conventional inlet BCs.Item Open 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, RolfThe 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.Item Open Access Simulation of transonic buffet with an automated zonal DES approach(2020) Ehrle, Maximilian; Waldmann, Andreas; Lutz, Thorsten; Krämer, EwaldA study of transonic buffet on the NASA Common Research Model at flight Reynolds numbers is presented. The ability of two different hybrid RANS/LES models as well as the URANS approach for resolving three-dimensional buffet motion was evaluated by means of spectral analysis. Automated Zonal DES and URANS simulations show similar results in terms of buffet frequency and spanwise propagation of buffet cells, whereas the Delayed Detached Eddy Simulation results indicate a strong interaction between flow separation and shock motion. The extracted characteristic frequencies which are associated with transonic buffet are located in a range of Sr = 0.2-0.65 for URANS and AZDES and are therefore in accordance with findings from related recent research. Furthermore, the simulation time series were investigated and a structure of spanwise moving buffet cells with varying convection speed and wavelength could be observed.Item Open Access A reinforcement learning based slope limiter for second‐order finite volume schemes(2023) Schwarz, Anna; Keim, Jens; Chiocchetti, Simone; Beck, AndreaHyperbolic 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.Item Open Access Entropy stable discontinuous Galerkin schemes on moving meshes for hyperbolic conservation laws(2020) Schnücke, Gero; Krais, Nico; Bolemann, Thomas; Gassner, Gregor J.This work is focused on the entropy analysis of a semi-discrete nodal discontinuous Galerkin spectral element method (DGSEM) on moving meshes for hyperbolic conservation laws. The DGSEM is constructed with a local tensor-product Lagrange-polynomial basis computed from Legendre-Gauss-Lobatto points. Furthermore, the collocation of interpolation and quadrature nodes is used in the spatial discretization. This approach leads to discrete derivative approximations in space that are summation-by-parts (SBP) operators. On a static mesh, the SBP property and suitable two-point flux functions, which satisfy the entropy condition from Tadmor, allow to mimic results from the continuous entropy analysis, if it is ensured that properties such as positivity preservation (of the water height, density or pressure) are satisfied on the discrete level. In this paper, Tadmor’s condition is extended to the moving mesh framework. We show that the volume terms in the semi-discrete moving mesh DGSEM do not contribute to the discrete entropy evolution when a two-point flux function that satisfies the moving mesh entropy condition is applied in the split form DG framework. The discrete entropy behavior then depends solely on the interface contributions and on the domain boundary contribution. The interface contributions are directly controlled by proper choice of the numerical element interface fluxes. If an entropy conserving two-point flux is chosen, the interface contributions vanish. To increase the robustness of the discretization we use so-called entropy stable two-point fluxes at the interfaces that are guaranteed entropy dissipative and thus give a bound on the interface contributions in the discrete entropy balance. The remaining boundary condition contributions depend on the type of the considered boundary condition. E.g. for periodic boundary conditions that are of entropy conserving type, our methodology with the entropy conserving interface fluxes is fully entropy conservative and with the entropy stable interface fluxes is guaranteed entropy stable. The presented proof does not require any exactness of quadrature in the spatial integrals of the variational forms. As it is the case for static meshes, these results rely on the assumption that additional properties like positivity preservation are satisfied on the discrete level. Besides the entropy stability, the time discretization of the moving mesh DGSEM will be investigated and it will be proven that the moving mesh DGSEM satisfies the free stream preservation property for an arbitrary s-stage Runge–Kutta method, when periodic boundary conditions are used. The theoretical properties of the moving mesh DGSEM will be validated by numerical experiments for the compressible Euler equations with periodic boundary conditions.Item Open Access Characterization of low levels of turbulence generated by grids in the settling chamber of a laminar wind tunnel(2022) Romblad, Jonas; Greiner, Michael; Guissart, Amandine; Würz, WernerAbstractWind tunnel investigations of how Natural Laminar Flow (NLF) airfoils respond to atmospheric turbulence require the generation of turbulence, whose relevant characteristics resemble those in the atmosphere. The lower, convective part of the atmospheric boundary layer is characterized by low to medium levels of turbulence. The current study focuses on the small scales of this turbulence. Detailed hot-wire measurements have been performed to characterize the properties of the turbulence generated by grids mounted in the settling chamber of the Laminar Wind Tunnel (LWT). In the test section, the very low base turbulence level of Tuu ≅ 0.02% (10 ≤ f ≤ 5000 Hz) is incrementally increased by the grids up to Tuu ≅ 0.5%. The turbulence spectrum in the u-direction shows the typical suppression of larger scales due to the contraction between grids and test section. Still, the generated turbulence provides a good mapping of the spectrum measured in flight for most of the frequency range 500 ≤ f ≤ 3000 Hz, where Tollmien-Schlichting (TS)-amplification occurs for typical NLF airfoils. The spectra in v and w-direction exhibit distinct inertial subranges with slopes being less steep compared to the - 5/3 slope of the Kolmogorov spectrum. The normalized spectra in u-direction collapse together well for all grids, whereas in v- and w-directions the inertial- and dissipative subranges are more clearly distinguished for the coarser grids. It is demonstrated that the dissipation rate ε is a suitable parameter for comparing the wind tunnel turbulence with the atmospheric turbulence in the frequency range of interest. By employing the grids, turbulence in the range 4.4 × 10-7 ≤ ε ≤ 0.40 m2/s3 at free-stream velocity U∞ = 40 m/s can be generated in the LWT, which covers representative dissipation rates of free flight NLF applications. In the x-direction, the spectra of the v and w-components develop progressively more pronounced inertial- and dissipative subranges, and the energy below f ≈ 400 Hz decreases. In contrast, the spectral energy of the u-component increases across the whole frequency range, when moving downstream. This behavior can be explained by the combination of energy transport along the Kolmogorov cascade and the incipient return to an isotropic state.Graphic Abstract