NUMElUCAL ANALYSIS OF FOI{CU> CONVECTION INFLI IJ<:N('E ON Tli E COUHSE OF SOLI I>IFICATION PROCESS
|
|
- Lawrence Richards
- 5 years ago
- Views:
Transcription
1 Solitliłkation of i\ldah a111l..\lln.d KrL e pn i ę cte: M e talt t S tor>ńw. IR 1'1. ISSN H.WX - 'J.\Xfl NUMElUCAL ANALYSIS OF FOI{CU> CONVECTION INFLI IJ<:N('E ON Tli E COUHSE OF SOLI I>IFICATION PROCESS jan Szajnar Fot/1/ilrv 111.\lifllft'. Silnilfil /'edtnic lll l i ni\'1'1'\ ll\', (.'!i ll'ic. l'olrllul ABSTRACT In this paper the mathcmatical model and numerical computations cor1eerning the pure metal (Al 99.7) solidification in rotational magnctic field are prcscnted. This ticld causes the convectional movemcnts of molten metal (forced convection) and in sub-arca corrcsponding to liquid stale t he heat transfer conditions change in essential way. This process has unquestionable influence on the course of solidification and temperature field in casting domain. The aim of researeh was to eonstruci the computer program simulating heat transfer in considered area at the same time the convectional movements of molten metal should be taken into account. The problem was solved on the basis of eontroi volume method (in enthalpy convention). In this paper the macroscepie approach is prcscnted it means trat boundary-initial problem described by Fourier's equations and adequate conditions i s analyzed. In the fina! part of the paper the results of computations are shown. The problems concerning the crystallization of metal (macroscopic-microscopic scale) will be a subject of future research works. ANALIZA NUMERYCZNA WPLYWU WYMUSZONEJ KONWEKCJI W CIEKłYM METALU NA PRZEBIEG PROCESU KRZF.PNIĘCIA STRESZCZENIE W pracy przedstawiono model matematyczny i obliczenia numeryczne dotycz4ce krzepnięcia czystego metalu (Al 99. 7) w obecności wirującego pola magnetycznego. Wirujące pole magnetyczne powoduje ruchy konwekcyjne fazy cieklej w krzepnącym odlewie (konwekcja wymuszona). Tak więc w podobszarze zajmowanym przez ciekły metal zmieniają si ~ istotnie warunki przepływu ciepła co ma niewątpliwy wpływ na kinetykę krzepnięcia i pole temperatury w krzepnącym odlewie. Celem pracy była konstrukcja programu komputerowego symulującego przepływ ciepła w odlewie przy uwzględnieniu wpływu przemieszczania się cieklego metalu w rozważanym obszarze. Zaclanie rozwiązano wykorzystując metodę bilansów elementarnych (control volume methocl) w ujęciu entalpowym. Wyniki oblicze11 zostały przedstawione w kork owej C7c~ci artykułu. W pracy prezentowane jest ujęcie makroskopowe bazujące na równaniach energii 117. upełnionych odpowiednimi warunkami brzegowo - początkowymi. Problemy dotyojcc krv,talizacji mctal11 w polu magnetycznym będą przedmiotem dalszych badań.
2 170 MATI-IEMATIC:\1. l >I.'>C'Rli"IION A non -stla!lly tempc r;~tu re lido in the domain of solidified and cooling metal is described by a system of tw11 d,fferential cquat1ons conceming a liquid and solid parts of casting, namely at".(x, t) ] c (T )p (T) +w gradt.(x,t) = [ "' "' "' '" at (l) = div[a",(t".)gradt.,(x, t) j, m=l, 2 where c,., p,., t..m are the thermophysical parameters of liquid and solid state, wis a velocity field in considered object, T, X, r are the temperature, spatial co-ordinates and time. The expression in b racket s on t he left s ide o f eqn (l) i s called a substantial derivative. The velocity field i s generated by action o f rotational magnetic field and w=o for XE D 2 (t), w;ćo for XED 1 (t). In this paper it is assumed that intensive stirring of molten metal caused by magnetic field action takes place in the all sub-area of liquid zone except of boundary layer close to interface r (comp. Fig. l). Because it is rather difficult to determine in convincing way the thickness of this layer therefore the some values of this quantity in numerical s ~ mulation has been taken into account. flow o f me l t D 0 (t) f* (t) l l l l / mol d Fig. l. The ronsidert"-"j suh-areas
3 l i l The adequate differential equations and bottndary - initial conditions have been transformed to so-called enthalpy convention. The physical enthalpy of metal related to <111 unit of volume is defined as T H (X. t) = f c ( 11) p (li) d li + 11 (T) L o (2) where T< T' ll(t) = l o' l ' T< T' T~ T' (3) L is a latent heat related to an unit of volume, T' - solidification point (Fig. 2). [J /m 3] L Fig. 2. Entalphy distribulion The energy equations expressed in cnthalpy convention are of t he form at the same time a,.=>...,! p.., c,.
4 172 Along the surface r (t) the Stefan boundary condition is accepted J - a 1 n gradh 1 (X, t) = - a 2 n gradh 2 (X, t)+ Lv" l Al = Az + L (5) The. sense of A 1 and A 2 symbols is shown in Fig. 2, v. is a solidification rate in normai to interface r direction. Along the contact surface between casting and mold the 1st kind condition resulting from Schwarz' solution is assumed (6) where He is the enthalpy corresponding to contact temperature found on the basis of Schwarz' solution [1]. In conventionally assumed limits of considered object the adiabatic conditions in the form XEr.: n gradh(x, t) = O (7) have been taken into account. Initial condition describing enthalpy field at the moment t=o is of the form t= O: H(X, O) = HP. (8) The physical experiments show that molten metal stirring caused by rotational magnetic field is very strong and the temperature of liquid core is practically homogenous. In numerical reali zation this effect can be achieved by the acceptation for this sub-area the thermophysical parameter called effective thermal conductivity [2]. This one is a multiple of real thermal conductivity A 1 and in practice the values of range ~=7>.. 1 are assumed. The introducing of ~ in mathematical description allows to neglect the substantial derivative in eqns (l) and (4) however it is necessary to consider three sub-domains, namely solid state D 2 (t), boundary layer D 1 (t) and shifting melt area D 0 (t). Finally one obtains the system of equations ah (X, t) "" = div[a grad H (X, t)j, m =0, l, 2 at "" "' where index O identilies the convectional movements sub-area. Along the surface between D 0 (t) and D 1 (t) the conditions of temperature and heat flux continuity should be assumed. CONTROI. VOLUME METHOD In this papcr the object shown in Fig. l has been considered. lt was assumed that solidificat.ion process procecds in normai to mold surface direction and the problem was treated as a one d i mensionaj - at the same time in adequate energy bajances the considered form o f liquid - solid in terface h as been taken in to account (moditication o f so-called shape functions in adequate difference equations). Control volume method depends on the (9)
5 11\ discrt'tintion of consrdern l an a 011, Jn ncnr ary, nltllrlcs <t ihi lll.tkillg up of encr)!y halancc.> for distinguishcd vollltnc ~. J'ltc SILpol.\1 rcsulrin g from l ill it' dr s,,, rization is 1 nnsidcted. For elementary volume.1 V,, one obtai11\ thc fpllov.ing c qu a t~o tl\ (lo) where T(ll) is an inver.\c fun ction to fun CIIOil // -= /1(7') (comp. Fig. 2), 1 1, t 1 " - the poinls being the heginning and end of time interval tj./,.,e" distinguished dircction of heat flow, R, -!herma! resistances in dircction' e. <t>, - sltape functions which are the ratio of a surface area!'j.s, Iimiting element tj. V,, in.. e" direction to volume tj. V., [3, 4], m - the number of considered directions. Above bajances written for all eontroi volumc\ eonstilule a lincar sysrem of equations from which on the basis of initial and next pscudo-initial condition (cnthalpy field at the moment t f) one can directly fin d t he cnthalpy field for lcvcl o f t i me. Presented algorithm belongs to the group of explicit schemcs it means that main matrix of the system has only one non zero diagonal This schcme must fulfili the stability condition [5] Iimiting acceptable interval of time!'j.t. As it was mcntioncd in chapter l the algorithm of heat transfer computations in considcrcd ohjcct has bccn used to numerical simulation of several geomctrical variants (thickness of hnundary layer o). THE RESULTS OF NUMERICAL SlMULATIONS The solidification process in area of aluminium casting has hccn considered. The following thermophysical parameters for Dm(r) have hcen acccptcd: c 0 =C1, 'A 0 =7>...1, p 0 =p 1, c 1 =1050 [J/kgK], >...,=100 [W/mK], p 1 =2700 [k g/m 3 ], c 2 = 1290, >... 2 =250, p 2 =2380, L=29000(J-p 2 [J/m 3 ], T'=660 C, Tr=700, Tc=650. The function descrihing the physical enthalpy [J/m 3 ] is of the form T< 660 T> 660 ( 11 ) whereas an inverse function T(H) H H< I , H lo H E < , > H > (12) The step of mesh in normai direction to mol d smface i s h = 0. 12[ mm 1 at t he same t i me t he shape of eontroi volume corresponding to r has hcen takcn in to account. In Fig. 3 the kinetics of so1idification for o,-hx>, o2=0.7.'i and 03=0.25 [mml is shown. One can notice that the solidification processin immovahle volume o f molten metal <h -+oo)
6 174 proceeds slowly in comparison with solidification rate for the case o< oo and this rate increases with the rate acceleration of molten metal it means with reduction of boundary layer thickness. The very interesting is a fact that in initial stage of duration of the process (comp. Fig. 3) the phenomena proceeds conversely. It results from very intensive heat abstraction (superheating enthalpy) from liquid area to inteńace r and then the crystallization rate is smaller than in typical solidification conditions. After the essential reduction of temperature in D 0 (t) the process proceeds in expected way. ~[mm] t------ł ł----ł----' t [s] Fig. 3. Kinetics o f solidification l: o,... oo, 2:.5 2 =0. 75mm, 3:.5 3 =0.25mm Very distinct are the qualitative and quantitative differences of temperature fields in considered area (comp. Figures 4, 5, 6). As one can anticipate the temperature in D 0 (t) quickly tends to homogenous field and its value only slightly surpasses the solidification point. Presented above model comprises the numerous simplifications. They as a ruje could be discarded but only by essential extension of algorithm and adequate computer program. It should be pointed that from the other hand the experimental verification in fuli confirmed the effectiveness of proposed numerical model.
7 17 ~ Fig. 4. Temperature field (ó -+oo) e11 ee7 &ej / 1/ l l ~ 1--- e--? V V l-----' ~ - s 1.5s 3s if; V i 655 o ~[mm] Fig. 5. Temperature field (ó=0.75 mm)
8 176 l s l.ss o:: ~6 ~[mm) Fig. 6. Temperature field (ó=0.25 mm) ACKNOWLEDGEMENT This research work has been supported by KBN (Grant No Ol). The author wish to thank for prof. B. Mochnacki for the help in the range of this paper preparing. REFERENCES [l] W.Longa. Solidification of Castings., Śląsk, Katowice 1985 [2] R.Grzymkowski, B.Mochnacki, The Analysis of Solidification Process in Continuous Steel Casting Domain, Solidification of Metais and Alloys, Ossolineum, 2, (1980), pp [3] J.Szargut, Thermal Computations in lndustry Furnaces, Śląsk, Katowice 1977 [ 4] B. Mochnacki, J.S.Suchy, Modeliing and Numerical Simulation of Casting Solidification, PWN Warszawa 1993 [5] A.A.Samarski, Teoria raznostnych schem, Nauka, Moskwa 1983
A R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume Issue 1 DOI: /amm
A R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume 60 205 Issue DOI: 0.55/amm-205-0033 M. RYWOTYCKI, Z. MALINOWSKI, K. MIŁKOWSKA-PISZCZEK, A. GOŁDASZ, B. HADAŁA COMPARISON OF FOUR MODELS
More informationA R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume Issue 1 DOI: /v
A R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume 57 2012 Issue 1 DOI: 10.2478/v10172-012-0030-7 Z. MALINOWSKI, T. TELEJKO, B. HADAŁA INFLUENCE OF HEAT TRANSFER BOUNDARY CONDITIONS
More informationA R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume Issue 2 DOI: /amm
A R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume 59 2014 Issue 2 DOI: 10.2478/amm-2014-0080 M. RYWOTYCKI, Z. MALINOWSKI, J. GIEŁŻECKI, A. GOŁDASZ MODELLING LIQUID STEEL MOTION CAUSED
More informationMODEL OF CAST IRON SOLIDIFICATION USING THE ARTIFICIAL MUSHY ZONE APPROACH
42/22 Archives of Foundry, Year 26, Volume 6, 22 Archiwum Odlewnictwa, Ro 26, Roczni 6, Nr 22 PAN Katowice PL ISSN 1642-538 MODEL OF AST IRON SOLIDIFIATION USING THE ARTIFIIAL MUSHY ZONE APPROAH S.LARA
More informationMODELING OF STEADY AND TRANSIENT TEMPERATURE DISTRIBUTION IN THE DEVICE FOR MEASURING THE THERMAL CONDUCTIVITY
TECHNICAL TRANSACTIONS MECHANICS CZASOPISMO TECHNICZNE MECHANIKA 1-M/2013 PIOTR DUDA*, ROMAN DUDA* MODELING OF STEADY AND TRANSIENT TEMPERATURE DISTRIBUTION IN THE DEVICE FOR MEASURING THE THERMAL CONDUCTIVITY
More informationEFFECT OF MASS MATRIX FORMS ON NUMERICAL SIMULATIONS RESULTS IN HEAT CONDUCTION MODELING. Maria Zych
Journal of Applied Mathematics and Computational Mechanics 2015, 14(3), 149-156 www.amcm.pcz.pl p-issn 2299-9965 DOI: 10.17512/jamcm.2015.3.16 e-issn 2353-0588 EFFECT OF MASS MATRIX FORMS ON NUMERICAL
More informationA R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume Issue 1 DOI: /v
A R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume 58 2013 Issue 1 DOI: 10.2478/v10172-012-0174-5 A. SMALCERZ ASPECTS OF APPLICATION OF INDUSTRIAL ROBOTS IN METALLURGICAL PROCESSES
More informationMODELING OF HEAT FLOW THROUGH MULTILAYER INTERNAL SUPPORTS OF CRYOGENIC VESSELS
WOJCIECH CZYŻYCKI * MODELING OF HEAT FLOW THROUGH MULTILAYER INTERNAL SUPPORTS OF CRYOGENIC VESSELS MODELOWANIE PRZEPŁYWU CIEPŁA PRZEZ WIELOWARSTWOWE PODPORY MIĘDZYPŁASZCZOWE W ZBIORNIKACH KRIOGENICZNYCH
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 informationTable of Contents. Foreword... Introduction...
Table of Contents Foreword.... Introduction.... xi xiii Chapter 1. Fundamentals of Heat Transfer... 1 1.1. Introduction... 1 1.2. A review of the principal modes of heat transfer... 1 1.2.1. Diffusion...
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
More informationNUMERICAL STUDY OF MELTING OF TIN WITHIN A RECTANGULAR CAVITY INCLUDING CONVECTIVE EFFECTS
NUMERICAL STUDY OF MELTING OF TIN WITHIN A RECTANGULAR CAVITY INCLUDING CONVECTIVE EFFECTS Christiano Garcia da Silva Santim, chrisoff22@yahoo.com.br Luiz Fernando Milanez, milanez@fem.unicamp.br Universidade
More informationThermal Analysis measurements
Thermal Analysis measurements R W McCallum Ames Laboratory And Materials Science and Engineering Phase vs Phase Field phase set of states of a macroscopic physical system that have relatively uniform chemical
More informationA. K. M. B. Rashid Professor, Department of MME BUET, Dhaka
A. K. M. B. Rashid Professor, Department of MME BUET, Dhaka 1 Introduction to thermodynamics Power and limitations of thermodynamics The structure of thermodynamics MME6701: Brief lecture format The term
More informationA R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume Issue 4 DOI: /v
A R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume 57 2012 Issue 4 DOI: 10.2478/v10172-012-0124-2 M. MAJ, W. OLIFERUK, ANALYSIS OF PLASTIC STRAIN LOCALIZATION ON THE BASIS OF STRAIN
More informationTHe continuous casting of metals, alloys, semiconductor
A reproduction of boundary conditions in three-dimensional continuous casting problem Iwona Nowak, Jacek Smolka, and Andrzej J. Nowak Abstract The paper discusses a 3D numerical solution of the inverse
More informationVIBRATION EFFECTS ON HEAT TRANSFER DURING SOLIDIFICATION OF PARAFFIN
HEFAT 2012 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics July 16-18, 2012 Malta VIBRATION EFFECTS ON HEAT TRANSFER DURING SOIDIFICATION OF PARAFFIN Johnathan Vadasz,
More informationRelationship to Thermodynamics. Chapter One Section 1.3
Relationship to Thermodynamics Chapter One Section 1.3 Alternative Formulations Alternative Formulations Time Basis: CONSERVATION OF ENERGY (FIRST LAW OF THERMODYNAMICS) An important tool in heat transfer
More informationfirst law of ThermodyNamics
first law of ThermodyNamics First law of thermodynamics - Principle of conservation of energy - Energy can be neither created nor destroyed Basic statement When any closed system is taken through a cycle,
More informationSTUDY OF HEAT TRANSFER MECHANISMS DURING THE LENS TM PROCESS
STUDY OF HEAT TRANSFER MECHANISMS DURING THE LENS TM PROCESS Liang Wang 1 and Sergio Felicelli 1. Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, MS 3976, USA; email:
More informationNumerical Study of a High Temperature Latent Heat Storage ( C) Using NaNO 3 -KNO 3 Binary Mixture
1 Presented at the COMSOL Conference 2010 Paris Numerical Study of a High Temperature Latent Heat Storage (200-300 0 C) Using NaNO 3 -KNO 3 Binary Mixture Foong Chee Woh, 17-11-2010 2 Background 3 Background
More informationCoupled Thermomechanical Contact Problems
Coupled Thermomechanical Contact Problems Computational Modeling of Solidification Processes C. Agelet de Saracibar, M. Chiumenti, M. Cervera ETS Ingenieros de Caminos, Canales y Puertos, Barcelona, UPC
More informationLaser processing of materials. Temperature distributions
Laser processing of materials Temperature distributions Prof. Dr. Frank Mücklich Dr. Andrés Lasagni Lehrstuhl für Funktionswerkstoffe Sommersemester 7 Contents: Temperature distributions 1. Definitions.
More informationChapter 10: Steady Heat Conduction
Chapter 0: Steady Heat Conduction In thermodynamics, we considered the amount of heat transfer as a system undergoes a process from one equilibrium state to another hermodynamics gives no indication of
More informationCHEMISTRY 12 UNIT II EQUILIBRIUM D Learning Goals
CHEMISTRY 12 UNIT II EQUILIBRIUM D Learning Goals 1. Chemical equilibrium is said to by dynamic because a. The reaction proceeds quickly b. The mass of the reactants is decreasing c. The macroscopic properties
More informationReview of Mathematics
Review for Exam #1 Review of Mathematics 2 Weighted Mean A certain property of material 1 is P 1 and that of material 2 is P 2 If x 1 amount (weight or volume) of material 1 is mixed with x 2 amount of
More informationVolume Issue 3 DOI: /amm Effect of turbulence modelling in numerical analysis of melting process in an induction furnace
A R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume 60 2015 Issue 3 DOI: 10.1515/amm-2015-0275 P. BULIŃSKI*,#, J. SMOŁKA*, S. GOLAK**, R. PRZYŁUCKI**, L. BLACHA***, R. BIAŁECKI*, M.
More information6. Transport phenomena. Basel, 2008
6. Transport phenomena Basel, 2008 6. Transport phenomena 1. Introduction 2. Phenomenological equations for transport properties 3. Transport properties of a perfect gas 4. Diffusion in a fluid 5. Measurement
More informationFigure Q3. Boundary conditions for a tubular reactor.
1. Please illustrate or explain the meaning the following questions: (a) (10%) Please proof the five dimensionless numbers in fluid transport by Buckingham method as Eu=f(Re, Fr, We, Ma). (b) (10%) Please
More informationFLOW MODELLING OF SLURRY ICE IN A CONTROL VALVE MODELOWANIE PRZEPŁYWU ZAWIESINY LODOWEJ W ZAWORZE REGULACYJNYM
TECHNICAL TRANSACTIONS MECHANICS CZASOPISMO TECHNICZNE MECHANIKA 1-M/2015 ŁUKASZ MIKA*, BOGUSŁAW GÓRSKI**, RYSZARD KANTOR* FLOW MODELLING OF SLURRY ICE IN A CONTROL VALVE MODELOWANIE PRZEPŁYWU ZAWIESINY
More informationPareto - ABC Analysis of Temperature Field in High Voltage Three-Phase Cable Systems
Janusz TYKOCKI 1, Andrzej JORDAN 2, Dariusz SUROWIK 3 Lomza State University of Applied Sciences (1)(2), Bialystok University (3) doi:10.15199/48.2015.05.26 Pareto - ABC Analysis of Temperature Field in
More information3. First Law of Thermodynamics and Energy Equation
3. First Law of Thermodynamics and Energy Equation 3. The First Law of Thermodynamics for a ontrol Mass Undergoing a ycle The first law for a control mass undergoing a cycle can be written as Q W Q net(cycle)
More informationSENSITIVITY OR TEMPERATURE FIELD IN THE SYSTEM PROTECTIVE CLOTHING FOREARM WITH RESPECT TO PERTURBATIONS OF EXTERNAL HEATING CONDITIONS
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationUnits. specified to be meaningful. between different scales. - All physical quantities must have their units. - We will always used SI units.
Units - All physical quantities must have their units specified to be meaningful. - We will always used SI units. - We must be comfortable with conversions between different scales. 2 General View of the
More informationHeat Transfer Modeling using ANSYS FLUENT
Lecture 1 - Introduction 14.5 Release Heat Transfer Modeling using ANSYS FLUENT 2013 ANSYS, Inc. March 28, 2013 1 Release 14.5 Outline Modes of Heat Transfer Basic Heat Transfer Phenomena Conduction Convection
More information- Apply closed system energy balances, observe sign convention for work and heat transfer.
CHAPTER : ENERGY AND THE FIRST LAW OF THERMODYNAMICS Objectives: - In this chapter we discuss energy and develop equations for applying the principle of conservation of energy. Learning Outcomes: - Demonstrate
More information1) The unit of the coefficient of viscosity in SI system is: a) m/kg-s b) m-s/kg 2 c) kg/m-s 2 d) kg/m-s
1) The unit of the coefficient of viscosity in SI system is: a) m/kg-s b) m-s/kg 2 c) kg/m-s 2 d) kg/m-s 2) Joule x see is unit of a) energy b) momentum c) angular momentum d) Power 3) The unit of surface
More informationCorresponding author.
EFFECTS OF BEAM SIZE AND PULSE DURATION ON THE LASER DRILLING PROCESS Nazia Afrin, Pengfei Ji, J. K. Chen and Yuwen Zhang 1 Department of Mechanical and Aerospace Engineering University of Missouri Columbia,
More informationNUMERICAL AND EXPERIMENTAL COLLAPSE ANALYSIS OF TUBULAR MULTI-MEMBER ENERGY ABSORBERS UNDER LATERAL COMPRESSION. Notation
JOURNAL OF THEORETICAL AND APPLIED MECHANICS 42, 4, pp. 841-857, Warsaw 2004 NUMERICAL AND EXPERIMENTAL COLLAPSE ANALYSIS OF TUBULAR MULTI-MEMBER ENERGY ABSORBERS UNDER LATERAL COMPRESSION Sebastian Lipa
More informationSOLIDIFICATION SURFACE SPEED CONTROL OF FERROMAGNETIC PIECES USING EDDY CURRENT HEATING
SOLIDIFICATION SURFACE SPEED CONTROL OF FERROMAGNETIC PIECES USING EDDY CURRENT HEATING MIHAI MARICARU, MARILENA STĂNCULESCU, 1 VALERIU ŞTEFAN MINCULETE, 1 FLOREA IOAN HĂNŢILĂ 11 Key words: Coupled eddy
More informationThe experimental determination of the thermal conductivity of melting chocolate: thermal resistance analogies and free convection boundary conditions
Advanced Computational Methods and Experiments in Heat Transfer XIII 505 The experimental determination of the thermal conductivity of melting chocolate: thermal resistance analogies and free convection
More information3.012 PS Issued: Fall 2004 Due: pm
3.012 PS 2 3.012 Issued: 09.15.04 Fall 2004 Due: 09.22.04 5pm Graded problems: 1. In discussing coordination numbers and deriving the permitted range of radius ratio, R A / R B, allowed for each ( where
More informationCOEFFICIENTS OF VISCOSITY FOR A FLUID IN A MAGNETIC FIELD OR IN A ROTATING SYSTEM
Hooyman, G.J. Physica XXI Mazur, P. " 355-359 Groot, S. R. de 1955 COEFFICIENTS OF VISCOSITY FOR A FLUID IN A MAGNETIC FIELD OR IN A ROTATING SYSTEM by G. J. HOOYMAN, P. MAZUR and S. R. DE GROOT Instituut
More information4.1. Physics Module Form 4 Chapter 4 - Heat GCKL UNDERSTANDING THERMAL EQUILIBRIUM. What is thermal equilibrium?
4.1 4 UNDERSTANDING THERMAL EQUILIBRIUM What is thermal equilibrium? 1. ( Heat, Temperature ) is a form of energy that flows from a hot body to a cold body. 2. The SI unit for ( heat, temperature) is Joule,
More informationSimilarities and differences:
How does the system reach equilibrium? I./9 Chemical equilibrium I./ Equilibrium electrochemistry III./ Molecules in motion physical processes, non-reactive systems III./5-7 Reaction rate, mechanism, molecular
More informationEffectiveness NTU performance of finned PCM storage unit K.A.R. Ismail, M.M. Gongalves
Effectiveness NTU performance of finned PCM storage unit K.A.R. Ismail, M.M. Gongalves Campinas, SP, Brazil ABSTRACT This paper presents a mathematical model based upon two dimensional formulation of the
More informationThe effect of natural convection on solidification in tall tapered feeders
ANZIAM J. 44 (E) ppc496 C511, 2003 C496 The effect of natural convection on solidification in tall tapered feeders C. H. Li D. R. Jenkins (Received 30 September 2002) Abstract Tall tapered feeders (ttfs)
More informationENSC 388. Assignment #8
ENSC 388 Assignment #8 Assignment date: Wednesday Nov. 11, 2009 Due date: Wednesday Nov. 18, 2009 Problem 1 A 3-mm-thick panel of aluminum alloy (k = 177 W/m K, c = 875 J/kg K, and ρ = 2770 kg/m³) is finished
More informationModelling of Solidification and Melting in a Latent Heat Storage
Modelling of Solidification and Melting in a Latent Heat Storage A Quasi-Stationary Approach Felix Eckl, Simon Maranda, Anastasia Stamatiou, Ludger Fischer, Jörg Worlitschek Lucerne University of Applied
More informationAutumn 2005 THERMODYNAMICS. Time: 3 Hours
CORK INSTITUTE OF TECHNOOGY Bachelor of Engineering (Honours) in Mechanical Engineering Stage 3 (Bachelor of Engineering in Mechanical Engineering Stage 3) (NFQ evel 8) Autumn 2005 THERMODYNAMICS Time:
More informationNUMERICAL SIMULATION OF ELECTROWETTING USING ELECTRIC BODY FORCE DENSITY AND SURFACE TENSION
Tan Il SUNG NUMERICAL SIMULATION OF ELECTROWETTING USING ELECTRIC BODY FORCE DENSITY AND SURFACE TENSION ABSTRACT The electrowetting is a phenomenon in which the shape of a droplet is deformed by external
More informationDevelopment of an evaporation boundary condition for DSMC method with application to meteoroid entry
Development of an evaporation boundary condition for DSMC method with application to meteoroid entry F. Bariselli, S. Boccelli, A. Frezzotti, A. Hubin, T. Magin Annual METRO meeting 29th November 2016
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 informationSensible Heat and Enthalpy Calculations
Sensible Heat and Enthalpy Calculations Sensible Heat - The amount of heat that must be added when a substance undergoes a change in temperature from 298 K to an elevated temperature without a change in
More informationSensible Heat and Enthalpy Calculations
* Sensible Heat and Enthalpy Calculations Sensible Heat - The amount of heat that must be added when a substance undergoes a change in temperature from 298 K to an elevated temperature without a change
More informationENTROPY OF MELTING AND STRUCTURAL CONDITION OF A MELT 12,, S. PROKHORENKO 2 V. PROKHORENKO. Rzeszów University, Poland Rzeszów 2
22/6 ARCHIWUM ODLEWNICTWA Rok 2002, Rocznik 2, Nr 6 Archives of Foundry Year 2002, Volume 2, Book 6 PAN - Katowice PL ISSN 1642-5308 ENTROPY OF MELTING AND STRUCTURAL CONDITION OF A MELT V. PROKHORENKO
More informationHEAT TRANSFER UNDER AN AIR-WATER COOLING JET
HEAT TRANSFER UNDER AN AIR-WATER COOLING JET F. Kavicka 1, J. Stetina 1, B. Sekanina 1, J. Heger 2, B. Velicka 3, L. Čamek 4, V. Gontarev 5, B. Kosec 5 1 VUT - FSI, Brno, Technická 2, Brno 616 69, Česká
More information3-D Turbulent Multiphase Modeling of Molten Steel Flow and Heat Transfer in a Continuous Slab Caster. Abstract
Presented at FX User's onference, Wilmington, DE, Oct. 1, 1998. 3-D Turbulent Multiphase Modeling of Molten Steel Flow and Heat Transfer in a ontinuous Slab aster D.T. reech and B.G. Thomas, University
More informationFinite Element Analysis of the Heat Transfer in a Copper Mould during Continuous Casting of Steel Slabs. 14 May 2005
Finite Element Analysis of the Heat Transfer in a Copper Mould during Continuous Casting of Steel Slabs 14 May 2005 D. Hodgson, Sami Vapalahti, and B.G. Thomas Department of Mechanical and Industrial Engineering
More informationA R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume Issue 3
A R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume 55 2010 Issue 3 M. KIELOCH, Ł. PIECHOWICZ, J. BORYCA, A. KLOS NUMERICAL ANALYSIS OF CORRELATION BETWEEN HEAT CONSUMPTION AND THE STEEL
More informationThe Structure of Thermodynamics
MME6701: Lecture 02 The Structure of Thermodynamics A. K. M. B. Rashid Professor, Department of MME BUET, Dhaka Topics to discuss q Thermodynamic Systems Classification of thermodynamic systems q Thermodynamic
More information39th International Physics Olympiad - Hanoi - Vietnam Experimental Problem
DIFFERENTIAL THERMOMETRIC METHOD In this problem, we use the differential thermometric method to fulfill the two following tasks: 1. Finding the temperature of solidification of a crystalline solid substance.
More informationConduction Heat Transfer. Fourier Law of Heat Conduction. x=l Q x+ Dx. insulated x+ Dx. x x. x=0 Q x A
Conduction Heat Transfer Reading Problems 10-1 10-6 10-20, 10-48, 10-59, 10-70, 10-75, 10-92 10-117, 10-123, 10-151, 10-156, 10-162 11-1 11-2 11-14, 11-20, 11-36, 11-41, 11-46, 11-53, 11-104 Fourier Law
More informationAnalysis of Forced Convection Heat Transfer in Microencapsulated Phase Change Material Suspensions
JOURNAL OF THERMOPHYSICS AND HEAT TRANSFER Vol. 9, No. 4, October-December 1995 Analysis of Forced Convection Heat Transfer in Microencapsulated Phase Change Material Suspensions Yuwen Zhangh and Amir
More informationAdvanced Subsidiary Unit 1: The Core Principles of Chemistry
Write your name here Surname Other names Pearson Edexcel International Advanced Level Centre Number Chemistry Advanced Subsidiary Unit 1: The Core Principles of Chemistry Candidate Number Thursday 9 January
More informationThermal Capacity Measurement of Engineering Alloys in Dependence on Temperature.
Thermal Capacity Measurement of Engineering Alloys in Dependence on Temperature. Z. Jedlicka, 1 I. Herzogova 2 ABSTRACT Thermophysical laboratory at Department of Thermal Engineering of Technical University
More informationTransactions on Modelling and Simulation vol 10, 1995 WIT Press, ISSN X
Numerical analysis of melting process in a rectangular enclosure J. Goscik Department of Mechanical Engineering, Technical University of Bialystok, 15-351 Bialystok, Wiejska 45C, Poland Abstract In the
More informationASSUMPTIONS: (1) One-dimensional, radial conduction, (2) Constant properties.
PROBLEM 5.5 KNOWN: Diameter and radial temperature of AISI 00 carbon steel shaft. Convection coefficient and temperature of furnace gases. FIND: me required for shaft centerline to reach a prescribed temperature.
More informationModule 16. Diffusion in solids II. Lecture 16. Diffusion in solids II
Module 16 Diffusion in solids II Lecture 16 Diffusion in solids II 1 NPTEL Phase II : IIT Kharagpur : Prof. R. N. Ghosh, Dept of Metallurgical and Materials Engineering Keywords: Micro mechanisms of diffusion,
More informationLaw of Heat Transfer
Law of Heat Transfer The Fundamental Laws which are used in broad area of applications are: 1. The law of conversion of mass 2. Newton s second law of motion 3. First and second laws of thermodynamics
More information1 One-Dimensional, Steady-State Conduction
1 One-Dimensional, Steady-State Conduction 1.1 Conduction Heat Transfer 1.1.1 Introduction Thermodynamics defines heat as a transfer of energy across the boundary of a system as a result of a temperature
More informationTHERMAL ANALYSIS OF A SPENT FUEL TRANSPORTATION CASK
Excerpt from the Proceedings of the COMSOL Conference 2009 Bangalore THERMAL ANALYSIS OF A SPENT FUEL TRANSPORTATION CASK P. Goyal*, Vishnu Verma, R.K. Singh & A.K. Ghosh Reactor Safety Division Bhabha
More informationMinistry of Higher Education And Scientific Research. University Of Technology Chemical Engineering Department. Heat Transfer
Ministry of Higher Education And Scientific Research University Of Technology Heat Transfer Third Year By Dr.Jamal Al-Rubeai 2008-2009 Heat Transfer 1. Modes of Heat Transfer: Conduction, Convection and
More informationTHEORY OF DISTRIBUTIONS
THEORY OF DISTRIBUTIONS THE SEQUENTIAL APPROACH by PIOTR ANTOSIK Special Research Centre of the Polish Academy of Sciences in Katowice JAN MIKUSltfSKI Special Research Centre of the Polish Academy of Sciences
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 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 informationBasic Thermodynamics Prof. S. K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur. Lecture No 16
Basic Thermodynamics Prof. S. K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No 16 Properties of Pure Substances-I Good afternoon. In the last class, we were
More informationPTT 277/3 APPLIED THERMODYNAMICS SEM 1 (2013/2014)
PTT 77/3 APPLIED THERMODYNAMICS SEM 1 (013/014) 1 Energy can exist in numerous forms: Thermal Mechanical Kinetic Potential Electric Magnetic Chemical Nuclear The total energy of a system on a unit mass:
More informationTHE INFLUENCE OF GLAZING TYPE, FRAME PROFILES SHAPE AND SIZE OF THE WINODOW OF THE FINAL VALUE OF WINDOW THERMAL TRANSMITTANCE U
MAŁGORZATA ROJEWSKA-WARCHAŁ * THE INFLUENCE OF GLAZING TYPE, FRAME PROFILES SHAPE AND SIZE OF THE WINODOW OF THE FINAL VALUE OF WINDOW THERMAL TRANSMITTANCE U ANALIZA WPŁYWU PRZESZKLENIA, RAM ORAZ KSZTAŁTU
More informationA Mixed Finite Element Formulation for Solving Phase Change Problems with Convection
A Mixed Finite Element Formulation for Solving Phase Change Problems with Convection Youssef Belhamadia 1, Abdoulaye S. Kane 2, and André Fortin 3 1 University of Alberta, Campus Saint-Jean and Department
More informationIdentification of model parameters of a sandwich beam incorporating magnetorheological fluid
XXIV Symposium Vibrations in Physical Systems, Poznan Bedlewo, May 12-15, 2010 Identification of model parameters of a sandwich beam incorporating magnetorheological fluid Bogdan SAPIŃSKI AGH University
More informationChapter 1 INTRODUCTION AND BASIC CONCEPTS
Heat and Mass Transfer: Fundamentals & Applications 5th Edition in SI Units Yunus A. Çengel, Afshin J. Ghajar McGraw-Hill, 2015 Chapter 1 INTRODUCTION AND BASIC CONCEPTS Mehmet Kanoglu University of Gaziantep
More informationCombining Multiphysics Modeling and Solution Thermodynamics Using M4Dlib, an External Library
Combining Multiphysics Modeling and Solution Thermodynamics Using M4Dlib, an External Library Tanai L. Marin-Alvarado *1 1 M4Dynamics Inc. *Corresponding author: 1120 Finch Ave. W. Suite 701, Toronto,
More informationMathematical Modelling of Ceramic Block Heat Transfer Properties
Proceedings of the 3 RD INTERNATIONAL CONFERENCE ADVANCED CONSTRUCTION 18-19 October, 2012, Kaunas, Lithuania Kaunas University of Technology, Faculty of Civil Engineering and Architecture Studentu st.
More informationHEAT- I Part - A C D A B. Te m p. Heat input
e m p HE- I Part -. solid material is supplied with heat at a constant rate. he temperature of the material is changing with heat input as shown in the graph. Study the graph carefully and answer the following
More informationheat transfer process where a liquid undergoes a phase change into a vapor (gas)
Two-Phase: Overview Two-Phase two-phase heat transfer describes phenomena where a change of phase (liquid/gas) occurs during and/or due to the heat transfer process two-phase heat transfer generally considers
More informationMultiscale modelling of microstructure formation in polymer casting
TECHNISCHE MECHANIK, 30, 1-3, (2010), 259 268 submitted: October 31, 2009 Multiscale modelling of microstructure formation in polymer casting T. Pitkänen, S. Majaniemi, T. Ala-Nissila A data bank approach
More informationGOZO COLLEGE BOYS SECONDARY SCHOOL
GOZO COLLEGE BOYS SECONDARY SCHOOL Embracing Diversity Half Yearly Exams 2013-2014 FORM 4 CHEMISTRY TIME: 1h 30min Name: Class: Useful Data: One Faraday is equivalent to 96500 C. 1 2 3 4 5 6 7 0 1 H 1
More informationChapter 2 Heat, Temperature and the First Law of Thermodynamics
Chapter 2 Heat, Temperature and the First Law of Thermodynamics 2.1. Temperature and the Zeroth Law of Thermodynamics 2.2. Thermal Expansion 2.3. Heat and the Absorption of Heat by Solids and Liquids 2.4.
More informationChapter 17 Temperature and heat
Chapter 17 Temperature and heat 1 Temperature and Thermal Equilibrium When we speak of objects being hot and cold, we need to quantify this by some scientific method that is quantifiable and reproducible.
More informationA R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume Issue 3
A R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume 52 2007 Issue 3 A. FEDORYSZYN P. ZYZAK CHARACTERISTICS OF THE ABRASIVE STREAM DEPENDING ON DESIGN FEATURES AND OPERATIONAL PARAMETERS
More informationTemperature and Heat. Ken Intriligator s week 4 lectures, Oct 21, 2013
Temperature and Heat Ken Intriligator s week 4 lectures, Oct 21, 2013 This week s subjects: Temperature and Heat chapter of text. All sections. Thermal properties of Matter chapter. Omit the section there
More informationLiquid. T > Tm Liquid has. Solid T < Tm Solid has. the lower free energy T. Demo. the lower free energy. Solutions.
Just to be clear about Free Energy Super Cooled or Super Heated G = H - TS straight line assumes that H and S are independent of temperature Slope is given by S Liquid has a larger entropy and therefore
More informationAP PHYSICS 2 WHS-CH-14 Heat Show all your work, equations used, and box in your answers! 1 108kg
AP PHYSICS 2 WHS-CH-4 Heat Show all your work, equations used, and box in your answers! James Prescott Joule (88 889) James Prescott Joule studied the nature of heat, and discovered its relationship to
More informationNumerical analysis of natural convection in a latent heat thermal energy storage system containing rectangular enclosures
EUROTHERM99-01-074 Numerical analysis of natural convection in a latent heat thermal energy storage system containing rectangular enclosures Julian Vogel, Maike Johnson, Markus Eck, Dörte Laing German
More informationDO NOT TURN OVER UNTIL TOLD TO BEGIN
UIVERSITY OF LODO MSci EXAMIATIO 007 For Students of the University of London DO OT UTIL TOLD TO BEGI PH411A: STATISTICAL MECHAICS Time Allowed: TWO AD A HALF hours Answer THREE questions only. o credit
More informationIntroduction to the Thermodynamics of Materials
Solutions Manual to accompany Introduction to the hermodynamics of Materials Sixth Edition David R. Gaskell School of Materials Engineering Purdue University West Lafayette, IN David E. Laughlin ALCOA
More informationINSTRUCTIONS: CHEM Exam I. September 13, 1994 Lab Section
CHEM 1314.05 Exam I John I. Gelder September 13, 1994 Name TA's Name Lab Section Please sign your name below to give permission to post, by the last 4 digits of your student I.D. number, your course scores
More informationA R C H I V E S O F M E T A L L U R G Y A N D M A T E R I A L S Volume Issue 4 DOI: /v
A R C H I V E S O F M E A L L U R G Y A N D M A E R I A L S Volume 57 2012 Issue 4 DOI: 10.2478/v10172-012-0133-1. SKRZYPCZAK SHARP INERFACE NUMERICAL MODELING OF SOLIDIFICAION PROCESS OF PURE MEAL SPOSÓB
More informationHeat Transfer Benchmark Problems Verification of Finite Volume Particle (FVP) Method-based Code
PROCEEDING OF 3 RD INTERNATIONAL CONFERENCE ON RESEARCH, IMPLEMENTATION AND EDUCATION OF MATHEMATICS AND SCIENCE YOGYAKARTA, 16 17 MAY 2016 Heat Transfer Benchmark Problems Verification of Finite Volume
More information