FLOW STRUCTURE IN THE TURBULENT BOUNDARY LAYER OVER DIRECTIONAL RIBLETS SURFACES

本文采用时间解析的二维粒子图像测速技术,对零压力梯度光滑以及汇聚和发散沟槽表面平板湍流边界层统计特性和流动结构进行了研究.结果表明在垂直于汇聚和发散沟槽表面的对称平面内,相对于光滑壁面,发散沟槽壁面使当地边界层厚度、壁面摩擦阻力、湍流脉动、雷诺应力等明显减小;而汇聚沟槽壁面对湍流边界层特性和流动结构的影响正好相反,汇聚沟槽使壁面流体有远离壁面向上运动的趋势,因而导致边界层厚度增加了约43%;同时,在汇聚沟槽表面情况下流向大尺度相干结构更容易形成,这对减阻是不利的.此外,顺向涡数量在湍流边界层的对数区均存在一个极大值,发散沟槽表面所对应的极大值位置更靠近沟槽壁面,而在汇聚沟槽表面则有远离壁面的趋势,由顺向涡诱导产生的较强的喷射和扫掠运动会在湍流边界层中产生较强的剪切作用,顺向涡数量的减少是发散沟槽壁面当地摩擦阻力降低的主要原因.     

Layered structure of turbulent plane wall jet

This paper investigates the layered structure of a turbulent plane wall jet at a distance from the nozzle exit. Based on the force balances in the mean momentum equation, the turbulent plane wall jet is divided into three regions: a boundary layer-like region (BLR) adjacent to the wall, a half free jet-like region (HJR) away from the wall, and a plug flow-like region (PFR) in between. In the PFR, the mean streamwise velocity is essentially the maximum velocity, and the simplified mean continuity and mean momentum equations result in a linear variation of the mean wall-normal velocity and Reynolds shear stress. In the HJR, as in a turbulent free jet, a proper scale for the mean wall-normal flow is the mean wall-normal velocity far from the wall and a proper scale for the Reynolds shear stress is the product of the maximum mean streamwise velocity and the velocity scale for the mean wall-normal flow. The BLR region can be divided into four sub-layers, similar to those in a canonical pressure-driven turbulent channel flow or shear-driven turbulent boundary layer flow. Building on the log-law for the mean streamwise velocity in the BLR, a new skin friction law is proposed for a turbulent wall jet. The new prediction agrees well with the correlation of Bradshaw and Gee (1960) over moderate Reynolds numbers, but gives larger skin frictions at higher Reynolds numbers.     

Near wall characterization of the flow over a two-dimensional steep smooth hill

Near-wall data for the strongly perturbed flow in a neutrally stable boundary layer encountering a steep, smooth, two-dimensional hill are presented. Observations were made on the centerplane of a water channel at thirteen stations relative to the hill by laser Doppler anemometry. The large reverse flow region that is formed on the lee of the hill was particularly scrutinized through seven measuring stations. Results are presented for the mean and turbulent properties of the flow. Wall shear stress was evaluated through fitting procedures that resorted to the near wall behavior of the velocity profile. Logarithmic fits as well as predictions through the Reynolds stress profiles are also presented.     

An experimental study of the dynamic characteristics in boundary layer in the flow past a two-dimensional circular cylinder at subcritical reynolds numbers

An experimental study of the dynamic characteristics of flow past a two-dimensional circular cylinder is described. The fluctuationsoof wall shear stress, surface-pressure and velocity of the flow are measured with hot-film, hot-wire and pressure transducer. The frequency feature of fluctuations of wall shear stress is given. The cross-correlation functions of these fluctuations at any two points are calculated. The experimental results reveal that there is an overall syncronous fluctuation, at the shedding frequency, in boundary layer in the flow past a two-dimensional circular cylinder at subcritical Reynolds number.     

Experiments in Turbulent Pipe Flow with Polymer Additives at Maximum Drag Reduction

In this paper we report on (two-component) LDV experiments in a fully developed turbulent pipe flow with a drag-reducing polymer (partially hydrolyzed polyacrylamide) dissolved in water. The Reynolds number based on the mean velocity, the pipe diameter and the local viscosity at the wall is approximately 10000. We have used polymer solutions with three different concentrations which have been chosen such that maximum drag reduction occurs. The amount of drag reduction found is 60–70%. Our experimental results are compared with results obtained with water and with a very dilute solution which exhibits only a small amount of drag reduction. We have focused on the observation of turbulence statistics (mean velocities and turbulence intensities) and on the various contributions to the total shear stress. The latter consists of a turbulent, a solvent (viscous) and a polymeric part. The polymers are found to contribute significantly to the total stress. With respect to the mean velocity profile we find a thickening of the buffer layer and an increase in the slope of the logarithmic profile. With respect to the turbulence statistics we find for the streamwise velocity fluctuations an increase of the root mean square at low polymer concentration but a return to values comparable to those for water at higher concentrations. The root mean square of the normal velocity fluctuations shows a strong decrease. Also the Reynolds (turbulent) shear stress and the correlation coefficient between the stream wise and the normal components are drastically reduced over the entire pipe diameter. In all cases the Reynolds stress stays definitely non-zero at maximum drag reduction. The consequence of the drop of the Reynolds stress is a large polymer stress, which can be 60% of the total stress. The kinetic-energy balance of the mean flow shows a large transfer of energy directly to the polymers instead of the route by turbulence. The kinetic energy of the turbulence suggests a possibly negative polymeric dissipation of turbulent energy. This revised version was published online in July 2006 with corrections to the Cover Date.     

BOUNDARY TREATMENT AND AN EFFICIENT PRESSURE ALGORITHM FOR INTERNAL TURBULENT FLOWS USING THE PDF METHOD

The generalized Langevin model, which is used to model the motion of stochastic particles in the velocity–composition joint probability density function (PDF) method for reacting turbulent flows, has been extended to incorporate solid wall effects. Anisotropy of Reynolds stresses in the near-wall region has been addressed. Numerical experiments have been performed to demonstrate that the forces in the near-wall region of a turbulent flow cause the stochastic particles approachi ng a solid wall to reverse their direction of motion normal to the wall and thereby, leave the near-wall layer. This new boundary treatment has subsequently been implemented in a full-scale problem to prove its validity. The test problem considered here is that of an isothermal, non-reacting turbulent flow in a two-dimensional channel with plug inflow and a fixed back-pressure. An efficient pressure correction method, developed in the spirit of the PISO algorithm, has been implemented. The pressure correction strategy is easy to implement and is completely consistent with the time- marching scheme used for the solution of the Lagrangian momentum equations. The results show remarkable agreement with both k–ϵ and algebraic Reynolds stress model calculations for the primary velocity. The secondary flow velocity and the turbulent moments are in better agreement with the algebraic Reynolds stress model predictions than the k– ϵ predictions. © 1997 by John Wiley & Sons, Ltd.     

Reynolds Stress Budgets in Couette and Boundary Layer Flows

Reynolds stress budgets for both Couette and boundary layer flows are evaluated and presented. Data are taken from direct numerical simulations of rotating and non-rotating plane turbulent Couette flow and turbulent boundary layer with and without adverse pressure gradient. Comparison of the total shear stress for the two types of flows suggests that the Couette case may be regarded as the high Reynolds number limit for the boundary layer flow close to the wall. The limit values of turbulence statistics close to the wall for the boundary layer for increasing Reynolds number approach the corresponding Couette flow values. The direction of rotation is chosen so that it has a stabilizing effect, whereas the adverse pressure gradient is destabilizing. The pressure-strain rate tensor in the Couette flow case is presented for a split into slow, rapid and Stokes terms. Most of the influence from rotation is located to the region close to the wall, and both the slow and rapid parts are affected. The anisotropy for the boundary layer decreases for higher Reynolds number, reflecting the larger separation of scales, and becomes close to that for Couette flow. The adverse pressure gradient has a strong weakening effect on the anisotropy. All of the data presented here are available on the web [36]. This revised version was published online in July 2006 with corrections to the Cover Date.     

Riblet flow calculation with a low Reynolds number κ - ε model

A low Reynolds number κ - ε model has been used to calculate the turbulent boundary layer over riblets. Calculated mean velocity, Reynolds shear stress and kinetic energy distributions are generally in good agreement with available experimental data. The comparison between these distributions and those in a corner flow points to a significant difference between the two flows and the unlikelihood of counter-rotating vortices within the riblet grooves. One shortcoming of the present κ - ε model is the relatively slow return to a two-dimensional turbulence state as the distance from the riblet surface increases.     

Reynolds number effects on a turbulent boundary layer with separation,reattachment, and recovery

The present paper addresses experimental studies of Reynolds number effects on a turbulent boundary layer with separation, reattachment, and recovery. A momentum thickness Reynolds number varies from 1,100 to 20,100 with a wind tunnel enclosed in a pressure vessel by varying the air density and wind tunnel speed. A custom-built, high-resolution laser Doppler anemometer provides fully resolved turbulence measurements over the full Reynolds number range. The experiments show that the mean flow is at most a very weak function of Reynolds number while turbulence quantities strongly depend on Reynolds number. Roller vortices are generated in the separated shear layer caused by the Kelvin–Helmholtz instability. Empirical Reynolds number scalings for the mean velocity and Reynolds stresses are proposed for the upstream boundary layer, the separated region, and the recovery region. The inflectional instability plays a critical role in the scaling in the separated region. The near-wall flow recovers quickly downstream of reattachment even if the outer layer is far from an equilibrium state. As a result, a stress equilibrium layer where a flat-plate boundary layer scaling is valid develops in the recovery region and grows outward moving downstream.     

Numerical and experimental investigation of turbulent characteristics in a drag-reducing flow with surfactant additives

Cetyltrimethyl ammonium chloride (CTAC) surfactant additives, because of their long-life characteristics, can be used as promising drag-reducers in district heating and cooling systems. In the present study we performed both numerical and experimental tests for a 75 ppm CTAC surfactant drag-reducing channel flow. A two-component PIV system was used to measure the instantaneous streamwise and wall-normal velocity components. A Giesekus constitutive equation was adopted to model the extra stress due to the surfactant additives, with the constitutive parameters being determined by well-fitting apparent shear viscosities, as measured by an Advanced Rheometric Expansion System (ARES) rheometer. In the numerical study, we connected the realistic rheological properties with the drag-reduction rate. This is different from previous numerical studies in which the model parameters were set artificially. By performing consistent comparisons between numerical and experimental results, we have obtained an insight into the mechanism of the additive-induced drag-reduction phenomena.

Our simulation showed that the addition of surfactant additives introduces several changes in turbulent flow characteristics: (1) In the viscous sublayer, the mean velocity gradient becomes gentler due to the viscoelastic forces introduced by the additives. The buffer layer becomes expanded and the slope of the velocity profile in the logarithmic layer increases. (2) The locations where the streamwise velocity fluctuation and Reynolds shear stress attain their maximum value shifted from the wall region to the bulk flow region. (3) The root-mean-square velocity fluctuations in the wall-normal direction decrease for the drag-reducing flow. (4) The Reynolds shear stress decreases dramatically and the deficit of the Reynolds shear stress is mainly compensated by the viscoelastic shear stress. (5) The turbulent production becomes much smaller and its peak-value position moves toward the bulk flow region. All of these findings agree qualitatively with experimental measurements.

Regarding flow visualization, the violent streamwise vortices in the near wall region become dramatically suppressed, indicating that the additives weaken the ejection and sweeping motion, and thereby inhibit the generation of turbulence. The reduction in turbulence is accomplished by additive-introduced viscoelastic stress. Surfactant additives have dual effects on frictional drag: (1) introduce viscoelastic shear stress, which increases frictional drag; and (2) dampen the turbulent vortical structures, decrease the turbulent shear stress, and then decrease the frictional drag. Since the second effect is greater than the first one, drag-reduction occurs.     

Anisotropy measurements in the boundary layer over a flat plate with suction

Laser Doppler velocity measurements are carried out in a turbulent boundary layer subjected to concentrated wall suction (through a porous strip). The measurements are taken over a longitudinal distance of 9× the incoming boundary layer thickness ahead of the suction strip. The mean and rms velocity profiles are affected substantially by suction. Two-point measurements show that the streamwise and wall-normal autocorrelations of the streamwise velocity are reduced by suction. It is found that suction alters the redistribution of the turbulent kinetic energy k between its components. Relative to the no-suction case, the longitudinal Reynolds stress contributes more to k than the other two normal Reynolds stresses; in the outer region, its contribution is reduced which suggests structural changes in the boundary layer. This is observed in the anisotropy of the Reynolds stresses, which depart from the non-disturbed boundary layer. With suction, the anisotropy level in the near-wall region appears to be stronger than that of the undisturbed layer. It is argued that the mean shear induced by suction on the flow is responsible for the alteration of the anisotropy. The variation of the anisotropy of the layer will make the development of a turbulence model quite difficult for the flow behind suction. In that respect, a turbulence model will need to reproduce well the effects of suction on the boundary layer, if the model is to capture the effect of suction on the anisotropy of the Reynolds stresses.     

Characteristic Regimes of Transitional Separation Bubbles in Unsteady Flow

This experimental investigation deals with transition phenomena of a separated boundary layer under unsteady inlet flow conditions. The main purpose of this investigation is to understand the influence of the rotor-stator interaction in turbomachinery on the subsequent, highly loaded boundary layer. The research project is divided into two phases. In the first phase, which has been completed recently, only the variation of mean velocity caused by upstream blades was simulated in the experiments while the free-stream turbulence intensity was retained at a constant low level. The experiments are carried out in an Eifel-type wind tunnel to investigate the laminar separated boundary layer of a flat plate under oscillating inlet conditions. The adverse pressure gradient, similar to that of turbomachines, is generated by the contoured upper wall. The unsteadiness is produced by a rotating flap located downstream of the test section. The reduced frequency, the amplitude and the mean Reynolds number are varied to simulate the conditions prevailing in turbomachines. In addition to the Kelvin–Helmholtz instability of the separated shear layer, a lower frequency instability was observed. This is frequently referred to as `free shear layer flapping' and results in two distinctly different ways of re-attachment, depending primarily on the Reynolds number. For low momentum thickness Reynolds numbers at the separation point, large-scale vortices locked to the frequency of the unsteady main flow are identified. They originate nearly at the top of the separation bubble and are ejected downstream. A fully turbulent boundary layer develops after these vortices mix out. For higher Reynolds numbers, transition is completed within a short length of the free shear layer and there-attachment region. The characteristic momentum thickness Reynolds number separating these two regimes in unsteady flow is about 125. The Strouhal number (reduced frequency) does not appear to have any significant effect. Based on the experimental results, this behaviour is discussed in some detail. This revised version was published online in July 2006 with corrections to the Cover Date.     

Three-dimensional particle-tracking velocimetry measurement of turbulence statistics and energy budget in a backward-facing step flow

Using a three-dimensional (3-D) particle-tracking velocimeter, detailed turbulent flow measurements were made in a plane channel with a one-sided 50% abrupt expansion, which acted as a backward-facing step. The turbulent channel flow reached a fully developed state well upstream of the step. The Reynolds number based on the upstream centerline velocity and the step height H was 5540. With the mean reattachment point located at 6.51H downstream of the step, the measurement region ranged from −2H upstream to 12H downstream of the step. Various turbulent statistics and the energy budget were calculated from numerous instantaneous vector distributions. As in previous experimental investigations, the Reynolds normal and shear stresses had maximum values upstream of the reattachment. The stress anisotropy tensor revealed a peculiar phenomenon near the reattachment wall, wherein the spanwise normal stress was the largest among the three normal stresses. The triple velocity correlations indicated large values in the separating shear layer, and hence the turbulent diffusion was a major term in the energy budget. Comparison was made between the present results and those of the direct numerical simulation (DNS) of Le et al. (1993), and it was found that the mean and fluctuating velocities, the Reynolds shear stress, and the turbulent energy budget were in excellent agreement, although there was a considerable difference in the inflow conditions.     

Instantaneous structure and statistical feature of unsteady flow in a channel obstructed by a square rod

A two-dimensional numerical computation has been made for an unsteady flow in a channel obstructed by an inserted square rod. The results of the computation made for the flow with a parabolic inlet velocity profile at a specific value of channel Reynolds number are analyzed in detail. The obtained results reveal that momentum transfer is enhanced due to the apparent shear stress resulting from the nonzero value of cross-correlation between the streamwise and normal components of fluctuating velocity, , just as in turbulent shear flows, although the studied flow is quite different from turbulent flows in the sense that it is highly periodical and therefore free from randomness. This periodicity leads to a quick recovery of the velocity defect in some region of the wake of the rod. Special attention is paid to the time variation of flow structure. The crisscross motion of the Karman vortex previously found to occur is discussed again, and how it appears is explained in terms of the interaction between the Karman vortex and the disturbed wall shear layer. In the discussion of this relationship, wavering motion of the separation vorticity layers formed on both sides of the rod and the periodic formation of an isolated vortex island from the lifted tip of the wall vorticity layer are analyzed. The vortex island is found to play an important role not only for the occurrence of the crisscross motion of Karman vortex but also for the generation of the nonzero value of .     

Flow field investigation in rotating rib-roughened channel by means of particle image velocimetry

The turbulent velocity field over the rib-roughened wall of an orthogonally rotating channel is investigated by means of two-dimensional particle image velocimetry (PIV). The flow direction is outward, with a bulk Reynolds number of 1.5 × 10⁴ and a rotation number ranging from 0.3 to 0.38. The measurements are obtained along the wall-normal/streamwise plane at mid-span. The PIV system rotates with the channel, allowing to measure directly the relative flow velocity with high spatial resolution. Coriolis forces affect the stability of the boundary layer and free shear layer. Due to the different levels of shear layer entrainment, the reattachment point is moved downstream (upstream) under stabilizing (destabilizing) rotation, with respect to the stationary case. Further increase in rotation number pushes further the reattachment point in stabilizing rotation, but does not change the recirculation length in destabilizing rotation. Turbulent activity is inhibited along the leading wall, both in the boundary layer and in the separated shear layer; the opposite is true along the trailing wall. Coriolis forces affect indirectly the production of turbulent kinetic energy via the Reynolds shear stresses and the mean shear. Two-point correlation is used to characterize the coherent motion of the separated shear layer. Destabilizing rotation is found to promote large-scale coherent motions and accordingly leads to larger integral length scales; on the other hand, the spanwise vortices created in the separating shear layer downstream of the rib are less organized and tend to be disrupted by the three-dimensional turbulence promoted by the rotation. The latter observation is consistent with the distributions of span-wise vortices detected in instantaneous flow realizations.     

Turbulent Heat and Fluid Flow over a Two-Dimensional Hill

Experimental investigation has been made on the flow and thermal fields over a heated two-dimensional hill with a cosine-squared shape. The detailed turbulent characteristics are measured by a backscatter-type two-component LDV, a PIV system, a fine thermocouple and a cold-wire probe. In the reverse-flow region on the leeward side of the hill, the turbulence intensities and the Reynolds shear stress show much larger values than in a canonical wall-bounded shear flow. The mean temperature maintains a relatively high value below the location where the horizontal mean velocity rapidly decreases. At the outer edge of the reverse-flow region, there exists a second maximum intensity of temperature fluctuations. The instantaneous temperature waveforms near the heated surface show very large amplitude consisting of high-frequency fluctuations superimposed on the low-frequency motions. Simultaneous measurement of velocity and temperature is also done using a combination of a two-component LDV and a fine-wire thermocouple together with digital response compensation. In the downstream region of the hill, the horizontal and vertical turbulent heat fluxes become maximum at the hill-top height, and tend to diffuse in the vertical direction. On the other hand, in the reverse-flow region formed behind the hill, both of the heat fluxes decrease remarkably. In particular, it is worth noting that the horizontal turbulent heat flux near the surface becomes opposite in sign to that in the forward flow region. This is mainly due to the reversal of the mean flow direction.     

Numerical analysis of turbulent flow separation in a rectangular duct with a sharp 180‐degree turn by algebraic Reynolds stress model

Turbulent flow in a rectangular duct with a sharp 180‐degree turn is difficult to predict numerically because the flow behavior is influenced by several types of forces, including centrifugal force, pressure‐driven force, and shear stress generated by anisotropic turbulence. In particular, this type of flow is characterized by a large‐scale separated flow, and it is difficult to predict the reattachment point of a separated flow. Numerical analysis has been performed for a turbulent flow in a rectangular duct with a sharp 180‐degree turn using the algebraic Reynolds stress model. A boundary‐fitted coordinate system is introduced as a method for coordinate transformation to set the boundary conditions next to complicated shapes. The calculated results are compared with the experimental data, as measured by a laser‐Doppler anemometer, in order to examine the validity of the proposed numerical method and turbulent model. In addition, the possibility of improving the wall function method in the separated flow region is examined by replacing the log‐law velocity profile for a smooth wall with that for a rough wall. The analysis results indicated that the proposed algebraic Reynolds stress model can be used to reasonably predict the turbulent flow in a rectangular duct with a sharp 180‐degree turn. In particular, the calculated reattachment point of a separated flow, which is difficult to predict in a turbulent flow, agrees well with the experimental results. In addition, the calculation results suggest that the wall function method using the log‐law velocity profile for a rough wall over a separated flow region has some potential for improving the prediction accuracy. Copyright © 2007 John Wiley & Sons, Ltd.     

On the effects of dilute polymers on driven cavity turbulent flows

Effects of dilute polymer solutions on a lid-driven cubical cavity turbulent flow are studied via particle image velocimetry (PIV). This canonical flow is a combination of a bounded shear flow, driven at constant velocity and vortices that change their spatial distribution as a function of the lid velocity. From the two-dimensional PIV data we estimate the time averaged spatial fields of key turbulent quantities. We evaluate a component of the vorticity–velocity correlation, namely 〈ω₃v〉, which shows much weaker correlation, along with the reduced correlation of the fluctuating velocity components, u and v. There are two contributions to the reduced turbulent kinetic energy production −〈u v〉S_uv, namely the reduced Reynolds stresses, −〈u v〉, and strongly modified pointwise correlation of the Reynolds stress and the mean rate-of-strain field, S_uv. The Reynolds stresses are shown to be affected because of the derivatives of the Reynolds stresses, ∂〈u v〉/∂y that are strongly reduced in the same regions as the vorticity–velocity correlation. The results, combined with the existing evidence, support the phenomenological model of polymer effects propagating from the polymer scale to the velocity derivatives and through the mixed-type correlations and Reynolds stress derivatives up to the turbulent velocity fields. The effects are shown to be qualitatively similar in different flows regardless of forcing type, homogeneity or presence of liquid–solid boundaries.     

Turbulence measurements in the inlet plane of a centrifugal compressor vaneless diffuser

Detailed flow measurements at the inlet of a centrifugal compressor vaneless diffuser are presented. The mean 3-d velocities and six Reynolds stress components tensor are used to determine the turbulence production terms which lead to total pressure loss. High levels of turbulence kinetic energy were observed in both the blade and passage wakes, but these were only associated with high Reynolds stresses in the blade wakes. For this reason the blade wakes mixed out rapidly, whereas the passage wake maintained its size, but was redistributed across the full length of the shroud wall. Peak levels of Reynolds stress occurred in regions of high velocity shear and streamline curvature which would tend to destabilize the shear gradient. Four regions in the flow are identified as potential sources of loss - the blade wake, the shear layers between passage wake and jet, the thickened hub boundary layer and the interaction region between the secondary flow within the blade wake and the passage vortex. The blade wakes generate most turbulence, with smaller contributions from the hub boundary layer and secondary flows, but no significant contribution is apparent from the passage wake shear layers.     

倾斜吹吸控制下湍流边界层减阻的直接数值模拟

雷诺切应力是壁湍流高摩擦阻力的重要来源, 有理论认为可以通过壁面生成负雷诺应力(数值上为正)的方式来削弱湍流流场中雷诺应力的分布, 以此获得流动减阻. 而通过对雷诺平均运动方程的法向二次积分, 可以发现壁面生成正雷诺应力(数值上为负)对壁面摩擦阻力系数才有负贡献. 文中在湍流边界层流动的控制区域下边界设置一系列倾斜狭缝, 利用该装置通过周期性吹吸的方法产生壁面生成正(负)雷诺应力, 并采用直接数值模拟方法考察和验证上文提到的减阻理论. 文中采用的湍流边界层流动模型, 其流动雷诺数(基于外流速度及动量损失厚度)从300 发展到860. 文中通过多组数值模拟算例, 考察了射流强度和频率对壁面摩擦阻力系数的影响, 并对比了壁面生成正或负雷诺应力对流动的影响. 研究表明, 壁面生成正雷诺应力控制的减阻率能达到3.26, 而壁面生成负雷诺应力控制的减阻效果较壁面生成正雷诺应力控制的要差; 壁面生成的正雷诺应力对壁面摩擦阻力有负贡献, 而壁面生成的负雷诺应力对壁面摩擦阻力有正贡献; 通过考察控制的收支比, 发现控制方案不能获得能量净收益.