No Arabic abstract
We conduct depth-resolved three-dimensional Direct Numerical Simulations (DNS) of bi-disperse turbidity currents interacting with complex bottom topography in the form of a Gaussian bump. Several flow characteristics such as suspended particle mass, instantaneous wall shear stress, transient deposit height are shown via videos. Furthermore, we investigate the influence of the obstacle on the vortical structure and sedimentation of particles by comparing the results against the same setup and but with a flat bottom surface. We observe that the obstacle influences the deposition of coarse particles mainly in the vicinity of the obstacle due to lateral deflection, whereas for the sedimentation of fine particles the effects of topographical features are felt further downstream. The results shown in this fluid dynamics video help us develop a fundamental understanding of the dynamics of turbidity currents interacting with complex seafloor topography.
The propagation of surface water waves interacting with a current and an uneven bottom is studied. Such a situation is typical for ocean waves where the winds generate currents in the top layer of the ocean. The role of the bottom topography is taken into account since it also influences the local wave and current patterns. Specific scaling of the variables is selected which leads to approximations of Boussinesq and KdV types. The arising KdV equation with variable coefficients, dependent on the bottom topography, is studied numerically when the initial condition is in the form of the one soliton solution for the initial depth. Emergence of new solitons is observed as a result of the wave interaction with the uneven bottom.
We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply-connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface. We prove that the resulting second-kind Fredholm integral equations are invertible. In the velocity potential formulation, invertibility is achieved after a physically motivated finite-rank correction. The integral equations for the two methods are closely related, one being the adjoint of the other after modifying it to evaluate the layer potentials on the opposite side of each interface. In addition to a background flow, both formulations allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. The proposed boundary integral methods are compatible with graph-based or angle-arclength parameterizations of the free surface. In the latter case, we show how to avoid curve reconstruction errors in interior Runge-Kutta stages due to incompatibility of the angle-arclength representation with spatial periodicity. The proposed methods are used to study gravity-capillary waves generated by flow over three elliptical obstacles with different choices of the circulation parameters. In each case, the free surface forms a structure resembling a Crapper wave that narrows and eventually self intersects in a splash singularity.
We present a comparison of different particles velocity and acceleration statistics in two paradigmatic turbulent swirling flows: the von Karman flow in a laboratory experiment, and the Taylor-Green flow in direct numerical simulations. Tracers, as well as inertial particles, are considered. Results indicate that, in spite of the differences in boundary conditions and forcing mechanisms, scaling properties and statistical quantities reveal similarities between both flows, pointing to new methods to calibrate and compare models for particles dynamics in numerical simulations, as well as to characterize the dynamics of particles in simulations and experiments.
A general, two-way coupled, point-particle formulation that accounts for the disturbance created by the dispersed particles in obtaining the undisturbed fluid flow field needed for accurate computation of the force closure models is presented. Specifically, equations for the disturbance field created by the presence of particles are first derived based on the inter-phase momentum coupling force in a finite-volume formulation. Solution to the disturbance field is obtained using two approaches: (i) direct computation of the disturbance velocity and pressure using the reaction force due to particles at computational control volumes, and (ii) a linearized, approximate computation of the disturbance velocity field, specifically applicable for low Reynolds number flows. In both approaches, the computed disturbance field is used to obtain the undisturbed fluid velocity necessary to model the aerodynamic forces on the particle. The two approaches are thoroughly evaluated for a single particle in an unbounded and wall-bounded flow on uniform, anisotropic, as well as unstructured grids to show accurate computation of the particle motion and inter-phase coupling. The approach is straightforward and can be applied to any numerical formulation for particle-laden flows including Euler-Lagrange as well as Euler-Euler formulations.
A hybrid parallel approach for fully resolved simulations of particle-laden flows in sediment transport is proposed. To overcome the challenges of load imbalance in the traditional domain decomposition method when encountering highly uneven distributions of particles in space, we develop a hybrid parallel approach adopting the domain decomposition method for the carrier phase and a mirror domain technique for the disperse phase. We modify the mirror domain technique originally developed for point particles to fully resolved particle simulations, which are more challenging since a finite-sized particle may be split into different subdomains; thus, more complex treatments of particle-fluid interactions are needed. By utilizing the mirror domain technique, in which each processor handles nearly the same number of particles regardless of the particle spatial distribution, excellent load balance is achieved. The present hybrid parallel approach also shows strong scalability and high parallel efficiency in a test of a fully resolved simulation case of sediment transport. Furthermore, a novel memory optimization method is proposed for spherical particles of equal size, which can substantially reduce the memory cost and enable the simulation of millions of fully resolved particles on a common highly parallel computing platform. Our code is validated by several benchmark cases, and the results show good agreement with experimental and computational data in the literature.