next up previous
Next: Gaussova kvadratura Up: Metoda robnih elementov - Previous: Robni elementi

Robni pogoji

V nasi nalogi imamo izracun toka idelanega fluida. Definirano je hitrostno polje s potencialom tex2html_wrap_inline1168 in hitrostmi na robu q = v. Ce skozi rob ne tece fluid, je to ovira, za fluid. Pogoj za obtekanje je torej podan z normalno hitrostjo na povrsini:
 equation405
Telesa, ki jih fluid obteka imajo na ograji q = 0 in smer obhoda v smeri urinega kazalca.

  figure408
Figure 2: Obtekanje telesa

Primer obtekanja telesa ze slika 2. Ograja, v kateri zajamemo fluid pa ima smer krivulje naproti obtekajocem trikotniku. Ce zelimo videti gibanje fluida, moramo na ograjo podati hitrost q na elementih v katerih fluid priteka in odteka, s tem, da upostevamo enacbo 13. Za elemente skozi pretoka ni podamo nicno hitrost tex2html_wrap_inline1236. Laplaceovo enacbo 10 resujemo z znanimi hitrostmi na vseh elementih. Sistem enacb 30 se tako poenostavi na
 equation429
kjer je tex2html_wrap_inline1238 vektor izracunan z enostavnim mnozenjem matrike tex2html_wrap_inline1220 in vektorja hitrosti na ograji tex2html_wrap_inline1242.



Leon Kos
Thu Nov 20 12:56:40 CET 1997