No Arabic abstract
The influence of turbulent effects on a fluid flow through a (pseudo) porous media is studied by numerically solving the set of Reynolds-averaged Navier-Stokes equations with the $kappa$-$epsilon$ model for turbulence. The spatial domains are two-dimensional rectangular grids with different {it porosities} obtained by the random placing of rigid obstacles. The objective of the simulations is to access the behavior of the generalized friction factor with varying Reynolds number. A good agreement with the Forchheimers equation is observed. The flow distribution at both low and high Reynolds conditions is also analyzed.
Huge data advent in high-performance computing (HPC) applications such as fluid flow simulations usually hinders the interactive processing and exploration of simulation results. Such an interactive data exploration not only allows scientiest to play with their data but also to visualise huge (distributed) data sets in both an efficient and easy way. Therefore, we propose an HPC data exploration service based on a sliding window concept, that enables researches to access remote data (available on a supercomputer or cluster) during simulation runtime without exceeding any bandwidth limitations between the HPC back-end and the user front-end.
We investigate the elastoviscoplastic flow through porous media by numerical simulations. We solve the Navier-Stokes equations combined with the elastoviscoplastic model proposed by Saramito for the stress tensor evolution. In this model, the material behaves as a viscoelastic solid when unyielded, and as a viscoelastic Oldroyd-B fluid for stresses higher than the yield stress. The porous media is made of a symmetric array of cylinders, and we solve the flow in one periodic cell. We find that the solution is time-dependent even at low Reynolds numbers as we observe oscillations in time of the unyielded region especially at high Bingham numbers. The volume of the unyielded region slightly decreases with the Reynolds number and strongly increases with the Bingham number; up to 70% of the total volume is unyielded for the highest Bingham numbers considered here. The flow is mainly shear dominated in the yielded region, while shear and elongational flow are equally distributed in the unyielded region. We compute the relation between the pressure drop and the flow rate in the porous medium and present an empirical closure as function of the Bingham and Reynolds numbers. The apparent permeability, normalized with the case of Newtonian fluids, is shown to be greater than 1 at low Bingham numbers, corresponding to lower pressure drops due to the flow elasticity, and smaller than 1 for high Bingham numbers, indicating larger dissipation in the flow owing to the presence of the yielded regions. Finally we investigate the effect of the Weissenberg number on the distribution of the unyielded regions and on the pressure gradient.
Direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a face-centered cubic porous unit cell on the transport, clustering, and deposition of fine particles at different Stokes numbers ($St = 0.01, 0.1, 0.5, 1, 2$) and at a pore Reynolds number of 500. Particles are advanced using one-way coupling and collision of particles with pore walls is modeled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyze the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with Stokes number. However, unlike homogeneous, isotropic turbulence, formation of large number of very small volumes was observed at all Stokes numbers and is attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density showed peak of clustering shifts towards larger scales with increase in Stokes number. Scale-dependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers.
We study the gravity-driven flow of two fluid phases in a one-dimensional homogeneous porous column when history dependence of the pressure difference between the phases (capillary pressure) is taken into account. In the hyperbolic limit, solutions of such systems satisfy the Buckley-Leverett equation with a non-monotone flux function. However, solutions for the hysteretic case do not converge to the classical solutions in the hyperbolic limit in a wide range of situations. In particular, with Riemann data as initial condition, stationary shocks become possible in addition to classical components such as shocks, rarefaction waves, and constant states. We derive an admissibility criterion for the stationary shocks and outline all admissible shocks. Depending on the capillary pressure functions, flux function, and the Riemann data, two cases are identified a priori for which the solution consists of a stationary shock. In the first case, the shock remains at the point where the initial condition is discontinuous. In the second case, the solution is frozen in time in at least one semi-infinite half. The predictions are verified using numerical results.
Transport of liquid mixtures through porous membranes is central to processes such as desalination, chemical separations and energy harvesting, with ultrathin membranes made from novel 2D nanomaterials showing exceptional promise. Here we derive, for the first time, general equations for the solution and solute fluxes through a circular pore in an ultrathin planar membrane induced by a solute concentration gradient. We show that the equations accurately capture the fluid fluxes measured in finite-element numerical simulations for weak solute-membrane interactions. We also derive scaling laws for these fluxes as a function of the pore size and the strength and range of solute-membrane interactions. These scaling relationships differ markedly from those for concentration-gradient-driven flow through a long cylindrical pore or for flow induced by a pressure gradient or electric field through a pore in an ultrathin membrane. These results have broad implications for transport of liquid mixtures through membranes with a thickness on the order of the characteristic pore size.