Explore the vast range of titles available.

The mass transfer in nanopores is understood respectively from the multiscale flow regime and the nanoscale non-continuum flow regime according to the size of the nanopore. If no interfacial slippage occurs, the mass flow rate through the nanopore is normally far smaller than the classical Hagen–Poiseuille equation calculation especially for small nanopores owing to the fluid-pore wall interaction which results in the significant effects of the viscosity enhancement and the non-continuum property of the very thin adhering layer in the nanopore. If the interfacial slippage occurs, in the multiscale flow regime, the adhering layer-pore wall interfacial slippage is advantageous over the adhering layer-continuum fluid interfacial slippage especially for small nanopores because of generating much greater flow rates through the nanopore. A hydrophobic nanopore wall is thus preferential. In the nanoscale non-continuum flow regime, the wall slippage more easily occurs, and its effect is determined by both the power loss on the nanopore and the intrinsic parameter of the nanopore.

A method for solving the Boltzmann equation is presented that makes it possible to calculate gas flows in the continuum flow regime described by the Navier–Stokes equations. Progress into the region of continuum flows was achieved by applying the conservative projection method for calculating the Boltzmann collision integral, which preserves the leading term of the Enskog–Chapman asymptotics. Optimization of this method that made it possible to considerably decrease the amount of computations is described. Examples of the longitudinal subsonic flow around a flat plate for the case of the Knudsen numbers are discussed.

The authors study small disturbances to the periodic, plane Couette flow in the 3D incompressible Navier-Stokes equations at high Reynolds number Re. They prove that for sufficiently regular initial data of size $\\epsilon \\leq c_0\\mathbf {Re}^-1$ for some universal $c_0 > 0$, the solution is global, remains within $O(c_0)$ of the Couette flow in $L^2$, and returns to the Couette flow as $t \\rightarrow \\infty $. For times $t \\gtrsim \\mathbf {Re}^1/3$, the streamwise dependence is damped by a mixing-enhanced dissipation effect and the solution is rapidly attracted to the class of \"2.5 dimensional\" streamwise-independent solutions referred to as streaks.

Selective solvent and solute transport across nanopores is fundamental to membrane separations, yet it remains poorly understood, especially for non-aqueous systems. Here, we design a chemically robust nanoporous graphene membrane and study molecular transport in various organic liquids under subnanometre confinement. We show that the nature of the solvent can modulate solute diffusion across graphene nanopores, and that breakdown of continuum flow occurs when pore size approaches the solvent’s smallest molecular cross-section. By holistically engineering membrane support, modelling pore creation and defect management, high rejection and ultrafast organic solvent nanofiltration of dye molecules and separation of hexane isomers are achieved. The membranes exhibit stable fluxes across a range of solvents, consistent with flow across rigid pores whose size is independent of the solvent. These results demonstrate that nanoporous graphene is a rich materials system for controlling subcontinuum flow that could enable new membranes for a range of challenging separation needs. A study of molecular transport in various organic liquids under subnanometre confinement shows that the nature of the solvent can modulate solute diffusion across graphene nanopores, and that breakdown of continuum flow occurs when pore size approaches the solvent’s smallest molecular cross-section.

The thermal conductivity of a molecular gas consists of the translational and internal parts. Although in continuum flows the total thermal conductivity itself is adequate to describe the heat transfer, in rarefied gas flows they need to be modelled separately, according to the relaxation rates of translational and internal heat fluxes in an homogeneous system. This paper is dedicated to quantifying how these relaxation rates affect rarefied gas dynamics. The kinetic model of Wu et al. (J. Fluid Mech., vol. 763, 2015, pp. 24–50) is adapted to recover the relaxation of heat fluxes, which is validated by the direct simulation Monte Carlo method. Then the model of Wu et al., which has the freedom to adjust the relaxation rates, is used to investigate the rate effects of thermal relaxation in problems such as the normal shock wave, creep flow driven by Maxwell's demon and thermal transpiration. It is found that the relaxation rates of heat flux affect rarefied gas flows significantly, even when the total thermal conductivity is fixed.

Following the stability analysis method in classic fluid dynamics, a linear stability equation (LSE) suitable for rarefied flows is derived based on the Bhatnagar–Gross–Krook (BGK) equation. The global method and singular value decomposition method are used for modal and non-modal analysis, respectively. This approach is validated by results obtained from Navier–Stokes (NS) equations. The modal analysis shows that LSEs based on NS equations (NS-LSEs) begin to fail when the Knudsen number ($Kn$) increases past $\\sim$0.01, regardless of whether a slip model is used. When $Kn\\geq 0.01$, the growth rate of the least stable mode is generally underestimated by the NS-LSEs. Under a fixed wavenumber, the pattern (travelling or standing wave) of the least stable mode changes with $Kn$; when the mode presents the same pattern, the growth rate decreases almost linearly with increasing $Kn$; otherwise, rarefaction effects may not stabilize the flow. The characteristic lengths of the different modes are different, and the single-scale classic stability analysis method cannot predict multiple modes accurately, even when combined with a slip model and even for continuum flow. However, non-modal analysis shows that this error does not affect the transient growth because modes with small growth rates offer little contribution to the transient growth. In rarefied flow, as long as the Mach number ($Ma$) is large enough, transient growth will occur in some wavenumber ranges. The rarefaction effect plays a stabilizing role in transient growth. The NS-LSEs-based method always overestimates the maximum transient growth.

In geoenergy applications, mudrocks prevent fluids to leak from temporary (H 2 , CH 4 ) or permanent (CO 2 , radioactive waste) storage/disposal sites and serve as a source and reservoir for unconventional oil and gas. Understanding transport properties integrated with dominant fluid flow mechanisms in mudrocks is essential to better predict the performance of mudrocks within these applications. In this study, small-angle neutron scattering (SANS) experiments were conducted on 71 samples from 13 different sets of mudrocks across the globe to capture the pore structure of nearly the full pore size spectrum (2 nm–5 μm). We develop fractal models to predict transport properties (permeability and diffusivity) based on the SANS-derived pore size distributions. The results indicate that transport phenomena in mudrocks are intrinsically pore size-dependent. Depending on hydrostatic pore pressures, transition flow develops in micropores, slip flow in meso- and macropores, and continuum flow in larger macropores. Fluid flow regimes progress towards larger pore sizes during reservoir depletion or smaller pore sizes during fluid storage, so when pressure is decreased or increased, respectively. Capturing the heterogeneity of mudrocks by considering fractal dimension and tortuosity fractal dimension for defined pore size ranges, fractal models integrate apparent permeability with slip flow, Darcy permeability with continuum flow, and gas diffusivity with diffusion flow in the matrix. This new model of pore size-dependent transport and integrated transport properties using fractal models yields a systematic approach that can also inform multiscale multi-physics models to better understand fluid flow and transport phenomena in mudrocks on the reservoir and basin scale.

The characterization and management of karst systems is a challenging task due to their inherently heterogeneous nature and vulnerability with respect to contamination. Highly conductive features of the vadose zone (e.g., dissolution shafts and faults) induce flow channeling and preferential flow. This complicates any efforts to simulate rapid recharge dynamics in deep porous‐fractured vadose zones in the context of flood and contamination risk assessment. Therefore, a strong need for numerical modeling strategies arises that employ conceptually sound formulations of these dynamics based on physical processes. Here, we present a novel modeling strategy by extending the numerical discrete‐continuum flow model MODFLOW‐CFPv2 to allow the simultaneous computation of diffuse fluxes and film‐flow in the vadose zone, thus simulating infiltration via preferential pathways. We conduct a global sensitivity analysis of a synthetic karst system that addresses the importance of including such processes in karst modeling. While event‐averaged sensitivities are in alignment with commonly observed dominance of the phreatic zone properties, results of time‐dependent sensitivities suggest that during strong infiltration events the consideration of film‐flow and its controlling parameters, that is, the fracture facial area density and an applied upper threshold for its activation, can become important. Our distributed numerical method assists in the development of karst modeling strategies where a sufficiently large and developed vadose zone offers the capacity for preferential flow that may not be accurately reproduced by most bulk‐effective methods. Hence, it benefits the unique characterization of such systems and can be easily implemented in existing workflows such as CFPv2.

Aerothermodynamic characteristics of a sphere in a subsonic flow are calculated over a broad range of gas rarefaction by the direct simulation Monte Carlo method based on ab initio interatomic potentials and Cercignani–Lampis surface scattering kernel. Calculations of the drag and average energy transfer coefficients are performed for various noble gases in the range of Mach number from 0.1 to 1. The obtained results point out that the influence of the interatomic potential is weak in subsonic flows. A comparison of the present results with a linear theory shows that the numerical solutions at Mach number equal to 0.1 are close to those obtained from the linearized kinetic equation in the transitional and free-molecular regimes. In the near-continuum flow regime, the difference between the present solution and the linear theory is significant. To reveal the effects of the gas–surface accommodation, a few sets of the tangential momentum and normal energy accommodation coefficients are considered in simulations. It is shown that the effect of the accommodation coefficients on the sphere drag is not trivial, and, for non-diffuse scattering, the drag coefficient can be either larger or smaller than that for diffuse scattering. The effect of the sphere temperature is also investigated and the calculated values of the average energy transfer coefficient are used to find the Stanton number, recovery factor and adiabatic surface temperature. The numerical results for the sphere drag and energy transfer are compared with the semi-empirical fitting equations known from the literature.

This paper presents a comprehensive investigation into flow past a circular cylinder where compressibility and rarefaction effects play an important role. The study focuses on steady subsonic flow in the Reynolds-number range 0.1–45. Rarefaction, or non-equilibrium, effects in the slip and early transition regime are accounted for using the method of moments and results are compared to data from kinetic theory obtained from the direct simulation Monte Carlo method. Solutions obtained for incompressible continuum flow serve as a baseline to examine non-equilibrium effects on the flow features. For creeping flow, where the Reynolds number is less than unity, the drag coefficient predicted by the moment equations is in good agreement with kinetic theory for Knudsen numbers less than one. When flow separation occurs, we show that the effects of rarefaction and velocity slip delay flow separation and will reduce the size of the vortices downstream of the cylinder. When the Knudsen number is above 0.028, the vortex length shows an initial increase with the Reynolds number, as observed in the standard no-slip continuum regime. However, once the Reynolds number exceeds a critical value, the size of the downstream vortices decreases with increasing Reynolds number until they disappear. An existence criterion, which identifies the limits for the presence of the vortices, is proposed. The flow physics around the cylinder is further analysed in terms of velocity slip, pressure and skin friction coefficients, which highlights that viscous, rarefaction and compressibility effects all play a complex role. We also show that the local Knudsen number, which indicates the state of the gas around the cylinder, can differ significantly from its free-stream value and it is essential that computational studies of subsonic gas flows in the slip and early transition regime are able to account for these strong non-equilibrium effects.