Thermodynamic analysis of geothermal heat pump during the cold season

Similar documents
Studies on Nuclear Fuel Rod Thermal Performance

First Law of Thermodynamics. Control Mass (Closed System) Conservation of Mass. Conservation of Energy

Simulated Performance of Packed Bed Solar Energy Storage System having Storage Material Elements of Large Size - Part I

Entropy ISSN

Department of Mechanical Engineering ME 322 Mechanical Engineering Thermodynamics. Lecture 33. Psychrometric Properties of Moist Air

Thermal Performance of Electrocaloric Refrigeration using Thermal Switches of Fluid Motion and Changing Contact Conductance

Simulation of Fluid Flow and Heat Transfer in Porous Medium Using Lattice Boltzmann Method

A study of fluid flow simulation in convergentdivergent

Chapter 4 The second law of thermodynamics

Psychrometric Applications

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM

New Expansion and Infinite Series

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

CBE 291b - Computation And Optimization For Engineers

Study Guide Final Exam. Part A: Kinetic Theory, First Law of Thermodynamics, Heat Engines

CHAPTER 20: Second Law of Thermodynamics

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

A Brief Note on Quasi Static Thermal Stresses In A Thin Rectangular Plate With Internal Heat Generation

Consequently, the temperature must be the same at each point in the cross section at x. Let:

CONTRIBUTION TO THE EXTENDED DYNAMIC PLANE SOURCE METHOD

Travelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing

Freely propagating jet

- 5 - TEST 2. This test is on the final sections of this session's syllabus and. should be attempted by all students.

A5682: Introduction to Cosmology Course Notes. 4. Cosmic Dynamics: The Friedmann Equation. = GM s

HT Module 2 Paper solution. Module 2. Q6.Discuss Electrical analogy of combined heat conduction and convection in a composite wall.

Basic model for traffic interweave

1. a) Describe the principle characteristics and uses of the following types of PV cell: Single Crystal Silicon Poly Crystal Silicon

An inverse steady state thermal stresses in a thin clamped circular plate with internal heat generation

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

Effects of peripheral drilling moment on delamination using special drill bits

Energy Consideration


G. MATEESCU 1 A. MATEESCU 1 C. SAMOILĂ 2 D. URSUŢIU 2

The heat budget of the atmosphere and the greenhouse effect

The International Association for the Properties of Water and Steam. Release on the Ionization Constant of H 2 O

Module 2: Rate Law & Stoichiomtery (Chapter 3, Fogler)

Transient Aspects of Heat Flux Bifurcation in Porous Media: An Exact Solution

Optimum Design of a Horizontal Condenser

PART 1 MULTIPLE CHOICE Circle the appropriate response to each of the questions below. Each question has a value of 1 point.

Deteriorating Inventory Model for Waiting. Time Partial Backlogging

On the Linear Stability of Compound Capillary Jets

New data structures to reduce data size and search time

CHM Physical Chemistry I Chapter 1 - Supplementary Material

Some parameters of varicaps with gradient base area based on Shottky barrier

Topic 1 Notes Jeremy Orloff

1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation

The Thermodynamics of Aqueous Electrolyte Solutions

Lesson 8. Thermomechanical Measurements for Energy Systems (MENR) Measurements for Mechanical Systems and Production (MMER)

Minimum Energy State of Plasmas with an Internal Transport Barrier

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

Estimation of Global Solar Radiation at Onitsha with Regression Analysis and Artificial Neural Network Models

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011

Study on the Calculation of Magnetic Force Based on the Equivalent Magnetic Charge Method

In the next lecture... Tutorial on ideal cycles and component performance.

Measuring Electron Work Function in Metal

Chaos in drive systems

TANDEM QUEUE WITH THREE MULTISERVER UNITS AND BULK SERVICE WITH ACCESSIBLE AND NON ACCESSBLE BATCH IN UNIT III WITH VACATION

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

13.4 Work done by Constant Forces

MAC-solutions of the nonexistent solutions of mathematical physics

Ordinary differential equations

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018

Terminal Velocity and Raindrop Growth

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

SUPPLEMENTARY INFORMATION

JJMIE Jordan Journal of Mechanical and Industrial Engineering

Set up Invariable Axiom of Force Equilibrium and Solve Problems about Transformation of Force and Gravitational Mass

Improper Integrals, and Differential Equations

Effects of dry density on soil water characteristic curve of clay

Construction and Selection of Single Sampling Quick Switching Variables System for given Control Limits Involving Minimum Sum of Risks

Determination of the activation energy of silicone rubbers using different kinetic analysis methods

Application of Exp-Function Method to. a Huxley Equation with Variable Coefficient *

Sections 1.3, 7.1, and 9.2: Properties of Exponents and Radical Notation

APPLICATION OF WATER EQUIVALENTS METHOD TO CALCULATION OF CERAMIC HEAT EXCHANGERS

An Improved Model for Predicting Fluid Temperature in Deep Wells

A smoothed particle hydrodynamics method for evaporating. multiphase flows

LABYRINTH SEALS DYNAMIC COEFFICIENTS AND CRITICAL SPEED RESEARCH BASED ON CHILDS MODEL

Experimental and theoretical study of a humidification dehumidification water desalination system using solar energy

arxiv: v1 [math.gm] 30 Dec 2015

DIRECT CURRENT CIRCUITS

HOMEWORK SOLUTIONS MATH 1910 Sections 7.9, 8.1 Fall 2016

13: Diffusion in 2 Energy Groups

Numerical Study Of Coated Electrical Contacts

Factors affecting the phonation threshold pressure and frequency

df dx There is an infinite number of different paths from

Pressure Wave Analysis of a Cylindrical Drum

DISTRIBUTION OF SUB AND SUPER HARMONIC SOLUTION OF MATHIEU EQUATION WITHIN STABLE ZONES

Vadose Zone Hydrology

Predict Global Earth Temperature using Linier Regression

HYDROTHERMAL ANALYSIS OF THE ABSORBER TUBES USED IN LINEAR FRESNEL REFLECTOR SOLAR THERMAL SYSTEM

ScienceDirect. Thermal convection measurements inside aluminum foam and comparison to existing analytical solutions

Review of Calculus, cont d

Math 124A October 04, 2011

Conservation Law. Chapter Goal. 5.2 Theory

INVESTIGATION OF BURSA, ESKIKARAAGAC USING VERTICAL ELECTRICAL SOUNDING METHOD

Strategy: Use the Gibbs phase rule (Equation 5.3). How many components are present?

Modified Crout s method for an LU decomposition of an interval matrix

Newly Established Hydrometer Calibration Set Up at UM. Ümit Y. AKÇADAĞ, S. Eren SAN TÜBİTAK, Ulusal Metroloji Enstitüsü (UME) Gebze Kocaeli, Turkey

Overview. Before beginning this module, you should be able to: After completing this module, you should be able to:

Analysis on Parameters of Regeneration Subsystem in Liquid Desiccant Dehumidification Systems

Transcription:

IOP onference eries: Mterils cience nd ngineering PAPR OPN A hermodynmic nlysis of geotherml het pump during the cold seson o cite this rticle: G Dumitrcu et l 2016 IOP onf. er.: Mter. ci. ng. 147 012137 View the rticle online for updtes nd enhncements. Relted content - hermodynmic Anlysis on Vpor Phse pitxy of GAs by Gl 3 nd AsH 3 ystem Akinori oukitu, Fumio Hsegw nd Hisshi eki - hermodynmic Anlysis of In x G 1-x N Alloy omposition Grown by Metlorgnic Vpor Phse pitxy Akinori oukitu, Noyuki khshi, etsuy ki et l. - A hermodynmic Anlysis of the Vlidity of Wenzel nd ssie's qutions n hui-xi, Lu Xio-Ying, Li Wen et l. his content ws downloded from IP ddress 37.44.201.67 on 31/12/2017 t 16:47

7th Interntionl onference on Advnced oncepts in Mechnicl ngineering IOP onf. eries: Mterils cience nd ngineering 147 (2016) 012137 doi:10.1088/1757-899x/147/1/012137 hermodynmic nlysis of geotherml het pump during the cold seson G Dumitrşcu 1, A Dumencu 1, B Horbniuc 1 nd M V Atnsiu 1 1 Automotive nd Mechnicl ngineering Deprtment, Gheorghe Aschi echnicl University of Isi, Isi, Romni -mil: dumencu@yhoo.com* Abstrct. he pper is nlysing the performnces (OP, power nd, heting het rte function of time) for ground-coupled het pump tht is used to het spce during winter, for period of 180 dys. he nlysis purpose is to evlute the time bsed chnges in vlues of OP nd, energy trnsfers of geotherml het pump, considering scenrio for the vrition of the mbient temperture in time nd n nlyticl solution for the time dependence of the soil one. he tempertures nd the energy trnsfer rtes were determined on the bsis of the eversible entropy blnce eqution. 1. Introduction he ground-coupled het pump bsed energy systems re using geotherml underground energy for spce heting nd cooling. Becuse of lrge heting efficiency nd environmentl protection, the interest in these systems is incresing tody. here is wide vriety of these systems to be used in office buildings, houses nd historicl buildings. he study[1], regrding the office building heting systems using ground coupled het pumps with different borehole dimeters, showed tht the highest svings to investment rtio, during period of thirty yers, is 4.80. he pper [2]presents n open loop system tht uses for cooling wter lke het cpcity. For heting most het pump systems use closed loop. A generl review[3], presenting the most common het pumps bsed energy systems tht re used tody, is compring their performnces nd opertionl prmeters but for qusi-stedy stte opertion. he thermodynmic nlysis of ground-coupled het pump bsed energy systems emphsizes the importnt influence of underground soil temperture profile nd the geometry of the het exchngers extrcting the geotherml energy, e.g. either the superficil temperture profile of soil influenced by the nturl convection on the surfce [4], or tht fter depth of 15m, the initil underground temperture is lmost constnt. his pper perform thermodynmic nlysis of closed loop het pump bsed heting system coupling two externl het reservoirs with different vrible fetures. he het sink is the heted spce tht sks heting rte function of the environmentl temperture nd the het source is the soil with vrible temperture in time. he het pump cycle is n eversible one, nd the overll entropy blnce eqution llowed linking these two het reservoirs by choosing scenrios for the heting rte during the winter nd n nlyticl solution in time for the geotherml het rte. he het pump bsed heting system is presented in figure 1. ontent from this work my be used under the terms of the retive ommons Attribution 3.0 licence. Any further distribution of this work must mintin ttribution to the uthor(s) nd the title of the work, journl cittion nd DOI. Published under licence by Ltd 1

7th Interntionl onference on Advnced oncepts in Mechnicl ngineering IOP onf. eries: Mterils cience nd ngineering 147 (2016) 012137 doi:10.1088/1757-899x/147/1/012137 Figure 1. he heting system scheme.h heted spce; evportor; condenser; compressor; V throttling vlve; i nd o input nd output tempertures for ; i nd o input nd output tempertures for. 2. Restrictive conditions he heted spce temperture H = 295.15 = const. while outside mbient temperture is vryingfrom253.15 to 288.15 (from 20 to +15 ). For the mbient temperture ws dopted scenrio, respectively we considered continuous heting during the winter (180 dys) with time dependence, see figure 2 where t is the time in hours. t 288.15 35 sin 4320 (1) Figure 2. = (t). onstnt initil soil temperture i = 283.15. For the soil temperture vrition in time ws imposed n nlyticl dependence, see figure 3: 2

7th Interntionl onference on Advnced oncepts in Mechnicl ngineering IOP onf. eries: Mterils cience nd ngineering 147 (2016) 012137 doi:10.1088/1757-899x/147/1/012137 where t is the time in hours. 0.375 283.15 exp 0.0025 t (2) Figure 3. = (t). Het crriers were ssumed to hve constnt het cpcities, = 500 W/ = const. nd = 1000 W/ = const. he heting system uses the ir s het crrier flowing in loop between the condenser nd the heted spce. Het pump het exchngers (, nd soil borehole) were presumed to hve constnt effectiveness (ε = 0.75 = const., ε = 0.75 = const., nd ε = 0.95 = const.). 3. Mthemticl lgorithm 3.1. he heted spce nd the condenser he heted spce needs to be mintined t constnt temperture, H = 295.15, while outside mbient temperture is vryingfrom253.15 to 288.15. hus the heting rte supplied by equlizes the lost het rte, see figure 4. Figure 4. herml energy blnce of the heted spce. he heted spce nd condenser het blnce equtions re: lost H H h A (3) lost lost N h h H H H A A H H H H N o N i N k k o on i in i. (4) 3

When the mbient temperture fluctutes between +15 nd 20 o, the ir input/output tempertures re: i = in = H = 295.15 = const. nd o = o ( ) function of the heting system design tempertures. We imposed on = 313.15 t N = 253.15. he dependences nd, = ( ) nd, o = o ( ) re obtined on the bsis of the het blnce equtions1, 3 nd 4, see figure 5: in on N H H i o 7 3.64286 421, [] (5) i o i, [] (6) i i o, [W]. (7) () (b) (c) Figure 5. he dependences () o = o ( ), (b) = ( ), (c) for 253.15 288.15. 3.2. he soil nd the evportor he het blnce equtions relted to the therml connection evportor soil re: i o i, [W] (8) o, [W]. (9) From equtions 2, 8, nd 9 we derived: i, [] (10) o, [] (11), [W]. (12) 7th Interntionl onference on Advnced oncepts in Mechnicl ngineering IOP onf. eries: Mterils cience nd ngineering 147 (2016) 012137 doi:10.1088/1757-899x/147/1/012137 4

7th Interntionl onference on Advnced oncepts in Mechnicl ngineering IOP onf. eries: Mterils cience nd ngineering 147 (2016) 012137 doi:10.1088/1757-899x/147/1/012137 3.3. Het pump eversibility he het pump eversibility (Irr) considered tht the isentropic efficiency of the compressor is 0.95, nd pressures drops in evportor nd condenser, on the suction nd dischrge nd liquid lines re of 0.1. Figure 6[5] present bsic eversible reverse cycle, consisting of processes: 1 2r eversible dibtic compression, 2r 3r eversible cooling, 3r 4r eversible dibtic process (expnsion or throttling), 4r 1 eversible heting. he internl eversibility is lso emphsized by the entropy vritions of the working fluid during the het exchnges with externl het reservoirs. Hence, the entropy vritions re: s s is the entropy vrition long the heting process 4r-1t (4-1) s s4 1 1t 4r 12 23 34 41 23 s2r s3t s41 s s s s s41 1 s41 N (13) s 41 is the entropy vrition long the cooling process 2r-3t (2-3). 2t 1 2 s 2r s 3t 3r 4 1 s 2 3 s 4t 4r 1t 1 3 4 s s Figure 6. cheme of n eversible reverse cycle, temperture entropy digrm. he eversible entropy internlly generted by flow with friction is stored in the entropy vrition during the cooling process, nd ccordingly we cn identify it in the first lw efficiency (OP): OP 23 W 23 23 41 m 2r3t mq m s 2r3t mq 23 s m 23 4r1t mq s 41 2r3t mq 4r1t mq 2r3t mq 4r1t mq N N 1 Irr (14) Irr 1 where: mq is the men thermodynmic temperture during the heting nd cooling processes, nd nd re the condensing nd evporting ones in nd. 2r3t s int mq N 1 nd Irr N, Irr re the numbers of internl s 4r1t 41 mq eversibility relted to the chosen reference tempertures of processes 2r-3t nd 4r-1t. For n idel reverse cycle, both numbers of eversibility re equl to one, i.e. rnot. hus the 5

7th Interntionl onference on Advnced oncepts in Mechnicl ngineering IOP onf. eries: Mterils cience nd ngineering 147 (2016) 012137 doi:10.1088/1757-899x/147/1/012137 entropy blnce eqution becomes [6]: Irr Irr OP 1 OP nd, by using oolpck we get for R717 (mmoni): or simplified: OPR (15) 0, where 4 1 nd, 2 3 (16) OPR 1.5672 2.1225 1.56564 (17) OPR 568338 0.60229 1., (errors: +1.1%, 1.8%). (18) olving the entropy blnce eqution it is obtined the reference temperture (t). his issue llows to evlute i, o nd from equtions (14), (15) nd (16). 2 3.4. Power blnce nd, first lw efficiency W (19) OP W. (20) 4. Numericl results his simultion ws developed for period of 6 months, round 180 dys. ffectiveness of het exchngers, condenser nd, evportor nd, soil het exchnger were ssumed to be ε = 0.75 = const., ε = 0.75 = const., nd ε = 0.95 = const. Initil soil temperture for this simultion is 10, i = 283.15. Underground het exchnger is considered to be borehole50m depth nd rdius of 0.2m. Overll het trnsfer surfce for extrcting het by convection is A = 62.8m 2. olving the problem strts by using equtions from evportor temperture section, simulting different time steps. Results ccording to time for soil temperture ( ), mbient temperture ( ), output temperture from condenser ( o ), temperture inside condenser ( ), temperture input to evportor ( i ), temperture output from evportor ( o ) nd temperture inside evportor re presented in tble 1. he grph in figure 7 indictes the tempertures ccording to time dependences. ble 2 presents the het rte inside condenser, evportor nd the power of the het pump. Het energy rtes function of time, re presented in figure 8. heir dependences with mbient temperture/time show tht het rtes in condenser, evportor nd, the power hve their peks during the lowest mbient temperture. In figure 9 is shown the OP function of time. 6

7th Interntionl onference on Advnced oncepts in Mechnicl ngineering IOP onf. eries: Mterils cience nd ngineering 147 (2016) 012137 doi:10.1088/1757-899x/147/1/012137 ble 1. empertures in reltion to time. ime [h] [] [] o [] [] i [] o [] [] 1 282.44 288.12 298.16 299.16 282.36 280.98 280.51 50 280.09 286.87 298.69 299.87 280.01 278.40 277.86 100 279.19 285.60 299.23 300.60 279.10 277.25 276.64 200 278.03 283.07 300.32 302.04 277.91 275.61 274.84 300 277.20 280.57 301.39 303.47 277.05 274.31 273.40 400 276.53 278.11 302.45 304.88 276.36 273.20 272.15 500 275.96 275.70 303.48 306.26 275.77 272.20 271.01 1000 273.86 264.88 308.12 312.44 273.58 268.31 266.55 1500 272.37 257.10 311.45 316.89 272.03 265.65 263.53 1750 271.74 254.69 312.48 318.26 271.39 264.69 262.46 2160 270.82 253.15 313.15 319.15 270.46 263.59 261.30 2320 270.49 253.38 313.04 319.01 270.13 263.31 261.04 2820 269.56 257.10 311.45 316.89 269.23 262.95 260.86 3320 268.72 264.88 308.12 312.44 268.45 263.33 261.62 3820 267.96 275.70 303.48 306.26 267.78 264.36 263.22 3920 267.82 278.11 302.45 304.88 267.66 264.63 263.63 4020 267.68 280.57 301.39 303.47 267.54 264.93 264.06 4120 267.54 283.07 300.32 302.04 267.42 265.24 264.51 4220 267.40 285.60 299.23 300.60 267.31 265.57 264.99 4270 267.33 286.87 298.69 299.87 267.25 265.74 265.23 4320 267.27 288.15 298.15 299.15 267.20 265.91 265.48 Figure 7. Opertionl tempertures(,, o,, i, o, )function of time. 7

7th Interntionl onference on Advnced oncepts in Mechnicl ngineering IOP onf. eries: Mterils cience nd ngineering 147 (2016) 012137 doi:10.1088/1757-899x/147/1/012137 ble 2. Het energy rtes functions of time. ime [h] [W] [W] W [W] OP 1 1505.45 1388.39 117.05 12.86 50 1772.64 1611.37 161.27 10.99 100 2044.93 1843.18 201.75 10.13 200 2586.98 2298.51 288.47 8.96 300 3123.29 2739.16 384.12 8.13 400 3651.02 3162.88 488.14 7.47 500 4167.38 3567.95 599.42 6.95 1000 6485.96 5271.48 1214.48 5.34 1500 8152.58 6375.40 1777.17 4.58 1750 8669.08 6691.37 1977.71 4.38 2160 9000 6868.32 2131.67 4.22 2320 8949.28 6823.24 2126.04 4.21 2820 8152.58 6275.03 1877.55 4.34 3320 6485.96 5125.65 1360.31 4.76 3820 4167.38 3423.89 743.48 5.60 3920 3651.02 3026.01 625 5.84 4020 3123.29 2611.89 511.40 6.10 4120 2586.98 2183.14 403.84 6.40 4220 2044.93 1741.59 303.33 6.74 4270 1772.64 1516.65 255.99 6.92 4320 1500 1289.28 210.71 7.12 Figure 8. Het rte in evportor, condenser nd power in time. 8

7th Interntionl onference on Advnced oncepts in Mechnicl ngineering IOP onf. eries: Mterils cience nd ngineering 147 (2016) 012137 doi:10.1088/1757-899x/147/1/012137 Figure 9. Het pump OP function of time 5. onclusions he pper dels with entropy blnce eqution in order to evlute the time dependent performnces of heting system, during the winter seson, involving geotherml het pump. he thermodynmic ssessment considered restrictive conditions describing imposed time bsed evolution of het pump externl het reservoirs. he entropy blnce eqution of the eversible het pump llowed complete linking of ll vrible in time prmeters, tempertures, energy trnsfer rtes nd OP. he mthemticl model of this ssessment might be pplied to other rel time bsed restrictive conditions. he geotherml het pump bsed heting systems seem to be very competitive to cogenertion bsed ones. hus, for this ssessment, the men vlue of the OP is round 5.46. If we evlute the first lw efficiency relted to the fossil bsed het energy operting the engine delivering the power for the het pump, we obtin: FLHP, men OPmen * 5.46* 0.35 to 0.45 1.911 to 2.457. Here η is the first lw efficiency of the engine delivering the power to the het pump. Nomenclture temperture [] h convection het coefficient [W/(m 2 )] hetrte [W] W power consumed [W] A surfce [m 2 ] ε het exchnger efficiency [ - ] m mss flow [kg/h] het cpcity of het crrier [W/] s entropy [J/] Irr eversibility OP coefficient of performnce ubscripts: mbient evportor condenser o output i input N nominl hc het crrier soil H heted spce 6. References [1] Luo J, Rohn J, Byer M nd Priess A 2013 herml performnce nd economic evlution of double U-tube borehole het exchnger with three different borehole dimeters nergy nd Buildings 67 pp 217-224 9

7th Interntionl onference on Advnced oncepts in Mechnicl ngineering IOP onf. eries: Mterils cience nd ngineering 147 (2016) 012137 doi:10.1088/1757-899x/147/1/012137 [2] Lo Russo nd ivit M V 2009 Open-loop groundwter het pumps development for lrge buildings: A cse study Geothermics 38 pp 335 345 [3] rbu I nd ebrchievici 2014 Generl review of ground-source het pump systems for heting nd cooling of buildings nergy nd Buildings 70 pp 441-454 [4] Ouzzne M, slmi-nejd P, Aidoun Z nd Lmrche L 2014 Anlysis of the convective het exchnge effect on the undisturbed ground temperture olr nergy 108 pp 340-347 [5] Dumitrscu G 2008 he wy to optimize the eversible cycles ermotehnic 2/2008 pp 18-22 [6] Dumitrscu G nd Horbniuc B 2010 Optimizre exergoeconomic Politehnium 10