Exact coherent structures (ECS), unstable three-dimensional solutions of the Navier–Stokes equations, play a fundamental role in transitional and turbulent wall flows. Dempsey et al. (J. Fluid Mech., vol. 791, 2016, pp. 97–121) demonstrate that at large Reynolds number reduced equations can be derived that simplify the computation and facilitate mechanistic understanding of these solutions. Their analysis shows that ECS in plane Poiseuille flow can be sustained by a novel inner–outer interaction between oblique near-wall Tollmien–Schlichting waves and interior streamwise vortices.

In the analysis of velocity fields in turbulent boundary layers, the traditional Reynolds decomposition is universally employed to calculate the fluctuating component of streamwise velocity. Here, we demonstrate the perils of such a determination of the fluctuating velocity in the context of structural analysis of turbulence when applied in the outer region where the flow is intermittently turbulent at a given wall distance. A new decomposition is postulated that ensures non-turbulent regions in the flow do not contaminate the fluctuating velocity components in the turbulent regions. Through this new decomposition, some of the typical statistics concerning the scale and structure of turbulent boundary layers are revisited.

Oceanic waves registered by satellite observations often have curvilinear fronts and propagate over various currents. In this paper we study long linear and weakly nonlinear ring waves in a stratified fluid in the presence of a depth-dependent horizontal shear flow. It is shown that, despite the clashing geometries of the waves and the shear flow, there exists a linear modal decomposition (different from the known decomposition in Cartesian geometry), which can be used to describe distortion of the wavefronts of surface and internal waves, and systematically derive a $2+1$-dimensional cylindrical Korteweg–de Vries-type equation for the amplitudes of the waves. The general theory is applied to the case of the waves in a two-layer fluid with a piecewise-constant current, with an emphasis on the effect of the shear flow on the geometry of the wavefronts. The distortion of the wavefronts is described by the singular solution (envelope of the general solution) of the nonlinear first-order differential equation, constituting generalisation of the dispersion relation in this curvilinear geometry. There exists a striking difference in the shapes of the wavefronts of surface and interfacial waves propagating over the same shear flow.

The dissolution process of small (initial (equivalent) radius $R_{0}<1$  mm) long-chain alcohol (of various types) sessile droplets in water is studied, disentangling diffusive and convective contributions. The latter can arise for high solubilities of the alcohol, as the density of the alcohol–water mixture is then considerably less than that of pure water, giving rise to buoyancy-driven convection. The convective flow around the droplets is measured, using micro-particle image velocimetry (${\rm\mu}$PIV) and the schlieren technique. When non-dimensionalizing the system, we find a universal $Sh\sim Ra^{1/4}$ scaling relation for all alcohols (of different solubilities) and all droplets in the convective regime. Here $Sh$ is the Sherwood number (dimensionless mass flux) and $Ra$ is the Rayleigh number (dimensionless density difference between clean and alcohol-saturated water). This scaling implies the scaling relation ${\it\tau}_{c}\propto R_{0}^{5/4}$ of the convective dissolution time ${\it\tau}_{c}$, which is found to agree with experimental data. We show that in the convective regime the plume Reynolds number (the dimensionless velocity) of the detaching alcohol-saturated plume follows $Re_{p}\sim Sc^{-1}Ra^{5/8}$, which is confirmed by the ${\rm\mu}$PIV data. Here, $Sc$ is the Schmidt number. The convective regime exists when $Ra>Ra_{t}$, where $Ra_{t}=12$ is the transition $Ra$ number as extracted from the data. For $Ra\leqslant Ra_{t}$ and smaller, convective transport is progressively overtaken by diffusion and the above scaling relations break down.

This paper is concerned with the rather broad issue of the impact of abrupt changes (such as isolated roughness, gaps and local suctions) on boundary-layer transition. To fix the idea, we consider the influence of a two-dimensional localized hump (or indentation) on an oncoming Tollmien–Schlichting (T–S) wave. We show that when the length scale of the former is comparable with the characteristic wavelength of the latter, the key physical mechanism to affect transition is through scattering of T–S waves by the roughness-induced mean-flow distortion. An appropriate mathematical theory, consisting of the boundary-value problem governing the local scattering, is formulated based on triple deck formalism. The transmission coefficient, defined as the ratio of the amplitude of the T–S wave downstream of the roughness to that upstream, serves to characterize the impact on transition. The transmission coefficient appears as the eigenvalue of the discretized boundary-value problem. The latter is solved numerically, and the dependence of the eigenvalue on the height and width of the roughness and the frequency of the T–S wave is investigated. For a roughness element without causing separation, the transmission coefficient is found to be approximately 1.5 for typical frequencies, indicating a moderate but appreciable destabilizing effect. For a roughness causing incipient separation, the transmission coefficient can be as large as $O(10)$, suggesting that immediate transition may take place at the roughness site. A roughness element with a fixed height produces the strongest impact when its width is comparable with the T–S wavelength, in which case the traditional linear stability theory is invalid. The latter however holds approximately when the roughness width is sufficiently large. By studying the two hump case, a criterion when two roughness elements can be regarded as being isolated is suggested. The transmission coefficient can be converted to an equivalent $N$-factor increment, by making use of which the $\text{e}^{N}$-method can be extended to predict transition in the presence of multiple roughness elements.

Flow generated by an oscillating sphere in a quiescent fluid is a classical problem in fluid mechanics whose solution is used ubiquitously. Miniaturisation of mechanical devices to small scales and their operation at high frequencies in fluid, which is common in modern nanomechanical systems, can preclude the use of the unsteady Stokes equation for continuum flow. Here, we explore the combined effects of gas rarefaction and unsteady motion of a sphere, within the framework of the unsteady linearised Boltzmann–BGK (Bhatnagar–Gross–Krook) equation. This equation is solved using the method of characteristics, and the resulting solution is valid for any oscillation frequency and arbitrary degrees of gas rarefaction. The resulting force provides the non-continuum counterpart to the (continuum) unsteady Stokes drag on a sphere. In contrast to the Stokes solution, where the flow is isothermal, non-continuum effects lead to a temperature jump at the sphere surface and non-isothermal flow. Unsteady effects and heat transport are found to mix strongly, leading to marked differences relative to the steady case. The solution to this canonical flow problem is expected to be of significant practical value in many applications, including the optical trapping of nanoparticles and the design and application of nanoelectromechanical systems. It also provides a benchmark for computational and approximate methods of solution for the Boltzmann equation.

We consider convective instability in a deep porous medium cooled from above with a linearised thermal exchange at the upper surface, thus determining the impact of using a Robin boundary condition, in contrast to previous studies using a Dirichlet boundary condition. With the linearised surface exchange, the thermal flux out of the porous layer depends linearly on the temperature difference between the effective temperature of a heat sink at the upper boundary and the temperature at the surface of the porous layer. The rate of this exchange is characterised by a dimensionless Biot number, $\mathit{Bi}$, determined by the effective thermal conductivity of exchange with the heat sink relative to the physical thermal conductivity of the porous layer. For a given temperature difference between the heat sink at the upper boundary and deep in the porous medium, we find that imperfectly cooled layers with finite Biot numbers are more stable to convective instabilities than perfectly cooled layers which have large, effectively infinite Biot numbers. Two regimes of behaviour were determined with contrasting stability behaviour and characteristic scales. When the Biot number is large the near-perfect heat transfer produces small corrections of order $1/\mathit{Bi}$ to the perfectly conducting behaviour found when the Biot number is infinite. In the insulating limit as the Biot number approaches zero, a different behaviour was found with significantly larger scales for the critical wavelength and depth of convection both scaling proportional to $1/\sqrt{\mathit{Bi}}$.

The steady-state nearly resonant water waves with time-independent spectrum in deep water are obtained from the full wave equations for inviscid, incompressible gravity waves in the absence of surface tension by means of a analytic approximation approach based on the homotopy analysis method (HAM). Our strategy is to mathematically transfer the steady-state nearly resonant wave problem into the steady-state exactly resonant ones. By means of choosing a generalized auxiliary linear operator that is a little different from the linear part of the original wave equations, the small divisor, which is unavoidable for nearly resonant waves in the frame of perturbation methods, is avoided, or moved far away from low wave frequency to rather high wave frequency with physically negligible wave energy. It is found that the steady-state nearly resonant waves have nothing fundamentally different from the steady-state exactly resonant ones, from physical and numerical viewpoints. In addition, the validity of this HAM-based analytic approximation approach for the full wave equations in deep water is numerically verified by means of the Zakharov’s equation. A thought experiment is discussed, which suggests that the essence of the so-called ‘wave resonance’ should be reconsidered carefully from both of physical and mathematical viewpoints.

For the fabrication of P-OLED displays, using inkjet printing, it is important to control the final shape resulting from evaporation of droplets containing polymer. Due to peripheral pinning and consequent outward capillary flow, a ring-like final shape is typically observed. This is often undesirable, with a spatially uniform film usually required. Several experimental studies have shown that binary liquid inks can prevent ring formation. There is no consensus of opinion on the mechanism behind this improvement. We have developed a model for the drying of thin, binary liquid droplets, based on thin-film lubrication theory, and we solve the governing equations to predict the final shape. White-light interferometry experiments are conducted to verify the findings. In addition, we present the results of a linear stability analysis that identifies the onset of an instability driven by a difference in surface tension. If the more volatile liquid is more abundant, an instability becomes increasingly likely.

Kolmogorov introduced dissipative scales based on the mean dissipation $\langle {\it\varepsilon}\rangle$ and the viscosity ${\it\nu}$, namely the Kolmogorov length ${\it\eta}=({\it\nu}^{3}/\langle {\it\varepsilon}\rangle )^{1/4}$ and the velocity $u_{{\it\eta}}=({\it\nu}\langle {\it\varepsilon}\rangle )^{1/4}$. However, the existence of smaller scales has been discussed in the literature based on phenomenological intermittency models. Here, we introduce exact dissipative scales for the even-order longitudinal structure functions. The derivation is based on exact relations between even-order moments of the longitudinal velocity gradient $(\partial u_{1}/\partial x_{1})^{2m}$ and the dissipation $\langle {\it\varepsilon}^{m}\rangle$. We then find a new length scale ${\it\eta}_{C,m}=({\it\nu}^{3}/\langle {\it\varepsilon}^{m/2}\rangle ^{2/m})^{1/4}$ and $u_{C,m}=({\it\nu}\langle {\it\varepsilon}^{m/2}\rangle ^{2/m})^{1/4}$, i.e. the dissipative scales depend rather on the moments of the dissipation $\langle {\it\varepsilon}^{m/2}\rangle$ and thus the full probability density function (p.d.f.) $P({\it\varepsilon})$ instead of powers of the mean $\langle {\it\varepsilon}\rangle ^{m/2}$. The results presented here are exact for longitudinal even-ordered structure functions under the assumptions of (local) isotropy, (local) homogeneity and incompressibility, and we find them to hold empirically also for the mixed and transverse as well as odd orders. We use direct numerical simulations (DNS) with Reynolds numbers from $Re_{{\it\lambda}}=88$ up to $Re_{{\it\lambda}}=754$ to compare the different scalings. We find that indeed $P({\it\varepsilon})$ or, more precisely, the scaling of $\langle {\it\varepsilon}^{m/2}\rangle /\langle {\it\varepsilon}\rangle ^{m/2}$ as a function of the Reynolds number is a key parameter, as it determines the ratio ${\it\eta}_{C,m}/{\it\eta}$ as well as the scaling of the moments of the velocity gradient p.d.f. As ${\it\eta}_{C,m}$ is smaller than ${\it\eta}$, this leads to a modification of the estimate of grid points required for DNS.

The force-driven Poiseuille flow of dense gases between two parallel plates is investigated through the numerical solution of the generalized Enskog equation for two-dimensional hard discs. We focus on the competing effects of the mean free path ${\it\lambda}$, the channel width $L$ and the disc diameter ${\it\sigma}$. For elastic collisions between hard discs, the normalized mass flow rate in the hydrodynamic limit increases with $L/{\it\sigma}$ for a fixed Knudsen number (defined as $Kn={\it\lambda}/L$), but is always smaller than that predicted by the Boltzmann equation. Also, for a fixed $L/{\it\sigma}$, the mass flow rate in the hydrodynamic flow regime is not a monotonically decreasing function of $Kn$ but has a maximum when the solid fraction is approximately 0.3. Under ultra-tight confinement, the famous Knudsen minimum disappears, and the mass flow rate increases with $Kn$, and is larger than that predicted by the Boltzmann equation in the free-molecular flow regime; for a fixed $Kn$, the smaller $L/{\it\sigma}$ is, the larger the mass flow rate. In the transitional flow regime, however, the variation of the mass flow rate with $L/{\it\sigma}$ is not monotonic for a fixed $Kn$: the minimum mass flow rate occurs at $L/{\it\sigma}\approx 2{-}3$. For inelastic collisions, the energy dissipation between the hard discs always enhances the mass flow rate. Anomalous slip velocity is also found, which decreases with increasing Knudsen number. The mechanism for these exotic behaviours is analysed.

The effect of entrained air turbulence on dispersion of droplets (with Stokes number based on the Kolmogorov time scale, $St_{{\it\eta}}$, of the order of 1) in a polydispersed spray is experimentally studied through simultaneous and planar measurements of droplet size, velocity and gas flow velocity (Hardalupas et al., Exp. Fluids, vol. 49, 2010, pp. 417–434). The preferential accumulation of droplets at various measurement locations in the spray was examined by two independent methods viz. counting droplets on images by dividing the image in to boxes of different sizes, and by estimating the radial distribution function (RDF). The dimension of droplet clusters (obtained by both approaches) was of the order of Kolmogorov’s length scale of the fluid flow, implying the significant influence of viscous scales of the fluid flow on cluster formation. The RDF of different size classes indicated an increase in cluster dimension for larger droplets (higher $St_{{\it\eta}}$). The length scales of droplet clusters increased towards the outer spray regions, where the gravitational influence on droplets is stronger compared to the central spray locations. The correlation between fluctuations of droplet concentration and droplet and gas velocities were estimated and found to be negative near the spray edge, while it was close to zero at other locations. The probability density function of slip between fluctuating droplet velocity and gas velocity ‘seen’ by the droplets signified presence of considerable instantaneous slip velocity, which is crucial for droplet–gas momentum exchange. In order to investigate different mechanisms of turbulence modulation of the carrier phase, the three correlation terms in the turbulent kinetic energy equation for particle-laden flows (Chen & Wood, Can. J. Chem. Engng, vol. 65, 1985, pp. 349–360) are evaluated conditional on droplet size classes. Based on the comparison of the correlation terms, it is recognized that although the interphase energy transfer due to fluctuations of droplet concentration is low compared to the energy exchange only due to droplet drag (the magnitude of which is controlled by average droplet mass loading), the former cannot be considered negligible, and should be accounted in two phase flow modelling.

Nonlinear interactions between sea waves and the sea bottom are a major mechanism for energy transfer between the different wave frequencies in the near-shore region. Nevertheless, it is difficult to account for this phenomenon in stochastic wave forecasting models due to its mathematical complexity, which mostly consists of computing either the bispectral evolution or non-local shoaling coefficients. In this work, quasi-two-dimensional stochastic energy evolution equations are derived for dispersive water waves up to quadratic nonlinearity. The bispectral evolution equations are formulated using stochastic closure. They are solved analytically and substituted into the energy evolution equations to construct a stochastic model with non-local shoaling coefficients, which includes nonlinear dissipative effects and slow time evolution. The nonlinear shoaling mechanism is investigated and shown to present two different behaviour types. The first consists of a rapidly oscillating behaviour transferring energy back and forth between wave harmonics in deep water. Owing to the contribution of bottom components for closing the class III Bragg resonance conditions, this behaviour includes mean energy transfer when waves reach intermediate water depths. The second behaviour relates to one-dimensional shoaling effects in shallow water depths. In contrast to the behaviour in intermediate water depths, it is shown that the nonlinear shoaling coefficients refrain from their oscillatory nature while presenting an exponential energy transfer. This is explained through the one-dimensional satisfaction of the Bragg resonance conditions by wave triads due to the non-dispersive propagation in this region even without depth changes. The energy evolution model is localized using a matching approach to account for both these behaviour types. The model is evaluated with respect to deterministic ensembles, field measurements and laboratory experiments while performing well in modelling monochromatic superharmonic self-interactions and infra-gravity wave generation from bichromatic waves and a realistic wave spectrum evolution. This lays physical and mathematical grounds for the validity of unexplained simplifications in former works and the capability to construct a formulation that consistently accounts for nonlinear energy transfers from deep to shallow water.

We examine the sensitivity of Saffman–Taylor fingers to controlled variations in channel depth by investigating the effects of centred, rectangular occlusions in Hele-Shaw channels. For large occlusions, the geometry is known to support symmetric, asymmetric and oscillatory propagation states when air displaces a more viscous fluid from within the channel. A previously developed depth-averaged model is found to be in quantitative agreement with laboratory experiments once the aspect ratio (width/height) of the tube’s cross-section reaches a value of 40. We find that the multiplicity of solutions at finite occlusion heights arises through interactions of the single stable and multiple unstable solutions already present in the absence of the occlusion: the classic Saffman–Taylor viscous fingering problem. The sequence of interactions that occurs with increasing occlusion height is the same for all aspect ratios investigated, but the occlusion height required for each interaction decreases with increasing aspect ratio. Thus, the system becomes more sensitive as the aspect ratio increases in the sense that multiple solutions are provoked for smaller relative depth changes. We estimate that the required depth changes become of the same order as the typical roughnesses of the experimental system ($1~{\rm\mu}\text{m}$) for aspect ratios beyond 155, which we conjecture underlies the extreme sensitivity of experiments conducted in such Hele-Shaw channels.

We investigate the statistical properties, based on numerical simulations and analytical calculations, of a recently proposed stochastic model for the velocity field (Chevillard et al., Europhys. Lett., vol. 89, 2010, 54002) of an incompressible, homogeneous, isotropic and fully developed turbulent flow. A key step in the construction of this model is the introduction of some aspects of the vorticity stretching mechanism that governs the dynamics of fluid particles along their trajectories. An additional further phenomenological step aimed at including the long range correlated nature of turbulence makes this model dependent on a single free parameter, ${\it\gamma}$, that can be estimated from experimental measurements. We confirm the realism of the model regarding the geometry of the velocity gradient tensor, the power-law behaviour of the moments of velocity increments (i.e. the structure functions) including the intermittent corrections and the existence of energy transfer across scales. We quantify the dependence of these basic properties of turbulent flows on the free parameter ${\it\gamma}$ and derive analytically the spectrum of exponents of the structure functions in a simplified non-dissipative case. A perturbative expansion in power of ${\it\gamma}$ shows that energy transfer, at leading order, indeed take place, justifying the dissipative nature of this random field.

We investigate numerically the nonlinear interactions between hetons. Hetons are baroclinic structures consisting of two vortices of opposite sign lying at different depths. Hetons are long-lived. They most often translate (they can sometimes rotate) and therefore they can noticeably contribute to the transport of scalar properties in the oceans. Heton interactions can interrupt this translation and thus this transport, by inducing a reconfiguration of interacting hetons into more complex baroclinic multipoles. More specifically, we study here the general case of two hetons, which collide with an offset between their translation axes. For this purpose, we use the point vortex theory, the ellipsoidal vortex model and direct simulations in the three-dimensional quasi-geostrophic contour surgery model. More specifically, this paper shows that there are in general three regimes for the interaction. For small horizontal offsets between the hetons, their vortices recombine as same-depth dipoles which escape at an angle. The angle depends in particular on the horizontal offset. It is a right angle for no offset, and the angle is shallower for small but finite offsets. The second limiting regime is for large horizontal offsets where the two hetons remain the same hetonic structures but are deflected by the weaker mutual interaction. Finally, the intermediate regime is for moderate offsets. This is the regime where the formation of a metastable quadrupole is possible. The formation of this quadrupole greatly restrains transport. Indeed, it constrains the vortices to reside in a closed area. It is shown that the formation of such structures is enhanced by the quasi-periodic deformation of the vortices. Indeed, these structures are nearly unobtainable for singular vortices (point vortices) but may be obtained using deformable, finite-core vortices.

Tornadoes are one type of violent flow phenomenon and occur in many places in the world. There are many research methods that aim to reduce the loss of human lives and material damage caused by tornadoes. One effective method is numerical simulation such as that in Ishihara et al. (J. Wind Engng Ind. Aerodyn., vol. 99, 2011, pp. 239–248). The swirling structure of the Navier–Stokes flow is significant for both the mathematical analysis and numerical simulations of tornadoes. In this paper, we try to clarify the swirling structure. More precisely, we performed numerical computations on axisymmetric Navier–Stokes flows with a no-slip flat boundary. We compared a hyperbolic flow with swirl and one without swirl, and observed that the following phenomenon occurs only in the swirl case: the distance between the point with the maximum magnitude of velocity $|\boldsymbol{v}|$ and the $z$-axis changed drastically at a specific time (which we call the turning point). Besides, an ‘increasing velocity phenomenon’ occurred near the boundary, and the maximum value of $|\boldsymbol{v}|$ was obtained near the axis of symmetry and the boundary when the time was close to the turning point in the swirl case.

Debris and pyroclastic flows often have bouldery flow fronts, which act as a natural dam resisting further advance. Counter intuitively, these resistive fronts can lead to enhanced run-out, because they can be shouldered aside to form static levees that self-channelise the flow. At the heart of this behaviour is the inherent process of size segregation, with different sized particles readily separating into distinct vertical layers through a combination of kinetic sieving and squeeze expulsion. The result is an upward coarsening of the size distribution with the largest grains collecting at the top of the flow, where the flow velocity is greatest, allowing them to be preferentially transported to the front. Here, the large grains may be overrun, resegregated towards the surface and recirculated before being shouldered aside into lateral levees. A key element of this recirculation mechanism is the formation of a breaking size-segregation wave, which allows large particles that have been overrun to rise up into the faster moving parts of the flow as small particles are sheared over the top. Observations from experiments and discrete particle simulations in a moving-bed flume indicate that, whilst most large particles recirculate quickly at the front, a few recirculate very slowly through regions of many small particles at the rear. This behaviour is modelled in this paper using asymmetric segregation flux functions. Exact non-diffuse solutions are derived for the steady wave structure using the method of characteristics with a cubic segregation flux. Three different structures emerge, dependent on the degree of asymmetry and the non-convexity of the segregation flux function. In particular, a novel ‘lens-tail’ solution is found for segregation fluxes that have a large amount of non-convexity, with an additional expansion fan and compression wave forming a ‘tail’ upstream of the ‘lens’ region. Analysis of exact solutions for the particle motion shows that the large particle motion through the ‘lens-tail’ is fundamentally different to the classical ‘lens’ solutions. A few large particles starting near the bottom of the breaking wave pass through the ‘tail’, where they travel in a region of many small particles with a very small vertical velocity, and take significantly longer to recirculate.

We present new observations from an experimental investigation of the classical problem of the crown splash and sealing phenomena observed during the impact of spheres onto quiescent liquid pools. In the experiments, a 6 m tall vacuum chamber was used to provide the required ambient conditions from atmospheric pressure down to $1/16\text{th}$ of an atmosphere, whilst high-speed videography was exploited to focus primarily on the above-surface crown formation and ensuing dynamics, paying particular attention to the moments just prior to the surface seal. In doing so, we have observed a buckling-type azimuthal instability of the crown. This instability is characterised by vertical striations along the crown, between which thin films form that are more susceptible to the air flow and thus are drawn into the closing cavity, where they atomize to form a fine spray within the cavity. To elucidate to the primary mechanisms and forces at play, we varied the sphere diameter, liquid properties and ambient pressure. Furthermore, a comparison between the entry of room-temperature spheres, where the contact line pins around the equator, and Leidenfrost spheres (i.e. an immersed superheated sphere encompassed by a vapour layer), where there is no contact line, indicates that the buckling instability appears in all crown sealing events, but is intensified by the presence of a pinned contact line.

Invariant solutions of shear flows have recently been extended from spatially periodic solutions in minimal flow units to spatially localized solutions on extended domains. One set of spanwise-localized solutions of plane Couette flow exhibits homoclinic snaking, a process by which steady-state solutions grow additional structure smoothly at their fronts when continued parametrically. Homoclinic snaking is well understood mathematically in the context of the one-dimensional Swift–Hohenberg equation. Consequently, the snaking solutions of plane Couette flow form a promising connection between the largely phenomenological study of laminar–turbulent patterns in viscous shear flows and the mathematically well-developed field of pattern-formation theory. In this paper we present a numerical study of the snaking solutions of plane Couette flow, generalizing beyond the fixed streamwise wavelength of previous studies. We find a number of new solution features, including bending, skewing and finite-size effects. We establish the parameter regions over which snaking occurs and show that the finite-size effects of the travelling wave solution are due to a coupling between its fronts and interior that results from its shift-reflect symmetry. A new winding solution of plane Couette flow is derived from a strongly skewed localized equilibrium.