No Arabic abstract
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.
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.
A homogenization approach is proposed for the treatment of porous wall boundary conditions in the computation of compressible viscous flows. Like any other homogenization approach, it eliminates the need for pore-resolved fluid meshes and therefore enables practical flow simulations in computational fluid domains with porous wall boundaries. Unlike alternative approaches however, it does not require prescribing a mass flow rate and does not introduce in the computational model a heuristic discharge coefficient. Instead, it models the inviscid flux through a porous wall surrounded by the flow as a weighted average of the inviscid flux at an impermeable surface and that through pores. It also introduces a body force term in the governing equations to account for friction loss along the pore boundaries. The source term depends on the thickness of the porous wall and the concept of an equivalent single pore. The feasibility of the latter concept is demonstrated using low-speed permeability test data for the fabric of the Mars Science Laboratory parachute canopy. The overall homogenization approach is illustrated with a series of supersonic flow computations through the same fabric and verified using supersonic, pore-resolved numerical simulations.
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.
A flow vessel with an elastic wall can deform significantly due to viscous fluid flow within it, even at vanishing Reynolds number (no fluid inertia). Deformation leads to an enhancement of throughput due to the change in cross-sectional area. The latter gives rise to a non-constant pressure gradient in the flow-wise direction and, hence, to a nonlinear flow rate--pressure drop relation (unlike the Hagen--Poiseuille law for a rigid tube). Many biofluids are non-Newtonian, and are well approximated by generalized Newtonian (say, power-law) rheological models. Consequently, we analyze the problem of steady low Reynolds number flow of a generalized Newtonian fluid through a slender elastic tube by coupling fluid lubrication theory to a structural problem posed in terms of Donnell shell theory. A perturbative approach (in the slenderness parameter) yields analytical solutions for both the flow and the deformation. Using matched asymptotics, we obtain a uniformly valid solution for the tubes radial displacement, which features both a boundary layer and a corner layer caused by localized bending near the clamped ends. In doing so, we obtain a ``generalized Hagen--Poiseuille law for soft microtubes. We benchmark the mathematical predictions against three-dimensional two-way coupled direct numerical simulations (DNS) of flow and deformation performed using the commercial computational engineering platform by ANSYS. The simulations show good agreement and establish the range of validity of the theory. Finally, we discuss the implications of the theory on the problem of the flow-induced deformation of a blood vessel, which is featured in some textbooks.
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.