Model v elektromagnetismu
The Havriliak – Negami relaxace je empirická modifikace Debye relaxace model v elektromagnetismu. Na rozdíl od modelu Debye je relaxace Havriliak – Negami odpovědná za asymetrie a šíři dielektrická disperze křivka. Model byl poprvé použit k popisu dielektrické relaxace některých polymery ,[1] přidáním dvou exponenciální parametry k Debye rovnici:
ε ^ ( ω ) = ε ∞ + Δ ε ( 1 + ( i ω τ ) α ) β , { displaystyle { hat { varepsilon}} ( omega) = varepsilon _ { infty} + { frac { Delta varepsilon} {(1+ (i omega tau) ^ { alpha}) ^ { beta}}},} kde ε ∞ { displaystyle varepsilon _ { infty}} je permitivita na vysokofrekvenčním limitu, Δ ε = ε s − ε ∞ { displaystyle Delta varepsilon = varepsilon _ {s} - varepsilon _ { infty}} kde ε s { displaystyle varepsilon _ {s}} je statická nízkofrekvenční permitivita a τ { displaystyle tau} je charakteristika čas na odpočinek média. Exponenty α { displaystyle alpha} a β { displaystyle beta} popsat asymetrii a šíři odpovídajících spekter.
V závislosti na aplikaci je Fourierova transformace natažená exponenciální funkce může být životaschopnou alternativou, která má o jeden parametr méně.
Pro β = 1 { displaystyle beta = 1} rovnice Havriliak – Negami se redukuje na Cole – Coleova rovnice , pro α = 1 { displaystyle alpha = 1} do Cole – Davidsonova rovnice .
Matematické vlastnosti Skutečné a imaginární části Úložná část ε ′ { displaystyle varepsilon '} a ztrátová část ε ″ { displaystyle varepsilon ''} permitivity (zde: ε ^ ( ω ) = ε ′ ( ω ) − i ε ″ ( ω ) { displaystyle { hat { varepsilon}} ( omega) = varepsilon '( omega) -i varepsilon' '( omega)} ) lze vypočítat jako
ε ′ ( ω ) = ε ∞ + Δ ε ( 1 + 2 ( ω τ ) α cos ( π α / 2 ) + ( ω τ ) 2 α ) − β / 2 cos ( β ϕ ) { displaystyle varepsilon '( omega) = varepsilon _ { infty} + Delta varepsilon left (1 + 2 ( omega tau) ^ { alpha} cos ( pi alpha / 2) + ( omega tau) ^ {2 alpha} vpravo) ^ {- beta / 2} cos ( beta phi)} a
ε ″ ( ω ) = Δ ε ( 1 + 2 ( ω τ ) α cos ( π α / 2 ) + ( ω τ ) 2 α ) − β / 2 hřích ( β ϕ ) { displaystyle varepsilon '' ( omega) = Delta varepsilon left (1 + 2 ( omega tau) ^ { alpha} cos ( pi alpha / 2) + ( omega tau) ^ {2 alpha} right) ^ {- beta / 2} sin ( beta phi)} s
ϕ = arktan ( ( ω τ ) α hřích ( π α / 2 ) 1 + ( ω τ ) α cos ( π α / 2 ) ) { displaystyle phi = arctan left ({( omega tau) ^ { alpha} sin ( pi alpha / 2) nad 1 + ( omega tau) ^ { alpha} cos ( pi alpha / 2)} vpravo)} Ztráta ztráty Maximum ztrátové části leží na
ω m A X = ( hřích ( π α 2 ( β + 1 ) ) hřích ( π α β 2 ( β + 1 ) ) ) 1 / α τ − 1 { displaystyle omega _ { rm {max}} = left ({ sin left ({ pi alpha over 2 ( beta +1)} right) over sin left ({ pi alpha beta over 2 ( beta +1)} right)} right) ^ {1 / alpha} tau ^ {- 1}} Superpozice Lorentzianů Relaxaci Havriliak – Negami lze vyjádřit jako superpozici jednotlivých relaxací Debye
ε ^ ( ω ) − ϵ ∞ Δ ε = ∫ τ D = 0 ∞ 1 1 + i ω τ D G ( ln τ D ) d ln τ D { displaystyle {{ hat { varepsilon}} ( omega) - epsilon _ { infty} přes Delta varepsilon} = int _ { tau _ {D} = 0} ^ { infty} {1 nad 1 + i omega tau _ {D}} g ( ln tau _ {D}) d ln tau _ {D}} s distribuční funkcí
G ( ln τ D ) = 1 π ( τ D / τ ) α β hřích ( β θ ) ( ( τ D / τ ) 2 α + 2 ( τ D / τ ) α cos ( π α ) + 1 ) β / 2 { displaystyle g ( ln tau _ {D}) = {1 over pi} {( tau _ {D} / tau) ^ { alpha beta} sin ( beta theta) nad (( tau _ {D} / tau) ^ {2 alpha} +2 ( tau _ {D} / tau) ^ { alpha} cos ( pi alfa) +1) ^ { beta / 2}}} kde
θ = arktan ( hřích ( π α ) ( τ D / τ ) α + cos ( π α ) ) { displaystyle theta = arctan left ({ sin ( pi alpha) over ( tau _ {D} / tau) ^ { alpha} + cos ( pi alpha)} vpravo )} pokud je argument arkustangensu pozitivní, jinak[2]
θ = arktan ( hřích ( π α ) ( τ D / τ ) α + cos ( π α ) ) + π { displaystyle theta = arctan left ({ sin ( pi alpha) over ( tau _ {D} / tau) ^ { alpha} + cos ( pi alpha)} vpravo ) + pi} Logaritmické momenty První logaritmický moment tohoto rozdělení, průměrná logaritmická relaxační doba je
⟨ ln τ D ⟩ = ln τ + Ψ ( β ) + E u α { displaystyle langle ln tau _ {D} rangle = ln tau + { Psi ( beta) + { rm {Eu}} přes alpha}} kde Ψ { displaystyle Psi} je funkce digamma a E u { displaystyle { rm {Eu}}} the Eulerova konstanta .[3]
Inverzní Fourierova transformace Inverzní Fourierovu transformaci funkce Havriliak-Negami (odpovídající relaxační funkce v časové oblasti) lze vypočítat numericky.[4] Lze ukázat, že zapojené série expanzí jsou speciální případy Funkce Fox – Wright .[5] Zejména v časové doméně odpovídající ε ^ ( ω ) { displaystyle { hat { varepsilon}} ( omega)} lze reprezentovat jako
X ( t ) = ε ∞ δ ( t ) + Δ ε τ ( t τ ) α β − 1 E α , α β β ( − ( t / τ ) α ) , { displaystyle X (t) = varepsilon _ { infty} delta (t) + { frac { Delta varepsilon} { tau}} vlevo ({ frac {t} { tau}} vpravo) ^ { alpha beta -1} E _ { alpha, alpha beta} ^ { beta} (- (t / tau) ^ { alpha}),} kde δ ( t ) { displaystyle delta (t)} je delta funkce Dirac a
E α , β y ( z ) = 1 Γ ( y ) ∑ k = 0 ∞ Γ ( y + k ) z k k ! Γ ( α k + β ) { displaystyle E _ { alpha, beta} ^ { gamma} (z) = { frac {1} { Gamma ( gamma)}} součet _ {k = 0} ^ { infty} { frac { Gamma ( gamma + k) z ^ {k}} {k! Gamma ( alpha k + beta)}}} je speciální instance Funkce Fox – Wright a přesně jsou to tři parametry Funkce Mittag-Leffler [6] také známá jako funkce Prabhakar. Funkce E α , β y ( z ) { displaystyle E _ { alpha, beta} ^ { gamma} (z)} lze numericky vyhodnotit například pomocí kódu Matlab.[7]
Reference ^ Havriliak, S .; Negami, S. (1967). "Složitá rovinová reprezentace dielektrických a mechanických relaxačních procesů v některých polymerech". Polymer . 8 : 161–210. doi :10.1016/0032-3861(67)90021-3 . ^ Zorn, R. (1999). „Použitelnost distribučních funkcí pro spektrální funkci Havriliak – Negami“. Journal of Polymer Science Part B . 37 (10): 1043–1044. Bibcode :1999JPoSB..37.1043Z . doi :10.1002 / (SICI) 1099-0488 (19990515) 37:10 <1043 :: AID-POLB9> 3.3.CO; 2-8 . ^ Zorn, R. (2002). „Logaritmické momenty rozdělení časů relaxace“ (PDF) . Journal of Chemical Physics . 116 (8): 3204–3209. Bibcode :2002JChPh.116.3204Z . doi :10.1063/1.1446035 . ^ Schönhals, A. (1991). "Rychlý výpočet časově závislé dielektrické permitivity pro funkci Havriliak-Negami". Acta Polymerica . 42 : 149–151. ^ Hilfer, J. (2002). "H -funkční reprezentace pro napjatou exponenciální relaxaci a non-Debyeovu náchylnost ve skleněných systémech ". Fyzický přehled E . 65 : 061510. Bibcode :2002PhRvE..65f1510H . doi :10.1103 / physreve.65.061510 . ^ Gorenflo, Rudolf; Kilbas, Anatoly A .; Mainardi, Francesco; Rogosin, Sergei V. (2014). Springer (ed.). Funkce Mittag-Leffler, související témata a aplikace . ISBN 978-3-662-43929-6 . ^ Garrappa, Roberto. "Funkce Mittag-Leffler" . Citováno 3. listopadu 2014 . Viz také