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

Sự chuyển dịch của ngành công nghiệp ô tô toàn cầu từ các dòng xe sử dụng động cơ đốt trong sang xe điện (EV) đặt pin Lithium-ion (Li-ion) vào vị trí công nghệ lưu trữ năng lượng cốt lõi. Hiệu năng, phạm vi hoạt động, mức độ an toàn và chi phí sản xuất của xe điện phụ thuộc trực tiếp vào các thành phần bên trong cell pin, trong đó lớp màng ngăn (separator) và dung dịch chất điện phân (electrolyte) đóng vai trò quyết định đến cơ chế truyền khối và dẫn truyền ion Lithium ($Li^+$). Để tối ưu hóa mật độ năng lượng, rút ngắn thời gian sạc, kéo dài tuổi thọ và giảm thiểu nguy cơ mất kiểm soát nhiệt, việc nghiên cứu thông qua mô hình hóa số và mô phỏng trên máy tính trở thành phương thức thiết yếu nhằm xác định các thông số vật lý và điện hóa tối ưu mà không đòi hỏi chi phí thử nghiệm thực nghiệm tốn kém.

Về tình hình nghiên cứu ngoài nước, các công trình quốc tế đã tập trung khai thác sâu các khía cạnh cấu trúc và hóa lý của cell pin:

  • Yifu Li, Zhongchao Tan và Wei Zhang (2021) khi phát triển mô hình số cho pin $LiFePO_4$ 38120 đã chỉ ra rằng việc tăng độ dày màng ngăn từ $5\ \mu\text{m}$ lên $100\ \mu\text{m}$ làm giảm mật độ năng lượng từ $148{,}8\text{ Wh/kg}$ xuống $110{,}6\text{ Wh/kg}$, giảm điện áp ổn định $0{,}02\text{ V}$ và dung lượng xả giảm từ $11{,}1\text{ Ah}$ xuống $6{,}3\text{ Ah}$, đồng thời xác định độ xốp tối ưu $80%$ giúp cải thiện truyền khối.
  • Dhrupad Parikh, Jeffrey Allen và Robert Kee (2019) so sánh hai loại màng thương mại Celgard 2325 và Celgard 2500, chứng minh Celgard 2500 có điện trở thấp hơn ($2{,}23\ \Omega/\text{cm}^2$) và khả năng thấm ướt tốt hơn, nâng cao dung lượng thêm $57%$ ở tốc độ $2\text{C}$ và $47%$ ở tốc độ $3\text{C}$.
  • Sebastian Klick, Martin Schulze và Hans Müller (2020) chứng minh pin có lượng chất điện phân thấp ($2{,}1\text{ mL}$, hệ số lấp đầy $0{,}9$) bị suy giảm dung lượng từ $60–85%$ sau 4 tuần do hiện tượng khô cạn điện phân làm gia tăng trở kháng ion, so với mức suy giảm chỉ $10%$ ở lượng chất điện phân $3{,}5\text{ mL}$.
  • Rahifa Ranom, Ahmad Zaki và Siti Noor (2022) ứng dụng phương pháp đường thẳng (MOL) trên pin $LiFePO_4$, ghi nhận độ khuếch tán ion Lithium đạt $6 \times 10^{-11}\text{ m}^2/\text{s}$ làm tăng hiệu suất phóng điện thêm $30%$, đồng thời khi tăng nồng độ ion ban đầu từ $800\text{ mol/m}^3$ lên $1200\text{ mol/m}^3$ thì hiệu suất tăng thêm $14%$.
  • Các nghiên cứu khác như Chris Manzie et al. (2021) phát triển mô hình hạt đơn từ 17 phương trình đạo hàm riêng; Manuel Stephan (2006) khảo sát chất điện phân polymer gel PVdF-HFP đạt độ dẫn ion $10^{-3}\text{ S/cm}$; Costa et al. (2013) chế tạo màng PVDF-HFP/MgO đạt độ dẫn $5 \times 10^{-2}\text{ S/cm}$; Hun Lee et al. (2014) và Meltem Yanilmaz et al. (2014) phát triển các dạng màng không dệt và màng lai nhằm tối ưu hóa khả năng hấp thụ dung dịch và độ ổn định nhiệt.

Tại Việt Nam, nghiên cứu về pin cho xe điện đang được đẩy mạnh nhằm đáp ứng điều kiện vận hành nhiệt đới:

  • Đỗ Nguyễn Trung Kiên và Nguyễn Văn Dũng (2025) tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh đã xây dựng mô hình điện hóa một chiều kết hợp phần mềm MATLAB/Simulink và phương pháp thể tích hữu hạn để phân tích điện áp, điện thế và nồng độ ion theo thời gian cũng như tốc độ sạc/xả.
  • Nguyễn Văn Hùng, Trần Thị Minh và Lê Văn Tuấn (2023) tại Đại học Bách khoa Hà Nội khảo sát màng ngăn polyolefin trong pin $LiFePO_4$, xác định màng ngăn dày $20\ \mu\text{m}$ với độ xốp $50%$ giúp cải thiện hiệu suất phóng điện $12%$ và hạn chế quá nhiệt ở dải nhiệt độ $35–45^\circ\text{C}$.
  • Phạm Thị Lan, Hoàng Văn Nam và Nguyễn Thị Hồng (2024) tại Đại học Khoa học Tự nhiên TP.HCM nghiên cứu dung dịch điện phân lỏng $LiPF_6$ (độ dẫn ion $1–2\text{ S/cm}$) duy trì tính ổn định ở $45^\circ\text{C}$, hỗ trợ trạm sạc xe điện tích hợp năng lượng mặt trời.

Trên cơ sở đó, đồ án tốt nghiệp của sinh viên Phùng Hoàng Nam và Nguyễn Minh Trí tập trung giải quyết bài toán mô phỏng ảnh hưởng của màng ngăn và nồng độ chất điện phân đến hoạt động của pin Lithium Graphite.

Nhiệm vụ nghiên cứu của đề tài gồm 3 nội dung:

  1. Tìm hiểu lý thuyết cấu tạo của pin, quá trình hoạt động nạp, xả của pin Lithium.
  2. Nghiên cứu, phân tích và tính toán kích thước màng ngăn, nồng độ chất điện phân ảnh hưởng tới quá trình hoạt động của pin Lithium Graphite.
  3. Mô phỏng, phân tích, dự đoán sự thay đổi của quá trình hoạt động của pin sử dụng phần mềm MATLAB/Simulink.

Mục tiêu nghiên cứu cụ thể được xác định gồm:

  1. Nghiên cứu, tìm hiểu và nắm rõ cấu tạo, nguyên lý hoạt động của pin Lithium Graphite.
  2. Nghiên cứu, phân tích sự thay đổi các thông số của pin Lithium Graphite trong việc thay đổi các thông số màng ngăn và nồng độ chất điện phân qua quá trình hoạt động của pin.
  3. Mô phỏng và so sánh kết quả với kết quả ban đầu của một viên pin sau khi đã thay đổi thông số như nồng độ chất điện phân và độ dày màng ngăn trong quá trình mô phỏng.
  4. Nhận xét và đưa ra thông số tối ưu để pin có hiệu suất tốt nhất.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Pin Lithium Graphite (mẫu cell pin Lithium-ion NMC, dung lượng định mức $6\text{ Ah}$).
  • Phạm vi nghiên cứu: Khảo sát quá trình hoạt động của pin Lithium Graphite bằng phần mềm mô phỏng MATLAB/Simulink trong điều kiện vận hành lý tưởng.

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

Đồ án ứng dụng mô hình điện hóa tựa hai chiều P2D (Pseudo-two-dimensional) được thiết lập trên hệ trục tọa độ Descartes. Không gian tính toán của cell pin được chia thành ba miền liên tục: điện cực âm (anode than chì - Graphite, độ dày $L_N$), lớp màng ngăn cách điện (separator, độ dày $L_S$), và điện cực dương (cathode NMC, độ dày $L_P$).

Mô hình nhiệt động học và động học điện hóa dựa trên các hệ phương trình cơ bản:

  • Động học phản ứng trao đổi bề mặt Butler-Volmer: Mô tả mối quan hệ giữa mật độ dòng điện chuyển dịch $j_{Li}$ với quá thế kích hoạt $\eta$: $$j_{Li} = i_0 \left[ \exp\left(\frac{\alpha_a F \eta}{R T}\right) - \exp\left(-\frac{\alpha_c F \eta}{R T}\right) \right]$$ Trong đó: $i_0$ là mật độ dòng trao đổi tại bề mặt hạt; $\alpha_a, \alpha_c$ là các hệ số chuyển giao điện tích; $F$ là hằng số Faraday ($96485\text{ C/mol}$); $R$ là hằng số khí lý tưởng; $T$ là nhiệt độ tuyệt đối; $\eta = \varphi_s - \varphi_e - U$ với $\varphi_s$ là điện thế pha rắn, $\varphi_e$ là điện thế pha lỏng, và $U$ là điện thế cân bằng nhiệt động học phụ thuộc vào tỷ lệ nồng độ Lithium ở điện cực âm ($U_N$) và cực dương ($U_P$).
  • Phương trình bảo toàn điện tích ở pha rắn (điện cực) và pha lỏng (chất điện phân): Xác định sự phân bố điện thế $\varphi_s$ trong chất điện cực và $\varphi_e$ trong dung dịch điện môi.
  • Quan hệ cấu trúc vi mô theo định luật Bruggeman: Độ dẫn ion hiệu quả $\kappa_{eff}$ và hệ số khuếch tán hiệu quả $D_e^{eff}$ trong môi trường xốp được hiệu chỉnh thông qua độ xốp $\varepsilon_e$ và số mũ Bruggeman $\gamma$: $$\kappa_{eff} = \kappa \cdot \varepsilon_e^\gamma, \quad D_e^{eff} = D_e \cdot \varepsilon_e^\gamma$$
  • Xác định nội trở màng ngăn ($R_{sep}$) và nội trở dung dịch điện phân ($R_{electrolyte}$): Dựa trên kích thước hình học (độ dày $L_S, L_N, L_P$, diện tích điện cực $A$) và các giá trị độ dẫn ion hiệu quả tương ứng ($\kappa_{eff,anode}, \kappa_{eff,sep}, \kappa_{eff,cathode}$).
  • Mô hình sinh nhiệt và biến thiên nhiệt độ: Xác định tổn hao nhiệt dựa trên hiệu ứng Joule của nội trở tổng $R_{Internal}$ kết hợp cùng quá thế kích hoạt điện hóa khi có dòng tải $I$ chạy qua.

Về phương pháp giải số, đề tài sử dụng phương pháp thể tích hữu hạn (Finite Volume Method - FVM):

  • Miền tính toán một chiều dọc theo bề dày cell pin được chia thành các thể tích kiểm soát con (control volumes) với các điểm nút lưới có khoảng cách $\Delta x$ và $\delta x$. Điện cực âm được rời rạc hóa thành $N_n$ miền thể tích con, điện cực dương gồm $N_p$ miền thể tích con, xen giữa là miền màng ngăn.
  • Các phương trình vi phân từng phần phi tuyến mô tả điện thế pha rắn $\varphi_s$ được tuyến tính hóa thành phần nguồn ($S$).
  • Thuật toán giải ma trận ba đường chéo (Tri-Diagonal Matrix Algorithm - TDMA) được lập trình để giải hệ phương trình đại số tuyến tính sau rời rạc hóa, cho phép xác định chính xác giá trị $\varphi_s$ và $\varphi_e$ tại từng điểm nút lưới theo các bước thời gian $t$.
  • Toàn bộ thuật toán FVM-TDMA và các phương trình điện hóa được cấu hình và chạy mô phỏng trên nền tảng phần mềm MATLAB/Simulink. Nguồn dữ liệu thông số đầu vào được trích xuất từ mẫu cell thực tế Lithium-ion NMC dung lượng $6\text{ Ah}$.

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

Chương 1: Tổng quan

Chương này trình bày lý do lựa chọn đề tài xuất phát từ nhu cầu tối ưu hóa hệ thống pin cho xe điện trong quá trình chuyển đổi giao thông xanh. Tác giả tổng thuật 20 công trình quốc tế tiêu biểu về cấu trúc màng cách điện, hệ dung dịch điện phân lỏng/gel, và các nghiên cứu mô hình hóa trong nước (Đỗ Nguyễn Trung Kiên & Nguyễn Văn Dũng, Nguyễn Văn Hùng et al., Phạm Thị Lan et al.). Từ đó, chương xác lập rõ ràng 4 mục tiêu, 3 nhiệm vụ đề tài, xác định đối tượng nghiên cứu là pin Lithium Graphite và giới hạn phạm vi mô phỏng số trên phần mềm MATLAB/Simulink.

Chương 2: Tổng quan về pin điện Lithium trên ô tô

Chương 2 cung cấp hệ thống kiến thức toàn diện về công nghệ pin lưu trữ trên xe điện:

  • Lịch sử phát triển: Từ phát minh pin axit-chì của Gaston Planté (1859), chiếc xe điện 3 bánh của Siemens hợp tác cùng James Starley (1880), xe điện Thomas Parker (1884), kỷ lục tốc độ $105{,}88\text{ km/h}$ của Camille Jenatzy trên xe Jamais Contente (1899), cho đến khi Sony thương mại hóa thành công pin Lithium-ion vào năm 1991. So sánh hiệu suất chỉ ra pin Lithium-ion duy trì sự ổn định dung lượng vượt trội so với sự sụt giảm nhanh của pin axit-chì.
  • Cấu tạo và nguyên lý: Mô tả chi tiết các thành phần chính (anode Graphite, cathode, màng ngăn, chất điện phân), cơ chế xen kẽ và tách chiết của ion $Li^+$ giữa hai điện cực trong hai chu trình sạc và xả.
  • Ưu nhược điểm và công tác quản lý an toàn: Đánh giá mật độ năng lượng, tuổi thọ chu kỳ, các nguy cơ thoát nhiệt; vai trò của hệ thống quản lý pin (BMS) trong việc giám sát điện áp theo ngưỡng an toàn, kiểm soát nhiệt độ vận hành và duy trì cân bằng điện áp giữa các cell pin.
  • Phân tích các thông số kỹ thuật cốt lõi: Dung lượng cell $Q$ ($6\text{ Ah}$), năng lượng dự trữ $E$ ($Wh$), điện áp hở mạch $U_{OCV}$, điện áp làm việc $U_{working}$, nội trở bên trong $R_{Internal}$, tốc độ dòng nạp/xả C-rate, công suất tức thời $P$, và chu kỳ vòng đời đi kèm hệ số suy hao $\delta$.

Chương 3: Cơ sở lý thuyết

Chương 3 thiết lập nền tảng toán học và vật lý mô tả các quá trình bên trong cell pin:

  • Quy ước phân bố nồng độ ion Lithium: Xác định nồng độ trong pha rắn $C_s$, nồng độ tại bề mặt hạt cầu $C_{s,surface}$, và nồng độ cực đại lý thuyết $C_{s,max}$.
  • Hiệu suất dẫn truyền ion: Áp dụng hệ số khuếch tán $D_S, D_E$ và công thức Bruggeman để tính toán độ dẫn ion hiệu quả $\kappa_{eff}$ qua môi trường xốp của màng ngăn và điện cực.
  • Nội trở vật lý: Xây dựng công thức tính nội trở riêng biệt của màng ngăn $R_{sep}$ và chất điện phân $R_{electrolyte}$.
  • Cân bằng nhiệt: Thiết lập phương trình cân bằng năng lượng để tính toán tổng nhiệt lượng sinh ra do tổn hao thuần trở Joule kết hợp phản ứng hóa học, từ đó dự báo mức độ tăng nhiệt của cell.
  • Các phương trình điện hóa: Trình bày phương trình bảo toàn điện tích, phương trình động học điện hóa Butler-Volmer và phương trình tính thế cân bằng nhiệt động học $U$ tại cực âm ($U_N$) và cực dương ($U_P$).

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

Chương 4 mô tả quy trình hiện thực hóa mô hình toán học thành công cụ mô phỏng số:

  • Áp dụng phương pháp thể tích hữu hạn (FVM): Thiết lập hệ lưới không gian trên tọa độ Descartes, chia cực âm thành $N_n$ phân đoạn và cực dương thành $N_p$ phân đoạn, gắn các biến trạng thái vào các điểm nút lưới.
  • Xử lý hệ phương trình: Tuyến tính hóa thành phần nguồn $S$ và triển khai thuật toán giải ma trận ba đường chéo (TDMA) cho phương trình vi phân từng phần để tìm nghiệm điện thế pha rắn $\varphi_s$ theo tọa độ và thời gian $t$.
  • Thiết kế kịch bản khảo sát:
    • Lựa chọn và tính toán thay đổi thông số độ dày màng ngăn $L_S$ ở 3 mức: $10\ \mu\text{m}$, $20\ \mu\text{m}$, và $30\ \mu\text{m}$.
    • Lựa chọn và tính toán thay đổi nồng độ chất điện phân ban đầu $C_e$ ở 3 mức: $0{,}5\text{ M}$, $1{,}2\text{ M}$, và $2\text{ M}$.
    • Thiết lập chương trình mô phỏng sự biến thiên của $\varphi_s$ và điện áp làm việc $U_{working}$ trong suốt các chu kỳ nạp và xả tương ứng.

Chương 5: Kết quả mô phỏng

Chương này công bố dữ liệu mô phỏng thu được trên mẫu pin Lithium-ion NMC dung lượng $6\text{ Ah}$:

  • Điện áp làm việc ($U_{working}$): Đưa ra các bảng số liệu chi tiết và các đường cong so sánh điện áp cell khi thay đổi độ dày màng ngăn ($10, 20, 30\ \mu\text{m}$) và thay đổi nồng độ chất điện phân ($0{,}5; 1{,}2; 2\text{ M}$) ở cả hai quá trình sạc và xả.
  • Điện thế chất điện cực ($\varphi_s$): Đồ thị và số liệu phân bố điện thế $\varphi_s$ tại các điểm nút dọc theo bề dày cực âm và cực dương ứng với từng trường hợp $L_s$ ($10, 20, 30\ \mu\text{m}$) và $C_e$ ($0{,}5; 1{,}2; 2\text{ M}$).
  • Nhiệt học cell pin: Phân tích bảng tính toán nhiệt lượng tỏa ra và mức độ gia tăng nhiệt độ của cell pin tương ứng với từng điều kiện màng ngăn và nồng độ chất điện phân.

Chương 6: Kết luận

Chương cuối cùng tổng kết các kết quả đạt được của đề tài, đánh giá mức độ tương thích giữa mô hình mô phỏng số FVM-TDMA trên MATLAB/Simulink với lý thuyết động học điện hóa, đồng thời đưa ra các nhận định về thông số tối ưu cho pin Lithium Graphite.

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

Đồ án đã lượng hóa ảnh hưởng của hai thông số cấu trúc (độ dày màng ngăn $L_S$ và nồng độ dung dịch chất điện phân $C_e$) đến hoạt động của cell pin Lithium NMC dung lượng $6\text{ Ah}$ thông qua các dữ liệu mô phỏng:

Kịch bản khảo sát Giá trị khảo sát Ảnh hưởng đến $U_{working}$ và nội trở Ảnh hưởng đến phân bố điện thế $\varphi_s$ và nhiệt sinh
Độ dày màng ngăn ($L_S$) $10\ \mu\text{m}$ Giảm nội trở màng ngăn $R_{sep}$, giảm sụt áp dung dịch, duy trì mức điện áp làm việc $U_{working}$ cao hơn khi xả. Phân bố điện thế $\varphi_s$ đồng đều hơn do khoảng cách giữa hai điện cực thu hẹp; lượng nhiệt sinh do trở kháng giảm.
$20\ \mu\text{m}$ (Giá trị gốc) Đảm bảo tính ổn định cơ bản về nội trở và điện áp định danh của cell chuẩn $6\text{ Ah}$. Phân bố thế ổn định, cân bằng giữa tốc độ chuyển dời ion và khoảng cách an toàn điện học.
$30\ \mu\text{m}$ Tăng nội trở $R_{sep}$, cản trở tốc độ di chuyển của ion $Li^+$, gây sụt áp làm việc lớn hơn trong chu trình xả. Độ dốc phân bố $\varphi_s$ tăng, gia tăng tổn hao Joule thuần trở dẫn đến tăng nhiệt lượng tỏa ra của cell pin.
Nồng độ chất điện phân ($C_e$) $0{,}5\text{ M}$ Thiếu hụt lượng hạt mang điện dẫn đến giảm độ dẫn ion hiệu quả $\kappa_{eff}$, làm tăng nội trở $R_{electrolyte}$, suy giảm điện áp làm việc. Tăng độ phân cực điện hóa, chênh lệch điện thế $\varphi_s$ tại bề mặt điện cực lớn, tốc độ sinh nhiệt tăng cao.
$1{,}2\text{ M}$ (Giá trị chuẩn) Đạt độ dẫn ion hiệu quả tối ưu cho hệ điện cực Graphite/NMC, duy trì đường cong phóng/nạp ổn định. Đường phân bố $\varphi_s$ tại cả cực âm và cực dương duy trì gradient ổn định, hạn chế suy thoái nhiệt.
$2\text{ M}$ Nồng độ ion Lithium cao hỗ trợ quá trình khuếch tán ở tốc độ lớn, tuy nhiên nồng độ tăng cao có thể làm tăng độ nhớt dung dịch và ảnh hưởng đến hệ số khuếch tán $D_e^{eff}$. Phân bố $\varphi_s$ tập trung tại vùng sát bề mặt tiếp xúc màng ngăn; nhiệt độ cell duy trì ở mức kiểm soát được.

Đóng góp khoa học và kỹ thuật của đồ án:

  • Xây dựng hoàn chỉnh mô hình tính toán số áp dụng phương pháp thể tích hữu hạn (FVM) kết hợp thuật toán TDMA để giải hệ phương trình vi phân từng phần phi tuyến mô tả phân bố điện thế $\varphi_s$ của pin Lithium trên phần mềm MATLAB/Simulink.
  • Cung cấp bộ dữ liệu so sánh định lượng về phân bố điện thế, điện áp làm việc và các giá trị nhiệt lượng tỏa ra ứng với các biến thiên thực tế của màng ngăn và chất điện phân, làm cơ sở khoa học để lựa chọn thông số kết cấu phù hợp trong thiết kế cell pin xe điện.

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

Hạn chế của nghiên cứu:

  • Giới hạn mô phỏng ở điều kiện lý tưởng: Đề tài tiến hành khảo sát mô hình khi pin hoạt động trong trạng thái lý tưởng, chưa tính đến ảnh hưởng của các yếu tố ngoại cảnh khắc nghiệt như dải nhiệt độ môi trường thay đổi đột ngột hoặc rung động cơ học trên xe.
  • Giả định mô hình một chiều: Miền tính toán FVM được xây dựng dọc theo một trục tọa độ Descartes theo phương bề dày cell pin, chưa phản ánh đầy đủ hiệu ứng gradient nhiệt và mật độ dòng điện phân bố không đồng đều theo ba chiều không gian (3D) của cell thực tế.
  • Thiếu kiểm chứng thực nghiệm vật lý: Do hạn chế về thời gian thực hiện đồ án tốt nghiệp, toàn bộ kết quả dừng lại ở mức mô phỏng số trên phần mềm MATLAB/Simulink mà chưa được kiểm chứng đối sánh trực tiếp với cell pin vật lý thực tế trên các băng thử nghiệm sạc/xả.

Hướng nghiên cứu tiếp:

  • Mở rộng sang mô hình nhiệt - điện hóa ghép nối đa chiều (coupled 3D electro-thermal modeling) để đánh giá chính xác trường nhiệt độ trong các cell pin cỡ lớn.
  • Tích hợp các cơ chế thoái hóa dài hạn như quá trình hình thành và phát triển của lớp màng thụ động bề mặt (SEI) và hiện tượng mạ Lithium (Lithium plating) qua hàng ngàn chu kỳ sạc nhanh.
  • Thiết lập quy trình đo đạc thực nghiệm trên cell pin thực tế trong phòng thí nghiệm để tinh chỉnh và kiểm chứng các hệ số động học của mô hình.

Giá trị tham khảo

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

  • Sinh viên, học viên và giảng viên chuyên ngành Công nghệ Kỹ thuật Ô tô, Kỹ thuật Cơ khí Động lực, Kỹ thuật Năng lượng và Kỹ thuật Điều khiển - Tự động hóa đang nghiên cứu về hệ thống truyền động điện và công nghệ lưu trữ năng lượng trên ô tô điện.
  • Kỹ sư phát triển phần mềm mô phỏng hệ thống pin, kỹ sư thiết kế hệ thống quản lý pin (BMS) cần thuật toán và cơ sở toán học để ước lượng trạng thái bên trong cell.

Các nội dung có giá trị tham khảo chuyên môn cao nhất gồm:

  • Toàn bộ nội dung Chương 3 và Chương 4: Cung cấp phương pháp chi tiết về cách rời rạc hóa các phương trình điện hóa Butler-Volmer và phương trình bảo toàn điện tích bằng phương pháp thể tích hữu hạn (FVM) kết hợp thuật toán TDMA trên môi trường MATLAB/Simulink.
  • Bảng số liệu và đồ thị trong Chương 5: Cung cấp nguồn thông số đối chiếu trực tiếp về biến thiên điện áp và phân bố điện thế $\varphi_s$ khi thay đổi kích thước màng ngăn ($10, 20, 30\ \mu\text{m}$) và nồng độ dung dịch ($0{,}5; 1{,}2; 2\text{ M}$) trên cell pin NMC dung lượng $6\text{ Ah}$.

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

1. Đề tài khảo sát mẫu cell pin cụ thể nào và có dung lượng định mức bao nhiêu?
Đề tài thực hiện khảo sát và mô phỏng trên mẫu cell pin Lithium-ion sử dụng cực âm than chì (Graphite) kết hợp cực dương NMC (Lithium Nickel Manganese Cobalt Oxide), có dung lượng định mức là $6\text{ Ah}$.

2. Độ dày màng ngăn và nồng độ chất điện phân được lựa chọn ở những mức giá trị nào để khảo sát?
Trong Chương 4 và Chương 5, độ dày màng ngăn $L_S$ được khảo sát tại 3 giá trị: $10\ \mu\text{m}$, $20\ \mu\text{m}$ (thông số chuẩn), và $30\ \mu\text{m}$. Nồng độ chất điện phân $C_e$ được khảo sát tại 3 mức: $0{,}5\text{ M}$, $1{,}2\text{ M}$ (thông số chuẩn), và $2\text{ M}$.

3. Phương pháp toán học nào được nhóm tác giả sử dụng để giải các phương trình vi phân trong mô hình?
Tác giả sử dụng phương pháp thể tích hữu hạn (FVM) để rời rạc hóa không gian miền tính toán trên hệ tọa độ Descartes, kết hợp tuyến tính hóa thành phần nguồn ($S$) và ứng dụng thuật toán ma trận ba đường chéo (TDMA) để giải hệ phương trình vi phân từng phần xác định điện thế pha rắn $\varphi_s$.

4. Đại lượng nào được tập trung tính toán để đánh giá trạng thái điện học tại các điện cực trong quá trình sạc và xả?
Đại lượng trọng tâm được tính toán và vẽ đồ thị phân bố theo không gian và thời gian là điện thế tại chất điện cực $\varphi_s$, được theo dõi độc lập tại miền thể tích cực âm ($N_n$ miền con) và cực dương ($N_p$ miền con), bên cạnh thông số điện áp làm việc tổng thể của cell ($U_{working}$).

5. Động học phản ứng tại bề mặt tiếp xúc giữa hạt điện cực và chất điện phân được mô tả bằng phương trình nào?
Động học phản ứng trao đổi ion được mô tả thông qua phương trình Butler-Volmer, xác định mật độ dòng chuyển dời $j_{Li}$ phụ thuộc vào mật độ dòng trao đổi $i_0$, các hệ số chuyển giao điện tích $\alpha_a, \alpha_c$, và quá thế kích hoạt $\eta = \varphi_s - \varphi_e - U$.

Kết luận

Đồ án tốt nghiệp của sinh viên Phùng Hoàng Nam và Nguyễn Minh Trí đã xây dựng thành công mô hình toán học và công cụ mô phỏng số trên nền tảng MATLAB/Simulink nhằm phân tích hoạt động của cell pin Lithium-ion Graphite/NMC dung lượng $6\text{ Ah}$. Thông qua phương pháp thể tích hữu hạn kết hợp thuật toán TDMA, đề tài đã làm rõ mối quan hệ giữa độ dày màng ngăn ($10–30\ \mu\text{m}$), nồng độ dung dịch điện phân ($0{,}5–2\text{ M}$) đối với điện áp làm việc, phân bố điện thế pha rắn và lượng nhiệt sinh của cell. Kết quả nghiên cứu cung cấp cơ sở dữ liệu học thuật và kỹ thuật tin cậy, đóng góp vào việc lựa chọn thông số kết cấu tối ưu trong công nghệ thiết kế pin phục vụ phương tiện giao thông chạy điện.