Project
In Construction

Battery Cooling Testing

How this project demonstrates each skill

First-Principles Problem Solving
Skill

Discharged cells, built a test rig, collected thermocouple data, and derived an h-correlation from it.

Hands-On & Shop Experience
Skill

Built the physical airflow test rig. Has a detailed writeup already (link to project page).

Physical Testing & Instrumentation
Skill

Chose experiment over CFD because the geometry was too irregular to trust a sim. Derived h-coefficient from raw thermocouple and power data. De-risked accumulator overheating during endurance.


The Summary

Decision: 6 W server fan operating point. Delivers , holding cells under 60 °C up to 39 °C ambient, which covers nearly all FSAE competition conditions. Enough margin to downgrade from a 10k to a 7k RPM fan without risk.

What it eliminated: Whether our pack needed active cooling at all, how much fan power, and which fan topology. Server fans outperforming blowers at equal power was the unexpected result that changed the hardware selection.

Evidence boundary: The h-correlation itself is experimentally validated across four power levels. Max safe ambient estimates depend on mass flow rate, which was back-calculated from h rather than measured directly.

Motivation

Over the course of endurance (22 km, ~25 minutes) our battery pack has to maintain temperatures under 60 °C (FSAE rules for maximum cell temperature). My cell and race profile modeling tools showed either cooling or power limiting must be done to keep our pack within temperatures.

Action

It was determined that doing physical airflow testing on our battery would produce higher quality results faster than CFD. The result of this experiment was creating a correlation between fan power and the heat transfer coefficient (cooling ability).

Result

Fan power vs heat transfer coefficient correlation
Fan power vs heat transfer coefficient — the key deliverable from physical testing

The correlation allows us to predict cooling performance at any fan power setting. At 10 W fan power, we achieve . However, the h-coefficient alone doesn't tell the whole story. As air flows past each row of cells, it absorbs heat and warms up, reducing cooling effectiveness for downstream rows. Our analysis accounts for this cumulative air heating effect to determine realistic “max safe ambient” temperatures, giving us high confidence that our cooling system can maintain cell temperatures below 60 °C even in unusually hot conditions.

Why Physical Testing?

CAD view of Enepaq module showing irregular airflow path
Front view of our module design (23 Enepaq bricks in a module)
Single Enepaq brick structure
Single Enepaq brick structure

The irregular geometry of Enepaq modules makes analytical heat transfer calculations (from array cooling textbook examples) unreliable. Physical testing was the fastest path to valid results.


Understanding

The fundamental energy balance for cell cooling:

During cooldown testing, we simplify by eliminating heat generation and assuming convection dominates:

Which gives us Newton's law of cooling with thermal mass:

CFD vs Experimentation

This project originally started because I set out to create realistic cooling estimates for our battery. I had a brief stint of trying to set up thermal CFD but quickly realized that I was neither confident in heat transfer fundamentals (started this in 2024), nor understood CFD solvers. I ended up asking the FSAE discord and Ethan Perrin (Tesla Battery Engineer, personal communication) pointed me toward IRL physical testing being much more likely to give valid results much faster.

Heat Transfer First Principles and MATLAB Files

The start of this project was in the semester before I had taken my university's heat transfer course. Before I designed any test setups or plans I worked on a couple thermal models in MATLAB.

For the algebra behind how scales with velocity and fan power, see the Nusselt correlation derivation — it predicts roughly 16× fan power to double cooling, which matches what this experiment measured.

Remaining Unknowns

The MATLAB models gave me fundamental understanding, but several parameters remained uncertain. I have ideas to look deeper into these, but more pressing priorities take up my time for the amount of benefit these would give.

  • Airflow velocity— fan specs don't translate directly to air speed over cells
  • Turbulence effects— how does fan-induced turbulence affect heat transfer?
  • Pressure drop— what's the system impedance across the module?
  • Cumulative error— how do uncertainties compound in detailed models?

These unknowns motivated physical testing over analytical approaches.


Designing

The goal: create a mock Enepaq brick that could be heated and have air blown past it. Success criteria:

  • Build quickly
  • Approximate thermal characteristics of real batteries
  • Test quickly and repeatably

Aluminum turned out to be ideal. It has similar density and heat capacity to lithium-ion cells at room temperature. The problem was I could only source imperial stock in Canada (5/8" vs 18 mm cells).

I wrote a dynamic similarity script (see MATLAB dropdown above) to check if imperial stock would skew results. Verdict: close enough for valid data, but 18 mm stock would make the mock modules dimensionally identical to real ones. Worth the wait.

The AliExpress Saga

University purchasing delays pushed this project back over a month. I ended up buying the aluminum myself from AliExpress. I designed two mock brick versions:

  • 3D printed: Cheap, fast. Mainly for validating system airflow impedance
  • Laser cut: Wood + aluminum rods. Enables oven preheating to test progressive air heating

The designs worked surprisingly well. The downside was by the time everything arrived and laser cutting was done, finals season was approaching.

Mock Brick Design

Laser-cut mock Enepaq brick with aluminum rods on desk
Laser-cut mock Enepaq brick with aluminum rods. Similar thermal mass to real 18650 cells

Testing

The physical testing can be broken up into two categories:

  • Aluminum mock Enepaq brick (heated externally by heat gun)
  • Real Enepaq brick (heated by electrical discharge)

Scrappy Test Setup

Temperature measurement used K-type thermocouples with thermal paste, taped to the aluminum surface. Data was logged via a Chinese TC logger to CSV files.

Test setup showing thermocouple placement and airflow configuration
Test setup showing thermocouple placement and airflow configuration. This was the first day everything was printed out and the DAQ was working. This is NOT the setup used, it's my best photo that displays everything in one view.
Messy MATLAB plot from initial aluminum brick testing
Early testing with aluminum 'dummy' bricks — these initial experiments gave us unusually high h-coefficients due to unreliable thermocouple adherence and didn't account for the 'hottest cell rule' or progressive air heating across rows

Real Enepaq Brick Testing

Full battery cooling test setup with Enepaq brick, airflow segment, power supply, and data logging
Full test setup overview — Enepaq brick, airflow segment, power supply, and data logging

The Enepaq brick's built-in thermistors provided comparison data for the live cell tests. These internal sensors are shielded by the plastic casing, which means they lag behind actual cell surface temperature. Important context when comparing to thermocouple readings.

Test Parameters

ParameterValue
Cell rows6
Cells per row4
Segment bricks1 (of 23 in full pack)
Discharge current35 A constant
Start voltage4.185 V
End voltage2.8 V (most tests)
RunFan TypePowerMedian h*
Run 110k Server Fan2.75 W13.4 W/m²K
Run 210k Server Fan6.0 W21.0 W/m²K
Run 310k Server Fan10.0 W23.6 W/m²K
Run 410k Server Fan16.0 W26.2 W/m²K
Run 55015 Blower2.75 W9.3 W/m²K

*h-coefficient calculated from the hottest cell temperature — this is the conservative approach since we care about the worst-case cell.

Why a 10k Server Fan?

A server fan is massive overkill for this application, but it was chosen deliberately to give us flexibility during testing. We could easily characterize performance across a wide power range and have confidence we'd find a workable operating point and better correlation.

Representative Test Analysis

Run 3 four-axis analysis plot
Run 3 (10k server fan at 10 W): Four-axis plot showing h-coefficient stabilization, heat loss, temperature decay, and battery (recovery) voltage over time

The h-coefficient calculation uses Newton's law of cooling:

Where:

  • = thermal mass of cells (5 cells × 49 g × 900 J/kg·K)
  • = exposed surface area (90% of cylindrical surface)
  • = temperature difference between cell and ambient

Comparison Across Runs

Overlay of h-coefficient vs time for all test runs
Overlay of h-coefficient vs time for all test runs. Faster cooling = shorter test time. Plots truncated to keep clean data. Median h-coefficient is taken from these graphs to avoid the effect of h rising due to thermal inertia.

Creating Takeaways

Purposeful Pessimism

This is UVic FSAE's first year with this battery architecture. We're applying conservative safety factors because:

  1. Overheating = DNF— exceeding 60 °C during endurance means we fail
  2. Data quality will improve— once the car is running, we'll get real-world validation
  3. We have the budget — our LV power budget can accommodate oversized cooling if needed

Better to be overly cautious now and dial back later.

The Correlation

Power law fit for h-coefficient vs fan power
Power law fit: h = 9.97 × P^0.462 for the 10k server fan configuration. Blower immediately seen as less effective even at its max nominal power draw

This correlation is the primary deliverable. It allows predicting heat transfer coefficient for any fan power setting without additional physical testing.

Golden Table (with High Safety Margins)

RunFan ConfigPowerMedian hCell ΔTAir HeatingMax Safe Ambient
Run 110k Server2.75 W13.4332.5 °C+11.0 °C16.5 °C
Run 210k Server6.0 W21.0015.3 °C+5.4 °C39.3 °C
Run 310k Server10.0 W23.5713.9 °C+4.4 °C41.8 °C
Run 410k Server16.0 W26.2312.4 °C+3.5 °C44.1 °C
Run 5Blower2.75 W9.2928.8 °C+15.0 °C16.3 °C

Max Safe Ambient = 60 °C limit − Cell ΔT − Air Heating offset.
Air heating estimates depend heavily on mass flow rate. Mass flow wasn't measured directly — it was back-calculated from the h-coefficient correlation to airspeed. This table will be updated once I measure actual volumetric flow rate with my “low-backpressure trash bag air catcher” (yes, it's exactly what it sounds like).

Air Temperature Gradient

Estimated air temperature rise across 6 cell rows
Estimated air temperature rise across 6 cell rows — higher fan power reduces cumulative heating

Mass Flow Measurement Pending

These curves rely on estimated mass flow. Direct measurement will improve accuracy significantly.

Design Implications

  1. Fan selection: Server fans actually outperform blowers at equivalent power — this was unexpected. The 10k server fan handles both static pressure and high airflow well.
  2. Operating point: 6 W fan power provides , sufficient for almost all race conditions (ambient < 39 °C).
  3. Blowers inadequate: Same power input yields ~30% lower h-coefficient, unexpected result.
  4. Room to optimize: Very confident we can downgrade to a 7k RPM fan and still meet thermal requirements with margin.

Analysis Code

Processing pipeline:

  • Smoothing for noise reduction
  • Newton's law of cooling for h-coefficient extraction
  • Power law curve fitting for correlation
  • Row-by-row air temperature gradient simulation

Future Work

  • Weather statistics to make decisions on ambient conditions at competition venues
  • Effect of additional heating from sun exposure and hot pavement
  • Pre-cooling strategies (ice intake air) for extreme conditions
  • Trash-bag air measurement test