Optimum form of lifting bodies at hypersonic speeds

The variational problem of the form of bodies with minimum drag for given lift force, volume, and other constraints in general leads to a second-order partial differential equation even for the simplest methods of drag calculation (Newton law and averaged friction coefficient). The solution of this equation is not justified; in its place an approximate solution is suggested which consists of: a) selection of a scheme characterized by certain parameters which are determined from the solution of the extremal problem, b) determination of the optimal surface form for the selected scheme with the aid of the system of ordinary Euler equations. This paper presents a comparison of the body schemes with minimum drag and maximum L/D and presents the solution of several variational problems.At the present time we have quite complete information on the form of minimum-drag bodies for zero lift (nonlifting bodies), and both approximate and quite rigorous methods are known for solving the corresponding variational problems. This cannot be said at all of the form of lifting bodies, for which the requirements are numerous, differing essentially for vehicles of different application, and are generally not limited to a single flight regime. Account for all the mandatory requirements in solving the variational problems is not possible; therefore in the majority of cases these solutions do not yield answers which are directly suitable in practice; rather they yield limiting estimates.The natural tendency to utilize for lifting bodies the axisymmetric form which is customary for nonlifting bodies leads to the study of axisymmetric bodies at angle of attack, axisymmetric bodies with skewed base, sections of axisymmetric bodies cut by planes, etc. In order to obtain a broader view of the optimum forms of lifting bodies we must, obviously, drop the limitation to axisymmetric bodies and bodies with similar cross sections. However, in the case of an arbitrary extremal surface the Euler equation is a second-order partial differential equation, and its simple solution is difficult. In practice it seems wise to solve those variational problems whose Euler equation may be reduced to a system of ordinary differential equations. Therefore, we propose the following method for selecting the optimum forms: a) we select a scheme, a form, which is formed by a set of planes and cylindrical, conical, spherical surfaces and which is defined by parameters that are found from the solution of the extremal problem; b) for the selected scheme the generators of the scheme surface are found from the solution of the variational problem.For the calculation of the air pressure on the body surface we use the empirical Newton law, which yields in the majority of cases results which are very close to the results of the more rigorous methods.It is assumed that the pressure may vanish only at the trailing edge of the body. The frictional drag coefficient, averaged over the body surface, is assumed to be independent of the body shape. In the case of a body of simple form the hottest portion is the frontal portion and account for the thermal protection requirements reduces to the selection of suitable dimensions of this portion of the body. In the general case the problem is stated as follows: find the form of the minimum-drag body for a given lift force, volume, length, and other conditions. To the particular case of the body with maximum L/D corresponds the value of the Lagrange multiplier =–1/k.All the results of calculations presented in the paper are intended only to illustrate the method. After the present paper was submitted for publication, another study [3] appeared which also proposes a method for determining the optimal parameters.     

Axisymmetric jets impinging on porous walls

The flow of axisymmetric turbulent jets impinging on porous walls has been studied experimentally. It is shown how the overall flow structure depends on the porosity of the surface. For high porosities (open area ratios, β, in excess of around 40% say) the porous wall, or screen, leads to a sudden increase in jet width and decrease in mean and fluctuating velocities, a direct consequence of the momentum flux extracted because of the screen drag. Lower porosities can lead to the appearance of radial wall jets on the upstream side of the screen but, in contrast to the corresponding case of planar jet impingement (Cant et al. in Exp Fluids 32:16–26, 2002), such wall jets never occur on the downstream side. The axial downstream velocities thus remain positive for all porosities. Jet growth rates for are initially increased by the screen, but once β≤0.4 momentum extraction by the screen is virtually complete, so that velocities become very small. Again, unlike in the corresponding planar case (for β≈0.4), recirculating regions upstream of the screen never occur. A simple argument is suggested to explain the fundamental differences in flow behaviour between planar and axisymmetric jet impingement onto porous screens and it is concluded that in the latter case the effects of the screen are generally more benign and unsurprising. Nonetheless, these axisymmetric flows, like the corresponding planar ones, provide a serious challenge for computational modelling.     

A variational approach to non-steady non-newtonian flow in a circular pipe

The transient response of a non-Newtonian power-law fluid to several assumed forms of pressure pulse in a circular tube is analysed by the semi-direct variational method of Kanntovorich. Velocity profiles are shown for several power-law indices, and by comparing the results for the Newtonian case with the exact solution given by Szymanski, it is observed that the results are good to 5%. More accurate solutions have been found for the case involving Newtonian fluid flow. New results are reported concerning the effect of a triangular pressure pulse on the development and transient response of the flow field of a non-Newtonian fluid.     

Fountain flow of pseudoplastic and viscoplastic fluids

Numerical simulations have been undertaken for the benchmark problem of fountain flow present in injection-mold filling. The Finite Element Method (FEM) is used to provide numerical results for both cases of planar and axisymmetric domains under steady-state conditions. The Herschel–Bulkley model of viscoplasticity is used, which reduces with appropriate modifications to the Bingham, power-law and Newtonian models. The present results extend previous ones regarding the shape of the front, which is essential in correctly capturing the flow field. In particular the centreline front position is found as a function of the dimensionless power-law index (in the case of pseudoplasticity) and the dimensionless yield stress (in the case of viscoplasticity). The pressures from the simulations have been used to compute the excess pressure losses in the system (front pressure correction or exit correction). Both shear-thinning and shear-thickening lead to more extended front positions relative to the Newtonian values, which are 0.895 for the planar case and 0.835 for the axisymmetric one. Viscoplasticity leads also to more extended front positions as the dimensionless yield stress goes from zero (Newtonian behaviour) to higher values of the yield stress. In both cases of non-Newtonian behaviour, the front tends to follow the development of the fully developed Poiseuille velocity profile, which tends towards a plug-like profile at the extreme cases of non-Newtonianness. The front pressure (exit) correction increases monotonically with the decrease in the power-law index and the increase in the dimensionless yield stress.     

Effects of Entropic Wave in Vibroacoustic Problem Using Boundary Element Analysis

A variational formulation for a vibroacoustic problem of a membrane and a viscothermal fluid is investigated in this paper. The formulation combines a variational formulation by integral equations of the fluid, that takes into account the acoustic and entropic waves coupling, with a variational formulation of the membrane. The formulation has been implemented numerically for the problems with axisymmetric geometry. The numerical results are compared to the analytical solution for a circular membrane coupled to a cylindrical cavity filled with air. These results show the validity of numerical implementation and illustrate the thermal effects of air on the membrane-cavity system modes in the micro cavities cases.     

Macroscopic strain potentials in nonlinear porous materials

By taking a hollow sphere as a representative volume element (RVE), the macroscopic strain potentials of porous materials with power-law incompressible matrix are studied in this paper. According to the principles of the minimum potential energy in nonlinear elasticity and the variational procedure, static admissible stress fields and kinematic admissible displacement fields are constructed, and hence the upper and the lower bounds of the macroscopic strain potential are obtained. The bounds given in the present paper differ so slightly that they both provide perfect approximations of the exact strain potential of the studied porous materials. It is also found that the upper bound proposed by previous authors is much higher than the present one, and the lower bounds given by Cocks is much lower. Moreover, the present calculation is also compared with the variational lower bound of Ponte Castañeda for statistically isotropic porous materials. Finally, the validity of the hollow spherical RVE for the studied nonlinear porous material is discussed by the difference between the present numerical results and the Cocks bound.     

The free energy of mixing for n-variant martensitic phase transformations using quasi-convex analysis

The construction of effective models for materials that undergo martensitic phase transformations requires usable and accurate functional representations for the free energy density. The general representation of this energy is known to be highly non-convex; it even lacks the property of quasi-convexity. A quasi-convex relaxation, however, does permit one to make certain estimates and powerful conclusions regarding phase transformation. The general expression for the relaxed free energy is however not known in the n-variant case. Analytic solutions are known only for up to 3 variants, whereas cases of practical interests involve 7-13 variants. In this study we examine the n-variant case utilizing relaxation theory and produce a seemingly obvious but very powerful observation regarding a lower bound to the quasi-convex relaxation that makes practical evolutionary computations possible. We also examine in detail the 4-variant case where we explicitly show the relation between three different forms of the free energy of mixing: upper bound by lamination, the Reuß lower bound, and a lower estimate of the -measure bound. A discussion of the bounds and their utility is provided; sample computations are presented for illustrative purposes.     

Features of flow in the region of the stagnation point in supersonic flow past power-law bodies

Bodies of power-law shape with a generator of the form y=x, 0.5<<1, have at the apex a vertical tangent (like blunt bodies) and infinite curvature (like sharp ones). Supersonic flow past axisymmetric bodies with power-law longitudinal profile has been studied by Radvogin in [1, 2], where on the basis of the analysis of a large quantity of numerical calculations certain empirical laws of similarity were formulated. It follows from these relations that the position and form of the subsonic sector of the shock wave are not determined by the singularity in the body's profile, but by the position on the body of the sonic point and its bluntness in relation to the cone of the critical opening angle. In O. V. Titov's work the results obtained in [1, 2] are confirmed analytically, but here it is assumed that the curvature of the shock wave and the second derivative of the curvature with respect to the longitudinal curvilinear coordinate are finite at the apex. This assumption imposes a limitation on the contour of the body; if it is satisfied, the curvature of the profile past which the flow takes place is also finite. Therefore it is natural to consider the case when this assumption is not made. In this paper we study the flow past bodies of power-law shape with shock waves of infinite curvature at the apex and finite curvature but an infinite second derivative of the curvature with respect to the longitudinal coordinate.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 138–142, May–June, 1985.     

On the nonlinear dynamic buckling mechanism of autonomous dissipative/nondissipative discrete structural systems

Summary Nonlinear dynamic buckling of nonlinearly elastic dissipative/nondissipative multi-mass systems, mainly under step load of infinite duration, is studied in detail. These systems, under the same loading applied statically, experience a limit point instability. The analysis can be readily extended to the case of dynamic buckling under impact loading. Energy, topological and geometrical aspects for the total potential energyV, which is constrained to lie in a region of phase-space whereV0, allow conclusions to be drawn directly regarding dynamic buckling. Criteria leading to very good, approximate and lower/upper bound dynamic buckling estimates are readily established without solving the highly nonlinear set of equations of motion. The theory is illustrated with several analyses of a two-degree-of-freedom model.     

The non-newtonian creeping flow around a rotating spheroidal spindle

The rotational flow in a region bounded by two axisymmetric closed surfaces has been studied by variational methods. The approach is described in detail for the global power-law representation of a viscosity function. Numerical, as well as simple, analytical results applicable to immersion rotational viscometry are given for the spheroidal geometry of bounding surfaces.     

Analysis of surface crack on forward extruding bar during axisymmetric cup-bar combined extrusion process

In this paper, the kinematically admissible velocity field with surface crack on forward extruding bar is put forward during the axisymmetric cup-bar combined extrusion process, in accordance with the results of model experiments.On the basis of velocity field, the necessary condition for surface crack formation on the forward extruding bar is derived, with the help of upper bound theorem and the minimum energy principle. Meanwhile, the relationships between surface crack formation and combination of reduction in area for the part of forward and backward extursions (ε_b,ε_f) relative residual thickness of billet (T/R₀),frictional factor (m) or relative land length of ram and chamber (l_b|R₀,l_f|R₀) are calculated during the extrusion process. Therefore, whether the surface crack on forward exturding bar occurs can be predicted before extruding the lower-plasticity metals for axisymmetric cup-bar combined extrusion process.The analytical results agree very well with experimental results of aluminium alloy LY12 (ASTM 2024) and LC4 (ASTM 7075).     

Approximate integration op nonstationary equations of a diffusion or thermal boundary layer

An approximate method is proposed for integrating the nonstationary equations of a diffusion or thermal boundary layer using the known steady solution in the planar or axisymmetric case. It is shown that the proposed method is exact in problems involving mass or heat transfer of reacting drops and bubbles in a laminar flow of a viscous incompressible fluid and also particles moving in an ideal fluid. An integral equation is obtained for the local diffusion or heat flux in the case of abrupt activation of a reaction on the surface of a particle.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 87–92, September–October, 1982.     

Study of convection in a liquid layer with vibrating forces

We consider convection in a gravity force field in a liquid enclosed in a vessel which is vibrating along the vertical axis according to the lawa/ sin t ( ).Time-averaged convection equations describing the basic motion in first approximation are derived in [1]. In addition, criteria are introduced in [1] which determine the initiation of convection in this case, and a model problem is considered: the case of spatially periodic disturbances. It is found that in this case high frequency vibrations stabilize the state of relative rest.In the present paper we consider convection in a liquid layer between two horizontal planes whose equations are z=±l/2.The temperatures T₁ and T₂ are specified on the upper and lower surfaces, respectively. It is shown (§1) that for small values of the vibrational criterion the principle of stability variation is satisfied. In this case there is a variational principle (§2) from which it follows that high frequency vibrations prevent the occurrence of convection in the horizontal liquid layer. With the aid of the variational principle a calculation is made of the dependence of the values of Rayleigh criterion on the vibrational criterion.The author wishes to thank I. B. Sinomenko and V. I. Yudovich for their continued interest in this study.     

Optimization of the shape of a body with respect to two criteria: Radiation heat flux and wave drag

The numerical solution of the two-criteria variational problem of the body contour with minimum radiation heat flux and wave drag is obtained in the class of axisymmetric and plane slender bodies in hypersonic flow. Solutions obtained using the Pareto, ideal point and minimax methods are compared. It is shown that in the class of axisymmetric slender bodies the optimum body gives a decrease in the radiation heat flux as compared with a cone of up to 15% for the Pareto method, up to 13% for the ideal point method, and up to 5% for the minimax method. A solution is also obtained in the subclass of power-law slender bodies and it is shown that the optimum power-law bodies are inferior, as compared with the optimum bodies from the general class of such bodies, in reducing both radiation heating and resistance.     

Steady flow of a power-law non-Newtonian fluid across an unconfined square cylinder

A two-dimensional flow of a non-Newtonian power-law fluid directed normally to a horizontal cylinder with a square cross section is considered in the present paper. The problem is investigated numerically with a finite volume method by using the commercial code Ansys Fluent with a very large computational domain so that the flow could be considered unbounded. The investigation covers the power-law index from 0.1 to 2.0 and the Reynolds number range from 0.001 to 45.000. It is found that the drag coefficient for low Reynolds numbers and low power-law index (n ≤ 0.5) obeys the relationship C_D = A/Re. An equation for the quantity A as a function of the power-law index is derived. The drag coefficient becomes almost independent of the power-law index at high Reynolds numbers and the wake length changes nonlinearly with the Reynolds number and power-law index.     

The axisymmetric problem with initial conditions for the equations of motion of an ideal incompressible fluid

The present paper presents a proof of the existence and uniqueness theorem for the solution of the axisymmetric problem with initial conditions for the Euler equations in the case of an incompressible fluid. We consider the case of the nonporous wall, and also the transpiration problem in the formulation given in [1]. Global unique solvability is proved for assumptions only on the smoothness of the conditions and for all values of the time t. The existence theorem for a small time segment in the case of a nonporous wall has been proved for the general three-dimensional problem in [2, 3]. For the proof we use a method analogous to that developed in [1] for planar flows. The a priori estimate of the vorticity which is used in the present study was obtained previously in [4],The author wishes to thank V. I. Yudovich for continued interest in the study and many valuable suggestions.     

Variational Proof of the Existence of the Super-Eight Orbit in the Four-Body Problem

Using the variational method, Chenciner and Montgomery (Ann Math 152:881–901, 2000) proved the existence of an eight-shaped periodic solution of the planar three-body problem with equal masses. Just after the discovery, Gerver numerically found a similar periodic solution called “super-eight” in the planar four-body problem with equal mass. In this paper we prove the existence of the super-eight orbit by using the variational method. The difficulty of the proof is to eliminate the possibility of collisions. In order to solve it, we apply the scaling technique established by Tanaka (Ann Inst H Poincaré Anal Non Linéaire 10:215–238, 1993), (Proc Am Math Soc 122:275–284, 1994) and investigate the asymptotic behavior of a binary collision.     

Variational problem of optimum side wall design for a narrow three-dimensional nozzle

The optimum design of the side walls of the supersonic section of a three-dimensional nozzle with two planes of symmetry is considered in the narrow channel model approximation, which reduces three-dimensional to two-dimensional flow. This nozzle realizes maximum thrust for given sonic or supersonic inlet flow, upper and lower walls, maximum permissible length and pressure outside the nozzle. In general, an approximate solution of the variational problem can be obtained by the indeterminate control contour method [1]. For nozzles with nonexpanding end sections of the upper and lower walls this is a rigorous solution. Numerical algorithms, based on the method of characteristics, for constructing the optimum, side walls and calculating the flow in narrow channels are developed in the formulation adopted using the optimality conditions found, which generalize the wellknown conditions for plane and axisymmetric configurations [1]. In addition, the three-dimensional supersonic flow in the nozzles thus designed has been calculated in accordance with a shock-capturing marching scheme [2], which for the uniform grids employed in the calculations gives a second-order approximation. A rather complex relation is established between the thrust of the optimum configurations constructed and the shape of their inlet cross sections.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.2, pp. 102–112, March–April, 1992.The authors are grateful to L. E. Sternin for drawing their attention to the problem and to V. A. Vostretsova for assisting with the work.     

The end correction for power-law fluids in the capillary rheometer

The finite element scheme developed by Nickell, Tanner and Caswell is used to compute the entry and exit losses for creeping flow of power-law fluids in a capillary rheometer. The predicted entry losses for a Newtonian fluid agree well with available experimental and theoretical results. The entry losses for inelastic power-law fluids increased with decreasing flow behaviour index and show an increasing deviation from available upper bound results as the flow behaviour index in the power-law decreases.The exit losses are found to be finite for inelastic power-law fluids and increase as the flow behaviour index decreases. The predicted die swell for Newtonian fluids agrees well with the available experimental data while the influence of shear thinning is to reduce the die swell.The end correction which is the sum of the entry and exit losses relative to twice the viscometric wall shear stress varies from 0.834 for n = 1 to 2.917 for n = 1/6. This figure reaches a very high value as n tends to zero. The experimental variation in the Couette correction factor in capillary rheometry is explained in terms of the shear thinning characteristics of the fluid. It is concluded that the exit flow is not viscometric, contrary to a common assumption.     

Equilibrium shapes and stability of rotating drops

Direct variational methods are used to obtain nonaxisymmetric equilibrium shapes of a rotating drop. The results are given of experiments on the interaction of freely floating drops of viscous fluid. It is shown that if the dimensionless angular momentum of the system is > 3.4, then a sequence of three-dimensional shapes similar to ellipsoids bifurcates from the sequence of axisymmetric shapes The stability of these shapes is studied in Poincaré's scheme.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 13–20, July–August, 1982.