Aerodynamic drag of bodies in transonic flow. Theory and applications to computational aerodynamics

Formulas for all the components of the aerodynamic drag (total, friction, inductive, wave, pressure, and heat-transfer) are uniformly derived as applied to flows governed by the Navier-Stokes and Reynolds equations. For flows of this type the definition of the aerodynamic drag components is refined and the physical basis of the chosen method of breaking up the total drag into components is discussed. Ways of calculating the aerodynamic drag components using the methods of computational aerodynamics are considered. On the basis of the refined formulas the drag components are calculated for flows around airfoils and a high-aspect-ratio wing in transonic flow.     

Mechanism of drag reduction by spanwise oscillating Lorentz force in turbulent channel flow

Numerical simulations and experimental research are both carried out to investigate the controlled effect of spanwise oscillating Lorentz force on a turbulent channel flow. The variations of the streaks and the skin friction drag are obtained through the PIV system and the drag measurement system, respectively. The flow field in the near-wall region is shown through direct numerical simulations utilizing spectral method. The experimental results are consistent with the numerical simulation results qualitatively, and both the results indicate that the streaks are tilted into the spanwise direction and the drag reduction utilizing spanwise oscillating Lorentz forces can be realized. The numerical simulation results reveal more detail of the drag reduction mechanism which can be explained, since the spanwise vorticity generated from the interaction between the induced Stokes layer and intrinsic turbulent flow in the near-wall region can make the longitudinal vortices tilt and oscillate, and leads to turbulence suppression and drag reduction.     

Nonisothermal flow of a viscous incompressible fluid in a two-dimensional diffuser

The velocity and temperature distributions in a viscous incompressible fluid flow in a two-dimensional diffuser are analyzed. Fully developed flow is considered, i.e., the influence of the entrant section is disregarded. It is assumed that the diffuser walls are maintained at a temperature depending on the polar radius. The dynamic viscosity is considered to be an exponential function of the temperature. The problem is reduced to the solution of a system of ordinary differential equations, which is solved by the method of successive approximations. The convergence of the iterative scheme is proved.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 40–48, July–August, 1973.The author is indebted to L.A. Galin and N. N. Gvozdkov for assistance with the study.     

An analytic approximation of the drag coefficient for the viscous flow past a sphere

We give an analytic solution at the 10th order of approximation for the steady-state laminar viscous flows past a sphere in a uniform stream governed by the exact, fully non-linear Navier-Stokes equations. A new kind of analytic technique, namely the homotopy analysis method, is applied, by means of which Whitehead's paradox can be easily avoided and reasonably explained. Different from all previous perturbation approximations, our analytic approximations are valid in the whole field of flow, because we use the same approximations to express the flows near and far from the sphere. Our drag coefficient formula at the 10th order of approximation agrees better with experimental data in a region of Reynolds number R_d<30, which is considerably larger than that (R_d<5) of all previous theoretical ones.     

超声速钝体逆向喷流减阻的数值模拟研究

为研究逆向喷流技术对超声速钝体减阻的影响,采用标准k-ε湍流模型,通过求解二维Navier-Stokes方程对超声速球头体逆向冷喷流流场进行了数值模拟,并着重分析了喷口总压、喷口尺寸对流场模态和减阻效果的影响。计算结果显示:随着喷流总压的变化,流场可出现两种流动模态,即长射流穿透模态和短射流穿透模态;喷流能使球头体受到的阻力明显减小;存在最大减阻临界喷流总压值(在所研究参数范围内最大减阻可达51.1%);在其它喷流物理参数不变时,随着喷口尺寸的增大,同一流动模态下的减阻效果下降。本文的研究对超声速钝体减阻技术在工程上的应用具有一定的参考价值。     

Heat transfer associated with viscous two-dimensional channel flow in a transverse acoustic field

An analytic solution is proposed for the problem of the effect of a transverse resonant acoustic field on the heat transfer process in laminar two-dimensional channel flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 41–49, May–June, 1985.The author wishes to thank V. E. Nakoryakov for formulating the problem and discussing the results.     

The role of small scale coherent structures in turbulent drag reduction

With a new visualization technique, a strong suppression of small scale coherent structures (filaments) in turbulence is observed when a drag reducing polymer is added. Simultaneously, the rate of formation of large scale structures (eddies) out of these filaments decreases, which may be an important observation in the explanation of turbulent drag reduction.     

Quasi-steady deformation of a two-dimensional bubble placed within a potential viscous flow

The shape evolution of a two-dimensional bubble, bounded by a simple closed curve, which is initially placed within a potential viscous flow, is analysed. It is assumed that the influence of gravity and inertia forces is negligible, so the quasi-steady approximation can be applied. Reformulating the problem for Stokes equations with relevant boundary conditions at the free surface in terms of the bianalytic stress-stream function, and using the time-dependent conformal mappingz(,t) of a unit disk onto an unbounded flow domain sought, an infinite system of ordinary differential equations for the Laurent coefficients ofz(,t) is derived. A class of exact solutions is found for the case when the principal part of the complex velocity of the dominant flow at infinity is a polynomial, and the problem of formation of a pointed bubble is discussed.

---

Sommario E analizzata l'evoluzione di una bolla bi-dimensionale, limitata da una curva chiusa semplice, inizialmente posta in un flusso potenziale viscoso. Si assume che l'influenza della gravità e delle forze inerziali sia trascurabile, cosicchè si può applicare l'approssimazione quasi-stazionaria. Riformulando il problema per le equazioni di Stokes con le opportune condizionial contorno sulla superficie libera in termini della funzione stress-stream bianalitica, e usando la tecnica delle trasformazioni conformi dipendenti dal tempoz(,t) di un disco unitario su un dominio di flusso non limitato incognito, viene derivato un sistema infinito di equazioni differenziali ordinarie per i coefficienti di Laurent diz(,t). Viene trovata una classe di soluzioni esatte per il caso in cui la parte principale della velocità complessa del flusso dominante all'infinito è una polinomiale, ed è discusso il problema della formazione di una bolla lenticolare.

     

The oseen flow of finite clusters of spheres and the wake effect on the drag factor

In this paper, the wake effect on drag factor in the axisymmetric Oseen flow of the finite clusters of equally spaced spheres with same size is studied. Putting the Oseen lets on the centres of all the spheres, the series solution of the problem is obtained. By truncating the infinite series and applying the collocation method to solve a set of the linear algebraic equations, the approximate solution of the Oseen flow of finite clusters of spheres and the drag factor for each sphere are presented. The effect of the sphere number and spacing on the drag factor of each sphere under different Reynolds numbers are calculated and the wake effect as well as the shielding effect and the end effect are revealed. The influence of various parameters on the effects is considered and compared with the corresponding results of the Stokes flow. The convergence of the method is also studied numerically in this paper.     

平板大攻角绕流升力和阻力系数的计算

二维平板或二维对称薄翼型大攻角绕流升力和阻力系数与攻角之间存在的函数关系一般用数据表格的形式给出。本文根据垂直平板绕流阻力实验数据和对称薄翼型全攻角绕流实验数据,分析得到了平板大攻角绕流总压力及其升力分量和阻力分量系数的近似计算公式。结果表明：平板总压力系数约等于攻角正弦值的2倍;总压力的阻力分量系数约等于攻角正弦值平方的2倍;升力分量系数约为攻角2倍的正弦值。计算结果与两组试验数据具有较好的一致性。     

液浮转子微陀螺筒间流动阻力矩参数影响

建立简化计算模型,将一种具有空间指数收敛特性的谱元-Fourier法引入复杂的液浮转子微陀螺筒间流动的直接数值模拟。对于柱坐标系下不可压缩Navier-Stokes方程,在子午面的径向和轴向采用Galerkin谱元法,进行标准化的Gauss-Lobatto-Legendre节点展开,在周向采用Fourier级数展开法。设定基准算例,计算分析了高径比、旋转速度和黏性系数等参数对内筒侧面和内筒端面阻力矩的影响。结果表明,两种阻力矩都随这些参数的增大而增大,且对总阻力矩的贡献量级相当。     

Energy-based decomposition of friction drag in turbulent square-duct flows

The generation of friction drag in turbulent duct flows has direct connection with statistical quantities and corresponding turbulence dynamics in the duct cross-section. In this study, we generalize the RD identity (Renard and Deck, 2016) to a ‘two-dimensional’ form which we exploit to decompose the mean friction drag in turbulent square-duct flows into contributions associated with viscosity, turbulence and cross-stream convection. The friction Reynolds number of the duct flows ranges from 220 to 2000. The scaling, spatial distribution and local normalization of the contributions to friction are investigated and compared with those in pipe and channel flows. As in other canonical flows, we find logarithmic growth of the turbulent contribution in contrast to the viscous one, the former thus becoming dominant at high enough Reynolds numbers. Whereas cross-stream convection has no net effect on friction, its contribution may be locally comparable to the other two, hence may be responsible for redistribution of friction along the duct perimeter.     

Optimal flow control for Navier–Stokes equations: drag minimization

Optimal control and shape optimization techniques have an increasing role in Fluid Dynamics problems governed by partial differential equations (PDEs). In this paper, we consider the problem of drag minimization for a body in relative motion in a fluid by controlling the velocity through the body boundary. With this aim, we handle with an optimal control approach applied to the steady incompressible Navier–Stokes equations. We use the Lagrangian functional approach and we consider the Lagrangian multiplier method for the treatment of the Dirichlet boundary conditions, which include the control function itself. Moreover, we express the drag coefficient, which is the functional to be minimized, through the variational form of the Navier–Stokes equations. In this way, we can derive, in a straightforward manner, the adjoint and sensitivity equations associated with the optimal control problem, even in the presence of Dirichlet control functions. The problem is solved numerically by an iterative optimization procedure applied to state and adjoint PDEs which we approximate by the finite element method. Copyright © 2007 John Wiley & Sons, Ltd.     

Analytical solution in the problem of constructing axisymmetric noses with minimum wave drag

An analytical procedure for determining the axisymmetric nose shapes that ensure minimum wave drag is developed. The solution is constructed as an improving variation on the conical shape. The target function is built up in an approximate form on the basis of the assumption that the relationship between the geometric and gasdynamic parameters is local in nature. It is shown that the optimal bodies are truncated power-law bodies with an exponent equal to 2/3. The bodies thus obtained are compared with those constructed in accordance with an exact formulation of the problem.     

An experimental study of drag on a single tube and on a tube in an array under two-phase cross flow

In this work, the drag coefficient and the void fraction around a tube subjected to two-phase cross flow were studied for a single tube and for a tube placed in an array. The drag coefficients were determined by measuring the pressure distribution around the perimeter of the tube. Single tube drag data were taken when the tube was held both rigidly and flexibly. The test tube was made of acrylic and was 2.2 cm in diameter and 20 cm in length. In the experiments, liquid Reynolds number ranged from 430 to 21,900 for the single tube and liquid gap Reynolds number ranged from 32,900 and 61,600 for the tube placed in a triangular array. Free stream void fraction was varied from 0 to 0.4. At low Reynolds numbers, the ratio of two-phase to single-phase drag coefficient is found to be a strong function of εGr/Re2. However, at high Reynolds numbers only void fraction is the important parameter. Empirical correlations have been developed for the ratio of two-phase drag on a single tube and on a tube placed in an array.     

A boundary integral method for two-dimensional (non)-Newtonian drops in slow viscous flow

A boundary integral method for the simulation of the time-dependent deformation of Newtonian or non-Newtonian drops suspended in a Newtonian fluid is developed. The boundary integral formulation for Stokes flow is used and the non-Newtonian stress is treated as a source term which yields an extra integral over the domain of the drop. The implementation of the boundary conditions is facilitated by rewriting the domain integral by means of the Gauss divergence theorem. To apply the divergence theorem smoothness assumptions are made concerning the non-Newtonian stress tensor. The correctness of these assumptions in actual simulations is checked with a numerical validation procedure. The method appears mathematically correct and the numerical algorithm is second order accurate. Besides this validation we present simulation results for a Newtonian drop and a drop consisting of an Oldroyd-B fluid. The results for Newtonian and non-Newtonian drops in two dimensions indicate that the steady state deformation is quite independent of the drop-fluid. The deformation process, however, appears to be strongly dependent on the drop-fluid. For the non-Newtonian drop a mechanical model is developed to describe the time-dependent deformation of the cylinder for small capillary numbers.