No Arabic abstract
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.
Discrete Boltzmann model (DBM) is a type of coarse-grained mesoscale kinetic model derived from the Boltzmann equation. Physically, it is roughly equivalent to a hydrodynamic model supplemented by a coarse-grained model for the relevant thermodynamic non-equilibrium (TNE) behaviours. The Navier-Stokes (NS) model is a traditional macroscopic hydrodynamic model based on continuity hypothesis and conservation laws. In this study, the two models are compared from two aspects, physical capability and computational cost, by simulating two kinds of flow problems including the thermal Couette flow and a Mach 3 step problem. In the cases where the TNE effects are weak, both the two models give accurate results for the hydrodynamic behaviour. Besides, DBM can provide more detailed non-equilibrium information, while the NS is more efficient if concern only the density, momentum, energy and their derived quantities. It is concluded that, if the TNE effects are strong or are to be investigated, the NS is insufficient while DBM is a good choice. While in the cases where the TNE effects are weak and only the macro flow fields are to be studied, the NS is more preferable.
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 study a correspondence between the multifractal model of turbulence and the Navier-Stokes equations in $d$ spatial dimensions by comparing their respective dissipation length scales. In Kolmogorovs 1941 theory the key parameter $h$, which is an exponent in the Navier-Stokes invariance scaling, is fixed at $h=1/3$ but is allowed a spectrum of values in multifractal theory. Taking into account all derivatives of the Navier-Stokes equations, it is found that for this correspondence to hold the multifractal spectrum $C(h)$ must be bounded from below such that $C(h) geq 1-3h$, which is consistent with the four-fifths law. Moreover, $h$ must also be bounded from below such that $h geq (1-d)/3$. When $d=3$ the allowed range of $h$ is given by $h geq -2/3$ thereby bounding $h$ away from $h=-1$. The implications of this are discussed.
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.