A partially implicit method for simulating viscous aerofoil flows

A partially implicit method for the unsteady compressible Navier-Stokes equations is developed. The method is based on an explicit treatment of streamwise fluxes and an implicit treatment of normal fluxes. This leads to a linear system which is generated by an efficient finite difference procedure and which is block pentadiagonal. The method is tested on a shock-induced oscillatory flow over an aerofoil. Parallel implementations of an explicit, fully implicit and partially implicit method are investigated.     

一种考虑界面相物性影响的界面单元法

界面是由复杂的界面相简化而成的,界面破坏实际是界面相材料的破坏。数值计算为了方便,如经典模型和内聚力模型等,都把很薄的界面相作无厚度化处理。导致只能考虑界面的面力,而无法考虑界面相内的应力(平行于界面方向的应力)。使界面失效准则先天性地排除了界面相内部应力的影响,从界面相材料失效机理的角度来看这是不够严谨的。本文将界面相材料等效为一种弹性连续体,由界面本构关系推导得到了一种新的界面单元。该单元具有界面参数易确定、对界面相物性可以进行等效描述等优点。通过商用有限元软件ABAQUS和用户子程序UEL实现了数值分析,并与直接物理模型的数值模拟结果进行对比,证明了本方法的简便及准确性。通过对不同界面相厚度结构的进一步分析,探讨了本文方法的可行范围。     

Non-linear argumental oscillators: Stability criterion and approximate implicit analytic solution

The behaviour of a space-modulated, so-called “argumental” oscillator, is studied. The oscillator is submitted to an external harmonic force, which is amplitude-modulated by the oscillator's position in space. An analytic expression of a stability criterion is given. Using the averaging method, an integrating factor and a Van der Pol representation in the (amplitude, phase)-space, an exact implicit analytic solution is given when there is no damping, and an approximate implicit analytic solution is given when there is damping, allowing the plotting of the separatrix curve. An attractor is identified.     

Interaction of an Ellipsoidal Inclusion with an Elliptic Crack in an Elastic Material under Triaxial Tension

The interaction of an elastic ellipsoidal inclusion with an elliptic crack in an infinite elastic medium under triaxial loading is analyzed. The stress state in the elastic space is represented as a superposition of the principal state and perturbed states, which are due to the presence and interaction of the inclusion and the crack. The analytical solution of the problem is found using the method of equivalent inclusion, the potential of an inhomogeneous ellipsoid, and a system of harmonic functions for an elliptic crack. The effect of triaxial loading on the stress intensity factors is analyzed     

可压缩流体的Rayleigh-Taylor和Kelvin-Helmholtz不稳定性

在状态方程为压力是密度的任意单值函数形式情况下，运用小扰动分析和奇异摄动法，给出了流体微扰方程渐近解和界面不稳定性的色散关系. 分析表明：对Rayleigh-Taylor不稳定性，在重力场作用下流体可压缩性形成的密度分布是致稳因素；而扰动流体的膨胀收缩效应助长不稳定性的发展；上层重流体的可压缩性是稳定因素，下层轻流体可压缩性是失稳因素. 而对Kelvin-Helmholtz不稳定性，流体可压缩性助长扰动的发展，是不稳定因素.     

Deformation of an ellipsoidal ferrogel sample in a uniform magnetic field

The elongation of a ferroelastic material sample (whose initial shape is a sphere or an ellipsoid of revolution) under the action of an external magnetic field is studied in an in approximation of small strains. For a sphere, there is a classical estimate obtained under the assumption that elongating in the direction of the field, it becomes a spheroid and the stress and strain fields remain uniform. In the present calculation, it is assumed that the body is an ellipsoid (a sphere in a particular case) only in the absence of an external field; the shape of the sample in the presence of a field is not specified in advance but is found from the condition of balance of surface forces (elastic and magnetic). For the spherical case, the problem is solved exactly: it is shown, that the contour of the deformed body is described by a third-order algebraic equation. The case where the initial configuration is an ellipsoid of revolution is studied numerically. It is shown that in all versions, the refined solution leads to an appreciable increase in the elongation of the sample compared to the classical estimate. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 153–164, May–June, 2005.     

Experimental Investigation of the Interaction between an Electric Arc and a Gas Flow

The results of investigating the shape of an electric arc carrying an electric current of less than 100 A and exposed to an air, nitrogen, or argon stream at a stream velocity of from zero to 20 m/s, when the interelectrode gap is less than 20 mm, are reported. The shape of the arc is qualitatively investigated as a function of the gas type, blowing velocity, and interelectrode gap. It is shown that an arc burning in a channel with insulated walls is much more stable under these conditions than an arc with the same parameters burning in free space.     

Solving axisymmetric dynamic problems for reinforced shells of revolution on an elastic foundation

The dynamic behavior of reinforced shells of revolution in an elastic medium is modeled. Pasternak’s model is used. A problem of vibration of discretely reinforced shells of revolution is formulated and a numerical algorithm is developed to solve it. Results from an analysis of the dynamic behavior of a reinforced spherical shell on an elastic foundation are presented as an example Translated from Prikladnaya Mekhanika, Vol. 45, No. 2, pp. 99–106, February 2009.     

Deformations of an elastic pipe submitted to gravity and internal fluid flow

This paper describes the deformation of an elastic pipe submitted to gravity and to an internal fluid flow. The pipe is clamped horizontally at one end and free at the other end. As the fluid velocity increases, the shape changes from an elastic beam deflected by its own weight towards an horizontal position. The shape of the pipe is characterized experimentally and is compared with a theoretical model based on the Euler–Bernoulli approximation and the conservation of the fluid momentum. We study how the determination of the pipe deformation provides an estimation of the conveyed fluid flow. Finally, the vertical force produced by the conveyed fluid to lift off a mass is deduced.     

Effects of an axial flow component on the Honji instability

This paper is concerned with a numerical study of the three-dimensional Honji instability that can arise in an oscillatory flow impinging on a circular cylinder. It is well known that when the fluid motion far from the cylinder is perpendicular to its axis then the flow is liable to a three-dimensional breakdown via this instability which initially appears as an axially periodic mushroom-like structure attached to the surface of the cylinder. Here the focus is on examining the Honji instability under an oblique inflow. The obliqueness of the free stream is represented by an angle of attack through introducing an axial flow component. It is found that the Honji mode is suppressed by increasing the axial flow component, and when this component is sufficiently large the instability mechanism is no longer operative so that all that remains is a featureless two-dimensional columnar flow. At smaller values of angle of attack, though the Honji structure remains, it is deformed by the axial flow component. The developed two-layer near-cylinder vortical structures can be related to the energy and momentum transfer between the two layers.     

Self-Similar Unsteady Viscous Flow

Four examples of self-similar flows of a viscous fluid are considered: separated flow over an expanding plate immersed in an unbounded unsteady viscous flow, the evolution of the velocity field induced by a vortex-source, the flow near an unsteadily moving permeable flat plate, and the flow near an unsteadily rotating disc. For the first example, a numerical solution is constructed. For the next two examples, an analytical solution is found, while the solution of the last problem is reduced to a system of ordinary differential equations.     

Uniform difference scheme for a singularly perturbed linear 2nd order hyperbolic problem with zeroth order reduced equation

In this paper a singularly perturbed linear second order hyperbolic problem with zeroth order reduced equation is discussed. Firstly, an energy inequality of the solution and an estimate of the remainder term of the asymptotic solution are given. Then an exponentially fitted difference scheme is developed in an equidistant mesh. Finally, uniform convergence in small parameter is proved in the sense of discrete energy norm.     

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.     

基于实测挠度、转角和曲率的细长梁分段抗弯刚度识别研究

静载试验是获得结构真卖特性的重要方法和途径，具有直观、可靠的特点。本文基于静载试验实测的挠度、转角和曲率，建立目标函数，通过优化的方法，实现结构分段截面抗弯刚度参数的识别。研究了优化目标函数的建立、优化算法的选取与实现；通过数值算例和一片简支梁静载试验实例分析，对方法的正确性与可行性进行了验证，所得结果可用于结构损伤识别和安全性能评估。     

Unsteady flow of a second-order fluid between concentric cylinders

The unsteady motion of an incompressible second-order fluid contained between two finite coaxial cylinders is examined when the outer cylinder is held fixed while the inner cylinder is constrained to execute an arbitrary angular velocity. A solution is obtained in closed form with the aid of transforms and an expression is obtained for the couple experienced. The particular case of a periodic angular velocity is then examined and some numerical work done. There is a marked difference between the results obtained and their classical counterparts.     

MICROMECHANICAL DAMAGE MODEL FOR ROCKS AND CONCRETES SUBJECTED TO COUPLED TENSILE AND SHEAR STRESSES

Based on the analysis of the deformation in an infinite isotropic elastic matrix with an embedded elliptic crack under far field coupled tensile and shear stresses, the energy release rate and a mixed fracture criterion are obtained using an energy balance approach. The additional compliance tensor induced by a single opening elliptic microcrack in a representative volume element is derived, and the effect of microcracks with random orientations is analyzed with the Taylor＇s scheme by introducing an appropriate probability density function. A micromechanical damage model for rocks and concretes is obtained and is verified with experimental results.     

On a new electromechanical switch using the reversible wavy elastic response of metallic glass ribbons

For the first time, practical applications as an alarm device and automatic filling of an aquarium using an electromechanical switch manufactured from metallic glass (MG) ribbon is proposed. The elastic response of an initial arc-shaped MG ribbon-based Fe_90.65B_3.9Cr_2.75Si_2.7 is studied and exploited. Under the applied load F, the amorphous material exhibits a reversible elastic wavy response. During the elastic deformation and multiplication of harmonic undulations, a perfect linear contact between the waves and support is established. This contact position is the same for the pair waves, and can be employed to ensure the passage of an electric current, since the ribbon is Fe-based. The reversible elastic wavy response of MG ribbon can be used as an electromechanical switch. The lifespan of the ribbon used as a switch is also considered.     

多自由度非线性动力学方程的能量校准算法

针对非线性动力学方程,通过Taylor展开和Duhamel积分,得到一个具有待定参数的逐步积分求解公式;通过数学变换,将原动力学方程转换为一个能确定待定参数的能量校准方程;最后将该参数回代入逐步积分公式,得到数值解.数值算例的结果说明了该方法的有效性,可以消除算法阻尼和抑制数值解发散,同时,在大步长的条件下也得到了非常准确而且稳定的结果,可以对系统长期性态进行仿真.     

Convective self-oscillations near an air-bubble surface in a horizontal rectangular channel

The concentration convection in an isothermal fluid near an air bubble clamped between the vertical walls of a horizontal channel with a rectangular cross-section is studied experimentally and numerically. The channel is filled with an aqueous solution of a surfactant with a nonuniform concentration. As a result of the competition between the gravitational convection in the cavity volume and the Marangoni convection near the bubble surface, an oscillation flow regime is established. This regime is observed experimentally over several hours. In the numerical experiment, the oscillations are obtained in the presence of an initial horizontal surfactant concentration gradient. Against the background of gravitational convection, short bursts of Marangoni convection with ten times greater intensity are observed. The convective flow patterns and the oscillation periods obtained experimentally and numerically are in fairly good agreement.     

SH波对浅埋弹性圆柱及裂纹的散射与地震动

采用Green函数、复变函数和多极坐标等方法研究含圆柱形弹性夹杂的弹性半空间中任意位置、任意方位有限长度裂纹对SH波的散射与地震动. 构造了含圆柱形弹性夹杂的半空间对SH波的散射波,并求解了适合本问题Green函数,即含有圆柱形弹性夹杂的半空间内(表面)任意一点承受时间谐和的出平面线源载荷作用时位移函数的基本解答. 利用裂纹``切割'方法在任意位置构造任意方位的裂纹,可以得到基体中圆柱形弹性夹杂和裂纹同时存在条件下的位移场与应力场. 通过数值算例,讨论各种参数对夹杂上方地表位移的影响.