Tesla/EV Battery Thermal Runaway Simulation

 Simulating Thermal Runaway in electric vehicle (EV) batteries within Ansys Fluent is primarily handled through the Thermal Abuse Model found within the Battery Module. This model allows you to simulate the violent exothermic reactions that occur when a battery reaches critical temperatures due to internal shorts, overcharging, or external heating.

1. Core Kinetics Models for Thermal Runaway

Ansys Fluent provides three levels of complexity for modeling the chemical kinetics of thermal runaway:

• One-Equation Kinetics Model: A simplified model where the abuse reaction rate can be defined using an Arrhenius formulation, a Table Data lookup (defining \(dT/dt\) vs. temperature), or a User-Defined Function (UDF).
• Four-Equation Kinetics Model: A more detailed physics-based model that explicitly accounts for four key decomposition reactions:
  1. SEI Decomposition: The breakdown of the solid-electrolyte interface.
  2. Negative-Solvent Reaction: Reaction between the anode and the electrolyte.
  3. Positive-Solvent Reaction: Reaction between the cathode and the electrolyte.
  4. Electrolyte Decomposition: The final breakdown of the electrolyte fluid.
• General N-Equation Model: Allows for an arbitrary number of abuse reactions, each defined by eight rate parameters, including specific heat release (\(H\)), activation energy (\(E\)), and pre-exponential factors (\(A\)).

2. Setup Procedure in Ansys Fluent

To initiate a thermal runaway simulation, follow these technical steps:

1. Enable the Battery Model: Navigate to Setup → Models → Battery Model and select Enable Battery Model.
2. Activate Thermal Abuse: In the Advanced Options tab, check the Thermal Abuse Model box.
3. Select Kinetics Model: Choose between the One-Equation, Four-Equation, or N-Equation models based on the available data for your specific battery chemistry (e.g., LiCoO2, LiFePO4).
4. Import Material Data: For the Four-Equation model, you can use the Material Database for Abuse Kinetics to copy standard kinetics parameters directly into your case.
5. Internal Short Circuit (Optional): If the runaway is triggered by a short, enable Internal Short Heat. You must specify a trigger temperature (\(T_{trigger}\)) and the heat of the short-circuit reaction.
6. Venting Model (Optional): To simulate the release of gases during runaway, enable the Battery Venting Model in the Advanced Options. You can define the venting gas composition (e.g., \(H_2, CO_2, CH_4\)) and the release faces.

3. Advanced Simulation Techniques

• Standalone Preview: Before running a full 3D CFD simulation, use the Preview Thermal Abuse Model button. This opens a standalone solver to compute and plot the temperature profile and heating rate based on your parameters, allowing for rapid calibration against experimental data (like ARC tests).
• Reduced Order Models (ROM): For large battery packs (like those in a Tesla), a full transient 3D simulation may be computationally prohibitive. You can use the Battery ROM Tool Kit to generate LTI-ROM or SVD-ROM training data, which provides high-fidelity results at a fraction of the time.
• Coupled Solver: Ansys Fluent automatically switches to the transient solver when the battery model is enabled, which is necessary to capture the rapid temperature spikes characteristic of thermal runaway.

4. Technical Best Practices

• Energy Source Option: Ensure Enable E-Chem Heat Source is active to account for the electrochemical heat alongside the abuse heat.
• Precision: Always use the Double Precision solver, as the high energy gradients during runaway can lead to numerical instability in single precision.
• Time-Stepping: Use a small time step (often in the range of milliseconds) during the runaway phase to ensure convergence of the highly non-linear Arrhenius terms.

To provide you with the most effective workflow, here are additional advanced tips for Ansys Fluent and Mechanical APDL, ranging from mesh optimization to solver stability and specialized physics setup.

1. Mesh Quality and Geometry Tips

• Polyhedral Conversion: For tetrahedral meshes, converting the entire domain to a polyhedral mesh generally results in a lower cell count and faster convergence, which can save significant computational expense.
• Node-Based Gradients: Use node-based gradients when working with unstructured tetrahedral meshes. This scheme is notably more accurate than the default cell-based scheme for these mesh types.
• Boundary Alignment: Define your inlets and outlets before meshing to ensure they are prioritized. Avoid using the "Normal to Boundary" velocity method if the resulting mesh faces are misaligned with the original geometry, as this can impair convergence.
• Mesh Repair: If you encounter high aspect ratio cells, use the text user interface (TUI) command mesh/repair-improve/repair-face-handedness to improve the mesh quality.
• Local Remeshing: Instead of standard smoothing or collapsing, you can use the hotkey Ctrl+Shift+R for a Local Remesh to remove specific skewed faces.

2. Solver Stability and Convergence

• Energy Residuals: For pressure-based solvers, the scaled energy residual should decrease to at least \(10^{-6}\). For high-precision CHT or radiation tasks, aiming for \(10^{-8}\) is recommended to ensure the energy balance is fully closed.
• Steady-State Divergence: If you face convergence issues in a steady-state flow using the SIMPLE scheme, try reducing the under-relaxation factors to 0.2 for Pressure and 0.4 for Momentum and Turbulence.
• Transient Accuracy: If using second-order time accuracy, avoid changing the time step size by more than 20% during the calculation to maintain temporal accuracy.
• Residual Smoothing: In density-based explicit solvers, you can use implicit residual smoothing (with a factor of 0.5 and 2 passes) to allow the Courant (CFL) number to be doubled.
• PISO for Transient: The PISO algorithm with neighbor correction is highly recommended for transient calculations where you wish to use larger time steps.

3. Thermal and Radiation Specialized Tips

• Coupled Boundaries: In Conjugate Heat Transfer (CHT) problems, always maintain the Coupled setting for walls or shadows to ensure field continuity.
• Radiation Discretization: For the Discrete Ordinates (DO) model, the default of 2 for Theta and Phi divisions is often too coarse. Increasing these to 3 or up to 5 will provide much more reliable results for complex geometries.
• View Factor Accuracy: In Mechanical APDL, you can increase the accuracy of 3D view factor calculations by increasing the HEMIOPT resolution from the default of 10.
• Closed Systems: Periodically check that the summation of view factors for any radiating surface in a closed system equals 1.0; a value higher than this indicates a calculation error, often due to blocking effects.

4. Advanced Physics and Initialization

• Hypersonic Initialization: For Mach 10+ flows, Full Multigrid (FMG) Initialization is essential. It is recommended to lower the CFL to 0.5 or 0.25 during this initialization to enhance robustness.
• Turbulence Initialization: Always start from a suitable state of turbulence. Typical values for initialization should be between 1-10% for Turbulent Intensity and 1-100 for the Viscous Ratio.
• Incompressible Ideal Gas: When using the incompressible ideal gas law, always double-check the Operating Density setting, as it directly impacts the body-force term in the momentum equations.
• Phase Change Analysis: In Mechanical APDL, use lower-order thermal elements (such as PLANE55 or SOLID70) and enable automatic time stepping to allow the program to adjust for the sudden changes during phase transformations.

5. UI and Productivity Shortcuts (Fluent)

• Verify Normals: Use Ctrl+Shift+C followed by Ctrl+N to color the mesh by normal direction, allowing you to visually verify if they are oriented correctly for prism growth or boundary setup.
• Range Selection: In any multiple-selection list within the Fluent UI, you can select a range of items by selecting the first and then Shift-clicking the last.
• Double Buffering: Enable Double Buffering in graphics settings to dramatically reduce screen flicker during visual updates.