Εξισώσεις της Μαγνητοϋδροδυναμικής

Εξισώσεις της Μαγνητοϋδροδυναμικής (MHD) ονομάζουμε τις μη γραμμικές διαφορικές εξισώσεις με μερικές παραγώγους, οι οποίες προκύπτουν από την σύζευξη των εξισώσεων του Maxwell για το Η/Μ πεδίο ,σε συνδυασμό με τις εξισώσεις Νάβιερ -Στόουκς , όταν εφαρμοστούν σε αγώγιμο ιξώδες ρευστό το οποίο κινείται μέσα σε ηλεκτρομαγνητικό πεδίο, ιδιαίτερα όταν τα ρεύματα που δημιουργούνται από την κίνηση του ρευστού μεταβάλλουν το μαγνητικό πεδίο, έτσι ώστε το ηλεκτρομαγνητικό πεδίο και οι εξισώσεις της δυναμικής να συζεύγνυνται. Οι εξισώσεις MHD εφαρμόζονται στα αγώγιμα ρευστά , είτε είναι υγρά είτε αέρια, στα οποία γίνονται συγκεκριμένες θεωρητικές παραδοχές (π.χ. το ρεύμα μετατόπισης του Maxwell θεωρείται αμελητέο) και το ρευστό μπορεί να μελετηθεί ως ένα συνεχές μέσο χωρίς φαινόμενα μέσης ελεύθερης διαδρομής. Οι εξισώσεις αυτές εφαρμόστηκαν αρχικά σε αστροφυσικά και γεωφυσικά προβλήματα όπου ακόμα εφαρμόζονται με επιτυχία στην περιγραφή της κίνησης ισχυρά ιονισμένων ρευστών ( π.χ. στην εξήγηση της κίνησης των δίσκων προσαύξησης γύρω από υπερμαζικά αστρικά σώματα όπως μαύρες τρύπες, αστέρες νετρονίων , στη περιγραφή του Πλάσματος στο ηλιακό Στέμμα (Αστρονομία) κ.λ.π.) Από τη δεκαετία του 50 με την κατασκευή των θερμοπυρηνικών αντιδραστήρων σύντηξης (ΤΟΚΑΜΑΚ ) για την παραγωγή ισχύος, οι εξισώσεις της Μαγνητοϋδροδυναμικής εφαρμόστηκαν για την μαθηματική μοντελοποίηση το υ περιορισμού του υπέρθερμου πλάσματος με ηλεκτρομαγνητικές δυνάμεις, έτσι ώστε να αποφεύγεται η επαφή του με τα τοιχώματα του αντιδραστήρα.

 Οι εξισώσεις της μαγνητοϋδροδυναμικής περιλαμβάνουν τις εξισώσεις του Maxwell ∇ × B → = μ 0 j {\displaystyle \nabla \times {\overrightarrow {B}}=\mu _{0}j} (νόμος Ampere Maxwell) (το ρεύμα μετατόπισης θεωρείται μηδενικό) ∇ × E → = − ϑ B → ϑ t {\displaystyle \nabla \times {\overrightarrow {E}}=-{\vartheta {\overrightarrow {B}} \over \vartheta t}} (νόμος της επαγωγής Faraday-Henry) ∇ B → = 0 {\displaystyle \nabla {\overrightarrow {B}}=0} (δεν υπάρχουν μονόπολα) και την πυκνότητα της ηλεκτρομαγνητικής δύναμης Lorentz την οποία την προσθέτουμε σαν δύναμη χώρου στην εξίσωση Navier-Stokes για ένα ασυμπίεστο Νευτώνιο ρευστό: d u d t = − ∇ ( p / ρ ) + g + ν ∇ 2 u {\displaystyle {du \over dt}=-\nabla (p/\rho )+g+\nu \nabla ^{2}u} (Navier-Stokes) (όπου ν= ο συντελεστής ιξώδους του ρευστού μετρημένος σε m2/sec) Δύναμη Lorentz Η δύναμη Lorentz J → × B → {\displaystyle {\overrightarrow {J}}\times {\overrightarrow {B}}} είναι δυνατόν εάν θέσουμε J → = ∇ × B → μ 0 {\displaystyle {\overrightarrow {J}}={\operatorname {\nabla \times {\overrightarrow {B}}} \over \operatorname {\mu _{0}} }} και χρησιμοποιώντας την αντίστοιχη μαθηματική ταυτότητα της διανυσματικής ανάλυσης, να γραφεί στη μορφή: J → × B → = ( B → ⋅ ∇ ) B μ 0 → − ∇ B 2 2 μ 0 {\displaystyle {\overrightarrow {J}}\times {\overrightarrow {B}}=({\overrightarrow {B}}\cdot \nabla ){\overrightarrow {\operatorname {B} \over \operatorname {\mu _{0}} }}-\nabla {\operatorname {B^{2}} \over \operatorname {2\mu _{0}} }} όπου ο δεύτερος όρος είναι η μαγνητική πίεση και ο πρώτος ονομάζεται τάση Faraday ή τάση Maxwell των δυναμικών γραμμών, εξαιτίας του γεγονότος ότι η δύναμη η οποία προκύπτει d F → = 1 μ 0 B → B → d S → {\displaystyle d{\overrightarrow {F}}={\operatorname {1} \over \operatorname {\mu _{0}} }{\overrightarrow {B}}{\overrightarrow {B}}d{\overrightarrow {S}}} είναι μηδενική στα τοιχώματα των ροϊκών σωλήνων. Στα 'άκρα' των δυναμικών γραμμών αυτές οι τάσεις παραμόρφωσης προκαλούν μια μη μηδενική δύναμη τάσης πάνω στους ροϊκούς σωλήνες. Σύμφωνα με τον νόμο του Gauss η δύναμη χώρου πάνω σε μια στοιχειώδη μάζα του ρευστού μπορεί να γίνει αντιληπτή ως το αποτέλεσμα της επιφανειακής ολοκλήρωσης αυτών των τάσεων και της πίεσης. Επιστρέφοντας στην έννοια μιάς δύναμης χώρου πάνω στις στοιχειώδεις μάζες του ρευστού μπορεί κανείς να αναλύσει το ανισότροπο κομμάτι του διανύσματος J → × B → {\displaystyle {\overrightarrow {J}}\times {\overrightarrow {B}}} σε εφαπτομενική και ορθοκανονική συνιστώσα της καμπύλης γραμμής η οποία είναι προσδεμένη πάνω στη δυναμική γραμμή: ( B → ⋅ ∇ ) B → = B ∂ B ∂ s e ^ t − B 2 R e ^ n {\displaystyle ({\overrightarrow {B}}\cdot \nabla ){\overrightarrow {B}}=B{\partial B \over \partial s}{\widehat {e}}_{t}-{\operatorname {B^{2}} \over \operatorname {R} }{\widehat {e}}_{n}} όπου R είναι η ακτίνα της τοπικής καμπυλότητας της δυναμικής γραμμής. Σημειώνουμε ότι αυτή η ανάλυση είναι ισοδύναμη με αυτό που μπορεί να γίνει με τον μη γραμμικό όρο ( u → ⋅ ∇ ) u → {\displaystyle ({\overrightarrow {u}}\cdot \nabla ){\overrightarrow {u}}} στην εξίσωση Navier - Stokes , εκτός από μια σημαντική διαφορά στο πρόσημο. Ο πρώτος όρος σε αυτή την περίπτωση είναι μια δύναμη η οποία αντιτίθεται στην επιτάχυνση η οποία σχετίζεται με την αδράνεια, δηλ. την επιτάχυνση μιάς στοιχειώδους μάζας του ρευστού. Ο δεύτερος όρος είναι τότε μια φυγόκεντρος δύναμη η οποία κατευθύνεται προς τα έξω. Οι όροι της εξίσωσης αντιπροσωπεύουν αντίστοιχα μια επιταχύνουσα δύναμη οποτεδήποτε το μαγνητικό πεδίο αυξάνει κατά μήκος μιάς δυναμικής γραμμής και μια δύναμη προς τα μέσα (κεντρομόλος) όταν οι δυναμικές γραμμές καμπυλώνονται. Αυτές οι δυνάμεις έχουν σημαντικές επιπτώσεις για την ευστάθεια ενός πλάσματος. Ενώ η αδρανειακή δύναμη αντιστέκεται στην επιτάχυνση , ο πρώτος όρος της εξίσωσης προκαλεί επιτάχυνση στην κατεύθυνση της αύξησης της έντασης του μαγνητικού πεδίου. Η δεύτερη συνιστώσα ενεργεί για να ευθυγραμμίσει τις μαγνητικές δυναμικές γραμμές, σύμφωνα με την ιδέα των δυνάμεων τάσης οι οποίες ενεργούν στα άκρα των δυναμικών γραμμών. Σε όρους αυτών των τάσεων η δύναμη στην κατεύθυνση της αύξησης του μαγνητικού πεδίου προκύπτει εξαιτίας της διαφοράς του μέτρου μεταξύ των αντίθετων τάσεων στα άκρα των δυναμικών γραμμών. Συμπερασματικά , η δύναμη Lorentz πάνω στο ρευστό μπορεί να θεωρηθεί ότι προκύπτει από την τάση που ασκείται όταν παραμορφώνονται οι δυναμικές γραμμές του μαγνητικού πεδίου μέσα στο ρευστό. Καθώς αυτές οι μαγνητικές γραμμές είναι σε κάποιο βαθμό 'παγωμένες' μέσα στο ρευστό αυτές οι τάσεις ως εκ τούτου ενεργούν ως δυνάμεις χώρου πάνω στο ρευστό. Σαν αποτέλεσμα οι δυναμικές γραμμές μπορούν να δρούν παρόμοια με μια χορδή κάτω από την επίδραση τάσης και προκαλούν μια κυματική κίνηση η οποία ονομάζεται 'κύματα Αλφεν' , όπου η απαραίτητη αδράνεια προσφέρεται από το ρευστό. Η χαρακτηριστική ταχύτητα με την οποία διαδίδονται τα μαγνητοακουστικά κύματα ονομάζεται ταχύτητα Alfven και ισούται με: V A = B 2 μ 0 ρ {\displaystyle V_{A}={\operatorname {B} \over {\sqrt {\operatorname {2\mu _{0}\rho } }}}} (Ταχύτητα Alfven )Στην MHD η μαγνητική επαγωγή Β και η ταχύτητα πλάσματος v προσδιορίζονται από την εξίσωση κίνησης: ρ d v d t = ρ ϑ v ϑ t + ρ ( v ∇ ) v = − ∇ p + j → × B → + ρ g + n ∇ 2 v {\displaystyle \rho {dv \over dt}=\rho {\vartheta v \over \vartheta t}+\rho (v\nabla \ )v=-\nabla \ p+{\overrightarrow {j}}\times {\overrightarrow {B}}+\rho g+n\nabla ^{2}v} όπου: − ∇ p = β α θ μ ι δ α π ι ϵ σ η ς {\textstyle -\nabla p=\beta \alpha \theta \mu \iota \delta \alpha \ \pi \iota \epsilon \sigma \eta \varsigma } (N/m3) J → × B → = H / M δ υ ν α μ η L o r e n t z ( N / m 3 ) {\displaystyle {\overrightarrow {J}}\times {\overrightarrow {B}}=\mathrm {H} /\mathrm {M} ~\delta \upsilon \nu \alpha \mu \eta ~Lorentz~(N/m^{3})} (πυκνότητα δύναμης) ρ g=δύναμη βαρύτητας (Ν/m3) n ∇ 2 v = δ υ ν α μ η α π τ o ι ξ ω δ ϵ ς {\displaystyle n\nabla ^{2}v=\delta \upsilon \nu \alpha \mu \eta \ \alpha \pi ~\tau o~\iota \xi \omega \delta \epsilon \varsigma } του ρευστού (N/m3) n= o συντελεστής ιξώδους του ρευστού με μονάδα μέτρησης ( kg /m sec) Επίσης άλλες εξισώσεις που περιλαμβάνονται στον φορμαλισμό της Μαγνητοϋδροδυναμικής είναι οι εξής: Η εξίσωση της συνέχειας για την μάζα και το φορτίο: ∂ ρ ∂ t + ∇ ( ρ V ) = 0 {\displaystyle {\partial \rho \over \partial t}+\nabla (\rho V)=0} εξίσωση συνέχειας για την μάζα όπου ρ = πυκνότητα του ρευστού ∂ ρ ∂ t + ∇ J → = 0 {\displaystyle {\partial \rho \over \partial t}+\nabla {\overrightarrow {J}}=0} εξίσωση διατήρησης για το φορτίο όπου ρ=πυκνότητα φορτίου και J=πυκνότητα ρεύματος Η Αδιαβατική εξίσωση ∂ p ∂ t + ( V → ⋅ ∇ ) p + γ p ∇ ⋅ V → = 0 {\displaystyle {\partial p \over \partial t}+({\overrightarrow {V}}\cdot \nabla )p+\gamma p\nabla \cdot {\overrightarrow {V}}=0} όπου γ=cp /cvH εξίσωση διατήρησης της ενέργειας : ∇ ⋅ ( E → × H → ) + ∂ ∂ t ( ρ m V 2 2 + p γ - 1 + B 2 2 μ 0 ) + n J 2 + ∇ ⋅ ( ρ m V 2 2 + p γ - 1 + B 2 2 μ 0 ) V → = 0 {\displaystyle \nabla \cdot ({\overrightarrow {E}}\times {\overrightarrow {H}})+{\partial \over \partial t}({\operatorname {\rho _{m}V^{2}} \over \operatorname {2} }+{\operatorname {p} \over \operatorname {\gamma -1} }+{B^{2} \over 2\mu _{0}})+nJ^{2}+\nabla \cdot ({\operatorname {\rho _{m}V^{2}} \over \operatorname {2} }+{\operatorname {p} \over \operatorname {\gamma -1} }+{B^{2} \over 2\mu _{0}}){\overrightarrow {V}}=0} Αδιάστατοι αριθμοί Προσθέτοντας μια δύναμη χώρου στην εξίσωση Navier - Stokes , προσθετει κανείς ένα επιπλέον αδιάστατο αριθμό στον ήδη γνωστό αριθμό Reynolds. η χαρακτηριστική αναλογία μεταξύ του όρου της δύναμης Lorentz και της αδράνειας δίνεται από μια 'παράμετρο αλληλεπίδρασης' η οποία ονομάζεται αριθμός Stuart και ορίζεται από την σχέση: N = | J → × B → | | ρ u ⋅ ∇ u → ⁡ | = σ B 2 l ρ u = H a 2 Re {\displaystyle N={\operatorname {|{\overrightarrow {J}}\times {\overrightarrow {B}}|} \over \operatorname {|\rho u\cdot \nabla {\overrightarrow {u}}} |}={\operatorname {\sigma B^{2}l} \over \operatorname {\rho u} }={\operatorname {Ha^{2}} \over \operatorname {Re} }} Εδώ ο αδιάστατος αριθμός Hartman δίνει εναλλακτικά την χαρακτηριστική αναλογία μεταξύ των ηλεκτρομαγνητικών δυνάμεων Lorentz και των δυνάμεων του ιξώδους. H a = | J → × B → | | ν ∇ 2 u | = B l σ ρ ν = N R e {\displaystyle Ha={\sqrt {{|{\overrightarrow {J}}\times {\overrightarrow {B}}|} \over |\nu \nabla ^{2}u|}}=Bl{\sqrt {{\sigma } \over \rho \nu }}={\sqrt {NRe}}} =αριθμός Hartman Στην εξίσωση της επαγωγής από την άλλη μεριά μπορεί να οριστεί ένας 'μαγνητικός αριθμός Reynolds' ως η χαρακτηριστική αναλογία μεταξύ του όρου της μεταφοράς και της διάχυσης. R e m = | ∇ × ( u → × B → ) | | λ ∇ 2 B → | {\displaystyle Re_{m}={\sqrt {{|\nabla \times ({\overrightarrow {u}}\times {\overrightarrow {B}})|} \over |\lambda \nabla ^{2}{{\overrightarrow {B}}|}}}} = ul λ = μ 0 σ u l {\displaystyle ={\operatorname {ul} \over \operatorname {\lambda } }=\mu _{0}\sigma ul} = μαγνητικός αριθμός Reynolds Χρησιμοποιώντας τις τιμές: ρ~4x10-7kg/m3 (n~1020m-3) , σ~2.4x109 S/m (~40 φορές περισσότερο από του χαλκού) , l~4m , B~4T , u~105m/s(M~0.15) και v~1m2/s οι οποίες είναι τυπικές τιμές για έναν ΤΟΚΑΜΑΚ τύπου ITER προκύπτει για τους αδιάστατους αριθμούς: Ν~4x1012 , Re~4x105 , Ha~109 και Rem ~109 (λ~3x10-4 m2 /sec ).Βλέπουμε έτσι ότι οι δυνάμεις Lorentz είναι πολύ μεγαλύτερες από τις αδρανειακές δυνάμεις(Ν>>1) και ότι οι αδρανειακές δυνάμεις είναι με την σειρά τους πολύ μεγαλύτερες από τις δυνάμεις ιξώδους (Re>>1) . Στην εξίσωση της επαγωγής τα μη ιδανικά φαινόμενα μπορούν εύκολα να αγνοηθούν(Rem >>1) . Θεώρημα Alfven Για ένα μαγνητικό πεδίο B → = B y ( x , t ) j ^ {\displaystyle {\overrightarrow {B}}=B_{y}(x,t){\widehat {j}}} και μια ροή ρευστού με ταχύτητα V → = v k ( x , t ) i ^ {\displaystyle {\overrightarrow {V}}=v_{k}(x,t){\widehat {i}}} η εξίσωση της επαγωγής ∂ B ∂ t = ∇ × ( V → × B → ) {\displaystyle {\partial B \over \partial t}=\nabla \times ({\overrightarrow {V}}\times {\overrightarrow {B}})} γίνεται ∂ B y ∂ t = − ∂ ∂ x ( v x B y ) ⇒ ∂ B y ∂ t + ∂ ∂ x ( v x B y ) = 0 {\displaystyle {\partial B_{y} \over \partial t}=-{\partial \over \partial x}(v_{x}B_{y})\Rightarrow {\partial B_{y} \over \partial t}+{\partial \over \partial x}(v_{x}B_{y})=0} (1) . Τώρα θεωρείστε την μαγνητική ροή μεταξύ δύο γραμμών οι οποίες είναι κάθετες στον άξονα x και περνούν από τα σημεία x=α και x=b, τα οποία κινούνται μαζί με το πλάσμα. Ετσι α=α(t) και b=b(t) έτσι ώστε d α dt = v x ( α , t ) , db dt = v x ( α , t ) ) {\displaystyle {\operatorname {d\alpha } \over \operatorname {dt} }=v_{x}(\alpha ,t),{\operatorname {db} \over \operatorname {dt} }=v_{x}(\alpha ,t))} H ροή είναι τότε F l u x ϕ = ∫ α ( t ) b ( t ) B y ( x , t ) d x {\displaystyle Flux~~\phi =\textstyle \int \limits _{\alpha (t)}^{b(t)}\displaystyle B_{y}(x,t)dx} . Σημειώστε ότι η ροή είναι μια συνάρτηση με μεταβλητή μόνο το χρόνο. Η χωρική εξάρτηση έχει ολοκληρωθεί. Η μαγνητική ροή μεταβάλλεται με το χρόνο χάρις σε δύο φαινόμενα. Δηλαδή επειδή μεταβάλλεται το Βy με το χρόνο και επειδή τα ακραία σημεία κινούνται με ταχύτητες υα και υb αυξάνοντας ή ελαττώνοντας τα όρια της ολοκλήρωσης. Διαφορίζοντας σε σχέση με το χρόνο πρέπει να θυμηθούμε να δίαφορίσουμε όχι μόνο το ολοκλήρωμα αλλά επίσης και τα όρια του ολοκληρώματος. Ετσι d dt ∫ α b B y d x = ∫ α b ∂ B y ∂ t d x + db dt B y ( b , t ) − d α dt B y ( α , t ) {\displaystyle {\operatorname {d} \over \operatorname {dt} }\textstyle \int \limits _{\alpha }^{b}\displaystyle B_{y}dx=\textstyle \int \limits _{\alpha }^{b}\displaystyle {\partial B_{y} \over \partial t}dx+{\operatorname {db} \over \operatorname {dt} }B_{y}(b,t)-{\operatorname {d\alpha } \over \operatorname {dt} }B_{y}(\alpha ,t)} Ως εκ τούτου μπορούμε να γράψουμε: d dt ∫ α b B y d x = ∫ α b ∂ B y ∂ t d x + [ v x B y ] α b = ∫ α b [ ∂ B y ∂ t + ∂ ∂ x ( v x B y ) ] d x = 0 ⇒ {\displaystyle {\operatorname {d} \over \operatorname {dt} }\textstyle \int \limits _{\alpha }^{b}\displaystyle B_{y}dx=\textstyle \int \limits _{\alpha }^{b}\displaystyle {\partial {B_{y}} \over \partial {t}}dx+[v_{x}B_{y}]_{\alpha }^{b}=\textstyle \int \limits _{\alpha }^{b}\displaystyle [{{\partial B_{y}} \over {\partial t}}+{{\partial } \over {\partial x}}(v_{x}B_{y})]dx=0\Rightarrow } d dt ∫ α b B y d x = 0 {\displaystyle {\operatorname {d} \over \operatorname {dt} }\textstyle \int \limits _{\alpha }^{b}\displaystyle B_{y}dx=0} Με άλλα λόγια , το ποσό της μαγνητικής ροής παραμένη σταθερό , και λέμε ότι η μαγνητική ροή είναι παγωμένη μέσα στο πλάσμα. Αυτό το αποτέλεσμα αποτελεί το θεώρημα του Alfven ή αλλιώς το θεώρημα της " παγωμένης μαγνητικής ροής " το οποίο απέδειξε ο Σουηδός Φυσικός Χάνες Αλφβέν το (1943), και το οποίο διατυπώνεται ως εξής: "Σε ένα τέλεια αγώγιμο ρευστό (Rm ⟶ {\displaystyle \longrightarrow } ∞ {\displaystyle \infty } ) , οι μαγνητικές δυναμικές γραμμές κινούνται μαζί με το ρευστό, ή αλλιώς : οι δυναμικές γραμμές είναι παγωμένες μέσα στο πλάσμα" Στην Ηλιακή Μαγνητοϋδροδυναμική απαλείφουμε το ηλεκτρικό πεδίο , Ε, και την πυκνότητα του ηλεκτρικού ρεύματος ,J, και εργαζόμαστε με την πρωταρχική μεταβλητή Β. Στην πραγματικότητα η σταθερά διάχυσης του μαγνητικού πεδίου δίνεται από τον Braginsky ως n=109T-3/2 m2s-1 ≈ {\displaystyle \approx } 1 στο ηλιακό στέμμα (Εάν θέσουμε Τ=106 Κ=θερμοκρασία στέμματος =>109(106)-3/2=10910-9=1). Ο τύπος L2/nτ είναι ένας αδιάστατος αριθμός , Rm ο οποίος λέγεται μαγνητικός αριθμός του Reynolds R m = L 2 n τ = L V n {\displaystyle R_{m}={L^{2} \over n\tau }={LV \over \ n}} όπου v είναι η τυπική ταχύτητα του πλάσματος. Στον Ηλιο ο μαγνητικός αριθμός Reynolds είναι φυσιολογικά πολύ μεγάλος επειδή τα L και V είναι τυπικά μεγάλα. Ο μαγνητικός αριθμός Reynolds είναι ένα μέτρο του μεγέθους του όρου της μαγνητικής μεταφοράς ∇ × ( v → × B → ) {\displaystyle \nabla \times ({\overrightarrow {v}}\times {\overrightarrow {B}})} συγκρινόμενος με το μέγεθος του όρου διάχυσης n ∇ 2 B → {\displaystyle n\nabla ^{2}{\overrightarrow {B}}} . Στις ηλιακές κηλίδες , οι οποίες είναι ψυχρότερες από το πλάσμα του ηλιακού στέμματος με ένα παράγοντα περίπου 100 , n=103m2s-1 , L=104km , V=1km s-1 =103ms-1 . Ως εκ τούτου ο μαγνητικός αριθμός Reynolds είναι ίσος με: R m = 10 7 ∙ 10 3 10 3 = 10 7 >> 1 {\displaystyle R_{m}={10^{7}\bullet 10^{3} \over 10^{3}}=10^{7}>>1} Εάν Rm >>1 τότε η εξίσωση του νόμου της επαγωγής προσεγγίζεται από την σχέση: ∂ B → ∂ t = ∇ × ( v → × B → ) {\displaystyle {\partial {\overrightarrow {B}} \over \partial t}=\nabla \times ({\overrightarrow {v}}\times {\overrightarrow {B}})} Εάν Rm <<1 τότε η εξίσωση του νόμου της επαγωγής προσεγγίζεται από τη σχέση: ∂ B → ∂ t = n ∇ 2 B → {\displaystyle {\partial {\overrightarrow {B}} \over \partial t}=n\nabla ^{2}{\overrightarrow {B}}} Μολονότι υπάρχουν πολύ παράμετροι που προσδιορίζουν έναν αντιδραστήρα θερμοπυρηνικής σύντηξης, υπάρχουν επίσης πάρα πολλές σχέσεις και περιορισμοί μεταξύ τους. Σε ένα θερμοπυρηνικό αντιδραστήρα σύντηξης η θερμοκρασία του πλάσματος φτάνει σε εκατομύρια βαθμούς Κelvin. Για να διατηρήσουμε ένα θερμό πλάσμα, πρέπει να το περιορίσουμε και να το κρατήσουμε μακριά από το τοίχωμα του θαλάμου κενού που το περιέχει. Η περισσότερα υποσχόμενη μέθοδος για ένα τέτοιο περιορισμό ενός υπέρθερμου πλάσματος είναι η χρήση κατάλληλων μαγνητικών πεδίων . Για τέτοια συστήματα μαγνητικού περιορισμού πρέπει να ισχύει μια συνθήκη ισορροπίας. Ισορροπία Πίεσης Όταν ένα πλάσμα βρίσκεται σε μόνιμη κατάσταση η μαγνητοϋδροδυναμική εξίσωση κίνησης δίνει: ∇ p = J → × B → {\displaystyle \nabla p={\overrightarrow {J}}\times {\overrightarrow {B}}} (1) και ∇ × B → = μ 0 J → {\displaystyle \nabla \times {\overrightarrow {B}}=\mu _{0}{\overrightarrow {J}}} (2) ∇ B → = 0 {\displaystyle \nabla {\overrightarrow {B}}=0} (3) ∇ J → = 0 {\displaystyle \nabla {\overrightarrow {J}}=0} (4). Από την εξίσωση ισορροπίας (1) προκύπτει B → ∇ p = B → ( J → × B → ) = 0 {\displaystyle {\overrightarrow {B}}\nabla p={\overrightarrow {B}}({\overrightarrow {J}}\times {\overrightarrow {B}})=0} (5) J → ∇ p = J → ( J → × B → ) = 0 {\displaystyle {\overrightarrow {J}}\nabla p={\overrightarrow {J}}({\overrightarrow {J}}\times {\overrightarrow {B}})=0} (6) H εξίσωση (5) δείχνει ότι τα διανύσματα B → {\displaystyle {\overrightarrow {B}}} και ∇ p {\displaystyle \nabla p} είναι ορθογώνια και οι επιφάνειες σταθερής πίεσης συμπίπτουν με τις μαγνητικές επιφάνειες. Η εξίσωση (6) δείχνει ότι το διάνυσμα πυκνότητας ρεύματος J → {\displaystyle {\overrightarrow {J}}} είναι παντού παράλληλο προς τις επιφάνειες σταθερής πίεσης. Με αντικατάσταση της (2) στην (1) προκύπτει: ∇ ( p + B 2 2 μ 0 ) = ( B → ⋅ ∇ ) B → μ 0 = B 2 μ 0 ( − 1 R n ^ + ∂ B / ∂ l B b ^ ) {\displaystyle \nabla (p+{\operatorname {B^{2}} \over \operatorname {2\mu _{0}} })=({\overrightarrow {B}}\cdot \nabla ){\operatorname {\overrightarrow {B}} \over \operatorname {\mu _{0}} }={\operatorname {B^{2}} \over \operatorname {\mu _{0}} }(-{\operatorname {1} \over \operatorname {R} }{\widehat {n}}+{\operatorname {\partial {B}/\partial {l}} \over \operatorname {B} }{\widehat {b}})} (7) Οπου χρησιμοποιήσαμε τις ακόλουθες σχέσεις: B → × ( ∇ × B → ) + ( B → ⋅ ∇ ) B → = ∇ ( B 2 / 2 ) , ( B → ⋅ ∇ ) B → = B 2 [ ( b ^ ∇ ) b ^ + b ^ ( ( b ^ ∇ B → ) B ) ] {\displaystyle {\overrightarrow {B}}\times (\nabla \times {\overrightarrow {B}})+({\overrightarrow {B}}\cdot \nabla ){\overrightarrow {B}}=\nabla (B^{2}/2),~~({\overrightarrow {B}}\cdot \nabla ){\overrightarrow {B}}=B^{2}[({\widehat {b}}\nabla ){\widehat {b}}+{\widehat {b}}({\operatorname {({\widehat {b}}\nabla {\overrightarrow {B}})} \over \operatorname {B} })]} Οπου R είναι η ακτίνα καμπυλότητας της μαγνητικής δυναμικής γραμμής και n ^ {\displaystyle {\widehat {n}}} είναι το μοναδιαίο διάνυσμα το οποίο κατευθύνεται προς ένα σημείο πάνω στη δυναμική γραμμή από το κέντρο της καμπύλης. (l) είναι το μήκος της καμπύλης στη διεύθυνση της δυναμικής γραμμής. Βρίσκουμε ότι το δεξί σκέλος της εξ. (7) μπορεί να αγνοηθεί όταν η ακτίνα καμπυλότητας είναι πολύ μεγαλύτερη από το μήκος στο οποίο το μέτρο της p μεταβάλλεται σημαντικά , δηλαδή το μέγεθος του πλάσματος , και η μεταβολή του Β κατά μήκος της δυναμικής γραμμής είναι πολύ μικρότερη από τη μεταβολή του B → {\displaystyle {\overrightarrow {B}}} στην κάθετη διεύθυνση. Τότε η (7) γίνεται: p + B 2 2 μ 0 ∽ B 0 2 2 μ 0 {\displaystyle p+{\operatorname {B^{2}} \over \operatorname {2\mu _{0}} }\backsim {\operatorname {B_{0}^{2}} \over \operatorname {2\mu _{0}} }} όπου Β0 είναι η τιμή του μαγνητικού πεδίου στο όριο του πλάσματος(p=0). Όταν το σύστημα είναι αξονικά συμμετρικό και ∂ ∂ z = 0 {\displaystyle {\partial \over \partial z}=0} η (7) συμπτύσσεται στην ∂ ∂ r ( P + B z 2 + B θ 2 2 μ 0 ) = − B θ 2 r μ 0 {\displaystyle {\partial \over \partial r}(P+{\operatorname {B_{z}^{2}+B_{\theta }^{2}} \over \operatorname {2\mu _{0}} })=-{\operatorname {B_{\theta }^{2}} \over \operatorname {r\mu _{0}} }} (8) Πολλαπλασιάζοντας την (8) επί r2 και ολοκληρώνοντας κατά μέλη παίρνουμε: ( P + B z 2 + B θ 2 2 μ 0 ) v = α = 1 π α 2 ∫ 0 α ( P + B z 2 2 μ 0 ) 2 π r d r {\displaystyle (P+{\operatorname {B_{z}^{2}+B_{\theta }^{2}} \over \operatorname {2\mu _{0}} })_{v=\alpha }={\operatorname {1} \over \operatorname {\pi \alpha ^{2}} }\int _{0}^{\alpha }(P+{\operatorname {B_{z}^{2}} \over \operatorname {2\mu _{0}} })2\pi rdr} δηλαδή < P > + < B z 2 > 2 μ 0 = P α + B z 2 ( α ) + B θ 2 2 μ 0 {\displaystyle <P>+{\operatorname {<B_{z}^{2}>} \over \operatorname {2\mu _{0}} }=P_{\alpha }+{\operatorname {B_{z}^{2}(\alpha )+B_{\theta }^{2}} \over \operatorname {2\mu _{0}} }} (9) όπου το σύμβολο < > σημαίνει μέση τιμή. Καθώς το B 2 2 μ 0 {\displaystyle {\operatorname {B^{2}} \over \operatorname {2\mu _{0}} }} είναι η πίεση του μαγνητικού πεδίου Β0 η εξ. (9) είναι η εξίσωση ισορροπίας της πίεσης. Ο λόγος της πίεσης του πλάσματος προς την πίεση του εξωτερικού μαγνητικού πεδίου Β0 β = P B 0 2 / 2 μ 0 = n ( T e + T i ) B 0 2 ⁡ / 2 μ 0 {\displaystyle \beta ={\operatorname {P} \over \operatorname {B_{0}^{2}/2\mu _{0}} }={\operatorname {n(T_{e}+T_{i})} \over \operatorname {\operatorname {B_{0}^{2}} /\operatorname {2\mu _{0}} } }} ονομάζεται λόγος βήτα . Για ένα περιορισμένο πλάσμα , το β είναι πάντοτε μικρότερο από τη μονάδα και χρησιμοποιείται σαν ένα μέτρο της αξίας του μαγνητικού πεδίου περιορισμού πλάσματος. Το γεγονός ότι το εσωτερικό μαγνητικό πεδίο είναι μικρότερο από το εξωτερικό δείχνει το διαμαγνητισμό του πλάσματος. Μαγνητοϋδροδυναμική Αστάθεια Πλάσματος Για μιά πίεση 1Atm απαιτείται ένα μαγνητικό πεδίο 5044 Gauss για τον μαγνητικό περιορισμό του πλάσματος. Ένα σωληνοειδές θα έπρεπε να έχει 4000 Αμπεροστροφές/εκατοστό για να παραγάγει ένα τέτοιο αξονικό πεδίο, πράγμα όχι πολύ εύκολο. Το σύνηθες φαινόμενο "pinch" (=στένωση των δυναμικών γραμμών για τον περιορισμό του πλάσματος), δημιουργείται από την ίδια τη ροή του πλάσματος μεσα στο θάλαμο του τοροειδούς. Το μαγνητικό πεδίο έχει την μορφή δακτύλιων γύρω από το ρεύμα. Εάν το πλάσμα είναι περιορισμένο σε μια ακτίνα R τότε η μέση πίεση είναι P ∼ ( I R ) 2 [ P = B 2 2 μ 0 , B = μ 0 I 2 π r ] {\displaystyle P\sim ({\operatorname {I} \over \operatorname {R} })^{2}[P={\operatorname {B^{2}} \over \operatorname {2\mu _{0}} },~B={\operatorname {\mu _{0}I} \over \operatorname {2\pi r} }]} Εάν θεωρήσουμε το ρεύμα να περιορίζεται σε ένα μικρότερο εμβαδό διατομής , μια "στένωση" σε ένα σημείο, τότε μικρότερη ακτίνα σημαίνει μεγαλύτερο μαγνητικό πεδίο . Αυτό με τη σειρά του σημαίνει μεγαλύτερη μαγνητική πίεση, η οποία δημιουργεί ακόμα μεγαλύτερη στένωση. Εάν το ρεύμα αποκλίνει από την ευθεία σε μια "τεθλασμένη" , τότε οι γραμμές της επαγωγής πλησιάζουν μεταξύ τους ( έτσι το Β αυξάνει) στο εσωτερικό τμήμα της τεθλασμένης, ενώ το αντίθετο συμβαίνει στην εξωτερική πλευρά της τεθλασμένης. Αυτό δημιουργεί παραπέρα τεθλασμένες τροχιές , και έτσι δημιουργείται μια αστάθεια στον περιορισμό του πλάσματος. Davidson, P.A. (2001). An Introduction to Magnetohydrodynamics. The Pitt Building , Trumpington Street,Cambridge UK: Cambridge University Press. ISBN 0 521-79149-9. Jackson, John David (1962). Classical Electrodynamics. John Wiley & Sons. ISBN 0 471 43131 1. Kenro Miyamoto (2011).'Fundamentals of Plasma Physics and Controlled Fusion'. Iwanami Book Service Center. (ISBN 4-900491-11-X)