No Arabic abstract
To study the dynamics of particles in turbulence when their sizes are comparable to the smallest eddies in the flow, the Kolmogorov length scale, efficient and accurate numerical models for the particle-fluid interaction are still missing. Therefore, we here extend the treatment of the particle feedback on the fluid based on the volume-averaged fluid equations (VA simulation) in the previous study of the present authors, by estimating the fluid force correlated with the disturbed flow. We validate the model against interface-resolved simulations using the immersed-boundary method. Simulations of single particles show that the history effect is well captured by the present estimation method based on the disturbed flow. Similarly, the simulation of the flow around a rotating particle demonstrates that the lift force is also well captured by the proposed method. We also consider the interaction between non-negligible size particles and an array of Taylor-Green vortices. For density ratios $rho_d/rho_cgeq$ 10, the results show that the particle motion captured by the VA approach is closer to that of the fully-resolved simulations than that obtained with a traditional two-way coupling simulation. The flow disturbance is also well represented by the VA simulation. In particular, it is found that history effects enhance the curvature of the trajectory in vortices and this enhancement increases with the particle size. Furthermore, the flow field generated by a neighboring particle at distances of around ten particle diameters significantly influences particle trajectories. The computational cost of the VA simulation proposed here is considerably lower than that of the interface-resolved simulation.
A dynamic procedure for the Lagrangian Averaged Navier-Stokes-$alpha$ (LANS-$alpha$) equations is developed where the variation in the parameter $alpha$ in the direction of anisotropy is determined in a self-consistent way from data contained in the simulation itself. The dynamic model is initially tested in forced and decaying isotropic turbulent flows where $alpha$ is constant in space but it is allowed to vary in time. It is observed that by using the dynamic LANS-$alpha$ procedure a more accurate simulation of the isotropic homogeneous turbulence is achieved. The energy spectra and the total kinetic energy decay are captured more accurately as compared with the LANS-$alpha$ simulations using a fixed $alpha$. In order to evaluate the applicability of the dynamic LANS-$alpha$ model in anisotropic turbulence, a priori test of a turbulent channel flow is performed. It is found that the parameter $alpha$ changes in the wall normal direction. Near a solid wall, the length scale $alpha$ is seen to depend on the distance from the wall with a vanishing value at the wall. On the other hand, away from the wall, where the turbulence is more isotropic, $alpha$ approaches an almost constant value. Furthermore, the behavior of the subgrid scale stresses in the near wall region is captured accurately by the dynamic LANS-$alpha$ model. The dynamic LANS-$alpha$ model has the potential to extend the applicability of the LANS-$alpha$ equations to more complicated anisotropic flows.
This paper has been withdrawn by the authors for adding some results.
An investigation of optimal feedback controllers performance and robustness is carried out for vortex shedding behind a 2D cylinder at low Reynolds numbers. To facilitate controller design, we present an efficient modelling approach in which we utilise the resolvent operator to recast the linearised Navier-Stokes equations into an input-output form from which frequency responses can be computed. The difficulty of applying modern control design techniques to complex, high-dimensional flow systems is thus overcome by using low-order models identified from these frequency responses. The low-order models are used to design optimal control laws using $mathcal{H}_{infty}$ loop shaping. Two distinct control arrangements are considered, both of which employ a single-input and a single-output. In the first control arrangement, a velocity sensor located in the wake drives a pair of body forces near the cylinder. Complete suppression of shedding is observed up to a Reynolds number of $Re=110$. Due to the convective nature of vortex shedding and the corresponding time delays, we observe a fundamental trade-off: the sensor should be close enough to the cylinder to avoid any excessive time lag, but it should be kept sufficiently far from the cylinder to measure any unstable modes developing downstream. It is found that these two conflicting requirements become more difficult to satisfy for larger Reynolds numbers. In the second control arrangement, we consider a practical setup with a body-mounted force sensor and an actuator that oscillates the cylinder according to the lift measurement. It is shown that the system is stabilised only up to $Re=100$, and we demonstrate why the performance of the resulting feedback controllers deteriorates much more rapidly with increasing Reynolds number. The challenges of designing robust controllers for each control setup are also analysed and discussed.
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.
Three dimensional roll-type double-diffusive convection in a horizontally infinite layer of an uncompressible liquid is considered in the neighborhood of Hopf bifurcation points. A system of amplitude equations for the variations of convective rolls amplitude is derived by multiple-scaled method. An attention is paid to an interaction of convection and horizontal vortex. Different cases of the derived equations are discussed.