Using the dipole particles for simulation of 3d vortex flow of a viscous incompressible fluid

A fully lagrangian numerical method for simulation of 3D nonstationary flow of viscous and ideal incompressible fluid is developed in this work. This method is based on the representation of a vortex field as a set of dipole particles [1]. The introduced vector-function D describes density of dipole momentum. The equation for this function is in accordance with Navier-Stokes equations [2]. The vorticity is equal to curl of dipole momentum density. Thus vortex field is always solenoidal. The dipole particles are generated at a body surface and are moving interacting. The region where function D is essentially non-zero approximately coincides with the vortex region. Each dipole particle induces the velocity field which is equal to field of a point dipole at large distance from the particle. But near a particle the induced velocity field is another taking into account the particle volume and viscosity of the liquid. The method can be applied for simulation of an ideal and viscous flows

A Vortex Method for Bi-phasic Fluids Interacting with Rigid Bodies

We present an accurate Lagrangian method based on vortex particles, level-sets, and immersed boundary methods, for animating the interplay between two fluids and rigid solids. We show that a vortex method is a good choice for simulating bi-phase flow, such as liquid and gas, with a good level of realism. Vortex particles are localized at the interfaces between the two fluids and within the regions of high turbulence. We gain local precision and efficiency from the stable advection permitted by the vorticity formulation. Moreover, our numerical method straightforwardly solves the two-way coupling problem between the fluids and animated rigid solids. This new approach is validated through numerical comparisons with reference experiments from the computational fluid community. We also show that the visually appealing results obtained in the CG community can be reproduced with increased efficiency and an easier implementation

Wavelet transforms and their applications to MHD and plasma turbulence: a review

Wavelet analysis and compression tools are reviewed and different applications to study MHD and plasma turbulence are presented. We introduce the continuous and the orthogonal wavelet transform and detail several statistical diagnostics based on the wavelet coefficients. We then show how to extract coherent structures out of fully developed turbulent flows using wavelet-based denoising. Finally some multiscale numerical simulation schemes using wavelets are described. Several examples for analyzing, compressing and computing one, two and three dimensional turbulent MHD or plasma flows are presented.Comment: Journal of Plasma Physics, 201

Direct Numerical Simulation of decaying two-dimensional turbulence in a no-slip square box using Smoothed Particle Hydrodynamics

This paper explores the application of SPH to a Direct Numerical Simulation (DNS) of decaying turbulence in a two-dimensional no-slip wall-bounded domain. In this bounded domain, the inverse energy cascade, and a net torque exerted by the boundary, result in a spontaneous spin up of the fluid, leading to a typical end state of a large monopole vortex that fills the domain. The SPH simulations were compared against published results using a high accuracy pseudo-spectral code. Ensemble averages of the kinetic energy, enstrophy and average vortex wavenumber compared well against the pseudo-spectral results, as did the evolution of the total angular momentum of the fluid. However, while the pseudo-spectral results emphasised the importance of the no-slip boundaries as generators of long lived coherent vortices in the flow, no such generation was seen in the SPH results. Vorticity filaments produced at the boundary were always dissipated by the flow shortly after separating from the boundary layer. The kinetic energy spectrum of the SPH results was calculated using a SPH Fourier transform that operates directly on the disordered particles. The ensemble kinetic energy spectrum showed the expected k-3 scaling over most of the inertial range. However, the spectrum flattened at smaller length scales (initially less than 7.5 particle spacings and growing in size over time), indicating an excess of small-scale kinetic energy

How long do particles spend in vortical regions in turbulent flows?

We obtain the probability distribution functions (PDFs) of the time that a Lagrangian tracer or a heavy inertial particle spends in vortical or strain-dominated regions of a turbulent flow, by carrying out direct numerical simulation (DNS) of such particles advected by statistically steady, homogeneous and isotropic turbulence in the forced, three-dimensional, incompressible Navier-Stokes equation. We use the two invariants, $Q$ and $R$, of the velocity-gradient tensor to distinguish between vortical and strain-dominated regions of the flow and partition the $Q-R$ plane into four different regions depending on the topology of the flow; out of these four regions two correspond to vorticity-dominated regions of the flow and two correspond to strain-dominated ones. We obtain $Q$ and $R$ along the trajectories of tracers and heavy inertial particles and find out the time $\mathrm{t_{pers}}$ for which they remain in one of the four regions of the $Q-R$ plane. We find that the PDFs of $\mathrm{t_{pers}}$ display exponentially decaying tails for all four regions for tracers and heavy inertial particles. From these PDFs we extract characteristic times scales, which help us to quantify the time that such particles spend in vortical or strain-dominated regions of the flow

Extraction of coherent structures in a rotating turbulent flow experiment

The discrete wavelet packet transform (DWPT) and discrete wavelet transform (DWT) are used to extract and study the dynamics of coherent structures in a turbulent rotating fluid. Three-dimensional (3D) turbulence is generated by strong pumping through tubes at the bottom of a rotating tank (48.4 cm high, 39.4 cm diameter). This flow evolves toward two-dimensional (2D) turbulence with increasing height in the tank. Particle Image Velocimetry (PIV) measurements on the quasi-2D flow reveal many long-lived coherent vortices with a wide range of sizes. The vorticity fields exhibit vortex birth, merger, scattering, and destruction. We separate the flow into a low-entropy ``coherent'' and a high-entropy ``incoherent'' component by thresholding the coefficients of the DWPT and DWT of the vorticity fields. Similar thresholdings using the Fourier transform and JPEG compression together with the Okubo-Weiss criterion are also tested for comparison. We find that the DWPT and DWT yield similar results and are much more efficient at representing the total flow than a Fourier-based method. Only about 3% of the large-amplitude coefficients of the DWPT and DWT are necessary to represent the coherent component and preserve the vorticity probability density function, transport properties, and spatial and temporal correlations. The remaining small amplitude coefficients represent the incoherent component, which has near Gaussian vorticity PDF, contains no coherent structures, rapidly loses correlation in time, and does not contribute significantly to the transport properties of the flow. This suggests that one can describe and simulate such turbulent flow using a relatively small number of wavelet or wavelet packet modes.Comment: experimental work aprox 17 pages, 11 figures, accepted to appear in PRE, last few figures appear at the end. clarifications, added references, fixed typo

Small scale aspects of flows in proximity of the turbulent/non-turbulent interface

The work reported below is a first of its kind study of the properties of turbulent flow without strong mean shear in a Newtonian fluid in proximity of the turbulent/non-turbulent interface, with emphasis on the small scale aspects. The main tools used are a three-dimensional particle tracking system (3D-PTV) allowing to measure and follow in a Lagrangian manner the field of velocity derivatives and direct numerical simulations (DNS). The comparison of flow properties in the turbulent (A), intermediate (B) and non-turbulent (C) regions in the proximity of the interface allows for direct observation of the key physical processes underlying the entrainment phenomenon. The differences between small scale strain and enstrophy are striking and point to the definite scenario of turbulent entrainment via the viscous forces originating in strain.Comment: 4 pages, 4 figures, Phys. Fluid