Some characteristics of low-speed streaks under sheared air-water interfaces

The characteristics of low-speed fluid streaks occurring under sheared air-water interfaces were examined by means of hydrogen bubble visualization technique. A critical shear condition under which the streaky structure first appears was determined to beu _τ≈0.19 cm/s. The mean spanwise streak spacing increases with distance from the water surface owing to merging and bursting processes, and a linear relationship describing variation of non-dimensional spacing versusy ⁺ was found essentially independent of shear stress on the interface. Values of , however, are remarkably smaller than their counterparts in the near-wall region of turbulent boundary layers. Though low-speed streaks occur randomly in time and space, the streak spacing exhibits a lognormal probability distribution behavior. A tentative explanation concerning the formation of streaky structure is suggested, and the fact that takes rather smaller values than that in wall turbulence is briefly discussed. The project supported by the National Natural Science Foundation of China (19672070)     

Dimensional analysis of pore scale and field scale immiscible displacement

A basic re-examination of the traditional dimensional analysis of microscopic and macroscopic multiphase flow equations in porous media is presented. We introduce a macroscopic capillary number which differs from the usual microscopic capillary number Ca in that it depends on length scale, type of porous medium and saturation history. The macroscopic capillary number is defined as the ratio between the macroscopic viscous pressure drop and the macroscopic capillary pressure. can be related to the microscopic capillary number Ca and the LeverettJ-function. Previous dimensional analyses contain a tacit assumption which amounts to setting = 1. This fact has impeded quantitative upscaling in the past. Our definition for , however, allows for the first time a consistent comparison between macroscopic flow experiments on different length scales. Illustrative sample calculations are presented which show that the breakpoint in capillary desaturation curves for different porous media appears to occur at 1. The length scale related difference between the macroscopic capillary number for core floods and reservoir floods provides a possible explanation for the systematic difference between residual oil saturations measured in field floods as compared to laboratory experiment.     

Pseudo-molecular dynamics study of grain boundary segregation

The segregation of bismuth atoms on the [101] tilt copper grain boundaries Σ3 ( ) 70.53°, Σ33 ( ) 58.99°, Σ11 ( ) 50.48° and Σ9 ( ) 38.94° has been studied by pseudo-molecular dynamics using the empiricalN-body potentials. The relationship between bismuth segregation and grain boundary structure has been discussed in detail. The subject supported by the Chinese Academy of Sciences and National Natural Science Foundation of China     

On the dynamic behaviour of an elastically supported beam of infinite length, loaded by a concentrated force

Summary An elastically supported beam of infinite length, initially at rest, carries a variable concentrated force at a prescribed point A. General expressions are given for the deflection and the bending moment at A (6.3 and 6.4). Three special cases are considered; the first one is defined by =0 for and =K=const. for ; the second one by =0 for 0 > > , given function of for 0 ; the third one applies to problems in which, during the period of impact, itself is an unknown. The results given here may be of use in those railway-engineering problems in which a rail can be considered as a beam of infinite length, and in which the supporting ground has the required properties.     

Optimal cupolas of uniform strength: Spherical M-shells and axisymmetric T-shells

Summary The first part of this paper is concerned with the optimal design of spherical cupolas obeying the von Mises yield condition. Five different load combinations, which all include selfweight, are investigated. The second part of the paper deals with the optimal quadratic meridional shape of cupolas obeying the Tresca yield condition, considering selfweight plus the weight of a non-carrying uniform cover. It is established that at long spans some non-spherical Tresca cupolas are much more economical than spherical ones.

---

Optimale Kuppeln gleicher Festigkeit: Kugelschalen und axialsymmetrische Schalen

Übersicht Im ersten Teil dieser Arbeit wird der optimale Entwurf sphärischer Kuppeln behandelt, wobei die von Misessche Fließbewegung zugrunde gelegt wird. Fünf verschiedene Lastkombinationen werden untersucht. Der zweite Teil befaßt sich mit der optimalen quadratischen Form des Meridians von Kuppeln, die der Fließbedingung von Tresca folgen.

---

List of Symbols a_k, b_k, c_k, A_k, B_k, C_k coefficients used in series solutions - A, B constants in the nondimensional equation of the meridional curve - normal component of the load per unit area of the middle surface - meridional and circumferential forces per unit width - radial pressure per unit area of the middle surface, - skin weight per unit area of the middle surface, - vertical external load per unit horizontal area, - base radius, - R radius of convergence - s - cupola thickness, - u, w subsidiary functions for quadratic cupolas - vertical component of the load per unit area of middle surface - resultant vertical force on a cupola segment - structural weight of cupola, - combined weight of cupola and skin, - distance from the axis of rotation, - vertical distance from the shell apex, - z auxiliary variable in series solutions - specific weight of structural material of cupola - radius of the middle surface, - uniaxial yield stress - meridional stress, - circumferential stress, - _a, _b, _c, _d, _e subsidiary variables used in evaluating the meridional stress - auxiliary function used in series solutions This paper constitutes the third part of a study of shell optimization which was initiated and planned by the late Prof. W. Prager     

Molecular dynamics study of deformation and fracture of Bi-segregated copper bicrystals

The microprocesses of deformation and fracture of Bi-segregated copper bicrystals Σ33 ( ) 58.99°, Σ11 ( ) 50.48° and Σ9 ( ) 38.94° have been simulated by molecular dynamics in order to study the relationship between the grain boundary embrittlement (GBE) and grain boundary (GB) structure. It is shown that GBE is related to the segregated concentration and distribution of Bi atoms, while Bi segregation is related to the GB structure. Due to their different structures, the bicrystals Σ33, Σ11 and Σ9 show an increasing propensity for Bi segregated concentration. So under the action of external force, Σ33, Σ11 and Σ9 show transgranular ductile, intergranular tearing and intergranular brittle fracture modes, respectively. The subject supported by the Chinese Academy of Sciences and National Natural Science Foundation of China     

Generalized diagonalization of matrices over quaternion field

Ｉｎｒｅｃｅｎｔｙｅａｒｓ,ａｐｐｌｉｃａｔｉｏｎｓｏｆｑｕａｔｅｒｎｉｏｎｍａｔｒｉｃｅｓａｒｅｂｅｃｏｍｉｎｇｍｏｒｅａｎｄｍｏｒｅｉｍｐｏｒｔａｎｔａｎｄｅｘｔｅｎｓｉｖｅｉｎｒｉｇｉｄｍｅｃｈａｎｉｃｓ,ｑｕａｎｔｕｍｍｅｃｈａｎｉｃｓ,ｃｏｎｔｒｏｌｔｈｅｏｒｙａｎｄｈｅｌｉｃａｌｔｅｃｈｎｏｌｏｇｙ［１～３］．Ｗｉｔｈｔｈｅｒａｐｉｄｄｅｖｅｌｏｐｍｅｎｔｏｆｔｈｅａｂｏｖｅｄｉｓｃｉｐｌｉｎｅｓ,ｉｔｉｓｇｅｔｔｉｎｇｍｏｒｅａｎｄｍｏｒｅｎｅｃｅｓｓａｒｙｆｏｒｕｓｔｏｆｕｒｔｈ…     

A study of stationary crack-tip deformation fields in thin sheets by computer vision

The in-plane deformation fields near a stationary crack tip for thin, single edge-notched (SEN) specimens, made from Plexiglas, 3003 aluminum alloy and 304 stainless steel, have been successfully obtained by using computer vision. Results from the study indicate that (a) in-plane deformations ranging from elastic to fully plastic can be obtained accurately by the method, (b) for U, and , the size of the HRR dominant zone is much smaller than forV and , respectively. Since these results are in agreement with recent analytical work, suggesting that higher order terms will be needed to accurately predict trends in the data, it is clear that the region where the first term in the asymptotic solution is dominant is dependent on the component of the deformation field being studied, (c) the HRR solution can be used to quantity only in regions where theplastic strains strongly dominate the elastic strain components (i.e., when ); forV, the HRR zone appears to extend somewhat beyond this region, (d) the displacement componentU does not have the HRR singularity anywhere within the measurement region for either 3003 aluminum or 304 SS. However, the displacement componentV agrees with the HRR slope up to the plastic-zone boundary in 3003 aluminum ( ) and over most of the region where measurements were obtained ( ) in 304 SS and (e) the effects of end conditions must be included in any finite-element model of typical SEN specimen geometries to accurately calculate theJ integral and the crack-tip fields.Paper was presented at the 1992 SEM Spring Conference on Experimental Mechanics held in Las Vegas, NV on June 8–11.     

Transient response of a finite crack in an orthotropic strip under the in-plane shear impact

The transient response of a central crack in an orthotropic strip under the in-plane shear impact loading is studied by using the dual integral equation method proposed by Copson and Sih. The general formula for the shear stress intensity factor near the crack tip is derived. Numerical results of with in various cases are obtained by solving the second kind Fredholm integral equation and by performing the inverse Laplace transform.     

Druckabhängigkeit der Viskosität einiger Polyolefinschmelzen

Zusammenfassung Es wird die Druckabhängigkeit des nicht -Newtonschen Fließverhältens von Polyolefinschmelzen (Hochdruck-, Niederdruck-,Phillips-Polyäthylen und Polypropylen) experimentell untersucht und gefunden, daß der durch Gl. [1] definierte Druckkoeffizient mit zunehmender Deformationsgeschwindigkeit kleiner wird und dabei die (im einzelnen in der Tabelle 1 angeführten) Werte annimmt. Der Druckkoeffizient der Polyolefinschmelzen ist ebenso wie für vieleNewtonsche Flüssigkeiten bis 2000 kp cm^–2 unabhängig vom Druck, er wird mit zunehmender Temperatur kleiner und nimmt mit zunehmender Verzweigung zu. Die Meßergebnisse werden mit Hilfe eines Aufweitungsvolumens interpretiert. Es wird gezeigt, daß eine Deutung des Fließverhaltens von Polyäthylen durch das freie Volumen allein nicht möglich ist.

---

Summary The influence of pressure of the non-Newtonian flow behaviour of polyolefin melts (Low- and High density Polyethylene,Phillips-Polyethylene and Polypropylene) was investigated. The results are: the coefficient of pressure as defined by eq. [1], decreases with increasing shear rate and reaches the values given in table 1 . The pressure coefficient of polyolefin melts does not depend on pressure up to 2000 kp cm^–2. As observed with manyNewtonian fluids, decreases with increasing temperature and increases with the degree of branching. The experimental results are explained by means of a so called volume of expansion. It has been shown that it is impossible to explain the flow behaviour of polyethylene exclusively with the free volume.

---

---

Für die Diskussion und Förderung dieser Arbeit danke ich Herrn Professor Dr.K.-H. Hellwege und Herrn Dr.W. Knappe.     

Theoretical mean wave resistance of precursor solition generation in two-layer flow

A new concept of pseudo mean wave resistance is introduced to find theoretical mean wave resistances of the precursor soliton generation in two-layer flow over a localized topography at near-resonance in this paper. The pseudo mean wave resistance of the precursor soliton generation of two-layer flow is determined in terms of the AfKdV equation. From the theoretical results it is shown that the theoretical mean wave resistance is equal to the pseudo mean wave resistance times 1/m ₁, wherem ₁ is the coefficient of the fKdV equation. From the regional distribution of the energy of the precursor soliton generation at the resonant points, it is shown that ratios of the theoretical mean wave resistance and regional mean energy to the total mean energy are invariant constants, i.e. : (−1/2):1:1, in which and are the mean energy of the generating regions of the precursor solitons, of the depression and of the trailing, wavetrain at the resonant points respectively, and <D> are the total energy of the system and the theoretical mean wave resistance at the resonant points. A prediction of the theoretical mean wave resistance flow over the semicircular topography is carried out in terms of the theoretical results of the present paper. The comparison shows that the theoretical mean wave resistance is in good agreement with the numerical calculation. The project supported by the Foundation of The State Education Commission “The Dynamics of Upper Ocean” and the Grants of The Physical Oceanography Laboratory of Ocean University of Qingdao     

Plane stress dynamic fields near a propagating crack-tip in a power-law material

Based on the plastic-dynamic equations, the asymptotic behaviour of the near-tip fields for a plane stress tensile crack propagating in a power-law material has been studied in this paper. It is shown that the stress and strain singularities are, respectively, of the order and , whereA is a constant which is related to the size of plastic region,r is the distance to the crack tip,n is the power-law exponent. Projects sponsored by the National Science Foundation.     

Behaviour of the Reynolds stress on rough walls

Quadrant analysis is used to study the contributions, associated to the four quadrants of the (u,v) plane, to the production of the turbulent shear stress on rough walls. The measurements are described for a fully developed turbulent flow between two rough plates with varying the parameterλ _z(span/height ratio of roughness elements). The application of this technique indicates that the (Q ₂) events (ejections) and (Q ₄) events (sweeps) cause an intense production of the Reynolds stress —ρ as compared with the (Q ₁) and (Q ₃) quadrants. The (Q ₂) contribution to the Reynolds stress depends on the geometry factorλ _z. Variation of the parameterλ _zaffects the distributions of and forH⩽3, whereH is a particular threshold. Comparison with boundary layer flow shows that the region 0.2⩽y/h⩽0.7 is characterized by a Reynolds stress production, independent of the flow nature. The third moment of the longitudinal velocity fluctuations is found to be sensitive to the surface roughness.     

Interval stochastic matrices: A combinatorial lemma and the computation of invariant measures of dynamical systems

The concept of an interval stochastic matrix is introduced. We prove a combinatorial theorem which describes the network flow associated with an interval matrix. The semi-invariant vectors of are characterized in terms of eigenvectors with unit eigenvalue of stochastic matrices . These results are then applied to the approximation and machine computation of invariant measures of dynamical systems.Funded under Australian Research Council Grant A 8913 2609.     

About a condition for blow up of solutions of cauchy problem for a wave equation

ＣｏｎｓｉｄｅｒｔｈｅＣａｕｃｈｙｐｒｏｂｌｅｍｆｏｒｔｈｅｗａｖｅｅｑｕａｔｉｏｎｉｎＲＮ×Ｒ＋（Ｎ≥２）：２ｕ（ｘ,ｔ）ｔ２－ｘｉａｉｊ（ｘ）ｘｊｕ＝｜ｕ｜ｐ－１·ｕ　　（（ｘ,ｔ）∈ＲＮ×（０,Ｔ））,ｕ（ｘ,０）＝ｇ（ｘ）　（ｘ∈ＲＮ）,ｕｔ（ｘ,０）＝ｈ（ｘ）　（ｘ∈ＲＮ）,（１）ｗｈｅｒｅｕ（ｘ,ｔ）ｉｓｎｏｎｔｒｉｖｉａｌｓｏｌｕｔｉｏｎｗｉｔｈｆｉｎｉｔｅｓｐｅｅｄｏｆｐｒｏｐａｇａｔｉｏｎａｎｄｉｓｓｕｐｐｏｒｔｅｄｏｎａｆｏｒｗａｒｄｃｏｎｅ（ｘ,ｔ）·ｔ≥０,｜…     

Determining parameters in nonequilibrium flow of a gas with homogeneous condensation

Analytic expressions are obtained for the parameters which determine the number of particles (clusters) and their mean radius and the supercooling at the starting point of active homogeneous condensation (the Wilson point). With a given law p=f(T), such parameters are the temperature T_s at the dew point and the mean cooling rate = dT/dt of the gas from the dew point to the Wilson point. In this context, with T_s = const and = const, the processes of nonequilibrium homogeneous condensation in the P-T diagram corresponding to the various initial states after the Wilson point may be described approximately by a single curve. So the parameters T_s and are determining parameters for the whole process of nonequilibrium homogeneous condensation. Equations are given which enable us to estimate approximately the cooling rate for flows in nozzles and jets. Computation results are given for steady flow in a Laval nozzle and unsteady spherically symmetric expansion of a gaseous sphere into a vacuum; these results demonstrate the noted regularities. Despite the great amount of literature, both periodical and monograph, on the nonequilibrium condensation of vapors, the important question of the determining parameters for the process of a homogeneous nonequilibrium condensation has not been studied to any great extent [1].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 80–87, May–June, 1985.     

The flow regime in the turbulent near wake of a hemisphere

The work presented is a wind tunnel study of the near wake region behind a hemisphere immersed in three different turbulent boundary layers. In particular, the effect of different boundary layer profiles on the generation and distribution of near wake vorticity and on the mean recirculation region is examined. Visualization of the flow around a hemisphere has been undertaken, using models in a water channel, in order to obtain qualitative information concerning the wake structure.List of symbols C _p pressure coefficient, - D diameter of hemisphere - n vortex shedding frequency - p pressure on model surface - p ₀ static pressure - Re Reynolds number, - St Strouhal number, - U, V, W local mean velocity components - mean freestream velocity inX direction - U ^* shear velocity, - u, v, w velocity fluctuations inX, Y andZ directions - X Cartesian coordinate in longitudinal direction - Y Cartesian coordinate in lateral direction - Z Cartesian coordinate in direction perpendicular to the wall - it* boundary layer displacement thickness, - diameter of model surface roughness - elevation angleI - O boundary layer momentum thickness, - _w wall shearing stress - dynamic viscosity of fluid - density of fluid - streamfunction - _x longitudinal component of vorticity, - _y lateral component of vorticity, - _z vertical component of vorticity, This paper was presented at the Ninth symposium on turbulence, University of Missouri-Rolla, October 1–3, 1984     

The double velocity correlation function of homogeneous turbulence with constant mean velocity gradient

In this article, as the velocity gradient is taken as a constant value, we obtain the solutions of the equation of fluctuation velocity after Fourier transform. When the mean velocity gradient is small, they represent the picture of eddies, of which the homogeneous turbulence (both isotropic and nonisotropic) of the final period is composed. By using the eddies of these types at different times, we may compose the steady turbulent field with the constant velocity gradient and this field may represent the turbulent field in the central part of the channel flow or pipe flow approximately. Then we may obtain the double velocity correlation function of this turbulent field, which involves both longitudinal correlation coefficient and the transversal correlation coefficient . We compare these theoretical coefficients with the experimental data of these coefficients at initial period and final period of isotropic homogeneous turbulence. And then we obtain the relationship between the turbulent double velocity correlation coefficient and the mean velocity gradient. Finally, we get the expressions of the Reynolds stress and the eddy viscosity coefficient.     

Effect of a Short Roughness Strip on a Turbulent Boundary Layer

The response of a turbulent boundary layer to a short roughness strip is investigated using laser Doppler velocimetry (LDV) and laser induced fluorescence (LIF). Skin friction coefficients are inferred from accurate near-wall measurements. There is an undershoot in , where is the undisturbed smooth wall skin friction coefficient, immediately after the strip. Downstream of the strip, overshoots before relaxing back to unity in an oscillatory manner. The roughness strip has a major effect on the turbulent stresses ; these quantities increase, relative to the undisturbed smooth wall, in the region between the two internal layers originating at the upstream and downstream edges of the strip. The increase in the ratio suggests a decrease in near-wall anisotropy. From the flow visualizations, it is inferred that streamwise vortical structures are weakened immediately downstream of the strip. Consistently, streamwise length scales are also reduced; direct support for this is provided by measured two-point velocity correlations.     

Grundzüge einer Stabilitätstheorie für elastische Systeme unter nichtkonservativer Belastung

Zusammenfassung Um den Einfluß von nichtkonservativen Kräften zu erfassen, muß die Stabilitätstheorie nichtkonservativer elastischer Systeme auf dem Prinzip der virtuellen Verrückung aufgebaut werden, das von genügender Allgemeinheit ist. Als mathematische Methode steht nach Klärung von Konvergenzfragen das Galerkinsche Verfahren zur Verfügung, das dem ursprünglichen Problem der Stabilität elastischer Systeme das Problem der Stabilität der Bewegung eines mechanischen Systems von unendlich vielen Freiheitsgraden als gleichwertig zuordnet und die Algebraisierung des Stabilitätsproblems erlaubt. Damit wird die Ljapunowsche Theorie anwendbar. Als Wichtigstes ist die Berechtigung der Linearisierung zu untersuchen, wozu der Einfluß der Dämpfung festzustellen ist. Es ergibt sich, daß bei Berücksichtigung von Dämpfung die Stabilität stets asymptotisch ist, so daß die Linesarisierung gerechtfertigt wird. Darüber hinaus zeigt es sich, daß die Struktur der Matrizen und von entscheidender Bedeutung für die Zuständigkeit eines der in Frage kommenden Stabilitätskriterien ist: Ist symmetrisch und eine beliebige positive Diagonalmatrix bzw. schiefsymmetrisch, oder symmetrisierbar bzw. unsymmetrisch und noch gleichzeitig (die Dämpfungskoeffizienten alle gleich), so wirkt in diesen Fällen die Dämpfung günstig, und die Stabilitätskriterien des ungedämpften Vorganges, das statische bzw. das kinetische, bleiben daher hinreichend. Man darf dann auch ohne Dämpfung rechnen. Ist dagegen unsymmetrisch und eine positive Diagonalmatrix mit verschiedenen Elementen, so darf in solch einem Fall nicht ohne Dämpfung gerechnet werden, weil gegenüber dem ungedämpften Vorgang ein Sprung in den Stabilitätsbedingungen auftritt, durch welchen das Ergebnis der Rechnung ohne Dämpfung in Frage gestellt wird.