[Υπολογιστικός Ηλεκτρομαγνητισμός] Γενικές απορίες και επικαιρότητα 2020

[Thmmy.gr portal]

  Thmmy.gr portal
   Forum
   Downloads
   Ενεργ. Λογαριασμού
   Επικοινωνία
  
  Χρήσιμα links
   Σελίδα τμήματος
   Βιβλιοθήκη Τμήματος
   Φοιτητικά fora
   Πρόγραμμα Λέσχης
   Πρακτική Άσκηση
   Ηλεκτρονική Εξυπηρέτηση Φοιτητών
   Διανομή Συγγραμμάτων
   Ψηφιακό Καταθετήριο Διπλωματικών
   Πληροφορίες Καθηγητών
   Αντιστοίχηση Μαθημάτων ΝΠΣ με ΠΠΣ
   mTHMMY
  
  Φοιτητικές Ομάδες
   IEEE φοιτητικό παράρτημα ΑΠΘ
   Ομάδα Ρομποτικής
   BEST Thessaloniki
   IAESTE Thessaloniki
   Aristurtle
   Ariadne
   EESTEC LC Thessaloniki
  

THMMY.gr Web

Helios_MultiScribbles2SMF Default ThemeSMFone_Blue

[12:00-13:30]

Η γραμματεία είναι ανοιχτή καθημερινά 12:00-13:30

[gkokk]

Μάλιστα ναι δεν έκανα την σύνδεση, σ'ευχαριστώ που μου το είπες κιόλας να το ξέρω απο δω και πέρα.

[dtsaravas]

Καλησπέρα!

Σε σχέση με την εργασία φαίνεται να έχω κάνει σωστά την υλοποίηση του main grid. Όμως, όσες φορές και αν είδα την υλοποίηση του pml δεν μπορώ να βρω το λάθος μου. Συγκεκριμένα φαίνεται να προκύπτουν τεράστιες τιμές κοντά στο άπειρο (10^31) κτλ και το τρέξιμο σταματάει λόγω των πολύ μεγάλων αυτών τιμών όταν το κύμα φτάνει στο pml. Πατώντας να δω τον πίνακα με όλες τις τιμές εμφανίζονται ΝαΝ στα κελια του pml. Συνάντησε κανείς κάτι παραπλήσιο; Έχω τσεκάρει την υλοποίηση αρκετές φορές διορθώνοντας μερικά πράγματα, αλλά πλέον δεν πάει το μυαλό μου τι άλλο μπορεί να φταίει...
 

[lhmma_Jordan]

se=-e0*c*log(R)/(2^(O+2)*dx*Npml^(O+1));

από τις οδηγίες για το PML. Από που προκύπτει ο τύπος αυτός; Δεν μπόρεσα να βρω κάτι παρόμοιο στις σημειώσεις.

[ChampionsNet]

Πως μπορούμε να αξιολογήσουμε την ικανότητα απορρόφησης και πως μπορούμε να υπολογίσουμε τη φάση του κύματος σε κάθε σημείο;

[ChampionsNet]

Μέτρησα τη διαφορά φάσης αναπαριστώντας το Εz για ένα συγκεκριμένο y, μετρώντας τη διαφορά σε dx ανάμεσα σε δύο περιόδους και μου βγαίνει ταχύτητα διάδοσης 1.5c 

Τι πάει τόσο στραβά;

[lhmma_Jordan]

Μέτρησα τη διαφορά φάσης αναπαριστώντας το Εz για ένα συγκεκριμένο y, μετρώντας τη διαφορά σε dx ανάμεσα σε δύο περιόδους και μου βγαίνει ταχύτητα διάδοσης 1.5c 

Τι πάει τόσο στραβά;

Δεν το έκανα ακόμα αλλά θα πω μια ιδέα:

Αν έχεις το βιβλίο του Τσιμπούκη για το πεδίο (έχω την παλιά έκδοση, νμζ άλλαξε) έχει τον τύπο για
την φασική ταχύτητα. Γενικά σε μεταλλικούς κυματοδηγούς (όχι μόνο στις παράλληλες πλάκες) η φασική
ταχύτητα υp βγαίνει μεγαλύτερη από το c στον ελεύθερο χώρο (δες τον όρο στον παρονομαστή
για f>fc).

Επίσης, χωρίς να είμαι 100% σίγουρος, νομίζω πως αυτό δεν είναι παράλογο γιατί η ενέργεια διαδείδεται
με ug<c, όχι up που αφορά τα μέτωπα κύματος. Και πάλι δεν είμαι 100% βέβαιος για αυτό.

EDIT: και για του λόγου του αληθές, από την σχέση 15.75 (παλιά έκδοση πάντα) από το κεφάλαιο ΚΥΜΑΤΟΔΗΓΟΙ βγάζω για
f = 3.75Hz, fc = 3Hz up = 1.6666c, κοντά σε αυτό που λες αν σκεφτείς ότι περιμένουμε και
σφάλμα λόγω διακριτοποίησης.

[ChampionsNet]

EDIT: και για του λόγου του αληθές, από την σχέση 15.75 (παλιά έκδοση πάντα) από το κεφάλαιο ΚΥΜΑΤΟΔΗΓΟΙ βγάζω για
f = 3.75Hz, fc = 3Hz up = 1.6666c, κοντά σε αυτό που λες αν σκεφτείς ότι περιμένουμε και
σφάλμα λόγω διακριτοποίησης.

Ωωω τέλεια φίλε μου πάλι καλά γιατί θα έσκαγα. Είσαι ωραίος, ευχαριστώ!

[lhmma_Jordan]

Ωωω τέλεια φίλε μου πάλι καλά γιατί θα έσκαγα. Είσαι ωραίος, ευχαριστώ!

Ρίξε μια ματιά και εσύ να είσαι σίγουρος αλλά νομίζω το πέτυχα.

[lhmma_Jordan]

Σχετικά με το PML, πρέπει να απεικονίζεται και η διάδοση μέσα σε αυτό; Γενικά πως θα πρέπει να κάνουμε
την αξιολόγησή του; Αριθμητικά ή παρατηρώντας με το μάτι τι συμβαίνει για διάφορες τιμές;

[lhmma_Jordan]

Τι αλλαγές κάνατε στον κώδικα για το δεξί PML; Ενώ το αριστερό δουλεύει κανονικά, το δεξί μου δίνει όχι απλά ανάκλαση, τείνουν
στο άπειρο οι τιμές του Ez.  

Υ.Γ: Έκανα μια αλλαγή προσήμου και δούλεψε, αν και δεν ξέρω γιατί

[gkokk]

Και σε μένα είχε συμβεί το ίδιο αλλά με πρόλαβες στην απάντηση, το ίδιο είχα κάνει και δούλεψε  .Τώρα αν θες εξήγηση δεν ξέρω που άλλαξες εσύ πρόσημο, για μένα η αλλαγή προσήμου δεν ήταν τυχαία, προσωπικά εγώ το είχα λύσει με διόρθωση στην Hy_pml, η οποία κατα τα δεξιά έχει αντίθετο πρόσημο απτις τιμές κατά τα αριστερά και στην αρχή την είχα φτιάξει ωστέ να έχει ακριβώς τα ίδια πρόσημα

[lhmma_Jordan]

Και σε μένα είχε συμβεί το ίδιο αλλά με πρόλαβες στην απάντηση, το ίδιο είχα κάνει και δούλεψε  .Τώρα αν θες εξήγηση δεν ξέρω που άλλαξες εσύ πρόσημο, για μένα η αλλαγή προσήμου δεν ήταν τυχαία, προσωπικά εγώ το είχα λύσει με διόρθωση στην Hy_pml, η οποία κατα τα δεξιά έχει αντίθετο πρόσημο απτις τιμές κατά τα αριστερά και στην αρχή την είχα φτιάξει ωστέ να έχει ακριβώς τα ίδια πρόσημα

Thnx. Εγώ άλλαξα το πρόσημο του Hy που έχω ως είσοδο στο PML (στο Ez νμζ) , ακριβώς για τον λόγο που λες. Το θέμα είναι πως δεν έχω αναλυτική εξήγηση του γιατί θα αλλάξει έτσι η εξίσωση του PML απλά ότι δουλεύει γιατί στο PML θα έχεις το ίδιο ακριβώς πρόβλημα με αριστερά έτσι

[gkokk]

Ναι παρόλαυτα σκέψου και το αρχικό πρόβλημα χωρίς το PML, ή απλά οπτικοποίησε αντί της Ez την Hy. Βάσει των διαφορικών του Maxwell(και κατ'επέκταση της μεθόδου FDTD) φαίνεται ότι όταν διαδίδεται η Hy προς τα δεξιά, τα πρόσημα δεν ειναι "συμμετρικά" ως προς την Hy που διαδίδεται κατά τα αριστερά, όπως γίνεται π.χ. στην Ez και Hx. Άρα προφανώς όταν θα ορίσεις και τις PML συνθήκες πρέπει να έχεις αυτό στο μυαλό σου, και να μην προσπαθήσεις να τροποποιήσεις την προσήμανση "manually", απλά άσε τις εξισώσεις να κάνουν την δουλειά. Έγω είχα φανταστεί μια Hy_Pml που όντως λύνει το ίδιο πρόβλημα και απτις δυο πλευρές, μέχρι να συνειδητοποιήσω οτι στην διάδοση κατά τα δεξιά τα πρόσημα της Hy είναι αντίθετα. 

[lhmma_Jordan]

Ναι παρόλαυτα σκέψου και το αρχικό πρόβλημα χωρίς το PML, ή απλά οπτικοποίησε αντί της Ez την Hy. Βάσει των διαφορικών του Maxwell(και κατ'επέκταση της μεθόδου FDTD) φαίνεται ότι όταν διαδίδεται η Hy προς τα δεξιά, τα πρόσημα δεν ειναι "συμμετρικά" ως προς την Hy που διαδίδεται κατά τα αριστερά, όπως γίνεται π.χ. στην Ez και Hx. Άρα προφανώς όταν θα ορίσεις και τις PML συνθήκες πρέπει να έχεις αυτό στο μυαλό σου, και να μην προσπαθήσεις να τροποποιήσεις την προσήμανση "manually", απλά άσε τις εξισώσεις να κάνουν την δουλειά. Έγω είχα φανταστεί μια Hy_Pml που όντως λύνει το ίδιο πρόβλημα και απτις δυο πλευρές, μέχρι να συνειδητοποιήσω οτι στην διάδοση κατά τα δεξιά τα πρόσημα της Hy είναι αντίθετα. 

Ευχαριστώ, θα το κοιτάξω. Κάτι ακόμα, σου προέκυψε συμμετρικό αποτέλεσμα στις θύρες 3,4; Γιατί εμένα συμμετρικό μου βγαίνει και υποψιάζομαι ότι μπορεί
να μην είναι σωστό. Κατευθυντικό ζεύκτη μου θυμίζει αυτή η ιστορία με το άνοιγμα, που θα έπρεπε να δίνει μια θύρα απομονωμένη (μικρό ή καθόλου πεδίο) και μία διέλευσης,
όπου θα πήγαινε το μεγαλύτερο μέρος της ισχύος.

[gkokk]

δεν μέτρησα την θύρα 3, αλλα κοιτώντας όλο το φαινόμενο οπτικά ναι,αλλά γενικά θα έπρεπε θεωρώ, χωρίς να μπορώ να σκεφτώ κάποια μαθηματική απόδειξη βέβαια. Απλά και μόνο δομένου του ότι η οπή βρίσκεται στο κέντρο, είναι πάλι σαν να βλέπω φαινόμενο διάδοσης με διέγερση στο κέντρο, π.χ. όπως στη πρώτη περίπτωση. Για να επαναδιατυπώσω την ερώτηση σου, γιατί πέραν της υποψίας που έχεις να υπάρχει διαφορά στις θύρες?

edit: Να πω την αλήθεια όταν είδα το άνοιγμα και καμία άλλη αλλαγή εγώ σκέφτηκα ότι μοντελοποιούμε κάποιου είδους cross-talk ή περίπτωση που έχουμε κάποιο leakage στην μόνωση. Τώρα κουβέντα κάνουμε, αλλά οποιαδήποτε είδους μοντελοποίηση κατευθυντικού ζεύκτη θα θεωρούσα οτι θα έπρεπε να δίνεται ως στοιχείο στην εκφώνηση ώστε ακριβώς αν βλέπουμε ύποπτα αποτελέσματα να μπορούμε να έχουμε περεταίρω feedback.