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     

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.     

Thermal diffusivity of inhomogeneous systems

The possibility of analyzing the nonsteady temperature fields of inhomogeneous systems using the quasi-homogeneous-body model is investigated.Notation t, t_I, t_i temperature of quasi-homogeneous body inhomogeneous system, and i-th component of system - a, , c thermal diffusivity and conductivity and volume specific heat of quasi-homogeneous body - a_{i i}, c_i same quantities for the i-th component - q heat flux - S, V system surface and volume - x, y coordinates - macrodimension of system - dimensionless temperature Fo=a/² - Bi=/ Fourier and Biot numbers - N number of plates - =h/ ratio of micro- and macrodimensions - _V, volumeaveraged and mean-square error of dimensionless-temperature determination - time - m_i i-th component concentration Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 39, No. 1, pp. 126–133, July, 1980.     

Structure of thermal stresses arising in a viscoelastic half-space with thermal “memory”

We develop the structure of thermal stresses arising in a viscoelastic half-space owing to the thermal impact of a heat flux at the boundary.Notation z coordinate normal to the surface of the half-space - t time - =T–T_o temperature of the half-space - (t) relaxation function of the heat flux - (t) relaxation function of the internal energy - c_v specific heat at constant volume - q_o heat flux acting on the boundary of the half-space - _xx, _yy, _zz normal stresses - density of the material - _t coefficient of linear thermal expansion - u=[(0)/c_v]^1/2 heat-propagation velocity - _t coefficient of thermal conductivity - _r relaxation time of the heat flux - relaxation time of the stresses Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 37, No. 5, pp. 894–897, November, 1979.     

Magnetic field dependence of zero-sound attenuation close to the pair-breaking edge in3He-B due toJ=1−,J z =±1 collective modes

Our previous theory yielded for the Zeeman splitting of the imaginaryJ=1 collective mode in³He-B the result =2+0.25J _z ( is the effective Larmor frequency). In this paper we take into account the downward shift of the pair-breaking edge from 2 to 2₂– (₂ and ₁ are the longitudinal and transverse gap parameters). This leads to a complex Landé factor: the frequencies of theJ _z =±1 components become =2+0.39J _z , and the linewidths of these resonances become finite: =0.18. The coupling amplitudes of theJ _z =±1 components to density are found to be proportional to gap distortion, (₁–₂/(/)². Our results for the ultrasonic attenuation due to theJ _z =±1,J=1 modes are capable of explaining the field dependence of the attenuation close to the pair-breaking edge as observed by Dobbs, Saunders, et al. The observed peak is caused by theJ _z =–1 component: its height increases due to gap distortion as the field is increased, and the peak shifts downward in temperature and its width increases with the field due to the complex Landé factor. TheJ _z =+1 component gives rise to a corresponding dip relative to the continuum attenuation.     

Numerical method for calculating the filling of a mold and heat transfer of a melt under the action of centrifugal forces

This paper presents an explicit difference method for solving the conjugate problem of pouring molten metal into a casting mold and its solidification under the action of centrifugal forces with allowance for the free surface.Notation r, z transverse and longitudinal coordinates - u, v horizontal and vertical velocity components - V velocity vector - angular velocity of rotation - P pressure - P normalized pressure - T temperature - time - g free fall acceleration - coefficient of volumetric expansion - C heat capacity - thermal conductivity - density - L crystallization heat - v viscosity - Re=|V|h/v Reynolds grid number - h grid spacing - l mixing length in a turbulent flow - relaxation parameter - coefficient of convective heat transfer - coefficient defining the boundary conditions at the solid wall - D flow divergence - volumetric velocity of filling - emissivity - ₀ Stefan-Boltzmann constant - thickness of a layer Indices r, z, I, j numbers of grid nodes - n number of the integration step with respect to time - L, S temperatures of the liquidus and solidus, respectively - s temperature on the surface - med temperature of the medium - 0 initial state of the system - m metal - mol molten state of the metal - red reduced emissivity in the gap - rad radiant component of the heat transfer coefficient; g, gas-air gap - coat heat-insulating coating - fil filling Dneprodzerzhinsk Industrial Institute. Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 68, No. 4, pp. 678–686, July–August, 1995.     

Analysis of the effect of the decomposition of char-forming heat-insulating materials on the free discharge of a gas mixture

The article presents results of a numerical solution of a nonsteady problem on the free discharge of a mixture of gases from a hemispherical volume with allowance for thermal decomposition of heat-insulating materials.Notation V volume - S area - t - P p - T - u v - Q q, dimensional and dimensionless time, pressure, temperature, TIM decomposition rate, and heat flux - adiabatic exponent - R gas constant - density - H specific enthalpy - c specific heat - thermal conductivity - , , _s dimensionless complexes - coefficient expressing the radiative properties of the gas medium and the heat-transfer surface - Stefan-Boltzmann constant Indices 0 initial state and scale factors - s surface - coke - M TIM material - P pyrolysis front - A ablation front - v volatile degradation products - adiabatic conditions - c completion of discharge Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 54, No. 5, pp. 787–793, May, 1988.     

Determination of the thermophysical characteristics of alternatively freezing and thawing dispersed media by the method of solving inverse problems of heat conduction

The article explains an algorithm for determining the thermophysical characteristics of dispersed media with phase transitions based on the method of solving inverse problems of heat conduction.Notation r space coordinate - time - T temperature of the specimen - T₀ initial temperature - c_i, c_w, c_sk specific heat of ice, water, and of the organic-mineral skeleton, respectively - c_f, c_m, _f, _m specific heat and thermal conductivity in the frozen and melted zones, respectively - c effective heat capacity - thermal conductivity - p density - ₀, _sb bound and strongly bound moisture, respectively - (T) amount of nonfrozen water - R radius of the cylinder - q() heat flux - I functional - u₁(), U₂() measured temperatures of the specimen at the points r = 0 and r = R, respectively, at the instant - ₁, ₂ degree of confidence of the supplementary information - final instant of time - a, b, k, _s positive constants - L specific heat of melting - N number of grid nodes over space - n number of grid nodes over time - h grid step over space - grid step over time - solution of the conjugate system - s number of iteration Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 39, No. 2, pp. 292–297, August, 1980.     

Negatively isotone optimal policies for random walk type Markov decision processes

Summary This paper considers a random walk type Markov decision process in which the state spaceI is an integer subset of IR ^m , and the action spaceK is independent ofi I. The natural order, overI, and a quasi order,, overK, is assumed, together with aconditional convexity assumption on the returns {r _i ^k }, and certain other assumptions about these rewards and the transition probabilities in relationship to the orders and.A negatively isotone policy is one for whichi i(i))(i) (i.e.(i) (i) or(i) i)). It is shown that, under specified conditions, a negatively isotone optimal policy exists. Some consideration is given to computational implications in particular relationship to Howard's policy space method.

---

Zusammenfassung Wir betrachten einen Markovschen Entscheidungsprozeß vom random walk Typ. Der ZustandsraumI sei eine Teilmenge des IR^m, wobeii I ganzzahlige Komponenten habe. Die MengeK der zulässigen Aktionen ini I sei unabhängig voni I. Sei die natürliche Ordnung aufI und sei eine Quasiordnung aufK. Die Erträge {r _i ^k }seienbedingt konvex, darüberhinaus seien weitere Voraussetzungen über diese Erträge und die Übergangswahrscheinlichkeiten in Bezug auf die Ordnungen und erfüllt. Eine Politik heißt negativ isoton, falls ausi i folgti(i) (d. h.(i) (i) oder(i)(i)). Wir zeigen, daß unter gewissen Voraussetzungen einenegativ isotone optimale Politik existiert: Auch diskutieren wir einige Folgerungen für die Numerik, insbesondere hinsichtlich Howards Politikiteration.

     

Numerical calculation of a generalized thermal model of radio-electronic apparatus

A method is proposed for numerical calculation of the temperature field of a generalized model of electronic equipment with high component density.Notation x,y,z,x,y spatial coordinates, m - time, sec - L_x, L_v, L_z dimensions of heated zone, m - _x, _y, _z effective thermal-conductivity coefficients of heated zone, W/m·deg - ₂ thermal conductivity of chassis, W/m·deg - a _z thermal diffusivity of heated zone along z axis, m²/sec - c₁ effective specific heat of heated zone, J/kg·deg - ₁ effective density of heated zone, kg/m³ - c₃, ₃, c₂, ₂ thermophysical characteristics of cooling agent and chassis, J/kg·deg·kg/m³ - q_v(x, ), q(x, y) volume heat-source distribution, W/m³ - q_s (x) surface heat-source distribution, W/m² - p number of cooling agent channels - Fo Fourier number - Bi Biot number - Uⁱ coolant velocity in i-th channel, m/sec - T₁(x, ), T₂(x, ), T₃(x, ) temperature distribution of heated zone, chassis, and coolant, °K - T₃₀, T₁₀(x), T₂₀(x) initial temperatures, °K - T_3in coolant temperature at input to channel, °K - T_T(x) effective temperature distribution of heat loss elements, °K - T_C temperature of external medium, °K - dimensionless heated zone temperature - _v(x) local volume heat exchange coefficient, W/m³·deg - ₁₂(x), _1C(x), _1T(x) heat liberation coefficients - W/m²·sec; ₂₁(x, y), _2c(x, y), _2T(x, y) volume heat-exchange coefficients of chassis with heated zone, medium, and cooling elements, W/m³·deg Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 40, No. 5, pp. 876–882, May, 1981.     

Numerical method for solving heat-conduction problems for two-dimensional bodies of complex shape

A finite-difference scheme is described for a curvilinear orthogonal net which permits the use of a single algorithm for calculating bodies of various shapes.Notation x, y independent variables - u, v orthogonal coordinates - F(w)=F(u + iv) function of a complex variable - g(u,v)= F(w)/w Jacobian of transformation from (u,v) to (x,y) - thermal conductivity - c volumetric heat capacity - Q heat release per unit volume - T temperature - f value of temperature on boundary of region - time - L, L₁, L₂ differential operators - (u,v) solution of differential problem in canonical region - ^j, ₁ ^j , ₂ ^j , t^J, t ₁ ^j , t ₂ ^j network functions in canonical region - ^j, t^*j solutions of difference problems using rectangular and orthogonal nets respectively - {u_i, v_k} rectangular net in canonical region G - {x_i,k, y_i, k} orthogonal net in given region G^* - u_i, v_k dimensions of cell of rectangular net - u_i,v _i,k dimensions of cell of orthogonal net - h, maximum dimension of cell for rectangular and orthogonal nets respectively - ₁, ₂, difference operators for rectangular and orthogonal nets - A, B, C, D, A^*, B^*, C^*, D^* coefficients of difference scheme for rectangular net - D, Ã, B coefficients of difference scheme for orthogonal net Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 40, No. 3, pp. 503–509, March, 1981.     

Relationship between the energy supply parameters and the physicochemical l characteristics of polymer compositions during drying and thermal processing

The effect of the type of energy supply on the formation of temperature and concentration fields in the thermal processing of polymer compositions is considered.Notation T₀, T initial and current temperature of the coating - T_m temperature of the air - =(T-T_o)/(T_m-T₀) dimensionless temperature of the coating - a thermal diffusivity - A absorption power of the coating - D diffusion coefficient - thermal conductivity - c thermal capacity - density - _k convective heat transfer coefficient - i number of moles of reacting groups per unit volume of polymer - K₀ factor in front of the exponential - R gas constant - u concentration - Q thermal effect of the reaction - q_n density of the incident radiant flux - =x/ dimensionless coordinate over the thickness of the coating - Ki=Aq_n /(T_m-T₀) Kirpichev criterion characterizing the thermal effect of the reaction - Ki_p=Qi/c (T_m-T₀) analog of the Predvoditelev criterion, characterizing the rate of occurrence of a chemical excess in the system - Bu= Bouguer criterion - Lu=D/a Lykov number - Fo=a/² Fourier number - Bi= _k Biot number Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 39, No. 1, pp. 26–33, July, 1980.     

Theory of RF SQUIDs Operating in the Presence of Large Thermal Fluctuations

A comprehensive analytical theory is presented for non-hysteretic RF SQUIDs operating in the adiabatic mode in the presence of large thermal fluctuations. When 1 ( = 2LI_c/₀ is the hysteresis parameter, L is the SQUID inductance, I_c is the critical current of the Josephson junction, and ₀ is the flux quantum) the theory is applicable also for RF SQUIDs operating in the non-adiabatic mode. In contrast to previous theories in which the noise is treated perturbatively and which therefore are applicable only if the product 1 ( = 2k_BT/ ₀ I_c is the noise parameter, k_B is the Boltzmann constant, and T is the absolute temperature)—the case of small thermal fluctuations—the present theory is valid for around unity or higher. In the limit 0 the theory reproduces the results of small thermal fluctuations theories. It has been found that in the presence of large thermal fluctuations the screening current in the SQUID inductance is suppressed by a factor that increases with increasing . Taking into account this new basic fact, all SQUID characteristics (output signal, transfer function, noise spectral density and energy sensitivity) have been recalculated and a good agreement with experimental data has been obtained. It has been also found that RF SQUIDs can be operated with substantially higher values of the inductance and of the noise parameter than DC SQUIDs. These two aspects, which are of particular importance at liquid nitrogen temperature, make high T_c RF SQUIDs very attractive.     

The piezoresistance coefficients of copper and copper-nickel alloys

This paper attempts to further a better understanding of the piezoresistance coefficients by studying the piezoresistive effects in copper and copper-nickel alloys. The experimental evidence of isotropic piezoresistance coefficients (₁₁=₁₂) has been obtained for the annealed copper and copper-nickel alloys. The piezoresistance coefficients of the cold-worked copper and Cu₆₀Ni₄₀ alloy are of the tensor character (₁₁₁₂). A physical explanation has been given to the change of the ( _ij ) tensor.     

On a simple wave approximation of a set of linear dispersive wave equations

Summary The validity of an approximation ₀ of one of the solutions of a set of two linear coupled dispersive wave equations has been discussed. ₀ is the solution of a linear Korteweg-de Vries equation and satisfies the same initial condition as . It is shown that for square integrable solutions having a spectral range not exceeding [–, ] the approximation is useful if ^{5 2}t«1 in the sense that –₀(t)« (t)(L ₂ -norm). is a measure for the dispersion. The approximation fails in that sense ast . Some remarks to a similar nonlinear problem are made.     

Measurement of thermal diffusivity of diamond by the light-induced thermal lattice method

The thermal diffusivity coefficient of natural diamonds is measured by optical induction of thermal diffraction lattices.Notation gc thermal conductivity coefficient - thermal diffusivity coefficient - diffraction efficiency - I_dif diffracted radiation intensity - I_Probe probe radiation intensity - probe radiation wavelength - c specific heat - Q surface energy density - thermal lattice relaxation time - thermal lattice period Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 56, No. 5, pp. 745–748, May, 1989.     

Effect of nonlinearity of gas mixtures on the process of their partition by thermal diffusion

Under consideration is the effect of nonideality of the components in a gas mixture on the process of their separation by thermal diffusion. It is demonstrated that in the expressions for the heat flux and the mass flux, the thermodiffusion ratio and the characteristic of diffusional thermal conductivity the effect of nonideality appears in the heat of mixing.Notation p pressure - density - length of the mean free path for molecules during transport of particles - length of the mean free path for particles during a transfer of the mean velocity - n molecule concentration - M molecular weight - I particle flux - J mass flux - m mass of a molecule - t time - D_ij coefficient of interdiffusion for a binary mixture - D _i ^T coefficient of thermal diffusion - K_T thermodiffusion ratio - _T thermodiffusion constant - x_i molar fraction of the i-th component in the mixture (r), intermolecular interaction potential - r intermolecular distance - collision integrals - T temperature - T^* referred temperature - R universal gas constant - k Boltzmann constant - Ñ Avogadro's number - v mean velocity of molecules - ¯V diffusion rate - _{i, trans} thermal conductivity associated with translatory degrees of freedom - f_i(r, v, t) velocity distribution function of molecules - viscosity - _i chemical potential of the i-th component - c_i mass fraction - _o thermal conductivity at the initial instant of time - thermal conductivity in the steady state - _DT diffusional component of thermal conductivity - g and h molar thermodynamic functions - ¯g and ¯h specific thermodynamic functions - c_p specific heat - J_q heat flux - J_q reduced heat flux - B second virial coefficient - U^* transport energy - coefficient of thermal expansion - coefficient of isothermal compression - f_i activity coefficient for the i-th mixture component Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 40, No. 5, pp. 829–839, May, 1981.     

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     

Retention of a metastable bcc ε-Pu solid solution: ε-Pu(Ti)

A metastable -Pu (bcc) solid solution has been retained to room temperature by rapid quenching of Pu-rich Pu-Ti alloys from the liquid state. Until now, -Pu solid solutions were limited to high temperatures and had not been successfully quenched to room temperature without transformation. The apparatus used to quench the specimens was a modified gun-type splat-cooling unit, capable of producing extremely high cooling rates of from 10⁶ to 10⁸ ° C sec^–1. -Pu(Ti) was retained in the composition region from 20 to more than 45 at.% Ti, and extrapolation of the lattice parameter/composition curve yielded a value of a ₀ = 3.53₀ Å for -Pu at 20° C. This modification differs from the -Pu modification derived by extrapolating from high temperature to 20° C by a small valence increase of 0.1. Metastable -Pu (Ti) (fcc) solid solutions were also quenched-in with alloys containing lesser amounts of Ti, and evidence was found to indicate that the was probably a product of -Pu(Ti) solid-state decomposition.     

Kinetic interpretation ofα- andβ-Si3N4 formation from oxide-free high-purity silicon powder

A kinetic analysis of the isothermal nitridation of high-purity oxide-free silicon powder is described. The kinetic analysis suggests that the and polymorphs of Si₃N₄ are formed by separate and parallel reaction paths. This analysis provides for the decoupling and quantitative kinetic interpretation of- and-Si₃N₄ formation reactions. Consistent with existing microstructural and thermodynamic evidence, the-forming reaction is shown to obey a first-order rate law, whereas a phase-boundary controlled rate law describes the-forming reaction. A kinetic model employing these rate laws is developed and is used to predict the/ phase ratio as a function of isothermal reaction temperature and extent of reaction. The/ phase ratios so obtained are shown to be in good agreement with experimental observations made under a variety of reaction conditions.