Tuesday, August 23, 2016

Шийдлийн алгоритм дахь Шингэний эзлэхүүний арга (Volume of Fluid Method in SOLA)

Шийдлийн алгоритм дахь (SOLA)

Шингэний эзлэхүүний арга (VOF)

Шингэний чөлөөт гадаргууг ойролцоолох нь тооцон бодох гидравликийн нэгэн чухал асуудал байсаар ирсэн. Ойрын хэдэн жил байх ч болно. Чөлөөт гадаргууг тодорхойлох батлагдсан аргуудыг доорхи зургаас хараарай. Хамгийн анхны арга нь Маркер болон нүдний арга юм. Энэ аргад Лагранжийн координатаар явах маркар буюу цэгнүүд нь шингэний чөлөөт гадаргууг тодорхойлоход хэрэглэгддэг. Үүний дараа гарч ирсэн одооны ихэнх арилжааны програмуудын /Ansys, Solidworks, Xflow, 3Dflow гэх мэт/ ашигладаг арга бол шингэний эзлэхүүний арга /Volume of fluid/ юм. Энэ удаад хамгийн анх боловсруулагдсан шингэний эзлэхүүний аргын тухай авч үзье. Энд тооцооны төгсгөлөг нүдэнд байх шингэнийг эзлэхүүний хувиар илэрхийлж түүнийгээ шингэний фракц гэж авч үзнэ. Шингэний эзлэхүүний фракцийн хугацаан дахь өөрчлөлтийг энгийн адвекцийн тэгшитгэлээр шийднэ.
Гэхдээ энэ адвекцийн тэгшитгэлийн төгсгөлөг ялгаварын ойролцоолол нь нарийвчлал муутай байдаг учир давхар донор ба хүлээн авагчийн схемийг хэрэглэнэ.
Шингэний чөлөөт гадаргууг ойролцоолох тооцонгийн аргууд.

Wednesday, August 17, 2016

Эрдмийн ажлын зураг, графикийг Нюфлотоор дүрсжүүлэх /Scientific visualization and graphics with Gnuplot/

Эрдмийн ажлын зураг, графикийг Нюфлотоор дүрсжүүлэх, зурах
Хоёрдугаар хэсэг

Нэгдүгээр хэсгийг эндээс уншаарай.
Сэтгүүл болгон л өөрийн тодорхой шаардлагуудыг хэвлэгдэхээр илгээгдэн ирж байгаа өгүүллүүдэд тавьдаг. Хамгийн гол нь тэдгээр шаардлагыг сайн биелүүлэхийн тулд ямар ч асуудлыг шийдвэрлэх гарцтай байх хэрэгтэй. Бид графикийг ихэвчлэн Эксэл дээр бужигнуулдаг. Энэ ч амархан. Харин хэд хэдэн графикийг хооронд нь тэнхлэгээр нь нийлүүлэх гэх мэт ахисан түвшний ажлууд гараад ирвэл суугай майзаг болоод эхлэнэ. Нэг ашиглаад суричхвал Нюфлот арвин чадвартай. Энэ хэсгээр зарим нэг чухал үйлдлүүдийг хийж сурсанаар та бүгдтэй хуваалцая.

Муруйд тохирсон регрессийн тэгшитгэл зохиож зурах
Тоон өгөгдөлд тохирох шугаман болон өндөр эрэмбийн регрессийн тэгшитгэлийг эксэл дээр ч хийгээд коэффициентуудыг нь олж болно. Гэхдээ Нюфлот нь өөр дээрээ итераци хийж хамааралтай бүх хэмжигдэхүүнүүдийг тодорхойлдог. 
Энэ жишээнд postscript терминалыг ашигласан болно. Өндөр чанартай вектор зураглалуудыг ихэвчлэн postscript хийж болдог. 

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
# this script for compare remaining ice area
 reset
 set term postscript landscape enhanced font "Times-New-Roman, 30"  #size 700,500 #500,400  
 set output 'Compare.esp'

#
 set xlabel "{/=30 {/Times-Italic Time } (min)}"
 set ylabel "{/=30 {/Times-Italic Remaining ice area} (percent)}"
 set xrange [0:72]
 set yrange [0.5:1.01]
 set key top right 
 f(x) = a*x + b # fit line 1
 fit f(x) 'ice_example.dat' u 1:2 via a, b
 title_f(a,b) = sprintf('f(x) = %.2fx + %.2f', a, b)
 
 #f(x) = a*x **2 + b *x + e # fit line 1
 #fit f(x) 'ice_example.dat' u 1:2 via a, b, e
 #title_f(a,b,e) = sprintf('f(x) = %.2fx **2 + %.2fx + %.2f', a, b, e)
 g(x) = c*x + d # fit line 2
 fit g(x) 'ice_mass.dat' u 1:2 via c, d
 title_g(c,d) = sprintf('g(x) = %.2gx + %.2g', c, d)
 
 #f1(x) = E*exp(I*x**2) # fit line 3
 #fit f1(x) 'ice_mass.dat' u 1:2 via E, I
 #title_f1(E,I) = sprintf('f1(x) = %.2f1*exp(%.2f1x**2)', E, I)
 
 plot 'ice_example.dat' u 1:2 w lp pt 8 pi -40 lc rgb 'red' lw 7 title 'Free surface LBM' ,\
   'ice_mass.dat' u 1:2 w p pt 5 lc rgb 'black' lw 9 title 'Experiment',\
    f(x) w l ls 10 lw 4 title 'Trend line for LBM',\
    g(x) w l ls 9 lw 4 title 'and for experiment' 
#    f1(x) t title_f1(e,i), title_f(a,b),title_g(c,d)
   
 unset multiplot
 set term win
 unset output