Bối cảnh và vấn đề nghiên cứu

Trong xu thế phát triển giao thông vận tải bền vững, ô tô điện (Electric Vehicle - EV) đang trở thành giải pháp tối ưu nhằm giảm thiểu phát thải khí nhà kính và giải quyết bài toán khủng hoảng nhiên liệu hóa thạch. Trong hệ thống truyền động điện, nguồn tích trữ năng lượng đóng vai trò quyết định đến quãng đường di chuyển, công suất vận hành và tính an toàn của toàn xe. Pin Lithium-ion hiện là công nghệ lưu trữ năng lượng chủ đạo trên ô tô điện nhờ sở hữu mật độ năng lượng cao, tuổi thọ chu kỳ dài và khả năng nạp xả tái tạo hiệu quả. Mặc dù vậy, pin Lithium-ion vẫn tồn tại những rào cản kỹ thuật nghiêm trọng như sự suy giảm dung lượng sau quá trình làm việc liên tục và rủi ro mất kiểm soát nhiệt, cháy nổ dưới điều kiện nhiệt độ cao hoặc tải dòng lớn. Do đó, việc nghiên cứu các cơ chế động học điện hóa bên trong pin để mô hình hóa và mô phỏng chính xác hành vi làm việc trở thành yêu cầu cấp thiết, giúp tối ưu hóa thiết kế hình học và nâng cao hiệu quả quản lý năng lượng.

Đồ án tốt nghiệp "Nghiên cứu, mô hình hóa và mô phỏng pin Lithium - ion sử dụng cho ô tô điện bằng phần mềm Matlab/Simulink" của sinh viên Nguyễn Văn Dũng và Đỗ Nguyễn Trung Kiên (Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh, năm 2025) được thực hiện nhằm giải quyết các nhiệm vụ và mục tiêu cụ thể sau:

  1. Tìm hiểu và nắm vững nguyên lý hoạt động điện hóa của pin Lithium-ion trong cả hai chế độ nạp và phóng điện.
  2. Phân tích và dự đoán sự biến thiên của các thông số nội tại trong pin Lithium-ion sử dụng cực âm than chì thông qua việc thiết lập và giải toán hệ phương trình động học và động lực học vận chuyển ion.
  3. Xây dựng chương trình mô phỏng trên nền tảng số cho phép tùy biến thay đổi các thông số vật lý hình học (như độ dày các lớp điện cực), đối chiếu và kiểm chứng kết quả thu được với các công trình nghiên cứu khoa học đã công bố để xác định thông số tối ưu giúp pin đạt hiệu suất vận hành cao nhất.

Đối tượng nghiên cứu trọng tâm của đề tài là pin Lithium-ion sử dụng vật liệu cực âm là than chì (graphite - $\text{Li}_x\text{C}_6$) kết hợp với cực dương oxit kim loại chuyển tiếp Lithium Niken Mangan Coban ($\text{Li}_y\text{NiMnCoO}_2$ - NMC). Về phạm vi nghiên cứu, tác giả giới hạn khảo sát trên mô hình điện hóa một chiều (1D) với cấu trúc cực âm than chì, không mở rộng sang các hệ vật liệu cực âm khác. Không gian mô phỏng được thực hiện trên phần mềm tính toán kỹ thuật Matlab/Simulink trong điều kiện vận hành lý tưởng, đẳng nhiệt và không chịu tác động của các yếu tố nhiễu loạn môi trường bên ngoài.

Cơ sở lý thuyết và phương pháp

Khung lý thuyết của đồ án được xây dựng dựa trên mô hình điện hóa một chiều (1D electrochemical model) kế thừa từ lý thuyết điện cực xốp của Newman. Mô hình chia cấu trúc tế bào pin thành ba miền vật lý liên tiếp theo phương ngang: miền điện cực âm (anode), miền vách ngăn (separator) và miền điện cực dương (cathode). Các quá trình truyền chất và phản ứng bề mặt được mô tả thông qua hệ phương trình vi phân riêng phần (PDEs) phi tuyến:

  1. Phương trình bảo toàn khối lượng ion Lithium trong pha điện phân lỏng: $$\frac{\partial (\varepsilon_e C_e)}{\partial t} = \frac{\partial}{\partial x} \left( D_e^{\text{eff}} \frac{\partial C_e}{\partial x} \right) + \frac{1 - t_+}{F} j_{\text{Li}}$$ Trong đó, $C_e$ là nồng độ ion Lithium trong chất điện phân, $\varepsilon_e$ là phần thể tích chất điện phân, $D_e^{\text{eff}}$ là hệ số khuếch tán hiệu dụng, $t_+$ là hệ số chuyển dời của ion, $F$ là hằng số Faraday và $j_{\text{Li}}$ là mật độ dòng chuyển dời diện tích bề mặt.

  2. Phương trình bảo toàn điện tích trong pha rắn của điện cực: $$\frac{\partial}{\partial x} \left( \sigma^{\text{eff}} \frac{\partial \phi_s}{\partial x} động \right) = j_{\text{Li}}$$ Phương trình áp dụng cho hai miền điện cực rắn, với $\phi_s$ là điện thế pha rắn và $\sigma^{\text{eff}}$ là độ dẫn điện hiệu dụng của vật liệu điện cực.

  3. Phương trình bảo toàn điện tích trong pha lỏng của chất điện phân: $$\frac{\partial}{\partial x} \left( k^{\text{eff}} \frac{\partial \phi_e}{\partial x} \right) + \frac{\partial}{\partial x} \left( k_d^{\text{eff}} \frac{\partial \ln C_e}{\partial x} \right) = -j_{\text{Li}}$$ Phương trình mô tả sự phân bố điện thế pha lỏng $\phi_e$ trên toàn bộ miền pin, trong đó $k^{\text{eff}}$ là độ dẫn ion hiệu dụng và $k_d^{\text{eff}}$ là hệ số khuếch tán hiệu dụng.

  4. Phương trình động học bề mặt Butler-Volmer: Liên kết hai pha rắn - lỏng thông qua tốc độ phản ứng điện hóa tại ranh giới hạt: $$j_{\text{Li}} = a_s i_0 \left[ \exp\left( \frac{\alpha_a F}{RT} \eta \right) - \exp\left( -\frac{\alpha_c F}{RT} \eta \right) \right]$$ Với $a_s$ là diện tích bề mặt riêng của hạt điện cực, $i_0$ là mật độ dòng trao đổi, $\alpha_a$ và $\alpha_c$ là hệ số chuyển dời điện tích anode và cathode, $\eta$ là quá thế bề mặt được xác định bởi công thức $\eta = \phi_s - \phi_e - U$ ($U$ là điện thế cân bằng nhiệt động học).

  5. Phương trình điện áp làm việc và trạng thái sạc: Điện áp làm việc của pin được xác định qua hiệu điện thế pha rắn giữa hai bộ thu dòng kèm theo độ sụt áp tiếp xúc: $$V(t) = \phi_s(L, t) - \phi_s(0, t) - \frac{R_f}{A} I(t)$$ Trạng thái sạc (State of Charge - SOC) và độ xả sâu (Depth of Discharge - DOD) được xác định trên cơ sở nồng độ ion Lithium trung bình tích phân trong pha rắn và các tỷ lệ cân bằng hóa học bề mặt.

Phương pháp nghiên cứu thực nghiệm số của đồ án kết hợp ba phương pháp luận:

  • Phương pháp phân tích - tổng hợp tài liệu: Rút trích hệ thống tham số vật lý, hằng số động học điện hóa từ các công trình nghiên cứu khoa học chuyên sâu và thư viện mở (chẳng hạn như phần mềm Dandeliion).
  • Phương pháp rời rạc hóa thể tích hữu hạn (Finite Volume Method - FVM): Chia miền tính toán 1D thành các phân chia thể tích con (control volumes) trong hệ tọa độ Descartes, tích phân từng phương trình vi phân trên từng ô lưới để chuyển hệ phương trình vi phân liên tục thành hệ phương trình đại số.
  • Thuật toán giải ma trận ba đường chéo (Thomas Algorithm / TDMA): Tuyến tính hóa các thành phần nguồn phi tuyến $S = S_C + S_P \phi$, thiết lập ma trận ba đường chéo và giải theo hai bước khử tiến (forward substitution) và thế lùi (backward substitution).
  • Phương pháp so sánh, kiểm chứng: Đối chiếu trực tiếp các đường cong điện áp làm việc theo thời gian thu được từ mô hình Matlab với dữ liệu thực nghiệm đã công bố của mẫu pin NMC dung lượng 6Ah ở các tốc độ sạc/xả 1C, 2C và 10C.

Thiết kế và triển khai mô phỏng

Đồ án triển khai thuật toán tính toán trên phần mềm Matlab/Simulink thông qua chuỗi xử lý số học đồng bộ các phương trình vi phân từng phần:

Tham số / Phân vùng Miền cực âm (Anode) Miền vách ngăn (Separator) Miền cực dương (Cathode)
Kích thước độ dày $L_N$ $L_S$ $L_P$
Số miền thể tích con $n_N$ $n_S$ $n_P$
Vật liệu hoạt tính Than chì (Graphite - $\text{Li}_x\text{C}_6$) Màng Polymer vi sợi (PE/PP) Oxit $\text{Li}_y\text{NiMnCoO}_2$ (NMC)
Pha rắn $\phi_s$ Giải phương trình điện tích pha rắn Không giải ($\sigma^{\text{eff}} = 0$) Giải phương trình điện tích pha rắn
Pha lỏng $\phi_e$, $C_e$ Giải phương trình truyền chất/điện tích Giải phương trình truyền chất/điện tích Giải phương trình truyền chất/điện tích
Điều kiện biên thu dòng $-\sigma_-^{\text{eff}} \frac{\partial \phi_s}{\partial x}\Big|_{x=0} = \frac{I}{A}$ Cách điện pha rắn tại ranh giới $\sigma_+^{\text{eff}} \frac{\partial \phi_s}{\partial x}\Big|_{x=L} = \frac{I}{A}$
Điều kiện biên nồng độ $\frac{\partial C_e}{\partial x}\Big|_{x=0} = 0$ Liên tục tại ranh giới $\frac{\partial C_e}{\partial x}\Big|_{x=L} = 0$

Quy trình giải lặp tại mỗi bước thời gian $t \to t + \Delta t$ được tổ chức theo sơ đồ tuần tự:

  1. Tiếp nhận giá trị dòng điện nạp/xả $I(t)$ và trạng thái nồng độ hạt cầu bề mặt từ bước thời gian trước.
  2. Khai triển và giải hệ phương trình TDMA tìm điện thế pha rắn $\phi_s$ tại cực âm và cực dương.
  3. Tuyến tính hóa thành phần nguồn điện hóa và giải TDMA xác định điện thế chất điện phân $\phi_e$ trên toàn bộ ba miền ($L_N + L_S + L_P$).
  4. Giải phương trình khuếch tán bảo toàn khối lượng để tìm phân bố nồng độ ion $C_e$.
  5. Cập nhật mật độ dòng chuyển dời $j_{\text{Li}}$ theo Butler-Volmer, xác định điện thế cân bằng $U_N, U_P$, tính toán điện áp tế bào $V(t)$ và cập nhật trạng thái sạc SOC.

Nội dung chính theo từng chương

Chương 1. Giới thiệu đề tài

Chương mở đầu tập trung trình bày tính cấp thiết của việc phát triển xe điện và vị thế của pin Lithium-ion trong hệ thống lưu trữ năng lượng. Tác giả phân tích các hạn chế của pin hóa học khi làm việc ở chế độ tải biến thiên liên tục trên ô tô, nhấn mạnh sự cần thiết của các công cụ mô phỏng toán học nhằm cắt giảm chi phí thử nghiệm phá hủy trong phòng thí nghiệm. Chương này thiết lập tường minh các mục tiêu nghiên cứu, xác định rõ đối tượng khảo sát là pin Lithium-ion sử dụng cực âm than chì và công bố các kết quả dự kiến đạt được trên mô hình Matlab/Simulink.

Chương 2. Tổng quan về các loại pin sử dụng trên ô tô điện

Nội dung chương 2 khái quát toàn diện cơ sở kỹ thuật về nguồn điện hóa cho xe điện. Tác giả phân tích chuyên sâu các thông số hoạt động chính: dung lượng lý thuyết, dung lượng định mức, dung lượng thực tế; điện áp mạch hở ($U_{\text{OCV}}$), điện áp làm việc ($U_{\text{working}} = U_{\text{OCV}} - I \cdot R_{\text{Internal}}$), điện áp cắt xả (Cut-off Voltage); bản chất nội trở Ohmic và điện trở phân cực; tốc độ sạc/xả định mức ($C_{\text{rate}} = I / C$); mật độ công suất khối lượng ($P/m$); chu kỳ vòng đời (chuẩn suy giảm 80% dung lượng ban đầu).

Bên cạnh đó, chương 2 phân loại và so sánh cơ chế làm việc của ba thế hệ pin sạc chính:

  • Pin Niken-Cadmium (Ni-Cd): Cực dương $\text{NiOOH}$, cực âm $\text{Cd}$, dung dịch điện phân kiềm $\text{KOH}$. Pin có độ bền cơ học cao, chịu được môi trường khắc nghiệt nhưng nặng và chứa kim loại độc hại Cadmium gây ô nhiễm.
  • Pin Niken-Kim loại Hydrua (NiMH): Thay thế cực âm Cadmium bằng hợp kim hấp thụ Hydro ($MH$), thân thiện với môi trường hơn và dung lượng cao hơn Ni-Cd, song tỷ lệ tự phóng điện cao và dễ hư hỏng khi sạc quá mức.
  • Pin Lithium-ion: Cấu tạo bốn thành phần cốt lõi gồm cực âm than chì hoặc silic, cực dương oxit kim loại, màng ngăn polymer (PE/PP tráng gốm) và dung dịch điện phân chứa muối $\text{LiPF}_6$ hòa tan trong dung môi hữu cơ carbonat (EC, DMC, EMC). Tác giả so sánh chi tiết ưu nhược điểm của 5 họ pin Lithium thông dụng: $\text{LiCoO}_2$ (LCO), $\text{LiFePO}_4$ (LFP), $\text{LiMn}_2\text{O}_4$ (LMO), $\text{LiNiMnCoO}_2$ (NMC) và $\text{LiNiCoAlO}_2$ (NCA).
Đặc tính so sánh Pin Niken-Cadmium (Ni-Cd) Pin Niken-Kim loại Hydrua (NiMH) Pin Lithium-ion (NMC/LFP)
Vật liệu cực dương $\text{NiOOH}$ $\text{NiOOH}$ $\text{LiMO}_2$ (NMC, LCO) / $\text{LiFePO}_4$
Vật liệu cực âm Cadmium kim loại ($\text{Cd}$) Hợp kim kim loại tích trữ Hydro ($MH$) Than chì ($\text{C}$) hoặc Silic
Chất điện phân Dung dịch kiềm $\text{KOH}$ Dung dịch kiềm $\text{KOH}$ Muối $\text{LiPF}_6$ trong dung môi hữu cơ
Điện áp cắt xả $0,9 - 1,0\text{ V}$ Khoảng $0,9 - 1,0\text{ V}$ Khoảng $2,5 - 2,8\text{ V}$
Mật độ năng lượng Thấp Trung bình Rất cao
Tác động môi trường Độc hại cao (chứa Cadmium) Thân thiện hơn Ni-Cd Không chứa Cadmium, cần xử lý tái chế

Chương 3. Quy ước sự phân bố nồng độ ion Lithium ở điện cực

Chương 3 trình bày toán học hóa các hiện tượng truyền khối và điện hóa. Các giả định chính được thiết lập: nồng độ ion Lithium tại một vị trí tọa độ $x$ trên điện cực bằng nồng độ bề mặt của hạt hình cầu tại vị trí đó; các thông số vận chuyển hiệu dụng ($D_e^{\text{eff}}$, $k^{\text{eff}}$, $\sigma^{\text{eff}}$) được hiệu chỉnh thông qua độ xốp $\varepsilon$ theo tương quan Bruggeman. Tác giả dẫn xuất chi tiết phương trình bảo toàn khối lượng chất điện phân, hai phương trình bảo toàn điện tích, phương trình động học Butler-Volmer và các hàm điện thế cân bằng thực nghiệm $U_N$ và $U_P$ phụ thuộc vào tỷ lệ thành phần hóa học $x, y$ tại điện cực. Chương kết thúc bằng việc xây dựng mô hình toán xác định biến thiên SOC, DOD và điện áp đầu cực tế bào pin.

Chương 4. Phương pháp giải

Chương 4 tập trung giải quyết bài toán tính toán số cho hệ phương trình vi phân riêng phần phi tuyến:

  • Thiết lập lưới không gian 1D trong hệ tọa độ Descartes, chia các miền $L_N, L_S, L_P$ thành các thể tích kiểm soát hữu hạn có nút mạng $P$ ở tâm cùng hai nút lân cận Tây ($W$) và Đông ($E$).
  • Tích phân phương trình đạo hàm cấp hai dạng bảo toàn: $$\int_{x_W}^{x_E} \left[ \frac{d}{dx}\left( k \frac{d\phi}{dx} \right) + S \right] dx = 0 \Longleftrightarrow a_P \phi_P = b_E \phi_E + c_W \phi_W + d$$
  • Kỹ thuật tuyến tính hóa thành phần nguồn phi tuyến $S = S_C + S_P \phi$ để đảm bảo tính hội tụ đơn điệu của nghiệm số.
  • Triển khai thuật toán Thomas (TDMA): Khử tiến ma trận ba đường chéo về hệ phương trình truy hồi $\phi_{i-1} = P_{i-1} \phi_i + Q_{i-1}$, sau đó áp dụng phép thế lùi để tìm giá trị biến trạng thái tại tất cả các điểm lưới.
  • Xây dựng sơ đồ khối và thuật toán lập trình trong Matlab để lần lượt tính toán ba biến trường cốt lõi: điện thế pha rắn $\phi_s$, điện thế pha lỏng $\phi_e$ và nồng độ ion chất điện phân $C_e$ qua từng bước thời gian.

Chương 5. Kết quả mô phỏng và kết quả thu được

Chương 5 là phần thực nghiệm mô phỏng và kiểm chứng mô hình:

  • Đối tượng kiểm chứng là tế bào pin Lithium-ion NMC dung lượng thương mại 6Ah (được đặc trưng bởi Bảng thông số kỹ thuật 5.1 trong đồ án).
  • Thực hiện mô phỏng quá trình sạc và phóng điện đẳng dòng ở ba chế độ tải khác nhau: 1C (6A), 2C (12A) và 10C (60A).
  • Trích xuất và vẽ đồ thị so sánh đường cong điện áp làm việc (Cell voltage) theo thời gian; kết quả cho thấy sự tương thích chính xác với đường đặc tính tham chiếu từ các bài báo khoa học đã xuất bản.
  • Mô phỏng khảo sát độ nhạy hình học: Khảo sát ảnh hưởng của việc thay đổi độ dày điện cực ($L_N, L_P$) tới độ sụt áp và dung lượng khai thác tại 1C, 2C và 10C.
  • Phân tích chi tiết quy luật phân bố không gian và thời gian của điện thế pha rắn $\phi_s$, độ chênh lệch điện thế trong dung dịch điện phân $\phi_e$, cũng như gradien nồng độ ion Lithium $C_e$ trên toàn bộ chiều dài pin tại các thời điểm làm việc khác nhau.

Kết quả và đóng góp

Đồ án đã đạt được các kết quả kỹ thuật cụ thể, được ghi nhận qua các bảng số liệu và đồ thị trong báo cáo:

  • Xây dựng thành công bộ chương trình mô phỏng hoàn chỉnh trên Matlab tính toán đồng thời các biến trạng thái nội tại của pin Lithium-ion NMC 6Ah.
  • Kết quả mô phỏng đường cong điện áp làm việc của pin ở các tốc độ sạc và xả 1C, 2C và 10C hoàn toàn trùng khớp với dữ liệu thực nghiệm đã công bố trong các bài báo khoa học quốc tế dùng làm đối chứng. Khi tốc độ dòng xả tăng từ 1C lên 10C, độ sụt áp Ohmic và quá thế phân cực tăng mạnh, khiến điện áp hoạt động giảm nhanh chạm ngưỡng điện áp cắt xả trong thời gian ngắn hơn.
  • Kết quả khảo sát biến thiên thông số hình học chỉ ra rằng việc gia tăng độ dày lớp điện cực giúp gia tăng thể tích vật liệu hoạt tính tích trữ nhưng đồng thời làm gia tăng trở kháng khuếch tán nội tại, dẫn đến sự phân bố nồng độ ion Lithium $C_e$ kém đồng đều, làm tăng độ dốc suy giảm điện áp ở các chế độ xả dòng lớn (10C).
  • Trích xuất thành công sự phân bố không thời gian của điện thế pha rắn $\phi_s$ (tập trung độ dốc tại hai điện cực), điện thế pha lỏng $\phi_e$ (hình thành độ chênh lệch rõ rệt qua vách ngăn do trở kháng ion) và sự biến thiên gradien nồng độ $C_e$ giữa hai cực trong quá trình nạp và xả.

Đóng góp kỹ thuật của đề tài:

  • Tự thiết lập và làm chủ hoàn toàn thuật toán thể tích hữu hạn kết hợp giải thuật Thomas (TDMA) trên nền tảng lập trình Matlab, không phụ thuộc vào các module hộp đen thương mại đắt tiền.
  • Cung cấp một phương pháp luận tin cậy để đánh giá động thái điện - hóa học của tế bào pin xe điện, cho phép các kỹ sư ô tô tối ưu hóa kích thước hình học điện cực phục vụ việc chế tạo tế bào pin có tỷ số công suất/năng lượng tối ưu.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của đề tài:

  • Mô hình chỉ được giải trong điều kiện đẳng nhiệt lý tưởng (nhiệt độ phòng không đổi), chưa ghép nối với phương trình cân bằng nhiệt lượng để khảo sát hiện tượng sinh nhiệt do tổn hao Joule và nhiệt phản ứng điện hóa.
  • Đối tượng nghiên cứu giới hạn trên hệ vật liệu cực âm than chì và cực dương NMC dung lượng 6Ah; chưa mở rộng kiểm chứng cho các hóa chất pin khác như LFP hoặc cực âm hỗn hợp Si-Graphite.
  • Mô hình 1D giả định tính đồng nhất trên mặt phẳng tiết diện vuông góc, bỏ qua sự phân bố không đồng đều theo phương 2D/3D trong các cấu trúc cuộn (cylindrical) hoặc xếp lớp (prismatic/pouch).
  • Chưa xét đến các cơ chế lão hóa dài hạn theo chu kỳ như sự dày lên của lớp màng liên pha điện cực rắn (SEI) hay hiện tượng mạ Lithium (Lithium plating).

Hướng nghiên cứu tiếp được đề xuất:

  • Mở rộng mô hình điện hóa một chiều thành mô hình ghép nối điện hóa - nhiệt (coupled electrochemical-thermal model) 2D/3D nhằm dự báo chính xác trường nhiệt độ và nguy cơ thoát nhiệt trên bộ pin xe điện.
  • Tích hợp các phương trình mô tả cơ chế lão hóa (suy giảm dung lượng và tăng nội trở tiếp xúc $R_f, R_{\text{SEI}}$) theo thời gian sử dụng.
  • Ứng dụng mô hình để thiết kế các thuật toán ước lượng trạng thái sạc (SOC) và trạng thái sức khỏe (SOH) thời gian thực phục vụ hệ thống quản lý pin (BMS) trên ô tô điện.

Giá trị tham khảo

Đồ án tốt nghiệp là tài liệu tham khảo chuyên khảo hữu ích cho:

  • Sinh viên, học viên cao học và nghiên cứu sinh chuyên ngành Công nghệ Kỹ thuật Ô tô, Kỹ thuật Cơ điện tử, Kỹ thuật Năng lượng và Kỹ thuật Điện - Điện tử quan tâm đến lĩnh vực lưu trữ năng lượng và phương tiện giao thông điện hóa.
  • Các nhóm nghiên cứu làm đề tài về mô hình hóa pin: Phần trình bày tường minh các bước biến đổi toán học phương trình vi phân, kỹ thuật rời rạc thể tích hữu hạn và thuật toán ma trận ba đường chéo TDMA tại Chương 3 và Chương 4 có giá trị thực hành lập trình tính toán cao.
  • Các kỹ sư phát triển phần mềm nhúng BMS: Bảng thông số vật lý của pin NMC 6Ah và các phương trình đặc tính điện thế cân bằng cung cấp bộ dữ liệu đầu vào chuẩn xác để mô phỏng kiểm thử các thuật toán điều khiển sạc xả.

Câu hỏi thường gặp

1. Mô hình điện hóa một chiều trong đồ án giải quyết những phương trình vi phân cốt lõi nào? Mô hình tập trung giải quyết đồng thời bốn phương trình vi phân riêng phần cơ bản: phương trình bảo toàn khối lượng ion Lithium trong chất điện phân, phương trình bảo toàn điện tích trong điện cực rắn, phương trình bảo toàn điện tích trong chất điện phân và phương trình động học điện hóa Butler-Volmer liên kết mật độ dòng chuyển dời tại bề mặt phân chia pha.

2. Thuật toán ma trận ba đường chéo (TDMA) được ứng dụng trong đồ án nhằm mục đích gì? Sau khi tích phân các phương trình vi phân dạng bảo toàn trên các thể tích kiểm soát hữu hạn (FVM) và tuyến tính hóa thành phần nguồn, phương trình trạng thái tại mỗi điểm nút quy về dạng quan hệ phụ thuộc tuyến tính giữa ba nút liền kề. Thuật toán TDMA (thuật toán Thomas) được sử dụng để giải hệ phương trình tuyến tính này một cách nhanh chóng qua hai bước khử tiến và thế lùi mà không cần nghịch đảo toàn bộ ma trận hệ số, giúp tối ưu hóa tốc độ tính toán trong Matlab.

3. Vì sao tốc độ sạc/xả cao (10C) lại gây sụt giảm điện áp làm việc nhanh chóng hơn so với 1C? Ở tốc độ 10C, dòng điện $I$ chạy qua pin rất lớn làm tăng độ sụt thế thuần trở $I \cdot R_f$ và gây ra quá thế phân cực nồng độ cao trong dung dịch điện phân do ion Lithium không kịp khuếch tán bù đắp giữa các vùng điện cực. Sự phân cực điện hóa mạnh mẽ này làm điện áp làm việc của pin suy giảm nhanh về ngưỡng cắt xả, khiến dung lượng khả dụng suy giảm so với chế độ xả chậm.

4. Điện thế cân bằng của cực âm và cực dương được xác định như thế nào trong mô hình? Điện thế cân bằng ($U_N$ và $U_P$) là hàm số phi tuyến thực nghiệm phụ thuộc vào các hệ số cân bằng hóa học stoichiometry ($x$ và $y$), phản ánh tỷ lệ giữa nồng độ ion Lithium thực tế trên bề mặt hạt hình cầu với nồng độ ion Lithium bão hòa tối đa trong cấu trúc mạng tinh thể của vật liệu điện cực tương ứng.

Kết luận

Đồ án tốt nghiệp đã xây dựng thành công mô hình điện hóa một chiều hoàn chỉnh và công cụ tính toán số trên nền tảng Matlab/Simulink nhằm mô phỏng các hiện tượng truyền khối, phân bố điện thế và động học phản ứng của pin Lithium-ion NMC 6Ah ứng dụng cho ô tô điện. Việc áp dụng chuẩn xác phương pháp thể tích hữu hạn kết hợp thuật toán giải ma trận ba đường chéo TDMA cho phép mô hình dự báo chính xác đặc tính điện áp làm việc ở nhiều chế độ sạc/xả (1C, 2C, 10C) và đã được kiểm chứng tương thích với các công trình khoa học quốc tế. Công trình cung cấp nền tảng lý thuyết và công cụ kỹ thuật thực nghiệm số đáng tin cậy, đóng góp thiết thực cho công tác nghiên cứu, thiết kế tối ưu hóa cấu trúc pin và phát triển hệ thống điều khiển quản lý năng lượng cho xe điện hiện nay.