Bộ giải FEM 3D nguồn mở cho điện từ (Maxwell hài thời gian)

Tôi đã tự hỏi liệu có tồn tại bất kỳ bộ giải FEM mã nguồn mở tốt (chính xác / nhanh / dễ sử dụng) nào cho các phương trình Maxwell hài hòa thời gian 3D không. Tôi đang tìm cách mô phỏng các hệ thống một vài bước sóng lớn trong kích thước X / Y và hàng chục bước sóng trong kích thước Z, vì vậy các vấn đề kích thước tương đối nhỏ. Tất cả các điện môi của tôi là đồng nhất và hình chữ nhật hoàn hảo (hộp).

Tôi đã thực hiện một số tìm kiếm trực tuyến và tìm thấy các gói FEM mã nguồn mở trưởng thành khác nhau, nhưng dường như chúng không hỗ trợ Maxwell. Bất kỳ đề xuất sẽ được nhiều đánh giá cao!

— Chi phí là
nguồn

Bạn đã thử FEniCS?

— facetus

Câu trả lời:

deal.II (xem http://www.dealii.org/ ) không hỗ trợ các yếu tố Nedelec và, do đó, có thể giải quyết các vấn đề bạn quan tâm. (Từ chối trách nhiệm đầy đủ: Tôi là một trong những nhà phát triển chính của thỏa thuận .II.)

— Wolfgang Bangerth
nguồn

Cảm ơn, giáo sư Bangerth! Liệu deal.II có bất kỳ hỗ trợ tích hợp nào cho các tài liệu PML hay sẽ cần được mã hóa riêng?

— Costis

deal.II không thực hiện bất kỳ phương trình hoặc công thức cụ thể nào - nó chỉ cung cấp các công cụ để làm như vậy. Vì vậy, trong trường hợp của bạn, bạn sẽ cần cung cấp mẫu song tuyến bạn muốn giải quyết. Cho dù điều đó bao gồm PML hay bất kỳ điều kiện ranh giới hấp thụ nào khác là tùy thuộc vào bạn.

— Wolfgang Bangerth

Liệu deal.ii có hỗ trợ số phức ngày hôm nay không? Tôi nghĩ nó hữu ích vì việc giải một hệ thống tuyến tính phức tạp sẽ tự nhiên hơn nhiều so với việc giải quyết hệ thống thực tương đương có kích thước gấp đôi.

— Hui Zhang

Không, nó vẫn sử dụng hai hệ thống có giá trị thực. Thật khó để đảm bảo rằng người ta phân biệt giữa transpose và Hermiti, các sản phẩm chấm thông thường và phức tạp, v.v.

— Wolfgang Bangerth

Hypre có một số điều kiện tiên quyết tích hợp để giải phương trình Maxwell. Có một số gói giao diện với nó (bạn có thể sử dụng hypre từ PETSc) như một bộ giải cho các hệ đại số tuyến tính, nhưng nó cũng có một lưới có cấu trúc và giao diện phần tử hữu hạn.

— Daniel Shapero
nguồn