Model of Heat and Mass Transfer with Moving Boundary during Roasting of Meat in Convection-Oven

Size: px
Start display at page:

Download "Model of Heat and Mass Transfer with Moving Boundary during Roasting of Meat in Convection-Oven"

Transcription

1 Presented at the COMSOL Conference 2009 Milan Model of Heat and Mass Transfer with Moving Boundary during Roasting of Meat in Convection-Oven COMSOL Conference October 15, 2009, Milan, ITALY Aberham Hailu Feyissa, PhD student Jens Adler-Nissen, Prof, dr. techn. Krist V. Gernaey, Assoc.Prof, Ph.D Food Production Engineering

2 Outline Introduction Theoretical background Modelling Results Conclusion Governing model equations Boundary and initial condition Model Solution 2 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

3 Roasting process Model Better quantitative knowledge Prediction of T, C Quality, safety Minimize loss Upscale Mass and energy Process control Transport processes during oven convection roasting 3 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

4 Water transport Pure diffusion (often assumed) But doesn t capture pressure Convective T Denaturation WHC Shrinkage (protein network) Pressure gradient velocity 4 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

5 Change of microstructure Permeability(K) Elastic modulus (E) Raw surface o P i Case-1 center T- T- low high Roasted surface o T-high R-low P i Case-2 T- low R -high center Change of microstructure during roasting (cross-sectional view) 5 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

6 Assumptions Fat transport is negligible g (lean meat) Evaporation takes place at the surface No internal heat generation and no chemical reaction Dissolved matter lost with water is small Reduction of water holding capacity and shrinkage are considered. 6 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

7 , Governing equation of heat and mass transfer Heat transfer T ρ m cp m + ( km T) + ρwc p uw T = 0 t w Water transport ρ m = D = f c pm = C + Cu w = D C ( ) Emx t E ( T ) = E o + (1+ exp( E ( T E ) ) Velocity of water µ w = n f ( T ) f f ( C ) ( C ) ( T ) D u w = K µ w P P= EC ( Ceq(T) ) u w KE = C µ w ( C C ) eq Darcy s law C eq ( T ) = 0.75 (1 + 30exp( 0.25( T Tσ ))) 7 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

8 Shrinkage hi /interface velocity Proportional to the volume of liquid water removed Formation of air filled pores V = V 0 β V w, l β=1, no pore formation β = 0, no shrinkage 1-β= fraction of Vw,l replaced by air V = πr βv w, V0 1 V / 3 1/ 3 w, l βv w, l Z 0 1 V 0 V0 2/3 0β Vw, l d 1 ( Vw l ) V 0 V 0 dt βv 1 V l = m ( X X) ρ V ( 1 C ) dz Z β v z = =, dt 3 2/ 3 dr R0β βvw, l d v r = = 1, dt 3V 0 V 0 dt ( V ) w l = πr 2. Z V, = w l d 0 ρ w dvw, dt l = 0 0 ρ ρ0v 0 = ρ w 0 ( 1 C ) w 0 C 0 1 C C Cav 1 C av 2 av dc dt av 8 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

9 Boundary conditions and initial condition BC 1: Axial symmetry, vr=0 BC (2, 3, and 4): HT- heat flux MT- mass flux ALE- velocity IC: T C ( r, z ) = T 0 = const at t = 0 ( r, z ) = C 0 = const at t = 0 9 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

10 Solution COMSOL Multiphyics version3.5 2D cylindrical Radius of 20 mm and length of 54 mm Chemical Engineering module Transient heat transfer Transient mass transfer Moving mesh module (ALE) 10 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

11 Results Temperature distribution Water content distribution 11 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

12 t=500s Temperature profile across cylindrical sample (Z = 0) Iso-concentration lines C (kg/kg) at t = 500 s 12 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

13 Model prediction i with (-) MB vs. with FB (--) Temperature profile blue-center(0, 0) green-surface (r = 0.02, z = 0) Water content profile blue(0, 0), red (0.017, 0) green (0.019, 0), cyan (0.02, 0) 13 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

14 Radial shrinkage I II III Relative radius of cylinder as function of time (R/R0) 14 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

15 Conclusion A 2D model of HMT with MB is developed and models equations were solved using COMSOL for convection roasting process. Effect of WHC and shrinkage are incorporated. Better insight i is obtained. ACKNOWLEDGEMENT DTU for a Ph.D. Grant (Globalization Fund). 15 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

16 Thank you for your attention! 16 Food Production Engineering, DTU Fødevareinstituttet, Danmarks Tekniske Universitet

ASSUMPTIONS: (1) Homogeneous medium with constant properties, (2) Negligible radiation effects.

ASSUMPTIONS: (1) Homogeneous medium with constant properties, (2) Negligible radiation effects. PROBEM 5.88 KNOWN: Initial temperature of fire clay bric which is cooled by convection. FIND: Center and corner temperatures after 50 minutes of cooling. ASSUMPTIONS: () Homogeneous medium with constant

More information

I. Borsi. EMS SCHOOL ON INDUSTRIAL MATHEMATICS Bedlewo, October 11 18, 2010

I. Borsi. EMS SCHOOL ON INDUSTRIAL MATHEMATICS Bedlewo, October 11 18, 2010 : an : an (Joint work with A. Fasano) Dipartimento di Matematica U. Dini, Università di Firenze (Italy) borsi@math.unifi.it http://web.math.unifi.it/users/borsi porous EMS SCHOOL ON INDUSTRIAL MATHEMATICS

More information

Outline. Definition and mechanism Theory of diffusion Molecular diffusion in gases Molecular diffusion in liquid Mass transfer

Outline. Definition and mechanism Theory of diffusion Molecular diffusion in gases Molecular diffusion in liquid Mass transfer Diffusion 051333 Unit operation in gro-industry III Department of Biotechnology, Faculty of gro-industry Kasetsart University Lecturer: Kittipong Rattanaporn 1 Outline Definition and mechanism Theory of

More information

QUESTION ANSWER. . e. Fourier number:

QUESTION 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 information

q x = k T 1 T 2 Q = k T 1 T / 12

q x = k T 1 T 2 Q = k T 1 T / 12 Conductive oss through a Window Pane q T T 1 Examine the simple one-dimensional conduction problem as heat flow through a windowpane. The window glass thickness,, is 1/8 in. If this is the only window

More information

OPTIMAL Developing a decision support tool for optimization of process and product quality in fish

OPTIMAL Developing a decision support tool for optimization of process and product quality in fish Report 4/2017 Published February 2017 OPTIMAL Developing a decision support tool for optimization of process and product quality in fish Izumi Sone (Nofima), Marthe J. Blikra (Nofima), Jörn Schulz (University

More information

Modeling of the Bread Baking Process Using Moving Boundary and Arbitrary-Lagrangian-Eulerian (ALE) C. Anandharamakrishnan, N. Chhanwal, P.

Modeling of the Bread Baking Process Using Moving Boundary and Arbitrary-Lagrangian-Eulerian (ALE) C. Anandharamakrishnan, N. Chhanwal, P. Presented at the COMSOL Conference 2010 India Modeling of the Bread Baking Process Using Moving Boundary and Arbitrary-Lagrangian-Eulerian (ALE) Approaches C. Anandharamakrishnan, N. Chhanwal, P. Karthik,

More information

Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties. SIF file : phasechange solid-solid

Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties. SIF file : phasechange solid-solid Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties 3 6 1. Tb=1750 [K] 2 & 5. q=-10000 [W/m²] 0,1 1 Ω1 4 Ω2 7 3 & 6. α=15 [W/(m²K)] Text=300 [K] 4.

More information

Chapter 2: Heat Conduction. Dr Ali Jawarneh Department of Mechanical Engineering, Hashemite University

Chapter 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 information

Perforation Effect on a Rectangular Metal Hydride Tank at the Hydriding and Dehydriding Process

Perforation Effect on a Rectangular Metal Hydride Tank at the Hydriding and Dehydriding Process Comsol Conference 2012, Milan Perforation Effect on a Rectangular Metal Hydride Tank at the Hydriding and Dehydriding Process E.I Gkanas, S.S Makridis, E. S. Kikkinides and A.K. Stubos http://h2matters.weebly.com

More information

How are calculus, alchemy and forging coins related?

How are calculus, alchemy and forging coins related? BMOLE 452-689 Transport Chapter 8. Transport in Porous Media Text Book: Transport Phenomena in Biological Systems Authors: Truskey, Yuan, Katz Focus on what is presented in class and problems Dr. Corey

More information

Separation through Dialysis

Separation through Dialysis Separation through Dialysis SOLVED WITH COMSOL MULTIPHYSICS 3.5a COPYRIGHT 2008. All right reserved. No part of this documentation may be photocopied or reproduced in any form without prior written consent

More information

2. Modeling of shrinkage during first drying period

2. Modeling of shrinkage during first drying period 2. Modeling of shrinkage during first drying period In this chapter we propose and develop a mathematical model of to describe nonuniform shrinkage of porous medium during drying starting with several

More information

Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore

Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore Lecture No. # 26 Problem solving : Heterogeneous reactions Friends, in last few

More information

Modeling Flow and Deformation during Salt- Assisted Puffing of Single Rice Kernels. Tushar Gulati Cornell University 1

Modeling Flow and Deformation during Salt- Assisted Puffing of Single Rice Kernels. Tushar Gulati Cornell University 1 Modeling Flow and Deformation during Salt- Assisted Puffing of Single Rice Kernels Tushar Gulati Cornell University 1 Puffing 2 Puffing: Salient Features Quick process and a complex interplay of mass,

More information

Numerical Investigation on The Convective Heat Transfer Enhancement in Coiled Tubes

Numerical Investigation on The Convective Heat Transfer Enhancement in Coiled Tubes Numerical Investigation on The Convective Heat Transfer Enhancement in Coiled Tubes Luca Cattani Department of Industrial Engineering - University of Parma Excerpt from the Proceedings of the 2012 COMSOL

More information

Chapter 10: Steady Heat Conduction

Chapter 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 information

Chapter 5 Time-Dependent Conduction

Chapter 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 information

Fuel - Coolant Heat Transfer

Fuel - Coolant Heat Transfer Heat Transfer 5-1 Chapter 5 Fuel - Coolant Heat Transfer 5.1 Introduction The interface between the fuel and the coolant is centrally important to reactor design since it is here that the limit to power

More information

Flow toward Pumping Well, next to river = line source = constant head boundary

Flow toward Pumping Well, next to river = line source = constant head boundary Flow toward Pumping Well, next to river = line source = constant head boundary Plan view River Channel after Domenico & Schwartz (1990) Line Source Leonhard Euler 1707-1783 e i" +1 = 0 wikimedia.org Charles

More information

University of Rome Tor Vergata

University 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 information

Lecture 9 Laminar Diffusion Flame Configurations

Lecture 9 Laminar Diffusion Flame Configurations Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.

More information

Mathematical Modeling of Oil Shale Pyrolysis

Mathematical Modeling of Oil Shale Pyrolysis October, 19 th, 2011 Mathematical Modeling of Oil Shale Pyrolysis Pankaj Tiwari Jacob Bauman Milind Deo Department of Chemical Engineering University of Utah, Salt Lake City, Utah http://from50000feet.wordpress.com

More information

Solid Food Pasteurization by Ohmic Heating: Influence of Process Parameters

Solid Food Pasteurization by Ohmic Heating: Influence of Process Parameters Excerpt from the Proceedings of the COMSOL Conference 2008 Boston Solid Food Pasteurization by Ohmic Heating: Influence of Process Parameters Markus Zell 1, Denis A. Cronin 1, Desmond J. Morgan 1, Francesco

More information

TEMPERATURE DISTRIBUTION OF AN INFINITE SLAB UNDER POINT HEAT SOURCE

TEMPERATURE DISTRIBUTION OF AN INFINITE SLAB UNDER POINT HEAT SOURCE THERMAL SCIENCE, Year 14, Vol. 18, No. 5, pp. 1597-161 1597 TEMPERATURE DISTRIBUTION OF AN INFINITE SLAB UNDER POINT HEAT SOURCE by Zhao-Chun WU * and Dao-Lai CHENG School of Urban Construction and Safety

More information

UNIVERSITY OF SOUTHAMPTON VERY IMPORTANT NOTE. Section A answers MUST BE in a separate blue answer book. If any blue

UNIVERSITY OF SOUTHAMPTON VERY IMPORTANT NOTE. Section A answers MUST BE in a separate blue answer book. If any blue UNIVERSITY OF SOUTHAMPTON PHYS1013W1 SEMESTER 2 EXAMINATION 2011/12 ENERGY AND MATTER SOLUTIONS Duration: 120 MINS VERY IMPORTANT NOTE Section A answers MUST BE in a separate blue answer book. If any blue

More information

PROBLEM 1.2 ( ) 25 C 15 C dx L 0.30 m Ambient air temperature, T2 (C)

PROBLEM 1.2 ( ) 25 C 15 C dx L 0.30 m Ambient air temperature, T2 (C) PROBLEM 1.2 KNOWN: Inner surface temperature and thermal conductivity of a concrete wall. FIND: Heat loss by conduction through the wall as a function of ambient air temperatures ranging from -15 to 38

More information

Part I.

Part I. Part I bblee@unimp . Introduction to Mass Transfer and Diffusion 2. Molecular Diffusion in Gasses 3. Molecular Diffusion in Liquids Part I 4. Molecular Diffusion in Biological Solutions and Gels 5. Molecular

More information

Heat processes. Heat exchange

Heat 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 information

Clouds associated with cold and warm fronts. Whiteman (2000)

Clouds associated with cold and warm fronts. Whiteman (2000) Clouds associated with cold and warm fronts Whiteman (2000) Dalton s law of partial pressures! The total pressure exerted by a mixture of gases equals the sum of the partial pressure of the gases! Partial

More information

Study fluid dynamics. Understanding Bernoulli s Equation.

Study fluid dynamics. Understanding Bernoulli s Equation. Chapter Objectives Study fluid dynamics. Understanding Bernoulli s Equation. Chapter Outline 1. Fluid Flow. Bernoulli s Equation 3. Viscosity and Turbulence 1. Fluid Flow An ideal fluid is a fluid that

More information

Numerical Analysis And Experimental Verification of a Fire Resistant Overpack for Nuclear Waste

Numerical Analysis And Experimental Verification of a Fire Resistant Overpack for Nuclear Waste Numerical Analysis And Experimental Verification of a Fire Resistant Overpack for Nuclear Waste Author P. Geraldini 1, Co-author A. Lorenzo 2 1,2 Sogin S.p.A., Rome, Italy Via Marsala 51/C 00185 Rome 1

More information

A Simple Method for Thermal Characterization of Low-Melting Temperature Phase Change Materials (PCMs)

A Simple Method for Thermal Characterization of Low-Melting Temperature Phase Change Materials (PCMs) A Simple Method for hermal Characterization of Low-Melting emperature Phase Change Materials (PCMs) L. Salvador *, J. Hastanin, F. Novello, A. Orléans 3 and F. Sente 3 Centre Spatial de Liège, Belgium,

More information

Time-Dependent Conduction :

Time-Dependent Conduction : Time-Dependent Conduction : The Lumped Capacitance Method Chapter Five Sections 5.1 thru 5.3 Transient Conduction A heat transfer process for which the temperature varies with time, as well as location

More information

Mathematical Modelling of Heating and Evaporation of a Spheroidal Droplet

Mathematical Modelling of Heating and Evaporation of a Spheroidal Droplet ILASS Europe 2016, 27th Annual Conference on Liquid Atomization and Spray Systems, 4-7 September 2016, Brighton, UK Mathematical Modelling of Heating and Evaporation of a Spheroidal Droplet Vladimir Zubkov

More information

convection coefficient, h c = 18.1 W m K and the surrounding temperature to be 20 C.) (20 marks) Question 3 [35 marks]

convection coefficient, h c = 18.1 W m K and the surrounding temperature to be 20 C.) (20 marks) Question 3 [35 marks] COP 311 June Examination 18 June 005 Duration: 3 hours Starting time: 08:30 Internal examiners: Prof. T. Majozi Mnr. D.J. de Kock Mnr. A.T. Tolmay External examiner: Mnr. B. du Plessis Metallurgists: Questions

More information

Problem Solving 3: Calculating the Electric Field of Highly Symmetric Distributions of Charge Using Gauss s Law

Problem Solving 3: Calculating the Electric Field of Highly Symmetric Distributions of Charge Using Gauss s Law MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Problem Solving 3: Calculating the Electric Field of Highly Symmetric Distributions of Charge Using Gauss s Law REFERENCE: Section 4.2, 8.02

More information

Introduction to Heat and Mass Transfer. Week 5

Introduction 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 information

Lumped parameter thermal modelling

Lumped parameter thermal modelling Chapter 3 umped parameter thermal modelling This chapter explains the derivation of a thermal model for a PMSM by combining a lumped parameter (P) model and an analytical distributed model. The whole machine

More information

MAGNETOHYDRODYNAMICS

MAGNETOHYDRODYNAMICS Chapter 6 MAGNETOHYDRODYNAMICS 6.1 Introduction Magnetohydrodynamics is a branch of plasma physics dealing with dc or low frequency effects in fully ionized magnetized plasma. In this chapter we will study

More information

An Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation. Ahsan Munir, PhD Tom Spirka, PhD

An Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation. Ahsan Munir, PhD Tom Spirka, PhD An Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation Ahsan Munir, PhD Tom Spirka, PhD Agenda Provide an overview of COMSOL 4.3b Our products, solutions and applications Subsurface

More information

Specific heat capacity. Convective heat transfer coefficient. Thermal diffusivity. Lc ft, m Characteristic length (r for cylinder or sphere; for slab)

Specific 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 information

Modeling of Humidification in Comsol Multiphysics 4.4

Modeling of Humidification in Comsol Multiphysics 4.4 Modeling of Humidification in Comsol Multiphysics 4.4 Indrajit Wadgaonkar *1 and Suresh Arikapudi 1 1 Tata Motors Ltd. Pimpri, Pune, India, 411018. *Corresponding author: Indrajit Wadgaonkar, Tata Motors

More information

One dimensional steady state diffusion, with and without source. Effective transfer coefficients

One dimensional steady state diffusion, with and without source. Effective transfer coefficients One dimensional steady state diffusion, with and without source. Effective transfer coefficients 2 mars 207 For steady state situations t = 0) and if convection is not present or negligible the transport

More information

Analytical Modeling of Laser Moving Sources

Analytical Modeling of Laser Moving Sources Analytical Modeling of Laser Moving Sources Contains: Heat flow equation Analytic model in one dimensional heat flow Heat source modeling Point heat source Line heat source Plane heat source Surface heat

More information

Unit II Thermal Physics Introduction- Modes of Heat Transfer Normally there are three modes of transfer of heat from one place to another viz., conduction, convection and radiation. Conduction : Conduction

More information

ASSUMPTIONS: (1) One-dimensional, radial conduction, (2) Constant properties.

ASSUMPTIONS: (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 information

Diffusion and Reaction in Fe-Based Catalyst for Fischer- Tropsch Synthesis Using Micro Kinetic Rate Expressions

Diffusion and Reaction in Fe-Based Catalyst for Fischer- Tropsch Synthesis Using Micro Kinetic Rate Expressions Diffusion and Reaction in Fe-Based Catalyst for Fischer- Tropsch Synthesis Using Micro Kinetic Rate Expressions 3-D CFD Model for Shell & Tube Exchanger with 7 Tubes Ender Ozden and Ilker Tari (2010) Multitubular

More information

Aerodynamic Performance 1. Figure 1: Flowfield of a Wind Turbine and Actuator disc. Table 1: Properties of the actuator disk.

Aerodynamic Performance 1. Figure 1: Flowfield of a Wind Turbine and Actuator disc. Table 1: Properties of the actuator disk. Aerodynamic Performance 1 1 Momentum Theory Figure 1: Flowfield of a Wind Turbine and Actuator disc. Table 1: Properties of the actuator disk. 1. The flow is perfect fluid, steady, and incompressible.

More information

Moisture content in concrete

Moisture content in concrete Universal Moisture Transport Model in Concrete under Natural Environment Takumi Shimomura Nagaoka University of Technology Moisture content in concrete Shrinkage and creep of concrete Transport of aggressive

More information

Documentation of the Solutions to the SFPE Heat Transfer Verification Cases

Documentation 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 information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 03 Finned Surfaces Fausto Arpino f.arpino@unicas.it Outline Introduction Straight fin with constant circular cross section Long

More information

Models for measurement of thermophysical parameters

Models for measurement of thermophysical parameters Models for measurement of thermophysical parameters P. Dieška Department of Physics, STU FEI Ilkovičova 3, 812 19 Bratislava, Slovakia October 16, 2008 44 Model as a image of real experiment Heat equation:

More information

Coupled Electromagnetics- Multiphase Porous Media Model for Microwave Combination Heating

Coupled Electromagnetics- Multiphase Porous Media Model for Microwave Combination Heating Presented at the COMSOL Conference 2008 Boston Coupled Electromagnetics- Multiphase Porous Media Model for Microwave Combination Heating Vineet Rakesh and Ashim Datta Biological & Environmental Engineering

More information

Differential equations of mass transfer

Differential equations of mass transfer Differential equations of mass transfer Definition: The differential equations of mass transfer are general equations describing mass transfer in all directions and at all conditions. How is the differential

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree in Mechanical Engineering Numerical Heat and Mass Transfer 02-Transient Conduction Fausto Arpino f.arpino@unicas.it Outline Introduction Conduction ü Heat conduction equation ü Boundary conditions

More information

Taylor Dispersion Created by Robert P. Hesketh, Chemical Engineering, Rowan University Fall 2005

Taylor Dispersion Created by Robert P. Hesketh, Chemical Engineering, Rowan University Fall 2005 Taylor Dispersion Created by Robert P. Hesketh, Chemical Engineering, Rowan University Fall 005 In this problem you will simulate a tubular reactor with fluid flowing in laminar flow. The governing equations

More information

Supplementary Material

Supplementary Material Mangili et al. Supplementary Material 2 A. Evaluation of substrate Young modulus from AFM measurements 3 4 5 6 7 8 Using the experimental correlations between force and deformation from AFM measurements,

More information

Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay

Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Lecture No. 18 Forced Convection-1 Welcome. We now begin our study of forced convection

More information

Solutions to PS 2 Physics 201

Solutions to PS 2 Physics 201 Solutions to PS Physics 1 1. ke dq E = i (1) r = i = i k eλ = i k eλ = i k eλ k e λ xdx () (x x) (x x )dx (x x ) + x dx () (x x ) x ln + x x + x x (4) x + x ln + x (5) x + x To find the field for x, we

More information

MECH 375, Heat Transfer Handout #5: Unsteady Conduction

MECH 375, Heat Transfer Handout #5: Unsteady Conduction 1 MECH 375, Heat Transfer Handout #5: Unsteady Conduction Amir Maleki, Fall 2018 2 T H I S PA P E R P R O P O S E D A C A N C E R T R E AT M E N T T H AT U S E S N A N O PA R T I - C L E S W I T H T U

More information

PROCESS SYSTEMS ENGINEERING Dr.-Ing. Richard Hanke-Rauschenbach

PROCESS SYSTEMS ENGINEERING Dr.-Ing. Richard Hanke-Rauschenbach Otto-von-Guerice University Magdeburg PROCESS SYSTEMS ENGINEERING Dr.-Ing. Richard Hane-Rauschenbach Project wor No. 1, Winter term 2011/2012 Sample Solution Delivery of the project handout: Wednesday,

More information

The effect of Entry Region on Thermal Field

The effect of Entry Region on Thermal Field The effect of Entry Region on Thermal Field Flow Fractionation Nick Cox Supervised by Professor Bruce Finlayson University of Washington Department of Chemical Engineering June 6, 2007 Abstract Achieving

More information

Chapter 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 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 information

Mathematical Modelling and Geometry

Mathematical Modelling and Geometry http://mmg.tversu.ru ISSN 2311-1275 MMG Mathematical Modelling and Geometry Volume 4, No 1, p.34 40 (2016) The mathematical model of wax deposition thickness in a pipeline taking the aging of the deposits

More information

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?

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? 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 information

Table of Contents. Foreword... Introduction...

Table 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 information

Lecture 3a: Surface Energy Balance

Lecture 3a: Surface Energy Balance Lecture 3a: Surface Energy Balance Instructor: Prof. Johnny Luo http://www.sci.ccny.cuny.edu/~luo Surface Energy Balance 1. Factors affecting surface energy balance 2. Surface heat storage 3. Surface

More information

ON THE CONSERVATION OF MASS AND ENERGY IN HYGROTHERMAL NUMERICAL SIMULATION WITH COMSOL MULTIPHYSICS

ON THE CONSERVATION OF MASS AND ENERGY IN HYGROTHERMAL NUMERICAL SIMULATION WITH COMSOL MULTIPHYSICS ON THE CONSERVATION OF MASS AND ENERGY IN HYGROTHERMAL NUMERICAL SIMULATION WITH COMSOL MULTIPHYSICS Michele Bianchi Janetti, Fabian Ochs, and Wolfgang Feist,2 Unit for Energy Efficient Buildings, University

More information

Dynamic Simulation Using COMSOL Multiphysics for Heterogeneous Catalysis at Particle Scale

Dynamic Simulation Using COMSOL Multiphysics for Heterogeneous Catalysis at Particle Scale Dynamic Simulation Using COMSOL Multiphysics for Heterogeneous Catalysis at Particle Scale Ameer Khan Patan *1, Mallaiah Mekala 2, Sunil Kumar Thamida 3 Department of Chemical Engineering, National Institute

More information

Internal Forced Convection. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.

Internal Forced Convection. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Internal Forced Convection Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Introduction Pipe circular cross section. Duct noncircular cross section. Tubes small-diameter

More information

ECE309 INTRODUCTION TO THERMODYNAMICS & HEAT TRANSFER. 10 August 2005

ECE309 INTRODUCTION TO THERMODYNAMICS & HEAT TRANSFER. 10 August 2005 ECE309 INTRODUCTION TO THERMODYNAMICS & HEAT TRANSFER 0 August 2005 Final Examination R. Culham & M. Bahrami This is a 2 - /2 hour, closed-book examination. You are permitted to use one 8.5 in. in. crib

More information

Finite Element Modules for Enhancing Undergraduate Transport Courses: Application to Fuel Cell Fundamentals

Finite Element Modules for Enhancing Undergraduate Transport Courses: Application to Fuel Cell Fundamentals Finite Element Modules for Enhancing Undergraduate Transport Courses: Application to Fuel Cell Fundamentals Originally published in 007 American Society for Engineering Education Conference Proceedings

More information

Application of COMSOL Multiphysics in Nanoscale Electrokinetic Transport

Application of COMSOL Multiphysics in Nanoscale Electrokinetic Transport High Performance Computing Day Application of COMSOL Multiphysics in Nanoscale Electrokinetic Transport Selcuk Atalay, PhD Candidate Advisor: Prof. Shizhi Qian Institute of Micro/Nanotechnology Aerospace

More information

Ion Concentration and Electromechanical Actuation Simulations of Ionic Polymer-Metal Composites

Ion Concentration and Electromechanical Actuation Simulations of Ionic Polymer-Metal Composites October 5-7, 2016, Boston, Massachusetts, USA Ion Concentration and Electromechanical Actuation Simulations of Ionic Polymer-Metal Composites Tyler Stalbaum, Qi Shen, and Kwang J. Kim Active Materials

More information

Optimization of DPF Structures with a 3D-Unit Cell Model

Optimization of DPF Structures with a 3D-Unit Cell Model Optimization of DPF Structures with a 3D-Unit Cell Model Wieland Beckert, Marcel Dannowski, Lisabeth Wagner, Jörg Adler, Lars Mammitzsch Fraunhofer IKTS, Dresden, Germany *Corresponding author: FhG IKTS,

More information

Digital Holographic Measurement of Nanometric Optical Excitation on Soft Matter by Optical Pressure and Photothermal Interactions

Digital Holographic Measurement of Nanometric Optical Excitation on Soft Matter by Optical Pressure and Photothermal Interactions Ph.D. Dissertation Defense September 5, 2012 Digital Holographic Measurement of Nanometric Optical Excitation on Soft Matter by Optical Pressure and Photothermal Interactions David C. Clark Digital Holography

More information

ROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring

ROTATING 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 information

ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer

ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer Solution to Final Examination Date: May 4, 2011 1:00 3:00 PM Instructor: J. Murthy Open Book, Open Notes Total: 100 points Use the finite volume

More information

Internal Flow: Heat Transfer in Pipes

Internal Flow: Heat Transfer in Pipes Internal Flow: Heat Transfer in Pipes V.Vuorinen Aalto University School of Engineering Heat and Mass Transfer Course, Autumn 2016 November 15 th 2016, Otaniemi ville.vuorinen@aalto.fi First about the

More information

CONCURRENT DRYING OF SOYBEAN SEEDS: THE EFFECT OF THE RADIAL AIR PROFILE

CONCURRENT DRYING OF SOYBEAN SEEDS: THE EFFECT OF THE RADIAL AIR PROFILE Brazilian Journal of Chemical Engineering ISSN 0104-6632 Printed in Brazil www.abeq.org.br/bjche Vol. 22, No. 02, pp. 311-318, April - June, 2005 CONCURRENT DRYING OF SOYBEAN SEEDS: THE EFFECT OF THE RADIAL

More information

Turbulence Model Affect on Heat Exchange Characteristics Through the Beam Window for European Spallation Source

Turbulence Model Affect on Heat Exchange Characteristics Through the Beam Window for European Spallation Source International Scientific Colloquium Modelling for Material Processing Riga, September 16-17, 2010 Turbulence Model Affect on Heat Exchange Characteristics Through the Beam Window for European Spallation

More information

FORMULA SHEET. General formulas:

FORMULA SHEET. General formulas: FORMULA SHEET You may use this formula sheet during the Advanced Transport Phenomena course and it should contain all formulas you need during this course. Note that the weeks are numbered from 1.1 to

More information

2, where dp is the constant, R is the radius of

2, where dp is the constant, R is the radius of Dynamics of Viscous Flows (Lectures 8 to ) Q. Choose the correct answer (i) The average velocity of a one-dimensional incompressible fully developed viscous flow between two fixed parallel plates is m/s.

More information

THERMOCAPILLARY CONVECTION IN A LIQUID BRIDGE SUBJECTED TO INTERFACIAL COOLING

THERMOCAPILLARY CONVECTION IN A LIQUID BRIDGE SUBJECTED TO INTERFACIAL COOLING THERMOCAPILLARY CONVECTION IN A LIQUID BRIDGE SUBJECTED TO INTERFACIAL COOLING Melnikov D. E. and Shevtsova V. M. Abstract Influence of heat loss through interface on a supercritical three-dimensional

More information

Generating cloud drops from CCN. Wallace & Hobbs (1977)

Generating cloud drops from CCN. Wallace & Hobbs (1977) Generating cloud drops from CCN Wallace & Hobbs (1977) Cloud Drops and Equilibrium Considera3ons: review We discussed how to compute the equilibrium vapor pressure over a pure water drop, or a solu3on

More information

FIND: (a) Sketch temperature distribution, T(x,t), (b) Sketch the heat flux at the outer surface, q L,t as a function of time.

FIND: (a) Sketch temperature distribution, T(x,t), (b) Sketch the heat flux at the outer surface, q L,t as a function of time. PROBLEM 5.1 NOWN: Electrical heater attached to backside of plate while front surface is exposed to convection process (T,h); initially plate is at a uniform temperature of the ambient air and suddenly

More information

Advantages, Limitations and Error Estimation of Mixed Solid Axisymmetric Modeling

Advantages, Limitations and Error Estimation of Mixed Solid Axisymmetric Modeling Advantages, Limitations and Error Estimation of Mixed Solid Axisymmetric Modeling Sudeep Bosu TATA Consultancy Services Baskaran Sundaram TATA Consultancy Services, 185 Lloyds Road,Chennai-600086,INDIA

More information

Vortices in Superfluid MODD-Problems

Vortices in Superfluid MODD-Problems Vortices in Superfluid MODD-Problems May 5, 2017 A. Steady filament (0.75) Consider a cylindrical beaker (radius R 0 a) of superfluid helium and a straight vertical vortex filament in its center Fig. 2.

More information

Introduction to Heat and Mass Transfer. Week 7

Introduction to Heat and Mass Transfer. Week 7 Introduction to Heat and Mass Transfer Week 7 Example Solution Technique Using either finite difference method or finite volume method, we end up with a set of simultaneous algebraic equations in terms

More information

6.9 Applications of SolidLiquid Phase Change

6.9 Applications of SolidLiquid Phase Change 6.9 Applications of Solid-Liquid Phase Change 6.9.1 Latent Heat Thermal Energy Storage The theoretical model employed in this study is shown in Fig. 6.31. At the very beginning of the process (t =0), the

More information

PROBLEM 14.6 ( )( ) (b) Applying a species balance to a control volume about the hydrogen, dt 6 dt 6RAT dt 6RT dt

PROBLEM 14.6 ( )( ) (b) Applying a species balance to a control volume about the hydrogen, dt 6 dt 6RAT dt 6RT dt PROBLEM 14.6 KNOWN: Pressure and temperature of hydrogen stored in a spherical steel tank of prescribed diameter and thickness. FIND: (a) Initial rate of hydrogen mass loss from the tank, (b) Initial rate

More information

1. Thermal energy is transferred through the glass windows of a house mainly by. D. radiation and convection. (1)

1. Thermal energy is transferred through the glass windows of a house mainly by. D. radiation and convection. (1) 1. Thermal energy is transferred through the glass windows of a house mainly by A. conduction. B. radiation. C. conduction and convection. D. radiation and convection. 2. The specific latent heat of vaporization

More information

Chapter 17. Fundamentals of Atmospheric Modeling

Chapter 17. Fundamentals of Atmospheric Modeling Overhead Slides for Chapter 17 of Fundamentals of Atmospheric Modeling by Mark Z. Jacobson Department of Civil & Environmental Engineering Stanford University Stanford, CA 94305-4020 August 21, 1998 Mass

More information

Errata for SOLIDIFICATION (Second Edition, 2016)

Errata for SOLIDIFICATION (Second Edition, 2016) Errata for SOLIDIFICATION (Second Edition, 2016) J. A. Dantzig and M. Rappaz September 6, 2017 Nomenclature There are several minor typographical errors in this table. Please download the corrected version

More information

Chapter 24. Gauss s Law

Chapter 24. Gauss s Law Chapter 24 Gauss s Law Let s return to the field lines and consider the flux through a surface. The number of lines per unit area is proportional to the magnitude of the electric field. This means that

More information

Excerpt from the Proceedings of the COMSOL Users Conference 2007 Grenoble

Excerpt from the Proceedings of the COMSOL Users Conference 2007 Grenoble Excerpt from the Proceedings of the COSOL Users Conference 007 Grenoble Evaluation of the moderator temperature coefficient of reactivity in a PWR V. emoli *,, A. Cammi Politecnico di ilano, Department

More information

Thermal Diffusivity of Plastic

Thermal Diffusivity of Plastic Thermal Diffusivity of Plastic School of Physics and Astronomy University of Manchester Amir El-hamdy Sohail Mahmood November 19, 2015 Date Performed: October 6-13, 2015 Demonstrator: Henry Cox Abstract

More information

This section develops numerically and analytically the geometric optimisation of

This section develops numerically and analytically the geometric optimisation of 7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically

More information

Numerical Study of a High Temperature Latent Heat Storage ( C) Using NaNO 3 -KNO 3 Binary Mixture

Numerical 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 information