scholarly journals The effort of increasing Reynolds number in POD-Galerkin Reduced Order Methods: from laminar to turbulent flows

Author(s):  
Shafqat Ali ◽  
Saddam Hijazi ◽  
Sokratia Georgaka ◽  
Francesco Ballarin ◽  
Giovanni Stabile ◽  
...  

We present different strategies to be able to increase Reynolds number in Reduced Order Methods (ROMs), from laminar to turbulent flows, in the context of the incompressible parametrised Navier-Stokes equations. The proposed methodologies are based on different full order discretisation techniques: the finite element method and the finite volume method. For what concerns finite element full order discretisations which in this work aim to be used from low to moderate Reynolds numbers the

2007 ◽  
Vol 589 ◽  
pp. 57-81 ◽  
Author(s):  
G. GULITSKI ◽  
M. KHOLMYANSKY ◽  
W. KINZELBACH ◽  
B. LÜTHI ◽  
A. TSINOBER ◽  
...  

This is a report on a field experiment in an atmospheric surface layer at heights between 0.8 and 10m with the Taylor micro-scale Reynolds number in the range Reλ = 1.6−6.6 ×103. Explicit information is obtained on the full set of velocity and temperature derivatives both spatial and temporal, i.e. no use of Taylor hypothesis is made. The report consists of three parts. Part 1 is devoted to the description of facilities, methods and some general results. Certain results are similar to those reported before and give us confidence in both old and new data, since this is the first repetition of this kind of experiment at better data quality. Other results were not obtained before, the typical example being the so-called tear-drop R-Q plot and several others. Part 2 concerns accelerations and related matters. Part 3 is devoted to issues concerning temperature, with the emphasis on joint statistics of temperature and velocity derivatives. The results obtained in this work are similar to those obtained in experiments in laboratory turbulent grid flow and in direct numerical simulations of Navier–Stokes equations at much smaller Reynolds numbers Reλ ~ 102, and this similarity is not only qualitative, but to a large extent quantitative. This is true of such basic processes as enstrophy and strain production, geometrical statistics, the role of concentrated vorticity and strain, reduction of nonlinearity and non-local effects. The present experiments went far beyond the previous ones in two main respects. (i) All the data were obtained without invoking the Taylor hypothesis, and therefore a variety of results on fluid particle accelerations became possible. (ii) Simultaneous measurements of temperature and its gradients with the emphasis on joint statistics of temperature and velocity derivatives. These are reported in Parts 2 and 3.


2021 ◽  
Vol 930 ◽  
Author(s):  
Kartik P. Iyer ◽  
Katepalli R. Sreenivasan ◽  
P.K. Yeung

Using direct numerical simulations performed on periodic cubes of various sizes, the largest being $8192^3$ , we examine the nonlinear advection term in the Navier–Stokes equations generating fully developed turbulence. We find significant dissipation even in flow regions where nonlinearity is locally absent. With increasing Reynolds number, the Navier–Stokes dynamics amplifies the nonlinearity in a global sense. This nonlinear amplification with increasing Reynolds number renders the vortex stretching mechanism more intermittent, with the global suppression of nonlinearity, reported previously, restricted to low Reynolds numbers. In regions where vortex stretching is absent, the angle and the ratio between the convective vorticity and solenoidal advection in three-dimensional isotropic turbulence are statistically similar to those in the two-dimensional case, despite the fundamental differences between them.


2013 ◽  
Vol 721 ◽  
pp. 58-85 ◽  
Author(s):  
Kengo Deguchi ◽  
Philip Hall ◽  
Andrew Walton

AbstractThe recently understood relationship between high-Reynolds-number vortex–wave interaction theory and computationally generated self-sustaining processes provides a possible route to an understanding of some of the underlying structures of fully turbulent flows. Here vortex–wave interaction (VWI) theory is used in the long streamwise wavelength limit to continue the development found at order-one wavelengths by Hall & Sherwin (J. Fluid Mech., vol. 661, 2010, pp. 178–205). The asymptotic description given reduces the Navier–Stokes equations to the so-called boundary-region equations, for which we find equilibrium states describing the change in the VWI as the wavelength of the wave increases from $O(h)$ to $O(Rh)$, where $R$ is the Reynolds number and $2h$ is the depth of the channel. The reduced equations do not include the streamwise pressure gradient of the perturbation or the effect of streamwise diffusion of the wave–vortex states. The solutions we calculate have an asymptotic error proportional to ${R}^{- 2} $ when compared to the full Navier–Stokes equations. The results found correspond to the minimum drag configuration for VWI states and might therefore be of relevance to the control of turbulent flows. The key feature of the new states discussed here is the thickening of the critical layer structure associated with the wave part of the flow to completely fill the channel, so that the roll part of the flow is driven throughout the flow rather than as in Hall & Sherwin as a stress discontinuity across the critical layer. We identify a critical streamwise wavenumber scaling, which, when approached, causes the flow to localize and take on similarities with computationally generated or experimentally observed turbulent spots. In effect, the identification of this critical wavenumber for a given value of the assumed high Reynolds number fixes a minimum box length necessary for the emergence of localized structures. Whereas nonlinear equilibrium states of the Navier–Stokes equations are thought to form a backbone on which turbulent flows hang, our results suggest that the localized states found here might play a related role for turbulent spots.


1990 ◽  
Vol 220 ◽  
pp. 459-484 ◽  
Author(s):  
H. M. Badr ◽  
M. Coutanceau ◽  
S. C. R. Dennis ◽  
C. Ménard

The unsteady flow past a circular cylinder which starts translating and rotating impulsively from rest in a viscous fluid is investigated both theoretically and experimentally in the Reynolds number range 103 [les ] R [les ] 104 and for rotational to translational surface speed ratios between 0.5 and 3. The theoretical study is based on numerical solutions of the two-dimensional unsteady Navier–Stokes equations while the experimental investigation is based on visualization of the flow using very fine suspended particles. The object of the study is to examine the effect of increase of rotation on the flow structure. There is excellent agreement between the numerical and experimental results for all speed ratios considered, except in the case of the highest rotation rate. Here three-dimensional effects become more pronounced in the experiments and the laminar flow breaks down, while the calculated flow starts to approach a steady state. For lower rotation rates a periodic structure of vortex evolution and shedding develops in the calculations which is repeated exactly as time advances. Another feature of the calculations is the discrepancy in the lift and drag forces at high Reynolds numbers resulting from solving the boundary-layer limit of the equations of motion rather than the full Navier–Stokes equations. Typical results are given for selected values of the Reynolds number and rotation rate.


2015 ◽  
Vol 777 ◽  
pp. 219-244 ◽  
Author(s):  
Jesse T. Ault ◽  
Kevin K. Chen ◽  
Howard A. Stone

Direct numerical simulations were used to investigate the downstream decay of fully developed flow in a $180^{\circ }$ curved pipe that exits into a straight outlet. The flow is studied for a range of Reynolds numbers and pipe-to-curvature radius ratios. Velocity, pressure and vorticity fields are calculated to visualize the downstream decay process. Transition ‘decay’ lengths are calculated using the norm of the velocity perturbation from the Hagen–Poiseuille velocity profile, the wall-averaged shear stress, the integral of the magnitude of the vorticity, and the maximum value of the $Q$-criterion on a cross-section. Transition lengths to the fully developed Poiseuille distribution are found to have a linear dependence on the Reynolds number with no noticeable dependence on the pipe-to-curvature radius ratio, despite the flow’s dependence on both parameters. This linear dependence of Reynolds number on the transition length is explained by linearizing the Navier–Stokes equations about the Poiseuille flow, using the form of the fully developed Dean flow as an initial condition, and using appropriate scaling arguments. We extend our results by comparing this flow recovery downstream of a curved pipe to the flow recovery in the downstream outlets of a T-junction flow. Specifically, we compare the transition lengths between these flows and document how the transition lengths depend on the Reynolds number.


2013 ◽  
Vol 734 ◽  
pp. 275-316 ◽  
Author(s):  
Rashad Moarref ◽  
Ati S. Sharma ◽  
Joel A. Tropp ◽  
Beverley J. McKeon

AbstractWe study the Reynolds-number scaling and the geometric self-similarity of a gain-based, low-rank approximation to turbulent channel flows, determined by the resolvent formulation of McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382), in order to obtain a description of the streamwise turbulence intensity from direct consideration of the Navier–Stokes equations. Under this formulation, the velocity field is decomposed into propagating waves (with single streamwise and spanwise wavelengths and wave speed) whose wall-normal shapes are determined from the principal singular function of the corresponding resolvent operator. Using the accepted scalings of the mean velocity in wall-bounded turbulent flows, we establish that the resolvent operator admits three classes of wave parameters that induce universal behaviour with Reynolds number in the low-rank model, and which are consistent with scalings proposed throughout the wall turbulence literature. In addition, it is shown that a necessary condition for geometrically self-similar resolvent modes is the presence of a logarithmic turbulent mean velocity. Under the practical assumption that the mean velocity consists of a logarithmic region, we identify the scalings that constitute hierarchies of self-similar modes that are parameterized by the critical wall-normal location where the speed of the mode equals the local turbulent mean velocity. For the rank-1 model subject to broadband forcing, the integrated streamwise energy density takes a universal form which is consistent with the dominant near-wall turbulent motions. When the shape of the forcing is optimized to enforce matching with results from direct numerical simulations at low turbulent Reynolds numbers, further similarity appears. Representation of these weight functions using similarity laws enables prediction of the Reynolds number and wall-normal variations of the streamwise energy intensity at high Reynolds numbers (${Re}_{\tau } \approx 1{0}^{3} {\unicode{x2013}} 1{0}^{10} $). Results from this low-rank model of the Navier–Stokes equations compare favourably with experimental results in the literature.


2017 ◽  
Vol 25 (2) ◽  
Author(s):  
Xin He ◽  
Cornelis Vuik ◽  
Christiaan Klaij

Abstract The modified augmented Lagrangian preconditioner has attracted much attention in solving nondimensional Navier–Stokes equations discretized by the finite element method. In industrial applications the governing equations are often in dimensional form and discretized using the finite volume method. This paper assesses the capability of this preconditioner for dimensional Navier–Stokes equations in the context of the finite volume method. Two main contributions are made. First, this paper introduces a new dimensionless parameter that is involved in the modified augmented Lagrangian preconditioner. Second, with a number of academic test problems this paper reveals that the convergence of both nonlinear and linear iterations depend on this dimensionless parameter. A way to choose the optimal value of the dimensionless parameter is suggested and it is found that the optimal value is dependent of the Reynolds number, instead of the fluid’s properties, e.g., density and dynamic viscosity. The outcomes of this paper yield a potential rule to choose the optimal dimensionless parameter in practice, namely, correspondingly increasing with enlarging the Reynolds number.


Sign in / Sign up

Export Citation Format

Share Document