scholarly journals A FENE-P k–ε Viscoelastic Turbulence Model Valid up to High Drag Reduction without Friction Velocity Dependence

2020 ◽  
Vol 10 (22) ◽  
pp. 8140
Author(s):  
Michael McDermott ◽  
Pedro Resende ◽  
Thibaut Charpentier ◽  
Mark Wilson ◽  
Alexandre Afonso ◽  
...  

A viscoelastic turbulence model in a fully-developed drag reducing channel flow is improved, with turbulent eddies modelled under a k–ε representation, along with polymeric solutions described by the finitely extensible nonlinear elastic-Peterlin (FENE-P) constitutive model. The model performance is evaluated against a wide variety of direct numerical simulation data, described by different combinations of rheological parameters, which is able to predict all drag reduction (low, intermediate and high) regimes with good accuracy. Three main contributions are proposed: one with a simplified viscoelastic closure for the NLTij term (which accounts for the interactions between the fluctuating components of the conformation tensor and the velocity gradient tensor), by removing additional damping functions and reducing complexity compared with previous models; second through a reformulation for the closure of the viscoelastic destruction term, Eτp, which removes all friction velocity dependence; lastly by an improved modified damping function capable of predicting the reduction in the eddy viscosity and thus accurately capturing the turbulent kinetic energy throughout the channel. The main advantage is the capacity to predict all flow fields for low, intermediate and high friction Reynolds numbers, up to high drag reduction without friction velocity dependence.

2020 ◽  
Vol 142 (10) ◽  
Author(s):  
Xiao Xiong ◽  
Yan Zhang ◽  
Mohammad Azizur Rahman

Abstract Reynolds-averaged modeling is performed for polymer-induced drag reduction (DR) fluid at the fully developed turbulent regime in a concentric annulus by using the commercial code, ansys-fluent. The numerical approach adopted in this study relies on a modified k–ε–v2¯–f model to characterize the turbulence and the finitely extensible nonlinear elastic-Peterlin (FENE-P) constitutive model to represent the rheological behavior of the polymer solution. The near-wall axial velocity, Reynolds stress, and turbulent kinetic energy (TKE) budget near both walls of the annulus (fixed radius ratio of 0.4) are compared in detail at a constant Reynolds number (Re=10,587) and various rheological parameters (Weissenberg number We in the range of 1–7 and the maximum polymer elongation L = 30 and 100). Current simulation has predicted the redistributions of turbulent statistics in the annulus, where the two turbulent boundary layers (TBLs) of the DR flow differ more compared to those of its Newtonian counterpart. The difference is also found to be highly dependent on the rheological properties of the viscoelastic fluid.


2020 ◽  
Vol 36 (5) ◽  
pp. 691-698
Author(s):  
Che-Yu Lin ◽  
Chao-An Lin

ABSTRACTDirect numerical simulations have been applied to simulate flows with polymer additives. FENE-P (finite-extensible-nonlinear-elastic-Peterlin) dumbbell model solving for the conformation tensor is adopted to investigate the influence of the polymer on the flowfield. Boundary treatments of the conformation tensor on the flowfield are examined first, where boundary condition based on the linear extrapolation scheme provides more accurate results with second-order accurate error norms. Further simulations of the turbulent channel flow at different Weissenberg numbers are also conducted to investigate the influence on drag reduction. Drag reduction increases in tandem with the increase of Weissenberg number and the increase saturates at Weτ~200, where the drag reduction is close to the maximum drag reduction (MDR) limit. At the regime of y+ > 5, the viscous layer thickens with the increase of the Weissenberg number showing a departure from the traditional log-law profile, and the velocity profiles approach the MDR line at high Weissenberg number. The Reynolds stress decreases in tandem with the increase of Weτ, whereas the levels of laminar stress and polymer stress act adversely. However, as the Weissenberg number increases, the proportion of the laminar stress in the total stress increases, and this contributes to the drag reduction of the polymer flow.


2021 ◽  
Vol 90 ◽  
pp. 108799
Author(s):  
M. McDermott ◽  
P.R. Resende ◽  
M.C.T. Wilson ◽  
A.M. Afonso ◽  
D. Harbottle ◽  
...  

Author(s):  
Karsten Tawackolian ◽  
Martin Kriegel

AbstractThis study looks to find a suitable turbulence model for calculating pressure losses of ventilation components. In building ventilation, the most relevant Reynolds number range is between 3×104 and 6×105, depending on the duct dimensions and airflow rates. Pressure loss coefficients can increase considerably for some components at Reynolds numbers below 2×105. An initial survey of popular turbulence models was conducted for a selected test case of a bend with such a strong Reynolds number dependence. Most of the turbulence models failed in reproducing this dependence and predicted curve progressions that were too flat and only applicable for higher Reynolds numbers. Viscous effects near walls played an important role in the present simulations. In turbulence modelling, near-wall damping functions are used to account for this influence. A model that implements near-wall modelling is the lag elliptic blending k-ε model. This model gave reasonable predictions for pressure loss coefficients at lower Reynolds numbers. Another example is the low Reynolds number k-ε turbulence model of Wilcox (LRN). The modification uses damping functions and was initially developed for simulating profiles such as aircraft wings. It has not been widely used for internal flows such as air duct flows. Based on selected reference cases, the three closure coefficients of the LRN model were adapted in this work to simulate ventilation components. Improved predictions were obtained with new coefficients (LRNM model). This underlined that low Reynolds number effects are relevant in ventilation ductworks and give first insights for suitable turbulence models for this application. Both the lag elliptic blending model and the modified LRNM model predicted the pressure losses relatively well for the test case where the other tested models failed.


AIAA Journal ◽  
2001 ◽  
Vol 39 (12) ◽  
pp. 2402-2404 ◽  
Author(s):  
S. L. Woodruff ◽  
J. M. Seiner ◽  
M. Y. Hussaini ◽  
G. Erlebacher

Author(s):  
Stephen A Solovitz

Abstract Following volcanic eruptions, forecasters need accurate estimates of mass eruption rate (MER) to appropriately predict the downstream effects. Most analyses use simple correlations or models based on large eruptions at steady conditions, even though many volcanoes feature significant unsteadiness. To address this, a superposition model is developed based on a technique used for spray injection applications, which predicts plume height as a function of the time-varying exit velocity. This model can be inverted, providing estimates of MER using field observations of a plume. The model parameters are optimized using laboratory data for plumes with physically-relevant exit profiles and Reynolds numbers, resulting in predictions that agree to within 10% of measured exit velocities. The model performance is examined using a historic eruption from Stromboli with well-documented unsteadiness, again providing MER estimates of the correct order of magnitude. This method can provide a rapid alternative for real-time forecasting of small, unsteady eruptions.


2008 ◽  
Vol 18 (1) ◽  
pp. 12421-1-12421-5
Author(s):  
V. Mik ◽  
J. Myska ◽  
Z. Chara ◽  
P. Stern

AbstractEffectiveness of drag reduction by small addition of a surfactant in the turbulent flow of water depends on the structure and concentration of the additive, temperature of the solution and turbulence intensity, possible flow disturbance by a mechanical obstacle and the content of ions in water, but also on the age of the surfactant solution. We show how important aging effects are in connection with total surfactant concentration, in particular how rheological parameters of the drag reducing solution change with time.


2019 ◽  
Vol 874 ◽  
pp. 699-719 ◽  
Author(s):  
Jose M. Lopez ◽  
George H. Choueiri ◽  
Björn Hof

Polymer additives can substantially reduce the drag of turbulent flows and the upper limit, the so-called state of ‘maximum drag reduction’ (MDR), is to a good approximation independent of the type of polymer and solvent used. Until recently, the consensus was that, in this limit, flows are in a marginal state where only a minimal level of turbulence activity persists. Observations in direct numerical simulations at low Reynolds numbers ($Re$) using minimal sized channels appeared to support this view and reported long ‘hibernation’ periods where turbulence is marginalized. In simulations of pipe flow at $Re$ near transition we find that, indeed, with increasing Weissenberg number ($Wi$), turbulence expresses long periods of hibernation if the domain size is small. However, with increasing pipe length, the temporal hibernation continuously alters to spatio-temporal intermittency and here the flow consists of turbulent puffs surrounded by laminar flow. Moreover, upon an increase in $Wi$, the flow fully relaminarizes, in agreement with recent experiments. At even larger $Wi$, a different instability is encountered causing a drag increase towards MDR. Our findings hence link earlier minimal flow unit simulations with recent experiments and confirm that the addition of polymers initially suppresses Newtonian turbulence and leads to a reverse transition. The MDR state on the other hand results at these low$Re$ from a separate instability and the underlying dynamics corresponds to the recently proposed state of elasto-inertial turbulence.


Water ◽  
2021 ◽  
Vol 14 (1) ◽  
pp. 36
Author(s):  
Jiyang Qi ◽  
Yue Qi ◽  
Qunyan Chen ◽  
Fei Yan

In this study, the drag reduction effect is studied for a cylinder with different V-groove depths on its surface using a k-ω/SST (Shear Stress Transport) turbulence model of computational fluid dynamics (CFD), while a particle image velocimetry (PIV) system is employed to analyze the wake characteristics for a smooth cylinder and a cylinder with different V-groove depths on its surface at different Reynolds numbers. The study focuses on the characteristics of the different V-groove depths on lift coefficient, drag coefficient, the velocity distribution of flow field, pressure coefficient, vortex shedding, and vortex structure. In comparison with a smooth cylinder, the lift coefficient and drag coefficient can be reduced for a cylinder with different V-groove depths on its surface, and the maximum reduction rates of lift coefficient and drag coefficient are about 34.4% and 16%, respectively. Otherwise, the vortex structure presents a complete symmetry for the smooth cylinder, however, the symmetry of the vortex structure becomes insignificant for the V-shaped groove structure with different depths. This is also an important reason for the drag reduction effect of a cylinder with a V-groove surface.


Sign in / Sign up

Export Citation Format

Share Document