Radiologická fyzika Ultrazvuk – vlnové vlastnosti podzim 2011, devátá přednáška Ultrazvuk Ultrazvuk je zvukové vlnění s frekvencí f > 20 kHZ. Pro ultrazvukovou diagnostiku v medicíně se používají frekvence řádu jednotek až desítek MHz. λ vlastnosti prostředí vlastnosti zdroje Základní rovnice Základními rovnicemi jsou rovnice kontinuity a Eulerova rovnice Pro malé kmity ponecháme v rovnicích jen členy prvního řádu v malých veličinách δρ, δp a v Stejně jako každý pohyb v ideální tekutině je i šíření zvuku děj adiabatický. Proto můžeme psát takže dostáváme rovnici kontinuity a Eulerovu rovnici ve tvaru Rychlost pohybu elementu prostředí Elementárními úpravami dvou základních rovnic hydrodynamiky 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 totiž vektor rychlosti daného elementu jako gradient nějaké potenciálové funkce φ a místní změnu tlaku jako záporně vzatou derivaci podle času této funkce násobenou hustotou, je tato funkce řešením vlnové rovnice (Δ je Laplaceův operátor) Potenciálová funkce Obecný tvar vlnové rovnice pro potenciálovou funkci je V kartézských souřadnicích je φ=φ(x,y,z) a gradient a laplacián mají jednoduché vyjádření Pokud jde o rozměr potenciálové funkce φ, máme nebo také Rovinná vlna Vlnová rovnice pro potenciálovou funkci, která závisí jen na čase a jedné kartézské souřadnici (třeba x) je Nejjednodušším příkladem je rovinná monochromatická vlna postupující v kladném směru osy x kde konstanta φm je amplituda a konstanta α značí fázi při x=c.t. Frekvence a vlnová délka spolu souvisí vztahem Přímým derivováním se přesvědčíme, že řešením jsou libovolné funkce tvaru Kulová vlna Vlnová rovnice pro potenciálovou funkci, která závisí jen na čase a radiální sférické souřadnici (značíme ji r) je Nejjednodušším příkladem je kulová monochromatická vlna šířící se z počátku r=0 kde konstanta φm je amplituda a konstanta α značí fázi při x=c.t. Frekvence a vlnová délka spolu souvisí vztahem Přímým derivováním se přesvědčíme, že řešením jsou libovolné funkce tvaru Kmitající koule jako zdroj 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)=vm.sin(ωt). Potom je výraz pro potenciálovou funkci v obecném bodě tvořen relativně jednoduchým vztahem Při ultrazvukové diagnostice máme vždy ωR>>c (jinak zapsáno R>>λ), takže se vztah pro potenciálovou funkci zjednoduší na Časová střední hodnota Z různých důvodů není zajímavá a mnohdy 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 Pro rychlost pohybu elementu prostředí máme a pro místní změnu tlaku Pro intenzitu zvukové vlny (tok energie jednotkovou ploškou za jednotku času) dostáváme ve střední hodnotě Odraz a lom na rozhraní Pro dopadající, odraženou a prošlou vlnu máme x y c1 c2 Podmínky spojitosti Na rozhraní (x=0) musí být u všech vln stejná závislost na souřadnici y (vlastnosti prostředí se v tomto směru nemění), což vede k tomu, že úhel odrazu je roven úhlu dopadu a poměr sinů úhlů dopadu a lomu (úhel se měří od kolmice k rozhraní) je roven poměru rychlostí Dále se musí rovnat tlak a kolmá složka rychlosti na obou stranách rozhraní odkud Intenzita odražené a lomené vlny Řešení rovnic popisujících podmínky na rozhraní dává Energie dopadající na jednotkovou plošku rozhraní resp. energie odražená jednotkovou ploškou nebo propuštěná jednotkovou ploškou je Koeficienty propustnosti a odrazivosti jsou definovány jako Akustická impedance je definována jako součin hustoty prostředí a rychlosti zvuku v tomto prostředí Pro kolmý dopad na rozhraní dvou prostředí Akustická impedance Jednotkou akustické impedance je rayleigh: 1 Rayl = 1 kg m–2 s–1 .Tabulka uvádí potřebné hodnoty pro některá prostředí. prostředí 103ρ [kg m–3] c [m s–1] Z [Rayl] vzduch 0,0012 330 0,0004 voda 1 1430 1,43 měkká tkáň 1,1 1540 1,69 játra 1,05 1570 1,65 tuk 0,95 1450 1,38 kost 1,91 4080 7,8 Dopplerův jev I Do detektoru přicházejí vlny v časech Frekvence je tedy Dopplerův jev II Do detektoru přicházejí vlny v časech Frekvence je tedy Dopplerův jev III Absorpce Ultrazvukové vlnění je při šíření prostředím částečně absorbováno. Jako u jiných druhů vlnění, můžeme pro pokles intenzity psát exponenciální zákon – pokles intenzity je úměrný dráze, kterou paprsek vlnění prochází U zvuku je obvyklé charakterizovat útlum nikoliv lineárním koeficientem útlumu s rozměrem [m–1], ale koeficientem s rozměrem [dB.m–1]. Logaritmováním zákona absorpce dostáváme takže Koeficient útlumu prostředí μ [dB cm–1] pro f=1 Mhz vzduch 1,61 voda 0,0025 měkká tkáň 0,5 – 1,0 játra 1,1 tuk 0,6 kost 10,0 Zápis vln pomocí komplexních čísel Tento zápis přináší velké zjednodušení, Platí totiž Můžeme si ověřit například triviální příklad Pro zjednodušení zápisu vln je důležité, ž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 bodu zapisujeme jako Skládání 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í 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 Skládání vln – rovinný problém Pole φ je součtem polí od 2.n+1 zdrojů vzájemně fázově posunutých a v poloze (x=n.d,z=0) Fáze je Pro x=X dostáváme maximum součtu, takže volbou fázového posunutí jednotlivých zdrojů posouváme toto maximum. Skládání vln – rovinný problém Celkem 21 zdrojů s d=2, λ=2 a L=100 pro X=0, 2 a 4 (všechno v mm).