30.06.2013 Views

ดาวน์โหลด All Proceeding - AS Nida

ดาวน์โหลด All Proceeding - AS Nida

ดาวน์โหลด All Proceeding - AS Nida

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!