We report the results of a complete modal and nonmodal linear stability analysis of the electrohydrodynamic flow (EHD) for the problem of electroconvection in the strong injection region. Convective cells are formed by Coulomb force in an insulating
liquid residing between two plane electrodes subject to unipolar injection. Besides pure electroconvection, we also consider the case where a cross-flow is present, generated by a streamwise pressure gradient, in the form of a laminar Poiseuille flow. The effect of charge diffusion, often neglected in previous linear stability analyses, is included in the present study and a transient growth analysis, rarely considered in EHD, is carried out. In the case without cross-flow, a non-zero charge diffusion leads to a lower linear stability threshold and thus to a more unstable low. The transient growth, though enhanced by increasing charge diffusion, remains small and hence cannot fully account for the discrepancy of the linear stability threshold between theoretical and experimental results. When a cross-flow is present, increasing the strength of the electric field in the high-$Re$ Poiseuille flow yields a more unstable flow in both modal and nonmodal stability analyses. Even though the energy analysis and the input-output analysis both indicate that the energy growth directly related to the electric field is small, the electric effect enhances the lift-up mechanism. The symmetry of channel flow with respect to the centerline is broken due to the additional electric field acting in the wall-normal direction. As a result, the centers of the streamwise rolls are shifted towards the injector electrode, and the optimal spanwise wavenumber achieving maximum transient energy growth increases with the strength of the electric field.
The numerical simulation of a flow through a duct requires an externally specified forcing that makes the fluid flow against viscous friction. To this aim, it is customary to enforce a constant value for either the flow rate (CFR) or the pressure gra
dient (CPG). When comparing a laminar duct flow before and after a geometrical modification that induces a change of the viscous drag, both approaches (CFR and CPG) lead to a change of the power input across the comparison. Similarly, when carrying out the (DNS and LES) numerical simulation of unsteady turbulent flows, the power input is not constant over time. Carrying out a simulation at constant power input (CPI) is thus a further physically sound option, that becomes particularly appealing in the context of flow control, where a comparison between control-on and control-off conditions has to be made. We describe how to carry out a CPI simulation, and start with defining a new power-related Reynolds number, whose velocity scale is the bulk flow that can be attained with a given pumping power in the laminar regime. Under the CPI condition, we derive a relation that is equivalent to the Fukagata--Iwamoto--Kasagi relation valid for CFR (and to its extension valid for CPG), that presents the additional advantage of natively including the required control power. The implementation of the CPI approach is then exemplified in the standard case of a plane turbulent channel flow, and then further applied to a flow control case, where the spanwise-oscillating wall is used for skin friction drag reduction. For this low-Reynolds number flow, using 90% of the available power for the pumping system and the remaining 10% for the control system is found to be the optimum share that yields the largest increase of the flow rate above the reference case, where 100% of the power goes to the pump.
In order to generalize the well-known spanwise-oscillating-wall technique for drag reduction, non-sinusoidal oscillations of a solid wall are considered as a means to alter the skin-friction drag in a turbulent channel flow. A series of Direct Numeri
cal Simulations is conducted to evaluate the control performance of nine different temporal waveforms, in addition to the usual sinusoid, systematically changing the wave amplitude and the period for each waveform. The turbulent average spanwise motion is found to coincide with the laminar Stokes solution that is constructed, for the generic waveform, through harmonic superposition. This allows us to define and compute, for each waveform, a new penetration depth of the Stokes layer which correlates with the amount of turbulent drag reduction, and eventually to predict both turbulent drag reduction and net energy saving rate for arbitrary waveforms. Among the waveforms considered, the maximum net energy saving rate is obtained by the sinusoidal wave at its optimal amplitude and period. However, the sinusoid is not the best waveform at every point in the parameter space. Our predictive tool offers simple guidelines to design waveforms that outperform the sinusoid for given (suboptimal) amplitude and period of oscillation. This is potentially interesting in view of applications, where physical limitations often preclude the actuator to reach its optimal operating conditions.
Flow control with the goal of reducing the skin friction drag on the fluid-solid interface is an active fundamental research area, motivated by its potential for significant energy savings and reduced emissions in the transport sector. Customarily, t
he performance of drag reduction techniques in internal flows is evaluated under two alternative flow conditions, i.e. at constant mass flow rate or constant pressure gradient. Successful control leads to reduction of drag and pumping power within the former approach, whereas the latter leads to an increase of the mass flow rate and pumping power. In practical applications, however, money and time define the flow control challenge: a compromise between the energy expenditure (money) and the corresponding convenience (flow rate) achieved with that amount of energy has to be reached so as to accomplish a goal which in general depends on the specific application. Based on this idea, we derive two dimensionless parameters which quantify the total energy consumption and the required time (convenience) for transporting a given volume of fluid through a given duct. Performances of existing drag reduction strategies as well as the influence of wall roughness are re-evaluated within the present framework; how to achieve the (application-dependent) optimum balance between energy consumption and convenience is addressed. It is also shown that these considerations can be extended to external flows.
Harmonic oscillations of the walls of a turbulent plane channel flow are studied by direct numerical simulations to improve our understanding of the physical mechanism for skin-friction drag reduction. The simulations are carried out at constant pres
sure gradient in order to define an unambiguous inner scaling: in this case, drag reduction manifests itself as an increase of mass flow rate. Energy and enstrophy balances, carried out to emphasize the role of the oscillating spanwise shear layer, show that the viscous dissipations of the mean flow and of the turbulent fluctuations increase with the mass flow rate, and the relative importance of the latter decreases. We then focus on the turbulent enstrophy: through an analysis of the temporal evolution from the beginning of the wall motion, the dominant, oscillation-related term in the turbulent enstrophy is shown to cause the turbulent dissipation to be enhanced in absolute terms, before the slow drift towards the new quasi-equilibrium condition. This mechanism is found to be responsible for the increase in mass flow rate. We finally show that the time-average volume integral of the dominant term relates linearly to the drag reduction.
