Birds have to flap their wings to generate the needed thrust force, which powers them through the air. But how exactly do flapping wings create such force, and at what amplitude and frequency should they operate? These questions have been asked by many researchers. It turns out that much of the secret is hidden in the wake left behind the flapping wing. Exemplified by the study of Andersen et al. (J. Fluid Mech., vol. 812, 2017, R4), close examination of the flow pattern behind a flapping wing will inform us whether the wing is towed by an external force or able to generate a net thrust force by itself. Such studies are much like looking at the footprints of terrestrial animals as we infer their size and weight, figuring out their walking and running gaits. A map that displays the collection of flow patterns after a flapping wing, using flapping frequency and amplitude as the coordinates, offers a full picture of its flying ‘gaits’.

It is a new concept for porous media flow that a hydrodynamic lifting force is generated inside a highly compressible porous layer as a planing surface glides over it. The concept originated from the observation of the pop-out phenomena of red blood cells over the endothelial glycocalyx layer (EGL) lining the inner surface of our blood vessels (Feng & Weinbaum, J. Fluid Mech., vol. 422, 2000, pp. 282–317). In the current paper, we report an experimental study to examine this concept. A novel testing set-up was developed that consists of a running conveyer belt covered with a soft porous sheet, and a fully instrumented upper planar board, i.e. planing surface. The generation of pore pressure was observed and captured by pressure transducers when the planing surface glides over the porous sheet. Its distribution strongly depends on the relative velocity between the planing surface and the running belt, the mechanical and transport properties of the porous sheet as well as the compression ratios at the leading and trailing edges. The relative contribution of the transiently trapped air to the total lift was evaluated by comparing the pore pressure to the total lifting pressure measured by a load cell mounted between two adjacent pressure transducers. For a typical running condition with a polyester porous material ( $k=h_{2}/h_{1}=5$ , $\unicode[STIX]{x1D706}=h_{2}/h_{0}=1$ , $U=3.8~\text{m}~\text{s}^{-1}$ , where $h_{2}$ , $h_{1}$ , are the porous layer thickness at the leading and trailing edges, respectively; $h_{0}$ is the un-deformed porous layer thickness; and $U$ is the velocity of the running belt), over 68 % of the local lift is generated by the pore pressure. The results conclusively verified the validity of lift generation in a highly compressible porous layer as a planing surface glides over it. This study provides the foundation for the application of highly compressible porous media for soft lubrication with minimal frictional losses. It also sheds some light on the biophysics study of the EGL.

A semi-empirical model is presented that describes the development of a fully developed turbulent boundary layer in the presence of surface roughness with length scale $k_{s}$ that varies with streamwise distance $x$ . Interest is centred on flows for which all terms of the von Kármán integral relation, including the ratio of outer velocity to friction velocity $U_{\infty }^{+}\equiv U_{\infty }/u_{\unicode[STIX]{x1D70F}}$ , are streamwise constant. For $Re_{x}$ assumed large, use is made of a simple log-wake model of the local turbulent mean-velocity profile that contains a standard mean-velocity correction for the asymptotic fully rough regime and with assumed constant parameter values. It is then shown that, for a general power-law external velocity variation $U_{\infty }\sim x^{m}$ , all measures of the boundary-layer thickness must be proportional to $x$ and that the surface sand-grain roughness scale variation must be the linear form $k_{s}(x)=\unicode[STIX]{x1D6FC}x$ , where $x$ is the distance from the boundary layer of zero thickness and $\unicode[STIX]{x1D6FC}$ is a dimensionless constant. This is shown to give a two-parameter $(m,\unicode[STIX]{x1D6FC})$ family of solutions, for which $U_{\infty }^{+}$ (or equivalently $C_{f}$ ) and boundary-layer thicknesses can be simply calculated. These correspond to perfectly self-similar boundary-layer growth in the streamwise direction with similarity variable $z/(\unicode[STIX]{x1D6FC}x)$ , where $z$ is the wall-normal coordinate. Results from this model over a range of $\unicode[STIX]{x1D6FC}$ are discussed for several cases, including the zero-pressure-gradient ( $m=0$ ) and sink-flow ( $m=-1$ ) boundary layers. Trends observed in the model are supported by wall-modelled large-eddy simulation of the zero-pressure-gradient case for $Re_{x}$ in the range $10^{8}{-}10^{10}$ and for four values of $\unicode[STIX]{x1D6FC}$ . Linear streamwise growth of the displacement, momentum and nominal boundary-layer thicknesses is confirmed, while, for each $\unicode[STIX]{x1D6FC}$ , the mean-velocity profiles and streamwise turbulent variances are found to collapse reasonably well onto $z/(\unicode[STIX]{x1D6FC}x)$ . For given $\unicode[STIX]{x1D6FC}$ , calculations of $U_{\infty }^{+}$ obtained from large-eddy simulations are streamwise constant and independent of $Re_{x}$ when this is large. The present results suggest that, in the sense that $U_{\infty }^{+}(\unicode[STIX]{x1D6FC},m)$ is constant, these flows can be interpreted as the fully rough limit for boundary layers in the presence of small-scale linear roughness.

We experimentally investigate the behaviour of a line-source plume falling through a finite two-layer stratified ambient where the depth of the upper ambient layer increases in time. Laboratory observations suggest one of two possible flow regimes depending on the value of $\unicode[STIX]{x1D706}$ , which represents the relative loss of buoyancy experienced by the plume upon crossing the ambient interface. When $\unicode[STIX]{x1D706}>1$ , a classical filling-box-type flow is realized and plume fluid always reaches the bottom boundary. By contrast, when $\unicode[STIX]{x1D706}<1$ , we observe a transition by which an increasing fraction of plume fluid discharges along the interface. The approximate start time, $t_{v}$ , and end time, $t_{t}$ , of the transition process are well determined by $\unicode[STIX]{x1D706}$ . After transition, the ambient density evolves to form a three-layer fluid with an intermediate layer that grows in time. Measured densities of the intermediate layer in experiments with $\unicode[STIX]{x1D706}<1$ are well predicted using plume theory. We further characterize the horizontal speed of the intrusion that forms along the ambient interface, the mass of solute present in the intermediate layer at time $t_{t}$ and the rate of descent of the intrusion level for $t>t_{t}$ . The significance of our findings is discussed in the context of the ventilation of natural and hybrid ventilated buildings and of effluent discharge through marine outfall diffusers.

A gas lubricated bearing model is derived which is appropriate for a very small bearing face separation by including velocity slip boundary conditions and centrifugal inertia effects. The bearing dynamics is examined when an external harmonic force is imposed on the bearing due to the bearing being situated within a larger complex dynamical system. A compressible Reynolds equation is formulated for the gas film which is coupled to the bearing structure through an axial force balance where the rotor and stator correspond to spring–mass–damper systems. Surface slip boundary conditions are derived on the bearing faces, characterised by the slip length parameter. The coupled bearing system is analysed using a stroboscopic map solver with the modified Reynolds equation and structural equations solved simultaneously. For a sufficiently large forcing amplitude a flapping motion of the bearing faces is induced when the rotor and stator are in close proximity. The minimum bearing gap over the time period of the external forcing is examined for a range of bearing parameters.

The characteristic time of low-Reynolds-number fluid–structure interaction scales linearly with the ratio of fluid viscosity to solid Young’s modulus. For sufficiently large values of Young’s modulus, both time and length scales of the viscous-elastic dynamics may be similar to acoustic time and length scales. However, the requirement of dominant viscous effects limits the validity of such regimes to micro-configurations. We here study the dynamics of an acoustic plane wave impinging on the surface of a layered sphere, immersed within an inviscid fluid, and composed of an inner elastic sphere, a creeping fluid layer and an external elastic shell. We focus on configurations with similar viscous-elastic and acoustic time and length scales, where the viscous-elastic speed of interaction between the creeping layer and the elastic regions is similar to the speed of sound. By expanding the linearized spherical Reynolds equation into the relevant spectral series solution for the hyperbolic elastic regions, a global stiffness matrix of the layered elastic sphere was obtained. The maximal pressure difference induced by the acoustic wave on the creeping region was found to occur for identical viscous-elastic and acoustic length scales. Comparing an elastic sphere with an embedded creeping layer to a fully elastic sphere, a significant reduction in magnitude and fluctuations (with regard to wavelength) are observed for both the displacements of the solid and target strength of the sphere. This effect was most significant for identical viscous-elastic and acoustic time scales. This work relates viscous-elastic dynamics to acoustic scattering and may pave the way to the design of novel metamaterials with unique acoustic properties.

For Mach reflection in steady supersonic flow, the slipline and reflected shock wave from the triple point are disturbed by secondary Mach waves generated over the slipline and by the expansion fan from the rear wedge corner. Analytical expressions for the shape of the curved slipline and reflected shock wave are derived in this paper. It is found that, due to transmitted expansion waves from the expansion fan, the slipline has a slope discontinuity at the turning point, i.e., the intersection point of the slipline and the leading characteristics of the transmitted expansion wave. The hypothetical shock wave calculated by considering this slope discontinuity as flow deflection angle matches a similar wave observed in numerical results by computational fluid dynamics, suggesting the existence of a weak shock wave from this turning point. The effects of the secondary Mach waves upstream of the turning point and of the turning point weak shock wave mutually cancel out approximately so that the transmitted Mach waves can be approximated as straight characteristic lines. This simplification leads to a fast analytical model which can predict the Mach stem height and shape of the slipline and reflected shock wave with increasing accuracy for the decreasing deflection angle of the slipline at the triple point. The slipline slope discontinuity at the turning point and the hypothetical turning point weak shock wave are new phenomena found in this work.

Gas injection into a liquid cross-flow is examined for the case where the gas is injected beneath a horizontal flat surface. For moderate Froude numbers, the gas pocket that is formed will rise toward the flow boundary under the action of buoyancy, a condition that is conducive to the formation of gas layers for friction-drag reduction on the surface. At the location of gas injection, a plume whose geometry is related to the mass and momentum flux of the injected gas and liquid cross-flow is formed, and the influence of buoyancy is minimal. However, as the gas pocket convects downstream, buoyancy brings the gas back upward to the flow boundary, and leads to the bifurcation of the pocket into two distinct branches, forming a stable ‘V’-shape. Under some conditions, the flow between the two gas branches is almost entirely liquid, while for others there exists a bubbly flow or a continuous sheet of gas between the branches. The sweep angle and cross-sectional geometry of the gas branches are related to free-stream speed and boundary-layer thickness of the liquid cross-flow, the mass-injection rate of the gas, the diameter of the injection orifice and the gas outlet mean velocity and gas–jet angle. Data for a range of experimental conditions are used to scale the flow and results are compared to numerical computations of the flow, and these data are used to illustrate the underlying flow processes responsible leading to the formation the stable and straight gas branches. A simple model based on the balance of forces around a stable gas branch is presented and used to scale the observed data, and we use the results of this analysis and the computations to discuss how the process of gas injection may interact with the formation of the stable gas pockets farther downstream.

Time-resolved velocity measurements are obtained using a hot-wire in a nearly homogeneous and isotropic flow downstream of an active grid for a range of Taylor Reynolds numbers from $191$ to $659$ . It is found that the dimensionless dissipation rate, $C_{\unicode[STIX]{x1D716}}$ , is nearly a constant for sufficiently high values of Taylor Reynolds number, $R_{\unicode[STIX]{x1D706},u_{q}}$ , and is approximately equal to $0.87$ . This value is approximately $5\,\%$ less than the value reported by Bos et al. (Phys. Fluids, vol. 19 (4), 2007, 045101), which is obtained using DNS/LES (direct numerical simulation combined with large eddy simulation) for decaying homogeneous and isotropic turbulence, and is in excellent agreement with the active grid experiment of Thormann & Meneveau (Phys. Fluids, vol. 26 (2), 2014, 025112.). The results presented herein show that deviation from isotropy may cause inconsistencies in the computation of $C_{\unicode[STIX]{x1D716}}$ . As a result, it is suggested that the velocity scale be the square root of the turbulence kinetic energy. The integral length scale measurements obtained from the longitudinal velocity correlation are in close agreement with the integral length scale measured from the peak of the energy spectrum, $\unicode[STIX]{x1D705}E_{11}(\unicode[STIX]{x1D705})$ , where $\unicode[STIX]{x1D705}$ is the wavenumber and $E_{11}(\unicode[STIX]{x1D705})$ is the one-dimensional power spectrum of the downstream velocity.

A comprehensive study of the classical linear spin-down of a constant-density viscous fluid (kinematic viscosity $\unicode[STIX]{x1D708}$ ) rotating rapidly (angular velocity $\unicode[STIX]{x1D6FA}$ ) inside an axisymmetric cylindrical container (radius $L$ , height $H$ ) with rigid boundaries, which follows the instantaneous small change in the boundary angular velocity at small Ekman number $E=\unicode[STIX]{x1D708}/H^{2}\unicode[STIX]{x1D6FA}\ll 1$ , was provided by Greenspan & Howard (J. Fluid Mech., vol. 17, 1963, pp. 385–404). For that problem $E^{1/2}$ Ekman layers form quickly, triggering inertial waves together with the dominant spin-down of the quasi-geostrophic (QG) interior flow on the $O(E^{-1/2}\unicode[STIX]{x1D6FA}^{-1})$ time scale. On the longer lateral viscous diffusion time scale $O(L^{2}/\unicode[STIX]{x1D708})$ , the QG flow responds to the $E^{1/3}$ sidewall shear layers. In our variant, the sidewall and top boundaries are stress-free, a set-up motivated by the study of isolated atmospheric structures such as tropical cyclones or tornadoes. Relative to the unbounded plane layer case, spin-down is reduced (enhanced) by the presence of a slippery (rigid) sidewall. This is evidenced by the QG angular velocity, $\unicode[STIX]{x1D714}^{\star }$ , evolution on the $O(L^{2}/\unicode[STIX]{x1D708})$ time scale: spatially, $\unicode[STIX]{x1D714}^{\star }$ increases (decreases) outwards from the axis for a slippery (rigid) sidewall; temporally, the long-time $(\gg L^{2}/\unicode[STIX]{x1D708})$ behaviour is dominated by an eigensolution with a decay rate slightly slower (faster) than that for an unbounded layer. In our slippery sidewall case, the $E^{1/2}\times E^{1/2}$ corner region that forms at the sidewall intersection with the rigid base is responsible for a $\ln E$ singularity within the $E^{1/3}$ layer, causing our asymptotics to apply only at values of $E$ far smaller than can be reached by our direct numerical simulation (DNS) of the linear equations governing the entire spin-down process. Instead, we solve the $E^{1/3}$ boundary layer equations for given $E$ numerically. Our hybrid asymptotic–numerical approach yields results in excellent agreement with our DNS.

Steady incompressible flow down a slowly curving circular pipe is considered, analytically and numerically. Both real and complex solutions are investigated. Using high-order Hermite–Padé approximants, the Dean series solution is analytically continued outside its circle of convergence, where it predicts a complex solution branch for real positive Dean number, $K$ . This is confirmed by numerical solution. It is shown that other previously unknown solution branches exist for all $K>0$ , which are related to an unforced complex eigensolution. This non-uniqueness is believed to be generic to the Navier–Stokes equations in most geometries. By means of path continuation, numerical solutions are followed around the complex $K$ -plane. The standard Dean two-vortex solution is shown to lie on the same hypersurface as the eigensolution and the four-vortex solutions found in the literature. Elliptic pipes are considered and shown to exhibit similar behaviour to the circular case. There is an imaginary singularity limiting convergence of the Dean series, an unforced solution at $K=0$ and non-uniqueness for $K>0$ , culminating in a real bifurcation.

Three-dimensional (3D) wake transition for a circular cylinder placed near to a moving wall is investigated using direct numerical simulation (DNS). The study covers a parameter space spanning a gap ratio $(G/D)\geqslant 0.3$ and Reynolds number ( $Re$ ) up to 325. The wake transition regimes in the parameter space are mapped out. It is found that vortex dislocation associated with Mode A is completely suppressed at $G/D$ smaller than approximately 1.0. The suppression of vortex dislocation is believed to be due to the confinement of the Mode A streamwise vortices by the plane wall, which suppresses the excess growth and local dislocation of any Mode A vortex loop. Detailed wake transition is examined at $G/D=0.4$ , where the wake transition sequence is ‘two-dimensional (2D) $\rightarrow$ ordered Mode A $\rightarrow$ mode swapping (without dislocations) $\rightarrow$ Mode B’. Relatively strong three-dimensionality is found at $Re=160{-}220$ as the wake is dominated by large-scale structure of ordered Mode A, and also at $Re\geqslant 285$ , where Mode B becomes increasingly disordered. A local reduction in three-dimensionality is observed at $Re=225{-}275$ , where the wake is dominated by finer-scale structure of a mixture of ordered Modes A and B. Corresponding variations in the vortex shedding frequency and hydrodynamic forces are also investigated.

We present a new solution for the steady boundary-layer flow over the rotating sphere that also accounts for the eruption of the boundary layer at the equator and other higher-order viscous effects. Non-parallel corrections to the local type I and type II convective instability modes of this flow are also computed as a function of spin rate. Our instability results are associated with the previously observed spiral vortices and remarkable agreement between our predictions of the number of vortices and experimental observations is found. Vortices travelling at 70 %–80 % of the local surface speed are found to be the most amplified for sufficient spin rates, also consistent with prior experimental observations.

We perform a linearized local stability analysis for short-wavelength perturbations of a circular Couette flow with a radial temperature gradient. Axisymmetric and non-axisymmetric perturbations are considered and both the thermal diffusivity and the kinematic viscosity of the fluid are taken into account. The effect of asymmetry of the heating both on centrifugally unstable flows and on the onset of instabilities of centrifugally stable flows, including flows with a Keplerian shear profile, is thoroughly investigated. It is found that an inward temperature gradient destabilizes the Rayleigh-stable flow either via Hopf bifurcation if the liquid is a very good heat conductor or via steady state bifurcation if viscosity prevails over the thermal conductance.

The present contribution is a significant extension of the work by Kelbin et al. (J. Fluid Mech., vol. 721, 2013, pp. 340–366) as a new time-dependent helical coordinate system has been introduced. For this, Lie symmetry methods have been employed such that the spatial dependence of the originally three independent variables is reduced by one and the remaining variables are: the cylindrical radius $r$ and the time-dependent helical variable $\unicode[STIX]{x1D709}=(z/\unicode[STIX]{x1D6FC}(t))+b\unicode[STIX]{x1D711}$ , $b=\text{const.}$ and time $t$ . The variables $z$ and $\unicode[STIX]{x1D711}$ are the usual cylindrical coordinates and $\unicode[STIX]{x1D6FC}(t)$ is an arbitrary function of time $t$ . Assuming $\unicode[STIX]{x1D6FC}=\text{const.}$ , we retain the classical helically symmetric case. Using this, and imposing helical invariance onto the equation of motion, leads to a helically symmetric system of Euler and Navier–Stokes equations with a time-dependent pitch $\unicode[STIX]{x1D6FC}(t)$ , which may be varied arbitrarily and which is explicitly contained in all of the latter equations. This has been conducted both for primitive variables as well as for the vorticity formulation. Hence a significantly extended set of helically invariant flows may be considered, which may be altered by an external time-dependent strain along the axis of the helix. Finally, we sought new conservation laws which can be found from the helically invariant Euler and Navier–Stokes equations derived herein. Most of these new conservation laws are considerable extensions of existing conservation laws for helical flows at a constant pitch. Interestingly enough, certain classical conservation laws do not admit extensions in the new time-dependent coordinate system.

Most marine creatures exhibit remarkable flow sensing abilities. Their task of discerning hydrodynamic cues from local sensory information is particularly challenging because it relies on local and partial measurements to accurately characterize the ambient flow. This is in contrast to classical flow characterization methods, which invariably depend on the ability of an external observer to reconstruct the flow field globally and identify its topological structures. In this paper, we develop a mathematical framework in which a local sensory array is used to identify select flow features. Our approach consists of linearizing the flow field around the sensory array and providing a frame-independent parameterization of the velocity gradient tensor which reveals both the local flow ‘type’ and ‘intensity’. We show that a simple bioinspired sensory system that measures differences in flow velocities is capable of locally characterizing the flow type and intensity. We discuss the conditions under which such flow characterization is possible. Then, to demonstrate the effectiveness of this sensory system, we apply it in the canonical problem of a circular cylinder in uniform flow. We find excellent agreement between the sensed and actual flow properties. These findings will serve to direct future research on optimal sensory layouts and dynamic deployment of sensory arrays.

Typical flows contain internal flow barriers: specialised time-moving Lagrangian entities which demarcate distinct motions. Examples include the boundary between an oceanic eddy and a nearby jet, the edge of the Antarctic circumpolar vortex or the interface between two fluids which are to be mixed together in an microfluidic assay. The ability to control the locations of these barriers in a user-specified time-varying (unsteady) way can profoundly impact fluid transport between the coherent structures which are separated by the barriers. A condition on the unsteady Eulerian velocity required to achieve this objective is explicitly derived, thereby solving an ‘inverse Lagrangian coherent structure’ problem. This is an important first step in developing flow-barrier control in realistic flows, and in providing a postprocessing tool for observational/experimental velocity data. The excellent accuracy of the method is demonstrated using the Kelvin–Stuart cats-eyes flow and the unsteady double gyre, utilising finite-time Lyapunov exponents.

For a small sphere suspended in a background fluid flow near an obstacle, we calculate the hydrodynamic force on the sphere in the direction normal to the boundary of the obstacle. Using the Lorentz reciprocal theorem, we obtain analytical expressions for the normal force in the Stokes flow limit, valid for arbitrary separations of the particle from the obstacle, both for solid obstacles and those with free surfaces. The main effect of the boundary is to produce a normal force proportional to extensional flow gradients in the vicinity of the particle. The strength of this force is greatest when the separation between the surfaces of the particle and the obstacle is small relative to the particle size. While the magnitude of the force weakens for large separations between the sphere and the obstacle (decaying quadratically with separation distance), it can significantly modify Faxén’s law even at modestly large separation distances. In addition, we find a second force contribution due to the curvature of the background flow normal to the obstacle, which is also important when the sphere is close to the obstacle. The results of the theory are of importance to the dynamics of particles in confined geometries, whether bounded by a solid obstacle, the wall of a channel or a gas bubble.

This paper presents the results of a detailed experimental investigation into the effectiveness of sinusoidal leading edge serrations on aerofoils for the reduction of the noise generated by the interaction with turbulent flow. A detailed parametric study is performed to investigate the sensitivity of the noise reductions to the serration amplitude and wavelength. The study is primarily performed on flat plates in an idealized turbulent flow, which we demonstrate captures the same behaviour as when identical serrations are introduced onto three-dimensional aerofoils. The influence on the noise reduction of the turbulence integral length scale is also studied. An optimum serration wavelength is identified whereby maximum noise reductions are obtained, corresponding to when the transverse integral length scale is approximately one-fourth the serration wavelength. This paper proves that, at the optimum serration wavelength, adjacent valley sources are excited incoherently. One of the most important findings of this paper is that, at the optimum serration wavelength, the sound power radiation from the serrated aerofoil varies inversely proportional to the Strouhal number $St_{h}=fh/U$ , where $f$ , $h$ and $U$ are frequency, serration amplitude and flow speed, respectively. A simple model is proposed to explain this behaviour. Noise reductions are observed to generally increase with increasing frequency until the frequency at which aerofoil self-noise dominates the interaction noise. Leading edge serrations are also shown to reduce aerofoil self-noise. The mechanism for this phenomenon is explored through particle image velocimetry measurements. Finally, the lift and drag of the serrated aerofoil are obtained through direct measurement and compared against the straight edge baseline aerofoil. It is shown that aerodynamic performance is not substantially degraded by the introduction of the leading edge serrations on the aerofoil.

The electrification of particles embedded in a turbulent flow may cause hazards such as spark discharges but is also exploited in several industrial applications. Nonetheless, due to its complexity and sensitivity to the initial conditions, the process of build-up of particle charge is currently not well understood. In order to gain a deeper understanding of this phenomenon, we performed fully resolved numerical simulations of particle charging. More specifically, our study concerned the charging process of particles dispersed in a turbulent channel flow at a friction Reynolds number of $Re_{\unicode[STIX]{x1D70F}}=180$ . Emphasis was placed on the analysis of the interplay between the different physical mechanisms underlying particle electrification, such as fluid turbulence, particle dynamics and particle collisions. Further, we investigated the influence of some important physical parameters. According to our simulations the charge build-up depends strongly on the particle Stokes number, $Stk$ . In particular, at small Stokes numbers, $Stk=0.2$ , the turbopheretic drift inhibits particle charging. By contrast, at moderate Stokes numbers, $Stk=2$ , and low particle number densities, the electric charge builds up but cannot escape the viscous sublayer due to limited particle migration. However, in the case of high particle number densities, the charge is transported away from the wall via inter-particle charge diffusion. A further increase to $Stk=20$ leads to strong charging and particle-bound charge transport towards the bulk of the channel.