ترغب بنشر مسار تعليمي؟ اضغط هنا

Direct Numerical Simulation of a Turbulent Boundary Layer with Strong Pressure Gradients

67   0   0.0 ( 0 )
 نشر من قبل Riccardo Balin
 تاريخ النشر 2020
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

The turbulent boundary layer over a Gaussian shaped bump is computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. The two-dimensional bump causes a series of strong pressure gradients alternating in rapid succession. At the inflow, the momentum thickness Reynolds number is approximately 1,000 and the boundary layer thickness is 1/8 of the bump height. DNS results show that the strong favorable pressure gradient (FPG) causes the boundary layer to enter a relaminarization process. The near-wall turbulence is significantly weakened and becomes intermittent, however relaminarization does not complete. The streamwise velocity profiles deviate above the standard logarithmic law and the Reynolds shear stress is reduced. The strong acceleration also suppresses the wall-shear normalized turbulent kinetic energy production rate. At the bump peak, where the FPG switches to an adverse gradient (APG), the near-wall turbulence is suddenly enhanced through a partial retransition process. This results in a new highly energized internal layer which is more resilient to the strong APG and only produces incipient flow separation on the downstream side. In the strong FPG and APG regions, the inner and outer layers become largely independent of each other. The near-wall region responds to the pressure gradients and determines the skin friction. The outer layer behaves similarly to a free-shear layer subject to pressure gradients and mean streamline curvature effects. Results from a RANS simulation of the bump are also discussed and clearly show the lack of predictive capacity of the near-wall pressure gradient effects on the mean flow.



قيم البحث

اقرأ أيضاً

A study of large-scale motions from a new direct numerical simulation database of the turbulent boundary layer up to Re{theta} ~ 2500 is presented. The statistics of large-scale streamwise structures have been investigated using two-dimensional and t hree-dimensional extraction procedures. The large-scale structures are abstracted using a robust skeletonization method usually applied to other research domains to simplify complex 3D objects. Different structure parameters such as the length, shape or angle are investigated. The features of the detected structures are compared to their mean counterparts extracted from two-point correlations. Structures as large as 10 boundary layer thickness are observed. The streamwise length of these structures follows a -2 power law distribution, similar to the experimental findings at higher Reynolds numbers.
The turbulent boundary layer over a flat plate is computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations as a test bed for a synthetic turbulence generator (STG) inflow boundary condition. The inlet momentum thick ness Reynolds number is approximately 1,000. The study provides validation of the ability of the STG to develop accurate turbulence in 5 to 7 boundary layer thicknesses downstream of the boundary condition. Also tested was the effect of changes in the stabilization scheme on the development of the boundary layer. Moreover, the grid resolution required for both the development region and the downstream flow is investigated when using a stabilized finite element method.
The present numerical investigation uses well-resolved large-eddy simulations to study the low-frequency unsteady motions observed in shock-wave/turbulent-boundary-layer interactions. Details about the numerical aspects of the simulations and the sub sequent data analysis can be found in three papers by the authors (Theo. Comput. Fluid Dyn., 23:79--107 (2009); Shock Waves, 19(6):469--478 (2009) and J. of Fluid Mech. (2011)). The fluid dynamics video illustrates the complexity of the interaction between a Mach 2.3 supersonic turbulent boundary layer and an oblique shock wave generated by a 8-degree wedge angle. The first part of the video highlights the propagation of disturbances along the reflected shock due to the direct perturbation of the shock foot by turbulence structures from the upstream boundary layer. The second part of the video describes the observed block-like back-and-forth motions of the reflected shock, focusing on timescales about two orders of magnitude longer than the ones shown in the first part of video. This gives a visual impression of the broadband and energetically-significant peak in the wall-pressure spectrum at low frequencies. The background blue-white colouring represents the temperature field (with white corresponding to hot) and one can clearly appreciate why such low-frequency shock motions can lead to reduced fatigue lifetimes and is detrimental to aeronautical applications.
The numerical simulation of a flow through a duct requires an externally specified forcing that makes the fluid flow against viscous friction. To this aim, it is customary to enforce a constant value for either the flow rate (CFR) or the pressure gra dient (CPG). When comparing a laminar duct flow before and after a geometrical modification that induces a change of the viscous drag, both approaches (CFR and CPG) lead to a change of the power input across the comparison. Similarly, when carrying out the (DNS and LES) numerical simulation of unsteady turbulent flows, the power input is not constant over time. Carrying out a simulation at constant power input (CPI) is thus a further physically sound option, that becomes particularly appealing in the context of flow control, where a comparison between control-on and control-off conditions has to be made. We describe how to carry out a CPI simulation, and start with defining a new power-related Reynolds number, whose velocity scale is the bulk flow that can be attained with a given pumping power in the laminar regime. Under the CPI condition, we derive a relation that is equivalent to the Fukagata--Iwamoto--Kasagi relation valid for CFR (and to its extension valid for CPG), that presents the additional advantage of natively including the required control power. The implementation of the CPI approach is then exemplified in the standard case of a plane turbulent channel flow, and then further applied to a flow control case, where the spanwise-oscillating wall is used for skin friction drag reduction. For this low-Reynolds number flow, using 90% of the available power for the pumping system and the remaining 10% for the control system is found to be the optimum share that yields the largest increase of the flow rate above the reference case, where 100% of the power goes to the pump.
A direct numerical simulation (DNS) of a channel flow with one curved surface was performed at moderate Reynolds number (Re_tau = 395 at the inlet). The adverse pressure gradient was obtained by a wall curvature through a mathematical mapping from ph ysical coordinates to Cartesian ones. The code, using spectral spanwise and normal discretization, combines the advantage of a good accuracy with a fast integration procedure compared to standard numerical procedures for complex geometries. The turbulent flow slightly separates on the profile at the lower curved wall and is at the onset of separation at the opposite flat wall. The thin separation bubble is characterized with a reversal flow fraction. Intense vortices are generated near the separation line on the lower wall but also at the upper wall. Turbulent normal stresses and kinetic energy budget are investigated along the channel.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا