Numerical simulation of supersonic flow past reentry capsules

The flow fields over ARD (ESA's atmospheric reentry demonstrator), OREX (orbital reentry experiments) and spherically blunted cone-flare reentry configurations are numerically obtained by solving time-dependent, axisymmetric, compressible Navier–Stokes equations for freestream Mach numbers range of 1.2–6.0. The fluid dynamics are discretized in spatial coordinates employing a finite volume approach which reduces the governing equations to semi discretized ordinary differential equations. Temporal integration is performed using the multistage Runge–Kutta time-stepping scheme. A local time step is used to achieve steady-state solution. The numerical simulation is carried out on a structured grid. The flow-field features around the reentry capsule, such as bow shock wave, sonic line, expansion fan and recirculating flow in the base region are obtained. A good agreement is found between the calculated value of aerodynamic drag coefficient of the spherically blunted cone/fare reentry configuration with the experimental data. The effects of geometrical parameters, such as radius of the spherical cap, half cone angle, with sharp shoulder edge and with smooth shoulder edge on the flow-field have been numerically investigated for various reentry configuration which will be useful for optimization of the reentry capsule. PACS 47.11.Df, 47.40.Ki     

Internal-strain measurements of longitudinal pulses in conical bars

Longitudinal stress waves in a truncated 20-deg solid cone were investigated using embedded semiconductor strain gages. The cone, composed of an aluminum-filled epoxy, was struck normally at its small end with a 1/2-in.-diam steel ball traveling at a velocity of 170 jps. The results show the magnitude of the resulting stress wave to be nonuniform over a plane cross section perpendicular to the cone axis, the strain being greater at the center of the cone than near the surface, and the nonuniformity to increase with distance of travel from the impact end. The surface-strain measurements were compared with the one-dimensional theory of longitudinal waves in cones developed by Landon and Quinney as solved by Kenner and Goldsmith for a onehalf cycle sine-squared input pulse, and found to be in qualitative agreement with this theory, but to vary significantly in strain magnitude due to the strain nonuniformity over plane cross sections. The nonuniformity was compared with the Pochhammer-Chree theory for stress waves in cylindrical bars when that theory was evaluated for a cross section equivalent to the cone cross section. The trends of the deviations were similar, but the variations measured in the cone were consistently greater than that predicted by the theory.     

The determination of plough draught—Part I. prediction from soil and meteorological data with cone index as the soil strength parameter

Using cone index as an indication of soil strength, empirical equations are developed in accordance with soil mechanics theory to relate soil moisture content to plough draught. The plough draught equation comprises a quasi-static component dependent on cone index and a dynamic component which is a function of the soil specific weight, plough speed and mouldbard tail angle. It is further argued that the cohesive and frictional components of the cone penetration resistance can be predicted by means of a simple equation comprising a reciprocal function of the square of the soil moisture content and a linear function of the soil specific weight. The cone index equation explained 98% of the experimental data for threthree soils over a wide range of moisture contents. These empirical equations, together with a soil moisture model, provide a method of predicting plough draught directly from soil and meteorological data.     

Fluid inertia effects on a Controlled Stress Rheometer in its oscillatory mode

This paper is concerned with the effect of fluid inertia on experimental oscillatory stress data taken from a Controlled Stress Rheometer. A linear viscoelastic theory is developed which includes the effect of fluid inertia for cone and plate, parallel plate and concentric cylinder geometries. This theory is used to interpret dynamic data for both a slightly elastic and a highly viscoelastic fluid. It is shown that intertial effects are very small for both a cone and plate and parallel plate geometries. Inertial effects, however, can be important in the concentric cylinder geometry.     

On the integrability of the motion of a heavy particle on a tilted cone and the swinging Atwood machine

In the present note we study the motion of a particle under uniform gravity on a tilted smooth cone. We point out some integrable cases of that problem. In particular, we establish that motion on the cone with semi-vertical angle π/6 and one generator vertical is integrable. Moreover, we show that the swinging Atwood machine is equivalent to a special version of the motion on the tilted cone. This analogy makes it easier to analyze motion of SAM from a geometrical point of view.     

静力触探锥头阻力的近似理论与实验研究进展

锥头阻力在静力触探试验中扮演着十分重要的角色.从不同角度,对触探中锥头阻力的研究进行简要阐述,对承载力理论、空洞膨胀理论、应变路径法及运动点位错法等几种理论分析方法进行了回顾.另外,对数值分析和实验研究的进展情况进行了叙述.并对各种方法的适用性进行了比较.承载力理论虽然简单,但忽略了土的压缩性和探杆周围初始应力的增加,所以不能精确地模拟锥头的深层贯入.空洞膨胀理论提供了一个分析锥头阻力的简单而较精确的方法,它考虑了土的压缩性(或膨胀性)和锥头贯入过程中锥杆周围应力增加的影响.但这种方法是将锥头贯入与空洞膨胀之间做了一个等效模拟,所以不同的模拟方法,得到的结果差别较大.应变路径法能够有效解决饱和粘土中的不排水贯入,但不适用于砂土.运动点位错法因为考虑了部分排水,所以能较好地预测固结系数,但采用了线弹性分析,故位错法在其他方面的应用还需要大量的试验验证.有限元法在处理锥头贯入这类慢侵彻问题时缺乏一种很好的处理技术,导致它在进行破坏荷载计算时有显著的误差和数值计算困难.标定槽试验将在验证和建立锥头阻力与土的性能关系方面继续起到一个重要作用,但其结果需经过校正后才可应用到现场.最后对该领域的研究趋势进行了讨论.      

Features of the Application of Slender-Body Theory to the Calculation of Supersonic Water-Stream Cavitation Flow Past Cones

Features of the application of slender-body theory to the calculation of supersonic water-stream cavitation flows past cones are described. It is shown that the boundary condition on the cone edge in whose neighborhood Prandtl-Meyer flow develops cannot be satisfied within the framework of slender-body theory. The cone drag coefficients determined on the basis of slender-body theory are compared with the results of numerical calculations of conical flows.     

Bending analysis of functionally graded CNT reinforced doubly curved singly ruled truncated rhombic cone

The bending analysis of functionally graded carbon nanotube (CNT) reinforced doubly curved singly ruled truncated rhombic cone is investigated. In this study, a simple C⁰ isoparametric finite element formulation based on third order shear deformation theory is presented. To characterize the membrane-flexure behavior observed in a CNT reinforced truncated rhombic cone, a displacement field involving higher-order terms in in-plane fields is considered. The proposed kinematics field incorporates for transverse shear deformation and nonlinear variation of the in-plane displacement field through the thickness to predict the overall response of the CNT reinforced truncated rhombic cone in an accurate sense. The material properties of the CNT reinforced truncated rhombic cone are estimated according to the rule of mixture. The present model eliminates the need of shear correction factor and imposed zero-transverse shear strain at upper and lower surface of the truncated rhombic cone. The new feature in present model is simultaneous inclusion of twist curvature in strain field as well as curvature in displacement field that makes it suitable for moderately thick and deep truncated rhombic cone. The proposed new mathematical model is implemented in finite element code written in FORTRAN. The proposed model has been validated with analytical, experimental, and finite element results from the literature. This is first attempt to study bending response of CNT reinforced doubly curved singly ruled truncated rhombic cone. The effect of CNT distribution, boundary condition, loading pattern, and other geometric parameters are also examined.     

Second gradient electrodynamics: Green functions,wave propagation,regularization and self-force

In this work, the theory of second gradient electrodynamics, which is an important example of generalized electrodynamics, is proposed and investigated. Second gradient electrodynamics is a gradient field theory with up to second-order derivatives of the electromagnetic field strengths in the Lagrangian density. Second gradient electrodynamics possesses a weak nonlocality in space and time. In the framework of second gradient electrodynamics, the retarded Green functions, first-order derivatives of the retarded Green functions, retarded potentials, retarded electromagnetic field strengths, generalized Liénard–Wiechert potentials and the corresponding electromagnetic field strengths are derived for three, two and one spatial dimensions. The behaviour of the electromagnetic fields is investigated on the light cone. In particular, the retarded Green functions and their first-order derivatives show oscillations around the classical solutions inside the forward light cone and it is shown that they are singularity-free and regular on the light cone in three, two and one spatial dimensions. In second gradient electrodynamics, the self-force and the energy release rate are calculated and the equation of motion of a charged point particle, which is an integro-differential equation where the infamous third-order time-derivative of the position does not appear, is determined.     

On the accurate prediction of the wall-normal velocity in compressible boundary-layer flow

We consider a problem which arises in the numerical solution of the compressible two-dimensional or axisymmetric boundary-layer equations. Numerical methods for the compressible boundary-layer equations are facilitated by transformation from the physical (x, y) plane to a computational (ξ, η) plane in which the evolution of the flow is ‘slow’ in the time-like ξ direction. The commonly used Levy-Lees transformation results in a computationally well-behaved problem, but it complicates interpretation of the solution in physical space. Specifically, the transformation is inherently non-linear, and the physical wall-normal velocity is transformed out of the problem and is not readily recovered. Conventional methods extract the wall-normal velocity in physical space from the continuity equation, using finite-difference techniques and interpolation procedures. The present spectrally accurate method extracts the wall-normal velocity directly from the transformation itself, without interpolation, leaving the continuity equation free as a check on the quality of the solution. The present method for recovering wall-normal velocity, when used in conjunction with a highly accurate spectral collocation method for solving the compressible boundary-layer equations, results in a discrete solution which satisfies the continuity equation nearly to machine precision. As demonstration of the utility of the method, the boundary layers of three prototypical high-speed flows are investigated and compared: the flat plate, the hollow cylinder, and the cone. An important implication for classical linear stability theory is also briefly discussed.     

聚能装药对混凝土靶板的侵彻研究

系统开展了不同药型罩材料、不同锥角、不同壁厚的聚能装药在不同炸高下侵彻混凝土试验, 研究了罩材料、锥角、壁厚、炸高等结构参数对漏斗坑直径、侵彻孔洞直径、漏斗坑深度以及侵彻深度等参数的影响规律;应用空腔膨胀理论计算了混凝土靶体阻力, 采用改进的伯努利方程和两阶段空腔膨胀理论获得了混凝土靶板在侵彻体作用下的侵彻深度和孔洞直径, 理论结果与试验结果基本吻合;基于AUTODYN 软件平台, 采用与试验一致的聚能装药结构, 开展了57 种工况下侵彻体成形过程的数值模拟研究, 并对其中典型工况的侵彻混凝土过程进行了数值模拟, 计算所得孔洞直径和侵彻深度与试验结果吻合较好, 在此基础上深入探讨了聚能装药作用下混凝土漏斗坑的形成机理, 分析表明, 铝药型罩的开坑机理不同于钢和铜药型罩.      

Far field scattering from bodies of revolution

Summary By use of approximations based on physical reasoning radar cross-section results for bodies of revolution are found. In the Rayleigh region (wavelength large with respect to object dimensions) approximate solutions are found. Examples given include a finite cone, a lens, an elliptic ogive, a spindle and a finite cylinder. In the physical optics region (wavelength very small with respect to all radii of curvature) Kirchhoff theory and also geometric optics can be used. When the body dimensions are only moderately large with respect to the wavelength, Fock or Franz theory can be applied, and examples of the circular and elliptic cylinder are presented. In the region where some dimensions of the body are large with respect to the wavelength and other dimensions are small with respect to the wavelength, special techniques are used. One example, the finite cone, is solved by appropriate use of the wedgelike fields locally at the base. Another example is the use of traveling wave theory for obtaining approximate solutions for the prolate spheroid and the ogive. Other results are obtained for cones the base perimeter of which is of the order of a wavelength by using known results for rings of the same perimeter.This paper with minor revision is as the author presented it at the URSI XIIth General Assembly in Boulder, Colorado August 22–September 5, 1957.     

高压捕获翼位置设计方法研究

高压捕获翼构型是一种合理利用机体/上置翼(简称捕获翼)间的耦合关系提高飞行器升力,进而大幅提高升阻比的高速飞行器新概念构型.基于其设计原理,捕获翼的位置与机体压缩激波和自身二次压缩激波的位置均直接相关,一般难以利用理论方法直接获得.针对这一问题,本文运用均匀实验设计方法在设计空间内获取样本点并利用计算流体力学分析和迭代获得其设计位置,之后通过构造代理模型建立捕获翼位置与设计参数间的模拟映射关系,进而发展了一种捕获翼位置设计的有效方法.在方法研究基础上以锥体-捕获翼组合构型作为实例对其进行验证.结果表明,该方法可在较大设计空间范围内准确判定捕获翼的设计位置.此外,针对这一构型还开展了基于代理模型的设计参数单因素分析.发现在设计空间内,前缘压缩角、来流马赫数、和捕获翼钝化半径等3个关键参数均与捕获翼位置呈单调正比例关系.      

Plasma mitigation of shock wave: experiments and theory

Two types of plasma spikes, generated by on-board 60 Hz periodic and pulsed dc electric discharges in front of two slightly different wind tunnel models, were used to demonstrate the non-thermal plasma techniques for shock wave mitigation. The experiments were conducted in a Mach 2.5 wind tunnel. (1) In the periodic discharge case, the results show a transformation of the shock from a well-defined attached shock into a highly curved shock structure, which has increased shock angle and also appears in diffused form. As shown in a sequence with increasing discharge intensity, the shock in front of the model moves upstream to become detached with increasing standoff distance from the model and is eliminated near the peak of the discharge. The power measurements exclude the heating effect as a possible cause of the observed shock wave modification. A theory using a cone model as the shock wave generator is presented to explain the observed plasma effect on shock wave. The analysis shows that the plasma generated in front of the model can effectively deflect the incoming flow; such a flow deflection modifies the structure of the shock wave generated by the cone model, as shown by the numerical results, from a conic shape to a curved one. The shock front moves upstream with a larger shock angle, matching well with that observed in the experiment. (2) In the pulsed dc discharge case, hollow cone-shaped plasma that envelops the physical spike of a truncated cone model is produced in the discharge; consequently, the original bow shock is modified to a conical shock, equivalent to reinstating the model into a perfect cone and to increase the body aspect ratio by 70%. A significant wave drag reduction in each discharge is inferred from the pressure measurements; at the discharge maximum, the pressure on the frontal surface of the body decreases by more than 30%, the pressure on the cone surface increases by about 5%, whereas the pressure on the cylinder surface remains unchanged. The energy saving from drag reduction is estimated to make up two-thirds of the energy consumed in the electric discharge for the plasma generation. The measurements also show that the plasma effect on the shock structure lasts much longer than the discharge period.

---

     

Cone resistance in sand by incremental method

A new method has been developed for the determination of cone resistance under drained conditions. Numerical methods are used for the solution of the differential equations of plasticity theory for soils and for the determination of the stress states in the soil produced by the penetration of the cone. It is assumed that the stresses produced by the penetration of the cone remain ‘locked in’ the soil and constitute boundary conditions for further penetration. The computation starts with the cone base at the surface and is continued by successively incrementing the depth by a small amount. Charts are given for the computation of cone resistance in sands for various friction angles. The importance of the effect of the shear stresses generated at the surface of the cone and characterized by the interface friction angle, δ, is discussed in detail.     

基于微纳米压入法提取无铅焊料合金弹塑性力学参数的研究

本文通过微纳米压入法结合数值模拟研究了无铅焊料合金SnAg3.5 的弹塑性力学性能，分别采用圆柱形压头及两种不同锥角压头对无铅焊料合金进行压入测试：基于圆柱形压头测试过程中接触面积恒定的特点得到了无铅焊料的弹性模量，进一步采用塑性应变梯度理论对两种锥角压头的测试结果予以修正并通过数值模拟反分析得到相应的特征应力值，同时基于压入特征塑性应变与压头锥角的关系式得到两种不同锥角压头下的特征应变值，在此基础上经求解方程组得到焊料合金的初始屈服应力与应变强化因子，进而得到了焊料合金的幂强化弹塑性本构关系．该方法剔除了压入尺度效应的影响并保证了所得本构关系的唯一性，给出了一种通过原位压入测试表征金属材料弹塑性力学性能的有效方法．     

Technology showcase electronic penetrometer for field tests

The design and construction of a cone penetrometer for field use is described. The penetrometer is mounted on a frame to which two wheels are fitted for towing the device from one test site to the next. An electric motor drives the penetrometer during the test and simultaneous readings are performed of both cone depth and penetration resistance. The device weighs 70 kg and the maximum force on the cone typically exceeds 1.2 kN. The acquired data are transferred to a microcomputer and after each penetration the result is presented to the operator in graphic form. The microcomputer is equipped with a fairly large RAM memory which serves as storage medium for data during measuring work. Afterwards, data are easily moved to a mainframe computer for analysis.     

Stability analysis and transition prediction of hypersonic boundary layer over a blunt cone with small nose bluntness at zero angle of attack

Stability and transition prediction of hypersonic boundary layer on a blunt cone with small nose bluntness at zero angle of attack was investigated. The nose radius of the cone is 0.5 mm; the cone half-angle is 5°, and the Mach number of the oncoming flow is 6. The base flow of the blunt cone was obtained by direct numerical simulation. The linear stability theory was applied for the analysis of the first mode and the second mode unstable waves under both isothermal and adiabatic wall condition, and eN method was used for the prediction of transition location. The N factor was tentatively taken as 10, as no experimentally confirmed value was available. It is found that the wall temperature condition has a great effect on the transition location. For adiabatic wall, transition would take place more rearward than those for isothermal wall. And despite that for high Mach number flows, the maximum amplification rate of the second mode wave is far bigger than the maximum amplification rate of the first mode wave, the transition location of the boundary layer with adiabatic wall is controlled by the growth of first mode unstable waves. The methods employed in this paper are expected to be also applicable to the transition prediction for the three dimensional boundary layers on cones with angle of attack.     

NUMERICAL INVESTIGATION OF EVOLUTION OF DISTURBANCES IN SUPERSONIC SHARP CONE BOUNDARY LAYERS

The spatial evolution of 2-D disturbances in supersonic sharp cone boundary layers was investigated by direct numerical simulation (DNS) in high order compact difference scheme. The results suggested that, although the normal velocity in the sharp cone boundary layer was not small, the evolution of amplitude and phase for small amplitude disturbances would be well in accordance with the results obtained by the linear stability theory (LST) which supposes the flow was parallel. The evolution of some finite amplitude disturbances was also investigated, and the characteristic of the evolution was shown. Shocklets were also found when the amplitude of disturbances increased over some value.