Journal of Chemical and Pharmaceutical Research, 2013, 5(12): Research Article

Similar documents
The Open Petroleum Engineering Journal

Simulation and verification of a plate heat exchanger with a built-in tap water accumulator

RP 2.4. SEG/Houston 2005 Annual Meeting 1513

Carnot Factor of a Vapour Power Cycle with Regenerative Extraction

The Basics of Vacuum Technology

De-Coupler Design for an Interacting Tanks System

Study of Convective Heat Transfer through Micro Channels with Different Configurations

Lecture 38: Vapor-compression refrigeration systems

DYNAMIC MODELING OF ORGANIC RANKINE CYCLE (ORC) SYSTEM FOR FAULT DIAGNOSIS AND CONTROL SYSTEM DESIGN

Unsteady State Simulation of Vapor Compression Heat Pump Systems With Modular Analysis (This paper will not be presented.)

The entransy dissipation minimization principle under given heat duty and heat transfer area conditions

Chapter 3 Thermoelectric Coolers

Research on the Negative Permittivity Effect of the Thin Wires Array in Left-Handed Material by Transmission Line Theory

Chapters 19 & 20 Heat and the First Law of Thermodynamics

Lecture 10: Carnot theorem

Section A 01. (12 M) (s 2 s 3 ) = 313 s 2 = s 1, h 3 = h 4 (s 1 s 3 ) = kj/kgk. = kj/kgk. 313 (s 3 s 4f ) = ln

SENSITIVITY ANALYSIS FOR COUNTER FLOW COOLING TOWER- PART I, EXIT COLD WATER TEMPERATURE

A = h w (1) Error Analysis Physics 141

A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225

3.4 Worksheet: Proof of the Chain Rule NAME

WINKLER PLATES BY THE BOUNDARY KNOT METHOD

Analysis of Static and Dynamic Load on Hydrostatic Bearing with Variable Viscosity and Pressure

Distribution of reynolds shear stress in steady and unsteady flows

How to Find the Derivative of a Function: Calculus 1

Fundamentals of Heat Transfer Muhammad Rashid Usman

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households

. If lim. x 2 x 1. f(x+h) f(x)

CFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow Neale, A.; Derome, D.; Blocken, B.; Carmeliet, J.E.

Explicit Interleavers for a Repeat Accumulate Accumulate (RAA) code construction

Problem Solving. Problem Solving Process

Parshall Flume Discharge Relation under Free Flow Condition

Theoretical Analysis of Flow Characteristics and Bearing Load for Mass-produced External Gear Pump

lecture 26: Richardson extrapolation

Empirical models for estimating liquefaction-induced lateral spread displacement

Poisson Equation in Sobolev Spaces

Experimental Analysis of Heat Transfer Augmentation in Double Pipe Heat Exchanger using Tangential Entry of Fluid

Chapter 4 Optimal Design

International Power, Electronics and Materials Engineering Conference (IPEMEC 2015)

Velocity distribution in non-uniform/unsteady flows and the validity of log law

Effects of Radiation on Unsteady Couette Flow between Two Vertical Parallel Plates with Ramped Wall Temperature

The total error in numerical differentiation

Simulation and Analysis of Biogas operated Double Effect GAX Absorption Refrigeration System

MIMO decorrelation for visible light communication based on angle optimization

Sedimentation. Treatment Processes. Sedimentation. Sedimentation. Sedimentation. Sedimentation. CIVL 1112 Water Treatment - Sedimentation 1/7

CFD Analysis and Optimization of Heat Transfer in Double Pipe Heat Exchanger with Helical-Tap Inserts at Annulus of Inner Pipe

8-4 P 2. = 12 kw. AIR T = const. Therefore, Q &

A Reconsideration of Matter Waves

Computational Method of Structural Reliability Based on Integration Algorithms

2.3 Product and Quotient Rules

Precalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!

Chapter A 9.0-V battery is connected to a lightbulb, as shown below. 9.0-V Battery. a. How much power is delivered to the lightbulb?

Journal of Engineering Science and Technology Review 7 (4) (2014) 40-45

Desalination by vacuum membrane distillation: sensitivity analysis

CALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY

Comment on Experimental observations of saltwater up-coning

LEARNING FROM MISTAKES

CFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow. Annex 41 Kyoto, April 3 rd to 5 th, 2006

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.

PHYSICAL PROCESSES IN ANISOTROPIC THERMOELEMENT AND THEIR FEATURES

2016 PRELIM 2 PAPER 2 MARK SCHEME

Fast Computation of Capacitance Matrix and Potential Distribution for Multiconductor in Non-Homogenous Multilayered Dielectric Media

MVT and Rolle s Theorem

Analysis: The speed of the proton is much less than light speed, so we can use the

2.11 That s So Derivative

Effects of heat leak on the performance characteristics of Carnot like heat engines and heat pumps

Polynomial Interpolation

HOMEWORK HELP 2 FOR MATH 151

Determination of heat transfer intensity between free streaming water film and rigid surface using thermography

Non-linear Dynamics of Inlet Film Thickness during Unsteady Rolling Process

AN ANALYSIS OF AMPLITUDE AND PERIOD OF ALTERNATING ICE LOADS ON CONICAL STRUCTURES

ch (for some fixed positive number c) reaching c

Use of fin analysis for determination of thermal conductivity of material

The Laws of Thermodynamics

Chapter 5 FINITE DIFFERENCE METHOD (FDM)

Compressor 1. Evaporator. Condenser. Expansion valve. CHE 323, October 8, Chemical Engineering Thermodynamics. Tutorial problem 5.

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative

DEVELOPMENT AND PERFORMANCE EVALUATION OF THERMAL CONDUCTIVITY EQUIPMENT FOR LABORATORY USES

232 Calculus and Structures

Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 01 Aug 08, 2016.

Elmahdy, A.H.; Haddad, K. NRCC-43378

The Verlet Algorithm for Molecular Dynamics Simulations

Polynomial Interpolation

DEBONDING FAILURES OF RC BEAMS STRENGTHENED WITH EXTERNALLY BONDED FRP REINFORCEMENT: BEHAVIOUR AND MODELLING

Wind Turbine Micrositing: Comparison of Finite Difference Method and Computational Fluid Dynamics

COMPARISON OF FUZZY LOGIC CONTROLLERS FOR A MULTIVARIABLE PROCESS

Numerical Differentiation

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point

Performance Prediction of Commercial Thermoelectric Cooler. Modules using the Effective Material Properties

ADCP MEASUREMENTS OF VERTICAL FLOW STRUCTURE AND COEFFICIENTS OF FLOAT IN FLOOD FLOWS

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1

Bandwidth Selection in Nonparametric Kernel Testing

Taylor Series and the Mean Value Theorem of Derivatives

Given: Hot fluid oil, Cold fluid - water (T 1, T 2 ) (t 1, t 2 ) Water

Dedicated to the 70th birthday of Professor Lin Qun

f a h f a h h lim lim

Optimal Shape Design of a Two-dimensional Asymmetric Diffsuer in Turbulent Flow

A finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations

Transcription:

Available online.jocpr.com Journal of emical and Parmaceutical Researc, 013, 5(1):55-531 Researc Article ISSN : 0975-7384 ODEN(USA) : JPR5 Performance and empirical models of a eat pump ater eater system it eat recovery of aste ater Handong Wang Scool of Mecanical & Electrical Engineering, Senzen Polytecnic, Senzen, ina ABSTRAT In order to save energy, te autor developed a type of air-asteater source eat pump (AWSHP) ater eater system. Wat presented ere is about te experimental performance and empirical models of te AWSHP system. Experiments and analysis soed tat en te temperatures of aste ater and city ater varied in te ranges of 16.9~35.1 and 18.3~3.6, respectively, te AWSHP system could produce 41.~53.9 ot ater it flo rate of 3.53 10-5 ~9.31 10-5 m 3 /s. Te coefficient of eating performance (OP ) varied from 3.89 to 5.35. On te basis of energy balance, syntetic degree of grey incidence (SDGI) and grey-box metods, it as found tat te OP and ot ater flo rate depended on te ratio of inlet and outlet temperature differences of aste ater and ot ater, i.e. /. Te input poer of compressor as proportional to te outlet temperature of ot ater. Tus, te empirical correlation equations ere obtained. Keyords: eat recovery; eat pump; aste ater; empirical model; energy saving INTRODUTION In families and oter fields suc as artificial arm springs, sauna rooms and HVA engineering, ot ater is generally required above 40. But in tese fields, especially in families and sauna rooms, te aste ot ater generated during soer bat or sauna is directly drained into te environment. Tis causes bot energy-asting and termal pollution to te environment. Measurements soed tat te temperature of aste ater generated by soer bat as in te range of 30~36, ic is still valuable for utilizing. In order to recover te eat of soer aste ater and provide local air-conditioning, on te basis of researc results of [1 4], te autor developed a type of air-ater source eat pump (AWSHP). It uses air and aste ater, especially soer aste ater, as eat sources. Te AWSHP could produce enoug ot ater for bating by recovering te eat from air and soer aste ater, and soed satisfactory performance of energy saving. Wat presented in tis paper is about te experimental performance and te data-based empirical models of te AWSHP ater eater system, ic is a furter study on te performance of AWSHP presented by [5]. FLOW HART AND EXPERIMENTAL SET OF THE AWSHP Te basic flo cart diagram of te AWSHP ater eater is son in Figure 1. It can be installed in batroom or elseere. Tere is no need of ot ater tank and tis makes it more compact. Its basic principle is te same as common ater-source eat pump ater eater. In tis type of AWSHP, te air source evaporator (AE) and ater source evaporator (WE) are arranged in serial instead of parallel to maintain te different evaporating temperatures. 55

Handong Wang J. em. Parm. Res., 013, 5(1):55-531 Te AWSHP can efficiently recover eat from aste ater and air to save energy and avoid te eat pollution caused by directly drained ot ater. In te experimental set, te outlet temperature of ot ater is regulated by a condenser-pressure-dependent valve ic also controls te condenser pressure and temperature by regulating cooling ater floing troug te condenser. Te trottling part of te AWSHP is a set of capillary tubes. Te air source evaporator is a self-made eat excanger it a fan. A plate-type eat excanger and a sell-tube eat excanger are used as te condenser and aste ater evaporator, respectively. Te refrigerant is R and te input poer of aste ater pump is 5W. Figure 1. Diagram of flo cart and energy-excanging of AWSHP ater eater. omp compressor; ond condenser; ap capillary; AE air source evaporator; WE ater evaporator; WWHE ater-ater eat excanger; WWL aste ater collector. A data acquisition system based on LabVie virtual instrument tecnology is equipped to te experimental set. Te temperature sensors are T-type termocouples. Te flo rates of ot ater and aste ater are measured by intelligent electromagnetic flo meters it accuracy of ±1.0%. Te pressures are manually measured by refrigerant pressure gauges. Te input poer of te eat pump system is measured on line by a Fluke 43B poer quality analyzer. Te umidity of air floing troug te air source evaporator is measured by digital temperature-umidity meters it accuracy of ±3.0% (for relative umidity), and te velocity of air flo in te air source evaporator is measured by an anemometer it accuracy of ±.0%. ALULATION BASIS In te AWSHP ater eater, driven by input poer of compressor, te refrigerant only excanges energy it air, city ater and aste ater, as son in Figure 1. From te vie of energy balance, te performance of te AWSHP sould depend on input poer of compressor (N) and eat capacities of condenser ( ), air source evaporator ( AE ) and aste ater source evaporator ( WE ). Te coefficient of eating performance of te AWSHP is defined by Eq.1: OP = (1) N ere, N can be measured by te Fluke 43B poer quality analyzer. is also te eat transferred to ot ater. If te eat loss of condenser is neglected, can be calculated from te parameters of ot ater as Eq.: cv = ρ ( T, o T, m 3600 ) () 56

Handong Wang J. em. Parm. Res., 013, 5(1):55-531 ere T,m is te temperature of ot ater leaving off te WWHE,. T,m is determined by te eat excanging efficiency of WWHE and te initial temperature of city ater named as T,in (as son in Figure.1). Based on te energy balance of te system, if te eat loss of te system is neglected, tere is Eq.3: = N + + (3) AE WE Were, WE cv = ρ ( T, m T, o 3600 ) (4) AE = ρ V ) (5) A A ( A, in A, o In Eqs.4 and 5, T,m is te temperature of aste ater leaving off te WWHE, (as son in Figure 1). T,m is determined by te eat excanging efficiency of WWHE and te initial temperature of aste ater named as T,in. V A is te volume flo rate of air floing troug te AE, m 3 /s. V A =v F, v and F are te velocity of air flo (m/s) and area (m ) of te air supplying section of te air source evaporator, respectively. Eq.3 also can be used to judge te eat loss of te system. PERFORMANE AND EMPIRIAL MODELS OF THE AWSHP Tere are many teoretical models based on matematic equations of compressor, condenser, trottle valve (or capillary tube) and evaporator to simulate te performance of eat pump system, suc as [, 6 9]. But all of te teoretical models need many parameters and sopisticated computer calculation. Tey are suitable for simulating and designing te system but not convenient for engineering evaluation on te fields. Yu (1995) [10] establised an empirical OP model solely depending on te temperature difference of condensing and evaporating temperature, but it is found not suitable for evaluating te performance of te AWSHP son in Figure 1. In order to investigate te performance and establis te empirical models of te AWSHP, experiments are carried out at various conditions suc as different inlet temperatures of aste ater and city ater, it or itout aterater eat excanger, different outlet temperature of ot ater, different volume of carged refrigerant (R), and so on. Experimental performance Experiments and analysis so en te temperatures of aste ater and city ater are in te range of 16.9~35.1 and 18.3~3.6, respectively, te AWSHP can produce ot ater it temperature of 41.~53.9, and te volume flo rate of ot ater varies in te range of 0.13~0.34 m 3 / (3.53 10-5 ~ 9.31 10-5 m 3 /s). Te OP varies in te range of 3.89~5.35. Wen te input poer of pump and fan are taken into account, te overall energy efficiency ratio (EER) varies in 3.74~5.10. It is also found tat te WWHE deteriorates te OP, but it can efficiently increase te ot ater flo rate. Terefore, te WWHE is necessary to provide enoug ot ater en te city ater temperature is loer suc as in inter. Experiments so te temperature and relative umidity of air at outlet of te AE vary in te range of 1.3~0.8 and 83.1~93.4%, respectively. But te refrigeration capacity of te AE is just 5.~13.% of tat produced by te WE. Terefore, te influence of te AE on te performance of AWSHP can be neglected in te present analysis. It also indicates tat te AE sould be redesigned to supply enoug cold air for local air-conditioning. Empirical models Te empirical models are based on energy-excanging and te affecting factors of te AWSHP. It can be deduced From Eqs.1 3 tat te OP of te AWSHP is te function of N, AE and WE. As AE and WE depend on te temperatures and flo rates of air flo and aste ater flo, te OP can be expressed by te f 1 function of te parameters as Eq.6: 57

Handong Wang J. em. Parm. Res., 013, 5(1):55-531 OP =f 1 (N, A,in, A,o, V A, T,in, T,o, V ) (6) Because te experiments so tat te AE is very small and it takes little effect on OP, te parameters of te AE can be omitted by adjusting oter coefficients, e.g. T,in. It is also found tat te N is closely related to te outlet temperature of ot ater (T,o ). Tus, en te aste ater flo rate V is kept constant, Eq.6 can be simplified as Eq.7: OP =f (T,o, T,in, T,in, T,o ) (7) Furter analysis on experimental data sos tat OP of te AWSHP closely depends on te ratio of temperature difference / en V is constant. / can be calculated as Eq.8: T = T, in, o T T, o, in (8) Terefore, Eq.7 can be furter simplified as Eq.9: OP =f 3 ( / ) (9) Similar to OP, te influence factors of te performances of te AWSHP can be summarized as T,in, T,o, T,in, T,o,, and /. Based on te analysis metod of syntetic degree of grey incidence (SDGI) proposed by [11], te most important factors of te performances of te AWSHP can be determined as / and T,o. Tus, te empirical models can be proposed by coosing / or T,o as te variables. If te grey-box metod is used to determine te empirical models, en te aste ater flo rate is kept constant, te empirical regression equations of OP, V and N can be establised based on te experiment data and principle of te least square metod as Eqs.10 1: OP = A T + B + (10) V = D E (11) + N = m T, n (1) o + ere A, B,, D, E, m and n are constant coefficients determined by regression data. Te experiment data of OP, V and N are son in Figures 4. Te corresponding empirical regression equations of OP, V and N in Figures 4 are Eqs.13 15: OP.518 = + 5.704 + 3.1159 (13) = 0.454 + 0.0551 V (14) N 0.014, + 0.03 (15) = T o 58

Handong Wang J. em. Parm. Res., 013, 5(1):55-531 6 5 4 OP 3 1 0 0 0.1 0. 0.3 0.4 0.5 0.6 0.7 / Figure. Relationsip of OP and / 0.4 0.35 0.3 V (m 3 /) 0.5 0. 0.15 0.1 0.05 0 0 0.1 0. 0.3 0.4 0.5 0.6 0.7 / Figure 3. Relationsip of V and / N (kw) 1.6 1.4 1. 1 0.8 0.6 0.4 0. 0 40 45 50 55 T,o ( ) Figure 4. Relationsip of N and T,o Validation and discussion of empirical models In order to investigate te validity of models expressed by Eqs.10 1, experiments ere carried out on te AWSHP it reduced carged refrigerant (R). Te results soed te perfect accordance of Eqs.10 1 and only te constant coefficients ere different. For example, te OP at tis condition is son in Figure 5, and te corresponding empirical regression equation of OP is expressed by Eqs.16. 59

Handong Wang J. em. Parm. Res., 013, 5(1):55-531 OP 13.735 = + 11.895 +.3817 (16) It demonstrates te validity of empirical models expressed by Eqs.10 1. As te temperatures of ot ater and aste ater can be easily measured, tese models can be conveniently used to evaluate and modulate te performances of ater-source eat pump ater eaters on te field. It also indicates tat te OP as a maximum or optimal value at certain /. Tis provides an engineering guide to improve te eating performance of te AWSHP. Tus, designers sould keep in mind to set te / of AWSHP in close proximity to its optimal value to obtain te maximum OP. OP 5 4.9 4.8 4.7 4.6 4.5 4.4 4.3 0. 0.3 0.4 0.5 0.6 / Figure 5. Relationsip of OP and / in te AWSHP it reduced carged R. Above analysis sos tat te performance of te AWSHP ater eater can be determined using te inlet and outlet temperatures of asteater and ot ater, especially te ratio of temperature differences of asteater and ot ater (i.e., / ). Te variable of / makes te above empirical models different to te oters presented in literatures in most of ic te absolute, not te relative temperatures of fluids ere selected as model variables. Meanile, te empirical models presented in tis paper not only can be used to predict te OP, but also can be used to evaluate te ot ater flo rate and te compressor poer-consumption of te AWSHP ater eater. Tus, it te above empirical models, te performance of AWSHP ater eater can be investigated from different points of vie, but not only from te OP. Te above empirical models also provide an easy ay for engineers to analyze and evaluate te performance of te AWSHP ater eater because measuring te temperatures of asteater and ot ater is easier tan measuring refrigerant parameters. Te empirical models are also elpful for evaluating and optimizing operational conditions. Toug Eqs.10 1 are obtained from experiment data of te AWSHP, tey are also suitable for evaluating oter ater-source eat pump ater eaters because te influence of air-source evaporator as neglected. It sould be mentioned tat te constant coefficients in Eqs.10 1 depend on operating conditions, suc as carged volume of refrigerant, dimensions of capillary tube, it or itout ater-ater eat excanger, aste ater flo rate and so on. ONLUSION Experiments and analysis indicate tat te AWSHP ater eater as satisfactory energy saving performance. Tis system also provides a simple and efficient ay to obtain sustainable use of aste ater by recovering te lograde eat from it. Te main points of tis paper can be summarized as folloing: 1) Te coefficient of eating performance (OP ) varies in te range of 3.89~5.35. Toug te WWHE deteriorates te OP of AWSHP, it can efficiently increase te ot ater flo rate. But te AE in te system sould be redesigned to provide enoug cooling capacity. 530

Handong Wang J. em. Parm. Res., 013, 5(1):55-531 ) Based on te analysis of energy balance, SDGI and grey-box, empirical models are obtained and te models so tat te OP, V and N only depend on te temperatures of ot ater and aster ater, especially te ratio of temperature differences of asteater and ot ater (i.e., / ). 3) Te empirical models presented in tis paper can be used to predict te OP, te ot ater flo rate and te compressor poer-consumption of te AWSHP ater eater at te same time. It can be conveniently used to evaluate and modulate te performances of te AWSHP or oter ater-source eat pump ater eaters on te field. Furter researc on te performances of te AWSHP it improved air-source evaporator ill soon be carried out. Nomenclature c specific eat of ater, kj/(kg ) f 1, f, f 3 signs of functions specific entalpy, kj/(kg ) N input poer of compressor, kw eat capacity, kw T temperature, V volume flo rate, m 3 / Greek symbols temperature difference, ρ density, kg/m 3 Superscripts and subscripts A air A,in inlet air A,o outlet air c condenser ot ater ot ater,in inlet of ot ater,m middle status of ot ater,o outlet of ot ater ater WE aste source evaporator asteater,in inlet of asteater,m middle status of asteater,o outlet of asteater Acknoledgement Te autor ill be grateful to Senzen Baomei Ne-Energy Resource Ltd. ompany for supporting tis researc. REFERENES [1] GS en; F Li; Journal of Guangdong University of Tecnology, 00, 19(4), 47 49. [] GX Kou; H Wang; XL Wang; XL Li; HV&A, 004, 34(6), 43 45. [3] HP Tan; Y en; Fluid Macinery, 008, 36(6), 83 85. [4] HP Tan; Fluid Macinery, 009, 37(7), 65 68. [5] HD Wang; MAE011 Proceedings, IEEE Press, 011, 7, 5803 5806. [6] J Ji; G Pei; W He; J Dong; WP Zao; HV&A, 003, 33(), 19 3. [7] H Tian; YT Ma; Z Wen; Journal of Termal Science and Tecnology, 008, 7(4), 373 378. [8] GY Xu; SH Li; XS Zang; WB Wu; Journal of Harbin Institute of Tecnology, 009, 41(), 11 13. [9] J Li; S.M Jin; Fluid Macinery, 010, 38(8), 73 79. [10] L Yu; HV&A, 1995, 5(1), 1 14. [11] SF Liu; YG Dang; ZG Fang; Grey system teory and application, 3 rd Edition, Science Press, Beijing, 004. 531