ดาวน์โหลด All Proceeding - AS Nida
ดาวน์โหลด All Proceeding - AS Nida
ดาวน์โหลด All Proceeding - AS Nida
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
ดังนี้<br />
อัตราการลดลงของพลังงานจลนของคาความปนปวนเขียนได<br />
∂ ∂ ⎛⎛ µ ⎞ t ∂ε⎞<br />
( ρε ui<br />
) = ⎜⎜µ + ⎟ ⎟<br />
∂xi ∂x ⎜<br />
i σ k x ⎟<br />
⎝⎝ ⎠∂<br />
i ⎠<br />
่<br />
้<br />
⎛ ε ⎞ε<br />
ε<br />
+ ⎜Cε1 ( Gk + Cε3Gb) ⎟ −ρCε2<br />
⎝ k ⎠k<br />
k<br />
โดยที<br />
⎛∂uj⎞ Gk =−ρu′ i u ′<br />
j ⎜ ⎟<br />
⎝ ∂xi<br />
⎠<br />
(7)<br />
µ t<br />
Gb =−gi<br />
ρ PrT<br />
∂ρ<br />
∂xi<br />
(8)<br />
สําหรับคาคงที่ตางๆในสมการจะมีคาคงที่ของ<br />
Launder and<br />
Sharma [3] ดังนี<br />
C µ =0.09, 1 Cε =1.44, 2 Cε =1.92, 3 Cε =1, σ k =1, σ ε =1.3, σ T =<br />
0.9<br />
แบบจําลองความปนปวนชนิด<br />
Standard k − ω SST<br />
(SST, Shear Stress Transport Turbulence Model)<br />
แบบจําลอง k − ω SST ของ Menter [1] สมการพลังงาน<br />
จลนของความปนปวน<br />
(k) สามารถเขียนไดดังนี้<br />
∂ ∂ ⎛ ∂k⎞<br />
( ρ kui)<br />
= ⎜( µ + σ kµ t)<br />
⎟<br />
∂xi ∂xi ⎝ ∂xi<br />
⎠<br />
∂u<br />
j *<br />
−ρuu ′ ′<br />
i j −ρβ<br />
kω<br />
∂x<br />
i<br />
สมการอัตราการลดลงของพลังงานจลนของคาความปนปวน<br />
จําเพาะ ( ω) สามารถเขียนไดดังนี้<br />
∂ ∂ ⎛ ∂ω⎞<br />
( ρω ui)<br />
= ⎜( µ + σ ωµ<br />
t)<br />
⎟<br />
∂xi ∂xi ⎝ ∂xi<br />
⎠<br />
α ∂u<br />
j<br />
2<br />
− ρuu ′ ′<br />
i j −ρβω<br />
vt ∂xi<br />
1 ∂k∂ω + 21 ( −F1<br />
) ρσ ω,2<br />
ω ∂xj ∂xj<br />
โดยที่คา<br />
Eddy viscosity แสดงไดดังนี้<br />
k 1<br />
µ t = ρ<br />
ω ⎡ 1 ΩF2<br />
⎤<br />
max ⎢ , *<br />
⎣αa1ω⎥ ⎦<br />
(6)<br />
(9)<br />
(10)<br />
(11)<br />
โดยที่<br />
Ω= ΩijΩ ij<br />
(12)<br />
164<br />
1 ⎛∂u∂u⎞ i j<br />
Ω ij = ⎜ −<br />
2 ⎜<br />
⎟<br />
x j x ⎟<br />
⎝∂ ∂ i ⎠<br />
(13)<br />
F = tanh Φ (14)<br />
2 ( )<br />
2 2<br />
⎡ k 500µ<br />
⎤<br />
Φ 2 = max ⎢2 , 2 ⎥ (15)<br />
⎣ 0.09ωy<br />
ρωy<br />
⎦<br />
โดยที่พจนความเคนเรยโนลด<br />
ซึ่งสามารถเขียนเปนสมการ<br />
ตามสมมุติฐานของ Boussinesq [2] ดังสมการที่<br />
(3) ซึ่งคาคงที่ในสมการ<br />
*<br />
ไดแก β , σk, σωหาไดจากสมการที่ (16) โดยที่คา<br />
θ เปนคาคงที่<br />
ใดๆที่ใชในสมการ,<br />
θ1 เปนคาคงที่ที่มาจากแบบจําลองชนิด<br />
k − ω ,<br />
θ2 เปนคาคงที่ที่มาจากแบบจําลองชนิด<br />
k − ε ซึ่งคา<br />
( θ ) F ( θ) ( 1 F )( θ )<br />
่<br />
=<br />
โดยที<br />
1 + − 1 2<br />
(16)<br />
4 ( )<br />
F = tanh Φ (17)<br />
1 1<br />
⎡ ⎛ k 500µ ⎞⎛<br />
4ρk<br />
⎞⎤<br />
Φ 1 = min ⎢max ⎜ , 2 ⎟⎜ + 2 ⎟⎥(18)<br />
⎝0.09ωy ρy ω⎠ ⎜σω,2Dωy ⎟<br />
⎢⎣ ⎝ ⎠⎥⎦<br />
⎡<br />
+ 1 ∂k∂ω ⎤<br />
−20<br />
Dω<br />
= max ⎢2 ρ ,10 ⎥<br />
(19)<br />
⎢ σ ω ⎣<br />
∂x ,2 j ∂x<br />
ω<br />
j ⎥⎦<br />
สําหรับคาคงที่อื่นๆในสมการมีดังนี้<br />
β 1 =0.075, β 2 =0.0826, a 1 =0.31,<br />
*<br />
α =1<br />
2.2 วิธีดําเนินการศึกษา<br />
งานวิจัยนี้ไดดําเนินการศึกษาผลของการจําลองเชิงตัวเลขที่ได<br />
จากการวิเคราะหผลการไหลของน้ําที่ไหลผานเสน<br />
Profile ใบของลอ<br />
กังหันน้ําผลิตไฟฟาขนาดเล็กแบบแกนตั้ง<br />
โดยทําการประเมิน<br />
ประสิทธิภาพของความแมนยําและผลของการลูเขาของคําตอบที่รวดเร็ว<br />
และมีถูกตองของแบบจําลองความปนปวนที่เกิดจากสภาวะการไหลทั้ง<br />
สองชนิดที่เลือกใช<br />
ซึ่งในงานวิจัยนี้ไดเลือกชนิดของแบบจําลองไว<br />
2<br />
ชนิด คือ แบบจําลองชนิด k − ε Model และ แบบจําลองชนิด k − ω<br />
SST เพื่อนําผลการจําลองเชิงตัวเลขที่ไดมาทําการวิเคราะหและทําการ<br />
เปรียบเทียบผลที่ได<br />
จากผลการทดสอบจริงจากตนแบบที่ไดออกแบบ<br />
และสรางขึ้น<br />
ดังรูปที่<br />
1