DNS study of decaying homogeneous isotropic turbulence with polymer additives

2010 ◽  
Vol 665 ◽  
pp. 334-356 ◽  
Author(s):  
W.-H. CAI ◽  
F.-C. LI ◽  
H.-N. ZHANG

In order to investigate the turbulent drag reduction phenomenon and understand its mechanism, direct numerical simulation (DNS) was carried out on decaying homogeneous isotropic turbulence (DHIT) with and without polymer additives. We explored the polymer effect on DHIT from the energetic viewpoint, i.e. the decay of the total turbulent kinetic energy and energy distribution at each scale in Fourier space and from the phenomenological viewpoint, i.e. the alterations of vortex structures, the enstrophy and the strain. It was obtained that in DHIT with polymer additives the decay of the turbulent kinetic energy is faster than that in the Newtonian fluid case and a modification of the turbulent kinetic energy transfer process for the Newtonian fluid flow is observed due to the release of the polymer elastic energy into flow structures at certain small scales. Besides, we deduced the transport equations of the enstrophy and the strain, respectively, for DHIT with polymer additives. Based on the analyses of these transport equations, it was found that polymer additives depress both the enstrophy and the strain in DHIT as compared to the Newtonian fluid case, indicating the inhibition effect on small-scale vortex structures and turbulence intensity by polymers.

2017 ◽  
Vol 95 (12) ◽  
pp. 1271-1277 ◽  
Author(s):  
Yue Wang ◽  
Wei-Hua Cai ◽  
Xin Zheng ◽  
Hong-Na Zhang ◽  
Feng-Chen Li

In this paper, to study the viscoelastic effect on isotropic turbulence without wall effects, a two oscillating grid turbulence is built to investigate this phenomenon using particle image velocimetry. In the experiments, the classical drag-reducing additives are chosen: polyacrylamide (PAM) and cetyltrimethyl ammonium chloride (CTAC), which have shown remarkable drag-reducing effect in wall-bounded turbulent flows. The results show that the existence of drag-reducing additives makes velocity field more anisotropic and reduces turbulent kinetic energy. We propose an intuitive and natural definition for a reduction rate of turbulent kinetic energy to show viscoelastic effect. It suggests that there exists a critical concentration for the reduction rate of turbulent kinetic energy in the CTAC solution case. Also, the small-scale vortex structures are inhibited, which suggests the drag-reducing mechanism in grid turbulence without wall effect.


2014 ◽  
Vol 6 ◽  
pp. 193731 ◽  
Author(s):  
Jie Yan ◽  
Nan Gui ◽  
Gongnan Xie ◽  
Jinsen Gao

Two parallel swirling/rotating jets with a distance between them are termed biswirling jets here, which have important and complicated vortex structures different from the single swirling jet due to the negligible vortex-vortex interactions. The visualization of vortex-vortex interaction between the biswirling jets is accomplished by using direct numerical simulation. The evolution of vortex structures of the biswirling jets is found rather complicated. The turbulent kinetic energy and turbulence dissipation in the central convergence region are augmented locally and rather strongly. The modulation of turbulent kinetic energy by jet-jet interaction upon different scales of vortices is dominated by the swirling levels and the distance between the jets. The turbulent kinetic energy upon intermediate and small scale vortices in bijets with not very high swirling level and at a very close distance is smaller than that in single swirling jets, whereas the opposite is true under a far distance, and so forth.


Author(s):  
Lu Wang ◽  
Yue Wang ◽  
Wei-Hua Cai ◽  
Feng-Chen Li

The drag-reducing phenomenon in forced homogeneous isotropic turbulence (FHIT) with polymer additives was realized in large eddy simulation (LES) results, which causes the variation of turbulent characteristics. Study on intermittency of turbulence is of great importance in investigating turbulent drag-reduction mechanism, because the intermittency has close relationship with coherent structures (CSs) and the transfer of energy in turbulent flows. In the present work, the influences of polymers on intermittency in FHIT were analyzed in detail by extracting CSs, researching the flatness factor based on high-order correlation function of velocity derivative and wavelet transform, surveying local intermittence measure, and discussing four rotational invariants consisting of velocity-strain tensor and vorticity tensor. From the viewpoint of the results, it can be perceived that the intermittency occurs in both the Newtonian fluid and polymer solution flows; moreover polymer additives behave inhibitive effect on the intermittency in turbulent drag-reducing flows.


2019 ◽  
Vol 4 (10) ◽  
Author(s):  
Mohamad Ibrahim Cheikh ◽  
James Chen ◽  
Mingjun Wei

2020 ◽  
Author(s):  
Wojciech W. Grabowski ◽  
Lois Thomas

Abstract. Increase of the spectral width of initially monodisperse population of cloud droplets in homogeneous isotropic turbulence is investigated applying a finite-difference fluid flow model combined with either Eulerian bin microphysics or Lagrangian particle-based scheme. The turbulence is forced applying a variant of the so-called linear forcing method that maintains the mean turbulent kinetic energy (TKE) and the TKE partitioning between velocity components. The latter is important for maintaining the quasi-steady forcing of the supersaturation fluctuations that drive the increase of the spectral width. We apply a large computational domain, 643 m3, one of the domains considered in Thomas et al. (2020). The simulations apply 1 m grid length and are in the spirit of the implicit large eddy simulation (ILES), that is, with explicit small-scale dissipation provided by the model numerics. This is in contrast to the scaled-up direct numerical simulation (DNS) applied in Thomas et al. (2020). Two TKE intensities and three different droplet concentrations are considered. Analytic solutions derived in Sardina et al. (2015), valid for the case when the turbulence time scale is much larger than the droplet phase relaxation time scale, are used to guide the comparison between the two microphysics simulation techniques. The Lagrangian approach reproduces the scalings relatively well. Representing the spectral width increase in time is more challenging for the bin microphysics because appropriately high resolution in the bin space is needed. The bin width of 0.5 μm is only sufficient for the lowest droplet concentration, 26 cm−3. For the highest droplet concentration, 650 cm−3, even an order of magnitude smaller bin size is not sufficient. The scalings are not expected to be valid for the lowest droplet concentration and the high TKE case, and the two microphysics schemes represent similar departures. Finally, because the fluid flow is the same for all simulations featuring either low or high TKE, one can compare point-by-point simulation results. Such a comparison shows very close temperature and water vapor point-by-point values across the computational domain, and larger differences between simulated mean droplet radii and spectral width. The latter are explained by fundamental differences in the two simulation methodologies, numerical diffusion in the Eulerian bin approach and relatively small number of Lagrangian particles that are used in the particle-based microphysics.


1991 ◽  
pp. 422-434 ◽  
Author(s):  
J. C. R. Hunt ◽  
J. C. H. Fung ◽  
N. A. Malik ◽  
R. J. Perkins ◽  
J. C. Vassilicos ◽  
...  

2012 ◽  
Vol 7 (1) ◽  
pp. 53-69
Author(s):  
Vladimir Dulin ◽  
Yuriy Kozorezov ◽  
Dmitriy Markovich

The present paper reports PIV (Particle Image Velocimetry) measurements of turbulent velocity fluctuations statistics in development region of an axisymmetric free jet (Re = 28 000). To minimize measurement uncertainty, adaptive calibration, image processing and data post-processing algorithms were utilized. On the basis of theoretical analysis and direct measurements, the paper discusses effect of PIV spatial resolution on measured statistical characteristics of turbulent fluctuations. Underestimation of the second-order moments of velocity derivatives and of the turbulent kinetic energy dissipation rate due to a finite size of PIV interrogation area and finite thickness of laser sheet was analyzed from model spectra of turbulent velocity fluctuations. The results are in a good agreement with the measured experimental data. The paper also describes performance of possible ways to account for unresolved small-scale velocity fluctuations in PIV measurements of the dissipation rate. In particular, a turbulent viscosity model can be efficiently used to account for the unresolved pulsations in a free turbulent flow


Author(s):  
Y. Kagawa ◽  
B. Yu ◽  
Y. Kawaguchi

For the purpose of elucidating the mechanism of drag reduction by additives and finding a way to judge optimum drag-reducing additives through a simple rheological test, we performed DNS analysis of viscoelastic fluid turbulent flow in a two-dimensional channel. In this calculation, we employed the Giesekus constitutive equation to model the interaction between water-soluble polymer, or the elastic micellar network structure, and solvent. We calculated the fluid flow by varying the rheological parameters of the model. We examined the turbulent kinetic energy budget and studied the “viscoelastic contribution” term in the budget equation for turbulent intensity, which is not apparent in normal Newtonian fluid turbulence. Viscoelastic contribution has a characteristic effect on viscoelastic fluid turbulence. We concluded that the viscoelastic contribution plays a major role in turbulent frictional drag reduction. Dissipation and viscoelastic contribution serve as a key factor of turbulent kinetic energy loss in most areas of the channel. From the visualization of local and instantaneous eddy behavior, we discussed the relationship between viscoelastic contribution, elastic energy and turbulent production. We found that viscoelastic contribution serves as a direct local source of turbulent production, and that energy is stored in the elasticity.


Sign in / Sign up

Export Citation Format

Share Document