Two-Phase Models of Formation of Cavitating Spalls in a Liquid

The dynamics of the structure of a liquid layer structure (with microbubbles of a free gas) behind a rarefaction wave front is studied numerically using the two-phase Iordansky–Kogarko–van Wijngaarden model and the frozen mass-velocity field model. An analysis of the initial stage of cavitation by the Iordansky–Kogarko–van Wijngaarden model showed that tensile stresses behind the rarefaction wave front relax quickly and the mass-velocity field in the cavitation zone turns out to be frozen. This effect is used to describe the late stage of the development of the cavitation zone. These models were combined to study the formation of cavitating spalls in a free-surface liquid under shock-wave loading.     

The effect of rotation on conical wave beams in a stratified fluid

Experiments are conducted to test extant theory on the effect of uniform rotation on the angle of conical beam wave propagation excited by a sphere vertically oscillating at frequency in a density stratified fluid. The near-constant Brunt–Väisälä frequency stratification N produced in situ in a rotating cylindrical tank exhibits no effect of residual motion for the range of Froude numbers investigated. Good agreement between experiment and theory is found over the range of angles 15°<<65° using the synthetic schlieren visualization technique. In particular, the cut-off for wave propagation at =2, below which waves do not propagate, is clearly observed.     

Effect of radial and rotational lubricant inertia on the pressure distribution in externally pressurised thrust bearings

Expressions are obtained for the pressure distribution in an externally pressurised thrust bearing for the condition when one bearing surface is rotated. The influence of centripetal acceleration and the combined effect of rotational and radial inertia terms are included in the analysis. Rotation of the bearing causes the lubricant to have a velocity component in an axial direction towards the rotating surface as it spirals radially outwards between the bearing surfaces. This results in an increase in the pumping losses and a decrease in the load capacity of the bearing. A further loss in the performance of the bearing is found when the radial inertia term, in addition to the rotational inertia term is included in the analysis.Nomenclature r, z, cylindrical co-ordinates - V _r, V , V _z velocity components in the r, and z directions respectively - U, X, W representative velocities - coefficient of viscosity - p static pressure at radius r - p mean static pressure at radius r - Q volume flow per unit time - 2h lubricant film thickness - density of the lubricant - r ₂ outside radius of bearing = D/2 - angular velocity of bearing - R dimensionless radius = r/h - P dimensionless pressure = h ³ p/Q - Re channel Reynolds number = Q/h     

An application of non-singularity BEM to plate bending, buckling and natural vibration problems

In this paper the fundamental solution of the singular governing equation of plate static bending is taken as the Green's function, which can satisfy the governing equation precisely in the plate region. Based on the principle of superposition, let the function values on the plate boundary, induced by a set of the Green's function sources (including the known sources in the plate region and the unknown sources in the fictitious region), satisfy the prescribed conditions on specially chosen boundary matching points, and the corresponding semi-analytical and semi-numerical solution can be obtained, which is free from the restraint of boundary forms and boundary conditions. The more matching points there are on the boundary, the better the accuracy of results is. Finally, in static bending problems a set of linear algebraic equations has to be computed; in buckling problems the minimum value of buckling eigenvalue equation has to be found; in natural vibration problems the eigenvalues of the frequency equation have to be calculated. Numerical examples are given and the results are compared with those by the analytical method and other methods. It can be seen that they are very close to each other.     

Spatial Simple Waves on a Shear Flow

The paper studies simple waves of the shallowwater equations describing threedimensional wave motions of a rotational liquid in a freeboundary layer. Simple wave equations are derived for the general case. The existence of unsteady or steady simple waves adjacent continuously to a given steady shear flow along a characteristic surface is proved. Exact solutions of the equations describing steady simple waves were found. These solutions can be treated as extension of Prandtl–Mayer waves for sheared flows. For shearless flows, a general solution of the system of equations describing unsteady spatial simple waves was found.     

The problems of nonlinear bending for orthotropic rectangular plate with four clamped edges

（黄家寅）（秦圣立）ＴＨＥＰＲＯＢＬＥＭＳＯＦＮＯＮＬＩＮＥＡＲＢＥＮＤＩＮＧＦＯＲＯＲＴＨＯＴＲＯＰＩＣＲＥＣＴＡＮＧＵＬＡＲＰＬＡＴＥＷＩＴＨＦＯＵＲＣＬＡＭＰＥＤＥＤＧＥＳ￥ＨｕａｎｇＪｉａｙｉｎ;ＱｉｎＳｈｅｎｇｌｉ（ＱｕｆｕＮｏｒｍａｌＵｎ...     

Unique Periodic Orbits for Delayed Positive Feedback and the Global Attractor

The delay differential equation, (t)=–x(t)+f(x(t–1)), with >0 and a real function f satisfying f(0)=0 and f>0 models a system governed by delayed positive feedback and instantaneous damping. Recently the geometric, topological, and dynamical properties of a three-dimensional compact invariant set were described in the phase space C=C([–1, 0], ) of initial data for solutions of the equation. In this paper, for a set of and f which include examples from neural network theory, we show that this three-dimensional set is the global attractor, i.e., the compact invariant set which attracts all bounded subsets of C. The proof involves, among others, results on uniqueness and absence of periodic orbits.     

The fixed point and eigenelement of sum and product about concave and convex operators

In this paper, we obtain the sufficient conditions under which there exists the fixed point of sum and product about concave and – convex operators in the positive cone of linear semi-order space, and the iterative procedure and error estimate can be given. The relation between eigenvalue and eigenelement will also be studied in this paper.     

On the linear instability of the Hall-Stewartson vortex

It is demonstrated that the Hall-Stewartson leading-edge vortex is linearly unstable to viscous perturbations of the center-mode type. Center modes are found to occur in two reigons of Reynolds-number-wave-number space, in limits in which the axial wave number is large. The appropriate center-mode equations in these neighborhoods are established, and it emerges that the two sets are identical. The single system of equations, which depends on the azimuthal wave number m and a distance parameter only, is solved numerically for various values of m and . Highly unstable modes are found for large positive , and the results are shown to be in good agreement with proposed asymptotic expansions when >1. To lowest order, unstable modes have phase surfaces that rotate with the fluid: in addition constant phase surfaces propagate upstream but the group velocity is directed downstream. The growth rate of the instability increases faster than Reynolds number to the quarter power. This, together with the finding that the length scale of the unstable modes found goes to zero as the Reynolds number tends to infinity, makes this instability an unusual one.This work was supported by the Air Force Office of Scientific Research under contract AFOSR-89-0346 monitored by Dr. L. Sakell, and by the U.S. Army Research Office at the Mathematical Sciences Institute of Cornell University.     

Sources of Complexity in Human Systems

Complex is a special attribute we can give to many kinds of systems. Although it is used often as a synonym of difficult, it has a specific epistemological meaning, which is going to be shared by the incoming science of complexity. Difficult is an object which, by means of an adequate computational power, can be deterministically or stochastically predictable. On the contrary complex is an object which can not be predictable because of logical impossibility or because its predictability would require a computational power far beyond any physical feasibility, now and forever. For complexity refers to some observing system, it is always subjective, and thus it is defined as observed irreducible complexity. Human systems are affected by several sources of complexity, belonging to three classes, in order of descending restrictivity. Systems belonging to the first class are not predictable at all, those belonging to the second class are predictable only through an infinite computational capacity, and those belonging to the third class are predictable only through a trans-computational capacity. The first class has two sources of complexity: logical complexity, directly deriving from self-reference and Gödel's incompleteness theorems, and relational complexity, resulting in a sort of indeterminacy principle occurring in social systems. The second class has three sources of complexity: gnosiological complexity, which consists of the variety of possible perceptions; semiotic complexity, which represents the infinite possible interpretations of signs and facts; and chaotic complexity, which characterizes phenomena of nonlinear dynamic systems. The third class coincides with computational complexity, which basically coincides with the mathematical concept of intractability. Artificial, natural, biological and human systems are characterized by the influence of different sources of complexity, and the latter appear to be the most complex.     

The evolution laws of dilatational spherical and cylindrical weak nonlinear shock waves in elastic non-conductors

The theory of singular surfaces yields a set of coupled evolution equations for the shock amplitude and the amplitudes of the higher order discontinuities which accompany the shock. To solve these equations, we use perturbation methods with a perturbation parameter characterising the initial shock amplitude. It is shown that for decaying shock waves, if the accompanying second order discontinuity is of order one, the straightforward perturbation procedure yields uniformly valid solutions, but if the accompanying second order discontinuity is of order , the method of multiple scales is needed in order to render the perturbation solutions uniformly valid with respect to the distance of travel. We also construct shock wave solutions from modulated simple wave solutions which are obtained with the aid ofHunter & Keller's Weakly Nonlinear Geometrical Optics method. The two approaches give exactly the same results within their common range of validity. The explicit evolution laws thus obtained enable us to see clearly how weak nonlinear curved shock waves are attenuated because of the effects of geometry and material nonlinearity, and on what length scale these effects are most pronounced. Communicated by C. C. Wang     

Investigation of the flow past wings of infinite span with a blunt leading edge in the vorticity interaction regime

An asymptotic analysis of the Navier-Stokes equations is carried out for the case of hypersonic flow past wings of infinite span with a blunt leading edge when 0, Re , and M . Analytic solutions are obtained for an inviscid shock layer and inviscid boundary layer. The results of a numerical solution of the problems of vorticity interaction at the blunt edge and on the lateral surface of the wing are presented. These solutions are compared with the solution of the equations of a thin viscous shock layer and on the basis of this comparison the boundaries of the asymptotic regions are estimated.deceasedTranslated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 120–127, November–December, 1987.     

Flow over a body with development of a steady separation zone as re → ∞

This paper discusses formulation of the total problem of flow of an incompressible liquid over a body, with formation of a closed stationary separation zone as Re . The scheme used is based on the method of matched asymptotic expansions [1]. Following [1], it is postulated that the separated zone is developed (i.e., it is not infinitely fragmented and does not vanish as Re ), and the flow inside it has a definite degree of regularity with respect to Re. With these hypotheses we can use the Prandtl-Batchelor theorem [2], which states that, in the limit as Re , a region of circulating flow becomes vortex flow of an inviscid liquid with constant vorticity . Therefore, a basis for constructing matched asymptotic expansions is the vortex-potential problem (the problem of determining a stream function , satisfying the equation = 0 in the region of translational motion and the equation = in a certain region, unknowna priori, of circulating motion). In the general case the solution of the vortex-potential problem depends on two parameters: the total pressure p_o and the vorticity in the separated zone. These parameters appear in the condition for matching the solutions of the first and second boundary-layer approximations (at the boundary of the separated zone for the end Re values) with the corresponding solutions for the inviscid flow. It is shown in the present paper that the conditions for matching the cyclic boundary layer with the external translational flow are the same additional relations which allow us to close the total problem. Thus, in using the method of matched asymptotic expansions to solve the problem of flow over a body with closed stationary separation zones one must simultaneously consider no less than two approximations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 28–37, March–April, 1978.The authors thank G. Yu. Stepanov for discussion of the paper and valuable comments.     

Neck propagation as a shock wave

Neck propagation in the stretching of elastic solid filaments having a yield point was analyzed using the space one-dimensional thin filament governing equations developed previously by the authors and other researchers. Constitutive model for the filament was assumed to be expressible as engineering tensile stress(X) (tensile force) given as a function of elongational strain with the(X) curve having a yield point maxima followed by a minima and a breaking point greater than the yield point maxima. Also incorporated into the model is the hysteresis of irreversible plastic deformation. When inertia is taken into consideration, the thin filament equations were found to reduce to the nonlinear wave equation ² (X)/ ² =C ₁ ² X/ ² where is Lagrangean space coordinate, is time, andC ₁ is inertia coefficient. The above nonlinear wave equation yields a solutionX(, ) having a stepwise discontinuity inX which propagates along the axis. The zero speed limit of the step wave solution was found to describe the above neck propagation occurring in solid filaments. Furthermore, it was recognized that the nonlinear wave equation was known for many years to also govern the plastic shock wave which propagates axially within a metal rod subjected to a very strong impact on its end. The one-dimensional atmospheric shock wave also was known to be governed by the nonlinear wave equation upon making certain simplifying assumptions. The above and other evidences lead to the conclusion that neck propagation occurring in the extension of solid filament obeying the above(X) function can be formally described as a shock wave.     

Construction of exact numerical solutions of the stationary traveling wave type for viscous thin films

An effective numerical procedure, based on the Galerkin method, for finding solutions of the stationary traveling wave type in the complete formulation is proposed for the case of viscous liquid films. Examples of a viscous film flowing freely down a vertical surface have been calculated. The calculations have been made for various values of the dimensionless surface tension , including =0. The method makes it possible to predict a number of bifurcations that occur as decreases. The existence of numerous families of stationary traveling waves when 1 was demonstrated in [6]. The present study shows that as 1 all but one of these families of wave solutions disappear. The shape of the periodic and solitary waves and the pressure distribution in the film are found for various . When =0 and the wave number is fairly small, the periodic solution has a singularity, as predicted in [14]: at the crest of the wave a corner point appears; the angle between the tangents at this point =140–150. The method proposed can be used to calculate other wavy film flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 94–100, May–June, 1990.     

Shock Waves in a Fluid with Bubbles Containing Evaporating Drops of a Liquefied Gas

The propagation and reflection of one-dimensional plane unsteady waves and pulses in a mixture of a fluid with two-phase bubbles containing evaporating drops is investigated. A significant effect of unsteady evaporation of the drops in the zone ahead of the shock wave on the wave propagation is demonstrated. The evaporation of the drops results in a pressure increase ahead of the wave and the shock wave as it were climbs to increasing pressure level. In contrast to bubbly fluids with single-phase bubbles, in a fluid with two-phase bubbles, at a fixed phase volume fraction, a decrease in bubble size results in an increase rather than a decrease of the oscillation amplitude. The wave reflection from a solid wall is essentially nonlinear and the maximum pressure attained at the wall is several times greater than the incident-wave intensity.     

Asymptotic solution to the problem of hypersonic flow over bodies in the neighborhood of the separation point of a thin shock layer

A study is made of the problem of hypersonic flow of an inviscid perfect gas over a convex body with continuously varying curvature. The solution is sought in the framework of the asymptotic theory of a strongly compressed gas [1–4] in the limit M when the specific heat ratio tends to 1. Under these assumptions, the disturbed flow is situated in a thin shock layer between the body and the shock wave. At the point where the pressure found by the Newton-Buseman formula vanishes there is separation of the flow and formation of a free layer next to the shock wave [1–4]. The singularity of the asymptotic expansions with respect to the parameter ₁ = ( –1)/( + 1) associated with separation of the strongly compressed layer has been investigated previously by various methods [3–9]. Local solutions to the problem valid in the neighborhood of the singularity have been obtained for some simple bodies [3–7]. Other solutions [7, 9] eliminate the singularity but do not give the transition solution entirely. In the present paper, an asymptotic solution describing the transition from the attached to the free layer is constructed for a fairly large class of flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 99–105, January–February, 1982.     

Hypersonic flow past blunt-edged delta wings

A method is proposed for calculating hypersonic ideal-gas flow past blunt-edged delta wings with aspect ratios = 100–200. Systematic wing flow calculations are carried out on the intervals 6 M 20, 0 20, 60 80; the results are analyzed in terms of hypersonic similarity parameters.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 175–179, September–October, 1990.     

Asymptotic Modeling of Nonlinear Wave Processes in Shock‐Loaded Elastoplastic Materials

Nonlinear wave processes in shockloaded elastoplastic materials are modeled. A comparison of the results obtained with experimental data and numerical solutions of exact systems of dynamic equations shows that the model equations proposed qualitatively describe the stressdistribution evolution in both the elasticflow and plasticflow regions and can be used to solve one and twodimensional problems of pulsed deformation and fracture of elastoplastic media.