인쇄본 - 10졸업논문_배세욱_강내탄도에 사용되는 Ergun식에 대한 수치적 보정연구_초록 수정.hwp

Size: px

Start display at page:

Download "인쇄본 - 10졸업논문_배세욱_강내탄도에 사용되는 Ergun식에 대한 수치적 보정연구_초록 수정.hwp"

규철 황
5 years ago
Views:

4 - I -

5 - II -

6 max Maximum exponent of the velocity in pressure drop equation - III -

7 - IV -

8 - V -

9 - VI -

10 - VII -

11 - VIII -

12 - 1 -

13 하나하나의 입자에 대하여 Lagrangian 접근법을 사용하여 지배방정식을 풀게 된다. 이 때, 기상과 고상의 지배방정식의 물리적 정보를 결합시켜 줄 수 있는 관계식이 필요하게 된다. 대표적으로 Ergun[1]의 실험식이 사용되고 있으며, 그 식은 아래와 같다. (1) 여기서 U는 공탑속도(Superficial velocity), L은 이상유동 영역의 길이, 은 기 공률(Porosity), 는 입자의 평균 입자 지름이다. 그림 1-2 Packed bed 구성 그림 1-3 입자 배치 형태 그림 1-4 Discrete Particle Model 해석 방법 [8] - 2 -

14 [8] - 3 -

15 표 1-1 Reddy의 CFD 해석결과오차 레이놀즈수 오차율 (%)

16 표 2-1 이상유동수치적해석모델 [16] - 5 -

17 그림 2-1 이상유동해석기법 [16] - 6 -

18 - 7 -

19 그림 2-2 연구수행프로세스 - 8 -

표 2-2 CATIA 모델링설정 기공률 0.70 0.75 0.80 Tube 직경 [m] 0.702 입자크기 [m] 0.01 직경입자비 7.

20 표 2-2 CATIA 모델링설정 기공률 Tube 직경 [m] 입자크기 [m] 0.01 직경입자비 7.02 Inlet 길이 [m] 0.03 Bed 길이 [m] Outlet 길이 [m] 0.03 입자개수 그림 2-3 CATIA 모델링 - 9 -

21 표 2-3 Meshing 기법 격자형태세부설정특징 Volume meshing Polyhedral mesh Tetrahedral mesh 의결합형태 Extra meshing Prism layer mesh Extrude mesh Wall 근처에분포 Outlet 이후분포 그림 2-4 Polyhedral / Prism layer / Extrude mesh

22 그림 2-5 여러조건의난류유동 그림 2-6 무차원화된난류유동의구분

23 그림 2-7 wall 근처에서의유체해석법

그림 2-8 격자생성결과 표 2-4 격자생성결과정리 기공률 입자개수 격자수 0.

24 그림 2-8 격자생성결과 표 2-4 격자생성결과정리 기공률 입자개수 격자수 만개 만개 만개

25 표 2-5 수치해석모델설정 해석분류 Space Time Material Flow Equation of State 해석방법 Three Dimensional Steady Gas Segregated Flow Constant Density Turbulent Reynolds-Averaged Navier-Stokes Viscous Regime K-Epsilon Turbulence Realizable K-Epsilon Two-Layer Two-Layer All y+wall Treatment

26 - 15 -

27 그림 2-9 All y+ Treatment

28 보정연구를위한데이터수집방법으로 CFD 해석결과를사용하기위해서는실험결과와의비교를통한검증과정을거쳐신뢰성을갖출필요가있다. 따라서 R.K. Reddy의연구결과와의비교검증을통해서본연구에서사용하는 CFD 해석기법의신뢰성을확보하였다. Reddy는 0.1~10000의레이놀즈수범위에대하여 Ergun식과 CFD 해석및실험결과를비교하였다. 그결과, CFD 해석결과와실험결과는비교적일치한반면, Ergun식계산결과와 CFD 해석결과를비교하였을때차이가발생하였다. 레이놀즈수가 0.1~10 영역에서는비교적일치, 10~500 영역에서는매우일치, 500이상영역에서는오차가점차증가하는결과를얻었다. 따라서레이놀즈수가 10~500인영역에서 CFD 해석을수행하여 Ergun식계산결과와비교하여 CFD 해석의신뢰성을확인할수있다. Reddy와같은해석조건에서레이놀즈수가 300일때의해석을수행한결과, 표 2-6과같이오차율 3.31% 로매우일치함을확인하였다. 표 2-6 CFD 해석검증결과 해석방법및오차율 결과값 (Re = 300) Ergun식계산결과 Pa CFD 해석결과 Pa 오차율 3.31 %

29 CFD 해석에사용된컴퓨터는 HPC(High Performance Computing) 이가능한 워크스테이션급 PC 2 대이다. 주요구성사항은표 2-7 과같다. 표 2-7 하드웨어구성 순번 CPU 사양 Core 수 RAM 용량 1st PC 2nd PC Intel Zeon X5550 (2.66Ghz) AMD Phenom II-X6 (3.2Ghz) 8 (4core X 2cpu) 48 Gbyte 6 (6core X 1cpu) 16 Gbyte 만 만

30 그림 3-1 모델링구성

31 그림 3-2 직경입자비에따른압력강하그래프 ( 기공률 0.4 case)

32 Dimensional Pressure Drop Porosity 0.40, Re=10 [Eisfeld] Porosity 0.40, Re=100 [Eisfeld] Porosity 0.40, Re=1000 [Eisfeld] Porosity 0.73, Re=1000 Porosity 0.73, Re=10000 Porosity 0.73, Re= Diameter Particle Ratio 그림 3-4 직경입자비변화에따른레이놀즈수별압력강하비교 ( 기공률 0.73)

33 표 3-1 직경입자비에따른압력강하량비교 (Re=1,000) 직경입자비 무차원압력강하 입자수 격자수 만개 만개 만개 만개 13.0(Infinity) 만개 표 3-2 직경입자비에따른압력강하량비교 (Re=10,000) 직경입자비 무차원압력강하 입자수 격자수 만개 만개 만개 만개 13.0(Infinity) 만개 표 3-3 직경입자비에따른압력강하량비교 (Re=100,000) 직경입자비 무차원압력강하 입자수 격자수 만개 만개 만개 만개 13.0(Infinity) 만개

34 - 23 -

35 표 4-1 수치해석초기조건 해석초기조건 입구정압력 20 atm 입구정온도 300 K Re Velocity Inlet 입구속도 m/s m/s m/s m/s m/s Outlet 출구정압력 출구정온도 20 atm 300 K 유체밀도

36 그림 4-1 =0.70, Total pressure

37 그림 4-2 =0.70, Total pressure

38 그림 4-3 =0.75, Total pressure

39 그림 4-4 =0.75, Total pressure

40 그림 4-5 =0.70, Velocity

41 그림 4-6 =0.70, Velocity

42 그림 4-7 =0.75, Velocity

43 그림 4-8 =0.80, Velocity

44 그림 4-9 Ergun 식계산결과와 CFD 해석결과의압력강하비교

45 표 4-2 각 Case 별압력강하량 Porosity Re

46 표 4-3 각 Case 별오차율 Porosity Re % 212.7% 2.767% 3.523% 454.4% % 231.5% 302.6% 404.1% 501.4% % 251.5% 327.5% 433.5% 525.3%

47 max max max

48 표 5-1 전체 Ergun 식중 Laminar flow 항이차지하는비율 Porosity Re % 0.51% 0.26% 0.05% 0.02% % 0.43% 0.21% 0.04% 0.02% % 0.34% 0.17% 0.03% 0.02% 표 5-2 각 Case 별 CFD 해석을통해구해진 값 Porosity Re

49 그림 5-1 각 Case 별 Re 에따른 결과

50 max max max 그림 5-2 Re' 에따른 결과

51 max 그림 5-3 Re'' 에대한 결과

52 (5) where,,

53 그림 5-5 Ergun 식, 보정된 Ergun 식, CFD 해석결과비교 표 5-3 Ergun eqaution의보정전보정후의오차율비교 Re Porosity 0.70 보정전 120.3% 212.7% 2.767% 3.523% 454.4% 보정후 2.84% 1.36% 5.97% 6.13% 1.70% 0.75 보정전 131.7% 231.5% 302.6% 404.1% 501.4% 보정후 0.44% 0.17% 5.56% 1.98% 3.53% 0.80 보정전 139.0% 251.5% 327.5% 433.5% 525.3% 보정후 6.99% 3.08% 3.01% 4.09% 0.29%

54 그림 5-6 표 5-4 Ergun 식과 modified Ergun 식의오차율비교 Porosity Re 보정전 19.41% 23.44% 72.22% 148.6% 176.2% 보정후 6.5% 4.62% 10.27% 6.64% 1.46%

55 - 44 -

56 [1] Ergun, S. "Fluid flow through packed columns", Chem Eng Prog vol.48, pp.89-94, 1953 [2] FC Blake. The resistance of packing to fluid flow. Trans Am Inst Chem Eng 14:415, 1922 [3] J Kozeny. U"N ber Kapillare Leitung des Wassers im Boden. Sitzungsber Akad Wiss Wien 136:271, 1927 [4] PC Carman. Fluid flow through granular beds. Trans Inst Chem Eng 15:150, 1937 [5] SP Burke, WB Plummer. Gas flow through packed columns. Ind Eng Chem 20: , 1928 [6] WE Ranz. Friction and transfer coefficients for single particles and packed beds. Chem Eng Prog 48: , 1952 [7] CL Kusik, J Happel. Boundary layer mass transfer with heterogeneous catalysis. AIChE J 8:163, 1962 [8] TR Galloway, BH Sage. A model of the mechanism of transport in packed, distended, and fluidized beds. Chem Eng Sci 25: , 1970 [9] WH Gauvin, S Katta. Momentum transfer through packed beds of various particles in the turbulent flow regime. AIChE J 19: , 1973 [10] H Darcy. Les Fontaines publiques de la ville de Dizon. Paris, Dalmont, 1856 [11] LE Brownell, DL Katz. Flow of fluids through porous media I. Single homogeneous fluids. Chem Eng Prog 43: , October 1947 [12] L Gibilaro, RDI Delice, SP Waldram, PU Foscolo.Generalized friction factor and drag coefficient correlations for fluid.particle interactions. Chem Eng Sci 40: , 1985 [13] H.P.A. Calis et al, "CFD modelling and experimental validation of pressure drop and flow profile in a novel structured catalytic reactor packing", Chemical Engineering Science Vol.56 pp , 2001 [14] Theodoros Atmakidis et al, "CFD-based analysis of the wall effect on the pressure drop in packed beds with moderate tube/particle diameter ratios in the laminar flow regime", Chemical Engineering Journal vol.155, pp ,

57 [15] Rupesh K. Reddy et al, "CFD modeling of pressure drop and drag coefficient in fixed and expanded beds", chemical engineering research and design vol.86, pp , 2008 [16] M.A. van der Hoef et. al. Numerical Simulation of Dense Gas-Solid Fluidized Beds: A Multiscale Modeling Strategy, Annu. Rev. Fluid Mech. vol , [17] Eisfeld, B. and Schnitzlein, K., The influence of confining walls on the pressure drop in packed beds, Chem Eng Sci, 56: [18] J.A. Tallmadge, Packed bed pressure drop - an extension to higher Reynolds numbers, Alche J.16, 1092~1093, 1970 [19] Wen-Ching Yang, Handbook of Fluidization and Fluid-Particle systems, Marcel Dekker Inc, New York Basel, p , 2003 [20] STAR-CCM+Userguide , CD-adapco

58 - 47 -

82 제 1 발표장 (2 일금 ) CFD 응용 [V] 이남훈 1*, 류태광 2 NUMERICAL VERIFICATION OF SHAKE TABLE TEST FOR THE LIQUID STORAGE TANK N. Lee and T. Yoo 1.,.,,,.,. Baek e

82 제 1 발표장 (2 일금 ) CFD 응용 [V] 이남훈 1*, 류태광 2 NUMERICAL VERIFICATION OF SHAKE TABLE TEST FOR THE LIQUID STORAGE TANK N. Lee and T. Yoo 1.,.,,,.,. Baek e 82 제 1 발표장 (2 일금 ) 이남훈 1*, 류태광 2 NUMERICAL VERIFICATION OF SHAKE TABLE TEST FOR THE LIQUID STORAGE TANK N. Lee and T. Yoo 1.,.,,,.,. Baek et al.[1]. CFD FLUENT LS-DYNA. 2. 2.1 CFD FLUENT, LS-DYNA. FLUENT,