NONUNIQUENESS AND STABILITY FOR HEAT CONDUCTION THROUGH A DUPLEX HEAT EXCHANGER TUBE
|
|
- Felix Edwards
- 5 years ago
- Views:
Transcription
1 Journal of Thermal Stresses ISSN: (Print) X (Online) Journal homepage: NONUNIQUENESS AND STABILITY FOR HEAT CONDUCTION THROUGH A DUPLEX HEAT EXCHANGER TUBE J. R. Barber To cite this article: J. R. Barber (1986) NONUNIQUENESS AND STABILITY FOR HEAT CONDUCTION THROUGH A DUPLEX HEAT EXCHANGER TUBE, Journal of Thermal Stresses, 9:1, 69-78, DOI: / To link to this article: Published online: 21 May Submit your article to this journal Article views: 14 View related articles Citing articles: 15 View citing articles Full Terms & Conditions of access and use can be found at Download by: [University of Michigan] Date: 03 January 2018, At: 08:22
2 NONUNIQUENESS AND STABILITY FOR HEAT CONDUCTION THROUGH A DUPLEX HEAT EXCHANGER TUBE J. R. Barber Department ofmechanical Engineering and Applied Mechanics University of Michigan Ann Arbor, Michigan An analysis is given for radial heat flow throu8,h the wall ofa duplex heat exchanger tube with a pressure or gap-dependent contact resistance at the interface. Thermal expansion ofthe tube changes the value ofthis resistance, and it is shown that for certain thermal conditions there will be more than one steady-state solution. The stability ofthese multiple solutions is then investigated using a perturbation method, and it is found that there is always an odd number of solutions, alternately stable and unstable. This kind of behavior has been observed in other thermoelastic contact problems, but this is the first example in which it has been exhibited in a case ofcontact between similar materials. INTRODUCTION When heat is conducted across an interface between two solids, the thermal distortion of the solids will generally influence the contact conditions and hence the heat conduction process. If the materials of the two solids are different, thermal rectification occurs-i.e., the thermal resistance of the interface depends upon the direction of heat flow [1,2]. For a contact between two half-spaces, interface resistance increases with heat flow into the material with the higher distortivity 8, defined by a(1 + v) &=-- k (I) where a, v, and k: are, respectively, coefficient of thermal expansion, Poisson's ratio, and thermal conductivity. Mathematical difficulties were encountered in early attempts to treat problems of this kind. In particular, the conventional boundary conditions of perfect thermal The author gratefully acknowledges support from the U.S. Army under contract DAAG k 001. Journal of Thermal Stresses, 9:69-78, 1986 Copyright 1986 by Hemisphere Publishing Corporation 69
3 70 1. R. BARBER R contact seporotjon p- 0 _g Fig. 1 Typical variation of interface resistance, R with contact pressure p and gap g. contact (no resistance to heat flow, leading to continuity of temperature) in regions of mechanical contact and complete insulation (no heat flux) in regions of separation led to ill-posed boundary-value problems wherever the hotter solid had the higher distortivity [3-5]. An asymptotic analysis of the transitions between the various contact regions [5] suggests that the conventional boundary conditions can be safely applied when the heat flows in the opposite direction, but there is evidence that, in this case, the solution obtained is not necessarily unique [6]. In more recent papers [7,8], the author has probed this question by using more realistic boundary conditions that include a thermal contact resistance at the interface. It is anticipated that this resistance will increase monotonically with decreasing pressure in regions of contact [9] and be approximately linear with gap thickness in regions of separation (see Fig. 1). An idealization of this idea is the condition of "imperfect contact" introduced in [4] and subsequently used to solve a variety of thermoelastic contact problems [10-13]. Hitherto, it,has been assumed that these phenomena only arise in contact between two dissimilar materials, but the present paper will demonstrate their occurrence in an important heat transfer application involving similar materials. The author is grateful to Dr. M. G. Srinivasan for bringing this problem to his attention. STATEMENT OF THE PROBLEM We consider a duplex heat exchanger tube, as shown in Fig. 2, fabricated by shrinking one cylinder onto another, leaving the assembly in a state of prestress. Heat transfer occurs between the inner surface of the tube (radius r.) and a contained fluid Fig. 2 Geometry of the problem. The heat transfer coefficient is h, at radius r, and h, at radius r,. Both tubes are of the same material.
4 HEAT CONDUCTION THROUGH DUPLEX TUBE 71 at temperature T 1 through a coefficient hi and between the outer surface (radius '2) and surrounding fluid at temperature T 2 through a coefficient h 2 At the interface (radius '0) there is a contact resistance R that varies with pressure or gap in a monotonic but otherwise unprescribed way, as shown in Fig. l. It is assumed that T 1 and T 2 are axisymmetric and that the axial variation is sufficiently slow for the problem to be treated as one of two-dimensional plane strain. The initial state of the interface can be described by a radial interference do or an initial gap go. STEADY-STATE SOLUTION In the steady state, the temperature at radius, is where = T + ~ + qo'o log ~ '0 < r < '2 2 A 2h2 k r r, A. i = -, i = 1, 2 '0 and qo is the heat flux per unit area at the interface (radius '0), taken as positive when directed outwards. These equations imply the existence of a temperature jump.1 at the interface given by where (2) (3) (4) (5) (6) and (7) We can solve this equation for qo, obtaining
5 72 J. R. BARBER (8) The temperature field described in Eqs. (2) and (3) can also be used to calculate the free thermal expansion of the two cylinders, i.e., the radial displacement u, that would occur in the absence of interface or fluid pressure. Substituting Eqs. (2) and (3) into equation (c, 151) of Timoshenko and Goodier [14], we find (9) for the inner cylinder and.for the outer. If the steady-state solution involves separation at the interface, the gap must be where and (10) g = U,2 - UrI - do (11) Substituting from Eqs. (9) and (10), and eiiminating qo using Eq. (8), we find A~ A~ c = log A - log A 2 (A~ - 1) 2 (A~ - I) I (12) (13) (14) If Eq. (12) defines a negative value of g, interference is implied and we must superpose an elastic field with an interface pressure p chosen to reduce this interference to zero. Elementary calculations [14, 28] show that in plane strain an interface pressure p produces a differential radial displacement (15)
6 HEAT CONDUCTION TIlROUGH DUPLEX TUBE 73 and hence, from Eq. (12) and the condition of zero interference, p is given by Following [7,8], Eq. (16) can be subsumed under Eq. (12) by defining (-C3P) as a continuation of g into negative values. (16) GRAPHIC SOLUTION The geometry of the problem demands that r2 > ro > '1 and with this condition it can be shown that c\ > C2 > O. It follows thatf is a monotonically increasing function of R that tends to C2 - (>0) when R - 0 and to 1 when R C\ The function / is therefore monotonic with g as shown in Fig. 3, and we can obtain a graphic solution of Eq. (12) by writing it in the form The right-hand side of this equation can be plotted on Fig. 3 as a straight line passing through the point (-do, 0) with a slope of [0:(1 + v)ro(t 2 - TI))-I. Intersections of this line with with /(g) represent possible steady-state solutions of the problem. We can conclude from Fig. 3 that there will be one and only one such intersection if (T 2 - T I ) < 0 (i.e. if the straight line has negative slope) or if do < d.. where d, is the intercept obtained by drawing the tangent to the f(g) curve at the point of maximum slope and extending it to intersect the g axis. In the case of an initial gap go the hitter condition becomes -go < d.. U do > d, and T 2 > T.. there will be some finite range of (T 2 - TI) in which three solutions are obtained, but on either side of which there will be a unique solution. (17) 9+<1" ao+vjrj't,-t) fig) / / <1" 9<O,P>O; contact 0 0>0; separation Fig. 3 Graphical representation of Eq. (17). For the case illustrated there are threesolutions, A, B, and C, of which B is unstable and A and C are stable [condition given in Eq. (39)].
7 74 J. R. BARBER If there are several points of maximum slope giving intercepts d.. d 2, etc., a more complex behavior is obtained, with a possibility of more then three solutions in some ranges of (T 2 - T.). However, for given conditions, the number of possible solutions is always odd. We note that the conventional boundary conditions of perfect thermal contact (R = 0, p > 0) and perfect separation (R - 00, g > 0) correspond to the two line segments f = cdc.. p > 0; and f = 1, g > O. If there is a positive gap go under isothermal conditions, there will always be a range of (negative) values of (T 2 - for which the straight line in Fig. 3 does not intersect either of these line segments. As in previous problems, a solution can only be obtained in these cases if the idealized f curve is completed by the inclusion of the line segment g = 0, C2/C. < f < I, which corresponds to imperfect contact as defined in [4, 5, 7]. STABILITY CONSIDERATIONS Following [7,8], we investigate the stability of the preceding steady-state solutions by determining the conditions under which a small perturbation in the temperature field can grow exponentially with time t. The perturbation must satisfy the transient heat conduction equation where K is the thermal diffusivity. The exponentially growing solution of this equation is T = e D![A 1 / o(l3r) + B,Ko(l3r)], T = e D1 [A 2 / o(13r) + B 2Ko( 13r)], r. < r < ro ro < r < r2 where 13 2 = al«and la, Ko are modified Bessel functions of first and second kind, respectively. The constants A.. A 2, B.. and 8 2 are to be determined from the boundary conditions at r = r.. rs, and roo At r = r.. r2 we have radiation into a medium at zero temperature (since the fluid temperature is not perturbed) and hence T.) (18) (19) (20) (21) (22) giving the two equations (23)
8 HEAT CONDUCTION THROUGH DUPLEX TUBE 75 where At r = ro. we require continuity in heat flux and hence in atfar between the cylinders. giving (25) (26) LINEAR PERTURBATION OF THE INTERFACE RESISTANCE EQUATION The final boundary condition is obtained by linearizing the relation in Eq. (5) between 1:1, R. and- qo for small perturbations about the steady state. Differentiating, we obtain The perturbation in heat flux, &qo. is given by and the perturbation in 1:1 is at &qo = - k - (r == ro) = - - zae Q1 ar r«81:1 == Ttr = r(;) - T(r == rt) from Eqs. (19) and (20). 'Solving Eq. (13) for R. we find k (27) (28) (29) (30) ro [C1- C2 ] R==- ---c k 1 - f I (31) and differentiating (32) We can find the change in gap 8g by substituting Eqs. (19 and 20) into (c 151) of [14]. using /7 of [15] obtaining
9 76 J. R. BARBER _ [A)hlll(hIZ) - B)hIK1(h1Z) - A/t(z) + BtK1(Z)]} hi - 1 (33) CHARACTERISTIC EQUATION Finally, we substitute for R, ar, qo, 8qo, and lid from Eqs. (31), (32), (8), (28), and (30), respectively, into Eq. (27) and eliminate AI> BI> A 2, and B 2 using Eqs. (23), (24), and (26), to obtain the characteristic equation where A given steady-state solution of Eq. (12) will be unstable if arsmall perturbation grows exponentially, and this will be the case if Eq, (34) has one or more real roots other than zero. It can be shown that G(z) is an even function that varies monoton-. ically with z in the range z > 0 and that G(z) < 0 if z is sufficiently large, for all physically meaningful values of the parameters. It follows that a real positive root will occur if (34) Lt G(z) > Expanding the Bessel functions in Eq. (34) at small values of the argument, we find that the terms containing the factor (1 - f) cancel and (37) Hence, the solution will be unstable if
10 HEAT CONDUcnON THROUGH DUPLEX TUBE f' > a(i + v)(t 2 - T,)ro I 77 (39) since (CI - C2) > 0, (see above). Referring to Fig. 3, it is clear that this criterion describes those intersections at which the f(g) curve crosses (g + do)/a(l + v)ro(t 2 - T 1 ) from below with increasing g. Thus, solution B is unstable, while A and C are stable. More generally, there will always be an odd number of steady-state solutions that are alternately stable and unstable. If there is only one solution, it is always stable. We also note that as in [7,8] it is possible to define an "energy" function u(i + v)ro(t 2 - T 1 ) Jfdg - ~ (g + t4i V(g) = such that the condition for a steady-state solution is and the condition for stability is av -=0 ag Thus, in all respects, the problem treated in this paper exhibits the same behavior as the examples combined in [7,8], despite the fact that there is no difference in distortivity between the two contacting bodies. It is interesting to speculate whether the bodies being multiply-connected is an essential factor in causing this behavior with similar materials. REFERENCES (40) (41) (42) I. A. M. Clausing, Heat transfer at the interface of dissimilar metals-the influence of thermal strain, Int. J. Heat Mass Transfer vol. 9, pp , J. R. Barber, The effect of thermal distortion on constriction resistance, Int. J. Heat Mass Transfer. vol. 14, pp , J. R. Barber, Indentation of the semi-infinite elastic solid by a hot sphere, Int. J. Mech. Sci., vol. 15, pp , J. R. Barber, Contact problems involving a cooled punch, J. Elasticity. vol. 8, pp , M. Comninou and J. Dundurs, On the Barber boundary conditions for thermoelastic contact, J. App/. Mech., vol. 46, pp , M. Comninou and J. Dundurs, On lack of uniqueness in heat conduction through a solid to solid interface, ASME J. Heal Transfer. vol. 102, pp , 1980.
11 78 J. R. BARBER 7. J. R. Barber, J. Dundurs, and M. Comninou, Stability considerations in thermoelastic contact, J. AppJ. Mech., vol. 47, pp , J. R. Barber, Stability of thermoelastic contact for the Aldo model, J. Appl. Mech., vol. 48, pp , M. G. Cooper, B. B. Mikic, and M. M. Yovanovich, Thermal contact conductance, Int. J. Heal Mass Transfer, vol. 12, pp , J. R. Barber, Indentation of an elastic half-space by a cooled flat punch, Q. J. Mech. AppJ. Math., vol. 35, pp , M. Comninou, J. Dundurs, and J. R. Barber, Planar Hertz contact with heat conduction, J. Appl, Mech.. vol. 48, pp , M. Comninou, J. R. Barber, and J. Dundurs, Heat conduction through a flat punch, J. Appl. Mech.. vol. 48, pp , J. R. Barber and M. Comninou, The penny-shaped interface crack with heat flow: II Imperfect contact, J. Appl, Mech., vol. SO, pp , S. P. Timoshenko and J. N. Goodier, Theory of Elasticity, 3d ed., McGraw-Hill, New York, IS. M. Abramowitz and I. A. Stegun, eds., Handbook of mathematical functions, National Bureau of Standards, Washington, D.C Request reprints from J. R. Barber. Received February 20, 1985 Communicated by Donald E. Carlson
The distortion of a cylinder with non-uniform axial heat conduction
Journal of Elasticity 13 (1983) 441-447 1983 Martinus Nijhoff Publishers, The Hague. Printed in the Netherlands The distortion of a cylinder with non-uniform axial heat conduction J.R. BARBER Department
More informationEffect of cracks on the thermal resistance of aligned fiber composites
JOURNA OF APPIED PHYSICS VOUME 92, NUMBER 2 5 JUY 22 Effect of cracks on the thermal resistance of aligned fiber composites J. Dryden Department of Mechanical and Materials Engineering, University of Western
More informationThermoelastic Stresses in a Rod Subjected to Periodic Boundary Condition: An Analytical Treatment
Thermoelastic Stresses in a Rod Subjected to Periodic Boundary Condition: An Analytical Treatment Ahmet N. Eraslan Department of Engineering Sciences Middle East Technical University Ankara 06531, Turkey
More informationThe Subsurface Crack Under Conditions of Slip and Stick Caused by a Surface Normal Force
F.-K.Chang Maria Comninou Mem. ASME Sheri Sheppard Student Mem. ASME J. R. Barber Department of Mechanical Engineering and Applied Mechanics, University of Michigan, Ann Arbor, Mich. 48109 The Subsurface
More informationQuasi Static Thermal Stresses in A Limiting Thick Circular Plate with Internal Heat Generation Due To Axisymmetric Heat Supply
International Journal of Mathematics and Statistics Invention (IJMSI) E-ISSN: 2321 4767 P-ISSN: 2321-4759 Volume 1 Issue 2 ǁ December. 2013ǁ PP-56-63 Quasi Static Thermal Stresses in A Limiting Thick Circular
More informationExact Solution of an MHD Natural Convection Flow in Vertical Concentric Annulus with Heat Absorption
International Journal of Fluid Mechanics & Thermal Sciences 217; 3(5): 52-61 http://www.sciencepublishinggroup.com/j/ijfmts doi: 1.11648/j.ijfmts.21735.12 ISSN: 2469-815 (Print); ISSN: 2469-8113 (Online)
More informationChapter 4 TRANSIENT HEAT CONDUCTION
Heat and Mass Transfer: Fundamentals & Applications Fourth Edition Yunus A. Cengel, Afshin J. Ghajar McGraw-Hill, 2011 Chapter 4 TRANSIENT HEAT CONDUCTION LUMPED SYSTEM ANALYSIS Interior temperature of
More informationReview: Conduction. Breaking News
CH EN 3453 Heat Transfer Review: Conduction Breaking News No more homework (yay!) Final project reports due today by 8:00 PM Email PDF version to report@chen3453.com Review grading rubric on Project page
More informationSteam Generator Tubing Inspection
Steam Generator Tubing Inspection Analytical Determination of Critical Flaw Dimensions in Steam Generator Tubing I. Kadenko, N. Sakhno, R. Yermolenko, Nondestructive Examination Training and Certification
More informationModelling Built-up Steam Turbine Rotor Using Finite Difference Method
Modelling Built-up Steam Turbine Rotor Using Finite Difference Method Roger Li Dr. Adam Wittek School of Mechanical Engineering, University of Western Australia Mark Smart Verve Energy Abstract The low-pressure
More informationRing-shaped crack propagation in a cylinder under nonsteady cooling
High Performance Structures and Materials III 5 Ring-shaped crack propagation in a cylinder under nonsteady cooling V. A. Zhornik, Yu. A. Prokopenko, A. A. Rybinskaya & P. A. Savochka Department of Theoretical
More informationLaminated Beams of Isotropic or Orthotropic Materials Subjected to Temperature Change
United States Department of Agriculture Forest Service Forest Products Laboratory Research Paper FPL 375 June 1980 Laminated Beams of Isotropic or Orthotropic Materials Subjected to Temperature Change
More informationFINITE ELEMENT ANALYSIS OF FRICTIONALLY EXCITED THERMOELASTIC INSTABILITY
Journal of Thermal Stresses ISSN: 0149-5739 (Print) 1521-074X (Online) Journal homepage: http://www.tandfonline.com/loi/uths20 FINITE ELEMENT ANALYSIS OF FRICTIONALLY EXCITED THERMOELASTIC INSTABILITY
More informationContact Modeling of Rough Surfaces. Robert L. Jackson Mechanical Engineering Department Auburn University
Contact Modeling of Rough Surfaces Robert L. Jackson Mechanical Engineering Department Auburn University Background The modeling of surface asperities on the micro-scale is of great interest to those interested
More informationThe Penny-Shaped Interface Crack With Heat Flow
C.J. Martin-Moran Department of Mechanical Engineering, University of Newcastle-Upon-Tyne, NE1 7RU, England J. R. Barber Department of Mechanical Engineering and Applied Mechanics. M. Comninou Department
More informationREE Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics
REE 307 - Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics 1. Is the following flows physically possible, that is, satisfy the continuity equation? Substitute the expressions for
More informationMechanics of Materials and Structures
Journal of Mechanics of Materials and Structures ON TORSIONAL VIBRATIONS OF INFINITE HOLLOW POROELASTIC CYLINDERS M. Tajuddin and S. Ahmed Shah Volume 2, Nº 1 January 27 mathematical sciences publishers
More informationArbor, Michigan, USA Published online: 27 Apr 2007.
This article was downloaded by: [the Bodleian Libraries of the University of Oxford] On: 01 October 2013, At: 11:31 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number:
More informationQUESTION ANSWER. . e. Fourier number:
QUESTION 1. (0 pts) The Lumped Capacitance Method (a) List and describe the implications of the two major assumptions of the lumped capacitance method. (6 pts) (b) Define the Biot number by equations and
More informationModeling Contact between Rigid Sphere and Elastic Layer Bonded to Rigid Substrate
IEEE TRANSACTIONS ON COMPONENTS AND PACKAGING TECHNOLOGIES, VOL. 24, NO. 2, JUNE 2001 207 Modeling Contact between Rigid Sphere and Elastic Layer Bonded to Rigid Substrate Mirko Stevanović, M. Michael
More informationA FINITE ELEMENT STUDY OF ELASTIC-PLASTIC HEMISPHERICAL CONTACT BEHAVIOR AGAINST A RIGID FLAT UNDER VARYING MODULUS OF ELASTICITY AND SPHERE RADIUS
Proceedings of the International Conference on Mechanical Engineering 2009 (ICME2009) 26-28 December 2009, Dhaka, Bangladesh ICME09- A FINITE ELEMENT STUDY OF ELASTIC-PLASTIC HEMISPHERICAL CONTACT BEHAVIOR
More informationMicroelectronics Heat Transfer Laboratory
Microelectronics Heat Transfer Laboratory Department of Mechanical Engineering University of Waterloo Waterloo, Ontario, Canada http://www.mhtl.uwaterloo.ca Outline Personnel Capabilities Facilities Research
More informationTHERMOELASTIC ANALYSIS OF THICK-WALLED FINITE-LENGTH CYLINDERS OF FUNCTIONALLY GRADED MATERIALS
Journal of Thermal Stresses, 28: 391 408, 2005 Copyright # Taylor & Francis Inc. ISSN: 0149-5739 print/1521-074x online DOI: 10.1080/01495730590916623 THERMOELASTIC ANALYSIS OF THICK-WALLED FINITE-LENGTH
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationNonlinear Stability, Thermoelastic Contact, and the Barber Condition
J. A. Pelesko School of Mathematics, Georgia Institute of Technology, Atlanta, GA 3332-6 Nonlinear Stability, Thermoelastic Contact, and the Barber Condition The behavior of a one-dimensional thermoelastic
More informationROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring
ROTATIONAL STRESSES INTRODUCTION High centrifugal forces are developed in machine components rotating at a high angular speed of the order of 100 to 500 revolutions per second (rps). High centrifugal force
More informationJoule Heating Effect on the Coupling of Conduction with Magnetohydrodynamic Free Convection Flow from a Vertical Flat Plate
Nonlinear Analysis: Modelling and Control, 27, Vol. 12, No. 3, 37 316 Joule Heating Effect on the Coupling of Conduction with Magnetohydrodynamic Free Convection Flow from a Vertical Flat Plate M. A. Alim
More informationIf there is convective heat transfer from outer surface to fluid maintained at T W.
Heat Transfer 1. What are the different modes of heat transfer? Explain with examples. 2. State Fourier s Law of heat conduction? Write some of their applications. 3. State the effect of variation of temperature
More informationChapter 2 HEAT CONDUCTION EQUATION
Heat and Mass Transfer: Fundamentals & Applications 5th Edition in SI Units Yunus A. Çengel, Afshin J. Ghajar McGraw-Hill, 2015 Chapter 2 HEAT CONDUCTION EQUATION Mehmet Kanoglu University of Gaziantep
More informationAn Experimental Study of Stress Singularities at a Sharp Corner in a Contact Problem
An Experimental Study of Stress Singularities at a Sharp Corner in a Contact Problem by E.C. Miniatt, A.M. Waas and W.J. Anderson ABSTRACT--The photoelastic method was used to investigate the nature of
More informationChapter 2 HEAT CONDUCTION EQUATION
Heat and Mass Transfer: Fundamentals & Applications Fourth Edition Yunus A. Cengel, Afshin J. Ghajar McGraw-Hill, 2011 Chapter 2 HEAT CONDUCTION EQUATION Mehmet Kanoglu University of Gaziantep Copyright
More informationChapter 4: Transient Heat Conduction. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University
Chapter 4: Transient Heat Conduction Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University Objectives When you finish studying this chapter, you should be able to: Assess when the spatial
More informationSome anomalies of the curved bar problem in classical elasticity
Journal of Elasticity 26: 75-86, 1991. 1991 Kluwer Academic Publishers. Printed in the Netherlands. Some anomalies of the curved bar problem in classical elasticity J.R. BARBER Department of Mechanical
More informationUniversity of Rome Tor Vergata
University of Rome Tor Vergata Faculty of Engineering Department of Industrial Engineering THERMODYNAMIC AND HEAT TRANSFER HEAT TRANSFER dr. G. Bovesecchi gianluigi.bovesecchi@gmail.com 06-7259-727 (7249)
More informationOn the Barber Boundary Conditions for Thermoelastic Contact
Maria Comninou Associate Professor, Department of Applied Mechanics and Engineering Science, University of Michigan, Ann Arbor, Mich. 48109 Mem. ASME J. Dundurs Professor, Department of Civil Engineering)
More informationE. not enough information given to decide
Q22.1 A spherical Gaussian surface (#1) encloses and is centered on a point charge +q. A second spherical Gaussian surface (#2) of the same size also encloses the charge but is not centered on it. Compared
More information1. Nusselt number and Biot number are computed in a similar manner (=hd/k). What are the differences between them? When and why are each of them used?
1. Nusselt number and Biot number are computed in a similar manner (=hd/k). What are the differences between them? When and why are each of them used?. During unsteady state heat transfer, can the temperature
More informationThe effect of rigidity on torsional vibrations in a two layered poroelastic cylinder
Int. J. Adv. Appl. Math. and Mech. 3(1) (2015) 116 121 (ISSN: 2347-2529) Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics The effect of rigidity on
More informationA Finite Element Study of Elastic-Plastic Hemispherical Contact Behavior against a Rigid Flat under Varying Modulus of Elasticity and Sphere Radius
Engineering, 2010, 2, 205-211 doi:10.4236/eng.2010.24030 Published Online April 2010 (http://www. SciRP.org/journal/eng) 205 A Finite Element Study of Elastic-Plastic Hemispherical Contact Behavior against
More informationTransient Analysis of Disk Brake By using Ansys Software
Transient Analysis of Disk Brake By using Ansys Software G. Babukanth & M. Vimal Teja Department of Mechanical Engineering, Nimra College of Engineering & Technology, Ibrahimpatnam, Vijayawada E-mail :
More informationDocumentation of the Solutions to the SFPE Heat Transfer Verification Cases
Documentation of the Solutions to the SFPE Heat Transfer Verification Cases Prepared by a Task Group of the SFPE Standards Making Committee on Predicting the Thermal Performance of Fire Resistive Assemblies
More informationReview of Thermal Joint Resistance Models for Non-Conforming Rough Surfaces in a Vacuum
Review of Thermal Joint Resistance Models for Non-Conforming Rough Surfaces in a Vacuum M. Bahrami J. R. Culham M. M. Yovanovich G. E. Schneider Department of Mechanical Engineering Microelectronics Heat
More informationChapter 23. Gauss Law. Copyright 2014 John Wiley & Sons, Inc. All rights reserved.
Chapter 23 Gauss Law Copyright 23-1 Electric Flux Electric field vectors and field lines pierce an imaginary, spherical Gaussian surface that encloses a particle with charge +Q. Now the enclosed particle
More informationLightning Phenomenology Notes Note 23 8 Jan Lightning Responses on a Finite Cylindrical Enclosure
Lightning Phenomenology Notes Note 23 8 Jan 2014 Lightning Responses on a Finite Cylindrical Enclosure Kenneth C. Chen and Larry K. Warne Sandia National Laboratories, P. O. Box 5800, Albuquerque, NM 87185,
More informationApplication of Maximum Principal Strain Theory for Study of Coal Particle Disintegration when Subjected to Detonation Wave
Application of Maximum Principal Strain Theory for Study of Coal Particle Disintegration when Subjected to Detonation Wave Patadiya D. M. 1,1, Sheshadri T. S. 1,, Jaishankar S 1,3. 1,1 Department of Aerospace
More informationConvection When the radial flux of energy is carried by radiation, we derived an expression for the temperature gradient: dt dr = - 3
Convection When the radial flux of energy is carried by radiation, we derived an expression for the temperature gradient: dt dr = - 3 4ac kr L T 3 4pr 2 Large luminosity and / or a large opacity k implies
More informationOscillatory flow of a jeffrey fluid in an elastic tube of variable cross-section
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research 2012 3 (2):671-677 ISSN: 0976-8610 CODEN (USA): AASRFC Oscillatory flow of a jeffrey fluid in an elastic tube of
More informationAdvanced Structural Analysis EGF Cylinders Under Pressure
Advanced Structural Analysis EGF316 4. Cylinders Under Pressure 4.1 Introduction When a cylinder is subjected to pressure, three mutually perpendicular principal stresses will be set up within the walls
More informationThermal Analysis Contents - 1
Thermal Analysis Contents - 1 TABLE OF CONTENTS 1 THERMAL ANALYSIS 1.1 Introduction... 1-1 1.2 Mathematical Model Description... 1-3 1.2.1 Conventions and Definitions... 1-3 1.2.2 Conduction... 1-4 1.2.2.1
More informationExperimental Studies of Active Temperature Control in Solid Breeder Blankets
Experimental Studies of Active Temperature Control in Solid Breeder Blankets M. S. Tillack, A. R. Raffray, A. Y. Ying, M. A. Abdou, and P. Huemer Mechanical, Aerospace and Nuclear Engineering Department
More informationAnalytical Strip Method for Thin Isotropic Cylindrical Shells
IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 14, Issue 4 Ver. III (Jul. Aug. 2017), PP 24-38 www.iosrjournals.org Analytical Strip Method for
More informationPrediction of geometric dimensions for cold forgings using the finite element method
Journal of Materials Processing Technology 189 (2007) 459 465 Prediction of geometric dimensions for cold forgings using the finite element method B.Y. Jun a, S.M. Kang b, M.C. Lee c, R.H. Park b, M.S.
More informationHEAT CONDUCTION USING GREEN S FUNCTIONS
HEAT CONDUCTION USING GREEN S FUNCTIONS Preface to the first edition Preface to the second edition Author Biographies Nomenclature TABLE OF CONTENTS FOR SECOND EDITION December 2009 Page viii x xii xiii
More informationIntroduction to Heat and Mass Transfer. Week 5
Introduction to Heat and Mass Transfer Week 5 Critical Resistance Thermal resistances due to conduction and convection in radial systems behave differently Depending on application, we want to either maximize
More informationDAMPING OF GENERALIZED THERMO ELASTIC WAVES IN A HOMOGENEOUS ISOTROPIC PLATE
Materials Physics and Mechanics 4 () 64-73 Received: April 9 DAMPING OF GENERALIZED THERMO ELASTIC WAVES IN A HOMOGENEOUS ISOTROPIC PLATE R. Selvamani * P. Ponnusamy Department of Mathematics Karunya University
More informationSolutions for the Practice Final - Math 23B, 2016
olutions for the Practice Final - Math B, 6 a. True. The area of a surface is given by the expression d, and since we have a parametrization φ x, y x, y, f x, y with φ, this expands as d T x T y da xy
More informationFinite element solution to transient asymmetric heat conduction in multilayer annulus
Int. J. dv. ppl. Math. and Mech. 2(3) (2015) 119-125 (ISSN: 2347-2529) Journal homepage: www.ijaamm.com International Journal of dvances in pplied Mathematics and Mechanics Finite element solution to transient
More informationTransient Heat Transfer Experiment. ME 331 Introduction to Heat Transfer. June 1 st, 2017
Transient Heat Transfer Experiment ME 331 Introduction to Heat Transfer June 1 st, 2017 Abstract The lumped capacitance assumption for transient conduction was tested for three heated spheres; a gold plated
More informationHYDRODYNAMIC ANALOGY STRESS DTSTRIBUTION IN
RFJ RJNT 2 SFRIES J HYDRODYNAMIC ANALOGY STRESS DTSTRIBUTION IN BY M. A. BIOT Reprinted from JOURNAL OF APPLIED PHYSICS, Vol. 9, No. 1, January, 1938 Copyright 1937 by the American Institute of Physics
More informationA HYBRID SEMI-ANALYTICAL AND NUMERICAL METHOD FOR MODELING WELLBORE HEAT TRANSMISSION
PROCEEDINGS, Thirtieth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 31-February 2, 5 SGP-TR-176 A HYBRID SEMI-ANALYTICAL AND NUMERICAL METHOD FOR MODELING
More informationModule 8: BoiIing Lecture 34: Analysis of Bubble Formation on a Heated Wall. The Lecture Contains: Bubble Growth with Heat and Mass Transfer
The Lecture Contains: Bubble Growth with Heat and Mass Transfer Simple Analysis of Bubble Growth in a Uniform Temperature Field file:///d /Web%20Course%20(Ganesh%20Rana)/Dr.%20gautam%20biswas/Final/convective_heat_and_mass_transfer/lecture34/34_1.html[12/24/2014
More information1.1 GRAPHS AND LINEAR FUNCTIONS
MATHEMATICS EXTENSION 4 UNIT MATHEMATICS TOPIC 1: GRAPHS 1.1 GRAPHS AND LINEAR FUNCTIONS FUNCTIONS The concept of a function is already familiar to you. Since this concept is fundamental to mathematics,
More informationS.E. (Chemical) (Second Semester) EXAMINATION, 2012 HEAT TRANSFER (2008 PATTERN) Time : Three Hours Maximum Marks : 100
Total No. of Questions 12] [Total No. of Printed Pages 7 Seat No. [4162]-187 S.E. (Chemical) (Second Semester) EXAMINATION, 2012 HEAT TRANSFER (2008 PATTERN) Time : Three Hours Maximum Marks : 100 N.B.
More informationBESSEL FUNCTIONS APPENDIX D
APPENDIX D BESSEL FUNCTIONS D.1 INTRODUCTION Bessel functions are not classified as one of the elementary functions in mathematics; however, Bessel functions appear in the solution of many physical problems
More informationStress analysis of functionally graded discs under mechanical and thermal loads
Indian Journal of Engineering & Materials Sciences Vol. 8, April 0, pp. -8 Stress analysis of functionally graded discs under mechanical and thermal loads Hasan Çallioğlu, Metin Sayer* & Ersin Demir Department
More informationSOLVING THE HEAT EQUATION IN THE EDPS METHOD
SOLVING THE HEAT EQUATION IN THE EDPS METHOD Peter Krupa 1, Svetozár Malinarič 1, Peter Dieška 2 1 Department of Physics, Faculty of Natural Sciences, Constantine the Philosopher University in Nitra, Tr.
More informationMassachusetts Institute of Technology 22.68J/2.64J Superconducting Magnets. February 27, Lecture #4 Magnetic Forces and Stresses
Massachusetts Institute of Technology.68J/.64J Superconducting Magnets February 7, 003 Lecture #4 Magnetic Forces and Stresses 1 Forces For a solenoid, energy stored in the magnetic field acts equivalent
More informationIntroduction. Statement of Problem. The governing equations for porous materials with Darcy s law can be written in dimensionless form as:
Symbolic Calculation of Free Convection for Porous Material of Quadratic Heat Generation in a Circular Cavity Kamyar Mansour Amirkabir University of technology, Tehran, Iran, 15875-4413 mansour@aut.ac.ir
More informationStress Function of a Rotating Variable-Thickness Annular Disk Using Exact and Numerical Methods
ngineering 0 3 4-430 doi:0.436/eng.0.34048 Published Online April 0 (http://www.scip.org/ournal/eng) Stress Function of a otating Variable-Thicness Annular Dis Using xact and Numerical Methods Abstract
More informationPrinciples of Food and Bioprocess Engineering (FS 231) Exam 2 Part A -- Closed Book (50 points)
Principles of Food and Bioprocess Engineering (FS 231) Exam 2 Part A -- Closed Book (50 points) 1. Are the following statements true or false? (20 points) a. Thermal conductivity of a substance is a measure
More informationSpecific heat capacity. Convective heat transfer coefficient. Thermal diffusivity. Lc ft, m Characteristic length (r for cylinder or sphere; for slab)
Important Heat Transfer Parameters CBE 150A Midterm #3 Review Sheet General Parameters: q or or Heat transfer rate Heat flux (per unit area) Cp Specific heat capacity k Thermal conductivity h Convective
More informationChemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017
Chemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017 Objective: Text: To introduce the basic concepts of fluid mechanics and heat transfer necessary for solution of engineering
More informationThermo-Elastic Stress Analysis of the GHARR-1 Vessel during Reactor Operation Using ANSYS 13.0
Journal of Emerging Trends in Engineering and Applied Sciences (JETEAS) 6(5): 309-314 Scholarlink Research Institute Journals, 2015 (ISSN: 2141-7016) jeteas.scholarlinkresearch.com Journal of Emerging
More information( ) PROBLEM C 10 C 1 L m 1 50 C m K W. , the inner surface temperature is. 30 W m K
PROBLEM 3. KNOWN: Temperatures and convection coefficients associated with air at the inner and outer surfaces of a rear window. FIND: (a) Inner and outer window surface temperatures, T s,i and T s,o,
More informationChapter 5 Time-Dependent Conduction
Chapter 5 Time-Dependent Conduction 5.1 The Lumped Capacitance Method This method assumes spatially uniform solid temperature at any instant during the transient process. It is valid if the temperature
More information3.3 Unsteady State Heat Conduction
3.3 Unsteady State Heat Conduction For many applications, it is necessary to consider the variation of temperature with time. In this case, the energy equation for classical heat conduction, eq. (3.8),
More informationStresses Analysis of Petroleum Pipe Finite Element under Internal Pressure
ISSN : 48-96, Vol. 6, Issue 8, ( Part -4 August 06, pp.3-38 RESEARCH ARTICLE Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure Dr.Ragbe.M.Abdusslam Eng. Khaled.S.Bagar ABSTRACT
More informationEffective Scraping in a Scraped Surface Heat Exchanger: Some Fluid Flow Analysis
ICEF 9 003 Effective Scraping in a Scraped Surface Heat Exchanger: Some Fluid Flow Analysis D. L. Pyle (1), K.-H. Sun (1), M. E. M. Lee (), C. P. Please (), A. D. Fitt (), S. K. Wilson (3), B. R. Duffy
More informationCRACK-TIP DRIVING FORCE The model evaluates the eect of inhomogeneities by nding the dierence between the J-integral on two contours - one close to th
ICF 100244OR Inhomogeneity eects on crack growth N. K. Simha 1,F.D.Fischer 2 &O.Kolednik 3 1 Department ofmechanical Engineering, University of Miami, P.O. Box 248294, Coral Gables, FL 33124-0624, USA
More informationElectric Flux. To investigate this, we have to understand electric flux.
Problem 21.72 A charge q 1 = +5. nc is placed at the origin of an xy-coordinate system, and a charge q 2 = -2. nc is placed on the positive x-axis at x = 4. cm. (a) If a third charge q 3 = +6. nc is now
More informationFinite element Analysis of thermo mechanical stresses of two layered composite cylindrical pressure vessel
Finite element Analysis of thermo mechanical stresses of two layered composite cylindrical pressure vessel Arnab Choudhury *1, Samar Chandra Mondol #2, Susenjit Sarkar #3 *Faculty, Mechanical Engineering
More informationAvailable online at ScienceDirect. Procedia Engineering 144 (2016 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 144 (016 ) 170 177 1th International Conference on Vibration Problems, ICOVP 015 Propagation of Torsional Surface Wave in Anisotropic
More informationEffect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study
Effect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study Biplab Chatterjee, Prasanta Sahoo 1 Department of Mechanical Engineering, Jadavpur University
More informationBUCKLING COEFFICIENTS FOR SANDWICH CYLINDERS OF FINITE LENGTH UNDER UNIFORM EXTERNAL LATERAL PRESSURE
U.S. DEPARTMENT OF AGRICULTURE FOREST SERVICE FOREST PRODUCTS LABORATORY MADlSON, WIS. In Cooperation with the University of Wisconsin U. S. FOREST SERVICE RESEARCH NOTE FPL-0104 SEPTEMBER 1965 BUCKLING
More informationModeling Transient Conduction in Enclosed Regions Between Isothermal Boundaries of Arbitrary Shape
JOURNAL OF THERMOPHYSICS AND HEAT TRANSFER Vol. 19, No., July September 2005 Modeling Transient Conduction in Enclosed Regions Between Isothermal Boundaries of Arbitrary Shape Peter Teertstra, M. Michael
More informationCopyright 1988 by ASME
F. A. Sturla J. R. Barber Department of Mechanical Engineering and Applied Mechanics, The University of Michigan, Ann Arbor, Ml 48109 Thermal Stresses Due to a Plane Crack in General Anisotropic Material
More informationStress and Temperature Distribution Study in a Functionally Graded Brake Disk
Stress and Temperature Distribution Study in a Functionally Graded Brake Disk P. Hosseini Tehrani, M.Talebi2 Assistant professor, 2 MS c. Student, Center of excellence in railway transportation, School
More informationChapter 8. Conservation Laws. 8.3 Magnetic Forces Do No Work
Chapter 8. Conservation Laws 8.3 Magnetic Forces Do No Work 8.2 Momentum of EM fields 8.2.1 Newton's Third Law in Electrodynamics Consider two charges, q 1 and q 2, moving with speeds v 1 and v 2 magnetic
More informationChapter 3 NATURAL CONVECTION
Fundamentals of Thermal-Fluid Sciences, 3rd Edition Yunus A. Cengel, Robert H. Turner, John M. Cimbala McGraw-Hill, 2008 Chapter 3 NATURAL CONVECTION Mehmet Kanoglu Copyright The McGraw-Hill Companies,
More informationSimulation of a two-dimensional sheath over a flat wall with an insulatorõconductor interface exposed to a high density plasma
JOURNAL OF APPLIED PHYSICS VOLUME 94, NUMBER 5 1 SEPTEMBER 2003 Simulation of a two-dimensional sheath over a flat wall with an insulatorõconductor interface exposed to a high density plasma Doosik Kim
More informationCombined Effect of Magnetic field and Internal Heat Generation on the Onset of Marangoni Convection
International Journal of Fluid Mechanics & Thermal Sciences 17; 3(4): 41-45 http://www.sciencepublishinggroup.com/j/ijfmts doi: 1.11648/j.ijfmts.1734.1 ISSN: 469-815 (Print); ISSN: 469-8113 (Online) ombined
More informationFinite Element Modeling for Transient Thermal- Structural Coupled Field Analysis of a Pipe Joint
International Conference on Challenges and Opportunities in Mechanical Engineering, Industrial Engineering and Management Studies 88 Finite Element Modeling for Transient Thermal- Structural Coupled Field
More informationFRACTURE MECHANICS OF COMPOSITES WITH RESIDUAL STRESSES, TRACTION-LOADED CRACKS, AND IMPERFECT INTERFACES
Proc. 2 nd ESIS TC4 Conference on Polymers and Composites, in press, 1999 Author prepared reprint FRACTURE MECHANICS OF COMPOSITES WITH RESIDUAL STRESSES, TRACTION-LOADED CRACKS, AND IMPERFECT INTERFACES
More informationStudies on flow through and around a porous permeable sphere: II. Heat Transfer
Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India
More informationArbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact
Arbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact Philip P. Garland 1 and Robert J. Rogers 2 1 School of Biomedical Engineering, Dalhousie University, Canada 2 Department
More informationChapter 2: Heat Conduction. Dr Ali Jawarneh Department of Mechanical Engineering, Hashemite University
Chapter : Heat Conduction Equation Dr Ali Jawarneh Department of Mechanical Engineering, Hashemite University Objectives When you finish studying this chapter, you should be able to: Understand multidimensionality
More informationFinite Difference Solution of Unsteady Free Convection Heat and Mass Transfer Flow past a Vertical Plate
Daffodil International University Institutional Repository DIU Journal of Science and Technology Volume 1, Issue 1, January 17 17-1 Finite Difference Solution of Unsteady Free Convection Heat and Mass
More informationHeat processes. Heat exchange
Heat processes Heat exchange Heat energy transported across a surface from higher temperature side to lower temperature side; it is a macroscopic measure of transported energies of molecular motions Temperature
More informationQuestions Chapter 23 Gauss' Law
Questions Chapter 23 Gauss' Law 23-1 What is Physics? 23-2 Flux 23-3 Flux of an Electric Field 23-4 Gauss' Law 23-5 Gauss' Law and Coulomb's Law 23-6 A Charged Isolated Conductor 23-7 Applying Gauss' Law:
More informationp + µ 2 u =0= σ (1) and
Master SdI mention M2FA, Fluid Mechanics program Hydrodynamics. P.-Y. Lagrée and S. Zaleski. Test December 4, 2009 All documentation is authorized except for reference [1]. Internet access is not allowed.
More information