The Numerical Study of Transport Phenomena in Porous Media

Total Page:16

File Type:pdf, Size:1020Kb

The Numerical Study of Transport Phenomena in Porous Media A Numerical Study of Transport Phenomena in Porous Media by May-Fun Liou Submitted in partial fulfillment of the requirements For the degree of Doctor of Philosophy Thesis Adviser: Dr. Isaac Greber Department of Mechanical and Aerospace Engineering Case Western Reserve University August, 2005 CASE WESTERN RESERVE UNIVERSITY SCHOOL OF GRADUATE STUDIES We hereby approve the dissertation of ______________________________________________________ candidate for the Ph.D. degree *. (signed)_______________________________________________ (chair of the committee) ________________________________________________ ________________________________________________ ________________________________________________ ________________________________________________ ________________________________________________ (date) _______________________ *We also certify that written approval has been obtained for any proprietary material contained therein. For my beloved parents, Bow-Churng & Chung Wan ii Table of Contents List of Tables vii List of Figures viii Nomenclature xvi Abstract xx 1 Introduction 1 1.1 General Overview 2 1.2 Volume Averaging Method 8 1.3 Three Dimensional Navier-Stokes Equations 13 1.3.1 Fluid-phase Equations 13 1.3.2 Solid-phase Equation 15 1.4 Material under Studied 15 1.5 Outline of the Thesis 16 2 An Enabling Technique for Pore-scale Modeling 18 2.1 Introduction 18 2.2 Mesh-based Microstructure Representing Algorithm 22 2.2.1 Mesh Generator 28 2.2.2 Random Number Generator 32 2.2.3 Quality Control 34 2.2.4 Porosity 37 2.3 Formation of the Two - dimensional Geometry 38 iii 2.4 Formation of the Three - dimensional Geometry 40 2.5 Numerical Solutions 43 2.5.1 Governing equations - both fluid and solid phases 43 2.5.2 Boundary conditions 44 2.5.3 Numerical methods 46 2.6 Sensitivity Test 47 2.6.1 Geometric types 47 2.6.2 Grid dependence study 47 2.6.3 Convergence study in terms of residual 49 2.7 Discussion 50 2.8 Summary 51 3 Fluid Flow and Pressure Drop in Porous Media 53 3.1 Introduction 53 3.2 Validation Cases 56 3.2.1 Laminar flow over a square cylinder 56 3.2.2 Laminar flow over a circular cylinder 57 3.3 Numerical Experiments of 2D Porous Flow 60 3.4 Results and Discussion 61 3.4.1 Two-dimensional channel flow 62 3.4.1.1 Velocity field 63 3.4.1.2 Pressure drop and porosity 70 3.5 Summary 73 iv 4 Forced Convection Heat Transfer in Porous Media 74 4.1 Introduction 74 4.2 Validation Cases – forced convection in single plate-plate channel 84 4.2.1 Bottom plate as an inner plate 86 4.2.2 Bottom plate as an external plate 90 4.3 Numerical Setup for Forced Convective Heat Transfer 94 4.3.1 Two-dimensional channel flow 95 4.3.1.1 Results and discussion 97 4.3.1.2 Effective thermal conductivity 117 4.3.1.3 Thermal dispersion 122 4.3.2 Summary 137 4.4 Three-dimensional Rectangular Duct Flow 138 4.4.1 Flow with Heated Porous Section 138 4.4.2 Summary 150 4.5 Conclusion 150 5 Closure 5.1 Summary 152 5.2 Future Work 154 Appendix A Programs of Random number generating 156 v Appendix B Governing equations and numerical methods 160 Governing equations 160 Numerical method 162 Appendix C Fluid flow and pressure drop data 168 References 169 vi List of Tables Table 1 : List of thermophysical properties of solid materials under study……… 16 Table 2 : Results from grid-dependence study………………………………….. 48 Table 3 : Summary of previous works in evaluating viscous resistance (pressure drop) in porous media……………………………………………………….. 54 Table 4 : List of pressure drops for various porous systems…………………….. 61 Table 5 : Summary of numerical simulations of forced convective heat transfer in Porous media………………………………………………………….. 80 Table 6 : The effective heat conductivity of fluid phase at Re = 270…………… 120 Table 7 : The ratio of effective heat conductivity of fluid phase, f, at Re = 1,100. 121 Table 8 : The ratio of effective heat conductivity of fluid phase, f, at Re = 2,200.. 122 Table C-1: The pressure drop of a simplified electronic cooling system at various conditions………………………………………………………………. 168 vii List of Figures 1.1 Schematic of porous medium, flow variations in the representative Elementary Volume and the pore velocity vector Vp....................................................... 8 2.1 A Schematic explaining the procedures of generating a numerical porous sample………………………..................................................................... 28 2.2 A porous sample with non-uniform porosity……………………………… 34 2.3 Schematic of typical domain of a porous medium – two phase problem…. 37 2.4 Randomly formed solid matrix representing porous medium in 2-D domain (a) fibrous shape formed with triangular mesh, (b) small solids in circular shape formed with triangular mesh, (c) large solids in circular shape with hybrid mesh……………………… 39 2.5 Open-cell foam microstructure …………………………………………… 41 2.6 (a) the cross-section of a typical mesh-based porous medium is shown with solid balls in black and space taken by fluid in white……………. 41 2.6 (b) The cross-section of numerical generated open-cellular sample after converting all solid and fluid tags from Figure 2.8(a), into the oppose ones……………………………………………………………………... 42 2.6 (c) 3D Numerical fluid and solid mesh in cross cut planes along the axial direction – x direction. …………………..……………………………… 42 viii 2.7 Convergence history of residuals from solving Navier-Stokes equations for fluid and solving heat conduction equations for solid in the case of the channel flow heated from both top and bottom walls. …………………...………… 50 3.1 Streamlines of a low Reynolds number flow exhibits symmetry around the solid block ………………………………………………………………… 57 3.2 A hybrid mesh is used for studying unsteady laminar flow over a stationary circular cylinder. ………...………………………………………………… 58 3.3 Contours of u velocity indicate that the laminar flow is fully developed in the flow field …..……………………………………………………………… 59 3.4 Schematic numerical setup for measuring pressure drop …………………... 60 3.5(a) The velocity distribution along the distance from the bottom wall, y', with respect to sampling locations of ∆1, ∆2, ∆3, ∆4, ∆5 and at the exit of the channel …...………………………………………………………………. 64 3.5(b) The velocity distributions in the axial (X) and vertical direction (Y) respectively ..………… ..………………… ……………………………… 65 3.6(a) Velocity contours of the whole numerical domain for the inflow velocity 0.1m/s ………………………..……………………………………….... 66 3.6(b) Vortices contours of the whole numerical domain for the inflow velocity 0.1m/s …………………………………………………………………… 67 3.7 (a) The pressure distribution along the non-dimensional axial direction with respect to several y locations close to the wall……………………. …… 68 3.7 (b) The pressure distribution along the non-dimensional axial direction at three y locations, farther away from the wall………………………... 69 ix ∗ 3.8 Dimensionless pressure gradient parameter, P vs. Red, of various porous systems in Table 3 ………………………………………………………….. 70 ∗ 3.9 Non-dimensional pressure drop P vs. Red in a wide range of Reynolds numbers for two values of porosities.……………………………………… 71 ∗ 3.10 A close-up view of the non-dimensional pressure drop P vs. Red, exhibiting different behaviors caused by varying porosityε …………………………. 72 4.1 Schematic cooling of electronic board with a substrate between coolant and board. ………………………………………………………………………. 85 4.2 Silicon rubber substrate with isothermal walls at 900 K, (a) velocity contours, (b) temperature contours…………………………………….……………… 87 4.3 Local Nusselt number Nux versus dimensionless x of heat transfer at the contact surface of coolant and substrate of silicon rubber………………………. 88 4.4 Aluminum substrate with constant wall temperature, Tw=900K, (a) velocity contours, (b) temperature contours……………………………………… 89 4.5 Silicon rubber substrate with Tw = 900 K and the other parallel plate insulated (a) velocity contours, (b) temperature contours………………………… 90 4.6 Local Nusselt number Nux versus dimensionless x of heat transfer at the contact surface of coolant and substrate of aluminum………………… 92 4.7 Aluminum substrate with isothermal top wall at Tw = 900 K, and the lower parallel plate insulated: (a) velocity contours, (b) temperature contours. 93 4.8 Nux along longitudinal direction x at the contact surface of coolant and substrates of silicon rubber and aluminum…………………………….. 94 x 4.9 Schematic of a two-dimensional porous channel with heat source from below or above…………………………………………………………. 95 4.10(a) Temperature contours of the nonporous channel with the background showing the mesh used………………………………………………… 99 4.10(b) The invariant temperature profiles along the transverse direction in the channel at various sampling positions along the axial direction………. 99 4.10(c) Velocity profiles along the transverse direction at all sampling positions.100 4.10(d) Local Nusselt number Nux versus the longitudinal coordinate x of the baseline case…………………………………………………………… 100 4.11(a) Temperature contours of the porous sample in a channel, along with the mesh used and the porous region in which the solid matrices are shown in black…………………………………………………………………… 101 4.11(b) The invariant temperature profiles along the direction of y* in the channel are plotted for various sampling positions……………………………. 102 4.11(c) Profiles of non-dimensional x-component velocity along the transverse direction at various sampling locations………………………………… 104 4.11(d) Plot of local Nusselt number Nux versus the longitudinal coordinate x… 104 4.12 The ratio, f, of local heat conductivity of fluid phase to that at the inflow condition (300K and 1atm) in terms of the distance away from wall at various sampling locations……………………………………………… 105 4 .13 Distributions of local dimensionless heat transfer along the bottom wall of four solid phases, the upward spikes indicate the enhancement of heat transfer contributed from the physical blockages of porous media……. 108 xi 4.14 Vorticity distribution along the dimensionless transverse direction, y, on the mid-plane of the porous medium (∆3)…………………………………. 109 4.15 Vorticity distributions on the center plane of the hypothetical solid/air (system of θ = 1.0) under Case (b)……………………………………... 111 4.16 Local dimensionless heat transfer rate along the bottom plate at various Reynolds numbers.
Recommended publications
  • Guggenheim Aeronautical Laboratory~://; California
    R.ES~.LttCH REPCRT GUGGENHEIM AERONAUTICAL LABORATORY ~://; CALIFORNIA INSTITUTE OF TECHNOLOGY HYPERSC»JIC RESEARCH PROJECT Memorandum No. 47 December 15, 1958 AN EXPERIMENTAL INVESTIGATION OF THE EFFECT OF EJECTING A COOLANT GAS AT THE NOSE OF A BLUNT BODY by C. Hugh E. Warren ARMY ORDNANCE CONTRACT NO. DA-04-495-0rd-19 GUGGENHEIM AERONAUTICAL LABORATORY CALIFORNIA INSTITUTE OF TECHNOLOGY Pasadena1 California HYPERSONIC RESEARCH PROJECT Memorandum No. 4 7 December 151 1958 AN EXPERIMENTAL INVESTIGATION OF THE EFFECT OF EJECTING A COOLANT GAS AT THE NOSE OF A BLUNT BODY by C. Hugh E. Warren ~~C!a r~:Ml ia:nl bil'eCtOr Guggenheim Aeronautical Laboratory ARMY ORDNANCE CONTRACT NO. DA-04-495-0rd-19 Army Project No. 5B0306004 Ordnance Project No. TB3-0118 OOR Project No. 1600-PE ACKNOWLEDGMENTS The author wishes to express his appreciation to Professor Lester Lees for his assistance1 guidance and encouragement during this work1 to Mr. Howard McDonald and other members of the Aeronautics Machine Shop for fabricating and maintaining the models and other items of equipment1 to Mr. Paul Baloga and the staff of the Hypersonic Wind Tunnel for their assistance and advice during testing~ to Mrs. Betty Laue for much computational work1 and to Mrs. Gerry Van Gieson for typing the manuscript and piloting the report through after the author 1 s departure. Most of the work was done while the author was in receipt of a Harkness Fellowship of the Commonwealth Fund of New York1 to whom thanks are due for the opportunity to spend a year on such an instructive and interesting program in such congenial surroundings.
    [Show full text]
  • Application of Laser Doppler Velocimetry
    APPLICATION OF LASER DOPPLER VELOCIMETRY (LDV) TO STUDY THE STRUCTURE OF GRAVITY CURRENTS UNDER FIRE CONDITIONS B Moghtaderi† Process Safety and Environment Protection Group Discipline of Chemical Engineering, School of Engineering Faculty of Engineering & Built Environment, The University of Newcastle Callaghan, NSW 2308, Australia ABSTRACT Gravity currents are of considerable safety importance primarily because of their role in the spread and transport of smoke and hot gases in building fires. Despite recent progress in the field, relatively little is known about the structure of gravity currents under conditions pertinent to building fires. The present investigation is an attempt to address this shortcoming by studying the turbulent structure of gravity currents. For this purpose, a series of experiments was conducted in a rectangular tank with turbulent, sub-critical underflows. Laser-Doppler Velocimetry was employed to quantify the velocity field and associated turbulent flow parameters. Experimental results indicated that the mean flow within the head region primarily consisted of an undiluted large single vortex which rapidly mixed with the ambient flow in the wake region. Cases with isothermal wall boundary conditions showed three-dimensional effects whereas those with adiabatic walls exhibited two-dimensional behaviour. Turbulence was found to be highly heterogeneous and its distribution was governed by the location of large eddies. While all components of turbulence kinetic energy showed minima in the regions where velocity was maximum (i.e. low fluid shear), they reached their maximum in the shear layer at the upper boundary of the flow. KEY WORDS: gravity currents, laser Doppler velocimetry, turbulent structure INTRODUCTION Gravity currents are important physical phenomena which have direct implications in a wide variety of physical situations ranging from environmental phenomena‡ to enclosure fires.
    [Show full text]
  • Conduction Heat Transfer Notes for MECH 7210
    Conduction Heat Transfer Notes for MECH 7210 Daniel W. Mackowski Mechanical Engineering Department Auburn University 2 Preface The Notes on Conduction Heat Transfer are, as the name suggests, a compilation of lecture notes put together over 10 years of teaching the subject. The notes are not meant to be a comprehensive ∼ presentation of the subject of heat conduction, and the student is referred to the texts referenced below for such treatments. A goal of mine, in preparing the notes, has been to address an apparent shortcoming in many of the current texts, in that the texts present the mathematical formulation and analytical solution to a wide variety of conduction problems, yet they spend little if any time on discussing how numerical and graphical results can be obtained from the solutions. As will be seen, this task in itself is not trivial, and to this end mathematical software packages (in particular, the package Mathematica) will be used extensively in application of the analytical solutions. The notes were prepared using the LATEX typesetting program, which is freely available via internet download. I wish to thank my former students, who have (and continue) to catch the multitude of mistakes and typos in the notes. These notes are dedicated to the memory of Clifford Cremers, an outstanding teacher of heat transfer and a fine fly fisherman. References The ‘text’ in the course will consist of my lecture notes – which contain few if any literature citations. I will need to fix this if I ever expect to publish the notes as a book. The following reference texts were used to prepare the notes.
    [Show full text]
  • Thermodynamic Processes: the Limits of Possible
    Thermodynamic Processes: The Limits of Possible Thermodynamics put severe restrictions on what processes involving a change of the thermodynamic state of a system (U,V,N,…) are possible. In a quasi-static process system moves from one equilibrium state to another via a series of other equilibrium states . All quasi-static processes fall into two main classes: reversible and irreversible processes . Processes with decreasing total entropy of a thermodynamic system and its environment are prohibited by Postulate II Notes Graphic representation of a single thermodynamic system Phase space of extensive coordinates The fundamental relation S(1) =S(U (1) , X (1) ) of a thermodynamic system defines a hypersurface in the coordinate space S(1) S(1) U(1) U(1) X(1) X(1) S(1) – entropy of system 1 (1) (1) (1) (1) (1) U – energy of system 1 X = V , N 1 , …N m – coordinates of system 1 Notes Graphic representation of a composite thermodynamic system Phase space of extensive coordinates The fundamental relation of a composite thermodynamic system S = S (1) (U (1 ), X (1) ) + S (2) (U-U(1) ,X (2) ) (system 1 and system 2). defines a hyper-surface in the coordinate space of the composite system S(1+2) S(1+2) U (1,2) X = V, N 1, …N m – coordinates U of subsystems (1 and 2) X(1,2) (1,2) S – entropy of a composite system X U – energy of a composite system Notes Irreversible and reversible processes If we change constraints on some of the composite system coordinates, (e.g.
    [Show full text]
  • Heat Transfer in Nanofluid Boundary Layer Near Adiabatic Wall
    Copyright © 2018 by American Scientific Publishers Journal of Nanofluids All rights reserved. Vol. 7, pp. 1297–1302, 2018 Printed in the United States of America (www.aspbs.com/jon) Heat Transfer in Nanofluid Boundary Layer Near Adiabatic Wall D. Hopper, D. Jaganathan, J. L. Orr, J. Shi, F. Simeski, M. Yin, and J. T. C. Liu∗ School of Engineering (ENGN 2760) and Fluids@Brown, Brown University, Providence, Rhode Island 02912, USA Heat transfer in nanofluid boundary layers has been receiving significant attention in recent years due to widespread applications and ongoing research in both academia and industry. Cooling systems in batteries, semiconductors, and other components in the hardware technology industry are just a few areas where the limits of heat transfer are being pushed. The principal region of heat exchange between the wall and fluid is within the thermal boundary layer. Intriguingly, the heat transfer of particle-laden nanofluid is greatly affected by both the appearance of the particle and the nanoscale molecular effects; hence, a thorough characterization of the boundary layer behaviour and properties is needed. Prior research has shown that nanofluids enhance the heat transfer. Based on related works (J. T. C. Liu, Proc. Royal Society A 408, 2383 (2012); J. T. C. Liu, et al., Arch. Mech. 69, 75 (2017); C. J. B. de Castilho, et al., J. Heat Transfer Eng. 1 (2018)), convective heat transfer in nanofluid boundary layers is studied by a theoretical model with first-order perturbation of variables. Small nanoparticle volume concentration is assumed. A theoretical model of nanofluid flow over an adiabatic wall is ARTICLE presented, and the corresponding temperature and velocity profiles are obtained.
    [Show full text]
  • Thermodynamics
    CHAPTER TWELVE THERMODYNAMICS 12.1 INTRODUCTION In previous chapter we have studied thermal properties of matter. In this chapter we shall study laws that govern thermal energy. We shall study the processes where work is 12.1 Introduction converted into heat and vice versa. In winter, when we rub 12.2 Thermal equilibrium our palms together, we feel warmer; here work done in rubbing 12.3 Zeroth law of produces the ‘heat’. Conversely, in a steam engine, the ‘heat’ Thermodynamics of the steam is used to do useful work in moving the pistons, 12.4 Heat, internal energy and which in turn rotate the wheels of the train. work In physics, we need to define the notions of heat, 12.5 First law of temperature, work, etc. more carefully. Historically, it took a thermodynamics long time to arrive at the proper concept of ‘heat’. Before the 12.6 Specific heat capacity modern picture, heat was regarded as a fine invisible fluid 12.7 Thermodynamic state filling in the pores of a substance. On contact between a hot variables and equation of body and a cold body, the fluid (called caloric) flowed from state the colder to the hotter body! This is similar to what happens 12.8 Thermodynamic processes when a horizontal pipe connects two tanks containing water 12.9 Heat engines up to different heights. The flow continues until the levels of 12.10 Refrigerators and heat water in the two tanks are the same. Likewise, in the ‘caloric’ pumps picture of heat, heat flows until the ‘caloric levels’ (i.e., the 12.11 Second law of temperatures) equalise.
    [Show full text]
  • MODI Fled CROCCO- LEES M Lxi NG THEORY for SUPERSONIC SEPARATED and REATTACHING FLOWS by Herbert S
    RESEARCH REPORT ·~yl GUGGENHEIM AERONAUTICAL LABORATORY CALIFORNIA INSTITUTE OF TECHNOLOGY ~ HYPERSONIC RESEARCH PROJECT Memorandum No. 53 May 2, 1960 MODI FlED CROCCO- LEES M lXI NG THEORY FOR SUPERSONIC SEPARATED AND REATTACHING FLOWS by Herbert S . Glick AUG - 3 1960 ARMY ORDNANCE CONTRACT NO. DA-04-495-0rd-1960 - GUGGENHEIM AERONAUTICAL LABORATORY CALIFORNIA INSTITUTE OF TECHNOLOGY Pasadena, California HYPERSONIC RESEARCH PROJECT ERRATA FOR 1962 Memorandwn No. 53 .f~o ~ May 2, 1960 ~ £NGNG. SCI. t\Y.~~ "Modified Crocco-Lees Mixing Theory for Supersonic Separated and Reattaching Flows" by Herbert S. Glick page 74 -- line 12 -- (nwnerator) in the equation L = (denom1nator) the numerator of the equation should read: GUGGENHEIM AERONAUTICAL LABORATORY CALIFORNIA INSTITUTE OF TECHNOLOGY Pa eadena, California HYPERSONIC RESEARCH PROJECT Memorandum No. 53 May Z, 1960 MODIFIED CROCCO-LEES MIXING THEORY FOR SUPERSONIC SEPARATED AND REATTACHING FLOWS by Herbert S. Glick C~ ar 1 • n uector ARMY ORDNANCE CONTRACT NO. DA-04-495-0rd-1960 ACKNOWLEDGMENTS The author wishes to express his appreciation to Professor Lester Lees for his guidance throughout the course of the investigation. He also wishes to thank Mrs. Betty Laue for her enthusiastic and able assistance in carrying out most of the desk computations; Mrs. Truus van Harreveld, Miss Georgette A. Pauwels, and Mr. Andy Chapkis who assisted in the early computations; Mrs. Betty Wood for preparing the figures, and Mrs. Geraldine VanGieson for her typing of the manu­ script. The author acknowledges with gratitude the receipt of a sabbatical leave grant from the Cornell Aeronautical Laboratory, Inc. of Buffalo, New York for the year 1957-1958, and a fellowship from the Curtiss­ Wright Corporation for the year 1958-1959.
    [Show full text]
  • The Foundations of Classical Thermodynamics
    University of Wollongong Research Online Wollongong University College Bulletin Corporate Publications Archive 11-1959 The Foundations of Classical Thermodynamics L C. Woods Follow this and additional works at: https://ro.uow.edu.au/wucbull Recommended Citation Woods, L C., "The Foundations of Classical Thermodynamics" (1959). Wollongong University College Bulletin. 4. https://ro.uow.edu.au/wucbull/4 Research Online is the open access institutional repository for the University of Wollongong. For further information contact the UOW Library: [email protected] The Foundations of Classical Thermodynamics This serial is available at Research Online: https://ro.uow.edu.au/wucbull/4 BULLETIN NO. 2 THE UNIVERSITY OF NEW SOUTH WALES \ The Foundations of Classical Thermodynamics L. c. W O O D S School of Mechanical Engineering November, 1959 THF FOUNDATIONS of CLASSICAL THERMODYNAMICS by L. C. W O O D S (Nuffield Research Professor of Mechanical Engineering the University of New South Wales, Australia.) A course of lectures given in 1957. TKF FOUNDATION? of CL A ? ? T C AL Tr FRVO PYM AV- ICC 1, Temperature 1,1 Thermodynamic f-ystems As thermodynamics is a branch of physics, it is concerned with certain characteristics of a particular domain of space and matter, a domain of interest which we call a system,. Fverything outside a system which could have a direct influence on the fceha,viour of the system is termed the surroundings, A system can be defined by stating its various physical and chemical properties in fully but for many purposes
    [Show full text]
  • Section 1 Introduction to Classical Thermodynamics
    Section 1 Introduction to Classical Thermodynamics 1.1 Introduction – Macroscopic and microscopic descriptions Any system may be described macroscopically or microscopically (G pp. 1-4, F pp. xi-xii). The macroscopic description considers the system as a whole. The system is characterised, in equilibrium, by a relatively small number of variables e.g. pressure, volume, temperature, internal energy, entropy. The state of a system described in this way is referred to as a macrostate. The macrostate of a system can be changed in two principal ways: we can do work on the system and/or heat can flow in. (The macrostate would also change if particles flowed into or out of the system.) The general rules governing these macroscopic variables and their inter-relationships is the subject of thermodynamics. The microscopic description considers the detailed microscopic nature of the system. At this level the system is characterised by a large number of variable which specify the state of the microscopic entities that make it up. For example, classically the state of a gas is given at any instant by the position and momentum of each particle. So for N particles we need 6N variables. This is a very large number; for 1cm3 of gas at STP, N is of order 1019. The state of a system described in this way is referred to as a microstate. The positions and momenta of the particles, and thus the microstate, evolves with time according to the laws of mechanics. Practically, we could never hope to perform such a calculation! In the quantum mechanical description a microstate will correspond to the system having a given wave function.
    [Show full text]
  • The Portion of the Universe Under Study Surroundings
    Thermodynamics: a few definitions system: the portion of the universe under study surroundings: the rest of the universe boundary: separates the system from the surroundings closed system: matter cannot enter or leave (impermeable boundary) heat may enter or leave; work may be done open system: matter may enter or leave (permeable boundary) isolated system: boundary prevents any interaction with the system no heat or matter may enter or leave no work may be done by the system on the surroundings or by the surroundings on the system rigid, impermeable, adiabatic wall adiabatic system: no heat exchange between the system and the surroundings the boundary is an adiabatic wall an adiabatic system is a perfectly insulated system state: the state of a system is defined by specifying its properties properties with numerical values can be divided into two categories: intensive properties: magnitude does not depend upon the amount of substance temperature, pressure, concentration, density, molar volume extensive properties: magnitude depends on the amount of substance total value equals the sum of its values for all parts of the system volume, mass, change in enthalpy homogeneous system: intensive properties constant throughout (opposite: heterogeneous) a phase is a homogeneous part of a system A state variable (also called a state function, thermodynamic variable, thermodynamic function) depends only on the state of the system. state of equilibrium: (a) system shows no further tendency to change its properties with time (b) removal from contact with surroundings causes no further change in properties if (a) but not (b), system is in a steady state) can specify an equilibrium state by its state functions; in fact, only a few are needed (not as easy to specify a nonequilibrium state) one can further subdivide the concept of equilibrium: 1.
    [Show full text]
  • Investigation of Thermal Adiabatic Boundary Condition on Semitransparent Wall in Combined Radiation and Natural Convection
    Investigation of Thermal Adiabatic Boundary Condition on Semitransparent Wall in Combined Radiation and Natural Convection G Chanakya, Pradeep Kumar Numerical Experiment Laboratory (Radiation and Fluid Flow Physics) School of Engineering Indian Institute of Technology Mandi Mandi, Himachal Pradesh, India 175075 Abstract Two thermal adiabatic boundary conditions arise on the semitransparent win- dow owing to the fact that whether semitransparent window allows the energy to leave the system by radiation mode of heat transfer. It is assumed that being low conductivity of semitransparent material, energy does not leave by conduction mode of heat transfer. This does mean that the semitransparent window may behave as only conductively adiabatic (qc = 0) or combinedly con- ductively and radiatively adiabatic (qc +qr = 0). In the present work, the above two thermal adiabatic boundary conditions have been investigated in natural convection problem for the Rayleigh number (Ra) 105 and Prandtl number(Pr) 0.71 in a cavity, whose left vertical wall has been divided into upper and lower parts in the ratio of 4:6. The upper section is semitransparent window, while lower section is isothermal wall at a temperature of 296K. A collimated beam is irradiated with different values (0, 100, 500 and 1000 W=m2) on the semitrans- parent window at an angle of 450. The cavity is heated from the bottom by arXiv:2007.12484v1 [physics.flu-dyn] 23 Jul 2020 convective heating with free stream temperature of 305K and heat transfer coef- Email address: [email protected] (Pradeep Kumar) Collimated beam on Semitransparent Wall G chanakya, Pradeep Kumar ficient of 50 W=m2K while right wall is also isothermal at same temperature as of lower left wall and upper wall is adiabatic.
    [Show full text]
  • Hypersonic Laminar Boundary Layers with Adiabatic Wall Condition
    HYPERSONIC LAMINAR BOUNDARY LAYERS WITH ADIABATIC WALL CONDITION Paulo G. de P. Toro ([email protected]) Instituto de Aeronaútica e Espaço - IAE - Centro Técnico Aeroespacial - CTA São José dos Campos - SP 12228-904 - BRAZIL Zvi Rusak ([email protected]) Leik N. Myrabo ([email protected])and Henry T. Nagamatsu, ([email protected]) Department of Mechanical Engineering, Aeronautical Engineering, and Mechanics Rensselaer Polytechnic Institute, Troy, NY 12180-3590 - USA Abstract. A new set of self-similar solutions of a compressible laminar boundary layer is used for air as perfect gas and where the viscosity is a power function of temperature. Modified Levy-Mangler and Dorodnitsyn-Howarth transformations are presented to solve the flow in a thin laminar boundary layer with no external pressure gradients on a smooth flat plate. This result in an explicit relation between the stream function and the enthalpy fields described by a closed-coupled system of nonlinear ordinary differential equations. In the present work, boundary layer flows with external Mach numbers 4 and 10 over in adiabatic wall are studied. The present solution methodology provides a straightforward way of comparing results using the viscosity-temperature linear relation, Sutherland law, and the relation according to the kinetic theory. Also, the results may provide important data needed for the design of future hypersonic vehicles. Keywords: Self-similar equations, Compressible laminar boundary layer, Supersonic and hypersonic flows. 1. INTRODUCTION Supersonic and hypersonic flows with low enthalpy conditions can be considered and may be modeled by a calorically or thermally perfect gas equation of state. In a calorically perfect gas the specific heats, cp and cv, are considered as constant.
    [Show full text]