This study evaluates the thermal performance of a Z-type battery thermal management system (BTMS) designed for nine lithium-ion batteries discharged at a high rate of 5C, using Computational Fluid Dynamics (CFD) simulations. The investigation employs Response Surface Methodology (RSM) to optimize two critical thermal performance parameters: the maximum battery temperature (Tmax) and the maximum temperature difference between cells (ΔTmax). Various cooling strategies are explored to comprehensively assess the BTMS, including natural convection, forced convection, cooling fins, phase change material (PCM), and composite PCM. These methods are analyzed to determine their effectiveness in controlling the thermal behavior of the battery pack. The simulation results indicate that integrating different cooling techniques can significantly lower Tmax from 352.38 K to 309.14 K and reduce ΔTmax from 14.6 K to 3.31 K, depending on the method used. Under critical conditions, such as the failure of the active cooling system, the BTMS still maintained a Tmax of 310.64 K and a ΔTmax of 0.95 K, demonstrating its robustness and reliability. Further optimization identified the ideal configuration for the system, including an inlet air speed of 1.2 m/s, an inlet temperature of 297.15 K, and a PCM thickness of 3.8 mm, achieving optimal thermal performance with a Tmax of 303.97 K and ΔTmax of 3.17 K. This study offers valuable insights into the design and optimization of effective BTMS for enhanced battery safety and longevity.