A note on the turbulent energy as a parameter of turbulent flow

The art of modeling turbulence is a needed tool in the construction of computer codes for turbulent flows. The state to which this art has been developed is inadequate, and quotations from authoritative sources support this point of view. The energy contained in the turbulent fluctuations, i.e., the turbulent energy, is often used as a parameter in the modeling process. The present article attempts to examine this quantity as it is being created, transported, and dissipated. For this purpose experimental evidence from the author's own experiments (free jets), as well as theoretical conclusions from the elementary deductions of the basic equations, the concept of turbulent potential flow, and a general solution to the Navier-Stokes-Reynolds equations, is drawn to attention. Recirculating flow is given special attention. The paper concludes with recommendations for principles that must be satisfied if improved modeling is to be achieved. These principles are necessary; whether they are also sufficient is open to question.Nomenclature A ₀ Constant - b _1/2 Jet's half-width - b ₁ ^2/(0) Jet's half-width at z=z(in0) - E _z Kinetic energy contained in the jet's axial velocity at a given profile - E _r Kinetic energy contained in the jet's radial velocity at a given profile - f() Dimensionless velocity profile [f(0)=1] - F(), H() Defined functions - L _char Jet's characteristic length - m, n Exponents - p Pressure - q Kinetic energy in the turbulent fluctuations - Heat flux - q ² - r, , z Cylindrical coordinates - t Time - û Internal energy - u, v, w Velocity components - Mean velocity components - Mean velocity components - U ₀ Constant - U _plate Plate's velocity - Uskc/(0) Centerline velocity at z=z₀ - X, Y, Z Components of body force - W Total work done by surface stresses - W ₁ Recoverable work done by surface stresses - W ₂ Dissipated work - z ₀ Downstream distance from the nozzle beyond which self-similar velocity profiles occur - Fluid's kinematic viscosity - Fluid's density - Normal stresses - Shear stresses - Normal stresses with the pressure removed - Dimensionless Crossflow coordinate - ₀ Constant - Stress functions - Stress potential Paper dedicated to Professor Joseph Kestin.Definitions of symbols are given under Nomenclature.     

Bounding-surface plasticity: Stress space versus strain space

Summary A bounding-surface plasticity model is formulated in stress space in a general enough manner to accommodate a considerable range of hardening mechanisms. Conditions are then established under which this formulation can be made equivalent to its strain-space analogue. Special cases of the hardening law are discussed next, followed by a new criterion to ensure nesting. Finally, correlations with experimental data are investigated.Notation (a) centre of the stress-space (strain-space) loading surface; i.e., backstress (backstrain) - ^* (a ^*) centre of the stress-space (strain-space) bounding surface - (a ) target toward which the centre of the stress-space (strain-space) loading surface moves under purely image-point hardening - (b) parameter to describe how close the loading surface is to nesting with the bounding surface in stress (strain) space; see (H10) - (c) elastic compliance (stiffness) tensor - (d) parameter to describe how close the stress (strain) lies to its image point on the bounding surface; see (H10) - (D) generalised plastic modulus (plastic compliance); see (1) - function expressing the dependence of the generalised plastic modulus on (plastic complianceD ond) - ^* (D ^*) analogue to (D) for the bounding surface - function expressing the dependence of ^* on (D ^* ond) - () strain (stress) - ' (') deviatoric strain (stress) - ^P ( ^R ) plastic strain (stress relaxation); see Fig. 1 - () image point on the bounding surface corresponding to the current strain (stress) - _iso (f _iso) at the point of invoking consistency, the fraction of local loading-surface motion arising from a change of radius; i.e., fraction of isotropic hardening in the stress-space theory - _kin (f _kin) at the point of invoking consistency, the fraction of local loading-surface motion arising from a change in the backstress (backstrain); i.e., fraction of kinematic hardening in the stress-space theory - _nor (f _nor) at the point of invoking consistency, the fraction of backstress (backstrain) motion directed toward the image stress (strain); i.e., the image-point fraction of the kinematic hardening in the stress-space theory - _ima (f _ima) at the point of invoking consistency, the fraction of backstress (backstrain) motion directed toward the image stress (strain); i.e., the image-point fraction of the kinematic hardening in the stress-space theory - function relating _iso to , , and (f _iso tob,d, andl) - function relating _kin to , , and (f _kin onb,d, andl) - function relating _nor to , , and (f _nor onb,d, andl) - function relating _ima to , , and (f _ima onb,d, andl) - the fraction of outwardly normal bounding-surface motion at the Mróz image point which arises from a change of radius - the fraction of outwardly normal bounding-surface motion at the Mróz image point which arises from a change in the centre - function relating _iso ^* to (f _iso ^* tod) - function relating _kin ^* to (f _kin ^* tod) - (l) parameter to describe the full extent of plastic loading up to the present, giving the arc length of plastic strain (stress relaxation) trajectory; see (H10) - function relating the direction for image-point translation of the loading surface to various other tensorial directions associated with the current state; see (H5). With 6 Figures     

Higher approximations for supersonic flow past slowly oscillating bodies of revolution

Summary Supersonic flow past slowly oscillating pointed bodies of revolution is studied. Starting from the complete nonlinear potential equation an elementary linearized solution is discussed and it is shown how this solution together with the method of matched asymptotic expansions can be used to derive an elementary second-order slender body theory. This approach is further demonstrated for the oscillating cone and its range of validity is evaluated by comparison with other theoretical methods.

---

Zusammenfassung Es wird die Überschallströmung um langsam schwingende spitze Rotationskörper untersucht. Ausgehend von der vollständigen nichtlinearen Potentialgleichung wird zuerst eine elementare linearisierte Lösung besprochen und gezeigt, wie diese Lösung im Verein mit der Method of matched asymptotic expansions zur Herleitung einer elementaren Schlankkörpertheorie zweiter Ordnung verwendet werden kann. Die Theorie wird am Beispiel des schwingenden Kegels näher erläutert und mit anderen Methoden verglichen.

---

Symbols a Velocity of sound - c _N Normal force coefficient - Damping coefficient - F _(x) Dipole distribution - k Reduced frequency - M Mach number - R (x) Meridian profile - t Time - x, r, Cylindrical coordinates - - Ratio of specific heats - Amplitude of oscillation - Thickness ratio - Perturbation potential - Zero angle of attack potential - æ - Velocity potential - Out-of-phase potential - - In-phase potential - - Source coordinate With 4 Figures     

Heat flow model of rapid solidification with nonhomogeneous thermal contact

A heat flow model is presented of the solidification process of a thin melt layer on a heat conducting substrate. The model is based on the two-dimensional heat conduction equation, which was solved numerically. The effect of coexisting regions of good and bad thermal contact between foil and substrate is considered. The numerical results for thermal parameters of the Al-Cu eutectic alloy show considerable deviations from one-dimensional solidification models. Except for drastic differences in the magnitude of the solidification rate near the foil-substrate interface, the solidification direction deviates from being perpendicular to the substrate and large lateral temperature gradients occur. Interruption of the thermal contact may lead to back-melting effects. A new quantity, the effective diffusion length, is introduced which allows some conclusions to be drawn concerning the behaviour of the frozen microstructure during subsequent cooling.Nomenclature _i ,a _i Thermal diffusivity _i = _i /c _{i i} ,a _i = _i / ₁ - c _i Specific heat capacity - d Foil thickness - D Solid state diffusion coefficient - e_x, e_z Unit vectors - H Latent heat of fusion - h ,h Foil-substrate heat transfer coefficients - i Index: 1, melt; 2, solidified foil; 3, substrate - _i ,k _i Thermal conductivityk _i = _i / ₁ - n Normal unit vector - Nu ,Nu Nusselt numbers for regions of badNu(x,) and good thermal contact, respectivelyNu =h Nu _d / ₁,,Nu(x, )=h(x,)d/ ₁ - R Universal gas constant - , s Position of the liquid-solid interface ¯s/d=s=s _xe_x+s _ze_z - Local solidification rate /d = s =s _xe_x +s _ze_z - t Real time - T _i Temperature field - T ₀ Ambient temperature - T _f Melting temperature - u _i Dimensionless temperature fieldu _i (x, z,)=T _i (x,z,)/T _f - u ₀ Dimensionless ambient temperatureu ₀=T ₀/T _f - _i Local cooling rate within the foil _i = du _i /d - W Stefan numberW=H/c ₁ T _f - ,x Cartesian coordinate parallel to the foil-substrate interfacex= /d - ₀,x ₀ Lateral extension of foil sectionx ₀= ₀/d - ₁,x ₁ Lateral contact lengthx ₁= ₁/d - ,z Cartesian coordinate perpendicular to the foil-substrate interfacez= /d - ₀,z ₀ Substrate thicknessz ₀= ₀/d - E Activation energy of diffusion - T Initial superheat of the melt - u Dimensionless initial superheat u=T/T _f - (x) Step function - _eff Dimensionless effective diffusion length - _i Mass density - Dimensionless time=t ₁/d ² - _f, _f(x, z) Total and local dimensionless freezing time, respectively     

Angular dependence of the critical magnetic field in granular superconducting films

The critical magnetic fieldH _c () of granular Al films has been measured as a function of the angle between the field direction and the plane of the film at temperatures nearT _c0 .The film thicknessd is smaller than the temperature-dependent coherence length (T), the bulk electron mean free path1 is smaller than the BCS coherence length ₀, and 1 d. The experimental data onH _c () are well fitted by the Tinkham formula. However, the observed values ofH _c/H _care not always consistent with and increase with1/d. This fact suggests that the boundary scattering of electrons at the film surface enhancesH _c () and that the enhancement ofH _cis larger than that ofH _c.On leave from Department of Physics, Faculty of Science, Kyushu University, Fukuoka, Japan.     

On strong Dickson pseudoprimes

It is known that the Lucas sequenceV _n(,c)=aⁿ + bⁿ,a, b being the roots ofx ² – x + c=0 equals the Dickson polynomial .^n–2i Lidl, Müller and Oswald recently defined a number b to be a strong Dickson pseudoprime to the parameterc (shortlysDpp(c)) if [itg_n(b, c)b modn for all b. These numbers seem to be very appropriate for a fast probabilistic prime number test. In generalizing results of the above mentioned authors a criterion is derived for an odd composite number to be ansDpp(c) for fixedc. Furthermore the optimal parameterc for the prime number test is determined.     

Nonstationary process of heat transfer in a tube with longitudinal fins

An analytical solution to the problem of nonstationary thermal interaction of a flow of a heat-transfer agent and a thin-walled tube with longitudinal fins is constructed for variable parameters of heat transfer.Notation u, temperatures of the fins - ,w temperatures of the tube walls - temperature of the flow of the heat-transfer agent - _i ,i= coefficients of heat transfer from the ambient medium to the fins and the tube walls, respectively - _i ,i= temperature distributions for the ambient medium - coefficients of heat transfer from the flow of the heat-transfer agent to the tube walls - q _i density of the heat flux to the corresponding portions of the tube - heat capacity, thermal conductivity, density, and thickness of the fin and tube material - c _p , ,G, F heat capacity, density, and flow rate of the heat-transfer agent, cross-sectional area of the tube - dimensions of the tube Bauman Moscow State Technical University. Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 66, No. 6, pp. 673–680, June, 1994.     

An asymptotic solution for large Prandtl number free convection

Summary The set of ordinary differential equations governing free convection boundary layer flow past an isothermal semi-infinite vertical flat plate is solved for large Prandtl numbers by means of the method of matched asymptotic expansions. The analysis leads to an expression for heat transfer which contains the Prandtl number explicitly and which is very accurate for sufficiently large values of the Prandtl number. On the other hand the analysis also has qualitative assets. Before choosing the mathematical method of solution, the physical aspects of the large Prandtl number free convection boundary layer are investigated. The mathematical solution serves to enlarge our understanding of the physical implications of a free convection boundary layer in a large Prandtl number fluid.Nomenclature a_ij coefficient defined by - b_ij coefficient defined by F_j()=b_0j+b_1j +b_2j ²+.... - c coefficient defined by equation (3) - c_p specific heat - f non-dimensional stream function of inner expansion (7) - f_n n-th perturbation of f - F non-dimensional stream function of outer expansion (15) - g non-dimensional stream function (1) - ¯g acceleration due to gravity - Gr_x local Grashof number:g(T_w–T)x³/v² - h non-dimensional temperature (2) - k coefficient of heat conduction - Nu_x local Nusselt number: - T temperature - T_w wall-temperature - T ambient temperature - u longitudinal velocity - x co-ordinate measuring distance from the leading edge - y co-ordinate measuring distance normal to the plate Greek symbols coefficient of thermal expansion - _i expansion parameter (21) - expansion parameter (22) - _i expansion parameter (33) - expansion parameter (34) - expansion parameter: ^–1/2 - inner similarity co-ordinate (9) - non-dimensional temperature of inner expansion (8) - _n n-th perturbation of - non-dimensional temperature of outer expansion (16) - _n n-th perturbation of - similarity co-ordinate (3) - kinematic viscosity - outer similarity co-ordinate (17) - density - Prandtl number:c_p/k - stream function     

Hydrodynamics of a vortex ring

We considered the kinematics and dynamics of a vortex ring in an incompressible fluid in toroidal coordinates. We obtained the change in the pressure difference along the boundary between two flow regions in the case of a moving torus.Notation , , toroidal coordinates - (V ;V ;V ) velocity of a fluid particle and its projections in toroidal coordinates - g ,g ,g metric tensor components - the Jacobian of transition to curvilinear coordinates - V ₀ velocity at the center of a vortex ring on its symmetry axis - x, y, z Cartesian coordinates - z, y, cylindrical coordinates - a distance from the axis of a torus (V=0) to its axis of symmetry (Oz) - angle between the Oy axis and the line that connects a fluid particle on the streamline =const, which represents a circle [16], with the center of this circle - U _z,U _y velocities in the cylindrical system of coordinates - ₀ stream function of a stationary vortex ring - velocity circulation - U V ₁, velocity of a rectilinear flow at infinity - ₁ stream function of a rectilinear flow - = ₀ + ₁ superposition of two flows - n=k ⁴=V ₁/V ₀ velocity ratio coefficient - R radius of a vortical region - U velocity of fluid particles at the boundary in polar coordinates (r, ) with the center at the coordinate origin (point 0) - fluid density - p ₀,p pressure at infinity and at a certain point of flow - pressure difference Polotsk State University, Polotsk, Belarus. Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 68, No. 4, pp. 531–536, July–August, 1995.     

The block-transitive,point-imprimitive 2-(729, 8, 1) designs

The block-transitive point-imprimitive 2-(729,8,1) designs are classified. They all have full automorphism group of order 729.13 which is an extension of a groupN of order 729, acting regularly on points, by a group of order 13. There are, up to isomorphism, 27 designs withN elementary abelian, 13 designs withN=Z ₉ ³ and 427 designs withN the relatively free 3-generator, exponent 3, nilpotency class 2 group, a total of 467 designs. This classification completes the classification of block-transitive, point-imprimitive 2-(, k, 1) designs satisfying , which is the Delandtsheer-Doyen upper bound for the number of points of such designs. The only examples of block-transitive, point-imprimitive 2-(, k, 1) designs with are the 2-(729, 8, 1) designs constructed in this paper.The first three authors acknowledge the support of an Australian European Awards Program scholarship, a Deutsche Akademische Austauschdienst scholarship, and an Australian Research Council Research Fellowship, respectivelyThe authors wish to thank Brendan McKay for his independent verification of the non-isomorphism of the classes of designs found, and of their automorphism groups, using different, nauty techniques [6].     

Stresses in a cylindrical shell supported by a core of a different material

Summary The stress problem of a thin cylindrical shell supported by an elastic core of a different material and subjected to arbitrary loading on its curved surface is considered. The problem is solved by applying the three-dimensional theory of elasticity to the core and using membrane or bending solutions for the shell. Equilibrium and compatibility equations are satisfied at the junction of the shell and the core. It is pointed out that the procedure can easily be extended to the case of a hollow core with or without another shell of another material in it. Numerical results are presented to illustrate the effectiveness of even a weak core in reducing the shell stresses.

---

Zusammenfassung Gegenstand der Untersuchung ist eine dünne Kreiszylinderschale, die durch einen elastischen Kern aus einem anderen Werkstoff gestützt ist und eine beliebige Belastung trägt. Die Lösung verbindet die strenge, dreidimensionale Theorie des zylindrischen Kerns mit der Membran- oder Biegetheorie der Schale. An der Grenze zwischen beiden Teilen müssen die Verschiebungen und gewisse Spannungskomponenten stetig übergehen. Es wird darauf hingewiesen, daß die Lösung leicht auf den Fall ausgedehnt werden kann, daß der Kern ein Hohlzylinder ist, der möglicherweise auf der Innenseite mit einer zweiten Zylinderschale verbunden ist. Zahlenergebnisse zeigen, daß selbst ein verhältnismäsig nachgiebiger Kern einen großen (und günstigen) Einfluß auf die Spannungen in der Schale ausübt.

---

Principal Symbols a Radius of the middle surface of the shell - t Thickness of the shell - =1–t/2a - u _c,v _c,w _c Displacements respectively in the axial, circumferential and radial directions of a point in the core - X(x), (), (r/a) 3×3 square matrices - ,m Parameters - l Length of the cylinder - c A vector containing constantsc ₁,c ₂ andc ₃ - =r/a - =m+4(1–v _e) - E _c,v _e Elastic constants for the core material - Stresses at a point in the core - D _c - A vector containing _rx , _r and _r - (r/a) A 3×3 matrix - Displacements at the surfacer=a of the core - A vector containing - Amplitudes of displacements - A vector containing - =(x, ,a) - _ij Constants - A A square matrix containing constants _ij - Stress resultants in the shell as defined in reference [3] - p _x,p ,P _r Components of applied loading per unit area of shell's middle surface - () - ()· - u, v, w Displacements of a point on the middle surface of the shell - E _s,v _s Elastic constants for the shell material - D _s - K - k - p _xmn,p _mn,p _rmn Amplitudes of loadsp _x,p , p_r - u _mn, v_mn,w _mn Amplitudes of displacementsu, v, With 1 Figure     

Free convection from a heated vertical cylinder embedded in a fluid-saturated porous medium

Summary We consider the free convection boundary layer flow induced by a heated vertical cylinder which is embedded in a fluid-saturated porous medium. The surface of the cylinder is maintained at a temperature whose value above the ambient temperature of the surrounding fluid varies as then ^th power of the distance from the leading edge. Asymptotic analyses and numerical calculations are presented for the governing nonsimilar boundary layer equations and it is shown that, whenn<1, the asymptotic flowfield far from the leading edge of the cylinder takes on a multiple-layer structure. However, forn>1, only a simple single layer is present far downstream, but a multiple layer structure exists close to the cylinder leading edge. We have shown that the fully numerical and asymptotic calculations are in stisfactory agreement, especially for exponentsn close to zero. Comparisons of the present numerical solutions obtained using the Keller-box method with previous numerical solutions using local methods are also given.List of symbols a radius - scaled streamfunctions - f ₀,f ₁,f ₂ inner zone streamfunctions whenn<1 - leading order streamfunctions inn>1, 1 asymptotic solution - F ₀,F ₁ outer zone streamfunctions whenn<1 - G large parameter satisfyingG=X ² lnG - g gravitational acceleration - K permeability of the porous medium - n exponent in prescribed temperature law - r radial co-ordinate - r rescaled radial co-ordinate - R Darcy-Rayleigh number - T temperature of convective fluid - T _w temperature of cylinder at leading edge - T ambient temperature of fluid - u velocity in axial direction - v velocity in azimuthal direction - w velocity in radial direction - x axial co-ordinate - x escaled axial co-ordinate - X dimensionless axial co-ordinate - thermal diffusivity of the saturated medium - coefficient of thermal expansion - constant in the boundary conditions forF ₀ - dimensionless radial co-ordinate - co-ordinate for the outer zone in then<1 solution - scaled radial co-ordinates - scaled fluid temperature - similarity variable for then=1 problem - nondimensionalisation constant (Eq. (9)) - viscosity of fluid - scaled axial co-ordinates - density of fluid - co-ordinate for the inner zone in then<1 solution - azimuthal co-ordinate - similarity variables for then>1 problem - streamfunction     

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₀     

Die Entwicklung von Längswirbeln in zeitlich anwachsenden Grenzschichten an konkaven Wänden

Zusammenfassung Messungen des Anwachsens von Längswirbeln in zeitlich anwachsenden Grenzschichten an konkav gekrümmten Wänden (Görtler-Taylor-Wirbel) ergaben drei deutlich getrennte Bereiche: Es traten zunächst Wirbel mit der Wellenläge 0,9 auf (=Grenzschichtdicke, =Höhe einer Zelle, die zwei gegensinnig drehende Wirbel enthält). Je nach Größe der mit der Verdrängungsdicke ₁ der Grenzschicht gebildeten Reynolds-Zahl erschienen dann kurze Zeit später Wirbel mit 2,5, wenn war. Im Bereiche dagegen traten stattdessen bei den hier durchgeführten Versuchen immer Wirbel mit der Wellenlänge 6,5 auf. Bei werden die ersten Tollmien-Schlichting-Wellen mit der Wellenlänge _TS 6· angefacht. In ihren wandnahen Bereichen der Wellentäler könnten sich dann die oben genannten Längswirbel der Wellenlänge 6,5· ausbilden, die die zwei-in eine dreidimensionale Störung allseits gleicher Größenordnung verwandeln können.

---

The development of longitudinal vortices in boundary layers growing with time along concave walls

Summary Measurements of the growth of longitudinal vortices in boundary layers growing with time along concave walls (Görtler-Taylor vortices) rendered three distinctly separated regions. First, vortices with a wave-length 0.9 appeared (-boundary layer thicness, =height of a cell containing two counterrotating vortices). Then, depending on the Reynolds number R _{a 1}/v ₁=displacement thickness), vortices with 2.5 appeared shortly afterwards, provided . In the region , however, the wave-length was 6.5. For the first Tollmien-Schlichting waves with _TS 6 were excited. In the wave-throughs close to the wall the abovementioned longitudinal vortices with wave length 6.5 may then be formed. This might transform the two-dimensional into a three-dimensional flow of equal order of magnitude in all directions.

---

Zeichenerklärungen R _a Innenradius - Re _a Reynolds-Zahl gebildet mit dem InnenradiusR _a - Reynolds-Zahl gebildet mit der Verdrängungsdicke ₁ - kritische Taylor-Zahl - h Standhöhe der Flüssigkeit im Zylinder - t Zeit - z Anzahl - Steigungswinkel der Geraden - Grenzschichtdicke - ₁ Verdrängungsdicke - Wellenlänge (enthält ein gegensinnig rotierendes Längswirbelpaar) - v kinematische Zähigkeit - Winkelgeschwindigkeit Indizes K Knickpunkt der Geradensteigung - L unterhalb des Knickpunktes der Geradensteigung - TS Tollmien-Schlichting - e Einsatz der Wirbelentstehung     

First-order chemical reaction on flow past an impulsively started vertical plate with uniform heat and mass flux

Summary A finite-difference solution of the transient natural convection flow of an incompressible viscous fluid past an impulsively started semi-infinite plate with uniform heat and mass flux is presented here, taking into account the homogeneous chemical reaction of first order. The velocity profiles are compared with the available theoretical solution and are found to be in good agreement. The steady-state velocity, temperature and concentration profiles are shown graphically. It is observed that due to the presence of first order chemical reaction the velocity decreases with increasing values of the chemical reaction parameter. The local as well as average skin-friction, Nusselt number and Sherwood number are shown graphically.List of symbols C concentration - C species concentration in the fluid far away from the plate - C _w species concentration near the plate - C dimensionless concentration - D mass diffusion coefficient - Gc mass Grashof number - Gr thermal Grashof number - g acceleration due to gravity - j mass flux per unit area at the plate - K dimensionless chemical reaction parameter - K _l chemical reaction parameter - k thermal conductivity - Nu_x dimensionless local Nusselt number - dimensionless average Nusselt number - Pr Prandtl number - q heat flux per unit area at the plate - Sc Schmidt number - Sh_x dimensionless local Sherwood number - dimensionless average Sherwood number - T temperature - T temperature of the fluid far away from the plate - T _w temperature of the plate - T dimensionless temperature - t time - t dimensionless time - u ₀ velocity of the plate - U, V dimensionless velocity components inX,Y-directions, respectively - u, v velocity components inx, y-directions, respectively - X dimensionless spatial coordinate along the plate - x spatial coordinate along the plate - Y dimensionless spatial coordinate normal to the plate - y spatial coordinate normal to the plate - thermal diffusivity - volumetric coefficient of thermal expansion - * volumetric coefficient of expansion with concentration - coefficient of viscosity - kinematic viscosity - _x dimensionless local skin-friction - dimensionless average skin-friction     

Numerical solution to the solidification of aluminium in the presence of various fibres

In an attempt to understand the experimentally observed solidification microstructures in metal matrix composites, the influence of SiC, graphite and alumina fibres on the solidification of aluminium has been studied numerically. Irregular geometries of the composite material were mapped into simple rectangles through numerical conformal mapping techniques to analyse the influence of a single fibre or a row of fibres on a unidirectionally advancing planar solid-liquid interface. The fibres were assumed to be circular in cross-section and the direction of the interface movement was perpendicular to the length of the fibres. The study showed that for fibres with lower thermal conductivity than aluminium, the interface first goes through acceleration as it approaches and ascends the fibre and then deceleration as it descends the fibre. The acceleration and deceleration phenomena of the interface increases as the thermal conductivity ratio of fibre to liquid aluminium decreases. With low thermal conductivity ratios (K _f/K _L1), the interface is orthogonal to the fibre surface. When the conductivity of the fibre is lower than that of the melt, the interface becomes convex facing the fibre; this mode would lead to pushing of the fibre ahead if it was free to move, as has been experimentally observed in cast microstructures of metal matrix composites. The temperature versus solidification time plots of two points, one in the fibre and the other in aluminium, show that the fibre with a conductivity lower than the matrix is at a temperature higher than the melt; the temperature difference between the two points increases with increasing solidification rate for all the positions of the interface before it touches the fibre. The three-fibre study shows that as the number of fibres increases, the curvature of the interface increases upon approaching the subsequent fibres. The relationship between these numerical computations and experimental observations has been discussed.Nomenclature a reference length = diameter of the fibre - h - K thermal conductivity; in Equation 4 it is defined as K = (K + K _f)/2 for the common boundary between fibre and the freezing medium. For all the rest of the points K = K in Equation 4 - L latent heat of fusion - r non-dimensional variable in radial direction - S non-dimensional distance travelled by the interface - Ste Stefan number = - T non-dimensional temperature - t non-dimensional time - x a non-dimensional spatial coordinate of physical plane - y a non-dimensional spatial coordinate of physical plane - thermal diffusivity - non-dimensional axial coordinate of the mapped plane - non-dimensional vertical coordinate of the mapped plane - a polar coordinate - l liquid - m melting - s solid - O constant wall temperature - i initial - f fibre - * dimensional variables     

Decreasing of the initial shear modulus with increasing axial strain explained by means of a plastic-hypoelastic model

Summary The purpose of this work is to examine in detail the possibility to explain the decreasing of the initial shear modulus with increasing axial strain, observed first by Feigen, by means of the plastic-hypoelastic stress-strain relation suggested by Lehmann and by the author of the present paper.Notations _ij components of the infinitesimal strain tensor dilatation - strain deviator - _ij components of the stress tensor - spherical part of the stress tensor - stress deviator - ²= _{ij ij} second invariant of the stress deviator - = ₃₃ axial strain - e= ₁₃ shear component of the strain tensor - =2 ₁₃ shear strain - = ₃₃ axial stress - s= ₁₃ shear stress - T _ij components of Cauchy's stress tensor - F _ij components of the deformation gradient - L _ij components of the velocity gradient (Eulerian coordinates) - components of the rate of deformation tensor - components of the spin tensor - components of the rate of deformations deviator - components of Cauchy's stress deviator - T=T ₃₃ axial Cauchy's stress With 7 Figures     

Large deformation theory of coupled thermoplasticity including kinematic hardening

Summary Thermoplasticity is a topic central to important applications such as metalforming, ballistics and welding. The current investigation introduces a thermoplastic constitutive model accommodating the difficult issues of finite strain and kinematic hardening. Two potential functions are used. One is interpreted as the Helmholtz free energy. Its reversible portion describes elastic behavior, while its irreversible portion describes kinematic hardening. The second potential function describes dissipative effects and arises directly from the entropy production inequality. It is shown that the dissipation potential can be interpreted as a yield function. With two simplifying assumptions, the formulation leads to a simple energy equation, which is used to derive a rate variational principle. Together with the Principle of Virtual Work in rate form, finite element equations governing coupled thermal and mechanical effects are presented. Using a uniqueness argument, an inequality is derived which is interpreted as a finite strain thermoplastic counterpart to the classical inequality for stability in the small. A simple example is introduced using a von Mises yield function with linear kinematic hardening, linear isotropic hardening and linear thermal softening.Symbols D rate of deformation tensor - d VEC(D) - F deformation gradient tensor - h heat generation per unit mass - L velocity gradient tensor - q heat flux vector - workless internal variable - Lagrangian strain - e VEC() - E quasi-Eulerian strain - entropy - internal energy per unit mass - Helmholtz free energy - Cauchy stress tensor - Truesdell stress flux tensor - t VEC() - yield function - First Piola Kirchhoff stress - Second Piola Kirchhoff stress - s VEC() - s ^* backstress, center of the yield surface - Kronecker product symbol - VEC vectorization operator - tr(.) trace - DEV deviator of a tensor - TEN22 Kronecker tensor operator     

Die Herleitung von Schalengleichungen über ein erweitertes Variationsprinzip bei Berücksichtigung der Querschubverzerrungen

Ohne ZusammenfassungBezeichnungen L Bezugsgrößen für dimensionslose Koordinaten - L charakteristische Schalenabmessung - t Schalendicke - Schalenparameter - körperfeste, krummlinige, dimensionslose Koordinaten der Schalenmittelfläche - Dimensionslose Koordinate in Richtung der Schalennormalen - i, j,...=1,2,3 Indizierung des dreidimensionalen Euklidischen Raumes - ,,...=1,2 Indizierung des zweidimensionalen Riemannschen Raumes - (...), Partielle Differentiation nach der Koordinate - (...), Kovariante Differentiation für Tensorkomponenten des zweidimensionalen Raumes nach der Koordinate - (...)| Kovariante Differentiation für Tensorkomponenten des dreidimensionalen Raumes nach der Koordinate - Variationssymbol - a ,a ₃ Basisvektoren der Schalenmittelfläche - V Verschiebungsvektor - U ,U ₃ Verschiebungskomponenten des Schalenraumes - v ,w,w ,W Verschiebungskomponenten der Schalenmittelfläche - Verhältnis der Metriktensoren des Schalenraumes und der Schalenmittelfläche - _ik Verzerrungstensor des Raumes - _{(, )}, Symmetrische Verzerrungstensoren der Schalenmittelfläche - _{[, ]} Antimetrischer Term des Verzerrungsmaßes - ^, Spannungstensor - n ,m ,q Tensorkomponenten der Schnittgrößenvektoren - p ,p,c Tensorielle Lastkomponenten     

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