Search

In this work, we aim to deepen the understanding of inertial clustering and the role of sling events in high-Reynolds number ($Re$) particle-laden turbulence. To this end, we perform one-way coupled particle tracking in flow fields obtained from direct numerical simulations (DNS) of forced homogeneous isotropic turbulence. Additionally, we examine the impact of filtering utilized in large eddy simulations (LES) by applying a sharp spectral filter to the DNS fields. Our analysis reveals that while instantaneous clustering through the centrifuge mechanism explains clustering at early times, the path history effect--the sampling of fluid flow along particle trajectories--becomes important later on. The filtered fields expose small-scale fractal clustering that cannot be predicted by the instantaneous flow field. We show that there exists a filter-effective Stokes number that governs the degree of fractal clustering and preferential sampling, revealing scale-similarity in the spatial distributions and fractal dimensions. Sling events are prevalent throughout our simulations and impose prominent patterns on the particle fields. In pursuit of investigating the sling dynamics, we compute the relative velocity, ensemble-averaged over proximal neighboring particles, to identify particles undergoing caustics. As postulated in recent theories, we find that in fully resolved, high-$Re$ turbulence, sling events occur in thin sheets of high strain, situated between turbulent vortices. This behavior is driven by rare, extreme events of compressive straining, manifested by fluctuations of the flow velocity gradients that propagate back and forth the positive branch of the Vieillefosse line.

This work is concerned with the intricate interplay between node or pore pressures and connection or throat conductivities in flow or pore networks. A setting similar to pore networks is given by fracture networks. Recently, a non-local generalization of Darcy's law for flow and transport in porous media was presented in the context of unbounded or periodic pore networks. In this work, we first outline a robust method for the extraction of the hydraulic conductivity distribution, which is at the heart of the non-local Darcy formulation. Second, a theory for mean pressure and flow in bounded networks is outlined. Predictions of that theory are validated against numerical network results and it is demonstrated that the theory works well for networks with high connectivity involving pores with high coordination numbers. For other networks, improvements to the outlined theory are proposed and their accuracy is assessed.

Recently, distribution element trees (DETs) were introduced as an accurate and computationally efficient method for density estimation. In this work, we demonstrate that the DET formulation promotes an easy and inexpensive way to generate random samples similar to a smooth bootstrap. These samples can be generated unconditionally, but also, without further complications, conditionally utilizing available information about certain probability-space components., Comment: published online in the Journal of Computational and Graphical Statistics

In the context of stochastic two-phase flow in porous media, we introduce a novel and efficient method to estimate the probability distribution of the wetting saturation field under uncertain rock properties in highly heterogeneous porous systems, where streamline patterns are dominated by permeability heterogeneity, and for slow displacement processes (viscosity ratio close to unity). Our method, referred to as the frozen streamline distribution method (FROST), is based on a physical understanding of the stochastic problem. Indeed, we identify key random fields that guide the wetting saturation variability, namely fluid particle times of flight and injection times. By comparing saturation statistics against full-physics Monte Carlo simulations, we illustrate how this simple, yet accurate FROST method performs under the preliminary approximation of frozen streamlines. Further, we inspect the performance of an accelerated FROST variant that relies on a simplification about injection time statistics. Finally, we introduce how quantiles of saturation can be efficiently computed within the FROST framework, hence leading to robust uncertainty assessment., Comment: 34 pages

The estimation of probability densities based on available data is a central task in many statistical applications. Especially in the case of large ensembles with many samples or high-dimensional sample spaces, computationally efficient methods are needed. We propose a new method that is based on a decomposition of the unknown distribution in terms of so-called distribution elements (DEs). These elements enable an adaptive and hierarchical discretization of the sample space with small or large elements in regions with smoothly or highly variable densities, respectively. The novel refinement strategy that we propose is based on statistical goodness-of-fit and pair-wise (as an approximation to mutual) independence tests that evaluate the local approximation of the distribution in terms of DEs. The capabilities of our new method are inspected based on several examples of different dimensionality and successfully compared with other state-of-the-art density estimators.