Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDoes the choice of the forcing term affect flow statistics in DNS of turbulent channel flow?
128
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We seek possible statistical consequences of the way a forcing term is added to the Navier--Stokes equations in the Direct Numerical Simulation (DNS) of incompressible channel flow. Simulations driven by constant flow rate, constant pressure gradient and constant power input are used to build large databases, and in particular to store the complete temporal trace of the wall-shear stress for later analysis. As these approaches correspond to different dynamical systems, it can in principle be envisaged that these differences are reflect by certain statistics of the turbulent flow field. The instantaneous realizations of the flow in the various simulations are obviously different, but, as expected, the usual one-point, one-time statistics do not show any appreciable difference. However, the PDF for the fluctuations of the streamwise component of wall friction reveals that the simulation with constant flow rate presents lower probabilities for extreme events of large positive friction. The low probability value of such events explains their negligible contribution to the commonly computed statistics; however, the very existence of a difference in the PDF demonstrates that the forcing term is not entirely uninfluential. Other statistics for wall-based quantities (the two components of friction and pressure) are examined; in particular spatio-temporal autocorrelations show small differences at large temporal separations, where unfortunately the residual statistical uncertainty is still of the same order of the observed difference. Hence we suggest that the specific choice of the forcing term does not produce important statistical consequences, unless one is interested in the strongest events of high wall friction, that are underestimated by a simulation run at constant flow rate.
rate research
Read More
We present direct numerical simulations of turbulent channel flow with passive Lagrangian polymers. To understand the polymer behavior we investigate the behavior of infinitesimal line elements and calculate the probability distribution function (PDF) of finite-time Lyapunov exponents and from them the corresponding Cramers function for the channel flow. We study the statistics of polymer elongation for both the Oldroyd-B model (for Weissenberg number $Wi <1$) and the FENE model. We use the location of the minima of the Cramers function to define the Weissenberg number precisely such that we observe coil-stretch transition at $Wiapprox1$. We find agreement with earlier analytical predictions for PDF of polymer extensions made by Balkovsky, Fouxon and Lebedev [Phys. Rev. Lett., 84, 4765 (2000).] for linear polymers (Oldroyd-B model) with $Wi<1$ and by Chertkov [Phys. Rev. Lett., 84, 4761 (2000).] for nonlinear FENE-P model of polymers. For $Wi>1$ (FENE model) the polymer are significantly more stretched near the wall than at the center of the channel where the flow is closer to homogenous isotropic turbulence. Furthermore near the wall the polymers show a strong tendency to orient along the stream-wise direction of the flow but near the centerline the statistics of orientation of the polymers is consistent with analogous results obtained recently in homogeneous and isotropic flows.
The turbulent flow in an infinitely extended plane channel is analysed by solving the Navier-Stokes equations with a DNS approach. Solutions are obtained in a numerical solution domain of finite size in the streamwise as well as in the lateral direction setting periodic boundary conditions in both directions. Their impact on large scale structures in the turbulent flow field is analysed carefully in order to avoid their suppression. When this is done appropriately well known stripe patterns in these flows can be observed and analysed especially with respect to their relative motion compared to the mean flow velocity. Various details of this stripe pattern dominated velocity field are shown. Also global parameters like the friction factor in the flow field and the Nusselt number in the temperature field are determined based on the statistics of the flow and temperature data in a very large time period that guarantees fully developed turbulent flow and heat transfer.
An analysis of the statistics of the non-linear terms in resolvent analysis is performed in this work for turbulent Couette flow at low Reynolds number. Data from a direct numerical simulation of a minimal flow unit, at Reynolds number 400, is post-processed using Fourier analysis in both time and space, leading to the covariance matrix of the velocity. From the same data, we computed the non-linear terms of the Navier-Stokes equations (treated as forcing in the present formulation), which allowed us to compute the covariance matrix of the forcing for this case. The two covariances are related exactly by the resolvent operator; based on this, we explore the recovery of the velocity statistics from the statistics of the forcing as a function of the components of the forcing term. This is carried out for the dominant structures in this flow, which participate in the self-sustaining cycle of turbulence: (i) streamwise vortices and streaks, and (ii) spanwise coherent fluctuations of spanwise velocity. The present results show a dominance by four of the non-linear terms for the prediction of the full statistics of streamwise vortices and streaks; a single term is seen to be dominant for spanwise motions. A relevant feature observed in these cases is that forcing terms have significant coherence in space; moreover, different forcing components are also coherent between them. This leads to constructive and destructive interferences that greatly modify the flow response, and should thus be accounted for in modelling work.
The transitional regime of plane channel flow is investigated {above} the transitional point below which turbulence is not sustained, using direct numerical simulation in large domains. Statistics of laminar-turbulent spatio-temporal intermittency are reported. The geometry of the pattern is first characterized, including statistics for the angles of the laminar-turbulent stripes observed in this regime, with a comparison to experiments. High-order statistics of the local and instantaneous bulk velocity, wall shear stress and turbulent kinetic energy are then provided. The distributions of the two former quantities have non-trivial shapes, characterized by a large kurtosis and/or skewness. Interestingly, we observe a strong linear correlation between their kurtosis and their skewness squared, which is usually reported at much higher Reynolds number in the fully turbulent regime.
The cross-spectral density (CSD) of the non-linear forcing in resolvent analyses is here quantified for the first time for turbulent channel flows. Direct numerical simulations (DNS) at $Re_{tau} =179$ and $Re_{tau} =543$ are performed. The CSDs are computed for highly energetic structures typical of buffer-layer and large-scale motions, for different temporal frequencies. The CSD of the non-linear forcing is shown not to be uncorrelated (white) in space, which implies the forcing is structured. Since the non-linear forcing is non-solenoidal by construction and the velocity of an incompressible flow is affected only by the solenoidal part of the forcing, this solenoidal part is evaluated. It is shown that the solenoidal part of the non-linear forcing is the combination of oblique streamwise vortices and a streamwise component which counteract each other, as in a destructive interference. It is shown that a rank-2 approximation of the forcing, with only the most energetic SPOD (spectral proper orthogonal decomposition) modes, leads to the bulk of the response. The projections of the non-linear forcing onto the right-singular vectors of the resolvent are evaluated. The left-singular vectors of the resolvent associated with very low-magnitude singular values are non-negligible since the non-linear forcing term has a non-negligible projection onto the linear sub-optimals of resolvent analysis. The same projections are computed when the forcing is modelled with an eddy-viscosity approach. It is clarified that this modelling improves the accuracy of the prediction since the projections are closer to those associated with the non-linear forcing from DNS data.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
