Katedra elektromagnetického pole FEL ČVUT v Praze
Buňka – sférický rezonátor COMSOL Multiphysics
Biologická buňka jako sférický dutinový elektromagnetický resonátor Vzorová úloha v COMSOL Multiphysics Tento dokument byl vytvořen v rámci projektu FRVŠ 2300/2008
•
Úvod
Biologické buňky jsou různých tvarů a jejich membrány je oddělují od vnějšího prostředí. Rozdílný index lomu (odmocnina z permitivity) cytoplasmy buněk, buněčné membrány a vnějšího prostředí je fyzikálním předpokladem vzniku rezonancí – elektromagnetických vidů – uvnitř buňky po vybuzení vhodnou frekvencí elektromagnetického pole. Rozměry buněk jsou řádově v mikrometrech. Pro jednoduchost uvažujme model buňky reprezentovaný sférickou dutinou s dokonale vodivými stěnami vyplněný materiálem o relativní permitivitě 1.5. Ukážeme jak spočítat v COMSOL Multiphysics vidy dutinového elektromagnetického rezonátoru. • Nákres modelu Spusťte COMSOL Multiphysics, v Model navigator zvolte: 9 9
Space dimension : 3D RF Module, Electromagnetic Waves, Eigenfrequency analysis
Nákres sférické dutiny - buňky: 9 V Draw menu zvolte Sphere a nastavte následující parametry:
v. 1.0 Michal Cifra
1/5
Katedra elektromagnetického pole FEL ČVUT v Praze
•
Buňka – sférický rezonátor COMSOL Multiphysics
Nastavení subdomén
Tady nadefinujeme elektromagnetické charakteristiky modelovaného materiálu: 9 Ve Physics menu zvolte Subdomain Settings. Jediná doména je č. 1 je zvolena 9 Vyplňte kolonku εr hodnotou 1.5
•
Nastavení hranic
Tady nastavíme okrajové podmínky modelu. 9 Ve Physics menu zvolte Boundary Settings a ve výběru hranic zvolte zároveň všechny hranice 1-8 9 V Boundary condition nastavte PEC – dokonalý elektrický vodič (perfect electric conductor)
v. 1.0 Michal Cifra
2/5
Katedra elektromagnetického pole FEL ČVUT v Praze
•
Buňka – sférický rezonátor COMSOL Multiphysics
Nastavení mřížky
Pro spočítání veličin elektromagnetického pole je potřeba diskretizovat model pomocí vhodné mřížky. Základní mřížka se dá snadno nadefinovat následovně: 9 V Mesh menu zvolte Initialize Mesh •
Řešení problému
Pro vyřešení problému je důležité odhadnout, kde ve frekvenční oblasti má řešič rezonanční frekvence hledat. Odhadovanou frekvenci je potřeba nastavit nad mezní frekvenci dominantního vidu jinak může řešení divergovat. Čím více vidů chceme najít tím výše je potřeba odhadovanou frekvenci nastavit. 9 V Solve menu zvolte Solver parameters, tam pro nalezení prvních 20 vidů nastavte podle obrázku:
v. 1.0 Michal Cifra
3/5
Katedra elektromagnetického pole FEL ČVUT v Praze
9 9 9
Buňka – sférický rezonátor COMSOL Multiphysics
V Solve menu zvolte Solve Problem Výpočetní čas je 25.6 s s použitím Intel Core2 Duo E8400 3 GHz Počet stupňů volnosti (degrees of freedom) je 60824
•
Zpracování výsledků
•
Zobrazení vidů a jejich mezních frekvencí 9 V Postprocessing menu zvolte Plot Parameters 9 V záložce General sekce Solution to use si můžete prohlédnout nalezené frekvence 9 Pro zobrazení dominantního vidu zvolte první frekvenci (u sférického rezonátoru z důvodu symetrie 3x degenerovaný vid). POZOR: První frekvence v seznamu nemusí být dominantní vid, jestli jste nastavili frekvenci kolem které má řešič hledat moc vysoko nebo jste nastavili příliš malý počet hledaných vidů. V této úloze bylo zvoleno optimální nastavení, aby k nalezení dominantního vidu došlo.
9 V Plot Parameters zvolte kartu Slice a nastavte sekci Slice data a sekci slice
positioning podle obrázku:
v. 1.0 Michal Cifra
4/5
Katedra elektromagnetického pole FEL ČVUT v Praze
Buňka – sférický rezonátor COMSOL Multiphysics
9 Stisknete OK pro vykreslení následujícího obrázku – rozložení intenzity elektrického pole dominantního vidu (jedné z třech degenerovaných složek) sférického rezonátoru
v. 1.0 Michal Cifra
5/5
