Heat and mass transfer effects on flow past an impulsively started vertical plate

Summary An exact solution to the problem of lfow past an impulsively started infinite vertical plate in the presence of uniform heat and mass flux at the plate is presented by the Laplace-transform technique. The velocity, the temperature and the concentration profiles are shown graphically. The rate of heat transfer, the skin-friction, and the Sherwood number are also shown on graphs. The effect of different parameters like Grashof number, mass Grashof number, Prandtl number, and Schmidt number are discussed.List of symbols C species concentration near the plate - C species concentration in the fluid far away from the plate - C dimensionless concentration - C _p specific heat at constant pressure - D mass diffusion coefficient - g acceleration due to gravity - Gr thermal Grashof number - Gc mass Grashof number - j mass flux per unit area at the plate - K thermal conductivity of the fluid - Nu Nusselt number - Pr Prandtl number - q heat flux per unit area at the plate - Sc Schmidt number - t time - t dimensionless time - T temperature of the fluid near the plate - T temperature of the fluid far away from the plate - T _w temperature of the plate - u velocity of the fluid in thex-direction - u ₀ velocity of the plate - u dimensionless velocity - x coordinate axis along the plate - y coordinate axis normal to the plate - y dimensionless coordinate axis normal to the plate - volumetric coefficient of thermal expansion - ^* volumetric coefficient of expansion with concentration - coefficient of viscosity - kinematic viscosity - density - skin-friction - dimensionless skin-friction - dimensionless temperature - er fc complementary error function - similarity parameter     

Microstructure and mechanical properties of rapidly quenched L20 and L20+L12 alloys in Ni-Al-Fe and Ni-Al-Co systems

Ductile L2₀-type wires and+L1₂-type duplex wires with high strengths and large elongation in the Ni-Al-Fe and Ni-Al-Co ternary systems have been manufactured directly from the liquid state by an in-rotating-water spinning method. The wire diameter was in the range 80 to 180m and the average grain size was 2 to 4m for the wires and 0.2 to 1.0m for the+ wires. _y, _f and _p of the wires were found to be about 360 to 760 MPa, 560 to 960 MPa, and 0.2 to 5.5%, respectively, for the Ni-Al-Fe system, those of the+ wires were about 395 to 660 MPa, 670 to 1285 MPa, and 3.5 to 17%, respectively, for the Ni-Al-Fe system, and about 260 to 365 MPa, 600 to 870 MPa, and 4.0 to 7.0%, respectively, for the Ni-Al-Co system. Cold-drawing caused a significant increase in _y and _f and the values attained were about 1850 and 2500 MPa, respectively, for Ni-20Al-30Fe and Ni-25Al-30Co wires drawn to about 90% reduction in area. The high strengths, large elongation and good cold-workability of the melt-quenched and+ compound wires have been inferred to be due to the structural change into a low-degree ordered state containing a high density of phase boundaries, suppression of grain-boundary segregation and refinement of grain size.     

The influence of biaxiality of loading on the form of caustics in cracked plates

Summary The optical method of reflected caustics was applied up-to-now to problems of cracked plates under uniaxial loading. Only the problem of the biaxial tension of the plate has been considered for the particular case where the crack is transverse to the longitudinal axis of the plate which coincided with the loading axis. In this paper the influence of a biaxial loading of the plate on the form and orientation of the caustic was studied in connection with the orientation of the crack. New modified relations were given for the evaluation of the complex stress intensity factorK=K _I –iK _II in terms of the angle of the angular displacement of the caustic axis. For the accurate evaluation ofK _I andK _II nomograms of correction factors _y ^max , _x ^max and _x ^max were given in terms of the angle of inclination of the crack =(90–) and the biaxiality factork. Experimental evidence with PMMA internally cracked plates corroborated the results of theory.List of Symbols (z), (z) complex-stress function of Muskhelishvili - _xx , _yy , _xy crack tip stress referred to Cartesian coordinate system - r, polar coordinate system centered at crack tip - K _I ,K _II stress intensity factors for ModeI andII loading, respectively - angle of inclination of the crack - 90°- - k ratio of stresses at infinity - ₁, ₂ principal stresses at crack tip - a crack length - stress applied at infinity along the transverse boundaries of the plate - X _r,f ,Y _r,f parametric equations of the reflected caustics referred to the Cartesian systemOXY on the reference screen: (r) reflected caustics from rear face of the specimen and (f) reflected caustics from the front face of the specimen - r ₀ radius of the generatrix curve on the specimen around the crack tip (initial curve) - c _r,f optical constants of the material for reflections from the rear and front faces of the specimen respectively - _m magnification ratio of the optical set-up - z ₀ distance between the reference-screen and the middle plane of the specimen - z _i distance between the focus of the light beam and the middle plane of the specimen - d thickness of specimen - 2 for the reflected caustics from the rear face of the specimen and 1 for the reflected caustics from the front face of the specimen - C _r,f z ₀ dc _r,f / _m (2)^1/2 - v Poisson's ratio - E elastic modulus of the material - A (1+k)+(1–k) cos 2 - B (1–k) sin 2 - C 1+k ²+(1–k ²) cos 2 2 tan^–1 (B/A)=2 tan^–1 (K _II/K_I) - D _y ^max ,D _x ^max ,D _x ^min the maximum and the ninimum diameter of caustics along the axisOy andOx of the crack respectively - _y ^max , _x ^max , _x ^min the correction factors forD _y ^max ,D _x ^max andD _x ^min respectively - D _t ^max ,D _l ^max the maximum transverse and longitudinal diameters of the caustics respectively - _t ^max , _l ^max the correction factors forD _t ^max ,D _l ^max respectively With 13 Figures     

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.     

Combined heat and mass transfer in natural convection flow from a vertical wavy surface

Summary In the present paper, effects of combined buoyancy forces from mass and thermal diffusion by natural convection flow from a vertical wavy surface have been investigated using the implicit finite difference method. Here we have focused our attention on the evolution of the surface shear stress,f(0), rate of heat transfer,g(0), and surface concentration gradient,h(0) with effect of different values of the governing parameters, such as the Schmidt number Sc ranging from 7 to 1500 which are appropriate for different species concentration in water (Pr=7.0), the amplitude of the waviness of the surface ranging from 0.0 to 0.4 and the buoyancy parameter,w, ranging from 0.0 to 1.Notation C species concentration in the boundary layer - C species concentration of the ambient fluid - C _w species concentration at the surface - D chemical molecular diffusivity - f dimensionless stream function - g acceleration due to gravity - Gr _x local modified Grashof number - N ratio of the buoyancy forces due to the temperature difference and the concentration difference - p pressure of the fluid - T temperature of the fluid in the boundary layer - T temperature of the ambient fluid - T _w temperature at the surface - u, v thex- andy-components of the velocity field - x, y axis in the direction along and normal to the plate Greek symbols thermal diffusivity - _T volumetric coefficient of thermal expansion - _C volumetric coefficient of expansion with concentration - stream function - nondimensional similarity variable - x/L - density of the ambient fluid - v kinematic coefficient of viscosity - stream function - dimensionless skin friction - fluid viscosity     

A flat jet in a granular bed of finite height

The distribution of gas flows in the vicinity of the jet is discussed and the conditions of disruption of the static equilibrium of the bed, the formation and growth of a cavity, and the jet breakthrough of the bed are investigated qualitatively.Notation a, b functions calculated in [11] - C, C constants in (7) - F derivative of the complex potential - f function in (6) - G function defined in (19) - H dimensionless height of bed - h height of cavity - k coefficient introduced in (15) - p, po pressure inside bed and in cavity - p dimensionless pressure drop - Q, q dimensional and dimensionless jet flow rates - q₁, q₂ critical values - T dimensionless height of cavity - T₀, T₁ T₁, T₂ characteristic values of T - u,v filtration velocities - u, u_* initial filtration velocity in the bed and minimum fluidization velocity - u_o velocity scale introduced in (14) - u * velocity scale introduced in (14) - u* velocity of fictitious flow defined in (15) - U complex velocity - Z=X+iY, z=x+iy dimensionless coordinates - z=x+iy dimensional coordinates - coefficient of hydraulic resistance - parameter from (5) - specific weight of particles' material - porosity - =+i coordinates in the plane obtained from z=x+iy as a result a of conformai transformation - _m value of giving a minimum of the function G - f complex and real flow potentials - angle of internal friction - stream function - angle of inclination of boundaries of the region of plastic flow to the vertical Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 37, No. 5, pp. 804–812, November, 1979.     

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     

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.     

Viscous oscillating cascade aerodynamics and flutter by a locally analytical method

A mathematical model is developed to analyze the viscous aerodynamics of an harmonically oscillating flat plate airfoil cascade in an incompressible laminar flow. The steady flow field is described by the Navier-Stokes equations, with the unsteady viscous flow modeled as a small perturbation to this steady flow. Solutions for both the steady and the unsteady viscous flow fields are then obtained by developing locally analytical solutions. The significant effects of Reynolds number, elastic axis, interblade phase angle and incidence angle on the oscillating cascade unsteady aerodynamics and torsional flutter characteristics are then demonstrated.List of symbols C airfoil chord - C _M unsteady moment coefficient - k reduced frequency, U/U - Re Reynolds number, U C/v - S cascade spacing - U free-stream velocity magnitude - x _ea elastic axis location - x mean flow direction coordinate - y normal flow direction coordinate - y _m mean airfoil position - x x-direction step size - y y-direction step size - (x ₀, y ₀) center of grid element - ₀ mean flow incidence angle - amplitude of airfoil oscillation - interblade phase angle - nondimensional unsteady stream function - nondimensional steady stream function - nondimensional unsteady vorticity - nondimensional steady vorticity - cascade stagger angle     

High-frequency dielectric properties of Mg-Al-Si,Ca-Al-Si and Y-Al-Si oxynitride glasses

Recently developed coaxial line techniques [1] have been used to determine, at room temperature, the values of the real () and imaginary (') parts of the dielectric constants for some Mg-Al-Si, Ca-Al-Si and Y-Al-Si oxynitride glasses over the frequency range 500 MHz to 5 GHz. The frequency dependencies of and ' are consistent with the universal law of dielectric response in that (-_t8)^(n–1) and '^(n–1) for all glass compositions; the high experimental value of the exponent (n=1.0±0.1) suggests the limiting form of lattice loss [2] situation. In this frequency range, as previously reported [3] at longer wavelengths, the addition of nitrogen increases the dielectric constant, (); in both the oxide and oxynitride glasses is also influenced by the cation, being increased with cation type in the order magnesium, yttrium, calcium as at lower frequencies.     

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.     

Methodology for the determination of the interfacial properties of brittle matrix composites

The interfacial properties of a glass-ceramic matrix composite (SiC/CAS) were determined from single-fibre push-out tests using the interfacial test system. The coefficient of friction, , the residual clamping stress, _c, and fibre axial residual stress, _z , were extracted by fitting the experimental stress versus fibre-end displacement curves using the models of Hsueh, and Kerans and Parthasarathy. Using Hsueh's model, the intrinsic interfacial frictional stress (=_c) was found to be 11.1±3.2 MPa, whereas by using Kerans-Parthasarathy's model it was found to be 8.2±1.5 MPa. Comparisons between these models are included, together with a discussion of data analysis techniques.Nomenclature _z Axial fibre residual stress (Pa) - ^* Effective clamping stress (Pa) - _c Residual clamping stress (Pa) - _p Poisson's effect-induced clamping stress (Pa) - _d ⁰ Debond stress in the absence of residual stresses (Pa) - _d Experimental debond stress (Pa) - Compressive applied stress (Pa) - Interfacial shear stress (Pa) - u Fibre-end displacement (m) - h Debond length (m) - r Fibre radius (m) - E _f Fibre Young's modulus (Pa) - E _m Matrix Young's modulus (Pa) - v _f Fibre Poisson's ratio (dimensionless) - v _m Matrix Poisson's ratio (dimensionless) - f Fibre volume fraction (dimensionless) - k Parameter (dimensionless) - D Parameter (dimensionless) - Interfacial coefficient of friction (dimensionless) - G _i Interface toughness (J m^–2) - C _m Load-train compliance (m N^–1)     

Non-Darcy regime mixed convection on vertical plates in saturated porous media with lateral mass flux

Summary A unified treatment is presented of mixed convection on vertical plates embedded in fluid saturated porous media with prescribed variable plate temperature or surface heat flux for the case of non-Darcy limiting regime. The plates are permeable with lateral mass flux. By suitable similarity transformations, it is shown that the two problems of prescribed temperature and prescribed heat flux lead to identical differential equations with two common boundary conditions and third boundary condition differing in the two cases. The effect of lateral mass flux and the free stream on the Nusselt number and the energy transport by the boundary layer is investigated. The unified approach to the mixed convection problem includes free convection problem as a special case. Exact analytical solutions are obtained for two cases of free convection problem.Notation b inertial coefficient - c specific heat - D _p pore diameter - E rate of upward energy transport - Ê dimensionless rate of energy transport - f dimensionless stream function - f _w mass flux parameter - g acceleration due to gravity - k thermal conductivity - K permeability of the porous medium - m exponent in the variation of heat flux - M mixed convection parameter - Nu _x Nusselt number - Pe _x Peclet number - q _w surface heat flux - Ra _x local Rayleigh number - Ra _x ^* modified local Rayleigh number - T temperature - T _e ambient temperature - T _w plate temperature - T _w temperature difference=T _w-T_e - u velocity in thex-direction - u _e free stream velocity - v velocity in they-direction - v _w lateral velocity at the plate - x coordinate along the plate in the upward direction - y coordinate normal to the plate - equivalent thermal diffusivity =k/c _e - coefficient of thermal expansion - porosity - dimensionless variable - dimensionless temperature - viscosity - fluid density - _e ambient density - exponent in the variation of plate temperature - stream function     

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     

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     

The dielectric properties of glass-ceramic-on-metal substrates for microelectronics packaging

The dielectric properties of some glass-ceramic-on-metal substrates have been determined over the frequency range 500 Hz to 5 MHz using a.c. bridge techniques. The substrates consisted of cordierite-based glass-ceramics screen printed on molybdenum. For glass layers of thickness greater than 100 m both the permittivity, and the dielectric loss, , are frequency independent over this frequency range at room temperature giving the value of =6.5 and tan =8×10^–3; the room-temperature data are consistent with the universal law of dielectric response. The variation of permittivity with temperature has also been examined and, below 120 °C, the temperature coefficient [(–1) (+2)]^–1 (/T)_p, was found to be 1.3×10^–5 K^–1. The results are compared with those previously reported for Al₂O₃ and AIN substrates.     

Precipitate strengthening mechanisms in magnesium zinc alloy single crystals

The mechanical behaviour of Mg-5.1 wt % Zn alloy single crystals was studied in the 4.2 to 300° K temperature range. Quenched crystals have activation energies and volumes best associated with the cutting of small clusters of Zn atoms by dislocations. Fully hardened crystals contain fine ₁ and occasional ₂ precipitates with an average ₁ interparticle spacing of 330 to 660 Å. Strengthening in these crystals is mainly ascribed to the cutting of ₁ particles by dislocations. In the overaged condition ₁, ₂ and equilibrium particles are present and lead to a considerable temperature-dependence unusual for an overaged condition. Analysis of this temperature-dependence suggests that below 77° K the relatively easy cutting of ₁ particles by dislocations takes place in addition to the cutting of ₂ and particles. Above 77° K the difficult cutting of ₂ and particles alone controls the deformation, ₁ being more easily cut with the aid of thermal fluctuations.     

Temperature field in a plate containing a system of heat sources

The temperature field is determined in a circular plate with a system of thin extrinsic heat sources.Notation T temperature in the plate with the inclusions - r polar radius - polar angle - time - (r,) coefficient of thermal conductivity - (r,) heat transfer coefficient - C(r,) volume heat capacity - W(r,, ) specific intensity of the heat sources - half thickness of the plate - (x) Dirac's delta function - ¯T finite Fourier cosine transform of the temperature - p parameter for this transformation - T Laplace transform of the temperature - s its parameter - Iv(x) Bessel function with imaginary argument of order - K _v (x) the MacDonald function of order - and dimensionless temperature - Po Pomerantz number - Bi Biot number - Fo Fourier's number - dimensionless polar radius - b₁ ^* dimensionless radius of the circle on which the inclusions are placed - R^* dimensionless radius of the plate Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 40, No. 3, pp. 495–502, March, 1981.     

Heat transfer in a “pseudoturbulent” bed of dispersed material

A two-phase model is proposed for the steady heat exchange between a surface and a pseudoturbulent bed of dispersed material. Expressions are obtained for the temperature fields of the gaseous and solid phases.Notation _g effective thermal conductivity of gaseous phase - _s effective thermal conductivity of the mixed solid phase - porosity - _m molecular thermal conductivity - d particle diameter - temperature of dispersed bed at a large distance from heat source - , g gas temperature - _p particle temperature - _w wall temperature - x current coordinate in the direction perpendicular to the wall - l bed thickness - q heat flux - coefficient of heat exchange between wall and pseudoturbulent bed of dispersed material - ^* coefficient of interphase heat exchange - _g=_g/_w dimensionless gas temperature - _p = _p/_w dimensionless particle temperature - Y = x/d dimensionless coordinate - L =l/d dimensionless bed thickness - A_h dimensionless coefficient of interphase heat exchange - Nu_g = d/_s Nusselt number Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 41, No. 3, pp. 465–469, September, 1981.     

A gradient method of defining smooth solutions to inverse boundary-value problems in thermal conduction

An iterative algorithm is described for solving boundary-value inverse problems in thermal conduction by steepest descent, which utilizes information on the smoothness of the solution.Notation A, B linear operators - u element of solution space U - f exact reference data - f reference data uncertainty - value of reference data uncertainty - A^–1 inverse operator - u^(k)() k-th derivative of function u - _m length of observation interval - _i(t) polynomials of degree i–1 - A^*, B^*, L^* operators conjugate to the operators A, B, L - Jg discrepancy functional gradient - _n descent step along the discrepancy antigradient for the n-th iteration - K( –) kernel of integral equation - q() heat flux - T() measured temperature inside body Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 39, No. 2, pp. 259–263, August, 1980.