TL;DR

Lucas Åkerstedt, Harald Hultin (KTH & Saab Surveillance) và B. L. G. Jonsson vừa công bố trên IEEE Transactions on Antennas and Propagation số tháng 2/2026 một phương pháp Array Decomposition cho Method of Moments (MoM) dành riêng cho mảng anten hữu hạn có phần tử nối điện. Ý tưởng cốt lõi: biểu diễn toàn bộ mảng $N_x \times N_y$ bằng đúng 9 thành phần (1 phần tử trung tâm + 8 lề). Ma trận trở kháng thu được có cấu trúc multilevel block Toeplitz, nhờ đó bộ nhớ giảm từ $\mathcal{O}(N_x^2 N_y^2)$ xuống $\mathcal{O}(N_x N_y)$. Một mảng BoR 32×32 chỉ cần 3.97 GB thay vì 1.97 TB — giảm gần 500 lần.

Có gì mới

Mô phỏng full-wave mảng anten lớn bằng MoM cổ điển tốn bộ nhớ $\mathcal{O}(N^2)$ — với hàng nghìn phần tử, ngay cả cụm máy chủ cao cấp cũng đầu hàng. Đã có nhiều hướng tăng tốc (MLFMA, CBF, fast multipole), nhưng với mảng có phần tử nối điện — loại cấu hình quan trọng cho anten dải rộng, Vivaldi, tightly coupled — chọn lựa rất hạn chế. Phương pháp Higher-Order Array Decomposition Method (HO-ADM) của TICRA làm được, nhưng dùng Discontinuous Galerkin với basis functions bậc cao, không tương thích với các solver MoM dùng hàm cơ sở RWG chuẩn.

Bài báo mới lấp đúng khoảng trống đó: lần đầu có cách xây ma trận MoM multilevel block Toeplitz cho mảng nối điện mà vẫn dùng RWG chuẩn. Toolchain sẵn có — mesh, solver thương mại — có thể tận dụng trực tiếp.

Tại sao quan trọng

Các thế hệ mạng di động kế tiếp (6G) được kỳ vọng dựa trên mảng anten có hơn 1,000 phần tử để beam-forming / beam-focusing trong vùng trường gần. Không có công cụ mô phỏng đủ tầm, việc thiết kế phần cứng trở thành nút cổ chai. Một công cụ MoM chạy trên máy tính kỹ sư bình thường — không phải HPC — mở ra vòng lặp thiết kế nhanh hơn nhiều. Đồng tác giả từ Saab Surveillance cho thấy giá trị trực tiếp cho radar và aerospace; nghiên cứu được tài trợ một phần bởi chương trình SEE-6GIA của Swedish Research Council.

Số liệu kỹ thuật

BenchmarkPhương pháp mớiMoM trực tiếpTỉ lệ giảm
Ma trận trở kháng mảng BoR 32×32 (350,519 cạnh RWG)3.97 GB1.97 TB~496×
Hai mảng 9×9 ghép chặt (157,220 cạnh RWG)9.53 GB395 GB~41×
Scaling tiệm cận$\mathcal{O}(N_x N_y)$$\mathcal{O}(N_x^2 N_y^2)$

Với mảng 32×32 BoR: directivity broadside đạt 35 dBi; Total Active Reflection Coefficient (TARC) giữ dưới −10 dB trong phạm vi quét ±52.9° (E-plane)±52.5° (H-plane). Với cặp mảng 9×9 chạy băng 3–9 GHz: 36 điểm tần số tính chính xác bằng thuật toán Theta I, nội suy Loewner ra 500 điểm cách đều, sai số tương đối tối đa $10^{-2}$.

Thủ thuật 9 thành phần

Thay vì mesh toàn bộ mảng (vốn tạo ra hàng trăm nghìn cạnh RWG), phương pháp chỉ mesh một phần tử trung tâmtám thành phần lề (4 cạnh + 4 góc) đủ để dán biên. Toàn mảng được tái lập bằng dịch chuyển (translation) 9 thành phần này trên lưới. Thuật toán phân hoạch xác định các shared triangles tại biên giữa hai phần kề nhau, đảm bảo dòng mặt liên tục. Vì hàm Green trong MoM bất biến theo dịch chuyển, hai cặp phần cách nhau cùng một khoảng không gian cho ra cùng một submatrix — vậy là phần lớn ma trận trở kháng được tái sử dụng thay vì tính lại. Đây là nguồn gốc của cấu trúc block Toeplitz đa lớp, và cũng là nguồn gốc của khoản tiết kiệm bộ nhớ một bậc độ lớn.

So sánh với các phương pháp hiện có

Phương phápNối điện?BasisToeplitz đa lớp?RWG chuẩn?
MoM cổ điểnRWGKhôngCó (nhưng không scale)
MLFMA / Fast MultipoleRWGKhông
HO-ADM (TICRA)Có (qua DGM)Bậc caoKhông
Åkerstedt et al. 2026RWG

Use cases

  • Hạ tầng 6G: mô phỏng mảng massive MIMO trên 1,000 phần tử trên workstation thường.
  • Radar & aerospace: phân tích anten radar pha mảng với phần tử nối điện dải rộng.
  • Ghép cặp mảng dị chủng: bài báo minh hoạ mutual coupling giữa hai mảng 9×9 với phần tử khác nhau (T-slot dipole cạnh Vivaldi) cách nhau 225.5 mm — một kịch bản thực tế khi nhiều anten cùng tồn tại trên một platform.
  • Thiết kế & debug: vẽ dòng mặt liên tục toàn mảng, cực kỳ hữu ích khi chẩn đoán khớp biên và cộng hưởng.

Hạn chế & điều kiện

  • Yêu cầu đối xứng dịch chuyển hữu hạn — mảng phải nằm trên lưới Cartesian đều. Mảng bất quy tắc không áp dụng được.
  • Mesh tại biên nội bộ giữa các thành phần phải khớp tuần hoàn trên cả hai hướng.
  • Cấu trúc 9 thành phần (1 center + 8 margin) là bắt buộc — topology phức tạp hơn cần mở rộng thuật toán.
  • Preprint mở trên arXiv:2506.04710; bản IEEE TAP có phí.

Tiếp theo là gì

Bài báo hé lộ khả năng mở rộng 3D (sort shared triangles theo toạ độ z), mở đường cho mảng conformal và stack nhiều tầng. Kết hợp với HO-ADM (basis bậc cao) có thể cho ra solver vừa tiết kiệm bộ nhớ vừa giảm số bậc tự do. Hướng nghiên cứu tự nhiên: GPU-accelerated Toeplitz solver, differentiable simulation (TICRA cũng đang đi hướng này), và tích hợp vào các solver thương mại để kỹ sư anten dùng trực tiếp.

Nguồn: arXiv preprint, IEEE TAP vol 74 no 2, Feb 2026, TICRA HO-ADM.