No Arabic abstract
In a recent paper, Liu, Zhu and Wu (2015, {it J. Fluid Mech.} {bf 784}: 304) present a force theory for a body in a two-dimensional, viscous, compressible and steady flow. In this companion paper we do the same for three-dimensional flow. Using the fundamental solution of the linearized Navier-Stokes equations, we improve the force formula for incompressible flow originally derived by Goldstein in 1931 and summarized by Milne-Thomson in 1968, both being far from complete, to its perfect final form, which is further proved to be universally true from subsonic to supersonic flows. We call this result the textit{unified force theorem}, which states that the forces are always determined by the vector circulation $pGamma_phi$ of longitudinal velocity and the scalar inflow $Q_psi$ of transverse velocity. Since this theorem is not directly observable either experimentally or computationally, a testable version is also derived, which, however, holds only in the linear far field. We name this version the textit{testable unified force formula}. After that, a general principle to increase the lift-drag ratio is proposed.
In a recent paper, Liu et al. [``Lift and drag in three-dimensional steady viscous and compressible flow, Phys. Fluids 29, 116105 (2017)] obtained a universal theory for the aerodynamic force on a body in three-dimensional steady flow, effective from incompressible all the way to supersonic regimes. In this theory, the total aerodynamic force can be determined solely with the vorticity distribution on a single wake plane locating in the steady linear far field. Despite the vital importance of this result, its validity and performance in practice has not been investigated yet. In this paper, we performed Reynolds-averaged Navier-Stokes simulations of subsonic, transonic, and supersonic flows over a three-dimensional wing. The aerodynamic forces obtained from the universal force theory are compared with that from the standard wall-stress integrals. The agreement between these two formulas confirms for the first time the validity of the theory in three-dimensional steady viscous and compressible flow. The good performance of the universal formula is mainly due to the fact that the turbulent viscosity in the wake is much larger than the molecular viscosity therein, which can reduce significantly the distance of the steady linear far field from the body. To further confirm the correctness of the theory, comparisons are made for the flow structures on the wake plane obtained from the analytical results and numerical simulations. The underlying physics relevant to the universality of the theory is explained by identifying different sources of vorticity in the wake.
We present a series of three-dimensional discrete Boltzmann (DB) models for compressible flows in and out of equilibrium. The key formulating technique is the construction of discrete equilibrium distribution function through inversely solving the kinetic moment relations that it satisfies. The crucial physical requirement is that all the used kinetic moment relations must be consistent with the non-equilibrium statistical mechanics. The necessity of such a kinetic model is that, with increasing the complexity of flows, the dynamical characterization of non-equilibrium state and the understanding of the constitutive relations need higher order kinetic moments and their evolution. The DB models at the Euler and Navier-Stokes levels proposed by this scheme are validated by several well-known benchmarks, ranging from one-dimension to three-dimension. Particularly, when the local Mach number, temperature ratio, and pressure ratio are as large as $10^2$, $10^4$, and $10^5$, respectively, the simulation results are still in excellent agreement with the Riemann solutions. How to model deeper thermodynamic non-equilibrium flows by DB is indicated. Via the DB method, it convenient to simulate nonequilibrium flows without knowing exact form of the hydrodynamic equations.
Conflict between formation of a cyclonic vortex and isotropization in forced homogeneous rotating turbulence is numerically investigated. It is well known that a large rotation rate of the system induces columnar vortices to result in quasi-two-dimensional (Q2D) flow, while a small rotation rate allows turbulence to be three-dimensional (3D). It is found that the transition from the Q2D turbulent flow to the 3D turbulent flow and the reverse transition occur at different values of the rotation rates. At the intermediate rotation rates, bistability of these two statistically steady states is observed. Such hysteretic behavior is also observed for the variation of the amplitude of an external force.
In this chapter, we analyze the steady-state microscale fluid--structure interaction (FSI) between a generalized Newtonian fluid and a hyperelastic tube. Physiological flows, especially in hemodynamics, serve as primary examples of such FSI phenomena. The small scale of the physical system renders the flow field, under the power-law rheological model, amenable to a closed-form solution using the lubrication approximation. On the other hand, negligible shear stresses on the walls of a long vessel allow the structure to be treated as a pressure vessel. The constitutive equation for the microtube is prescribed via the strain energy functional for an incompressible, isotropic Mooney--Rivlin material. We employ both the thin- and thick-walled formulations of the pressure vessel theory, and derive the static relation between the pressure load and the deformation of the structure. We harness the latter to determine the flow rate--pressure drop relationship for non-Newtonian flow in thin- and thick-walled soft hyperelastic microtubes. Through illustrative examples, we discuss how a hyperelastic tube supports the same pressure load as a linearly elastic tube with smaller deformation, thus requiring a higher pressure drop across itself to maintain a fixed flow rate.
We study analytically the joint dispersion of Gaussian patches of salt and colloids in linear flows, and how salt gradients accelerate or delay colloid spreading by diffusiophoretic effects. Because these flows have constant gradients in space, the problem can be solved almost entirely for any set of parameters, leading to predictions of how the mixing time and the Batchelor scale are modified by diffusiophoresis. We observe that the evolution of global concentrations, defined as the inverse of the patches areas, are very similar to those obtained experimentally in chaotic advection. They are quantitatively explained by examining the area dilatation factor, in which diffusive and diffusiophoretic effects are shown to be additive and appear as the divergence of a diffusive contribution or of a drift velocity. An analysis based on compressibility is developed in the salt-attracting case, for which colloids are first compressed before dispersion, to predict the maximal colloid concentration as a function of the parameters. This maximum is found not to depend on the flow stretching rate nor on its topology (strain or shear flow), but only on the characteristics of salt and colloids (diffusion coefficients and diffusiophoretic constant) and the initial size of the patches.