Title 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue

Size: px

Start display at page:

Download "Title 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue"

ゆみかいりぐら
5 years ago
Views:

A2 (2017), 73 Issue Date 2017 URL http://hdl.handle.

1 Title 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue Date 2017 URL Right 2017 公益社団法人土木学会 Type Journal Article Textversion publisher Kyoto University

2 A2, Vol. 73, No. 2 Vol. 20, I_143-I_152, ( ) toriu.daisuke.8v@kyoto-u.ac.jp ushijima.satoru.3c@kyoto-u.ac.jp Key Words : mixture model, compressible fluid, moving solid, thermal fluid-solid interaction 1. β T β T 1 Yamamoto 1, 2) Qi 3) 4) 5) MICS 6) TCUP 7) 4) CFL (Courant-Friedrichs-Lewy) 8) 4) CFL 4) 4) 9) I_143

10) 4) Lee & Ha 11) CFL Ra Ra 10 6 β T β T 0.3 2.

3 10) 4) Lee & Ha 11) CFL Ra Ra 10 6 β T β T (1) 1 ( 1) 4) ) 1 1 T f1 T s T f1 T s 2 T f2 = T s 1 2 T f1 T f2 T T f1 T s 10) 7) ( ) 4) 4) (2) 1 ( ) ( ) 1 2 ( 1 ) 1 2 5) I_144

4 (ρu i ) (ρe) ρ + (ρu i) x i = 0 (1) + (ρu iu j ) + (ρeu j) = σ ij + ρf i (2) = σ ij u i (3) t x i f i x i ρ u i σ ij e (2) (3) 5) 1 (3) 1 2 (1) (2) (3) ρ σ ij u i e k w k,i 5) w k,i 0 w k,i k u k,i u i w k,i u k,i u i 5) ρ u i ρ = ϕ f1 ρ f1 + ϕ f2 ρ f2 (4) u i = ϕ f1ρ f1 u f1,i + ϕ f2 ρ f2 u f2,i ϕ f1 ρ f1 + ϕ f2 ρ f2 = 1 ρ (ϕ f1ρ f1 u f1,i + ϕ f2 ρ f2 u f2,i ) (5) f1 1 f2 2 ϕ f1 ϕ f ϕ f1, ϕ f2 1 ϕ f1 + ϕ f2 = 1 12) 12) 1 2 σ ij 4) σ ij = ϕ f1 σ f1,ij + ϕ f2 σ f2,ij pδ ij + τ ij (6) ( ui τ ij µ f + u ) j x i 2 3 µ u m f δ ij (7) x m δ ij p τ ij µ f e T f1 T f2 (= T ) e = Σ kϕ k ρ k C V,k T k Σ k ϕ k ρ k T k = C V,f Σ k ϕ k ρ k Σ k ϕ k ρ k C V,f T (8) k k k = f1, f2 1 2 C V,f1 = C V,f2 = C V,f 1 2 e T p T f1 T f2 (= T ) p = ϕ f1 p f1 + ϕ f2 p f2 = ϕ f1 ρ f1 (γ f1 1)C V,f1 T f1 +ϕ f2 ρ f2 (γ f2 1)C V,f2 T f2 (ϕ f1 ρ f1 + ϕ f2 ρ f2 )(γ f 1)C V,f T = ρ(γ f 1)C V,f T (9) 1 2 γ f1 = γ f2 = γ f (1) 1 1 (ρ f1 C V,f1 T f1 ) (ρ s C V,s T s ) = q f1,j (10) = q s,j (11) q f1,j q s,j 1 I_145

5 5) 1 (ρc V ) m T = q m,j (12) (ρc V ) m = ϕ f1 ρ f1 C V,f1 +ϕ s ρ s C V,s q m,j T f1 T f2 (= T ) q m,j T q m,j = ϕ f1 q f1,j + ϕ s q s,j = ϕ f1 λ f1 T f1 ϕ s λ s T s (ϕ f1 λ f1 + ϕ s λ s ) T = λ m T (13) λ m 1 λ m = ϕ f1 λ f1 + ϕ s λ s (3) 2. (1) 4) 4) 7) Q n Q Q Q n+1 a) 1 2 ρ + (ρu j) = 0 (14) Adv (ρu i ) + (ρu iu j ) = 0 (15) Adv (ρe) + (ρeu j) = 0 (16) Adv Euler MUSCL TVD 13) u i (15) (ρu i ) u i = 1 ρ m [ϕ f1 (ρu i ) + ϕ s (ρu i ) s ] (17) ρ m = ϕ f1 ρ + ϕ s ρ s 1 (9) p ( 2) ϕ s = ϕ f2 = 1 p = p n b) 1 2 (1) (2) (3) T 1 (12) T (1) (2) (3) ρ = 0 (18) Diff (ρu i ) = τ ij (19) Diff (ρe) u i = τ ij (20) Diff1 2 I_146

6 Euler u i T T T t 1 = ρ C V,f u i u i ρ t { (τ ij u i ) = τ ij (21) ρ (u 2 i u 2 2 t i ) } (22) 1 (12) T T T T 1 T = (λ ) t (ρc V ) m m (ρc V ) m (23) (ρc V ) m = ϕ f1 ρ f1c V,f1 + ϕ s ρ s C V,s (24) ρ f2 (24) ϕ f1 ρ f1 ϕ f1 ρ f1 = ρ ϕ f2 ρ f2 (25) (23) T T T = (1 ϕ sc )Θ + ϕ sc T sc (26) t Θ (23) ϕ sc T sc 1 2 7) k p k = γ k 1 ρk C P,k ρ k C (Tk Tk ) (27) P,kµ J,k + 1 p γ k k = f1 f2 C P,k k µ J,k k Joule-Thomson p = ϕ f1 p f1 + ϕ f2 p f2 1 2 µ J,k = 0 p p γ f 1 γ f (ϕ f1 ρ f1 + ϕ f2 ρ f2)c P,f (T T ) c) = γ f 1 γ f ρ C P,f (T T ) (28) 1 2 (1) (2) (3) (ρu i ) (ρe) ρ = 0 (29) Acous = p + ρf i (30) Acous x i = p u i (31) Acous x i 4) TCUP 7) 1 p n+1 p ρ a 2 = t ( 1 p n+1 ) x i ρ t + u i x i (32) a a = (γ f p )/ρ (32) CFL 4) 3. Lee & Ha 11) 11) 11) 11) 3 T h T c I_147

7 (a) t = 12 [s] (b) t = 60 [s] 3 (c) t = 72 [s] (d) 5 (Ra = 10 4 ) 4 11) T = T h T c = 5 [K] (β T = 0.017) 3 g x 2 ρ s /ρ f0 = 10 C s /C P,f = 10 λ s /λ f = 50 ρ f0 C P,f (T h +T c )/2 P r 0.70 L Ra Ra Ra = gβ T L3 ν f0 α f0 (33) ν f0 α f0 Ra = Ra = Lee & Ha 11) (a) t = 10 [s] (b) t = 20 [s] (c) t = 24 [s] (d) 6 (Ra = 10 6 ) Ra T/10 5(a) 6(a) I_148

8 1 Nu h U 1 U 2 Ra = 10 4 Nu h U 1 U (a) Ra = 10 3 (b) Ra = 10 4 Ra = 10 6 Nu h U 1 U ) Nu h (a) Ra = 10 5 (b) Ra = ( T/10, T = 5 [K]) T 7 T/ (1) 2. (2) T f1 T s u k,i u i (w k,i = 0) Nu h x 1 = L/2 U 1 x 2 = L/2 U 2 Nu h 1 Nu h = 0 θ X 2 1 (34) X2=0dX θ X i θ = (T T c )/ T X i = x i /L gβl T 1 Ra = Ra = 10 6 Nu h Ra Present Lee & Ha 11) Nu h U 1 U 2 CFL C a 14) C a C a,max { u1 + a C a = max t, u } 2 + a t x 1 x 2 (35) t x i C a,max Ra = CFL t Nu h Lee & Ha 11) 2 2 T = 100 [K] (β T = 0.286) Ra = T/10 8 T = 100 [K] I_149

9 8 ( T = 100 [K]) ( T/10) 3 T ρ/ ρ 0 T [K] Max Min T = 5 [K] T = 5 [K] Nu h 6.31 T = 5 [K] C a,max ρ ρ 0 ρ/ ρ 0 3 ρ/ ρ 0 3 T = 5 [K] 1% T = 100 [K] 15.7% T ( r h ) T h T c T h = 400 [K] T c = 300 [K] L [m] r t L/3 r h 3L/40 g 9.8 [m/s 2 ] ω ω ω = t t b π (0 t t b ) (36) ω = π (t b < t) (37) t b = 3.0 [s] t b 0 12) (Fe) ( ) I_150

10 4 ( ) ( ) Fluid (air) Solid (Fe) ρ [kg/m 3 ] λ [W/(m K)] C V [J/(kg K)] t [s] t = [s] Ra Ra = V ( u u2 2 ) Re t = [s] Re (T h T c )/ (b) (c) t = [s] ( 10 (c)) t = [s] C a ,000 7,000 x2 x1 x2 x1 x2 x1 (a) t = 0.12 [s] (b) t = 3.00 [s] (c) t = [s] ) 11) I_151

11 β T % 7,000 3, JSPS JP16K ) Yamamoto, S., Niiyama, D. and Shin, B. R.: A numerical method for natural convection and heat conduction around and in a horizontal circular pipe, Int. J. Heat Mass Transfer, Vol. 47, pp , ) Yamamoto, S.: Preconditioning method for condensate fluid and solid coupling problems in general curvilinear coordinates, J. Compt. Phys., Vol. 207, pp , ) Qi, S., Furusawa, T. and Yamamoto, S.: A numerical method applied to forced and natural convection flows over arbitrary geometry, Int. J. Heat Mass Transfer, Vol. 85, pp , ), CFL, A2( ), Vol. 71, No.2, pp. I 213 I 222, ), pp ,, ),, 3, B, Vol. 64, No.2, pp , ), ( CCUP -TCUP - ), (B ), Vol. 69, No.678, pp , ) P. J. Roache(, ), pp ,, ),,,,, Vol. 19, No. 6, pp , ),, A2( ), Vol. 72, No.2, pp. I 179 I 186, ) Lee, J. R. and Ha, M. Y.: Numerical simulation of natural convection in a horizontal enclosure with a heat-generating conducting body, Int. J. Heat Mass Transfer, Vol. 49, pp , ),, 3,, Vol. 51, pp , ) Yamamoto, S. and Daiguji, H.: Higher-order-accurate upwind schemes for solving the compressible Euler and Navier-Stokes equations, Computers & Fluids, Vol. 22, pp , ),, (B ), Vol. 69, No. 682, pp , ( ) COMPUTATIONAL METHOD FOR THERMAL INTERACTIONS BETWEEN COMPRESSIBLE FLUID AND MOVING SOLID BASED ON MIXTURE MODEL Daisuke TORIU and Satoru USHIJIMA In this paper, a computational method was proposed for thermal interactions between compressible fluids and moving solids with different physical properties based on the mixture model. In the present method, physical properties of solids are considered in computational stages for heat conduction and velocities of multiphase fields. The present method was applied to natural convection in the cavity containing a square solid and obtained averaged Nusselt numbers on the heated wall were in good agreement with reference results. In addition, applicability of the present method was confirmed for compressible high buoyancy flows and heat transfer around a rotating triangular solid. I_152

Title 鉛直配置された水平 2 円柱周りの自然対流に対する圧縮性流体と固体の熱連成計算手法の適用性 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (

Title 鉛直配置された水平 2 円柱周りの自然対流に対する圧縮性流体と固体の熱連成計算手法の適用性 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2016), 72 Issue Date 2016 URL http://hdl.handle.net/2433/229149