Radiologická fyzika Ultrazvuk – vlnové vlastnosti podzim 2010, osmá přednáška Poznámka Tato presentace se od předchozích odlišuje tím, že obsahuje matematicky náročnější části – pracuje se s funkcemi komplexní proměnné a integrály. Nebude proto asi pro většinu posluchačů zajímavá. Ultrazvuk Ultrazvuk je zvukové vlnění s vysokou frekvencí. Pro ultrazvukovou diagnostiku v medicíně se používají frekvence řádu jednotek až desítek MHz. λ vlastnosti prostředí vlastnosti zdroje Rychlost pohybu elementu prostředí Elementárními úpravami základnách rovnic hydrodynamiky (rovnice kontinuity a Eulerova rovnice) dospějeme k výrazu pro rychlost zvuku (čtverec rychlosti zvuku je dán změnou tlaku v daném elementu prostředí se změnou hustoty v tomto elementu, aniž si sousední elementy vyměňují teplo). Zapíšeme-li ještě vektor rychlosti daného elementu zapíšeme jako gradient nějaké potenciálové funkce φ, je tato funkce řešením vlnové rovnice (Δ je Laplaceův operátor) V kartézských souřadnicích φ=φ(x,y,z) a Kmitající koule Mějme kouli poloměru R vykonávající malé pulsace, tj. body na povrchu koule radiálně kmitají rychlostí u(t) R Obecná kulová vlna Mějme kouli poloměru R vykonávající malé pulsace, tj. body na povrchu koule radiálně kmitají rychlostí u(t). Řešení úlohy (tedy řešení vlnové rovnice) hledáme ve tvaru Na povrchu musí být rychlosti pohybu elementu budící koule i elementu okolního prostředí stejné takže řešením je Harmonická kulová vlna Předpokládejme teď, že u(t)=vRsin(ωt). Potom je výraz pro rychlost pohybu elementu v obecném bodě tvořen relativně jednoduchým vztahem Při ultrazvukové diagnostice máme vždy ωr>ωR>>c (jinak zapsáno r>R>>λ), takže se vztah pro rychlost pohybu elementu prostředí zjednoduší na Časová střední hodnota Z různých důvodů není zajímavá a mnohdu ani dobře měřitelná okamžitá hodnota fyzikální veličiny F(t), ale její střední hodnota za dobu T , tj. Tato doba bývá někdy velmi dlouhá (čekáme, až se nějaký jev ustálí a nezajímá nás přechodový jev). Pro veličinu s periodou ω=2πf je touto dobou přirozeně perioda T=1/f Zejména Intenzita harmonické kulové vlny Odvodili jsme výraz pro rychlost pohybu elementu prostředí Intenzita zvukové vlny je neboli Komplexní čísla Základem je třeba zavedení pojmu komplexní jednotky, čísla, které neleží na reálné ose Pomocí komplexní jednotky zavádíme komplexní čísla. Jsou-li a ,b, c, d reálná čísla, pro komplexní čísla w=a+i.b , z=c+i.d Každý si jistě pamatuje, že platí Pomocí komplexní čísel můžeme psát také Komplexní jednotka Jsou-li a ,b reálná čísla, definujeme ke komplexní číslu w číslo komplexně sdružené w* a modul |w| 1 i -1 -i cosα isinα α eiα Eulerův vztah Platí Zápis vln pomocí komplexních čísel Tento zápis přináší velké zjednodušení. Je například velmi výhodné, že nedochází při derivování (a integrování) k „přecházení“ od sinu ke kosinu Monochromatickou kulovou vlnu s frekvencí ω, která vychází z nějakého bodu zdroje zapisujeme jako Zdroj ultrazvukových vln není ale nikdy tak jednoduchý. Má jisté nezanedbatelné geometrické rozměry a budí vlny s daleko složitější časovou závislostí než jsou prosté harmonické kmity. Při odvozování obecně platných vztahů (a v mnoha dalších aplikacích) hraje velmi důležitou roli Fourierova transformace. Také její zápis pomocí komplexních proměnných je podstatně jednodušší. Skládání dvou harmonických vln – geometrie Dva zdroje v různých bodech kmitající ve fázi Vhodnou volbou souřadnicového systému (body A a B leží na ose x, počátek půlí jejich vzdálenost) můžeme psát A B O P Skládání dvou harmonických vln – aproximace Pole φ v bodě P je součtem polí φA a φB Závislost na d v čitateli zanedbáme, jde o pomalou změnu amplitudy, takže Totéž zanedbání nemůžeme ale provést v exponentu, protože jako součin velkého (k=2π/λ) a malého čísla nemusí být malé (ve srovnání s π, kdy při změně o tuto hodnotu se obrátí znaménko celého výrazu pro vlnu). Máme tedy Fourierova transformace Máme-li nějakou funkci času f(t), pak (za jistých matematických předpokladů) je možno ji zapsat jakou lineární kombinaci harmonických funkcí Funkce F(ω) je pak Říkáme, že funkce f(t) a F(ω) jsou navzájem zobrazeny pomocí Fourierovy transformace. Fourierovy transformace využijeme dále při hledání řešení vlnové rovnice pro ultrazvukové pole (například popsané lokální změnou tlaku) vyvolané nějakým kmitajícím zdrojem. Vlnová rovnice Vlnovou rovnici pro časově závislé tlakové pole převedeme Fourierovou transformací na rovnici pro harmonickou vlnu Známe-li rozložení tlakového pole na ploše zdroje (a při splnění některých dalších vcelku přirozených podmínek), získáme řešení rovnice pomocí Huygensova – Fresnelova principu. Řešení pro monochromatickou vlnu Huygensův – Fresnelův princip Integruje se po ploše detektoru S. (x,y,z) (xS,yS,zS) (0,0,0) Obecné řešení Pomocí Fourierovy transformace přejdeme od řešení pro monochromatickou vlnu k obecnému časově závislému řešení vlnové rovnice pro tlakové pole Volbou rozložení tlakového pole na ploše zdroje můžeme dosáhnout nejrozmanitějších profilů ultrazvukového svazku jak v prostoru tak v čase. Dříve se pracovalo pouze s nejjednodušším tvarem a v analogii se světelnou optikou byl ve svazku rozlišován Fresnelův nebo Fraunhoferův profil.