Deterministic Model of Homophase and Heterophase Fluctuations in a Liquid–Vapor System

Within the framework of the assembly-type catastrophe model, a nonlinear dynamic equation (DE) homogeneous in the parameter t with an aftereffect is constructed, in which t characterizes the deviation of the reduced density of a thin surface layer on the liquid-vapor interface from the mean density of the vapor-liquid system. This equation is used to treat a second-order nonlinear DE with a variable damping coefficient for a vapor-liquid system excited by periodic impacts (acts of evaporation and condensation of molecules). This DE is integrated over a finite time interval to find a two-dimensional mapping whose numerical solution describes the chaotic dynamics of the density in time, including homophase and heterophase fluctuations. For this system, the bifurcation diagrams are constructed and the Lyapunov exponents are found.     

Laminar flow of a viscous incompressible liquid over the surface of solids of revolution

The laminar flow of a viscous incompressible liquid over the surface of stationary solids of revolution is examined for the case of circular inflow of the stream.Notation v₁, v velocity projections - density of the liquid - kinematic viscosity coefficient of liquid - R radius of curvature of surface of the body in the flow - thickness of boundary layer - u₀ rate of inflow of stream - r current radius - Q total flow rate - r₀ radius of contact line (circle) - ₀ thickness of layer at contact line - h thickness of layer - r_c radius of stream - r_min radius at which layer thickness is minimal - h_min minimum thickness of layer - 2 angle at apex of conical packing - current angle - , t, parameters of the ellipse - r_s conjugate line radius of plane circular torus Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 39, No. 1, pp. 57–63, July, 1980.     

Optimal design of a thick-walled pipeline cross-section against creep rupture

Summary Cylinder under combined loadings (pressure, bending, axial force) is subject to non-linear creep described by Norton-Odqvist creep law. In view of bending a circularly-symmetric cross-section is no longer optimal in this case. Hence we optimize the shape of the cross-section; minimal area being the design objective under the constraint of creep rupture. Kachanov-Sdobyrev hypothesis of brittle creep rupture is applied. The solution is based on the perturbation method (expansions into double series of small parameters), adjusted to optimization problems.Notation A cross-sectional area - C, , creep rupture constants - K, n, _C , _C creep constants - F dimensionless creep modulus - M bending moment - N axial force - a(),b() internal and external radii of the cross-section - j creep modulus - p internal pressure - r, ,z cylindrical coordinates - s _r ,s ,s _z ,t _r dimensionless stresses - t _R time to rupture - stress function - , () dimensionless internal and external radii - _e effective strain rate - _kl strain rates - rate of curvature - rate of elongation of the central axis - dimensionless radius - _e effective stress - _I maximal principal stress - _S Sdobyrev's reduced stress - _r , , _z , _r components of the stress tensor - measure of material continuity - measure of deterioration With 7 Figures     

Breakup of an anomolously viscous liquid film in a centrifugal force field

An equation is obtained for the breakup radius with consideration of tipping moments and Laplacian pressure forces acting on the liquid ridge at the critical point.Notation K, n rhenological constants - density - surface tension - r current cup radius - R maximum cup radius - r_c critical radius for film breakup - ¯r=¯r=r/R dimensionless current radius - ¯r_c=r_c/R dimensionless critical radius - ₀, _c actual and critical film thicknesses - current thickness - R_r ridge radius - h₀ ridge height - h current ridge height - ₀ limiting wetting angle - current angle of tangent to ridge surface - angle between axis of rotation and tangent to cup surface - angular velocity of rotation - q volume liquid flow rate - v₁ and v meridional and tangential velocities - =4vv _lm/r,=4v_m/r dimensionless velocities - M moments of surface and centrifugal forces - M_v moment from velocity head - p_r pressure within ridge - P_vm pressure from velocity head - p_m, p_pm pressures from centrifugal force components tangent and normal to cup surface - deviation range of breakup radius from calculated value - ¯r_max, ¯r_min limiting deviations of breakup radius - _c angle of tangent to curve _c/°₀=f(¯r) at critical point - t random oscillation of ratio _c/_c Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 39, No. 1, pp. 51–56, July, 1980.     

Determination of the effective thermal conductivity of a plate in nonsteady experiment

The article describes a method of finding the effective thermal conductivity of multilayered structures at the nonsteady stage of the experiment and an experimental installation with automatic systems for setting the experimental regime and processing the experimental data.Notation thickness of the layer of heat insulation - x coordinate - t(x, ) temperature at the point of the material with the coordinate x at the instant - thermal conductivity - _ef effective thermal conductivity - c heat capacity - density of the material - C_d heat capacity of the disk per unit area of the base - m cooling rate of the disk - P_c thermal contact resistance - P_r thermal resistance of the investigated layer Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 55, No. 4, pp. 616–620, October, 1988.     

A superconducting vibrating reed applied to flux-line pinning. I. Theory

A theoretical treatment is given of a superconducting reed clamped at one end and performing flexural vibrations in a homogeneous longitudinal magnetic fieldB _a. When the flux lines are rigidly pinned the reed behaves like an ideal diamagnet whose bending distorts the external field. This generates a magnetic restoring force (line tension) B _a ² which is independent of the reed thicknessd, whereas the mechanical restoring force (stiffness) is d ³. Therefore, the resonance frequency /2 of a thin superconducting reed increases drastically when a fieldB _a is applied, or for a givenB _a, when the reed is cooled below its critical temperatureT _c. With decreasing pinning strength (characterized by Labusch's parameter ) the resonance frequency decreases, ²²– _pin ² where _pin ^{2 –1}, and an attenuation _v ^–2 occurs due to the viscous motion of flux lines. For larger vibration amplitudes an additional, amplitude-dependent damping _h ^–3 occurs due to the hysteretic losses caused by elastic instabilities during flux motion.On leave from Centro Atómico, Bariloche, Argentina.     

Thermal conductivity of hydrocarbons in the naphthene group under high pressures

The thermal conductivity of hydrocarbons in the naphthene group has been experimentally determined. An equation is now proposed for calculating the thermal conductivity over the given temperature and pressure ranges.Notation thermal conductivity - ₂₀ and ₃₀ values of the thermal conductivity at 20 and 30°C, respectively - t₀,P₀ thermal conductivity at t₀, p₀ - _t p thermal conductivity at temperature t and under pressure P - change in thermal conductivity - P pressure - P_melt melting pressure - P₀ atmospheric pressure - t₀ 20°C temperature - T, t temperature - T_cr critical temperature - temperature coefficient of thermal conductivity - ₂₀ temperature coefficient of density - density - ₂₀ density at 20°C - _cr critical density - M molar mass - =T/T_cr referred temperature - v specific volume - v₀ specific volume at 20°C - v change in specific volume - _{3 0} a coefficient - B (t) a function of the temperature - S a quadratic functional - W_i, weight of the i-th experimental point - _i error of the i-th experimental value of thermal conductivity - B _{y, =0.6} value of B (t) at T = 0.6T_cr - B = B (t)/B_{, =0.6} referred value of coefficient B (t) Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 41, No. 3, pp. 491–499, September, 1981.     

Response of a spinning liquid column to pitch excitation

Summary The response of a solidly rotating liquid bridge consisting of inviscid liquid is determined for pitch excitation about its undisturbed center of mass. Free liquid surface displacement and velocity distribution has been determined in the elliptic (>2₀) and hyperbolic (<2₀) excitation frequency range.List of symbols a radius of liquid column - h length of column - I ₁ modified Besselfunction of first kind and first order - J ₁ Besselfunction of first kind and first order - r, ,z cylindrical coordinates - t time - u, v, w velocity distribution in radial-, circumferential-and axial direction resp. - mass density of liquid - free surface displacement - velocity potential - ₀ rotational excitation angle - ₀ velocity of spin - forcing frequency - _1n natural frequency - surface tension - acceleration potential - for elliptic range >2₀ - for hyperbolic range >2₀     

Effectiveness of mixing coaxial flows swirled in opposite directions

The mixing of coaxial turbulent flows swirled in opposite directions is experimentally studied. The effectiveness of this mixing is compared with mixing after an agitating grid.Notation z, r, cylindrical coordinate system - r₁, r₂ inside and outside radius of annular channel - H=r₂–r₁ radial gap in annular channel - y=(r–r₁)/H dimensionless radial coordinate - V(v_z, v_r, v) mean velocity vector - v, v_z pulsative components of velocity in the direction of the mean velocity vector and in the axial direction - P_*, P total and static pressure - , loss coefficients Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 41, No. 3, pp. 407–413, September, 1981.     

A finite volume technique to simulate the flow of a viscoelastic fluid

Hydrodynamically developing flow of Oldroyd B fluid in the planar die entrance region has been investigated numerically using SIMPLER algorithm in a non-uniform staggered grid system. It has been shown that for constant values of the Reynolds number, the entrance length increases as the Weissenberg number increases. For small Reynolds number flows the center line velocity distribution exhibit overshoot near the inlet, which seems to be related to the occurrence of numerical breakdown at small values of the limiting Weissenberg number than those for large Reynolds number flows. The distributions of the first normal stress difference display clearly the development of the flow characteristics from extensional flow to shear flow.List of symbols D rate of strain tensor - L slit halfheight - P pressure, indeterminate part of the Cauchy stress tensor - R the Reynolds number - t time - U average velocity in the slit - u velocity vector - u,v velocity components - W the Weissenberg number based on the difference between stress relaxation time and retardation time - W ₁ the Weissenberg number based on stress relaxation time - x,y rectangular Cartesian coordinates - ratio of retardation time to stress relaxation time - zero-shear-rate viscosity, ₁ + ₂ - ₁ non-Newtonian contribution to - ₂ Newtonian contribution to - ₁ stress relaxation time - ₂ retardation time - density - (, , ) xx, yy and xy components of ₁, respectively - determinate part of the Cauchy stress tensor - ₁ non-Newtonian contribution to - ₂ Newtonian contribution to      

State space formulation to viscoelastic fluid flow of magnetohydrodynamic free convection through a porous medium

Summary In this work we formulate the state space approach for one-dimensional problems of viscoelastic magnetohydrodynamic unsteady free convection flow through a porous medium past an infinite vertical plate. Laplace transform techniques are used. The resulting formulation is applied to a thermal shock problem and to a problem for the flow between two parallel fixed plates both without heat sources. Also a problem with a distribution of heat sources is considered. A numerical method is employed for the inversion of the Laplace transforms. Numerical results are given and illustrated graphically for the problem considered.Notation C specific heat at constant pressure - g acceleration due to gravity - density - time - u velocity component parallel to the plate - H _x induced magnetic field - x, y coordinates system - T temperature distribution - T _o temperature of the plate - T temperature of the fluid away from the plate - ₀ limiting viscosity at small rates to shear - v _o ^* / - v _m magnetic diffusivity - Alfven velocity - ^* coefficient of volume expansion - thermal conductivity - ^* thermal diffusivity - G Grashof number - Pr Prandtl number - L some characteristic length - k _o the elastic constant - K permeability of the porous medium     

Evolution of transient thermal processes in layer counterflow apparatuses and heat exchangers

Results are given of an analytic investigation of transient processes inside counterflow apparatuses and heat exchangers with temperature disturbance in one of the heat carriers at the entry to the apparatus.Notation =(t–t₀)/(T₀–t₀),=(T–t₀)/(T₀ s-t₀) relative temperatures - t, T temperatures of material and gas respectively - t₀, T₀ same for the initial state - Z=[ _Vm₁/c(1–w/w_g)] [–(y₀–y)/w_g] dimensionless time - m₁=1/(1+Bi/) solidity coefficient - B₁=( _FR/) Biot number - _{F V} heat-exchange coefficients referred to 1 m² surface and 1 m³ layer - R depth of heat penetration in a portion - portion heat conductivity coefficient - shape coefficient (=0 for a plate,=1 for a cylinder,=2 for a sphere) - c, C_g heat capacities of material and gas respectively - , _g volumetric masses - w, W_g flow velocities of material and gas - y distance from the point of entry to the heating heat carrier - y₀ heat-exchanger length - Y= _Vm₁y/W_gC_{g g} dimensionless coordinate - m=cw/C_g– _gW_g water equivalent ratio Deceased.Translated from Inzhenerno-Fizicheskii Zhurnal, vol. 20, No. 5, pp. 832–840, May, 1971.     

Effect of molecular weight on the fracture surface energy of poly(methyl methacrylate) in cleavage

By the radiolysis of poly(methylmethacrylate) (PMMA), the fracture surface energy () was determined at room temperature as a function of viscosity average molecular weight (¯M _v). Using a modified parallel cleavage technique, results showed that decreased more than two orders of magnitude with decreasing molecular weight. In the high molecular weight region (¯M _v10⁵), (1×10⁵ erg cm^–2) was relatively insensitive to polymer chain length; whereas for 2.5×10⁴¯M _v 1×10⁵, was strongly dependent on molecular weight. A linear regression analysis in the range ¯M _v=2 to 2.25×10³ indicated that a truly glassy Griffith material was approached for which 750 erg cm^–2. The results confirm the sigmoidal dependence of on molecular weight tested in notched tension. The apparent independence which variations in crack velocity have on with decreasing ¯M _v is shown and explained in terms of the increasingly brittle character of PMMA. Problems associated with the measurement and interpretation of experimental data are considered, particularly with respect to the lower ¯M _v regions.     

Nonsteady thermocapillary mass transfer in the laser alloying of metals

This article examines convective mass transfer of an impurity in a shallow bath of molten metal with allowance for the motion of the fusion front during the laser alloying of metals.Notation r, z, cylindrical coordinates - t time - T_i temperature of the liquid (i=1) and solid (i=2) phases - q(r) absorbed energy flux - k concentration factor - T_m melting point - L heat of fusion - density - _i, _i thermal conductivity and diffusivity - T₀ initial temperature - , absolute and kinematic viscosities of the melt - v_r, v_z projections of the melt velocity on the coordinate axes r and z - p pressure - surface tension Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 56, No. 5, pp. 799–805, May, 1989.     

The tensile strength and ductility of continuous fibre composites

The plastic instability approach has been applied to the tensile behaviour of a continuous fibre composite. It is shown that the combination of two components with different strengths and degrees of work-hardening produces a new material with a new degree of work-hardening, which may be determined by the present analysis. Expressions for the elongation at rupture and the strength of a composite have been obtained and the results of the calculation are compared with some experimental data.List of symbols V _f volume fraction of fibres in composite - , , true strain of fibre, matrix and composite - s true stress - , , nominal stress on fibre, matrix and composite - _*, _*, _* critical stress of fibre, matrix and composite (ultimate tensile strength) - _*, _* critical strain of separate fibre and matrix - _* critical strain of composite - Q external load - A cross-sectional area - A ₀ initial value of area     

Transport properties of helium near the liquid-vapor critical point. II. Thermal conductivity of a3He-4He mixture

Measurements of the thermal conductivity are reported for an 80%³He-20%⁴He mixture above the critical point along several isotherms and near-critical isochores, using the same techniques and apparatus described for a study of³He. Using again the assumption that the observed conductivity can be decomposed into a sum of a regular and a singular contribution _reg and _sing, it is shown that along two near-critical isochores, _sing diverges. In particular, along the isochore showing the largest at the phase transition, the divergence is nearly the same as for³He and can be roughly characterized by a simple power law (T-T_c)^– with 0.58. This observation is contrary to predictions that foresee _sing0 asT _c is approached. The relaxation times characterizing the attainment of steady state conditions after switching the heat flux on and off show a similar behaviour as a function of reduced temperature as do those for pure³He. This result might indicate a substantial coupling between concentration and entropy diffusion. In the Appendix, the correlation length for³He nearT _c is calculated from heat conductivity, viscosity, and specific heat data.Work supported by Grant DMR 8024056 of the National Science Foundation.     

Thermal stresses in a thin circular disc of orthotropic material due to an instantaneous point heat source

Summary In this paper, we have analysed the transient plane thermal stress problem of a circular disc of orthotropic material with instantaneous point heat source. The variation of with time along different radius vectors is exhibited graphically and compared with that of the isotropic case.

---

Wärmespannungen in einer dünnen Kreisscheibe aus orthotropem Material zufolge einer punktförmigen instantanen Wärmequelle

Zusammenfassung In dieser Arbeit wurde das instationäre ebene Wärmespannungsproblem einer Kreisscheibe aus orthotropem Material zufolge einer punktförmigen instantanen Wärmequelle untersucht. Die Veränderung von über die Zeit für verschiedene Radien ist graphisch dargestellt und wird mit dem isotropen Fall verglichen.

---

Nomenclature r, polar coordinates - T temperature rise - ² ratio of conductivities - ² thermal diffusivity in -direction - J _n Bessel function ofn-th order - t time - p introduced in equation (2) - h heat transfer coefficient - a radius of circular disc - f(r, ) temperature distribution at initial state - r ₀, ₀ a point on the disc - T ₀ strength of point heat source - Dirac delta function - F stress function in two dimensions - ₁, ₂ coefficient of thermal expansion - a ₁₁,a ₁₂,a ₂₂,a ₆₆ elastic constants With 9 Figures     

Natural damped frequencies and axial response of a rotating finite viscous liquid column

Summary For a finite solidly rotating cylindrical liquid column the damped natural axisymmetric frequencies have been determined. The liquid was considered incompressible and viscous. The cases of freely slipping edges and that of anchored edges have been treated. It was found that instability appears in a purely aperiodic root for the spinning liquid bridge. This is in contrast to the instability appearing in the damped oscillatory natural frequency of a nonspinning liquid column at . The spinning viscous liquid column exhibits the same instability as the frictionless liquid. It appears at for axisymmetric oscillations.List of symbols a radius of liquid column - I _m modified Bessel function of first kind and orderm - s complex frequency ( ) - r, ,z polar cylindrical coordinates - p pressure - t time - u, v, w radial-, azimuthal- and axial velocities of liquid, respectively - Weber number - h height of liquid column - dynamic viscosity of liquid - v kinematic viscosity of liquid (v=/) - density of liquid - surface tension of liquid - _r , _rz shear stress - (r, z, t) circulation - (r, z, t) streamfunction - ₀ angular velocity of liquid column about the axis of symmetry - (,t) free surface displacement     

Equation of state of3He near its liquid-vapor critical point

We report high-resolution measurements of the pressure coefficient (P/T) for³He in both the one-phase and two-phase regions close to the critical point. These include data on 40 isochores over the intervals–0.1t+0.1 and–0.2+0.2, wheret=(T–T _c )/T _c and =(– _c )/ _c . We have determined the discontinuity (P/T) of (P/T) between the one-phase and the two-phase regions along the coexistence curve as a function of . The asymptotic behavior of (1/) (P/T) versus near the critical point gives a power law with an exponent (+–1)^–1=1.39±0.02 for0.010.2 or–1×10 ^–2t–10 ^–6 , from which we deduce =1.14±0.01, using =0.361 determined from the shape of the coexistence curve. An analysis of the discontinuity (P/T) with a correction-to-scaling term gives =1.17±0.02. The quoted errors are fromstatistics alone. Furthermore, we combine our data with heat capacity results by Brown and Meyer to calculate (/T) _c as a function oft. In the two-phase region the slope (²/T ²)_c is different from that in the one-phase region. These findings are discussed in the light of the predictions from simple scaling and more refined theories and model calculations. For the isochores 0 we form a scaling plot to test whether the data follow simple scaling, which assumes antisymmetry of – ( _c ,t) as a function of on both sides of the critical isochore. We find that indeed this plot shows that the assumption of simple scaling holds reasonably well for our data over the ranget0.1. A fit of our data to the linear model approximation is obtained for0.10 andt0.02, giving a value of =1.16±0.02. Beyond this range, deviations between the fit and the data are greater than the experimental scatter. Finally we discuss the (P/T) data analysis for ⁴ He by Kierstead. A power law plot of (1/) P/T) versus belowT _c leads to =1.13±0.10. An analysis with a correction-to-scaling term gives =1.06±0.02. In contrast to ³ He, the slopes (²/T ²)_c above and belowT _c are only marginally different.Work supported by a grant from the National Science Foundation.     

Response of differently axially excited spinning liquid columns

Summary A solidly rotating finite liquid column consisting of frictionless liquid is subjected to various axial excitation modes. The response of the free liquid surface displacement and velocity distribution has been determined in the elliptic (>2₀) and hyperbolic range (<2₀). Differences of the various cases are presented.List of symbols a radius of liquid bridge - h length of liquid bridge - I ₀,I ₁ modified Besselfunctions - J ₀,J ₁ Besselfunctions - p liquid pressure - r, ,z cylindrical polar coordinates - t time - u, v, w velocity distribution in rotating liquid - axial excitation amplitude - elliptic case (>2₀) - hyperbolic case (<2₀) - liquid density - surface tension - liquid surface displacement - _n damping factor - phase angle - acceleration potential - ₀ rotational speed - , ₁,₂ axial forcing frequency - natural frequency of rotating system - _0n natural frequency of harmonic axial response