19.1.4. MSMD Solution Method

The difficulty with modeling a lithium-ion (Li-ion) battery is due to its multi-domain, multi-physics nature. Vastly different length scales associated with different physics complicates the problem. When performing a thermal analysis, the goal is to determine the temperature distribution at the battery length scale. The physics governing the Li-ion transport occurs in the anode-separator-cathode sandwich layers (the electrode pair length scale). Li-ion transport in an active material occurs at the atomic length scale. The Multi-Scale Multi-Domain (MSMD) approach deals with different physics in different solution domains [[312]].

Battery thermal and electrical fields are solved in the CFD domain at the battery cell's scale using the following differential equations:

(19–3)

(19–4)

where and are the effective electric conductivities for the positive and negative electrodes, and are phase potentials for the positive and negative electrodes, and are the volumetric current transfer rate and the electrochemical reaction heat due to electrochemical reactions, respectively, and are the current transfer rate and heat generation rate due to battery internal short-circuit, respectively, and is the heat generation due to the thermal runaway reactions under the thermal abuse condition.

For normal operation, is set to zero. For more information, refer to Thermal Abuse Model. The source terms and are computed using an electrochemical submodel. If there is no internal short-circuit, and are equal to zero. For more information, refer to External and Internal Electric Short-Circuit Treatment.

A wide range of electrochemical models, from simple empirically-based to fundamental physics-based, is available in the open literature. In Ansys Fluent, the following electrochemical submodels are implemented:

  • Newman, Tiedemann, Gu and Kim (NTGK) or direct current internal resistance (DCIR) model

  • Equivalent Circuit Model (ECM) model

  • Newman Pseudo-2D (P2D) model

In addition, you can define your own electrochemical model via user-defined functions and hook it to the Fluent MSMD battery module.

To use the MSMD approach in an arbitrary finite volume of the cell composite, the following conditions must be met:

  • Battery micro layers must have the same orientation

  • Two potential fields, and , must be uniquely determined

These conditions are usually met for aligned stack cells or wound cells with extended foil-type continuous current tabs.

The MSMD solution method can resolve the potential distribution within a battery and, therefore, the detailed heat generation distribution within the battery. Since two extra transport equations are solved, the computational cost is higher compared to the previous methods.