1D FDTD mô phỏng trong C - ổn định vấn đề

[blah1]

Chào tất cả mọi người,

Tôi đã đọc về các FDTD khá một chút, cả hai đều thông qua các nguồn tài nguyên Tôi tìm thấy trên diễn đàn này và trong các thư viện địa phương.Tôi đã thực hiện một 1D FDTD thuật toán trong C, với một xung Gaussian như kích thích này.Mã này được đăng dưới đây.Nếu trên dòng 97 i đặt currentTime == 0, tôi nhận được các xung ban đầu như dự kiến.Tuy nhiên, ngay khi tôi cố gắng để có một bản chụp xung thêm xuống dòng thời gian này, nói rằng tại currentTime == 50, tôi có thể thấy rằng nó sẽ terribly không ổn định, chăm sóc, hướng tới vô cực.

Tôi nghĩ vấn đề của tôi là tôi không giới thiệu các nguồn xung một cách chính xác, nhưng cho cuộc sống của tôi không thể hiểu tại sao.

Bất kỳ manh mối để giúp tôi ra?

Cảm ơn!
Xin lỗi, nhưng bạn cần phải đăng nhập để xem tập tin đính kèm này

 

[blah1]

Ok, tôi đã nhận rằng vấn đề của sự mất ổn định được sắp xếp ra ngay bây giờ.Nó bật ra tôi đã quên một khung tại một trong những phương trình lặp.

Tuy nhiên, với cố định, tôi vẫn còn có một vấn đề.Làn sóng EM đang chuyển động quá "chậm".Tôi đã nhìn thấy giải pháp cho vấn đề này, và sau 50 thời gian các bước sóng nên chỉ ở phía trước của điện môi.Trong trường hợp của tôi, nó đã được phản ánh từ các chất điện môi mặc dù nó không phải đã đạt được chưa.Ngoài ra, việc tăng biên độ của sóng!Điều này, tất nhiên, là không thể kể từ khi làn sóng không thể đạt được năng lượng trong khi đi du lịch (trong trường hợp này ít nhất).Tôi thật sự khó khăn kể từ bây giờ tôi không biết những gì có thể làm điều này.

Nó sẽ là tuyệt vời nếu một đôi mắt có thể hãy xem, vì tôi tin rằng nó có thể là một sai lầm đơn giản một nơi nào đó, mà tôi chỉ cần giữ trên nhìn xuống.

Dù sao, mã mới được đính kèm, cũng như mô tả một áo sơ mi của vấn đề trong file doc..

Thanks in advance!

 

[drasko]

Quý vị có chăm sóc của các mối quan hệ giữa các bước không gian và bước thời gian?

FDTD là không ổn định nếu hai bước không thỏa mãn một điều kiện.

 

[blah1]

Ok, biến ra nó đã rối tung lên vì tôi xác định một trong những hằng số không chính xác.Mọi cố định bây giờ.Gửi bài mã, và mô tả các vấn đề trong tập tin đính kèm nếu có ai bao giờ cần nó.

Các mã này sẽ tạo một tập tin của các giá trị xy thích hợp cho âm mưu trong một cái gì đó như Gnuplot.

Chúc mừng!
Xin lỗi, nhưng bạn cần phải đăng nhập để xem tập tin đính kèm này