On surface waves in generalized thermoelasticity

Knowles' representation theorem for harmonically time-dependent free surface waves on a homogeneous, isotropic elastic half-space is extended to include harmonically time-dependent free processes for thermoelastic surface waves in generalized thermoelasticity of Lord and Shulman and of Green and Lindsay.r , , r , , .This work was done when author was unemployed.     

Interpolation formulas for maxwell viscosity of certain metals as a function of shear-strain intensity and temperature

The purpose of this study is the construction of interpolation formulas for the dependence of Maxwell viscosity, a quantity which is the reciprocal of shear-strain relaxation time , on shear-strain intensity and temperature for several metals: iron, aluminum, copper, and lead. This function was interpolated in various temperature and deformation velocity ranges in accordance with available experimental data for iron (0 10⁷ sec^–1, 200 ° T 1500 °); aluminum (0 10⁷ sec^–1, 300 ° T 900 °); copper (0 10⁵ sec^–1, 300 ° T 1300 °); lead (0 10⁶ sec^–1, 90 ° T 400 °); temperatures in °K.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 114–118, July–August, 1974.     

Extended Korn's inequalities and the associated best possible constants

We are concerned with the coerciveness of the strain energy E(u) (in linear elasticity) associated with a displacement vector u on the Sobolev space H¹ () or its subspaces, a domain in ⁿ representing an isotropic elastic body—certain specific cases are called Korn's inequalities. Sufficient (and necessary) conditions on the Lamé moduli for E(·) to be coercive (or uniformly positive) on such spaces are given, and the associated best possible constants are obtained for some cases.     

Asymptotic behavior of solutions of nonlinear elliptic equations

We study and obtain formulas for the asymptotic behavior as ¦x¦ of C ² solutions of the semilinear equation u=f(x, u), x (*) where is the complement of some ball in ⁿ and f is continuous and nonlinear in u. If, for large x, f is nearly radially symmetric in x, we give conditions under which each positive solution of (*) is asymptotic, as ¦x¦, to some radially symmetric function. Our results can also be useful when f is only bounded above or below by a function which is radially symmetric in x or when the solution oscillates in sign. Examples when f has power-like growth or exponential growth in the variables x and u usefully illustrate our results.     

Hypersonic flow over a narrow delta plate at large angles of attack

The problem of hypersonic flow over a flat delta plate with a high sweepback anglex at angles of attack close to /2 is solved using a numerical algorithm based on transition to the conical solution. The existence of conical flow at /2 with the velocity vector directed towards the apex of the plate is established. Values ofC _p/sin² and the thickness of the shock layer in the plane of symmetry of the plate are given as functions of the hypersonic similarity parameterk=tan tanx. A comparison of the calculated and experimental data shows that they are in good agreement.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.5, pp. 183–185, September–October, 1992.     

The contribution of internal degrees of freedom to the non-Newtonian rheology of model polymer fluids

We report non-equilibrium molecular dynamics simulations of rigid and non-rigid dumbbell fluids to determine the contribution of internal degrees of freedom to strain-rate-dependent shear viscosity. The model adopted for non-rigid molecules is a modification of the finitely extensible nonlinear elastic (FENE) dumbbell commonly used in kinetic theories of polymer solutions. We consider model polymer melts — that is, fluids composed of rigid dumbbells and of FENE dumbbells. We report the steady-state stress tensor and the transient stress response to an applied Couerte strain field for several strain rates. We find that the rheological properties of the rigid and FENE dumbbells are qualitatively and quantitatively similar. (The only exception to this is the zero strain rate shear viscosity.) Except at high strain rates, the average conformation of the FENE dumbbells in a Couette strain field is found to be very similar to that of FENE dumbbells in the absence of strain. The theological properties of the two dumbbell fluids are compared to those of a corresponding fluid of spheres which is shown to be the most non-Newtonian of the three fluids considered.Symbol Definition b dimensionless time constant relating vibration to other forms of motion - F force on center of mass of dumbbell - F _i force on bead i of dumbbell - F force between center of masses of dumbbells and - F _ij force between beads i and j - h vector connecting bead to center of mass of dumbbell - H dimensionless spring constant for dumbbells, in units of / ² - I moment of inertia of dumbbell - J general current induced by applied field - k _B Boltzmann's constant - L angular momentum - m mass of bead, (= m/2) - M mass of dumbbell, g - N number of dumbbells in simulation cell - P translational momentum of center of mass of dumbbell - P pressure tensor - P _xy xy component of pressure tensor - Q separation of beads in dumbbell - Q _eq equilibrium extension of FENE dumbbell and fixed extension of rigid dumbbell - Q ₀ maximum extension of dumbbell - r _ij vector connecting beads i and j - r position vector of center of mass dumbbell - R vector connecting centers of mass of two dumbbells - t time - t ^* dimensionless time, in units of m/ - T ^* dimensionless temperature, in units of /k - u potential energy - u velocity vector of flow field - u _x x component of velocity vector - V volume of simulation cell - X general applied field - strain rate, s^–1 - ^* dimensionless shear rate, in units of /m ² - general transport property - Lennard-Jones potential well depth - friction factor for Gaussian thermostat - shear viscosity, g/cms - ^* dimensionless shear viscosity, in units of m/ ² - ^* dimensionless number density, in units of ^–3 - Lennard-Jones separation of minimum energy - relaxation time of a fluid - angular velocity of dumbbell - orientation angle of dumbbell      

Thermal free convection from a solid sphere

Summary Thermal free convection from a sphere has been studied by melting solid benzene spheres in excess liquid benzene (Pr=8,3; 10⁸<Gr<10⁹). Overall heat transfer as well as local heat transfer were investigated. For the effect of cold liquid produced by the melting a correction has been applied. Results are compared with those obtained by other workers who used alternative experimental methods.Nomenclature coefficient of heat transfer - d characteristic length, here diameter of sphere - thermal conductivity - g acceleration of free fall - cubic expansion coefficient - T temperature difference between wall and fluid at infinity - kinematic viscosity - density - c specific heat capacity - a thermal diffusivity (=/c) - D diffusion coefficient - Nu dimensionless Nusselt number (=d/) - Nu* the analogous number for mass transfer (=kd/D) - mean value of Nusselt number - Gr dimensionless Grashof number (=gd ³T/ ²) - Gr* the analogous number for mass transfer (=gd ³x/ ²) - Pr dimensionless Prandtl number (=/a) - Sc dimensionless Schmidt number (=/D)     

Mean and turbulent velocity measurements of supersonic mixing layers

The behavior of supersonic mixing layers under three conditions has been examined by schlieren photography and laser Doppler velocimetry. In the schlieren photographs, some large-scale, repetitive patterns were observed within the mixing layer; however, these structures do not appear to dominate the mixing layer character under the present flow conditions. It was found that higher levels of secondary freestream turbulence did not increase the peak turbulence intensity observed within the mixing layer, but slightly increased the growth rate. Higher levels of freestream turbulence also reduced the axial distance required for development of the mean velocity. At higher convective Mach numbers, the mixing layer growth rate was found to be smaller than that of an incompressible mixing layer at the same velocity and freestream density ratio. The increase in convective Mach number also caused a decrease in the turbulence intensity ( _u/U).List of symbols a speed of sound - b total mixing layer thickness between U ₁ – 0.1 U and U ₂ + 0.1 U - f normalized third moment of u-velocity, f u³/(U)³ - g normalized triple product of u² , g u²/(U)³ - h normalized triple product of u ², h u ²/(U)³ - l _u axial distance for similarity in the mean velocity - l _u axial distance for similarity in the turbulence intensity - M Mach number - M _c convective Mach number (for ₁ = ₂), M _c (U ₁ – U ₂)/(a ₁ + a ₂) - P static pressure - r freestream velocity ratio, r U ₂/U ₁ - Re unit Reynolds number, Re U/ - s freestream density ratio, s ₂/₁ - T _t total temperature - u instantaneous streamwise velocity - u deviation of u-velocity, uu – U - U local mean streamwise velocity - U ₁ primary freestream velocity - U ₂ secondary freestream velocity - average of freestream velocities, (U ₁ + U ₂)/2 - U freestream velocity difference, U U ₁ – U ₂ - instantaneous transverse velocity - v deviation of -velocity, – V - V local mean transverse velocity - x streamwise coordinate - y transverse coordinate - y ₀ transverse location of the mixing layer centerline - ensemble average - ratio of specific heats - boundary layer thickness (y-location at 99.5% of free-stream velocity) - similarity coordinate, (y – y ₀)/b - compressible boundary layer momentum thickness - viscosity - density - standard deviation - dimensionless velocity, (U – U ₂)/U - 1 primary stream - 2 secondary stream A version of this paper was presented at the 11th Symposium on Turbulence, October 17–19, 1988, University of Missouri-Rolla     

Certain yawed magnetogasdynamic flows

Certain steady yawed magnetogasdynamic flows, in which the magnetic field is everywhere parallel to the velocity field, are related to certain reduced three-dimensional compressible gas flows having zero magnetic field. Under a restriction, the reduced flows are linked, by certain reciprocal relations, to a four parameter class of plane gas flows. In the instance of constant entropy an approximation method is suggested for obtaining magnetogasdynamic flows from the corresponding plane, irrotational gasdynamic flows and examples are given.

Nomenclature

magnetogasdynamic flow variables H magnetic intensity - q fluid velocity - fluid density - p pressure - s entropy - Q _t, H _t component of q, H in the x–y plane - w , h component of q, H perpendicular to the x–y plane reduced gasdynamic flow factor of proportionality - q* fluid velocity - * fluid density - p* pressure - Q _t ^* =u*î+v*, w* components of q* - l arbitrary constant - A _v Alfvén speed - Q _t, , p fluid velocity, density, pressure of the reciprocal gas dynamic flow - L, n, k, arbitrary constants - , velocity potential, stream function - angle made by Q _t, Q _t ^* , and V with the x-axis - adiabatic gas constant - a ²=(–1)/2 constant - M Mach number - W constant value of w* - E approximate constant value of g(p) - * modified potential function - modified velocity coordinate - +i - complex potential of the irrotational flow - B arbitrary constant - V incompressible flow velocity - V modified fluid velocity - X _p, Y _p points on the profile     

Transport in ordered and disordered porous media I: The cellular average and the use of weighting functions

In this work we consider transport in ordered and disordered porous media using singlephase flow in rigid porous mediaas an example. We defineorder anddisorder in terms of geometrical integrals that arise naturally in the method of volume averaging, and we show that dependent variables for ordered media must generally be defined in terms of thecellular average. The cellular average can be constructed by means of a weighting function, thus transport processes in both ordered and disordered media can be treated with a single theory based on weighted averages. Part I provides some basic ideas associated with ordered and disordered media, weighted averages, and the theory of distributions. In Part II a generalized averaging procedure is presented and in Part III the closure problem is developed and the theory is compared with experiment. Parts IV and V provide some geometrical results for computer generated porous media.Roman Letters A interfacial area of the- interface contained within the macroscopic region, m² - A_e area of entrances and exits for the-phase contained within the macroscopic system, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor, m² - L general characteristic length for volume averaged quantities, m - characteristic length (pore scale) for the-phase - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - N unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p0 reference pressure in the-phase, N/m² - p traditional intrinsic volume averaged pressure, N/m² - r0 radius of a spherical averaging volume, m - r position vector, m - r position vector locating points in the-phase, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - v velocity vector in the-phase, m/s - v traditional superficial volume averaged velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V/V, volume average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

The theory of a vibrating-rod viscometer

The paper presents a complete theory for a viscometer based upon the principle of a circular-section rod, immersed in a fluid, performing transverse oscillations perpendicular to its axis. The theory is established as a result of a detailed analysis of the fluid flow around the rod and is subject to a number of criteria which subsequently constrain the design of an instrument. Using water as an example it is shown that a practical instrument can be designed so as to enable viscosity measurement with an accuracy of ±0.1%, although it is noted that many earlier instruments failed to satisfy one or more of the newly-established constraints.Nomenclature A, D constants in equation (46) - A _m , B _m , C _m , D _m constants in equations (50) and (51) - A _j , B _j constants in equation (14) - a _j ⁺ , a _j ^– wavenumbers given by equation (15) - C _f drag coefficient defined in equation (53) - c speed of sound - D _b drag force of fluid b - D ₀ coefficient of internal damping - E extensional modulus - f(z) initial deformation of rod - f(), F _m () functions of defined in equation (41) - F force in the rod - force per unit length near t=0 - F dimensionless force per unit length near t=0 - g _m amplitude of transient force - G modulus of rigidity - h, h* functions defined by equations (71) and (72) - H functions defined by equation (69) and (70) - I second moment of area - I _0,1, J _0,1, K _0,1 modified Bessel functions - k, k functions defined in equations (2) - L half-length of oscillator - Ma Mach number - m _b added mass per unit length of fluid b - m _s mass per unit length of solid - n _j eigenvalue defined in equations (15) and (16) - R radius of rod - R _c radius of container - r radial coordinate - T tension - T _visc temperature rise due to heat generation by viscous dissipation - t time - v _r , v radial and angular velocity components - y lateral displacement - y ₀ initial lateral displacement - y ₁, y ₂ successive maximum lateral displacement - z axial coordinate - dimensionless tension - dimensionless mass of fluid - dimensionless drag of fluid - amplification factor - logarithmic decrement in a fluid - _a , _b logarithmic decrement in fluids a and b - ₀ logarithmic decrement in vacuo - _j logarithmic decrement in mode j in a fluid - spatial resolution of amplitude - _v voltage resolution - r, , , _s, , increments in R, , , _s , , - dimensionless amplitude of oscillation - dimensionless axial coordinate - angular coordinate - _f thermal conductivity of fluid - viscosity of fluid - viscosity of fluid calculated on assumption that * - _a , _b viscosity of fluids a and b - _m constants in equation (10) - dimensionless displacement - _j j the component of - density of fluid - _a , _b density of fluids a and b - _s density of tube or rod material - dimensionless radial coordinate - * dimensionless radius of container - dimensionless times - spatial component of defined in equation (11) - _j , _tm jth, mth component of - dimensionless streamfunction - ⁰, ¹ components of in series expansion in powers of - streamfunction - dimensionless frequency (based on ) - angular frequency - ₀ angular frequency in absence of fluid and internal damping - _j angular frequency in mode j in a fluid - _a , _b frequencies in fluids a and b     

Effect of suspended particles on the stability of plane poiseuille flow

The behavior of the neutral stability curves is investigated for various values of the particle relaxation time and mass concentration 0 100 and 0 f 0.1. It is shown that as increases from zero the flow is at first destabilized and then at >6 becomes stable, while at >40 the stabilizing effect of the dispersed phase grows weaker. It is found that there is a certain interval 10< <40 on which the flow is most stable.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 46–53, January–February, 1986.     

Positively invariant regions for a problem in phase transitions

Positively invariant regions for the system v _t + p(W) _x = V _xx , W _t –V _x = W _xx are constructed where p < 0, w < , w > , p(w) = 0, w , > 0. Such a choice of p is motivated by the Maxwell construction for a van der Waals fluid. The method of an analysis is a modification of earlier ideas of Chueh, Conley, & Smoller [1]. The results given here provide independent L bounds on the solution (w, v).Dedicated to Professor James Serrin on the occasion of his sixtieth birthday     

Representation of the rheological properties of polymer melts in terms of complex fluidity

In dynamic rheological experiments melt behavior is usually expressed in terms of complex viscosity ^* () or complex modulusG ^* (). In contrast, we attempted to use the complex fluidity ^* () = 1/µ ^* () to represent this behavior. The main interest is to simplify the complex-plane diagram and to simplify the determination of fundamental parameters such as the Newtonian viscosity or the parameter of relaxation-time distribution when a Cole-Cole type distribution can be applied. ^* () complex shear viscosity - () real part of the complex viscosity - () imaginary part of the complex viscosity - G ^* () complex shear modulus - G() storage modulus in shear - G() loss modulus in shear - J ^* () complex shear compliance - J() storage compliance in shear - J() loss compliance in shear - shear strain - rate of strain - angular frequency (rad/s) - shear stress - loss angle - ^* () complex shear fluidity - () real part of the complex fluidity - () imaginary part of the complex fluidity - ₀ zero-viscosity - ₀ average relaxation time - h parameter of relaxation-time distribution     

Filament stretching equation in Lagrangian coordinates

The general integral of the very simple equation ^21/n/() was found to describe the cross sectional area of filaments of isothermal power law fluids while in transient stretching where is time and is the initial location of fluid molecules at time = 0 given as the distance from a reference point fixed in space. Any such stretching transient given as a solution of the above equation is physically realizable subject to the restrictions > 0 and/ < 0.     

Zur Variationsformulierung nichtlinearer Randwertprobleme

Übersicht Es werden verschiedene Bedingungen aufgestellt, die es erlauben, die durch die beiden (Systeme von) nichtlinearen DifferentialgleichungenA (u, ) = q, B (u, ) = und Randbedingungen zusammen mit den nichtlinearen algebraischen Relationenq = C(u, ), = D(u, ) beschriebene Aufgabe durch äquivalente Variationsprobleme zu ersetzen. Dabei zeigt sich ein enger Zusammenhang mit den in der Festkörpermechanik wohlbekannten Prinzipien der virtuellen Verschiebungen und der virtuellen Kräfte. Die auf systematischem Weg konstruierten Variationsfunktionale enthalten viele in der Physik bekannte Funktionale als Sonderfälle, insbesondere jene, die in der Elastomechanik nach Green, Castigliano, Hellinger, Reißner, Hu und Washizu benannt werden.

---

Summary In this paper there are established various conditions which allow a variational formulation of the problem described by the two (systems of) nonlinear differential equationsA(u, ) = q, B(u, ) = and boundary conditions together with the nonlinear algebraic relationsq = C(u, ), = D(u, ). Besides a close relationship is revealed to the principles of virtual displacements and virtual forces which are wellknown in solid mechanics. The systematically constructed variational functional contain many functionals in physics as special cases, mainly those of Green, Castigliano, Hellinger, Reißner, Hu and Washizu in elastomechanics.

     

Die Torsionsspannungen in Wellen mit tiefen Längsnuten

Übersicht MitF(x, y) als Spannungsfunktion einer Welle ohne Nut und(, y) als Potentialfunktion des Quelle-Senke-Systems erhält man Spannungsfunktionen(, y) =F(x, y) –(, y) für Wellen mit tiefen Längsnuten. Es wird gezeigt, daß sich damit die Schubspannungen in den Läufern von Schraubenverdichtern ermitteln lassen.

---

Shearing stresses in shafts with deep longitudinal grooves

Summary The stress functions(, y) of shafts with deep longitudinal grooves may be represented by(, y) =F(x, y) –(, y) whereF(x, y) is the stress function of a cylindrical shaft without grooves and(, y) denotes the potential function of the source-sink system. It is shown that the shearing stresses in rotors of screw-compressors may be obtained in this way.

     

Free convection heat transfer from an isothermal plate with arbitrary inclination

This paper presents a new formulation for the laminar free convection from an arbitrarily inclined isothermal plate to fluids of any Prandtl number between 0.001 and infinity. A novel inclination parameter is proposed such that all cases of the horizontal, inclined and vertical plates can be described by a single set of transformed equations. Moreover, the self-similar equations for the limiting cases of the horizontal and vertical plates are recovered from the transformed equations by setting=0 and=1, respectively. Heated upward-facing plates with positive and negative inclination angles are investigated. A very accurate correlation equation of the local Nusselt number is developed for arbitrary inclination angle and for 0.001 Pr .

---

Wärmeübertragung bei freier Konvektion an einer isothermen Platte mit beliebiger Neigung

Zusammenfasssung Diese Untersuchung stellt eine neue Formulierung der laminaren freien Konvektion von Flüssigkeiten mit einer Prandtl-Zahl zwischen 0,001 und unendlich an einer beliebig schräggestellten isothermen Platte dar. Ein neuer Neigungsparameter wird eingeführt, so daß alle Fälle der horizontalen, geneigten oder vertikalen Platte von einem einzigen Satz transformierter Gleichungen beschrieben werden können. Die unabhängigen Gleichungen für die beiden Fälle der horizontalen and vertikalen Platte wurden für=0 und=1 aus den transformierten Gleichungen wieder abgeleitet. Es wurden erwärmte aufwärtsgerichtete Platten mit positiven und negativen Neigungswinkeln untersucht. Eine sehr genaue Gleichung wurde für die lokale Nusselt-Zahl bei beliebigen Neigungswinkeln und für 0,001 Pr entwickelt.

---

Nomenclature C _p specific heat - f reduced stream function - g gravitational acceleration - Gr local Grashof number,g(T _w –T _w ) x³/v² - h local heat transfer coefficient - k thermal conductivity - n constant exponent - Nu local Nusselt number,hx/k - p pressure - Pr Prandtl number, v/ - Ra local Rayleigh number,g(T _w –T )J x³/v - T fluid temperature - T _w wall temperature - T temperature of ambient fluid - u velocity component in x-direction - v velocity component in y-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - (Ra¦sin¦)^1/4/( Ra cos()^1/5 - pseudo-similarity variable, (y/) - dimensionless temperature, (T–T )/(T _w–T ) - ( Ra cos)^1/5+(Rasin)^1/4 - v kinematic viscosity - 1/[1 +(Ra cos)^1/5/( Ra¦sin)^1/4] - density of fluid - Pr/(1+Pr) - _w wall shear stress - angle of plate inclination measured from the horizontal - stream function - dimensionless dynamic pressure     

Die Eigenschaften von Wasser und Wasserdampf nach „The 1968 IFC Formulation“

Zusammenfassung Diese Arbeit enthält Druck-Temperatur-Diagramme für 6 spezifische Zustandsgrößen und 16 erste Ableitungen und zusammengesetzte Größen von Wasser und Wasserdampf, die nach einem Gleichungssystem berechnet wurden, das unter dem Namen The 1968 IFC Formulation for Scientific and General Use von der 6. Internationalen Konferenz für die Eigenschaften des Wasserdampfes angenommen wurde. Einige Konsequenzen der thermodynamischen Konsistenz, das Verhalten im kritischen Gebiet und bei sehr kleinen Drücken werden diskutiert. Ferner werden die kinematische Viskosität und die Temperaturleitfähigkeit, sowie eine Beziehung zwischen dynamischer Viskosität und isenthalpem Drosselkoeffizienten angegeben.

---

This paper contains pressure-temperature diagrams for 6 properties and 16 first derivatives and combined terms for water and steam. These were calculated from a system of equations accepted by the 6^th International Conference on the Properties of Steam, and called The 1968 IFC Formulation For Scientific and General Use. Some consequences of thermodynamic consistency, and the behaviour in the critical region and at very small pressures are discussed. Further, the kinematic viscosity and the thermal diffusivity and a relation between the dynamic viscosity and the throttling coefficient at constant enthalpy are given.

---

Bezeichnungen (s. auch Tabelle 1) k Temperaturleitfähigkeit:k=/c _p - p Druck - r spezifische Verdampfungsenthalpie:r=h –h - T thermodynamische oder Kelvin-Temperatur - t Celsius-Temperatur - dynamische Viskosität - Wärmeleitfähigkeit - v kinematische Viskosität:=/ - Dichte:=1/v Indices und Sättigungswerte des Dampfes und der Flüssigkeit Differenz der Sättigungswerte, z. B. h=h –h     

Indentation of the semi-infinite elastic solid by a concave rigid punch

A solution is obtained for the relationship between load, displacement and inner contact radius for an axisymmetric, spherically concave, rigid punch, indenting an elastic half-space. Analytic approximations are developed for the limiting cases in which the ratio of the inner and outer radii of the annular contact region is respectively small and close to unity. These approximations overlap well at intermediate values. The same method is applied to the conically concave punch and to a punch with a central hole. , , . , . . .