DOCSLIB.ORG

Natural Convection in Porous Media: a Numerical Study of Brinkman Model

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Natural Convection in Porous Media: a Numerical Study of Brinkman Model Transactions on Modelling and Simulation vol 24, © 1999 WIT Press, www.witpress.com, ISSN 1743-355X Natural convection in porous media: a numerical study of Brinkman model R. Jecl?>, L. Skergef\ E. Petresin^ ^ Faculty of Civil Engineering, University ofMaribor, 2000 Maribor, Slovenia EMail: [email protected] Faculty of Mechanical Engineering, University ofMaribor, 2000 Maribor, Slovenia EMail: [email protected] Abstract The problem of natural convection in porous media is being investigated utilising a Boundary Domain Integral Method. The Brinkman equation (Brinkman- extended Darcy formulation with inertia! term included), which allows the no-slip boundary condition to be satisfied, is used as the starting momentum equation. With the inclusion of the continuity and energy equations the obtained set of governed non-linear partial differential equations allows the transport phenomena of natural convection in homogeneous isotropic non-deformable porous media, saturated with single phase Newtonian fluid, to be conveniently studied. The purpose of this work is to present flow and heat transfer characteristics of the fluid: a) in a vertical porous cavity with vertical walls maintained at different constant temperature where the top and the bottom walls are adiabatic, as well as b) in a horizontal porous layer heated from below. 1 Introduction Natural convection is one of the most frequently studied transport phenomena in porous media. Due to the general complexity of this phenomena, our studies are based on simplified mathematical model where we assume that: Transactions on Modelling and Simulation vol 24, © 1999 WIT Press, www.witpress.com, ISSN 1743-355X 146 Boundary Elements - porous media is a material consisting of a solid and a fluid phases. The solid phase is homogeneous, isotropic and non-deformable while the fluid is single phase, Newtonian and the fluid density is taken not to depend on pressure variations, but only on variations of temperature, - porous media is fluid saturated, - the solid and the fluid phase are in thermal equilibrium (7 = 7} =7^) and consequently the thermal behaviour of the porous media is described only by a single equation in terms of the average temperature T. When the temperature of the saturating fluid phase in a porous media is not uniform, some flows induced by buoyancy effects may occur. Commonly called free or natural convection movements, these flows depend on density differences due to temperature gradients and boundary conditions (Bear [1]). 2 Governing Equations The system of partial differential equations governing the phenomena of natural convection in porous media represents the basic conservation balances of mass, momentum and energy and consist of: - continuity equation • momentum equation 1 dv, 1 <?v/v, Darcy Law Brinkman extension - energy equation <5T dv T _ &T d [ „ dT where v- is filtration velocity; ^, /?, y stands for porosity, fluid density, kinematic viscosity of fluid, respectively; K is permeability of porous media, dPjdXi pressure gradient in the flow direction and s.j is the strain rate tensor defined as s^ = 1/2 (dv^fffxj +dvj/dx}. The kinematic viscosity is partitioned into constant and perturbated parts so that y = y + / ; T stands for temperature and O is thermal diffiisivity of the porous media defined as a^ -^/pCp , Transactions on Modelling and Simulation vol 24, © 1999 WIT Press, www.witpress.com, ISSN 1743-355X Boundary Elements 147 where ^ is heat conductivity of porous media /^ = (1 - ^)A +^A with 4 being heat conductivity of the solid part. The thermal diffiisivity a^ is also partitioned into constant and perturbated parts, i.e., a^ = a^ + a^ . The coefficient $) is so called heat capacity ratio, where p is density of fluid and solid phase and c^ is specific heat capacity. Brinkman extension, in our case written in a way to allow a changeable viscosity (as our attempt in the future is to deal also with non-Newtonian saturating fluids), expresses the viscous resistance or viscous drag force exerted by the solid phase on the flowing fluid at their contact surfaces. With the Brinkman equation we are able to satisfy the no-slip boundary conditions on an impermeable surfaces that bounds the porous media. It is important to stress out that the Brinkman equation is essentially an interpolation scheme between the Navier-Stokes and Darcy equation. It is well known that in the limit when the porosity approaches unity and consequently the permeability tends towards infinity, the Brinkman equation transforms into the classical Navier-Stokes equation for a pure fluid, meanwhile for the permeability converging to zero, the Brinkman term becomes negligible and the Darcy law is being recovered. With above formulated set of conservative equations named modified Navier- Stokes equations for porous media, we can solve in principal any transport phenomena in porous media if the appropriate hydrodynamic and thermal boundary conditions are precisely defined. 3 Velocity - Vorticity Formulation Governing equations may be written in different formulations as velocity- pressure, vorticity-stream function, velochy-vorticity, penalty formulations, ect. We have chosen to use the Velocity -Vorticity Formulation (WF) on account of its advantages when dealing with the Boundary Element Method (BEM). This formulation is well suited also for the Boundary Domain Integral Method (BDDM) which is an extension of the BEM. With the vorticity vector o)j = e^ dv^/dXj representing the curl of the velocity field, the computational scheme is partitioned into its kinematic and kinetic part so that the continuity and momentum equations are replaced by the equations of kinematics and kinetics (Skerget [2]). Applying the curl operator directly to the vorticity defined above, using the continuity equation and with addition of the relaxation parameter a , the kinematic can be formulated in the form of the parabolized kinematic equation written for the plane case of calculation as ^v, 1 <?V, do* dXdX a dt ** dx Transactions on Modelling and Simulation vol 24, © 1999 WIT Press, www.witpress.com, ISSN 1743-355X 148 Boundary Elements The kinetics is governed by the vorticity transport equation obtained as a curl of the momentum equation which reads: - 6v" - + — -,< .\ (5) ^ ^ ^ The new variable appears in vorticity equation, namely T^ = t/fi which is so called modified vorticity time step, introduced only as a necessary mathematical step allowing us to use the WF on our momentum equation. The last term of eqn. (5) represent a contribution from non-linear material properties defined as dx) dy dx ' dx dy By For the same reason as with vorticity kinetics we introduce new modified temperature time step also into energy equation ty = f/cr that permits us to write energy equation in the form 6T 4 Boundary Domain Integral Equations The obtained set of partial differential equations (4), (5), (6) have to be, using weighted residual techniques in combination with appropriate fundamental solutions - Green functions and completed with suitable boundary conditions, transformed into boundary domain integral equations. Integral representation can be derived by using integral representation of a parabolic diffusion-convective equation, rendering boundary domain integral equations for the plane (2D): - kinematics /9 ' - eg J a-j-dZl + plvj^ u <Kl describing the time dependent transport of velocity field v,, where w* is an elliptic modified Helmholtz fundamental solution which takes into account the effects of geometry, time step and material properties (Jecl [3]), - vorticity kinetics Transactions on Modelling and Simulation vol 24, © 1999 WIT Press, www.witpress.com, ISSN 1743-355X Boundary Elements 149 describing the time dependent transport of vorticity field CD in an integral formulation, where U* =y(/>u and u is the elliptic diffusion-convective fundamental solution of steady diffusion-convective PDE with first order reaction term (Skerget [2]), considering the effects of geometry, material properties, time step and velocity v, which has to be decomposed into an average constant vector and perturbation vector as v^ = v) + v," , - heat energy kinetics (9) "P a describing the time dependent transport of temperature in an integral formulation, where U* = au and %* is once again the elliptic diffusion-convective fundamental solution of steady diffiision-convective PDE with reaction term. 5 Discretized Boundary Domain Integral Equations For the numerical approximate solution of field functions (velocity, vorticity, temperature), the corresponding boundary domain integral representation is written in a discretized manner in which the integrals over the boundary and domain are approximated by a sum of integrals over E individual boundary elements and C internal cells (Hribersek [4]). The discretized equation for kinematic is written for all boundary nodes resulting in matrix system (10) to be solved for unknown boundary velocity components or their normal derivatives, while the computation of internal domain velocity components, if needed, is done afterwards in an explicit manner, point by point. Transactions on Modelling and Simulation vol 24, © 1999 WIT Press, www.witpress.com, ISSN 1743-355X 150 B oundary Elements The discretized equations for vorticity kinetic and for heat energy kinetic should be written for all boundary and internal nodes resulting in a matrix systems for vorticity kinetic yielding the solution of the unknown boundary vorticity flux values and unknown domain vorticity values. For the heat energy kinetic
Load more

Recommended publications
  • Chapter 6 Natural Convection in a Large Channel with Asymmetric
    Chapter 6 Natural Convection in a Large Channel with Asymmetric Radiative Coupled Isothermal Plates Main contents of this chapter have been submitted for publication as: J. Cadafalch, A. Oliva, G. van der Graaf and X. Albets. Natural Convection in a Large Channel with Asymmetric Radiative Coupled Isothermal Plates. Journal of Heat Transfer, July 2002. Abstract. Finite volume numerical computations have been carried out in order to obtain a correlation for the heat transfer in large air channels made up by an isothermal plate and an adiabatic plate, considering radiative heat transfer between the plates and different inclination angles. Numerical results presented are verified by means of a post-processing tool to estimate their uncertainty due to discretization. A final validation process has been done by comparing the numerical data to experimental fluid flow and heat transfer data obtained from an ad-hoc experimental set-up. 111 112 Chapter 6. Natural Convection in a Large Channel... 6.1 Introduction An important effort has already been done by many authors towards to study the nat- ural convection between parallel plates for electronic equipment ventilation purposes. In such situations, since channels are short and the driving temperatures are not high, the flow is usually laminar, and the physical phenomena involved can be studied in detail both by means of experimental and numerical techniques. Therefore, a large experience and much information is available [1][2][3]. In fact, in vertical channels with isothermal or isoflux walls and for laminar flow, the fluid flow and heat trans- fer can be described by simple equations arisen from the analytical solutions of the natural convection boundary layer in isolated vertical plates and the fully developed flow between two vertical plates.
    [Show full text]
  • Influence of Heat Source/Sink on Free Convection in Annular Porous Region
    International Journal of Heat and Technology Vol. 39, No. 3, June, 2021, pp. 841-850 Journal homepage: http://iieta.org/journals/ijht Influence of Heat Source/Sink on Free Convection in Annular Porous Region Anurag1, Shyam Lal Yadav2*, Ashok Kumar Singh1 1 Department of Mathematics, Institute of Science, Banaras Hindu University, Varanasi 221005, Uttar Pradesh, India 2 Department of Mathematics, Lok Nayak Jai Prakash Institute of Technology, Chhapra 841302, Bihar, India Corresponding Author Email: [email protected] https://doi.org/10.18280/ijht.390318 ABSTRACT Received: 4 January 2020 The significant interpretation of this model is to explore the influence of temperature- Accepted: 26 December 2020 dependent heat source/sink on laminar free-convective flow in an annular porous region such as petroleum engineering, thermal technique and groundwater hydrology. For a Keywords: unified solution of the Brinkman-Darcy model, the regulatory equations solved analytically fully developed flow, natural convection, by applying the variation of parameter technique in terms of Bessel's functions for the heat source and sink, isothermal and constant heat source and sink. Moreover, we have investigated the Variations of Darcy number, Heat flux, modified Bessel function source/sink and viscosity ratio in the presence of isothermal and constant heat flux sequentially. As a result, we received the critical value of the velocity for the radii ratio (R = 2.05 and 2.92) in both the cases of source and sink (S = 1.0 and Si = 0.1) respectively which is exhibited through the graphs. Further, the numerical outcomes present of the skin friction including volume flow with annular gap by the graphs as well as tables.
    [Show full text]
  • Forced and Natural Convection
    FORCED AND NATURAL CONVECTION Forced and natural convection ...................................................................................................................... 1 Curved boundary layers, and flow detachment ......................................................................................... 1 Forced flow around bodies .................................................................................................................... 3 Forced flow around a cylinder .............................................................................................................. 3 Forced flow around tube banks ............................................................................................................. 5 Forced flow around a sphere ................................................................................................................. 6 Pipe flow ................................................................................................................................................... 7 Entrance region ..................................................................................................................................... 7 Fully developed laminar flow ............................................................................................................... 8 Fully developed turbulent flow ........................................................................................................... 10 Reynolds analogy and Colburn-Chilton's analogy between friction and heat
    [Show full text]
  • Stability of Thermosolutal Free Convection on Electrically Conducting Superposed Fluid and Porous Layer
    International Journal of Applied Engineering Research ISSN 0973-4562 Volume 13, Number 17 (2018) pp. 13022-13030 © Research India Publications. http://www.ripublication.com Stability of Thermosolutal Free Convection on Electrically Conducting Superposed Fluid and Porous Layer K.Sumathi1*, T.Arunachalam2 and S.Aiswarya3 Department of Mathematics, 1,3PSGR Krishnammal College for Women, Coimbatore – 641 004, TamilNadu, India. 2Kumaraguru College of Technology, Coimbatore – 641 049, TamilNadu, India. (*Corresponding Author) Abstract investigated the onset of convection in a horizontal sparsely packed porous medium with imposed thermal and The effect of thermosolutal convection on linear stability of an concentration gradient at both the boundaries and found that incompressible fluid in a horizontal fluid/porous interface in Darcy number destabilizes the system in case of stationary the presence of a vertical magnetic field has been analyzed. and oscillatory convection. Asymptotic solution is obtained for rigid boundaries using normal mode analysis and the analysis is restricted to long Jena et al. [6] by employing Darcy-Brinkman-Forchheimer wave approximations. The influence of various non- model studied about the linear stability of onset of dimensional parameters such as Chandrasekhar number, thermosolutal convection in a rectangular annulus which is magnetic Prandtl number, Schmidt number, Prandtl number, highly heated and soluted from within and exposed to low porosity, wave number and depth ratio on stability thermal and solutal gradients from outside in a non-Darcian characteristics of flow field are represented numerically. region. Abdullah et al. [2] studied the effect of thermosolutal convection of a viscous incompressible fluid in the presence Keywords: Thermosolutal convection, magnetic field, fluid- of rotation on the stability of the system and observed that porous interface, stability analysis.
    [Show full text]
  • Convectional Heat Transfer from Heated Wires Sigurds Arajs Iowa State College
    Iowa State University Capstones, Theses and Retrospective Theses and Dissertations Dissertations 1957 Convectional heat transfer from heated wires Sigurds Arajs Iowa State College Follow this and additional works at: https://lib.dr.iastate.edu/rtd Part of the Physics Commons Recommended Citation Arajs, Sigurds, "Convectional heat transfer from heated wires " (1957). Retrospective Theses and Dissertations. 1934. https://lib.dr.iastate.edu/rtd/1934 This Dissertation is brought to you for free and open access by the Iowa State University Capstones, Theses and Dissertations at Iowa State University Digital Repository. It has been accepted for inclusion in Retrospective Theses and Dissertations by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. CONVECTIONAL BEAT TRANSFER FROM HEATED WIRES by Sigurds Arajs A Dissertation Submitted to the Graduate Faculty in Partial Fulfillment of The Requirements for the Degree of DOCTOR 01'' PHILOSOPHY Major Subject: Physics Signature was redacted for privacy. In Chapge of K or Work Signature was redacted for privacy. Head of Major Department Signature was redacted for privacy. Dean of Graduate College Iowa State College 1957 ii TABLE OF CONTENTS Page I. INTRODUCTION 1 II. THEORETICAL CONSIDERATIONS 3 A. Basic Principles of Heat Transfer in Fluids 3 B. Heat Transfer by Convection from Horizontal Cylinders 10 C. Influence of Electric Field on Heat Transfer from Horizontal Cylinders 17 D. Senftleben's Method for Determination of X , cp and ^ 22 III. APPARATUS 28 IV. PROCEDURE 36 V. RESULTS 4.0 A. Heat Transfer by Free Convection I4.0 B. Heat Transfer by Electrostrictive Convection 50 C.
    [Show full text]
  • Forced Convection Heat Transfer Convection Is the Mechanism of Heat Transfer Through a Fluid in the Presence of Bulk Fluid Motion
    Forced Convection Heat Transfer Convection is the mechanism of heat transfer through a fluid in the presence of bulk fluid motion. Convection is classified as natural (or free) and forced convection depending on how the fluid motion is initiated. In natural convection, any fluid motion is caused by natural means such as the buoyancy effect, i.e. the rise of warmer fluid and fall the cooler fluid. Whereas in forced convection, the fluid is forced to flow over a surface or in a tube by external means such as a pump or fan. Mechanism of Forced Convection Convection heat transfer is complicated since it involves fluid motion as well as heat conduction. The fluid motion enhances heat transfer (the higher the velocity the higher the heat transfer rate). The rate of convection heat transfer is expressed by Newton’s law of cooling: q hT T W / m 2 conv s Qconv hATs T W The convective heat transfer coefficient h strongly depends on the fluid properties and roughness of the solid surface, and the type of the fluid flow (laminar or turbulent). V∞ V∞ T∞ Zero velocity Qconv at the surface. Qcond Solid hot surface, Ts Fig. 1: Forced convection. It is assumed that the velocity of the fluid is zero at the wall, this assumption is called no‐ slip condition. As a result, the heat transfer from the solid surface to the fluid layer adjacent to the surface is by pure conduction, since the fluid is motionless. Thus, M. Bahrami ENSC 388 (F09) Forced Convection Heat Transfer 1 T T k fluid y qconv qcond k fluid y0 2 y h W / m .K y0 T T s qconv hTs T The convection heat transfer coefficient, in general, varies along the flow direction.
    [Show full text]
  • Chaotic Flow in a 2D Natural Convection Loop
    International Journal of Heat and Mass Transfer 61 (2013) 565–576 Contents lists available at SciVerse ScienceDirect International Journal of Heat and Mass Transfer journal homepage: www.elsevier.com/locate/ijhmt Chaotic flow in a 2D natural convection loop with heat flux boundaries ⇑ William F. Louisos a,b, , Darren L. Hitt a,b, Christopher M. Danforth a,c a College of Engineering & Mathematical Sciences, The University of Vermont, Votey Building, 33 Colchester Avenue, Burlington, VT 05405, United States b Mechanical Engineering Program, School of Engineering, The University of Vermont, Votey Building, 33 Colchester Avenue, Burlington, VT 05405, United States c Department of Mathematics & Statistics, Vermont Complex Systems Center, Vermont Advanced Computing Core, The University of Vermont, Farrell Hall, 210 Colchester Avenue, Burlington, VT 05405, United States article info abstract Article history: This computational study investigates the nonlinear dynamics of unstable convection in a 2D thermal Received 13 August 2012 convection loop (i.e., thermosyphon) with heat flux boundary conditions. The lower half of the thermosy- Received in revised form 1 February 2013 phon is subjected to a positive heat flux into the system while the upper half is cooled by an equal-but- Accepted 3 February 2013 opposite heat flux out of the system. Water is employed as the working fluid with fully temperature dependent thermophysical properties and the system of governing equations is solved using a finite vol- ume method. Numerical simulations are performed for varying levels of heat flux and varying strengths Keywords: of gravity to yield Rayleigh numbers ranging from 1.5 Â 102 to 2.8 Â 107.
    [Show full text]
  • Numerical Simulation of Williamson Combined Natural and Forced
    entropy Article Numerical Simulation of Williamson Combined Natural and Forced Convective Fluid Flow between Parallel Vertical Walls with Slip Effects and Radiative Heat Transfer in a Porous Medium Mohammad Yaghoub Abdollahzadeh Jamalabadi 1,*, Payam Hooshmand 2, Navid Bagheri 3, HamidReza KhakRah 4 and Majid Dousti 5 1 Department of Mechanical, Robotics and Energy Engineering, Dongguk Universit, Seoul 04620, Korea 2 Young Researchers and Elite Club, Mahabad Branch, Islamic Azad University, Mahabad 433-59135, Iran; [email protected] 3 Department of Energy Engineering, Graduate School of the Environment and Energy, Science and Research Branch, Islamic Azad University, Tehran 1477893855, Iran; [email protected] 4 Department of Mechanical Engineering, College of Engineering, Shiraz Branch, Islamic Azad University, Shiraz 71987-74731, Iran; [email protected] 5 Faculty of Engineering, Zanjan University, Zanjan 38791-45371, Iran; [email protected] * Correspondence: [email protected] or [email protected]; Tel.: +82-22-260-3073 Academic Editor: Milivoje M. Kostic Received: 29 February 2016; Accepted: 8 April 2016; Published: 18 April 2016 Abstract: Numerical study of the slip effects and radiative heat transfer on a steady state fully developed Williamson flow of an incompressible Newtonian fluid; between parallel vertical walls of a microchannel with isothermal walls in a porous medium is performed. The slip effects are considered at both boundary conditions. Radiative highly absorbing medium is modeled by the Rosseland approximation. The non-dimensional governing Navier–Stokes and energy coupled partial differential equations formed a boundary problem are solved numerically using the fourth order Runge–Kutta algorithm by means of a shooting method. Numerical outcomes for the skin friction coefficient, the rate of heat transfer represented by the local Nusselt number were presented even as the velocity and temperature profiles illustrated graphically and analyzed.
    [Show full text]
  • Finite Element Solutions for Non-Newtonian Pulsatile Flow in a Non-Darcian Porous Medium Conduit ∗
    Nonlinear Analysis: Modelling and Control, 2007, Vol. 12, No.3, 317–327 Finite Element Solutions for Non-Newtonian Pulsatile Flow in a Non-Darcian Porous Medium Conduit ∗ R. Bhargava1, H. S. Takhar2, S. Rawat1, Tasveer A. Bég3, O. Anwar Bég4 1Mathematics Department, Indian Institute of Technology Roorkee-247667, India [email protected] 2Engineering Department, Manchester Metropolitan University Oxford Rd., Manchester, M5 1GD, UK [email protected] 3Engineering Mechanics and Earthquake Dynamics Research 18 Milton Grove, Manchester, M16 OBP [email protected] 4Fire Engineering Sciences, Leeds College of Building / Leeds Metropolitan University North Street, Leeds, LS2 7QT, UK [email protected] Received: 17.05.2006 Revised: 08.11.2006 Published online: 31.08.2007 Abstract. The present analysis is motivated by the need to elucidate with more accuracy and sophistication the hydrodynamics of non-Newtonian flow via a channel containing a porous material under pulsating pressure gradient. A one-dimensional transient rheological model for pulsating flow through a Darcy-Forcheimmer porous channel is used. A modified Casson non-Newtonian constitutive model is employed for the transport fluid with a drag force formulation for the porous body force effects. The model is transformed and solved using a finite element numerical technique. Rheological effects are examined using a β parameter which vanishes in the limit (Newtonian flow). Velocity profiles are plotted for studying the influence of Reynolds number, Darcy number, Forchheimer number and the β (non-Newtonian) parameter. The channel considered is rigid with a pulsatile pressure applied via an appropriate pressure gradient term. The model finds applications in industrial filtration systems, pumping of polymeric fluids etc.
    [Show full text]
  • Enhancement of Natural Convection Heat Transfer Coefficient by Using V-Fin Array
    International Journal of Engineering Research and General Science Volume 3, Issue 2, March-April, 2015 ISSN 2091-2730 Enhancement of Natural convection heat transfer coefficient by using V-fin array Rameshwar B. Hagote, Sachin K. Dahake Student of mechanical Engg. Department, MET’s IOE, Adgaon, Nashik (Maharashtra,India). Email. Id: [email protected], Contact No. +919730342211 Abstract— Extended surfaces known as fins are, used to enhance convective heat transfer in a wide range of engineering applications, and offer an economical and trouble free solution in many situations demanding natural convection heat transfer. Fin arrays on horizontal, inclined and vertical surfaces are used in variety of engineering applications to dissipate heat to the surroundings. Studies of heat transfer and fluid flow associated with such arrays are therefore of considerable engineering significance. The main controlling variables generally available to the designer are the orientation and the geometry of the fin arrays. An experimental work on natural convection adjacent to a vertical heated plate with a multiple V- type partition plates (fins) in ambient air surrounding is already done. Boundary layer development makes vertical fins inefficient in the heat transfer enhancement. As compared to conventional vertical fins, this V-type partition plate works not only as extended surface but also as flow turbulator. This V-type partition plate is compact and hence highly economical. The numerical analysis of this technique is done using Computational Fluid Dynamics (CFD) software, Ansys CFX , for natural convection adjacent to a vertical heated plate in ambient air surrounding. In numerical analysis angle of V-fin is further optimized for maximum average heat transfer coefficient.
    [Show full text]
  • Chapter 7. Convection and Complexity
    Chapter 7 Convection and complexity ... if your theory is found to be convection has been taken as the classic exam­ against the second law of ple of thermal convection, and the hexagonal planform has been considered to be typical of thermodynamics, I can give you no convective patterns at the onset of thermal con­ hope; there is nothing for it but to vection. Fifty years went by before it was real­ collapse in deepest humiliation. ized that Benard's patterns were actually driven from above, by surface tension. not from below by Eddington an unstable thermal boundary layer. Experiments Contrary to current textbooks ... the showed the san1.e style of convection when the fluid was heated from above, cooled from below observed world does not proceed or when performed in the absence of gravity. This from lower to higher "degrees of confirmed the top-down surface-driven nature of disorder", since when all the convection which is now called Marangoni or B€mard-Marangoni convection. gravitationally-induced phenomena Although it is not generally recognized as are taken into account the emerging such, mantle convection is a branch of the newly result indicates a net decrease in the renamed science of complexity. Plate tectonics may be a self-drivenfarjrom-equilibrium system that orga­ "degrees of disorder", a greater nizes itself by dissipation in and between the "degree of structuring" ... classical plates, the mantle being a passive provider of equilibrium thermodynamics ... has energy and material. Far-from-equilibrium sys­ to be completed by a theory of tems, particularly those in a gravity field, can locally evolve toward a high degree of order.
    [Show full text]
  • Effect of Darcy, Reynolds, and Prandtl Numbers on the Performance of Two-Phase Flow Heat Exchanger Filled with Porous Media
    Heat Transfer—Asian Research, 43 (8), 2014 Effect of Darcy, Reynolds, and Prandtl Numbers on the Performance of Two-Phase Flow Heat Exchanger Filled with Porous Media M. Tarawneh,1 A.S. Alshiqirate,2 and A.M. Jawarneh1 1Mechanical Engineering, The Hashemite University, Zarqa, Jordan 2Mechanical Engineering, Alshoubak University College, Al-Balqa’ Applied University, Alsalt, Jordan Experimental investigation of two-phase laminar forced convection in a single porous tube heat exchanger is presented. The effect of Darcy, Reynolds, and Prandtl numbers on the performance of this heat exchanger during the condensation process of carbon dioxide at different test conditions were investigated. Gravel sand with different porosities is used as a porous medium. The flow in the porous medium is modeled using the Brinkman–Forchheimer-extended Darcy model. Parametric studies are also conducted to evaluate the effects of porosity and Reynolds and Prandtl numbers on the heat transfer coefficient and the friction factor. A dimensionless performance parameter is developed in order to be used in evaluating the porous tube heat exchanger based on both the heat transfer enhancement and the associated pressure. drop. The study covers a wide range of inlet pressure (Pin), mass flow rate ε ≤ ≤ (m), porosity of gravel. sand ( ), and Darcy number (Da) which ranged: 34.5 Pin 43 bars, 8*10–5 ≤ m ≤ 16*10–5 kg/s, 34.9% ≤ ε ≤ 44.5%, 1.6*10–6 ≤ Da ≤ 5*10–6, respectively. The study predicted the combined effect of the Reynolds number, Darcy number, porosity, and Prandtl number on the heat transfer and pressure drop of carbon dioxide during the condensation process in a porous medium.
    [Show full text]