We propose a theoretical framework to calculate capillary stresses in complex mesoporous materials, such as moist sand, nanoporous hydrates, and drying colloidal films. Molecular simulations are mapped onto a phase-field model of the liquid-vapor mixture, whose inhomogeneous stress tensor is integrated over Voronoi polyhedra in order to calculate equal and opposite forces between each pair of neighboring grains. The method is illustrated by simulations of moisture-induced forces in small clusters and random packings of spherical grains using lattice-gas Density Functional Theory. For a nano-granular model of cement hydrates, this approach reproduces the hysteretic water sorption/desorption isotherms and predicts drying shrinkage strain isotherm in good agreement with experiments. We show that capillary stress is an effective mechanism for internal stress relaxation in colloidal random packings, which contributes to the extraordinary durability of cement paste.

This letter presents a scaling theory of the coalescence of two viscous spherical droplets. An initial value problem was formulated and analytically solved for the evolution of the radius of a liquid neck formed upon droplet coalescence. Two asymptotic solutions of the initial value problem reproduce the well-known scaling relations in the viscous and inertial regimes. The viscous-to-inertial crossover experimentally observed by Paulsen et al. [Phys. Rev. Lett. 106, 114501 (2011)] manifests in the theory, and their fitting relation, which shows collapse of data of different viscosities onto a single curve, is an approximation to the general solution of the initial value problem.

Recent studies within the diffusive-thermal (constant-density) approximation have shown that, for premixed flames freely propagating in narrow adiabatic channels, the instabilities induced by differential diffusion may result in non-symmetric solutions and/or oscillating and rotating propagation modes. This has been shown in the context of lean mixtures, for which a single species transport equation with a single Lewis number (corresponding to the ratio of thermal to molecular diffusivity of the deficient reactant) can be used to describe the flame propagation problem. In the present work we extend the analysis to mixtures of any equivalence ratio. To this end, we consider a two-reactant model, where the different diffusivities of the two reactants introduce two different Lewis numbers. Steady-state computations and linear stability analysis are carried out for mixtures with large disparity between the Lewis number of the fuel (LeF) and the oxidizer (LeO), such as hydrogen-oxygen systems. It is shown that both differential diffusion and preferential diffusion have influence on the stability of the symmetric flame shape. For sufficiently lean and rich mixtures, the flame behaves as dictated by the Lewis number of the deficient reactant, i.e., the flame destabilizes toward non-symmetric solutions for large mass flow rates when the mentioned Lewis number is less than one. In near-stoichiometric mixtures the stability of the symmetric flame depends on a weighted average value of LeF and LeO. In particular, the symmetric solution is stable for large mass flow rates because of the difficulties found by the less diffusive reactant to reach the reactive zone of the flame.

Understanding rare events in turbulence provides a basis for science of extreme weather, for which the atmosphere is modeled by the Navier-Stokes equations (NSE). In solutions of the NSE for isotropic fluids, various quantities, such as fluid velocities, roughly follow Gaussian distributions, where rare events are promoted only in asymptotically small-scale quantities associated with the dissipation-dominating length scale. Using numerical simulations, here, we reveal another universal promotion mechanism at much larger scales, if three-dimensional fluids accompany strong two-dimensional anisotropies as in the atmosphere. The dimensional crossover between two and three dimensions generates prominent fat-tailed non-Gaussian distributions with intermittency accompanied by colossal large-scale chain-like structures with densely populated self-organized vortices. The promotion is caused by a sudden increase of the available phase space at the crossover. Since the discovered intermittency bears much larger energies than those in the conventional intermittency in small spatial scales, it would govern extreme events and chaotic unpredictability in the synoptic weather system on the earth.

This paper presents experimental and numerical analysis of grid generated turbulence with and without the effects of applied mean strain. We conduct a series of experiments on decaying grid generated turbulence and grid turbulence with mean strain. Experimental data of turbulence statistics including Reynolds stress anisotropies is collected, analyzed and then compared to the predictions of Reynolds Stress Models to assess their accuracy. The experimental data is used to evaluate the variability in the coefficients of the rate of dissipation model and the pressure strain correlation models used in Reynolds Stress Modeling. For both models we recommend optimal values of coefficients that should be used for experimental studies of grid generated turbulence.

Reynolds-averaged Navier-Stokes (RANS) simulations with turbulence models continue to play important roles in industrial flow simulations as high-fidelity simulations are prohibitively expensive for such flows. Commonly used linear eddy viscosity models are intrinsically unable to handle flows with non-equilibrium turbulence (e.g., flows with massive separation). Reynolds stress models, on the other hand, are plagued by their lack of robustness and stability. Recent studies found that even substituting Reynolds stresses from DNS databases (with errors below 0.5%) into RANS equations leads to grossly inaccurate velocities. Such an observation is not only disturbing for the recently emerging data-driven Reynolds stress models but also relevant for traditional, equation-based models. This observation cannot be explained by the global matrix condition number of the discretized RANS equations. In this work, we propose a metric based on local condition numbers for a priori evaluation of the stability of Reynolds stress models. Numerical tests on turbulent channel flows at various Reynolds numbers suggest that the proposed metric can adequately explain observations in previous studies, i.e., decreased model stability with increasing Reynolds number, and better stability of the implicit treatment of Reynolds stress compared to the explicit treatment.

We investigate the dynamics of hydrodynamically interacting motile and non-motile stress-generating swimmers or particles as they invade a surrounding viscous fluid. Colonies of aligned pusher particles are shown to elongate in the direction of particle orientation and undergo a cascade of transverse concentration instabilities. Colonies of aligned puller particles instead are found to elongate in the direction opposite the particle orientation and exhibit dramatic splay as the group moves into the bulk. A linear stability analysis of concentrated line distributions of particles is performed and growth rates are found, using an active slender-body approximation, to match the results of numerical simulations. Thin concentrated bands of aligned pusher particles are always unstable, while bands of aligned puller particles can either be stable (immotile particles) or unstable (motile particles) with a growth rate which is non-monotonic in the force dipole strength. We also prove a surprising "no-flow theorem": a distribution initially isotropic in orientation loses isotropy immediately but in such a way that results in no fluid flow anywhere at any time.

Many astrophysical systems encompass very large dynamical ranges in space and time, which are not accessible by direct numerical simulations. Thus, idealized subvolumes are often used to study small-scale effects including the dynamics of turbulence. These turbulent boxes require an artificial driving in order to mimic energy injection from large-scale processes. In this Letter, we show and quantify how the autocorrelation time of the driving and its normalization systematically change properties of an isothermal compressible magnetohydrodynamic flow in the sub- and supersonic regime and affect astrophysical observations such as Faraday rotation. For example, we find that $\delta$-in-time forcing with a constant energy injection leads to a steeper slope in kinetic energy spectrum and less efficient small-scale dynamo action. In general, we show that shorter autocorrelation times require more power in the acceleration field, which results in more power in compressive modes that weaken the anticorrelation between density and magnetic field strength. Thus, derived observables, such as the line-of-sight magnetic field from rotation measures, are systematically biased by the driving mechanism. We argue that $\delta$-in-time forcing is unrealistic and numerically unresolved, and conclude that special care needs to be taken in interpreting observational results based on the use of idealized simulations.