1999_The Vehicle Engine Cooling System Simulation Part 2 - Model Validation Using Transient Data
1999_The Vehicle Engine Cooling System Simulation Part 2 - Model Validation Using Transient Data
SAE TECHNICAL
PAPER SERIES 1999-01-0241
400 Commonwealth Drive, Warrendale, PA 15096-0001 U.S.A. Tel: (724) 776-4841 Fax: (724) 776-5760
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
The appearance of this ISSN code at the bottom of this page indicates SAE’s consent that copies of the
paper may be made for personal or internal use of specific clients. This consent is given on the condition,
however, that the copier pay a $7.00 per article copy fee through the Copyright Clearance Center, Inc.
Operations Center, 222 Rosewood Drive, Danvers, MA 01923 for copying beyond that permitted by Sec-
tions 107 or 108 of the U.S. Copyright Law. This consent does not extend to other kinds of copying such as
copying for general distribution, for advertising or promotional purposes, for creating new collective works,
or for resale.
SAE routinely stocks printed papers for a period of three years following date of publication. Direct your
orders to SAE Customer Sales and Satisfaction Department.
Quantity reprint rates can be obtained from the Customer Sales and Satisfaction Department.
To request permission to reprint a technical paper or permission to use copyrighted SAE publications in
other works, contact the SAE Publications Group.
No part of this publication may be reproduced in any form, in an electronic retrieval system or otherwise, without the prior written
permission of the publisher.
ISSN 0148-7191
Copyright 1999 Society of Automotive Engineers, Inc.
Positions and opinions advanced in this paper are those of the author(s) and not necessarily those of SAE. The author is solely
responsible for the content of the paper. A process is available by which discussions will be printed with the paper if it is published in
SAE Transactions. For permission to publish this paper in full or in part, contact the SAE Publications Group.
Persons wishing to submit papers to be considered for presentation or publication through SAE should send the manuscript or a 300
word abstract of a proposed manuscript to: Secretary, Engineering Meetings Board, SAE.
Printed in USA
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
1999-01-0241
1
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
Charge air cooler 01-23132-00 Allied Signal TEST MATRIX FOR VALIDATION
Fan + Fan drive DST variable speed clutch Kysor of Cadillac Table 2 gives the details of the tests used for validation. A
set of three tests at engine speeds of 1500 rpm, and
Coolant thermostat - Detroit Diesel Corp.
1800 rpm at different percent engine loads, and ambient
Oil Cooler - Harrison Division-GMC temperatures were used for stationary testing. Test no. 1
and 2 were with no airflow from the test cell boost fan
Oil thermostat - Detroit Diesel Corp. while test no. 3 was simulated a 25 m.p.h. RAM air veloc-
ity. For transient testing, three tests at 55%, 75% and
Turbocharger - Garrett Automotive
100% of full load on the vehicle were used over three dif-
ferent vehicle routes.
2
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
STATIONARY EXPERIMENTAL AND VECSS when the oil sump temperature first reaches 104°C (220
PREDICTED RESULTS °F). The oil thermostat and the oil sump temperature take
about 1800 seconds to stabilize. The predicted airflow
For all the stationary tests, the vehicle was stationary. A across the radiator depends on the fan speed which in
boost fan was used in the test cell to pressurize air over turn is a function of the engine speed and the radiator top
the vehicle to simulate airflow to the engine compartment tank temperature. As the top tank coolant temperature
from either wind speed or the vehicle speed. Load was reaches 88°C (190°F), the fan speed starts increasing
applied by means of a dynamometer. The coolant used and later reaches a steady speed of 1140 rpm. Table 4
was water. The air-conditioning was off and hence there gives the summary comparison of the VECSS predicted
was no energy release from the condenser. The cab parameters with the experimentally measured values
heater was also non-functional in all the stationary tests. when the system is at steady state condition. All the sim-
Initially, all the temperatures in the VECSS program were ulation runs were conducted on the PC with a 300 MHz
at ambient temperature. The “warming up” of the VECSS Pentium-II processor.
was done by using the DDEC-ECM measured engine
speed and fuel flow rate as input data to accurately repre- Table 2. Test Matrix for VECSS Validation
sent this transient operation of the engine / cooling sys-
tem from ambient temperature to steady state TEST TYPE No. TEST DESCRIPTION
In Figure 1 the first two time traces show the time traces Tests run with: 4. Cadillac to Spring Lake, Ambient: 18.9 deg.C (66
deg.F)
for the input data to the VECSS program. All the other * 50/50 Ethylene gly- 13609 kg (30,000 lbs) (55% of full load on vehicle)
time traces in that figure give the VECSS predicted col & water
5. Spring Lake to Plymouth, Ambient:1 deg.C ( 34
results. * Kysor DST variable deg.F)
speed clutch 18712 kg (41,250 lbs) (75% of full load on vehicle)
During the first 750 seconds, the fuel flow rate slowly
* Without shutters 6. Brighton to Spring Lake, Ambient: 11 deg.C (52
increases. This corresponds to the warm-up period for deg.F)
the engine test and the VECSS program. After 750 sec- 24949 kg (55,000 lbs) (100% of full load on vehi-
cle)
onds, the engine speed, and the fuel flow rate reach a
steady state value. The predicted engine power follows *% full load was calculated as the ratio of measured fuel flow to fuel flow from the
full load power curve for this engine
more or less the fuel flow rate curve. This predicted
engine power is the power obtained right at the engine
output shaft into the transmission and includes the power
used for the coolant pump, oil pump, and the fan. A repre- Table 3. Steady state input parameters for VECSS
sentative parameter for the engine operating conditions, program
i.e., the piston surface temperature, stabilizes after 1200
seconds. Time traces for intake air temperature and Input parameter Test No.1 Test No.2 Test No.3
mass flow rate are also shown in Figure 1. The top tank Engine speed (rpm) 1500 1800 1800
coolant temperature rises from the initial ambient temper-
Fuel flow rate (kg/min) 0.590 0.726 0.671
ature to a steady state value of 93°C (198°F). The time
trace for the percent open coolant thermostat shows that % Full load 60 80 75
it is first activated when the top tank coolant temperature
RAM air velocity (mph) 0 0 25
reaches 88°C (190°F). This results in flow of the coolant
through the radiator and hence a drop in the coolant tem- Average ambient temp. (deg.C) 23 21 17
perature. The ripples in top tank temperature time trace Total simulation real time (sec) 2800 2800 3750
are due to the opening and closing of the thermostat until
it reaches a steady state opening position. The oil ther-
mostat time trace involved modeling it with a 5% leak-
age, starting at time 0. The thermostat begins to activate
3
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
Figure 1. Time traces of select parameters for stationary Figure 2. Time traces of select parameters for stationary
Test No. 1 Test No. 1 (Continued)
4
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
As can be noted from Table 4 for the stationary tests, the TRANSIENT EXPERIMENTAL AND VECSS
oil sump temperature was predicted 0.5 to 1°C (1 to 2°F) PREDICTED RESULTS
lower than the experimental values. This can be attrib-
uted to the following reason. The oil sump temperature Both of the data systems were setup on the truck and the
was measured using only one thermocouple to represent truck was run over various routes in Michigan (refer Table
the entire oil sump, and might not be representative of 2). The tests were initiated in the starting cities and con-
the mass average temperature calculated by the VECSS ducted for the specified amount of time towards the desti-
since it could be close to an oil return from any of the nation city. Load on the truck was varied by changing the
auxiliaries (piston cooling, bearing, etc.). Furthermore, weight of load it carried. The cab-heater was off, just as
the differences between the actual and the predicted are in the case of the stationary tests. The coolant used was
only 0.5 to 1°C (1 to 2°F) low, which is well within the 50% water and 50% ethylene glycol (The VECSS pro-
experimental accuracy of the thermocouple used to mea- gram has been modeled to simulate either 50% water /
sure the oil sump temperature. 50% ethylene glycol or 100% water as the coolant [3]).
The predicted radiator top tank coolant temperature for Table 5 gives the input parameters for the transient tests.
each of the three tests agreed with the actual data to Most of these parameters do not remain constant with
within 0.5°C (1°F), which was within the experimental time, in contrast to the steady state tests, so either a
accuracy of the thermocouples used to measure the range or a time averaged value is given in this table. The
data. The predicted intake air temperature was within 8°C time averaged value is the average of about 2000 contin-
(14°F) and the power was within 4% of the experimen- uous data points. Figures 3,4 show sample time traces
tally recorded data, which also were within the accuracy for test no. 4 (Cadillac to Spring Lake) which are the time
of the experimental measurements. traces for key input and output variables for the VECSS
program.
Comparison between test no.1 and test no.2 from Tables
3,4 showed that with the ambient temperature being
Table 4. Comparison of predicted and experimental
close to the same values, the engine surface tempera-
stationary test parameters
tures for liner, head and piston increased with increase in
speed and percent full load of the engine. This also Test No.1 Test No.2 Test No.3
1500 RPM / 1800 RPM / 1800 RPM /
resulted in higher temperature rise for the coolant across 60% full 80% full 75% full
the engine. Due to this, the radiator top tank temperature Parameters
engine load / engine load engine load /
no RAM air / / no RAM air 25 m.p.h
was higher. The calculated and measured oil sump tem- 23 deg. C / 21 deg. C RAM air /
ambient ambient 17 deg. C
perature predictions for test no. 1 and 2 were also higher. ambient
The intake air temperature going into the engine was Pred. Expt. Pred. Expt. Pred. Expt.
higher at 80% load and 1800 rpm (test no. 2) as com- Engine Power (kW) 176 184 218 224 199 205
Exhaust temperature ( deg.C) 562 - 528 - 500 -
pared to 60% load and 1500 rpm (test no. 1).
Piston surface (deg.C) 313 - 328 - 312 -
The effects of change in ambient temperature and the Head surface ( deg.C) 253 - 261 - 249 -
Liner surface (deg.C) 155 - 161 - 155 -
RAM air can be studied by comparing test no. 3 and test
Coolant temperature rise across 3 - 3.5 - 3 -
no. 2 (Table 3,4). At the same engine speed and sightly engine (deg.C)
lower load, the engine liner, head and piston surface tem- Radiator top tank (deg. C) 93 93 94 94 91.5 91
peratures were lower at lower ambient temperatures. Radiator bottom tank (deg.C) 56 - 57 - 51 -
Thermostat flow through radiator (%) 9.9 - 10.5 - 9 -
This is because, for lower ambient temperatures (17°C Total time for which radiator thermostat 2163 - 2148 - 2827 -
(62°F) for test no. 3 and 21°C (70°F) for test no. 2) the was active (sec)
intake air temperature is lower for the same cooling air- Intake Manifold (deg. C) 38 46 50 55 44 49
Boost pressure (k Pa) 176 - 196 - 197 -
flow. The lower engine surface temperatures gave
Intake air mass (kg/min) 18.6 - 26 - 24 -
decreased coolant temperature rise across the engine A/F ratio 31.5 - 35 - 36 -
because of lower heat transfer to the coolant. Lower cool- Oil Sump (deg.C) 107 108 110 111 109 109
ant temperature results in a lower percentage opening of Oil thermostat factor open (%) 13.8 - 14.4 - 15 -
the thermostat (9 and 10.5%) and hence less coolant Actual air flow (ambient temp., m3/min) 125 - 152 - 141 -
Cooling air mass (kg/min) 102 - 126 - 120 -
flow through the radiator. Also, fan speed which is a func-
Fan speed (rpm) 1140 - 1480 - 1170 -
tion of the coolant temperature is lower for test no. 3. Time for which fan was fully engaged 0 - 0 - 0 -
Hence, although test no. 3 has RAM air of 25 m.p.h and (sec)
test no. 2 had no airflow, the cooling air mass remains Simulation real time (sec) - 2720 - 2640 - 3530
CPU time required (sec) 7800 - 7780 - 10,500 -
nearly the same. The oil sump temperature was also
slightly lower (109 and 110°C (228 and 230°F) from the
model) due to lower surface temperatures in the engine In Figure 3 the first three time traces show the engine
and higher heat rejection from the oil sump to the lower speed, vehicle velocity and the fuel flow rate which are
ambient air temperature. For the intake air temperature, input data to the VECSS program.
the lower ambient temperature results in lower intake air
temperature (44 and 50°C (111°F and 122°F)).
5
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
6
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
The predicted engine power follows the fuel flow rate The predicted top tank and oil sump temperatures agree
curve trend. One representative parameter for the engine closely with the experimental data both qualitatively and
operating conditions, the piston surface temperature, is quantitatively within 1°C (2°F) and 3.5°C (6°F) respec-
shown in Figure 3. Time traces for intake air temperature tively. Intake air temperature is predicted within 2.5°C
and mass flow rate are also shown. The ripples similar to (4.5°F).
the ones seen in the stationary test time traces (Figure 3)
Figures 6 gives the comparison for the oil sump tempera-
in the top tank temperature time trace are due to the
ture and Figure 7 gives the intake manifold air tempera-
opening and closing of the thermostat. The oil sump tem-
ture for the predicted and experimental data. The oil
perature slowly increases from the ambient temperature
temperature is initially at ambient temperature. The leak-
as the test time progresses. When it reaches 99°C
age through the oil thermostat is assumed at 5%. As the
(210°F), the oil thermostat begins to activate. The work-
oil temperature increases to 99°C (210°F), the oil thermo-
ing of the thermostat maintains the oil sump temperature
stat begins to open. The working of the oil thermostat
between 98 and 101°C (208 and 214°F). The predicted
maintains the average oil sump temperature at 97°C
airflow across the radiator follows the same trend as the
(207°F). The average temperature for the experimentally
vehicle velocity, since the fan speed which is a function of
measured data is 100°C (212°F). The trend in general
the coolant temperature and the engine speed, runs at a
lies within 3°C (5°F) qualitatively. The time averaged
nearly constant speed of 624 rpm as shown in Figure 4.
intake air temperature predicted by the VECSS program
is 27°C (81°F). The comparison of the trends for pre-
Table 5. Input parameters for the VECSS program for dicted and experimental data is shown in Figure 7. The
transient tests time averaged experimental data for intake air tempera-
Input parameter Test No.4 Test No.5 Test No.6
ture is 25.5°C (78°F), which is 1.5°C (3°F) less than the
predicted value.
Cadillac to Spring Lake Brighton to
Spring Lake to Plymouth Spring Lake
Comparison Plots for Test No 5 – The experimental
Engine speed (rpm) 1380 to 1450 1450 to 1500 1600 to 1650
test was run at an ambient temperature of 1°C (34°F)
Fuel flow rate range (kg/min) 0 to 1.2 0 to 1.2 0 to 1.2 with a load of 18712 kg (41,250 lbs) (75% of full load on
Average fuel flow (kg/min) 0.48 0.54 0.57
vehicle). The VECSS program took 11,400 secs of CPU
time to simulate 4200 secs of real time. For the top tank,
Cruising vehicle speed (m.p.h.) 57-60 60-65 60-65
oil sump, and the intake air temperatures, the time aver-
Ambient temperature (deg. C) 19 1 15.5 aged predicted and experimental data are 87 and 88, 97
% load on vehicle (24949.2 kg full 55 75 100
and 99.5, 11 and 13.5 deg. C(190 and 190, 207 and 211,
load) 51 and 56 deg. F) respectively.
Total simulation real time (sec) 5000 4000 4000
Comparison Plots for Test No 6 – A load of 24949 kg
(55,000 lbs) (100% full load on the vehicle) was carried
Comparison Plots for Test No 4 – For this test, the out at an ambient temperature of 15.5°C (60°F) for this
truck carried a load of 13609 kg (30,000 lbs), which was experimental test run. The VECSS program took 11,500
equal to 55% of the full load (24949 kg (55,000 lbs)) it secs of CPU time to simulate 4500 secs of real time. For
can carry on the vehicle. The VECSS program took the top tank, oil sump and the intake air temperatures,
14,250 secs of CPU time to simulate 5000 secs. of real the time averaged predicted and the experimental data
time. Figure 5 shows the comparison between the are 87 and 87, 97 and 100.5, and 22 and 24°C (189 and
VECSS predicted radiator top tank coolant temperature 189, 207 and 213, and 72 and 75°F) respectively.
and the experimental data. The time trace at the bottom
shows the percent opening of the coolant thermostat. As Table 6 gives the summary of the results obtained for the
the coolant temperature reaches 88°C (190°F), the ther- transient tests. This table presents predicted and experi-
mostat begins to open. The ripples in the coolant temper- mental data either in a range or as a time averaged value.
ature are due to initial opening and closing of the coolant With increase in the load on the vehicle, the average fuel
thermostat. The coolant temperature reaches an average flow rate increases, increasing the average power for test
temperature of 87°C (189°F). The small fluctuations in no. 4, test no. 5, and test no. 6 from 133 to 150 to 160
the temperature in the time range 1000-5000 seconds kW(181 to 205 to 218 hp). This conclusion can be made
show the continuous working of the thermostat. A similar because the load due to terrain effects (up and down hill
trend is seen in the experimental data. The average tem- effects) is probably averaged out at about the same value
perature for the experimental data is 88°C (190°F). for all the three tests although some of the differences in
power could be from terrain differences for the three runs.
7
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
In Table 6 although the load on the vehicle for test no. 5 is ing test no. 4 and 6 substantiate the above reasoning.
more than that of test no. 4, the piston, head, and the Test no. 6 is at 100% load on the vehicle with an ambient
liner temperatures are lower for test no. 5. This can be temperature of 15.5°C (60°F) for which the predicted
explained as follows. The ambient temperature is lower intake air temperature is 22°C (72°F). This is slightly
for test no. 5 resulting in lower intake air temperature. lower than the intake air temperature for test no. 4 and
Lower intake air temperatures even with increased power hence the engine surface temperatures are also slightly
resulted in lower engine surface temperatures. Compar- lower.
Total time for which radiator thermostat was active (sec) 4570 - 3600 - 4190 -
Actual cooling airflow range (at ambient temp. m3/min) 100-118 - 130-144 - 131-142 -
8
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
Figure 5. Comparison of top tank temperature and thermostat % open for Test No. 4
Figure 6. Comparison of sump temperature and thermostat % open for Test No. 4
9
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
Lower surface temperatures in the engine give lower tem- a high ambient temperature of 43°C (110°F) and at a low
perature rise in the coolant flowing across the engine ambient temperature of -29°C (-20°F). The driving cycle
because of lower heat transfer from the engine (refer fuel rate and engine speed used for these tests was
Table 6). At 75% load on the vehicle and lower ambient based on stationary test no. 2 which is at 21°C (70°F).
temperature (test no. 5), the radiator top tank remains the Test no. 2 used 0 m.p.h. ram air speed while these runs
same and the radiator bottom tank temperature is much used 58 m.p.h. to simulate the vehicle travel at this
lower. speed.
The oil sump temperature remains almost constant for all All the three tests were run with 50/50 ethylene glycol
the three tests. More than the load on the vehicle or heat and water as the coolant. The air-conditioning and the
rejection to the ambient temperature from the sump, the cab heater were in the off condition. No shutters were
oil thermostat plays a dominant role to maintain the tem- used for the tests. The fan used was variable speed
peratures constant. clutch.
For the tests number 4, 5, and 6, the fan operates at an Table 7 gives the input parameter values to the VECSS
average speed of 624, 660, and 659 rpm. respectively. program for the three runs.
Hence, the cooling airflow varies as the vehicle speed
Table 8 gives the summary of the results obtained from
changes. The vehicle speed and fan speed for test no. 4
the three tests.
are less than that for test no. 5, which results in lower air-
flow for test no. 4. But, for test no. 4 & 5, the average
vehicle speed and fan speed are the same. Since the Table 7. Input parameters for ambient temperature tests
ambient temperature is lower for test no. 5, the airflow is
Low Moderate High
higher as the air is more dense. The system resistance, Input parameter
Ambient Ambient Ambient
which is proportional to the air density decreases giving
higher volume flow rate. Engine speed (rpm) 1800 1800 1800
10
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
11
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
Figure 8. Trends observed from ambient tests Figure 9. Trends observed from ambient tests
(continued)
12
Downloaded from SAE International by Stony Brook Univ, Monday, August 13, 2018
ACKNOWLEDGEMENTS
13