ESTIMATION OF AUTOMOBILE COOLING LOADS FOR AIR CONDITIONING SYSTEM DESIGN

Similar documents
Cooling Load Calculation and Thermal Modeling for Vehicle by MATLAB

ME 476 Solar Energy UNIT THREE SOLAR RADIATION

CZECH TECHNICAL UNIVERSITY IN PRAGUE FACULTY OF MECHANICAL ENGINEERING DEPARTMENT OF ENVIRONMENTAL ENGINEERING

Chapter Seven. Solar Energy

ME 430 Fundamentals of Solar Energy Conversion for heating and Cooling Applications

SEM-2017(03HI MECHANICAL ENGINEERING. Paper II. Please read each of the following instructions carefully before attempting questions.

Time Series Model of Photovoltaic Generation for Distribution Planning Analysis. Jorge Valenzuela

Coolant. Circuits Chip

Draft for comment (613) Public Works Canada Design & Construction Technology DRAn' 1 July, Computer-Aided Design Centre (CAD) TSS COPY

(Refer Slide Time: 00:01:19 min)

Chapter 1: 20, 23, 35, 41, 68, 71, 76, 77, 80, 85, 90, 101, 103 and 104.

Orientation of Building

Greenhouse Steady State Energy Balance Model

Thermal behavior and Energetic Dispersals of the Human Body under Various Indoor Air Temperatures at 50% Relative Humidity

PROPOSAL OF SEVEN-DAY DESIGN WEATHER DATA FOR HVAC PEAK LOAD CALCULATION

ATMOSPHERIC ENERGY and GLOBAL TEMPERATURES. Physical Geography (Geog. 300) Prof. Hugh Howard American River College

Energy flows and modelling approaches

Heat and Mass Transfer Unit-1 Conduction

Interaction between building services and environment.

MDSRC Proceedings, December, 2017 Wah/Pakistan. Cooling load analysis by CLTD and RLF. Sher ali nawaz 1 Muzammil shahbaz 2

Double-Skin Facade in Low-Latitude: Study on the Absorptance, Reflectance, and Transmittance of Direct Solar Radiation

PERFORMANCE EVALUATION OF REFLECTIVE COATINGS ON ROOFTOP UNITS

Vertical Illuminance Measurement for Clear Skies in Tehran

Examination Heat Transfer

Autumn 2005 THERMODYNAMICS. Time: 3 Hours

SEASONAL AND DAILY TEMPERATURES

University of New Mexico Mechanical Engineering Spring 2012 PhD qualifying examination Heat Transfer

Chapter 2 Available Solar Radiation

SOLAR GEOMETRY (AND SOLAR RADIATION)

BACK TO BASICS: PIPE INSULATION

Chapter 1 Solar Radiation

The influence of solar radiation on the distribution of temperatures in historic masonry

DEVELOPMENT OF COOLING LOAD TEMPERATURE DIFFERENTIAL VALUES FOR BUILDING ENVELOPES IN THAILAND

R13 SET - 1 '' ''' '' ' '''' Code No RT21033

Level 7 Post Graduate Diploma in Engineering Heat and mass transfer

CONSTRUCTION AND CALIBRATION OF A LOCAL PYRANOMETER AND ITS USE IN THE MEASUREMENT OF INTENSITY OF SOLAR RADIATION

Why the Earth has seasons. Why the Earth has seasons 1/20/11

CAE 331/513 Building Science Fall 2016

TREES Training for Renovated Energy Efficient Social housing

MATHEMATICAL MODEL OF A VAPOUR ABSORPTION REFRIGERATION UNIT

Online homework #6 due on Tue March 24

Institut national des sciences appliquées de Strasbourg GENIE CLIMATIQUE ET ENERGETIQUE APPENDICES

Parametric Study of Rock Bed Thermal Regenerator for Space Heating

INVESTIGATION OF THE PERFORMANCE OF COOLING PANELS: CEILING AND FLOOR PANElS.

Thermodynamics Introduction and Basic Concepts

Analytical Design of Isolations for Cryogenic Tankers

Basic information about developed Calculator for solar / wind hybrid power supply

If there is convective heat transfer from outer surface to fluid maintained at T W.

Thermal Comfort; Operative Temperature in the Sun

Applied Thermal Engineering

Enerpool v Swimming Pool Simulation Simulation File: C:\ENERPOOL\SANDA.EPL Run on: Fri Feb 07 19:27:

The Numerical Psychrometric Analysis

Thermodynamics I Spring 1432/1433H (2011/2012H) Saturday, Wednesday 8:00am - 10:00am & Monday 8:00am - 9:00am MEP 261 Class ZA

CHAPTER 3. The sun and the seasons. Locating the position of the sun

Evaluation of monthly total global and diffuse solar radiation in Ibi, Taraba state, Nigeria

Winter Thermal Comfort in 19 th Century Traditional Buildings of the Town of Florina, in North-Western Greece

S.E. (Chemical Engineering) (Second Semester)EXAMINATION, 2012 THERMODYNAMICS-I (2008 PATTERN) Time : Three Hours Maximum Marks : 100

PROJECTING THE SOLAR RADIATION IN NASARAWA-NIGERIA USING REITVELD EQUATION

1. Basic state values of matter

Overall Heat Transfer Coefficient

Stochastic modeling of Extinction coefficients for solar power applications

Global Climate Change

SOLAR GEOMETRY (AND RADIATION)

Sunburned Products: Numerical Aging caused by Sun Exposition

Chapter 6. Solar Geometry. Contents

MAE 320 HW 7B. 1e. For an isolated system, please circle the parameter which will change with time. (a) Total energy;

Effect of moisture transfer on heat energy storage in simple layer walls

Principles of Energy Conversion Part 11A. Solar Energy Insolation

MC 405 MODEL TEST PAPER - 1 THERMAL SCIENCE & ENGINEERING. Time: Three Hours Maximum Marks: 100

EFFECT OF INTERNAL LONG WAVE RADIATION AND CONVECTION ON FENESTRATION SIMULATION

Solar radiation in Onitsha: A correlation with average temperature

Topic # 12 Natural Climate Processes

Common Terms, Definitions and Conversion Factors

True/False. Circle the correct answer. (1pt each, 7pts total) 3. Radiation doesn t occur in materials that are transparent such as gases.

Topic # 11 HOW CLIMATE WORKS continued (Part II) pp in Class Notes

Cooling of Electronics Lecture 2

Sunlight and its Properties II. EE 446/646 Y. Baghzouz

Calculating equation coefficients

Estimation of Hourly Solar Radiation on Horizontal and Inclined Surfaces in Western Himalayas

ApacheLoads Calculation Methodology as used in Revit MEP 2008 <Virtual Environment> 5.7

Radiation, Sensible Heat Flux and Evapotranspiration

AVAILABLE SOLAR RADIATION THEORETICAL BACKGROUND

UBMCC11 - THERMODYNAMICS. B.E (Marine Engineering) B 16 BASIC CONCEPTS AND FIRST LAW PART- A

Topic # 12 How Climate Works

A State Space Approach for the Dynamic Analysis of Automotive Air Conditioning System

Estimation of Hourly Global Solar Radiation for Composite Climate

Topic # 11 HOW CLIMATE WORKS PART II

S6. (a) State what is meant by an ideal gas...

Tick the box next to those resources for which the Sun is also the source of energy.

Div. 1 Div. 2 Div. 3 Div.4 8:30 am 9:30 pm 12:30 pm 3:30 pm Han Xu Ruan Pan

Eng Heat Transfer I 1

Average Monthly Solar Radiations At Various Places Of North East India

ANALYSIS OF GATE 2018*(Memory Based) Mechanical Engineering

UNIT TEST PRACTICE TEST

Outline. Stock Flow and temperature. Earth as a black body. Equation models for earth s temperature. Balancing earth s energy flows.

TankExampleNov2016. Table of contents. Layout

Section 3.5 Thermal Comfort and Heat Stress

FUNDAMENTALS OF HVAC

G109 Alternate Midterm Exam October, 2004 Instructor: Dr C.M. Brown

Chapter 3. Multiple Choice Questions

Transcription:

Nigerian Research Journal of Engineering and Environmental Sciences 596 Original Research Article ESTIMATION OF AUTOMOBILE COOLING LOADS FOR AIR CONDITIONING SYSTEM DESIGN *Omo-Oghogho, E., Aliu, S.A., Ebu-Nkamaodo, O.T. and Obanor, A. I. Department of Mechanical Engineering, Faculty of Engineering, University of Benin, Benin City, Nigeria *eghosa.omo-oghogho@uniben.edu ARTICLE INFORMATION Article history: Received 20 December 2017 Revised 25 December, 2017 Accepted 27 December, 2017 Available online 29 December, 2017 Keywords: Cooling loads Waste heat Automobile Air conditioning Heat transfer coefficient ABSTRACT Automobile air conditioning by vapour compression refrigeration system increases the amount of energy consumed. A vapour absorption refrigeration system which will utilize vehicle engine waste heat either from the exhaust gases or from the cooling system thereby reducing energy consumption is being proposed. Estimation of automobile cooling load was undertaken as part of an investigation to determine the suitability of a vapour absorption refrigeration system for automobile air conditioning system. Heat balance method was used to estimate the various thermal loads, radiation loads and cooling load of a car with four passengers in Lagos, Nigeria. The value of the cooling load obtained was 2.35kW and it was estimated under the hot weather condition of the sun to give a fairly accurate result for design purpose. 2017 RJEES. All rights reserved. 1. INTRODUCTION The first automobile air conditioning system introduced by Packer in 1939 used a vapour compression refrigeration system (Tiwari and Parishwad, 2012). In the last 70 years, automobile air conditioning systems have undergone gradual and continual improvements in performance and efficiency as a result of improvements in the individual components (McLaughlin, 2005). Alternative systems powered by heat which can be used in automobiles are solid adsorption cooling systems, thermo acoustic refrigeration, active magnetic regenerator system, thermo electric devices and vapour absorption refrigeration (VAR) systems (Zoontjens et al., 2005).

597 Air conditioning of a vehicle can be done by vapour compression refrigeration system (VCRS) or vapour absorption refrigeration system (VARS). Using a VARS will lead the automobile manufacturing industries into new technological ideas. For example, the compressor in a VCRS needs a pulley and a belt drive attached to the engine shaft before it can work, using a VARS eliminates that. A thorough cooling load calculation helps to minimize system over-sizing and the associated negative effects on thermal comfort and humidity control. The aim of this study was to properly estimate the various heat load of the air conditioning system of a vehicle. 2. MATERIALS AND METHOD 2.1. Materials A KIA PICANTO 2013 model with three cylinder engine was used for this study. The material properties used were in accordance with Fayazbakhsh and Bahrami (2013) and is shown in Table 1. Table 1: Material properties Property Glass Vehicle body Conductivity (W/mk) 1.05 0.2 Density (kg/m 3 ) 2500 1500 Transmisivity 0.5 0 Absorptivity 0.3 0.4 Specific heat (J/kgK) 840 1000 Thickness (mm) 3 10 2.2. Weather Data To effectively carry out this analysis weather data of the selected area of Lagos, Nigeria (Lagos State was choosing because the company that was used as a case study is situated at Lagos State, Kia picanto model was used because it was the model that was readily available and it carries only four persons), (Latitude: 6.5 o N, Longitude: 3.4 o E) were obtained. Dry bulb temperature (DBT) was obtained as 34 o C, while wet bulb temperature (WBT) was obtained as 26 o C. It is desired that the temperature be maintained at 24 o C DBT and relative humidity of 50%. 2.3. Heat Balance Method Heat balance method (HBM) was used to analyze the various thermal loads and radiation loads of the vehicle (Fayazbakhsh and Bahrami, 2013). The total thermal load on the vehicle cabin is given by Equation (1). The recommended heat-gain values for passengers and design specifications for the vehicle model are presented in Tables 2 and 3 respectively. = + + + + + (1) Where = Total load, = Ambient load, = Direct radiation load, = Engine load, = Exhaust load, = Reflected load and = Occupant load. 2.3.1. Heat gain due to vehicle occupants Heat gained due to vehicle occupants was determined for the driver and three passengers using Table 2. The sensible heat gain due to occupants (Q SCP) and latent heat gain due to occupants (Q LCP) were obtained from Equations (2) and (3) respectively. Q SCP = Q PS x Q AP (2)

598 Q LCP = Q AP - Q SCP (3) Where Q PS is the percentage sensible heat and Q AP is the adjusted male/female passenger load. Total heat gain due to occupants is given as: Q p = Q SCP + Q LCP (4) Total Sensible Heat is given as: Q SC = Q SCP + Q Dir + Q Ref + Q Amb + Q Exh + Q Eng (5) Sensible heat factor is given as: SHF = Table 2: Recommended heat-gain values for passengers and operators (ASHRAE, 2001) Application Adult male Adjusted % Sensible (w/person) male/female heat Inter-city and commuter car, seated 132 117 61 Inter-city and commuter car, standing 139 132 56 Subway and Light Rail Vehicle, Seated 139 132 56 Subway and Light Rail Vehicle, standing 161 147 50 Operators / drivers 132 132 56 Table 3: Kia Picanto 2013 model design specifications S/N Parameters Symbol Specification 1 Cabin surface area S C 3.8m 2 2 Surface area exposed to exhaust pipe S EXH 0.8774m 2 3 Surface area of cabin exposed to the engine S ENG 1.18m 2 4 Surface area of roof S R 0.963m 2 5 Side window area S SW 1.44m 2 6 Minimum speed of the vehicle V 5.56m 2 7 Petrol flow rate V P 0.12kg/s 8 Windscreen area S W 0.72m 2 2.3.2. Thermal and radiation loads i. Ambient load The ambient load which is the thermal load due to temperature difference between the ambient and cabin air is given by Equation (7). =! " #$ $ & (7) Where = Cabin surface area, $ = Temperature of surface and $ = Temperature of cabin ii Exhaust load Exhaust load is the load due to the high temperature of exhaust gases. The exhaust load was calculated using Equation (8). =! ' #$ $ & (8) = Surface area of exhaust, $ = Temperature of exhaust iii. Engine load This is the load due to the high temperature of the engine. It was calculated from Equation (9). =! ( )$ $ * (9) Where: (6)

599! " = " + (10), " = " - +. / + " 0 (11), " = Net thermal resistance for a unit surface area, h 2 = Outside convection coefficient. h = inside convection coefficient 3 = Surface thermal conductivity. 4= Surface element thickness. According to Fayazbakhsh and Bahrami (2013), the outside convection coefficient is given as: h 2 = 0.6 + 6.64 6 (12) where V= is the minimum vehicle speed $ = 0.138,78 17 (13) $ = -2 10-6,78 + 0.0355,78 + 77.5 = 1500 (14) According to Ding and Zito #2001&,h = 8.3 W/m 2 K! ' = " F,, ' =. / + " 0,! ( = " G,, ( =. / + " 0 iv Radiation load Radiation loads are obtained for both the glass surfaces and other surfaces of the car. Equations 15-28 are the commonly used models for air conditioning calculations as suggested by ASHRAE (2001). Direct radiation load = LM H I H J K cosp (15) J K = Q R (16) STUV W A is the apparent solar irradiation taken as 1230W/m 2 for months of January to December and 1080W/m 2 for mid-summer while B is the atmospheric extinction coefficient which has a values of 0.142 for the months of January to December. The declination, D is given as: X=23.47]^_ (`a #'bcdk& (`e Where day of the year starting from January, N= 21 for January 21, 2015 at Lagos, Nigeria Local Solar Time, LST = 14.00 (2.00pm) Hour angle, h = 15(LST - 12) (18) Tilt angle, f = 0 (horizontal surface) and 90 o (vertical surface), Wall azimuth, ε= 0 o (southern surface), ε= 90 o (western surface) and ε= 180 o (northern surface) The altitude angle, is given as: β= sin g" #cos#h&ij]#h&cos#k&+sin#h&sin#k&& (19) L = latitude, h = hour angle, d = angle of declination, I= roof element transmissivity Since surface element transmissivity of roof, wall and floor is zero, the direct radiation through those elements becomes zero. The incident angle is given as: P=ij] g" (sin#l&cos#f&+cos#l&cos#m&sin#f&& (20) The azimuth angle is given as: α = [n #o+f&] q (21) The solar azimuth is given as: (17)

600 γ = sin g" Q rst#u& tvw#& W (22) rst#x& F = 1 for afternoon, where n=180 Diffuse radiation load (Q DIF) = LM L{ I J (23) J = C J q ] #}/ ' ) (24) Where J = diffuse radiation and C, a constant = 0.058 for January. The view factor is given as: Fws = "d2 (25) ' Ɛ = tilt angle and its value will be 0 for horizontal surfaces and 90 for vertical surfaces. Reflected radiation loads = LM L{ I J (26) J = (J + J ) P G FwG (27) Ground reflectivity factor, P G = 0.2 Angle factor is given as: (FwG) = "g2 (28) ' FwG = 0, for horizontal surfaces and FwG = 0.5, for vertical surfaces. 3. RESULTS AND DISCUSSION The values obtained for different types of load which make the total cooling load of the vehicle are presented in Table 4. The total heat load, Q T was obtained as 2352.21W and the sensible heat factor, an important parameter for determining cooling coil capacity was 0.92. The various heat loads were gotten using various parameters which were gotten under the hot weather condition of the sun to get a fairly accurate heat load estimation so as to improve on the design and working conditions of the vehicle air conditioning system. Table 4: Total cooling load Q SCP 288.O3 Q LCP 194.97 Q AMB 163.78 Q EXH 705 83 Q ENG 852.04 Q DIR 115.2 Q DIF 7.11 Q REF 25.25 Q T 2352.21W 4. CONCLUSION The total heat load of the vehicle was estimated to be 2352.21W and as such the power to efficiently design and run the air conditioning system of the vehicle should be higher than the estimated heat load of the vehicle. This will help in proper and efficient design of the car air conditioning system.

601 5. CONFLICT OF INTEREST There is no conflict of interest associated with this work. REFERENCES ASHRAE (2001). Handbook of Fundamentals, American Society of Heating, Refrigerating and Air Conditioning, Atlanta, GA. Ding, Y. and Zito, R. (2001). Cabin heat transfer and air conditioning capacity (No. 2001-01-0284). SAE Technical Paper. Fayazbakhsh, M. A. and Bahrami, M. (2013). Comprehensive modeling of vehicle air conditioning loads using heat balance method (No. 2013-01-1507). SAE Technical Paper. McLaughlin, S. (2005). An Alternative Refrigeration System for Automotive Applications. M.S. Thesis, Mechanical Engineering Department, Mississippi State University. Tiwari, H. and Parishwad, G. V. (2012). Adsorption refrigeration system for cabin cooling of trucks. International journal of emerging technology and advanced engineering, 2(10), pp. 337-342. Zoontjens, L., Howard, C., Zander, A. and Cazzolato, B. (2005). Feasibility Study of an Automotive Thermo acoustics Refrigerator. In: Proceedings of ACOUSTICS 9-11 November 2005, Busselton, Western Australia.