Use of adjoint functions in investigations of heat conduction and transfer processes

An examination is made of the use of adjoint functions in heat conduction and convection theory. Formulas of perturbation theory are obtained for steady and unsteady cases, an interpretation of the physical meaning of adjoint temperature is given, and some applications of the theory are discussed.Notation (r,) thermal conductivity - t(r,) temperature - t *(r,) adjoint temperature - q_V(r,) density of heat release sources - p(r,) a parameter of adjoint equation - r generalized coordinate - time - (r_s, ) heat transfer coefficient - I linear functional of temperature - (r,;r₀,₀) and *(r,; r₀,₀) Green's function for t(r, ) and t *(r, ) - C(r,) volume specific heat - W(r, ) vector distribution of flow velocities - V, S volume and surface areas of body - R radius of HRE - r, radial and angular coordinates - F_in, F_out inlet and outlet flow areas of channel     

Motion of a thermal wave front in a nonlinear medium with absorption

We study the evolution of a thermal perturbation in a nonlinear medium whose thermal conductivity depends on the temperature and the temperature gradient according to a power law.Notation u temperature - k coefficient of thermal conductivity - t time - x spatial variable - x₊ a point on the thermal wave front - a ² generalized coefficient of thermal diffusivity - , , , and s parameters of the process - (x^s) Dirac delta-function - B[, ] a beta function - v(, x), (t) auxiliary functions - A, C, T_o, T_m, T*, R, r, p, and _m constants and parameters Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 39, No. 4, pp. 728–731, October, 1980.     

Unsteady motion of a fluid near a disk rotating in a magnetic field

The effect of a magnetic field on the velocity distribution in a fluid close to an unsteadily rotating disk is investigated.Notation r, , and z coordinates in the radial, circular, and axial directions - t time - u, v, and w radial, circular, and axial velocity components - u₀ radial velocity of external potential flux - v₀ circular velocity of the disk - (t) angular velocity of the disk - p pressure - density - v kinematic viscosity - B₀ characteristic of the applied magnetic field - electrical conductivity of fluid - R and Z dimensionless coordinates in the radial and axial directions - =Z/2 dimensionless coordinate - T dimensionless time - U, V, and W radial, circular, and axial components of dimensionless velocity - P dimensionless pressure - a, , and ₀ constants with dimensionality t^–1 - m, n, and positive numbers - k =a constant - = = B ₀ ² / parameter characterizing the magnetic field     

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     

The morphology of α′ martensite in a two-phase (α+γ) Fe-Cr-Ni stainless steel

Both lath-shaped and martensites are induced by tensile deformation within of a two-phase ( + ) Fe-Cr-Ni stainless steel, forms from the through the at an intersection of two crystals. These are observed both when is surrounded by and when borders . The amount of strain at which both and nucleate, increases with test temperatures in the range –196 to 50° C. Adjacent laths are either twin-related, or 5, 9, 15 or 19° off the twin relationship, as found by analysing electron diffraction patterns.     

Temperature and randomness dependence of the pair-breaking parameter in superconducting aluminum films

Magnetoconductance and excess conductance due to superconducting fluctuations in aluminum films are measured in order to study the temperature dependence of the pair-breaking parameter at temperatures nearT _c . The parameter _M is estimated from the relation =/8k _B T_in, where _in is the inelastic scattering time deduced from the analysis of the magnetoconductance. The parameter _F is determined by fitting theories to data on the excess conductance at zero magnetic field. It is shown that: (1) For films with a wide range of the sheet resistanceR , 12R 200 /, the temperature dependence of _M nearT _c agrees well with the theory of Brenig et al. (2) For clean films withR 100 /, the value of _F analyzed with theories including the correction term to the Maki-Thompson contribution shows almost the same temperature dependence as _M . In a film withR 200 /, however, a discrepancy between _M and _F remains.On leave from College of General Education, Kyushu University, Ropponmatsu, Fukuoka, Japan.     

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     

Free convective flow in an enclosure with a cooled inclined upper surface

Natural convective air flow in an enclosure with a horizontal lower wall, vertical side-walls and a straight inclined top wall has been numerically studied. The lower wall is at a uniform high temperature while the top wall is cooled to a uniform lower temperature. The temperature of the side-walls varies in a prescribed way between the temperatures of the bottom and top walls. This flow situation is related to that occurring in a proposed system for use in developing countries for drying crops such as corn and rice. It has been assumed that the flow is steady, laminar, and two-dimensional and that the fluid properties are constant except for the density change with temperature that gives rise to the buoyancy forces. The governing equations have been expressed in terms of stream function and vorticity and written in dimensionless form. The finite element method has been used to obtain the solution to these dimensionless equations. Results have been obtained for enclosures with aspect ratios of between 0.25 and 1 for top surface angles of inclination of between 0 and 45° and for Rayleigh numbers, based on the enclosure width, of between 1000 and 10⁷. The results have been used to study the effect of changes in the governing parameters on the flow pattern in the enclosure and on the mean heat transfer rate to the upper surface.List of symbols A aspect ratio, h/w - h mean height of cavity - k thermal conductivity - Nu local Nusselt number based on w - mean Nusselt number for top surface based on w - n n/w - n coordinate measured normal to boundary surface considered - p pressure - Pr Prandtl number - q local dimensionless heat transfer rate - mean heat transfer rate to top surface - Ra Rayleigh number based on w - T dimensionless temperature - T temperature - T _H temperature of bottom surface - T _C temperature of top surface - u velocity component in x direction - v velocity component in y direction - w width of cavity - x dimensionless x coordinate - x horizontal coordinate position - y dimensionless y coordinate - y vertical coordinate position - thermal diffusivity This work was supported by the Natural Sciences and Engineering Research Council (NSERC) of Canada.     

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₀     

Solidification behavior of Al-Si-Fe alloys and phase transformation of metastable intermetallic compound by heat treatment

Solidification behavior of Al-20 wt% Si-8 wt% Fe and Al-30 wt% Si-5 wt% Fe alloys during cooling with a cooling rate of 10 K/min has been studied using optical microscopy, X-ray diffractometry, differential thermal analysis, scanning electron microscopy and energy dispersive X-ray spectroscopy. In Al-20Si-8Fe alloy, metastable -Al₄FeSi₂ phase with tetragonal structure formed first from melt, followed by primary Si precipitation and then remaining liquid solidified finally into ternary eutectic of -Al, Si and phases. However, in Al-30Si-5Fe alloy, primary Si formed first, followed by the phase precipitation and then eutectic solidification. During isothermal heat treatment of as-solidified alloys, phase transformation from the phase to equilibrium phase began at the interface between phase and -Al matrix and progressed toward the inside of phase with co-precipitation of Si particles due to the difference in composition between -Al₄FeSi₂ and -Al₅FeSi phases.     

Effect of unsteady natural convection on the structure of a liquid flow in a horizontal mixing chamber

Results are presented of a numerical and experimental investigation of the effect of natural convection on the structure of a liquid flow in a horizontal mixing chamber with changes in the temperature of the liquid at the inlet.Notation t_o, t_in initial temperature and temperature at the inlet to the channel - dimensionless temperature - a heat conductivity - kinematic viscosity - coefficient of cubical expansion - density - P pressure - g acceleration due to gravity - d_equ equivalent diameter of porous body - time - v_o mean velocity at the inlet - X and dimensionless vertical and radial coordinates - U and V dimensionless vertical and horizontal components of velocity - H dimensionless height of channel - R radius of inlet to channel - W velocity of liquid - t temperature drop along channel height - Re = v_oR/ Reynolds number - Pe = v_oR/a Peclet number - Fr = v _o ² /g¦t_in–t_o R Froude number - Ho = v_o/R homochroneous number Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 39, No. 4, pp. 603–610, October, 1980.     

γ ? α2 Phase transformation in fractured high temperature stress rupture Ti-48Al-2Nb(at.%)

The ₂ phase transformation in fractured high temperature stress rupture Ti-48Al-2Nb(at.%) alloy has been studied by analytical electron microscopy. ₂ and phases were found at grain boundaries. ₂ layers that suspended in layers and interfacial ledge higher than 2d ₍₁₁₁₎ at /₂ interfaces were observed in the lamellar grains. These facts indicated that ₂ phase transformation and dynamic recrystallization have occurred during high temperature stress rupture deformation. It can be concluded that deformation induced ₂ phase transformation and dynamic recrystallization resulted in the presence of particles at grain boundaries. A structural and compositional transition area between deformation-induced ₂(or ) and its adjacently original (or ₂) phases was found by HREM and EDS and is suggested as a way to transform between and ₂ phase during high temperature stress rupture deformation. The transition area was formed by slide of partial dislocations on close-packed planes and diffusion of atoms.     

Temperature dependence of electrical resistivity of deformed and undeformed V-rich V3Si single crystals

The electrical resistivity (T) of V-rich V₃Si single crystals (T _c-11.4 K) was measured from 4.2 to 300 K along the directions of [1 0 0] and [1 1 1] before and after plastic deformation at 1573 K. Anisotropy of (T) was observed although V₃Si has the cubic A15 structure. Plastic deformation does not affect the normal-state (T) behaviour but changes the normal-superconducting transition width T_c. At low temperatures (T _c<T 40 K), (T) varies approximately as T ⁿ where n-2.5 and this behaviour does not contradict the (0)- phase-diagram plot proposed by Gurvitch, where is the electron-phonon coupling constant and (0) is the residual resistivity.     

Microstructural study of ε-iron-carbonitrides formed by a glow discharge technique

The microstructure and phase transformations occurring in-iron-carbonitrides have been studied by means of X-ray and electron diffraction, electron microscopy and Mössbauer spectrometry. Ordering of the interstitial atoms, N or C, results in a hexagonal unit cell for Fe₃(C, N) with parametersa=a3 andc=c wherea andc are the lattice parameters of the hexagonal close-packed (h c p) iron unit cell. Stacking faults on (0001) planes and partial dislocations with Burgers vector b=1/31 0 ¯1 0 are observed in quenched-Fe₄ (C, N). After quench-aging, the carbonitrides show a structural hardening due to the precipitation of a metastable phase. Slow cooling of-carbonitrides with less than 25 at% interstitials leads to the precipitation of-carbonitride and ferrite in-phase grains which allows the orientation relationships between the-,- and-phases to be defined and a model of the-phase--phase transformation to be proposed.     

Dielectric properties of thin solid films formed on silicon

We report on the measurement of the frequency-dependent complex permittivity, ()=()-i(), over the frequency range, 30 MHz to 6 GHz, of silicon wafers and of thin dielectric films formed on silicon. Measurements, as a function of temperature and time treatments, were obtained by means of an HP Network Analyzer and dielectric probe and the resulting ()and()plots for the silicon wafers are shown to have a Debye-type [1] profile, thereby indicating that the associated polarization mechanism is of the orientational variety.     

Effect of the weight of the admixture on the turbulence structure of a two-phase jet

The effect of gravity on the turbulence structure of an inclined two-phase jet is evaluated according to the Prandtl theory of mixing length.Notation C_x drag coefficient for a particle - D_p particle diameter - g_i components of the acceleration g due to gravity acting on a particle in the direction of jet flow (g_i=g sin ) and in the direction normal to it (g_i=g cos ) - V_poi ^±, V_goi ^± fluctuation components of the velocities of the particles and gas, respectively, at the end of a mole formation - V_fi free-fall velocity of a particle - l _u mixing length - m_p particle mass - t _p length of time of particle-mole interaction - V_pi ^±, V_gi ^± positive and negative fluctuation velocities of particles and of the gas respectively, with the components u_p ^±, u_g ^±, v_p ^±, v_g ^±, k=V_goi/V_fi - V_i ^± relative velocity of the gas - jet inclination angle relative to the earth's surface - empirical constant - _u, jet boundaries in terms of velocity and concentration - _u=y/ _u dimensionless velocity ordinate - =y/ dimensionless concentration ordinate - admixture concentration - u_m, _m velocity and the concentration of the admixture at the jet axis - _g dynamic viscosity of the gas - _s, _g densities of the particle material and of the gas - _g, _p shearing stresses in the gas and in the gas of particles - _m, ₀ shearing stresses in the mixture and in pure gas, respectively Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 40, No. 3, pp. 422–426, March, 1981.     

Pulse determination of thermal diffusivity and thermal conductivity for hemispherical specimens: nickel

A method is described for determining the thermal diffusivity, specific heat, and thermal conductivity in a hemispherical volume on the basis of duration of the reference signal.Notation r radius - R radius - r dimensionless coordinate - dimensionless temperature - time - _i duration of heat pulse - _1/2 time for temperature signal at r to attain half the maximum value - q_o amount of heat - a thermal diffusivity - thermal conductivity - density - c_p heat capacity Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 40, No. 5, pp. 864–869, May, 1981.     

Heat exchange in two-phase flow about a surface located in a rectangular channel

Similitude equations are obtained on the basis of the principle of superposition of separate effects to calculate heat exchange between surfaces with complexshaped cross sections located in a rectangular channel during their cooling by a two-phase flow.Notation T, q temperature and heat flux - T_w mean surface temperature - I, R current and electrical resistance - V volume of the material - , , anda heat-transfer coefficients, thermal conductivity, and linear expansion of the material - relative functions - =_m; ^* = _m ^{* *} _s ^* ; temperature factor - X relative weight content of liquid phase Indices w surface - f incoming flow - v volume - m two-phase flow - angle of attack - s shape of surface - * pertains to surface with swirl vanes Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 40, No. 5, pp. 780–786, May, 1980.     

Transient heat conduction in hollow spheres with a moving inner boundary

The finite integral transform method is used to obtain the solution of unsteady heat conduction problems for a hollow sphere with a moving internal boundary and various boundary conditions at the outer surface. For the solution of the problems of interest integral transform formulas are presented with kernels (16), (20), and (24) and the corresponding inversion formulas (18), (22), (26), (29) and characteristic equations (17), (21), (25), (28), (31), (33).Nomenclature a, thermal diffusivity and conductivity - t temperature of phase transformation - density - heat transfer coefficient - Q total quantity of heat passing through inner boundary - F latent heat of phase transformation - Fo(1,)=a/R ₁ ² , Fo(i,)=/r _i ² , Fo(i, _i)=a _i/r _i ² Fourier numbers - Bi₂=R₂/ Biot number     

Influence lines for bending under a ring load of a free shell,a shell embedded in a soft medium and a shell containing a soft core

Summary Theoretical expressions for stresses and displacements have been derived for bending under a ring load of a free shell, a shell embedded in a soft medium, and a shell containing a soft core. Numerical work has been done for typical cases with anElliot 803 Digital Computer and influence lines are drawn therefrom.

---

Einflußlinien für die Biegung einer freien Schale, einer Schale in einer weichen Bettung und einer Schale mit weichem Kern

Zusammenfassung Für die Biegung einer freien Schale, einer weich gebetteten Schale und einer Schale mit weichem Kern unter einer Ringlast werden Ausdrücke für die Spannungen und Verschiebungen hergeleitet. Die Ergebnisse wurden für einige typische Fälle mit einem DigitalrechnerElliot 803 numerisch ausgewertet. Die sich ergebenden Einflußlinien wurden graphisch dargestellt.

---

Nomenclature A (),B () Functions of - a, t Mean radius and thickness of the shell - E _S , _S Young's modulus andPoisson's ratio of the shell - G _c , _c Shear modulus andPoisson's ratio of the casting or core - I ₀ (r),I ₁ (r) ModifiedBessel functions of the first kind and order zero and one respectively - K ₀ (r),K ₁ (r) ModifiedBessel functions of the second kind and order zero and one respectively - p Ring load, lb/in - U, W Displacement components in the casing or core in thez andr direction - u, w Displacement components of a middle surface point in the shell - _r , _rz Radial and shearing stress components - Independent variable of infinite integrals - k [3(1–S ²)a ²/t ²]^1/4 With 13 Figures