Most of the existing theories on electrophoresis are based on the consideration of a weak applied electric field and ions as point charges, which create a mean electric potential and neglect ion–solvent interactions. These theories cannot demonstrate the dependence of electrophoretic mobility on the applied electric field (nonlinear electrophoresis), reversal in mobility with increasing ion concentration and/or surface charge density or counterion saturation in the electric double layer. In this study we consider a modified electrokinetic model to analyse nonlinear electrophoresis by taking into account the finite ion size effects and ion–ion electrostatic correlations. In this approach, the mean-field-based model is extended to capture the many-body phenomena by considering the non-local electrostatic contribution in the ion free energy functional and the ion–ion hydrodynamic steric interactions are incorporated through the volume exclusion effect in the electrochemical potential. The viscosity of the medium is considered to vary with the local ionic volume fraction. Stronger correlations for multivalent counterions create ion layering, charge density oscillation and mobility reversal. Such phenomena are captured by the present continuum model. The ion crowding attenuates the growth of the electrophoretic mobility with the electric field. At a higher range of the imposed electric field, the ion concentration in the electric double layer enhances, which modifies both the overscreening and ion crowding processes.

Experiments and numerical simulations of inertial particles in underexpanded jets are performed. The structure of the jet is controlled by varying the nozzle pressure ratio, while the influence of particles on emerging shocks and rarefaction patterns is controlled by varying the particle size and mass loading. Ultra-high-speed schlieren and Lagrangian particle tracking are used to experimentally determine the two-phase flow quantities. Three-dimensional simulations are performed using a high-order, low-dissipative discretization of the gas phase while particles are tracked individually in a Lagrangian manner. A simple two-way coupling strategy is proposed to handle interphase exchange in the vicinity of shocks. Velocity statistics of each phase are reported for a wide range of pressure ratios, particle sizes and volume fractions. An upstream shift of the Mach disk in the presence of particles reveals significant two-way coupling even at low mass loading. A semi-analytic model that predicts the extent of the Mach disk shift is presented based on a one-dimensional Fanno flow that takes into account volume displacement by particles and interphase exchange due to drag and heat transfer. The per cent shift in Mach disk is found to scale with the mass loading, nozzle pressure ratio and interphase slip velocity and inversely with the particle diameter.

Porous membranes are thin solid structures that allow the flow to pass through their tiny openings, called pores. Flow inertia may play a significant role in several filtration flows of natural and engineering interest. Here, we develop a predictive macroscopic model to describe solvent and solute flows past thin membranes for non-negligible inertia. We leverage homogenization theory to link the solvent velocity and solute concentration to the jumps of solvent stress and solute flux across the membrane. Within this framework, the membrane acts as a boundary separating two distinct fluid regions. These jump conditions rely on several coefficients, stemming from closure problems at the microscopic pore scale. Two approximations for the advective terms of Navier–Stokes and advection–diffusion equations are introduced to include inertia in the microscopic problem. The approximate inertial terms couple the micro- and macroscopic fields. Here, this coupling is solved numerically using an iterative fixed-point procedure. We compare the resulting models against full-scale simulations, with a good agreement both in terms of averaged values across the membrane and far-field values. Eventually, we develop a strategy based on unsupervised machine learning to improve the computational efficiency of the iterative procedure. The extension of homogenization towards weak-inertia flow configurations as well as the performed data-driven approximation may find application in preliminary analyses as well as optimization procedures towards the design of filtration systems, where inertia effects can be instrumental in broadening the spectrum of permeability and selectivity properties of these filters.

We report an experimental study of Rayleigh–Bénard convection of liquid metal GaInSn in a cuboid cell with an aspect ratio of 0.5 under the effect of a horizontal magnetic field. The Rayleigh number spans a range of $3.8\times 10^5 \leqslant Ra \leqslant 1.1\times 10^7$, while the magnetic field strength reaches up to 0.5 T, corresponding to a maximum Hartmann number to 2041. By combining temperature and velocity measurements, we identify several flow morphologies, including a novel cellular pattern characterized by four stacked vortices that periodically squeeze and induce velocity reversals. Based on the identified flow morphologies, we partition the entire ($Ra, Ha$) parameter space into five distinct flow regimes and systematically investigate the flow characteristics within each regime. The temperature gradient and oscillation frequency exhibit scaling relationships with the combined parameters $Ra$ and $Ha$. Notably, we observe a coupling between flow regime and global transport efficiencies, particularly in a regime dominated by the double-roll structure, which experiences a maximum 36 % decrease in heat transfer efficiency compared with the single-roll structure. The dependencies of heat and momentum transport on $Ra$ and $Ha$ follow scaling laws as $Nu \sim (Ha^{-2/3}RaPr^{-1})^{3/5}$ and $Re \sim (Ha^{-1}RaPr^{-1})^{4/3}$, respectively.

We employ direct numerical simulations to investigate the heat transfer and flow structures in turbulent Rayleigh–Bénard convection in both cylindrical cells and laterally periodic domains, spanning an unprecedentedly wide range of aspect ratios $0.075 \leqslant \varGamma \leqslant 32$. We focus on Prandtl number ${Pr}=1$ and Rayleigh numbers ${{Ra}}=2\times 10^7$ and ${{Ra}}=10^8$. In both cases, with increasing aspect ratio, the heat transfer first increases, then reaches a maximum (which is more pronounced for the cylindrical case due to confinement effects), and then slightly goes down again before it finally saturates at the large aspect ratio limit, which is achieved already at $\varGamma \approx 4$. Already for $\varGamma \gtrsim 0.75$, the heat transfers in both cylindrical and laterally periodic domains become identical. The large-$\varGamma$ limit for the volume-integrated Reynolds number and the boundary layer thicknesses are also reached at $\varGamma \approx 4$. However, while the integral flow properties converge at $\varGamma \approx 4$, the confinement of a cylindrical domain impacts the temperature and velocity variance distributions up to $\varGamma \approx 16$, as thermal superstructures cannot form close to the sidewall.

The dominant mode instability in hypersonic boundary-layer transition is the so-called second-mode instability, which has a peculiar nature strongly coupled with thermoacoustic phenomena. In linear stability theory, the unstable wave is associated with one of the two eigenvalues that originate from the acoustic branches, referred to as slow and fast modes. Interestingly, the unstable mode (slow or fast) reaches its maximum amplification as the other mode (fast or slow) attains a minimum. The phase velocity of the two modes is then very close, and this phenomenon is called synchronization. The aim of the present study is to unravel the physical mechanism that explains the second-mode growth. To that aim, second-order nonlinear equations are written for the disturbances given by linear stability. In this framework, entropy, kinetic energy and temperature energy budgets are obtained up to second order. The budgets are scrutinized for various Mach numbers and for adiabatic and cold-wall thermal conditions. Perturbation entropy budgets clearly show the process is a reversible one. An energy exchange between kinetic energy and temperature energy of the weakly nonlinear modes is driven by pressure–dilatation terms. As underlined in previous studies, the unstable mode experiences an alternate heating and cooling near the wall, which is shown to be a rather nonlinear process. The change in fluctuating thermal energy in the form of a dilatational wave is sustained by pumping disturbance kinetic energy through the pressure–dilatation term, the direction of the conversion being driven by the relative phase between pressure and dilatation. This process is similar for the slow and fast modes, the unstable mode being amplified and the other being damped. No change in the process has been noted at the location of the synchronization, suggesting that the modes have the same nature but evolve independently.