On the presence of spectral shortcut in the energy budget of an asymmetric jet–wake flow in a forward-curved centrifugal turbomachine as deduced from SPIV measurements

Viscous dissipation and its contribution to turbulent kinetic energy (TKE) budget are investigated in the asymmetric jet–wake flow of a forward-curved centrifugal turbomachine. Single-plane three-dimensional turbulent data are obtained using stereoscopic particle image velocimetry (SPIV). Viscous dissipation is indirectly estimated from subgrid-scale (SGS) dissipation (SGS energy flux) by filtering velocity field using a top-hat filter. The filter scale should be within the inertial sub-range and this is ensured by spectral analysis of the measured field. Reduction of turbulent energy flux for smaller filter scales plus underestimation of viscous dissipation as compared with other TKE terms both suggest the presence of spectral shortcut. This bypass energy transfer (from intermediate scales towards dissipative scales) works in parallel with direct SGS energy transfer and affects the classical energy cascade. Analysis of TKE budget in the rotor exit region shows significant radial/circumferential variations in the contributing terms. These variations are mainly due to jet–wake–volute interactions, circumferential asymmetry of volute area and expansion of flow toward the fan outlet.     

Assessment of subgrid-scale stress statistics in non-premixed turbulent wall-jet flames

We investigate the heat-release effects on the characteristics of the subgrid-scale (SGS) stress tensor and SGS dissipation of kinetic energy and enstrophy. Direct numerical simulation data of a non-premixed reacting turbulent wall-jet flow with and without substantial heat release is employed for the analysis. This study comprises, among others, an analysis of the eigenvalues of the resolved strain rate and SGS stress tensors, to identify the heat-release effects on their topology. An assessment of the alignment between the eigenvectors corresponding to the largest eigenvalues of these two tensors is also given to provide further information for modelling of the SGS stress tensor. To find out the heat-release effects on the dynamics of the turbulent kinetic energy and enstrophy dissipation, probability density functions (PDFs) and mean values are analysed. The mean SGS shear stress and turbulent kinetic energy both slightly increase in the buffer layer and substantially decrease further away from the wall, due to the heat-release effects. Contrary to the kinetic energy, heat release decreases the mean SGS dissipation of enstrophy in the near-wall region. Moreover, differences in the shapes of the PDFs between the isothermal and exothermic cases indicate changes in the intermittency level of both SGS dissipations. Heat release also increases the SGS stress anisotropy in the near-wall region. Although, the structure of the mean resolved strain-rate tensor only marginally differs between the isothermal and exothermic cases in the near-wall region, substantial differences are observed in the jet area, where compressibility effects are important and heat-release effects are found to promote compression states. The differences in the relative alignment between the SGS stress and resolved strain-rate tensors in the isothermal and exothermic cases are discussed in connection with the differences in the SGS dissipation of kinetic energy.     

Effects of variable specific heat on energy transfer in a high-temperature supersonic channel flow

An energy transfer mechanism in high-temperature supersonic turbulent flow for variable specific heat (VSH) condition through turbulent kinetic energy (TKE), mean kinetic energy (MKE), turbulent internal energy (TIE) and mean internal energy (MIE) is proposed. The similarities of energy budgets between VSH and constant specific heat (CSH) conditions are investigated by introducing a vibrational energy excited degree and considering the effects of fluctuating specific heat. Direct numerical simulation (DNS) of temporally evolving high-temperature supersonic turbulent channel flow is conducted at Mach number 3.0 and Reynolds number 4800 combined with a constant dimensional wall temperature 1192.60 K for VSH and CSH conditions to validate the proposed energy transfer mechanism. The differences between the terms in the two kinetic energy budgets for VSH and CSH conditions are small; however, the magnitude of molecular diffusion term for VSH condition is significantly smaller than that for CSH condition. The non-negligible energy transfer is obtained after neglecting several small terms of diffusion, dissipation and compressibility related. The non-negligible energy transfer involving TIE includes three processes, in which energy can be gained from TKE and MIE and lost to MIE. The same non-negligible energy transfer through TKE, MKE and MIE is observed for both the conditions.     

Similarity solutions for boundary layer flow and heat transfer of a FENE-P fluid with thermal radiation

This Letter presents a numerical study of the flow and heat transfer of an incompressible FENE-P fluid over a non-isothermal surface. The governing partial differential equations are converted into ordinary differential equations by a similarity transformation. The effects of the thermal radiation are considered in the energy equation, and the variations of dimensionless surface temperature and dimensionless surface temperature gradient, as well as the heat transfer characteristics with various physical parameters are graphed and tabulated. Two cases are studied, namely, (i) the sheet with prescribed surface temperature (PST case) and (ii) the sheet with prescribed heat flux (PHF case). Moreover, the mechanical characteristics of the corresponding flow are also presented.     

Influence of polymer additives on turbulent energy cascading in forced homogeneous isotropic turbulence studied by direct numerical simulations

Direct numerical simulations(DNS) were performed for the forced homogeneous isotropic turbulence(FHIT) with/without polymer additives in order to elaborate the characteristics of the turbulent energy cascading influenced by drag-reducing effects.The finite elastic non-linear extensibility-Peterlin model(FENE-P) was used as the conformation tensor equation for the viscoelastic polymer solution.Detailed analyses of DNS data were carried out in this paper for the turbulence scaling law and the topological dynamics of FHIT as well as the important turbulent parameters,including turbulent kinetic energy spectra,enstrophy and strain,velocity structure function,small-scale intermittency,etc.A natural and straightforward definition for the drag reduction rate was also proposed for the drag-reducing FHIT based on the decrease degree of the turbulent kinetic energy.It was found that the turbulent energy cascading in the FHIT was greatly modified by the drag-reducing polymer additives.The enstrophy and the strain fields in the FHIT of the polymer solution were remarkably weakened as compared with their Newtonian counterparts.The small-scale vortices and the small-scale intermittency were all inhibited by the viscoelastic effects in the FHIT of the polymer solution.However,the scaling law in a fashion of extended self-similarity for the FHIT of the polymer solution,within the presently simulated range of Weissenberg numbers,had no distinct differences compared with that of the Newtonian fluid case.     

Multi-scale properties of large eddy simulations: correlations between resolved-scale velocity-field increments and subgrid-scale quantities

We provide analytical and numerical results concerning multi-scale correlations between the resolved velocity field and the subgrid-scale (SGS) stress-tensor in large eddy simulations (LES). Following previous studies for Navier–Stokes equations, we derive the exact hierarchy of LES equations governing the spatio-temporal evolution of velocity structure functions of any order. The aim is to assess the influence of the subgrid model on the inertial range intermittency. We provide a series of predictions, within the multifractal theory, for the scaling of correlation involving the SGS stress and we compare them against numerical results from high-resolution Smagorinsky LES and from a-priori filtered data generated from direct numerical simulations (DNS). We find that LES data generally agree very well with filtered DNS results and with the multifractal prediction for all leading terms in the balance equations. Discrepancies are measured for some of the sub-leading terms involving cross-correlation between resolved velocity increments and the SGS tensor or the SGS energy transfer, suggesting that there must be room to improve the SGS modelisation to further extend the inertial range properties for any fixed LES resolution.     

充模过程中熔接痕的改进光滑粒子动力学方法模拟与预测

提出了一种黏弹性流体的改进光滑粒子动力学(SPH)方法以试探性地模拟和预测黏弹性FENE-P熔体充模过程中熔接痕的形态演化. 首先基于SPH方法建立了聚合物流动的宏微观耦合模型, 同时提出了黏弹性流体的改进SPH离散格式. 随后, 通过模拟一些基准算例验证了改进的SPH方法模拟聚合物宏微观耦合问题的有效性及收敛性, 以及所提出的黏弹性温度模型的有效性. 最后, 模拟了环型腔内的充模过程, 试探性地展示了充模过程中微观分子的变形过程. 同时采用顺序热流道技术模拟了多浇口C形腔内的充模过程, 并与其他数值结果做比较. 数值结果表明: 对于大制件多浇口充模过程, 顺序热流道技术能够改善甚至消除充模过程中的熔接痕.     

A new mixed subgrid-scale model for large eddy simulation of turbulent drag-reducing flows of viscoelastic fluids

A mixed subgrid-scale(SGS) model based on coherent structures and temporal approximate deconvolution(MCT) is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation(LES) of turbulent dragreducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence(FHIT) with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation(DNS) results. Compared with the LES results using the temporal approximate deconvolution model(TADM) for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number.For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives.     

A mixing controlled direct chemistry (MCDC) model for diesel engine combustion modelling using large eddy simulation

A mixing controlled direct chemistry (MCDC) combustion model with sub-grid scale (SGS) mixing effects and chemical kinetics has been evaluated for Large Eddy Simulation (LES) of diesel engine combustion. The mixing effect is modelled by a mixing timescale based on mixture fraction variance and sub-grid scalar dissipation rate. The SGS scalar dissipation rate is modelled using a similarity term and a scaling factor from the analysis of Direct Numerical Simulation (DNS) data. The chemical reaction progress is estimated from a kinetic timescale based on local internal energy change rate and equilibrium state internal energy. An optical research engine operating at conventional operating conditions and Low Temperature Combustion (LTC) conditions was used for evaluation of the combustion model. From the simulation results, the effect of SGS scalar mixing is evaluated at different stages of combustion. In the context of LES, the new approach provides improved engine modelling results compared to the Direct Chemistry Solver (DCS) combustion model.     

Large eddy simulation of orientation and rotation of ellipsoidal particles in isotropic turbulent flows

The rotational motion and orientational distribution of ellipsoidal particles in turbulent flows are of significance in environmental and engineering applications. Whereas the translational motion of an ellipsoidal particle is controlled by the turbulent motions at large scales, its rotational motion is determined by the fluid velocity gradient tensor at small scales, which raises a challenge when predicting the rotational dispersion of ellipsoidal particles using large eddy simulation (LES) method due to the lack of subgrid scale (SGS) fluid motions. We report the effects of the SGS fluid motions on the orientational and rotational statistics, such as the alignment between the long axis of ellipsoidal particles and the vorticity, the mean rotational energy at various aspect ratios against those obtained with direct numerical simulation (DNS) and filtered DNS. The performances of a stochastic differential equation (SDE) model for the SGS velocity gradient seen by the particles and the approximate deconvolution method (ADM) for LES are investigated. It is found that the missing SGS fluid motions in LES flow fields have significant effects on the rotational statistics of ellipsoidal particles. Alignment between the particles and the vorticity is weakened; and the rotational energy of the particles is reduced in LES. The SGS-SDE model leads to a large error in predicting the alignment between the particles and the vorticity and over-predicts the rotational energy of rod-like particles. The ADM significantly improves the rotational energy prediction of particles in LES.     

Elasticity of polymer solutions in Couette flow measured by fluorescence resonance energy transfer (FRET)

Polymer solutions are complex fluids that show elasticity and deformation in response to shear flows. A fluorescence resonance energy transfer (FRET) technique has been applied to measure the end-to-end distances of individual polymer molecules in Couette flow, using end-tagged reversible-addition fragmentation chain transfer (RAFT) polymerised poly(methyl methacrylate) (PMMA). Real-time rheofluorescence measurements on these polymers in solution above the critical overlap concentration are reported at several shear rates. The PMMA in Couette flow shows a systematic decrease in fluorescence, corresponding to a reduction in end-to-end distance of the polymer molecules with shear exposure. Full reversibility of the fluorescence signal is observed after the cessation of shear. These results show that polymer solution elasticity arises from compressive deformation of the polymer molecules in Couette flow. At polymer concentrations above the critical overlap, the polymer molecules are restricted by their neighbours and the net hydrodynamic forces are compressive rather than extensive.     

Generalised scale-by-scale energy-budget equations and large-eddy simulations of anisotropic scalar turbulence at various Schmidt numbers

A necessary condition for the accurate prediction of turbulent flows using large-eddy simulation (LES) is the correct representation of energy transfer between the different scales of turbulence in the LES. For scalar turbulence, transfer of energy between turbulent length scales is described by a transport equation for the second moment of the scalar increment. For homogeneous isotropic turbulence, the underlying equation is the well-known Yaglom equation. In the present work, we study the turbulent mixing of a passive scalar with an imposed mean gradient by homogeneous isotropic turbulence. Both direct numerical simulations (DNS) and LES are performed for this configuration at various Schmidt numbers, ranging from 0.11 to 5.56. As the assumptions made in the derivation of the Yaglom equation are violated for the case considered here, a generalised Yaglom equation accounting for anisotropic effects, induced by the mean gradient, is derived in this work. This equation can be interpreted as a scale-by-scale energy-budget equation, as it relates at a certain scale r terms representing the production, turbulent transport, diffusive transport and dissipation of scalar energy. The equation is evaluated for the conducted DNS, followed by a discussion of physical effects present at different scales for various Schmidt numbers. For an analysis of the energy transfer in LES, a generalised Yaglom equation for the second moment of the filtered scalar increment is derived. In this equation, new terms appear due to the interaction between resolved and unresolved scales. In an a-priori test, this filtered energy-budget equation is evaluated by means of explicitly filtered DNS data. In addition, LES calculations of the same configuration are performed, and the energy budget as well as the different terms are thereby analysed in an a-posteriori test. It is shown that LES using an eddy viscosity model is able to fulfil the generalised filtered Yaglom equation for the present configuration. Further, the dependence of the terms appearing in the filtered energy-budget equation on varying Schmidt numbers is discussed.     

Evolution of surfaces and mechanisms of contact electrification between metals and polymers

Contact electrification (CE) is a pretty common phenomenon, but still is poorly understood. The long-standing controversy over the mechanisms of CE related to polymers is particularly intense due to their complexity. In this paper, the CE between metals and polymers is systematically studied, which shows the evolution of surfaces is accompanied by variations of CE outputs. The variations of CE charge quantity are closely related to the creep and deformation of the polymer and metal surfaces. Then the relationship between CE and polymer structures is put forward, which is essentially determined by the electronegativity of elements and the functional groups in the polymers. The effects of load and contact frequency on the CE process and outputs are also investigated, indicating the increase of CE charge quantity with load and frequency. Material transfer from polymer to metal is observed during CE while electrons transfer from metal to polymer, both of which are believed to have an influence on each other. The findings advance our understanding of the mechanism of CE between metal and polymers, and provides insights into the performance of CE-based application in various conditions, which sheds light on the design and optimization of CE-based energy harvest and self-powered sensing devices.     

A-priori evaluations of subgrid-scale terms for large-eddy simulation of compressible turbulent flows

The subgrid-scale terms for different formulations of the energy equation are evaluated from a-priori tests using the direct numerical simulation (DNS) data of a compressible mixing layer at a moderate Mach number of M = 0.65. To extend the generality of the results, the simulations were performed with three different initial conditions for the velocity fields. To examine the impact of strong temperature variations on the subgrid scales, a non-isothermal mixing layer with lower to upper free-stream temperature ratio of 3 is also considered. For cold simulations, with equal free-stream temperatures, the total energy equation is shown to be the best choice in view of the accuracy and the subgrid-scale modelling requirements. For hot simulations, with the free-stream temperature ratio equal to 3, the total enthalpy equation is found to be the best formulation for the energy equation. Furthermore, it is shown that the subgrid-scale pressure dilatation term, which has been largely neglected so far, is of the same order of the subgrid-scale heat flux. Based on the present results, the contribution of the subgrid-scale pressure dilatation can be up to 46% of the total sugbrid-scale activity. Moreover, the time evolutions of the volume-average mean kinetic energy, turbulent kinetic energy, production, dissipation, and pressure dilatation terms are considered. Unlike the subgrid-scale pressure dilatation term, the volume-average pressure dilatation terms are negligible, and compressibility does not affect the large-scale evolutions of the mean and turbulent kinetic energies.     

Large-eddy simulations of a forced homogeneous isotropic turbulence with polymer additives

Large-eddy simulations （LES） based on the temporal approximate deconvolution model were performed for a forced homogeneous isotropic turbulence （FHIT） with polymer additives at moderate Taylor Reynolds number. Finitely extensible nonlinear elastic in the Peterlin approximation model was adopted as the constitutive equation for the filtered conformation tensor of the polymer molecules. The LES results were verified through comparisons with the direct numerical simulation results. Using the LES database of the FHIT in the Newtonian fluid and the polymer solution flows, the polymer effects on some important parameters such as strain, vorticity, drag reduction, and so forth were studied. By extracting the vortex structures and exploring the flatness factor through a high-order correlation function of velocity derivative and wavelet analysis, it can be found that the small-scale vortex structures and small-scale intermittency in the FHIT are all inhibited due to the existence of the polymers. The extended self-similarity scaling law in the polymer solution flow shows no apparent difference from that in the Newtonian fluid flow at the currently simulated ranges of Reynolds and Weissenberg numbers.     

Characteristics and generation of elastic turbulence in a three-dimensional parallel plate channel using direct numerical simulation

Direct numerical simulations(DNSs) of purely elastic turbulence in rectilinear shear flows in a three-dimensional(3D) parallel plate channel were carried out,by which numerical databases were established.Based on the numerical databases,the present paper analyzed the structural and statistical characteristics of the elastic turbulence including flow patterns,the wall effect on the turbulent kinetic energy spectrum,and the local relationship between the flow motion and the microstructures' behavior.Moreover,to address the underlying physical mechanism of elastic turbulence,its generation was presented in terms of the global energy budget.The results showed that the flow structures in elastic turbulence were 3D with spatial scales on the order of the geometrical characteristic length,and vortex tubes were more likely to be embedded in the regions where the polymers were strongly stretched.In addition,the patterns of microstructures' elongation behave like a filament.From the results of the turbulent kinetic energy budget,it was found that the continuous energy releasing from the polymers into the main flow was the main source of the generation and maintenance of the elastic turbulent status.     

Quenching of fluorescence of conjugated poly(p-phenylene) polymers by benzil and dimethylaminobenzil molecules

We have studied the mechanisms for quenching of the fluorescence of conjugated poly(p-phenylene) polymers by benzil and dimethylaminobenzil molecules. We have shown that molecules in the diketone series are quenching agents for the fluorescence of the indicated polymers, and can serve as singlet-triplet converters capable of populating the triplet state of the polymer. We have observed that the efficiency of quenching of the fluorescence of the studied polymers depends considerably on the presence of bulky side groups in the polymer or in the activator molecules. Based on analysis of the data obtained, we conclude that in the case of a rigid planar structure for the polymer, a significant contribution to quenching of its fluorescence comes from not only singlet-singlet energy transfer but also charge transfer, leading to formation of intermolecular complexes (exciplexes). __________ Translated from Zhurnal Prikladnoi Spektroskopii, Vol. 73, No. 6, pp. 756–762, November–December, 2006.     

Kinetic characteristics of the photoreduction of benzophenone in solid polymers in the presence of 4-halophenols

Benzophenone photoreduction in the presence of 4-halophenols (RC₆H₄OH; R = Cl, Br, and I) in a polymer glass is studied in terms of steady-state and nanosecond laser photolysis. The experimental data on the kinetics of the decay of the ketone triplet state are treated using a polychronous kinetic model in the assumption that two concurrent processes (hydrogen atom abstraction from the polymer and from the phenol) occur. The proton transfer rate constants averaged over the distribution and the parameter n that characterizes the distribution width were determined. The value of k _av for the halophenols is shown to be more than an order of magnitude higher. No heavy atom effect is observed. The reaction product composition is demonstrated to change upon addition of a halophenol. The photoreduction in glass polymers is controlled by hydrogen atom abstraction from the respective donor by triplet ketone molecules. The reaction occurs predominantly in a polymer cage, a kind of polymer nanoreactor.     

Large eddy simulation of channel flow with and without periodic constrictions using the explicit algebraic subgrid-scale model

We analyse the performance of the explicit algebraic subgrid-scale (SGS) stress model (EASSM) in large eddy simulation (LES) of plane channel flow and the flow in a channel with streamwise periodic hill-shaped constrictions (periodic hill flow) which induce separation. The LESs are performed with the Code_Saturne which is an unstructured collocated finite volume solver with a second-order spatial discretisation suitable for LES of incompressible flow in complex geometries. At first, performance of the EASSM in LES of plane channel flow at two different resolutions using the Code_Saturne and a pseudo-spectral method is analysed. It is observed that the EASSM predictions of the mean velocity and Reynolds stresses are more accurate than the conventional dynamic Smagorinsky model (DSM). The results with the pseudo-spectral method were, in general, more accurate. In the second step, LES with the EASSM of flow separation in the periodic hill flow is compared to LES with the DSM, no SGS model and a highly resolved LES data using the DSM. Results show that the mean velocity profiles, the friction and pressure coefficients, the length and shape of the recirculation bubble, as well as the Reynolds stresses are considerably better predicted by the EASSM than the DSM and the no SGS model simulations. It was also observed that in some parts of the domain, the resolved strain-rate and SGS shear stress have the same sign. The DSM cannot produce a correct SGS stress in this case, in contrast to the EASSM.     

A priori assessment of the subgrid scale viscous/scalar dissipation closures in compressible turbulence

An appraisal is made of several subgrid scale (SGS) viscous/scalar dissipation closures via a priori analysis of direct numerical simulation data in a temporally evolving compressible mixing layer. The effects of the filter width, the compressibility level and the Schmidt number are studied for several models. Based on the scaling of SGS kinetic energy, a new formulation for SGS viscous dissipation is proposed. This yields the best overall prediction of the SGS viscous dissipation within the inertial subrange. An SGS scalar dissipation model based on the proportionality of the turbulent time scale with the scalar mixing time scale also performs the best for the filter widths in the inertial subrange. Two dynamic methods are implemented for the determination of the model coefficients. The one based on the global equilibrium of dissipation and production is shown to be more satisfactory than the conventional dynamic model.