No Arabic abstract
We explore the utility of the recently proposed alpha equations in providing a subgrid model for fluid turbulence. Our principal results are comparisons of direct numerical simulations of fluid turbulence using several values of the parameter alpha, including the limiting case where the Navier-Stokes equations are recovered. Our studies show that the large scale features, including statistics and structures, are preserved by the alpha models, even at coarser resolutions where the fine scales are not fully resolved. We also describe the differences that appear in simulations. We provide a summary of the principal features of the alpha equations, and offer some explanation of the effectiveness of these equations used as a subgrid model for three-dimensional fluid turbulence.
We determine how the differences in the treatment of the subfilter-scale physics affect the properties of the flow for three closely related regularizations of Navier-Stokes. The consequences on the applicability of the regularizations as SGS models are also shown by examining their effects on superfilter-scale properties. Numerical solutions of the Clark-alpha model are compared to two previously employed regularizations, LANS-alpha and Leray-alpha (at Re ~ 3300, Taylor Re ~ 790) and to a DNS. We derive the Karman-Howarth equation for both the Clark-alpha and Leray-alpha models. We confirm one of two possible scalings resulting from this equation for Clark as well as its associated k^(-1) energy spectrum. At sub-filter scales, Clark-alpha possesses similar total dissipation and characteristic time to reach a statistical turbulent steady-state as Navier-Stokes, but exhibits greater intermittency. As a SGS model, Clark reproduces the energy spectrum and intermittency properties of the DNS. For the Leray model, increasing the filter width decreases the nonlinearity and the effective Re is substantially decreased. Even for the smallest value of alpha studied, Leray-alpha was inadequate as a SGS model. The LANS energy spectrum k^1, consistent with its so-called rigid bodies, precludes a reproduction of the large-scale energy spectrum of the DNS at high Re while achieving a large reduction in resolution. However, that this same feature reduces its intermittency compared to Clark-alpha (which shares a similar Karman-Howarth equation). Clark is found to be the best approximation for reproducing the total dissipation rate and the energy spectrum at scales larger than alpha, whereas high-order intermittency properties for larger values of alpha are best reproduced by LANS-alpha.
This paper has been withdrawn by the authors for adding some results.
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.
We introduce a model of interacting singularities of Navier-Stokes, named pin,cons. They follow a Hamiltonian dynamics, obtained by the condition that the velocity field around these singularities obeys locally Navier-Stokes equations. This model can be seen of a generalization of the vorton model of Novikov, that was derived for the Euler equations. When immersed in a regular field, the pin,cons are further transported and sheared by the regular field, while applying a stress onto the regular field, that becomes dominant at a scale that is smaller than the Kolmogorov length. We apply this model to compute the motion of a dipole of pin,cons. When the initial relative orientation of the dipole is inside the interval (0, pi/2), a dipole made of pin,con of same intensity exhibits a transient collapse stage, following a scaling with dipole radius tending to 0 like (tc - t) power 0.63. For long time, the dynamics of the dipole is however repulsive, with both components running away from each other to infinity.
We consider a computational model for complex-fluid-solid interaction based on a diffuse-interface model for the complex fluid and a hyperelastic-material model for the solid. The diffuse-interface complex-fluid model is described by the incompressible Navier-Stokes-Cahn-Hilliard equations with preferential-wetting boundary conditions at the fluid-solid interface. The corresponding fluid traction on the interface includes a capillary-stress contribution, and the dynamic interface condition comprises the traction exerted by the non-uniform fluid-solid surface tension. We present a weak formulation of the aggregated complex-fluid-solid-interaction problem, based on an Arbitrary-Lagrangian-Eulerian formulation of the Navier-Stokes-Cahn-Hilliard equations and a proper reformulation of the complex-fluid traction and the fluid-solid surface tension. To validate the presented complex-fluid-solid-interaction model, we present numerical results and conduct a comparison to experimental data for a droplet on a soft substrate.