이 연구에서 FLOW 3D 전산 유체 역학(CFD) 소프트웨어를 사용하여 파키스탄 Mirani 댐 방수로에 대한 에너지 소산 옵션으로 미국 매립지(USBR) 유형 II 및 USBR 유형 III 유역의 성능을 추정했습니다. 3D Reynolds 평균 Navier-Stokes 방정식이 해결되었으며, 여기에는 여수로 위의 자유 표면 흐름을 캡처하기 위해 공기 유입, 밀도 평가 및 드리프트-플럭스에 대한 하위 그리드 모델이 포함되었습니다. 본 연구에서는 5가지 모델을 고려하였다. 첫 번째 모델에는 길이가 39.5m인 USBR 유형 II 정수기가 있습니다. 두 번째 모델에는 길이가 44.2m인 USBR 유형 II 정수기가 있습니다. 3번째와 4 번째모델에는 길이가 각각 48.8m인 USBR 유형 II 정수조와 39.5m의 USBR 유형 III 정수조가 있습니다. 다섯 번째 모델은 네 번째 모델과 동일하지만 마찰 및 슈트 블록 높이가 0.3m 증가했습니다. 최상의 FLOW 3D 모델 조건을 설정하기 위해 메쉬 민감도 분석을 수행했으며 메쉬 크기 0.9m에서 최소 오차를 산출했습니다. 세 가지 경계 조건 세트가 테스트되었으며 최소 오류를 제공하는 세트가 사용되었습니다. 수치적 검증은 USBR 유형 II( L = 48.8m), USBR 유형 III( L = 35.5m) 및 USBR 유형 III 의 물리적 모델 에너지 소산을 0.3m 블록 단위로 비교하여 수행되었습니다( L= 35.5m). 통계 분석 결과 평균 오차는 2.5%, RMSE(제곱 평균 제곱근 오차) 지수는 3% 미만이었습니다. 수리학적 및 경제성 분석을 바탕으로 4 번째 모델이 최적화된 에너지 소산기로 밝혀졌습니다. 흡수된 에너지 백분율 측면에서 물리적 모델과 수치적 모델 간의 최대 차이는 5% 미만인 것으로 나타났습니다.
In this study, the FLOW 3D computational fluid dynamics (CFD) software was used to estimate the performance of the United States Bureau of Reclamation (USBR) type II and USBR type III stilling basins as energy dissipation options for the Mirani Dam spillway, Pakistan. The 3D Reynolds-averaged Navier–Stokes equations were solved, which included sub-grid models for air entrainment, density evaluation, and drift–flux, to capture free-surface flow over the spillway. Five models were considered in this research. The first model has a USBR type II stilling basin with a length of 39.5 m. The second model has a USBR type II stilling basin with a length of 44.2 m. The 3rd and 4th models have a USBR type II stilling basin with a length of 48.8 m and a 39.5 m USBR type III stilling basin, respectively. The fifth model is identical to the fourth, but the friction and chute block heights have been increased by 0.3 m. To set up the best FLOW 3D model conditions, mesh sensitivity analysis was performed, which yielded a minimum error at a mesh size of 0.9 m. Three sets of boundary conditions were tested and the set that gave the minimum error was employed. Numerical validation was done by comparing the physical model energy dissipation of USBR type II (L = 48.8 m), USBR type III (L =35.5 m), and USBR type III with 0.3-m increments in blocks (L = 35.5 m). The statistical analysis gave an average error of 2.5% and a RMSE (root mean square error) index of less than 3%. Based on hydraulics and economic analysis, the 4th model was found to be an optimized energy dissipator. The maximum difference between the physical and numerical models in terms of percentage energy absorbed was found to be less than 5%.
Abbasi S, Fatemi S, Ghaderi A, Di Francesco S (2021) The effect of geometric parameters of the antivortex on a triangular labyrinth side weir. Water (Switzerland) 13(1). https://doi.org/10.3390/w13010014
Amorim JCC, Amante RCR, Barbosa VD (2015) Experimental and numerical modeling of flow in a stilling basin. Proceedings of the 36th IAHR World Congress 28 June–3 July, the Hague, the Netherlands, 1, 1–6
Asaram D, Deepamkar G, Singh G, Vishal K, Akshay K (2016) Energy dissipation by using different slopes of ogee spillway. Int J Eng Res Gen Sci 4(3):18–22Google Scholar
Celik IB, Ghia U, Roache PJ, Freitas CJ, Coleman H, Raad PE (2008) Procedure for estimation and reporting of uncertainty due to discretization in CFD applications. J Fluids Eng Trans ASME 130(7):0780011–0780014. https://doi.org/10.1115/1.2960953ArticleGoogle Scholar
Chen Q, Dai G, Liu H (2002) Volume of fluid model for turbulence numerical simulation of stepped spillway overflow. J Hydraul Eng 128(7):683–688. 10.1061/共ASCE兲0733-9429共2002兲128:7共683兲 CE
Damiron R (2015) CFD modelling of dam spillway aerator. Lund University Sweden
Dunlop SL, Willig IA, Paul GE (2016) Cabinet Gorge Dam spillway modifications for TDG abatement – design evolution and field performance. 6th International Symposium on Hydraulic Structures: Hydraulic Structures and Water System Management, ISHS 2016, 3650628160, 460–470. 10.15142/T3650628160853
Frizell KW, Frizell KH (2015) Guidelines for hydraulic design of stepped spillways. Hydraulic Laboratory Report HL-2015-06, May
Ghaderi A, Abbasi S (2021) Experimental and numerical study of the effects of geometric appendance elements on energy dissipation over stepped spillway. Water (Switzerland) 13(7). https://doi.org/10.3390/w13070957
Ghaderi A, Dasineh M, Aristodemo F, Ghahramanzadeh A (2020) Characteristics of free and submerged hydraulic jumps over different macroroughnesses. J Hydroinform 22(6):1554–1572. https://doi.org/10.2166/HYDRO.2020.298ArticleGoogle Scholar
Nangare PB, Kote AS (2017) Experimental investigation of an ogee stepped spillway with plain and slotted roller bucket for energy dissipation. Int J Civ Eng Technol 8(8):1549–1555Google Scholar
Parsaie A, Moradinejad A, Haghiabi AH (2018) Numerical modeling of flow pattern in spillway approach channel. Jordan J Civ Eng 12(1):1–9Google Scholar
Pasbani Khiavi M, Ali Ghorbani M, Yusefi M (2021) Numerical investigation of the energy dissipation process in stepped spillways using finite volume method. J Irrig Water Eng 11(4):22–37Google Scholar
Peng Y, Zhang X, Yuan H, Li X, Xie C, Yang S, Bai Z (2019) Energy dissipation in stepped spillways with different horizontal face angles. Energies 12(23). https://doi.org/10.3390/en12234469
Raza A, Wan W, Mehmood K (2021) Stepped spillway slope effect on air entrainment and inception point location. Water (Switzerland) 13(10). https://doi.org/10.3390/w13101428
Rong Y, Zhang T, Peng L, Feng P (2019) Three-dimensional numerical simulation of dam discharge and flood routing in Wudu reservoir. Water (Switzerland) 11(10). https://doi.org/10.3390/w11102157
Saqib N, Akbar M, Pan H, Ou G, Mohsin M, Ali A, Amin A (2022) Numerical analysis of pressure profiles and energy dissipation across stepped spillways having curved risers. Appl Sci 12(448):1–18Google Scholar
Saqib N, Ansari K, Babar M (2021) Analysis of pressure profiles and energy dissipation across stepped spillways having curved treads using computational fluid dynamics. Intl Conf Adv Mech Eng :1–10
Saqib Nu, Akbar M, Huali P, Guoqiang O (2022) Numerical investigation of pressure profiles and energy dissipation across the stepped spillway having curved treads using FLOW 3D. Arab J Geosci 15(1):1363–1400. https://doi.org/10.1007/s12517-022-10505-8ArticleGoogle Scholar
Sarkardeh H, Marosi M, Roshan R (2015) Stepped spillway optimization through numerical and physical modeling. Int J Energy Environ 6(6):597–606Google Scholar
Valero D, Bung DB, Crookston BM, Matos J (2016) Numerical investigation of USBR type III stilling basin performance downstream of smooth and stepped spillways. 6th International Symposium on Hydraulic Structures: Hydraulic Structures and Water System Management, ISHS 2016, 3406281608, 635–646. https://doi.org/10.15142/T340628160853
슈트 여수로의 흐름 폭기는 캐비테이션 손상을 방지하는 가장 효과적이고 경제적인 방법 중 하나입니다. 수중 프리즘에 아주 작은 양의 공기가 흩어지면 표면 손상이 크게 줄어듭니다. 이를 위해 폭기 장치로 알려진 구조를 사용할 수 있습니다. 또한, 램프 각도는 폭기 효율에 영향을 미치는 요인 중 하나입니다. 이 연구에서는 Flow-3D 소프트웨어를 사용하여 3가지 다른 시나리오인 6, 8 및 10도의 램프 각도에서 Jarreh 댐의 방수로를 통해 흐름을 동반하는 공기의 값을 시뮬레이션했습니다. 6도의 경사각에서 유동 유체로 유입되는 공기의 결과를 검증하기 위해이란 TAMAB Company의 실험실에서 댐 방수로 물리적 모델의 관찰 결과를 사용했습니다. 결과에 따르면 램프 각도를 높이면 워터제트 기저귀로 유입되는 공기가 증가하고 10도 램프 각도는 최고의 폭기 효율을 제공합니다. Flow-3D 모델은 결과에 따라 여수로의 2단계 물-공기 흐름을 시뮬레이션할 수도 있습니다.
Flow aeration in chute spillway is one of the most effective and economic ways to prevent cavitation damage. Surface damage is significantly reduced when very small values of air are scattered in a water prism. A structure known as an aerator may be used for this purpose. Besides, ramp angle is one of the factors influencing aerator efficiency. In this research, the value of air entraining the flow through the Jarreh Dam’s spillway at the ramp angles of 6, 8 and 10 degrees, as three different scenarios, was simulated using the Flow-3D software. In order to validate the results of the inlet air into the flowing fluid at a ramp angle of 6 degrees, the observational results of the dam spillway physical model from the laboratory of TAMAB Company in Iran were used. According to the results, raising the ramp angle increases the inlet air to the water jet nappe, and a ten-degree ramp angle provides the best aeration efficiency. The Flow-3D model can also simulate the two-phase water-air flow on spillways, according to the results.
Fig. 1- Schematic of the general pattern of flow and aeration process in the aeratorsFig. 2- Experimental setup (Shamloo et al., 2012)Fig. 3- Results of numerical model validation in determining a) mean flow depth, b) mean velocity,
and c) static pressure in various discharges vs (Shamloo et al., 2012) research under a 6 degree ramp
angleFig. 4- Location of data extraction stations after aeration on a scale model of 1:50Fig.7- Changes in cavitation index in different discharges with changes in ramp angle: a) 6 degrees,
b) 8 degrees and c) 10 degrees
Baharvand, S., & Lashkar-Ara, B. (2021). 실험 모델과 CFD 모델을 결합한 수정 사행 C형 어로의 수력학적 설계기준. 생태 공학 , 164 . https://doi.org/10.1016/j.ecoleng.2021.106207
2- Bayon, A., Toro, JP, Bombardelli, FA, Matos, J., & López-Jiménez, PA(2018). VOF 기술, 난류 모델 및 이산화 방식이 계단식 배수로에서 폭기되지 않은 스키밍 흐름의 수치 시뮬레이션에 미치는 영향. 수력 환경 연구 저널 , 19 , 137–149. https://doi.org/10.1016/j.jher.2017.10.002
3- Brethour, JM, & Hirt, CW (2009). 2성분 흐름에 대한 드리프트 모델. Flow Science, Inc. , FSI – 09 – TN83Rev , 1–7.
4- Chanson, H. (1989). 공기 유입 및 폭기 장치 연구. 수력학 연구 저널 , 27 (3), 301–319. https://doi.org/10.1080/00221688909499166
5- Dong, Z., Wang, J., Vetsch, DF, Boes, RM, & Tan, G. (2019). 매우 높은 단위 배출에서 X자형 플레어링 게이트 교각 뒤의 계단식 배수로에서 공기-물 2상 흐름의 수치 시뮬레이션. 물(스위스) , 11 (10). https://doi.org/10.3390/w11101956
6- Flow-3D, V. 11. 2. (2017). 사용자 매뉴얼. Flow Science Inc.: Santa Fe, NM, USA;
7- Hirt, CW (2003). 자유 표면에서 공기의 난류 동반 모델링. Flow Science, Inc. , FSI – 03 – TN6 , 1–9.
8- Hirt, CW (2016). 드리프트 플럭스에 대한 동적 액적 크기. Flow Science, Inc. , 1–10.
9- Hirt, CW, & Nichols, BD (1981). 자유 경계의 역학에 대한 VOF(유체 체적) 방법. 전산 물리학 저널 , 39 (1), 201–225. https://doi.org/10.1016/0021-9991(81)90145-5
10- Kherbache, K., Chesneau, X., Zeghmati, B., Abide, S., & Benmamar, S. (2017). 계단식 배수로의 물 흐름에 대한 계단식 경사 및 공기 주입의 영향: 수치 연구. 유체 역학 저널 , 29 (2), 322–331. https://doi.org/10.1016/S1001-6058(16)60742-4
11- Kramer, M., & Chanson, H. (2019). 폭기된 여수로 흐름에서 광학 흐름 추정: 샘플링 매개변수에 대한 필터링 및 논의. 실험적 열 및 유체 과학 , 103 , 318–328. https://doi.org/10.1016/j.expthermflusci.2018.12.002
12- Mahmoudian, Z., Baharvand, S., & Lashkarara, B. (2019). Baffle Fishway Denil Type의 흐름 패턴 조사. 관개 과학 및 공학(JISE) , 42 (3), 179–196.
13- Meireles, IC, Bombardelli, FA 및 Matos, J. (2014). 가파른 계단식 배수로의 스키밍 흐름에서 공기 유입 시작: 분석. 수력학 연구 저널 , 52 (3). https://doi.org/10.1080/00221686.2013.878401
14- Parsaie, A., & Haghiabi, AH (2019). 1/4 원형 볏이 있는 계단식 배수로에서 흐름 폭기의 시작 지점. 유량 측정 및 계측 , 69 . https://doi.org/10.1016/j.flowmeasinst.2019.101618
15- Richardson, JF, & Zaki W N. (1979). 침전 및 유동화. 파트 1. 트랜스. Inst. 화학 영어 , 32 , 35–53.
16- Shamloo, H., Hoseini Ghafari, S., & Kavianpour, M. (2012). 슈트 여수로의 폭기에 대한 유입구 흐름의 영향에 대한 실험적 연구(사례 연구: 이란 Jare Dam). 제10차 토목 공학 발전에 관한 국제 회의, 중동 기술 대학, 앙카라, 터키 .
18- Wei, W., Deng, J., & Zhang, F. (2016). 초임계 슈트 흐름에 대한 자체 폭기 공정 개발. 다상 흐름의 국제 저널 , 79 , 172–180. https://doi.org/10.1016/j.ijmultiphaseflow.2015.11.003
20- Xu, Y., Wang, W., Yong, H., & Zhao, W. (2012). 슈트 폭기 장치에서 제트 흐름의 공동 역류에 대한 조사. 프로시디아 엔지니어링 , 31 , 51–56. https://doi.org/10.1016/j.proeng.2012.01.989
21- Yakhot, V., & Orszag, SA (1986). 난류의 재정규화 그룹 분석. I. 기본 이론. 과학 컴퓨팅 저널 , 1 (1), 3–51. https://doi.org/10.1007/BF01061452
22- Yang, J., Teng, P., & Lin, C. (2019). 넓은 여수로 폭기장치의 통풍구 배치 및 물-기류 거동. Theoretical and Applied Mechanics Letters , 9 (2), 130–143. https://doi.org/10.1016/j.taml.2019.02.009
23- Zhang, G., & Chanson, H. (2016). 자유 표면 폭기와 계단식 슈트의 총 압력 사이의 상호 작용. 실험적 열 및 유체 과학 , 74 , 368–381. https://doi.org/10.1016/j.expthermflusci.2015.12.011
많은 계단식 배수로 지오메트리 설계 지침이 평평한 단계를 위해 개발되었지만 통합 단계를 설계하는 것이 더 효율적으로 작동하는 배수로에 대한 적절한 대안이 될 수 있습니다.
이 논문은 POOL의 다른 높이에서 공기 연행과 보이드 비율의 시작점을 다루는 것을 목표로 합니다. 그 후, FLOW-3D 소프트웨어를 사용하여 POOL과 경사면의 높이를 다르게 하여 폭기된 지역과 폭기되지 않은 지역에서 압력 분포를 평가했습니다.
얻어진 수치 결과와 실험 결과의 비교는 본 연구에 사용된 모든 방류에 대해 잘 일치했습니다. POOL 높이는 시작 지점 위치에 미미한 영향을 미쳤습니다. 공극률의 값은 높은 방류에 비해 낮은 방전에서 더 많은 영향을 받았습니다.
여수로의 마루(통기되지 않은 지역)에서는 음압이 나타나지 않았으며 각 방류에서 마루를 따라 높이가 15cm인 수영장에서 최대 압력 값이 얻어졌습니다.
모든 사면에서 웅덩이 및 평평한 계단형 여수로의 계단층 부근에서는 음압이 형성되지 않았습니다. 그러나 평단식 여수로에 비해 평단식 여수로의 수직면 부근에서 음압이 더 많이 형성되어 평단식 슈트에서 캐비테이션 현상이 발생할 확률이 증가하였습니다.
Study of inception point, void fraction and pressure over pooled stWhile many stepped spillways geometry design guidelines were developed for flat steps, designing pooled steps might be an appropriate alternative to spillways working more efficiency. This paper aims to deal with the inception point of air-entrainment and void fraction in the different height of the pools. Following that, pressure distribution was evaluated in aerated and non-aerated regions under the effect of different heights of the pools and slopes through the use of the FLOW-3D software. Comparison of obtained numerical results with experimental ones was in good agreement for all discharges used in this study. Pools height had the insignificant effect on the inception point location. The value of void fraction was more affected in lower discharges in comparison with higher ones. Negative pressure was not seen over the crest of spillway (non-aerated region), and the maximum pressure values were obtained for pools with 15 cm height along the crest in each discharge. In all slopes, negative pressure was not formed near the step bed in the pooled and flat stepped spillways. However, negative pressure was formed in more area near the vertical face in the flat stepped spillway compared with the pooled stepped spillway which increases the probability of cavitation phenomenon in the flat stepped chute.
Design/methodology/approach
압력, 공극률 및 시작점을 평가하기 위해 POOL된 계단식 여수로가 사용되었습니다. 또한 POOL의 다른 높이가 사용되었습니다. 이 연구의 수치 시뮬레이션은 Flow-3D 소프트웨어를 통해 수행되었습니다. 얻어진 결과는 풀이 압력, 공극률 및 시작점을 포함한 2상 유동 특성에 영향을 미칠 수 있음을 나타냅니다.
Findings
마루 위에는 음압이 보이지 않았습니다. 압력 값은 사용된 모든 높이와 15cm 높이에서 얻은 최대 값에 대해 다릅니다. 또한, 풀링 스텝은 플랫 케이스에 비해 음압점 감소에 더 효과적인 역할을 하였습니다. 시작 지점 위치는 특히 9 및 15cm 높이에 대해 스키밍 흐름 영역과 비교하여 낮잠 및 전환 흐름 영역에서 더 많은 영향을 받았습니다.
Figure 1- Schematic diagram of pooled stepped spillway conducted by Felder et al. (2012A): Notes: h
step height (10 cm): w pool height (3.1 cm): l horizontal step length (20 cm): lw pool weir length (1.5 cm):
d’ is the water depth above the crest; y’ is the distance normal to the crest invertFigure 2- meshing domain and distribution of blocksFigure 3- Comparison of numerical simulation with experimental data by Felder et al. (2012A);
mesh convergence analysis; pooled stepped spillway (slope: 26.6 0
)Figure 4- Comparison of numerical simulation with experimental data by Felder et al. (2012A);
Flat stepped spillway (slope: 0
26 6. )Figure 5-Comparison of numerical simulation with experimental data by Felder et al. (2012B); pooled
and flat stepped spillways (slope: 0
9.8 )Figure 6- TKE distribution on steps 8, 9 and 10 for four different mesh numbers: 261252 (model 1),
288941 (model 2), 323578 (model 3) and 343154 (model 4)Figure 7- Comparison of obtained Void fraction distribution on step 10 in numerical simulation with
experimental work conducted by Felder et al. (2012A); (slope 26.60
)Figure 8- Results of inception point of air entrainment in different height of the pools: comparison with
empirical correlations (Eqs 8-9), experimental (Felder et al. (2012A)) and numerical dataFigure 9- Void fraction distribution for different pool heights on steps 10; slope 26.6 0Figure 10- Comparison of pressure distribution between numerical simulation and experimental work
conducted by Zhang and Chanson (2016); flat stepped spillway (slope: 0
45 )Figure 11- A comparison of the pressure distribution above the crest of the spillway; B comparison of the
free surface profile along the crest of the spillway.
Note: x’ indicates the longitudinal distance from the starting point of the crest.Figure 12- pressure distribution along crest of spillway in different discharges; slope 26.6Figure 13- Pressure distribution near the last step bed for different slopes and discharges: x” indicatesthe
longitudinal distance from the intersection of the horizontal and vertical faces of step 10; y” is the
distance from the intersection of the horizontal and vertical faces in the vertical directionFigure 14- Pressure distribution adjacent the vertical face of step 9 for different discharges and slopesTable1- Used discharges for assessments of mesh convergence analysis and hydraulic
characteristics
Conclusion
본 연구에서는 자유표면을 모사하기 위해 VOF 방법과 k -ε (RNG) 난류 모델을 활용하여 FLOW-3D 소프트웨어를 사용하였고, 계단식 배수로의 유동을 모사하기 위한 목적으로 난류 특성을 모사하였다. 얻은 결과는 수치 모델이 시작점 위치, 보이드 비율 및 압력을 적절하게 시뮬레이션했음을 나타냅니다. 풀의 높이는 공기 유입 위치에 미미한 영향을 미치므로 얻은 결과는 이 문서에서 제시된 상관 관계와 잘 일치했습니다. 즉, 사용 가능한 상관 관계를 서로 다른 풀 높이에 사용할 수 있습니다. 공극률의 결과는 스텝 풀 근처의 나프 유동 영역에서 공극율 값이 다른 배출보다 더 큰 것으로 나타났다. 더욱이 고방출량 .0 113m3/s에서 수영장 높이를 변경해도 수영장 표면 근처의 공극률 값에는 영향을 미치지 않았습니다.
낮잠 및 전환 체제의 압력 분포에 대한 0 및 3cm 높이의 수영장 효과는 많은 지점에서 대부분 유사했습니다. 더욱이 조사된 모든 높이에서 여수로의 마루를 따라 부압이 없었습니다. 여수로 끝단의 바닥 부근의 압력 결과는 평평하고 고인 경우 부압이 발생하지 않았음을 나타냅니다. 수직면 부근의 음압은 웅덩이에 비해 평평한 계단형 여수로의 깊이(w=0 cm)의 대부분에서 발생하였다. 또한 더 큰 사면에 대한 풀링 케이스에서 음압이 제거되었습니다. 평단식 여수로에서는 계단의 수직면에 인접한 더 넓은 지역에서 음압이 발생하였기 때문에 이 여수로에서는 고형단식여수로보다 캐비테이션 현상이 발생할 가능성이 더 큽니다.
In this study, the FLOW-3D software was used through utilizing the VOF method and k −ε (RNG) turbulence model in order to simulate free surface, and turbulence characteristics for the purpose of simulating flow over pooled stepped spillway. The results obtained indicated that the numerical model properly simulated the inception point location, void fraction, and pressure. The height of the pools has the insignificant effect on the location of air entrainment, so that obtained results were in good agreement with the correlations presented in this paper. In other words, available correlations can be used for different pool heights. The results of void fraction showed that the void fraction values in nappe flow regime near the step pool were more than the other discharges. Furthermore in high discharge, 0.113m3/s, altering pool height had no effect on the value of void fraction near the pool surface.
The effect of the pools with 0 and 3 cm heights over the pressure distribution in nappe and transition regimes was mostly similar in many points. Furthermore, in all examined heights there was no negative pressure along the crest of the spillway. The pressure results near the bed of the step at the end of the spillway indicated that negative pressure did not occur in the flat and pooled cases. Negative pressure near the vertical face occurred in the most part of the depth in the flat stepped spillway (w=0 cm) in comparison with the pooled case. Also, the negative pressure was eliminated in the pooled case for the larger slopes. Since negative pressure occurred in a larger area adjacent the vertical face of the steps in the flat stepped spillways, it is more likely that cavitation phenomenon occurs in this spillway rather than the pooled stepped spillways.
References
André, S. (2004), “High velocity aerated flows on stepped chutes with macro-roughness elements.” Ph.D. thesis, Laboratoire de Constructions Hydraulics (LCH), EPFL, Lausanne, Switzerland, 272 pages.
Attarian, A. Hosseini, Kh. Abdi, H and Hosseini, M. (2014), “The Effect of the Step Height on Energy Dissipation in Stepped Spillways Using Numerical Simulation”. Arabian Journal for Science and Engineering, 39(4), 2587-2594.
Bombardelli, F.A. Meireles. I. Matos, J. (2011), “Laboratory measurements and multi-block numerical simulations of the mean flow and turbulence in the non-aerated skimming flow region of steep stepped spillways”. Environmental fluid mechanics, 11(3) 263-288.
Chakib, B. (2013), “Numerical Computation of Inception Point Location for Flat-sloped Stepped Spillway”. International Journal of Hydraulic Engineering; 2(3): 47-52.
Chakib, B. Mohammed, H. (2015), “Numerical Simulation of Air Entrainment for Flat-Sloped Stepped Spillway. Journal of computational multiphase flows”, Volume 7. Number 1.
Chanson, H. Toombes, L. (2002), “Air–water flows down stepped chutes: turbulence and flow structure observations”. International Journal of Multiphase Flow, 28(11) 1737-1761
Chen, Q. Dai, G. Liu, H. (2002), “Volume of Fluid Model for Turbulence Numerical Simulation of Stepped Spillway Overflow”. DOI: 10.1061/(ASCE)0733-9429128:7(683).
Cheng, X. Chen, Y. Luo, L. (2006), “Numerical simulation of air-water two-phase flow over stepped spillways”. Science in China Series E: Technological Sciences, 49(6), 674-684.
Cheng, X. Luo, L. Zhao, W. (2004), “Study of aeration in the water flow over stepped spillway”. In: Proceedings of the world water congress.
Chinnarasri, Ch. Kositgittiwong, D. Julien, Y. (2013), “Model of flow over spillways by computational fluid dynamics”. Proceedings of the ICE – Water Management, Volume 167(3) 164 –175.
Dastgheib, A. Niksokhan, M.H. and Nowroozpour, A.R. (2012), “Comparing of Flow Pattern and Energy Dissipation over different forms of Stepped Spillway”. World Environmental and Water Resources Congress ASCE.
Eghbalzadeh, A. Javan, M. (2012), “Comparison of mixture and VOF models for numerical simulation of air entrainment in skimming flow over stepped spillway”. Procedia Engineering, 28. 657-660.
Felder, S, Chanson, H. (2012), “Free-surface Profiles, Velocity and Pressure Distributions on a Broad-Crested Weir: a Physical study “Free-surface Profiles, Velocity and Pressure Distributions on a Broad-Crested Weir: a Physical study
Felder, S. Fromm, Ch. Chanson, H. (2012B), “Air entrainment and energy dissipation on a 8.9 slope stepped spillway with flat and pooled steps”, School of Civil Engineering, The University of Queensland,. Brisbane, Australia.
Felder, S. Chanson, H. (2014A), Triple decomposition technique in air–water flows: application to instationary flows on a stepped spillway. International Journal of Multiphase Flow, 58, 139-153.
Felder, S. Chanson, H. (2014B), Effects of step pool porosity upon flow aeration and energy dissipation on pooled stepped spillways. Journal of Hydraulic Engineering, 140(4), 04014002.
Felder, S. Chanson, H. (2013A), “Air entrainment and energy dissipation on porous pooled stepped spillways”. Paper presented at the International Workshop on Hydraulic Design of Low-Head Structures.
Felder, S. Chanson, H. (2013B), “Aeration, flow instabilities, and residual energy on pooled stepped spillways of embankment dams”. Journal of irrigation and drainage engineering, 139(10) 880-887.
Felder, S. Guenther, Ph. Chanson, H. (2012A). “Air-water flow properties and energy dissipation on stepped spillways: a physical study of several pooled stepped configurations”, School of Civil Engineering, The University of Queensland,. Brisbane, Australia.
Flow Science, (2013). “FLOW-3D user’s manual”, version 10.1. Flow Science, Inc, Los Alamos.
Frizell, K.W. Renna, F.M. Matos, J. (2012), “Cavitation potential of flow on stepped spillways”. Journal of Hydraulic Engineering, 139(6), 630-636.
Gonzalez, C. (2005), “An experimental study of free-surface aeration on embankment stepped chutes”, department of civil engineering, Brisbane, Australia, Phd thesis.
Gonzalez, C.A. Chanson, H. (2008), “Turbulence manipulation in air–water flows on a stepped chute: An experimental study”. European Journal of Mechanics-B/Fluids, 27(4), 388-408.
Guenther, Ph.. Felder, S. Chanson, H. (2013), “Flow aeration, cavity processes and energy dissipation on flat and pooled stepped spillways for embankments”. Environmental fluid mechanics, 13(5) 503-525.
Hamedi, A. Mansoori, A. Malekmohamadi, I. Roshanaei, H. (2011), “Estimating Energy Dissipation in Stepped Spillways with Reverse Inclined Steps and End Sill”. World Environmental and Water Resources Congress, ASCE.
Hirt, C.W. (2003), “Modeling Turbulent Entrainment of Air at a Free Surface”. Flow Science Inc.
Hunt, S.L. Kadavy, K.C. (2013), “Inception point for enbankment dam stepped spillway”. J. Hydraul. Eng., 139(1), 60–64.
Hunt, S.L. Kadavy, K.C. (2010), “Inception Point Relationship for Flat-Sloped Stepped Spillways”. DOI: 10.1061/ASCEHY.1943-7900.0000297.
Matos, J. Quintela, A. (2000), “Air entrainment and safety against cavitation damage in stepped spillways over RCC dams. In: Proceeding Intl. Workshop on Hydraulics of Stepped Spillways”, VAW, ETH-Zurich, H.E. Minor and W.H. Hager. Balkema. 69–76.
Meireles, I. Matos, J. (2009), “Skimming flow in the nonaerated region of stepped spillways over embankment dams”. J. Hydraul. Eng., 135(8), 685–689.
Miang-liang, ZH. Yong-ming, SH. (2008), “Three dimentional simulation of meandering river basin on 3-D RNG k − ε turbulence model”. Journal of hydrodynamics, 20(4): 448-455.
Morovati, Kh. Eghbalzadeh, A. Javan, M. (2015), “Numerical investigation of the configuration of the pools on the flowPattern passing over pooled stepped spillway in skimming flow regime. Acta Mech, DOI 10.1007/s00707-015-1444-x
Morovati, Kh. Eghbalzadeh, A. Soori, S. (2016), “Numerical Study of Energy Dissipation of Pooled Stepped spillway”. Civil Engineering Journal. Vol. 2, No. 5.
Nikseresht, A.H. Talebbeydokhti, N. and Rezaei, M.J. (2013), “Numerical simulation of two-phase flow on steppool spillways”. Scientia Iranica, A 20 (2), 222–230.
Peyras, L. Royet, P. Degoutte, G. (1990), “Flow and energy dissipation over stepped gabion weirs”. ASCE Convention.
Qun, Ch. Guang-qing, D. Feu-qing, Zh. Qing, Y. (2004). “Three-dimensional turbulence numerical simulation of a stepped spillway overflow”. Journal of hydrodynamics, Ser. B, 1, 74-79.
Relvas, A. T. Pinheiro, A. N. (2008), Inception point and air concentration in flows on stepped chutes lined with wedge-shaped concrete blocks. Journal of Hydraulic Engineering, 134(8), 1042-1051
Sanchez, M. (2000), “Pressure field in skimming flow over a stepped spillways”. In: Proceeding Intl. Workshop on Hydraulics of Stepped Spillways, VAW, ETH-Zurich, H.E. Minor and W.H. Hager. Balkema, 137–146.
Sarfaraz, M. Attari, J. Pfister, M. (2012), “Numerical Computation of Inception Point Location for Steeply Sloping Stepped Spillways”. 9th International Congress on Civil Engineering, May 8-10. Isfahan University of Technology (IUT), Isfahan, Iran.
Savage, Bruce M. Michael C. Johnson. (2001), “Flow over ogee spillway: Physical and numerical model case study.” Journal of Hydraulic Engineering 127.8:640-649.
Shahhedari, H. Jafari Nodoshan, E. Barati, R. Azhdary moghadam, M. (2014). “Discharge coeficient and energy dissipation over stepped spillway under skimming flow regime”. KSCE Journal of Civil Engineering, DOI 10.1007/s12205-013-0749-3.
Tabbara, M. Chatila, J. Awwad, R. (2005), “Computational simulation of flow over stepped spillways”. Computers & structures, 83(27) 2215-2224.
Thorwarth, J. (2008), “Hydraulisches Verhalten der Treppengerinne mit eingetieften Stufen—Selbstinduzierte Abflussinstationaritäten und Energiedissipation” [Hydraulics of pooled stepped spillways— Self-induced unsteady flow and energy dissipation]. Ph.D. thesis, Univ. of Aachen, Aachen, Germany (in German).
WeiLin, XU. ShuJing, LUO, QiuWen, ZH. Jing, LUO. (2015), “Experimental study on pressure and aeration characteristics in stepped chute flows. SCIENCE CHINA. Vol.58 No.4: 720–726. doi: 10.1007/s11431-015- 5783-6.
Xiangju, Ch. Yongcan, C. Lin, L. (2006), “Numerical simulation of air-water two-phase flow over stepped spillways”. Science in China Series E: Technological Sciences, 49(6), 674-684.
Zare, K.H. Doering, J.C. (2012), “Inception Point of Air Entrainment and Training Wall Characteristics of Baffles and Sills on Stepped Spillways”. DOI: 10.1061/(ASCE)HY .1943-7900.0000630.
Zhan, J. Zhang, J. Gong, Y. (2016), “Numerical investigation of air-entrainment in skimming flow over stepped spillways”. Theoretical and Applied Mechanics Letters. Volume 6. Pages 139–142.
Zhang, G. Chanson, H. (2016), Hydraulics of the developing flow region of stepped spillways. II: Pressure and velocity fields. Journal of Hydraulic Engineering, 142(7).
Zhenwei, M. Zhiyan, Zh. Tao, Zh. (2012), “Numerical Simulation of 3-D Flow Field of Spillway based on VOF Method”. Procedia Engineering, 28, 808-812.
Zhi-yong, D. Hun-wei, L.J. (2006), “Numerical simulation of skimming flow over mild stepped channel”. Journal of Hydrodynamics, Ser. B, 18(3) 367-371.
ZhongDong, Q. XiaoQing, H. WenXin, H. António, A. (2009), “Numerical simulation and analysis of water flow over stepped spillways”. Science in China Series E: Technological Sciences, 52(7) 1958-1965.
Probabilistic investigation of cavitation occurrence in chute spillway based on the results of Flow-3D numerical modeling
Amin Hasanalipour Shahrabadi1*, Mehdi Azhdary Moghaddam2
1-University of Sistan and Baluchestan،amin.h.shahrabadi@gmail.com
2-University of Sistan and Baluchestan،Mazhdary@eng.usb.ac.ir
Abstract
Probabilistic designation is a powerful tool in hydraulic engineering. The uncertainty caused by random phenomenon in hydraulic design may be important. Uncertainty can be expressed in terms of probability density function, confidence interval, or statistical torques such as standard deviation or coefficient of variation of random parameters. Controlling cavitation occurrence is one of the most important factors in chute spillways designing due to the flow’s high velocity and the negative pressure (Azhdary Moghaddam & Hasanalipour Shahrabadi, ۲۰۲۰). By increasing dam’s height, overflow velocity increases on the weir and threats the structure and it may cause structural failure due to cavitation (Chanson, ۲۰۱۳). Cavitation occurs when the fluid pressure reaches its vapor pressure. Since high velocity and low pressure can cause cavitation, aeration has been recognized as one of the best ways to deal with cavitation (Pettersson, ۲۰۱۲). This study, considering the extracted results from the Flow-۳D numerical model of the chute spillway of Darian dam, investigates the probability of cavitation occurrence and examines its reliability. Hydraulic uncertainty in the design of this hydraulic structure can be attributed to the uncertainty of the hydraulic performance analysis. Therefore, knowing about the uncertainty characteristics of hydraulic engineering systems for assessing their reliability seems necessary (Yen et al., ۱۹۹۳). Hence, designation and operation of hydraulic engineering systems are always subject to uncertainties and probable failures. The reliability, ps, of a hydraulic engineering system is defined as the probability of safety in which the resistance, R, of the system exceeds the load, L, as follows (Chen, ۲۰۱۵): p_s=P(L≤R) (۱) Where P(۰) is probability. The failure probability, p_f, is a reliability complement and is expressed as follows: p_f=P[(L>R)]=۱- p_s (۲) Reliability development based on analytical methods of engineering applications has come in many references (Tung & Mays, ۱۹۸۰ and Yen & Tung, ۱۹۹۳). Therefore, based on reliability, in a control method, the probability of cavitation occurrence in the chute spillway can be investigated. In reliability analysis, the probabilistic calculations must be expressed in terms of a limited conditional function, W(X)=W(X_L ,X_R)as follows: p_s=P[W(X_L ,X_R)≥۰]= P[W(X)≥۰] (۳) Where X is the vector of basic random variables in load and resistance functions. In the reliability analysis, if W(X)> ۰, the system will be secure and in the W(X) <۰ system will fail. Accordingly, the eliability index, β, is used, which is defined as the ratio of the mean value, μ_W, to standard deviation, σ_W, the limited conditional function W(X) is defined as follows (Cornell, ۱۹۶۹): β=μ_W/σ_W (۴) The present study was carried out using the obtained results from the model developed by ۱:۵۰ scale plexiglass at the Water Research Institute of Iran. In this laboratory model, which consists of an inlet channel and a convergent thrower chute spillway, two aerators in the form of deflector were used at the intervals of ۲۱۱ and ۲۷۰ at the beginning of chute, in order to cope with cavitation phenomenon during the chute. An air duct was also used for air inlet on the left and right walls of the spillway. To measure the effective parameters in cavitation, seven discharges have been passed through spillway. As the pressure and average velocity are determined, the values of the cavitation index are calculated and compared with the values of the critical cavitation index, σ_cr. At any point when σ≤σ_cr, there is a danger of corrosion in that range (Chanson, ۱۹۹۳). In order to obtain uncertainty and calculate the reliability index of cavitation occurrence during a chute, it is needed to extract the limited conditional function. Therefore, for a constant flow between two points of flow, there would be the Bernoulli (energy) relation as follows (Falvey, ۱۹۹۰): σ= ( P_atm/γ- P_V/γ+h cosθ )/(〖V_۰〗^۲/۲g) (۵) Where P_atm is the atmospheric pressure, γ is the unit weight of the water volume, θ is the angle of the ramp to the horizon, r is the curvature radius of the vertical arc, and h cosθ is the flow depth perpendicular to the floor. Therefore, the limited conditional function can be written as follows: W(X)=(P_atm/γ- P_V/γ+h cosθ )/(〖V_۰〗^۲/۲g) -σ_cr (۶) Flow-۳D is a powerful software in fluid dynamics. One of the major capabilities of this software is to model free-surface flows using finite volume method for hydraulic analysis. The spillway was modeled in three modes, without using aerator, ramp aerator, and ramp combination with aeration duct as detailed in Flow-۳D software. For each of the mentioned modes, seven discharges were tested. According to Equation (۶), velocity and pressure play a decisive and important role in the cavitation occurrence phenomenon. Therefore, the reliability should be evaluated with FORM (First Order Reliable Method) based on the probability distribution functions For this purpose, the most suitable probability distribution function of random variables of velocity and pressure on a laboratory model was extracted in different sections using Easy fit software. Probability distribution function is also considered normal for the other variables in the limited conditional function. These values are estimated for the constant gravity at altitudes of ۵۰۰ to ۷۰۰۰ m above the sea level for the unit weight, and vapor pressure at ۵ to ۳۵° C. For the critical cavitation index variable, the standard deviation is considered as ۰.۰۱. According to the conducted tests, for the velocity random variable, GEV (Generalized Extreme Value) distribution function, and for the pressure random variable, Burr (۴P) distribution function were presented as the best distribution function. The important point is to not follow the normal distribution above the random variables. Therefore, in order to evaluate the reliability with the FORM method, according to the above distributions, they should be converted into normal variables based on the existing methods. To this end, the non-normal distributions are transformed into the normal distribution by the method of Rackwitz and Fiiessler so that the value of the cumulative distribution function is equivalent to the original abnormal distribution at the design point of x_(i*). This point has the least distance from the origin in the standardized space of the boundary plane or the same limited conditional function. The reliability index will be equal to ۰.۴۲۰۴ before installing the aerator. As a result, reliability, p_s, and failure probability, p_f, are ۰.۶۶۲۹ and ۰.۳۳۷۱, respectively. This number indicates a high percentage for cavitation occurrence. Therefore, the use of aerator is inevitable to prevent imminent damage from cavitation. To deal with cavitation as planned in the laboratory, two aerators with listed specifications are embedded in a location where the cavitation index is critical. In order to analyze the reliability of cavitation occurrence after the aerator installation, the steps of the Hasofer-Lind algorithm are repeated. The modeling of ramps was performed separately in Flow-۳D software in order to compare the performance of aeration ducts as well as the probability of failure between aeration by ramp and the combination of ramps and aeration ducts. Installing an aerator in combination with a ramp and aerator duct greatly reduces the probability of cavitation occurrence. By installing aerator, the probability of cavitation occurrence will decrease in to about ۴ %. However, in the case of aeration only through the ramp, the risk of failure is equal to ۱۰%.
확률적 지정은 수력 공학에서 강력한 도구입니다. 유압 설계에서 임의 현상으로 인한 불확실성이 중요할 수 있습니다. 불확실성은 확률 밀도 함수, 신뢰 구간 또는 표준 편차 또는 무작위 매개변수의 변동 계수와 같은 통계적 토크로 표현될 수 있습니다. 캐비테이션 발생을 제어하는 것은 흐름의 높은 속도와 음압으로 인해 슈트 여수로 설계에서 가장 중요한 요소 중 하나입니다(Azhdary Moghaddam & Hasanalipour Shahrabadi, ۲۰۲۰). 댐의 높이를 높이면 둑의 범람속도가 증가하여 구조물을 위협하고 캐비테이션으로 인한 구조물의 파손을 유발할 수 있다(Chanson, ۲۰۱۳). 캐비테이션은 유체 압력이 증기압에 도달할 때 발생합니다. 높은 속도와 낮은 압력은 캐비테이션을 유발할 수 있으므로, 통기는 캐비테이션을 처리하는 가장 좋은 방법 중 하나로 인식되어 왔습니다(Pettersson, ۲۰۱۲). 본 연구에서는 Darian 댐의 슈트 여수로의 Flow-۳D 수치모델에서 추출된 결과를 고려하여 캐비테이션 발생 확률을 조사하고 그 신뢰성을 조사하였다. 이 수력구조의 설계에서 수력학적 불확실성은 수력성능 해석의 불확실성에 기인할 수 있다. 따라서 신뢰성을 평가하기 위해서는 수력공학 시스템의 불확도 특성에 대한 지식이 필요해 보인다(Yen et al., ۱۹۹۳). 따라서 수력 공학 시스템의 지정 및 작동은 항상 불확실성과 가능한 고장의 영향을 받습니다. 유압 공학 시스템의 신뢰성 ps는 저항 R, 시스템의 부하 L은 다음과 같이 초과됩니다(Chen, ۲۰۱۵): p_s=P(L≤R)(۱) 여기서 P(۰)은 확률입니다. 고장 확률 p_f는 신뢰도 보완이며 다음과 같이 표현됩니다. Mays, ۱۹۸۰ 및 Yen & Tung, ۱۹۹۳). 따라서 신뢰성을 기반으로 제어 방법에서 슈트 여수로의 캐비테이션 발생 확률을 조사할 수 있습니다. 신뢰도 분석에서 확률적 계산은 제한된 조건부 함수 W(X)=W(X_L , X_R)은 다음과 같습니다. p_s=P[W(X_L,X_R)≥۰]= P[W(X)≥۰] (۳) 여기서 X는 부하 및 저항 함수의 기본 랜덤 변수 벡터입니다. 신뢰도 분석에서 W(X)> ۰이면 시스템은 안전하고 W(X) <۰에서는 시스템이 실패합니다. 따라서 표준편차 σ_W에 대한 평균값 μ_W의 비율로 정의되는 신뢰도 지수 β가 사용되며, 제한된 조건부 함수 W(X)는 다음과 같이 정의됩니다(Cornell, ۱۹۶۹). β= μ_W/σ_W (۴) 본 연구는 이란 물연구소의 ۱:۵۰ scale plexiglass로 개발된 모델로부터 얻은 결과를 이용하여 수행하였다. 이 실험 모델에서, 입구 수로와 수렴형 투수 슈트 여수로로 구성되며 슈트 중 캐비테이션 현상에 대처하기 위해 슈트 초기에 ۲۱۱과 ۲۷۰ 간격으로 편향기 형태의 2개의 에어레이터를 사용하였다. 여수로 좌우 벽의 공기 유입구에도 공기 덕트가 사용되었습니다. 캐비테이션의 효과적인 매개변수를 측정하기 위해 7번의 배출이 방수로를 통과했습니다. 압력과 평균 속도가 결정되면 캐비테이션 지수 값이 계산되고 임계 캐비테이션 지수 σ_cr 값과 비교됩니다. σ≤σ_cr일 때 그 범위에서 부식의 위험이 있다(Chanson, ۱۹۹۳). 슈트 중 캐비테이션 발생의 불확실성을 구하고 신뢰도 지수를 계산하기 위해서는 제한된 조건부 함수를 추출할 필요가 있다. 따라서 두 지점 사이의 일정한 흐름에 대해 다음과 같은 Bernoulli(에너지) 관계가 있습니다(Falvey, ۱۹۹۰). σ= ( P_atm/γ- P_V/γ+h cosθ )/(〖V_۰〗 ^۲/۲g) (۵) 여기서 P_atm은 대기압, γ는 물의 단위 중량, θ는 수평선에 대한 경사로의 각도, r은 수직 호의 곡률 반경, h cos θ는 바닥에 수직인 흐름 깊이입니다. 따라서 제한된 조건부 함수는 다음과 같이 쓸 수 있습니다. W(X)=(P_atm/γ- P_V/γ+h cosθ )/(〖V_۰〗^۲/۲g) -σ_cr (۶) Flow-۳D는 유체 역학의 강력한 소프트웨어. 이 소프트웨어의 주요 기능 중 하나는 수리학적 해석을 위해 유한 체적 방법을 사용하여 자유 표면 흐름을 모델링하는 것입니다. 방수로는 Flow-۳D 소프트웨어에 자세히 설명된 바와 같이 폭기 장치, 램프 폭기 장치 및 폭기 덕트가 있는 램프 조합을 사용하지 않고 세 가지 모드로 모델링되었습니다. 언급된 각 모드에 대해 7개의 방전이 테스트되었습니다. 식 (۶)에 따르면 속도와 압력은 캐비테이션 발생 현상에 결정적이고 중요한 역할을 합니다. 따라서 확률분포함수에 기반한 FORM(First Order Reliable Method)으로 신뢰도를 평가해야 한다 이를 위해 실험실 모델에 대한 속도와 압력의 확률변수 중 가장 적합한 확률분포함수를 Easy fit을 이용하여 구간별로 추출하였다. 소프트웨어. 확률 분포 함수는 제한된 조건부 함수의 다른 변수에 대해서도 정상으로 간주됩니다. 이 값은 단위 중량의 경우 해발 ۵۰۰ ~ ۷۰۰۰ m 고도에서의 일정한 중력과 ۵ ~ ۳۵ ° C에서의 증기압으로 추정됩니다. 임계 캐비테이션 지수 변수의 표준 편차는 ۰.۰۱으로 간주됩니다. . 수행된 시험에 따르면 속도 확률변수는 GEV(Generalized Extreme Value) 분포함수로, 압력변수는 Burr(۴P) 분포함수가 가장 좋은 분포함수로 제시되었다. 중요한 점은 확률 변수 위의 정규 분포를 따르지 않는 것입니다. 따라서 FORM 방법으로 신뢰도를 평가하기 위해서는 위의 분포에 따라 기존 방법을 기반으로 정규 변수로 변환해야 합니다. 이를 위해, 비정규분포를 Rackwitz와 Fiiessler의 방법에 의해 정규분포로 변환하여 누적분포함수의 값이 x_(i*)의 설계점에서 원래의 비정상분포와 같도록 한다. 이 점은 경계면의 표준화된 공간 또는 동일한 제한된 조건부 함수에서 원점으로부터 최소 거리를 갖습니다. 신뢰성 지수는 폭기 장치를 설치하기 전의 ۰.۴۲۰۴과 같습니다. 그 결과 신뢰도 p_s와 고장확률 p_f는 각각 ۰.۶۶۲۹과 ۰.۳۳۷۱이다. 이 숫자는 캐비테이션 발생의 높은 비율을 나타냅니다. 따라서 캐비테이션으로 인한 즉각적인 손상을 방지하기 위해 폭기 장치의 사용이 불가피합니다. 실험실에서 계획한 대로 캐비테이션을 처리하기 위해, 나열된 사양을 가진 두 개의 폭기 장치는 캐비테이션 지수가 중요한 위치에 내장되어 있습니다. 폭기장치 설치 후 캐비테이션 발생의 신뢰성을 분석하기 위해 Hasofer-Lind 알고리즘의 단계를 반복합니다. 경사로의 모델링은 폭기 덕트의 성능과 경사로에 의한 폭기 및 경사로와 폭기 덕트의 조합 사이의 실패 확률을 비교하기 위해 Flow-۳D 소프트웨어에서 별도로 수행되었습니다. 경사로 및 폭기 덕트와 함께 폭기 장치를 설치하면 캐비테이션 발생 가능성이 크게 줄어듭니다. 에어레이터를 설치하면 캐비테이션 발생 확률이 약 ۴%로 감소합니다. 그러나 램프를 통한 폭기의 경우 실패 위험은 ۱۰%와 같습니다. 폭기 설치 후 캐비테이션 발생의 신뢰성을 분석하기 위해 Hasofer-Lind 알고리즘의 단계를 반복합니다. 경사로의 모델링은 폭기 덕트의 성능과 경사로에 의한 폭기 및 경사로와 폭기 덕트의 조합 사이의 실패 확률을 비교하기 위해 Flow-۳D 소프트웨어에서 별도로 수행되었습니다. 경사로 및 폭기 덕트와 함께 폭기 장치를 설치하면 캐비테이션 발생 가능성이 크게 줄어듭니다. 에어레이터를 설치하면 캐비테이션 발생 확률이 약 ۴%로 감소합니다. 그러나 램프를 통한 폭기의 경우 실패 위험은 ۱۰%와 같습니다. 폭기장치 설치 후 캐비테이션 발생의 신뢰성을 분석하기 위해 Hasofer-Lind 알고리즘의 단계를 반복합니다. 경사로의 모델링은 폭기 덕트의 성능과 경사로에 의한 폭기 및 경사로와 폭기 덕트의 조합 사이의 실패 확률을 비교하기 위해 Flow-۳D 소프트웨어에서 별도로 수행되었습니다. 경사로 및 폭기 덕트와 함께 폭기 장치를 설치하면 캐비테이션 발생 가능성이 크게 줄어듭니다. 에어레이터를 설치하면 캐비테이션 발생 확률이 약 ۴%로 감소합니다. 그러나 램프를 통한 폭기의 경우 실패 위험은 ۱۰%와 같습니다. 경사로의 모델링은 폭기 덕트의 성능과 경사로에 의한 폭기 및 경사로와 폭기 덕트의 조합 사이의 실패 확률을 비교하기 위해 Flow-۳D 소프트웨어에서 별도로 수행되었습니다. 경사로 및 폭기 덕트와 함께 폭기 장치를 설치하면 캐비테이션 발생 가능성이 크게 줄어듭니다. 에어레이터를 설치하면 캐비테이션 발생 확률이 약 ۴%로 감소합니다. 그러나 램프를 통한 폭기의 경우 실패 위험은 ۱۰%와 같습니다. 경사로의 모델링은 폭기 덕트의 성능과 경사로에 의한 폭기 및 경사로와 폭기 덕트의 조합 사이의 실패 확률을 비교하기 위해 Flow-۳D 소프트웨어에서 별도로 수행되었습니다. 경사로 및 폭기 덕트와 함께 폭기 장치를 설치하면 캐비테이션 발생 가능성이 크게 줄어듭니다. 에어레이터를 설치하면 캐비테이션 발생 확률이 약 ۴%로 감소합니다. 그러나 램프를 통한 폭기의 경우 실패 위험은 ۱۰%와 같습니다. 에어레이터를 설치하면 캐비테이션 발생 확률이 약 ۴%로 감소합니다. 그러나 램프를 통한 폭기의 경우 실패 위험은 ۱۰%와 같습니다. 에어레이터를 설치하면 캐비테이션 발생 확률이 약 ۴%로 감소합니다. 그러나 램프를 통한 폭기의 경우 실패 위험은 ۱۰%와 같습니다.
Keywords
Aerator Probable Failure Reliability Method FORM Flow ۳D.
계단식 배수로에서 공기 유입 시뮬레이션에 대한 다양한 수치 매개 변수의 영향을 조사합니다.
엔지니어링 과제
쉽게 실험 데이터와 일치하도록 공기 유입 및 도움말 보정의 최초의 시뮬레이션에 수치 매개 변수의 영향을 연구에 사용할 수 있는 자동화된 워크 플로우 생성 1 .
연구할 매개 변수는 메쉬 크기, 난류 모델, 난류 길이 스케일 및 동적 대 고정 난류 길이 스케일입니다. 또한 FLOW-3D (x) 는 마지막 시간 단계에서 유입된 공기 농도의 이미지와 시뮬레이션에서 공기 유입의 진화를 보여주는 애니메이션을 생성합니다.
FLOW-3D (x) 워크 플로우
시뮬레이션은 동반된 공기의 양을 보고하기 위해 3, 4, 5 단계의 샘플링 볼륨으로 설정됩니다. FLOW-3D (x) 는 노드를 사용하여 자동화된 워크 플로를 구성합니다.
첫 번째 노드는 .csv 파일에서 시뮬레이션 매개 변수를 읽는데 사용됩니다. 그런 다음 매개 변수는 시뮬레이션을 실행 하기 위해 FLOW-3D 노드로 전송됩니다 . 후 처리 노드는 배수로의 각 단계에서 샘플링 볼륨에서 동반된 공기 볼륨을 추출하고, 마지막 시간 단계에서 동반된 공기의 이미지를 생성하고, 공기 동반 애니메이션을 생성합니다. 마지막 노드는 샘플링 볼륨에서 보고된 동반 공기 값을 .csv 파일에 씁니다.
매개 변수 정의 입력 파일에 18 개의 매개 변수 세트가 지정되어 있으므로 예산 또는 허용되는 반복 횟수가 18로 설정되었습니다. 단일 시뮬레이션의 런타임은 각 반복에서 사용되는 메시 크기에 따라 다릅니다.
매개 변수 연구 결과
사용 FLOW-3D (X) 의 데이터 분석 기능 및 자동 화상 생성하면 빠른 시각적 평가 결과의 검증을 허용합니다. 또한 각 시뮬레이션 실행에 대한 각 단계의 공기 유입 값은 보고된 .csv 파일에서 쉽게 액세스 할 수 있습니다. 최적화 연구 시간을 절약하기 위해 배치 실행이 사용되었습니다.
교정 전
0.01m의 메쉬 크기, k-ω 난류 모델 및 0.005m와 동일한 난류 길이 척도를 사용하는 시뮬레이션의 마지막 시간 단계에서 유입 공기
References
메쉬 크기 = 0.005m, k-ω 난류 모델 및 0.005m와 같은 난류 길이 척도의 시뮬레이션에서 마지막 시간 단계에서 혼입된 공기. 2 배 더 미세한 메쉬를 사용한 공기 혼입의 시작은 0.01m 메쉬보다 실험 결과와 유사하게 비교됩니다.
참고 문헌
1 Felder, Stefan (2013). Air-water flow properties on stepped spillways for embankment dams: Aeration, energy dissipation and turbulence on uniform, non-uniform and pooled stepped chutes. PhD Thesis, School of Civil Engineering, The University of Queensland.
시뮬레이션 비디오 갤러리에서 FLOW-3D 제품군으로 모델링 할 수 있는 다양한 가능성을 살펴보십시오 .
적층 제조 시뮬레이션 갤러리
FLOW-3D AM 은 레이저 파우더 베드 융합, 바인더 제트 및 직접 에너지 증착과 같은 적층 제조 공정을 시뮬레이션하고 분석합니다. FLOW-3D AM 의 다중 물리 기능은 공정 매개 변수의 분석 및 최적화를 위해 분말 확산 및 압축, 용융 풀 역학, L-PBF 및 DED에 대한 다공성 형성, 바인더 분사 공정을 위한 수지 침투 및 확산에 대한 매우 정확한 시뮬레이션을 제공합니다.
Multi-material Laser Powder Bed Fusion | FLOW-3D AM
Micro and meso scale simulations using FLOW-3D AM help us understand the mixing of different materials in the melt pool and the formation of potential defects such as lack of fusion and porosity. In this simulation, the stainless steel and aluminum powders have independently-defined temperature dependent material properties that FLOW-3D AM tracks to accurately capture the melt pool dynamics. Learn more about FLOW-3D AM’s mutiphysics simulation capabilities at https://www.flow3d.com/products/flow3…
FLOW-3D WELD 는 레이저 용접 공정에 대한 강력한 통찰력을 제공하여 공정 최적화를 달성합니다. 더 나은 공정 제어로 다공성, 열 영향 영역을 최소화하고 미세 구조 진화를 제어 할 수 있습니다. 레이저 용접 공정을 정확하게 시뮬레이션하기 위해 FLOW-3D WELD 는 레이저 열원, 레이저-재료 상호 작용, 유체 흐름, 열 전달, 표면 장력, 응고, 다중 레이저 반사 및 위상 변화를 특징으로 합니다.
FLOW-3D는 물고기 통로, 댐 파손, 배수로, 눈사태, 수력 발전 및 기타 수자원 및 환경 공학 과제 모델링을 포함하여 유압 산업에 대한 많은 응용 분야를 가지고 있습니다. 엔지니어는 수력 발전소의 기존 인프라 용량을 늘리고, 어류 통로, 수두 손실을 최소화하는 흡입구, 포 이베이 설계 및 테일 레이스 흐름을위한 개선 된 설계를 개발하고, 수세 및 퇴적 및 공기 유입을 분석 할 수 있습니다.
FLOW-3D CAST 에는 캐스팅을 위해 특별히 설계된 광범위하고 강력한 물리적 모델이 포함되어 있습니다. 이러한 특수 모델에는 lost foam casting, non-Newtonian fluids, and die cycling에 대한 알고리즘이 포함됩니다. FLOW-3D CAST 의 강력한 시뮬레이션 엔진과 결함 예측을 위한 새로운 도구는 설계주기를 단축하고 비용을 절감 할 수 있는 통찰력을 제공합니다.
HPDC |Comparison of slow shot profiles and entrained air during a filling simulation |FLOW-3D CAST
Shown is a video comparing two slow shot profiles. The graphs highlight the shot profiles through time and the difference in entrained air between the slow shots. Note the lack of air entrained in shot sleeve with calculated shot profile which yields a much better controlled flow within the shot sleeve.
FLOW-3D는 선박 설계, 슬로싱 다이내믹스, 파동 충격 및 환기 등 연안 및 해양 애플리케이션에 이상적인 소프트웨어입니다. 연안 애플리케이션의 경우 FLOW-3D는 연안 구조물에 심각한 폭풍과 쓰나미 파장의 세부 정보를 정확하게 예측하고 플래시 홍수 및 중요 구조물 홍수 및 손상 분석에 사용됩니다.
FLOW-3D POST 2023R2는 EXODUS II 형식을 기반으로 하는 완전히 새로운 결과 파일 형식을 도입하여 더 빠른 후처리를 가능하게 합니다. 이 새로운 파일 형식은 크고 복잡한 시뮬레이션의 후처리 작업에 소요되는 시간을 크게 줄이는 동시에(평균 최대 5배!) 다른 시각화 도구와의 연결성을 향상시킵니다.
FLOW-3D POST 2023R2 에서 사용자는 이제 선택한 데이터를 flsgrf 또는 EXODUS II 파일 형식으로 쓸 수 있습니다 . 새로운 EXODUS II 파일 형식은 각 개체에 대해 유한 요소 메시를 활용하므로 사용자는 다른 호환 가능한 포스트 프로세서 및 FEA 코드를 사용하여 FLOW-3D HYDRO 결과를 열 수도 있습니다. 새로운 워크플로우를 통해 사용자는 크고 복잡한 사례를 신속하게 시각화하고 임의 슬라이싱, 볼륨 렌더링 및 통계를 사용하여 보조 정보를 추출할 수 있습니다.
FLOW-3D POST 의 새로운 EXODUS II 파일 형식에서 볼륨 렌더링 기능을 사용하여 동반된 공기를 보여주는 예입니다 .
새로운 결과 파일 형식은 hydr3d 솔버의 성능을 저하시키지 않으면서 flsgrf 에 비해 시각화 작업 흐름에서 놀라운 속도 향상을 자랑합니다. 이 흥미로운 새로운 개발은 결과 분석의 속도와 유연성이 향상되어 원활한 시뮬레이션 경험을 제공합니다.
FLOW-3D HYDRO 2023R2는 2방정식(RANS) 난류 모델에 대한 동적 혼합 길이 계산을 크게 개선했습니다. 거의 층류 흐름 체계와 같은 특정 제한 사례에서는 이전 버전의 코드 계산 제한기가 때때로 과도하게 예측되어 사용자가 특정 혼합 길이를 수동으로 입력해야 할 수 있습니다.
새로운 동적 혼합 길이 계산은 이러한 상황에서 난류 길이와 시간 규모를 더 잘 설명하며, 이제 사용자는 고정(물리 기반) 혼합 길이를 설정하는 대신 더 넓은 범위의 흐름에 동적 모델을 적용할 수 있습니다.
적절한 고정 혼합 길이와 비교하여 접촉 탱크의 혼합 시뮬레이션을 위한 기존 동적 혼합 길이 모델과 새로운 동적 혼합 길이 모델 간의 비교
정수압 초기화
사용자가 미리 정의된 유체 영역에서 정수압을 초기화해야 하는 경우가 많습니다. 이전에는 대규모의 복잡한 시뮬레이션에서 정수압 솔버의 수렴 속도가 느려지는 경우가 있었습니다. FLOW-3D HYDRO 2023R2는 정수압 솔버의 성능을 크게 향상시켜 전처리 단계에서 최대 6배 빠르게 수렴할 수 있도록 해줍니다.
확장된 지형 표현 지원
GeoTIFF 지원
2023R2 릴리스에서 FLOW-3D HYDRO는 기본적으로 래스터 지형 및 수심 측량을 위한 GeoTIFF(.tif) 파일 형식을 지원합니다. 이제 사용자는 GeoTIFF 파일을 사용자 인터페이스로 직접 가져올 수 있습니다.
FLOW-3D HYDRO 에서 렌더링된 GeoTIFF(.tif) 래스터 파일의 예
LandXML 지원
측량 데이터가 균일하지 않거나 래스터 표면의 해상도가 충분하지 않은 경우 TIN 표면은 LandXML(.xml) 파일 형식을 통해 향상된 지형 지도를 제공합니다. FLOW-3D HYDRO 2023R2는 기본적으로 LandXML 파일을 가져옵니다.
래스터 파일과의 향상된 상호 작용
래스터 파일은 고해상도에서 넓은 지형 영역을 다루는 경우가 많으므로 사용자 인터페이스에서 3D 표현의 상호 작용 속도가 느려질 수 있습니다. 이제 사용자는 3D 표현의 품질을 제어하여 렌더링 시간을 크게 줄이고 상호 작용성을 크게 향상시킬 수 있습니다.
FLOW-3D HYDRO 2023R1 의 새로운 기능
FLOW-3D 소프트웨어 제품군의 모든 제품은 2023R1에서 IT 관련 개선 사항을 받았습니다.
FLOW-3D HYDRO 2023R1은 이제 Windows 11 및 RHEL 8을 지원합니다. Linux 설치 프로그램은 누락된 종속성을 보고하도록 개선되었으며 더 이상 루트 수준 권한이 필요하지 않으므로 설치가 더 쉽고 안전해집니다. 그리고 워크플로를 자동화한 분들을 위해 입력 파일 변환기에 명령줄 인터페이스를 추가하여 스크립트 환경에서도 워크플로가 업데이트된 입력 파일로 작동하는지 확인할 수 있습니다.
천수(shallow water) 난류 모델
난류는 물과 환경 흐름장의 주요 측면이며, 특히 천수(shallow water) 근사치로 모델링된 영역에서는 더욱 그렇습니다. 우리는 모델링 위험을 줄이고 더 나은 결과를 제공하기 위해 세 가지 새로운 난류 모델, 일정한 확산도, 혼합 길이 및 Smagorinsky 모델을 포함하도록 천수(shallow water) 모델의 난류 처리를 개선했습니다.
FLOW-3D HYDRO 2022R2 의 새로운 기능
FLOW-3D HYDRO 2022R2 출시로 Flow Science는 FLOW-3DHYDRO 의 워크스테이션과 HPC 버전을 통합하여 단일 노드 CPU 구성에서 다중 노드에 이르기까지 모든 유형의 하드웨어 아키텍처를 활용할 수 있는 단일 솔버 엔진을 제공했습니다. 병렬 고성능 컴퓨팅 실행. 추가 개발에는 향상된 공기 동반 기능과 물 및 환경 응용 분야에 대한 경계 조건 정의 개선이 포함됩니다.
통합 솔버
우리는 FLOW-3D 제품을 단일 통합 솔버로 마이그레이션하여 로컬 워크스테이션이나 고성능 컴퓨팅 하드웨어 환경에서 원활하게 실행했습니다.
많은 사용자가 노트북이나 로컬 워크스테이션에서 모델을 실행하지만, 고성능 컴퓨팅 클러스터에서도 더 큰 모델을 실행합니다. 2022R2 릴리스에서는 통합 솔버를 통해 사용자가 HPC 솔루션의 OpenMP/MPI 하이브리드 병렬화와 동일한 이점을 활용하여 워크스테이션과 노트북에서 실행할 수 있습니다.
증가하는 CPU 코어 수를 사용한 성능 확장의 예2소켓 워크스테이션에서 OpenMP/MPI 하이브리드 병렬화를 위한 메시 분해의 예
멀티 소켓 워크스테이션
다중 소켓 워크스테이션은 이제 매우 일반적이며 대규모 시뮬레이션을 실행할 수 있습니다. 새로운 통합 솔버를 사용하면 이러한 유형의 하드웨어를 사용하는 사용자는 일반적으로 HPC 클러스터 구성에서만 사용할 수 있었던 OpenMP/MPI 하이브리드 병렬화를 활용하여 모델을 실행할 수 있어 성능이 향상되는 것을 확인할 수 있습니다.
낮은 수준의 루틴으로 향상된 벡터화 및 메모리 액세스
대부분의 테스트 사례에서 10~20% 정도의 성능 향상이 관찰되었으며 일부 사례에서는 20%를 초과하는 런타임 이점이 나타났습니다.
정제된 체적 대류 안정성 한계
시간 단계 안정성 제한은 모델 런타임의 주요 동인이며, 2022R2에서는 새로운 시간 단계 안정성 제한인 3D 대류 안정성 제한을 숫자 위젯에서 사용할 수 있습니다. 실행 중이고 대류가 제한된(cx, cy 또는 cz 제한) 모델의 경우 새 옵션은 일반적인 속도 향상을 30% 정도 보여줍니다.
압력 솔버 프리컨디셔너
경우에 따라 까다로운 흐름 구성의 경우 과도한 압력 솔버 반복으로 인해 실행 시간이 길어질 수 있습니다. 이러한 어려운 경우 2022R2에서는 모델이 너무 많이 반복되면 FLOW-3D가 자동으로 새로운 사전 조절기를 활성화하여 압력 수렴을 돕습니다. 테스트의 런타임은 1.9에서 335까지 더 빨라졌습니다!
점탄성 유체에 대한 로그 형태 텐서 방법
점탄성 유체에 대한 새로운 솔버 옵션을 사용자가 사용할 수 있으며 특히 높은 Weissenberg 수에 효과적입니다.
FLOW-3D HYDRO 경계 조건 개선
FLOW-3D HYDRO 2022R2 에서는 물 적용 경계 조건에 대한 두 가지 개선 사항을 사용할 수 있습니다 . 천수(shallow water)의 유량 경계 조건이 개선되어 보다 현실적이고 공간적으로 변화하는 속도 프로파일을 생성하므로 사용자는 정확도를 잃지 않고 도메인 크기를 줄일 수 있습니다. 자연적인 입구 경계 조건의 경우 정격 곡선 완화 시간 옵션을 사용하여 과도 조건에 대한 응답을 향상시킬 수 있습니다.
입구 경계에서 흐름 방향으로 변하는 속도 프로파일의 예
향상된 공기 동반 기능
디퓨저 및 유사한 기포 흐름 응용 분야의 경우 이제 질량 공급원을 사용하여 물기둥에 공기를 유입할 수 있습니다. 또한, 동반 공기 및 용존 산소의 난류 확산에 대한 기본값이 업데이트되었습니다.
최근 FLOW Science, Inc에서는 토목 및 환경 엔지니어링 산업을위한 완벽한 CFD 모델링 솔루션인 FLOW-3D HYDRO 제품을 출시했습니다. 기존 FLOW-3D 사용자이거나 유압 엔지니어링 관행에 CFD 모델링 기능을 사용하시는 것에 관심이 있는 경우, 언제든지 아래 연락처로 연락주세요. 연락처 : 02-2026-0442 이메일 : flow3d@stikorea.co.kr
FLOW-3D HYDRO 는 더 높은 수준의 정확도와 모델 해상도를 제공하기 위해 3D 비 유압 모델링 기능이 필요한 경우 고급 모델링 도구로 사용할 수 있습니다. 일반적인 모델링 응용 분야는 소형 댐 / 인프라, 운송 수력학, 복잡한 3D 하천 수력학, 열 부력 연기, 배수구 및 오염 물질 수송과 관련됩니다.
FLOW-3D HYDRO의 핵심 기능은 전체 3D 모델과 동적으로 연결될 수있는 천수(shallow water) 모델입니다.
이 기능을 통해 사용자는 멀티 스케일 모델링 애플리케이션을위한 모델 도메인을 확장하여 필요한 모델 해상도로 계산 효율성을 극대화 할 수 있습니다. FLOW-3DHYDRO 또한 강 및 환경 응용 분야에 특화된 추가 기능과 고급 물리학을 포함합니다.
시뮬레이션 템플릿
FLOW-3D HYDRO 의 작업 공간 템플릿으로 시간을 절약하고 실수를 방지하며 일관된 모델을 실행하십시오 . 작업 공간 템플릿은 일반적인 응용 분야에 대한 유체 속성, 물리적 모델, 수치 설정 및 시뮬레이션 출력을 미리로드합니다.
작업 공간 템플릿은 7 가지 모델 클래스에 사용할 수 있습니다.
자유 표면 – TruVOF (기본값)
공기 유입
열 기둥
퇴적물 수송
천수(shallow water)
자유 표면 – 2 유체 VOF
자유 표면 없음
사전로드 된 예제 시뮬레이션
FLOW-3D HYDRO 의 40 개 이상의 사전로드 된 물 중심 예제 시뮬레이션 라이브러리는 애플리케이션 모델링을위한 훌륭한 시작점을 제공합니다. 사전로드 된 예제 시뮬레이션은 모델러에게 모델 설정 및 모범 사례의 로드맵뿐만 아니라 대부분의 애플리케이션에 대한 자세한 시작점을 제공합니다.이전다음
비디오 튜토리얼
비디오 자습서는 새로운 사용자가 다양한 응용 프로그램을 모델링하는 방법을 빠르게 배울 수있는 훌륭한 경로를 제공합니다. FLOW-3D HYDRO 비디오 튜토리얼 기능 :
광범위한 응용 및 물리학을위한 AZ 단계별 기록
“사용 방법”정보
모범 사례를위한 팁
CAD / GIS 데이터, 시뮬레이션 파일 및 후 처리 파일
고급 솔버 개발
Tailings Model
새로운 Tailings Model은 tailings dam failure로 인한 tailings runout을 시뮬레이션하기위한 고급 기능을 제공합니다. tailings정의에 대한 다층 접근 방식과 함께 미세하고 거친 입자 구성을 나타내는 이중 모드 점도 모델은 모든 방법으로 건설 된 tailings 댐의 모델링을 허용합니다.
천수(shallow water), 3D 및 하이브리드 3D / 천수(shallow water) 메싱을 포함한 유연한 메싱을 통해 얕은 지역에서 빠른 솔루션을 제공하면서 다층 tailings의 복잡성을 정확하게 모델링 할 수 있습니다. 점성 경계층의 정확한 표현을 위해 천수(shallow water) 메시에 2 층 Herschel-Bulkley 점도 모델을 사용할 수 있습니다.
모델 하이라이트
미세 입자 및 거친 입자 광미 조성물을위한 이중 모드 점도 모델
침전, 패킹 및 입자 종의 난류 확산을 포함한 Tailings 수송
천수(shallow water) 메시를위한 2 층 Herschel-Bulkley 점도 모델
3D, 천수(shallow water), 3D / 천수(shallow water) 하이브리드 메시를 포함한 유연한 메시 접근 방식
Multi-layer, variable composition tailings for general definition of tailings dam construction
Shallow Water
FLOW-3D HYDRO 의 천수(shallow water) 모델링 기능은 3D 메시를 천수(shallow water) 메시와 결합하여 탁월한 모델링 다양성을 제공하는 고유 한 하이브리드 메시를 사용합니다. 압력 솔버의 수치 개선으로 더 안정적이고 빠른 시뮬레이션이 가능합니다. 하이브리드 메쉬의 하단 전단 응력 계산이 크게 향상되어 정확도가 더욱 향상되었습니다. 지형에 거칠기를 적용하는 새로운 방법에는 Strickler, Chezy, Nikuradse, Colebrook-White, Haaland 및 Ramette 방정식이 포함됩니다.
Two-Fluid VOF Model
sharp 인터페이스가 있거나 없는 압축 가능 또는 비압축성 2 유체 모델은 항상 1 유체 자유 표면 모델과 함께 FLOW-3D 에서 사용할 수 있습니다 . 사실, sharp 인터페이스 처리는 TruVOF 기술을 자유 표면 모델과 공유하며 상용 CFD 소프트웨어에서 고유합니다. 최근 개발에는 2- 필드 온도 및 인터페이스 슬립 모델이 포함되었습니다. 이 모델은 오일 / 물, 액체 / 증기, 물 / 공기 및 기타 2 상 시스템에 성공적으로 적용되었습니다.
FLOW-3D HYDRO 는 2- 유체 솔루션의 정확성과 안정성에서 두 가지 중요한 발전을보고 있습니다. 운동량과 질량 보존 방정식의 강화 된 결합은 특히 액체 / 기체 흐름에서 계면에서 운동량 보존을 향상시킵니다. 연속성 방정식에서 제한된 압축성 항의 확장 된 근사값은 더 빠르고 안정적인 2 유체 압력 솔버를 만듭니다.
예를 들어, 터널 및 드롭 샤프트 설계와 같은 유압 응용 분야에서 공기가 종종 중요한 역할을 하기 때문에 두 개발 모두 FLOW-3DHYDRO 릴리스에 적시에 적용됩니다. 일반적으로 낮은 마하 수로 인해 이러한 경우 물과 공기에 제한된 압축성이 사용됩니다.
고성능 컴퓨팅 및 클라우드
일반 워크스테이션 또는 랩톱으로 많은 작업을 수행 할 수 있지만, 대형 시뮬레이션과 고화질 시뮬레이션은 더 많은 CPU 코어를 활용함으로써 엄청난 이점을 얻을 수 있습니다. FLOW-3D CLOUD 및 고성능 컴퓨팅은 더 빠르고 정확한 모델을 실행할 수있는 더 빠른 런타임과 더 많은 선택권을 제공합니다.
재료 비용을 줄이고 사이클 시간을 개선하기 위해 소비재 회사는 슬로 싱, 튀기 및 공기 혼입을 포함한 많은 자유 표면 유체 문제를 처리해야합니다.
Predicting Entrained Air in a Bottle Filling Example
혼입된 공기는 생산 라인에서 컨테이너가 채워질 때 액체의 부피를 증가시킬 수 있습니다. 아래 왼쪽 이미지는 높이가 약 20cm인 병에 1.2 초 동안 채우는 것을 보여줍니다. 색상 음영은 액체에서 공기의 부피 비율을 나타냅니다. 병에서 짧은 시간과 높은 수준의 혼합으로 인해 공기가 표면으로 올라와 빠져 나갈 시간이 없었습니다. 그러나 오른쪽 이미지에서 볼 수 있듯이 약 1.7 초의 추가 시간이 지나면 표면으로 상승하는 공기로 인한 액체 부피 감소가 명확하게 보입니다. FLOW-3D의 드리프트 플럭스 모델을 사용하면 액체의 기포와 같은 성분을 분리하여 분리 할 수 있습니다.
Air entrainment (left) and separation of air and liquid (right)
In by 9, out by 5 – Rapid evaluation of Tide® bottle filling
FLOW-3D를 사용하여 새로운 Tide 병 디자인의 채우기를 모델링하는 방법을 설명하는이 기사는 The Procter and Gamble Company의 기술 부문 책임자 인 John McKibben이 기고했습니다.
오전 9시에 긴급한 이메일을 받았다고 상상해보십시오.
새로운 Tide® 병 디자인 중 하나가 핸들을 채우고 충전 장비에 문제가있을 수 있음을 방금 깨달았습니다. 프로토 타입 병도없고 몇 주 동안도 없을 것입니다. 디자이너와 소비자는 디자인의 모양을 좋아하지만 그것이 채우는 방식은 우리 생산 시설의 쇼 스토퍼가 될 수 있습니다. 이 상황을 접했을 때 저는 3D 지오메트리 (그림 1)의 스테레오 리소그래피 (.stl) 파일을 요청하여 응답을 시작했고 제가 할 수있는 일을 확인했습니다. FLOW-3D는 .stl 파일을 사용하여 지오메트리를 입력 할 수 있으며 채우기에 대한 자유 표면 문제를 해결할 수 있어야한다는 것을 알고있었습니다. 나는 이것이 잠재적 인 문제에 대한 좋은 질적 이해를 제공 할 것으로 기대했지만,이 응용 프로그램에 대해 얼마나 정확한지에 대해서는 약간 불확실했습니다.
Setting up and Running the Simulation
오후 1 시경에 지오메트리 파일, 유량 및 유체 속성을 받았습니다. 몇 시간 내에 시뮬레이션이 실행되어 예비 결과를 제공했습니다. 저는 제 고객을 초대하여 결과를 간단히 살펴 보았고 그는 “보스의 상사”도 함께 살펴 보았습니다. 그래서 저녁 5 시까 지 예비 결과를보고 원래 우려 사항이 문제가 아니라고 판단했습니다.
그러나 결과는 몇 가지 다른 질문을 제기했습니다. 핸들을 채우면 유입되는 유체 분사가 많이 분리되었습니다. 나는 이것이 동반 된 공기와 거품의 양을 증가시킬 것이라는 것을 알고 있었다 (우리는 결국 세탁 세제를 채우고있다). FLOW-3D 공기 혼입 모델을 테스트하기로 결정했습니다. 이 모델은 원래 난류 제트 용으로 개발되었으며,이 층류 문제를 볼 때 얼마나 잘 수행 될지 확신 할 수 없었습니다.
Figure 2: Filled resultsExperimental comparison of bottle filling model with and without the air entrainment model, courtesy of The Procter & Gamble Company.
그림 2는 공기 유입 모델이 있거나없는 병 충전 모델의 결과를 보여줍니다. 혼입 된 공기가 포함되면 충전 레벨이 크게 증가합니다. 혼입 된 공기가 병 상단에서 액체를 밀어 내지는 않지만 공기 혼입 정확도를 확인해야 할만큼 충분히 가깝습니다.
그림 3은 몇 주 후에 실행 된 실험의 이미지와 공기 혼입 수준을 비교합니다 (시제품 병이 제공되었을 때). 제트 분리 및 충진 수준의 정 성적 일치는 우수하며 시뮬레이션이 병 설계를 선별하기에 충분히 정확하다는 것을 확인했습니다.
원심 주조 Workspace는 원심 주조 사용자에게 수평 및 수직 진정한 원심 주조, 부분 원심 주조 및 원심 주조 시뮬레이션을 위한 편리한 도구를 제공합니다. 새로운 원심 주조 Workspace를 사용하면 사용자가 프로세스를 모델링하고 설계 매개 변수를 최적화하는데 필요한 모든 도구를 찾을 수 있습니다. 금형을 고정시키고 회전하는 메쉬를 통해 사용자는 ladle 붓기를 포함하여 상상할 수 있는 모든 금형 모션을 모델링할 수 있는 유연성을 제공합니다.
원통형 메싱은 가능한 최고의 흐름 모델링 정확도를 제공하는 반면, 다중 블록 메싱은 주조물의 높은 전단 및 고온 구배 영역에서 훨씬 더 높은 정확도를 위한 효율적인 방법을 제공합니다. 이 솔루션은 적합하지 않은 금형 회전 속도에 따라 비처럼 떨어지는 것과 같은 흐름 관련 문제, 공기 유입 또는 응고 부위의 재용해과 같은 결함을 예측합니다. 몰드 예열 온도, 냉각 구성 및 금형 회전률과 같은 프로세스 매개변수는 모두 모델 설정의 일부가 될 수 있습니다.
FLOW-3D v12는 그래픽 사용자 인터페이스 (GUI)의 설계 및 기능에서 매우 큰 변화를 이룬 제품으로 모델 설정을 단순화하고 사용자 워크 플로를 향상시킵니다. 최첨단 Immersed Boundary Method(침수경계 방법)은 FLOW-3Dv12 솔루션의 정확성을 높여줍니다. 다른 주요 기능으로는 슬러지 침강 모델, 2-Fluid 2-Temperature 모델 및 Steady State Accelerator가 있으며,이를 통해 사용자는 자유 표면 흐름을 더욱 빠르게 모델링 할 수 있습니다.
Physical and Numerical Model
Immersed boundary method
힘과 에너지 손실에 대한 정확한 예측은 고체 주위의 흐름과 관련된 많은 엔지니어링 문제를 모델링하는 데 중요합니다. 새 릴리스 FLOW-3D v12에는 이러한 문제점 해결을 위해 설계된 새로운 고스트 셀 기반 Immersed Boundary Method (IBM)가 있습니다. IBM은 내 외부 흐름 해석을 위해, 벽 근처에서 보다 정확한 해를 제공하여 드래그 앤 리프트 힘의 계산을 향상시킵니다.힘과 에너지 손실의 정확한 예측은 고체 주위의 흐름을 포함하는 많은 공학적 문제를 모델링 하는데 중요합니다.
Two-field temperature for the two-fluid model
2 유체 열전달 모델은 각 유체에 대한 에너지 전달 방정식을 분리하기 위해 확장되었습니다. 각 유체는 이제 자체 온도 변수를 가지므로 인터페이스 근처의 열 및 물질 전달 솔루션의 정확도가 향상됩니다. 인터페이스에서의 열전달은 이제 시간의 표 함수가 될 수 있는 사용자 정의 열전달 계수에 의해 제어됩니다.
새로운 슬러지 정착 모델은 수처리 애플리케이션에 부가되어 사용자들이 수 처리 탱크와 클래리퍼의 고형 폐기물 역학을 모델링 할 수 있게 해 줍니다. 침전 속도가 분산상의 액적 크기의 함수 인 드리프트-플럭스 모델과 달리, 침전 속도는 슬러지 농도의 함수이며 기능 및 표 형식으로 입력 할 수 있습니다.
Void particles 가 기포 및 상 변화 모델에 추가되었습니다. Void particles는 붕괴 된 Void 영역을 나타내며, 항력 및 압력을 통해 유체와 상호 작용하는 작은 기포로 작용합니다. 주변 유체 압력에 따라 크기가 변하고 시뮬레이션이 끝날 때의 최종 위치는 공기 유입 가능성을 나타냅니다.
Sediment scour model
퇴적물 수송 및 침식 모델은 정확성과 안정성을 향상시키기 위해 정비되었습니다. 특히 퇴적물 종의 질량 보존이 크게 개선되었습니다.
고정 압력 경계 조건에는 압력 및 유체 분율을 제외한 모든 유량이 해당 경계의 상류의 유량 조건을 반영하는 ‘유출’옵션이 포함됩니다. 유출 압력 경계 조건은 고정 압력 및 연속 경계 조건의 하이브리드입니다.
Moving particle sources
시뮬레이션 중에 입자 소스를 이동할 수 있습니다. 시간에 따른 병진 및 회전 속도는 표 형식으로 정의됩니다. 입자 소스의 운동은 소스에서 방출 된 입자의 초기 속도에 추가됩니다.
Variable center of gravity
기변 무게중심은 중력 및 비관 성 기준 프레임 모델에서, 시간의 함수로서 무게 중심의 위치는 외부 파일에서 테이블로서 정의 될 수있다. 이 기능은 연료를 소비하고 분리 단계를 수행하는 로켓과 같은 모형을 모델링 할 때 유용합니다.
공기 유입 모델
가장 간단한 부피 기반 공기 유입 모델 옵션이 기존 질량 기반 모델로 대체되었습니다. 질량 기반 모델은 부피와 달리 주변 유체 압력에 따라 부피가 변화하는 동안 흡입된 공기량이 보존되기 때문에 물리학적 모델입니다.
Tracer diffusion
유동 표면에서 생성된 추적 물질은 분자 및 난류 확산 과정에 의해 확산될 수 있으며, 예를 들어 실제 오염 물질의 동작을 모방한다.
Model Setup
Simulation units
온도를 포함하여 단위 시스템은 완전히 정의해야하는데 표준 단위 시스템이 제공됩니다. 또한 사용자는 다양한 옵션 중에서 질량, 시간 및 길이 단위를 정의 할 수 있으므로 사용자 정의가 가능한 편리한 단위를 사용할 수 있습니다. 사용자는 압력이 게이지 또는 절대 단위로 정의되는지 여부도 지정해야합니다. 기본 시뮬레이션 단위는 기본 설정에서 설정할 수 있습니다. 단위를 완전히 정의하면 FLOW-3D 가 물리량의 기본값을 정의하고 범용 상수를 설정하여 사용자가 요구하는 작업량을 최소화 할 수 있습니다.
Shallow water model
Manning’s roughness in shallow water model
Manning의 거칠기 계수는 지형 표면의 전단 응력 평가를 위해 천수(shallow water) 모델에서 구현되었습니다. 표면 결함의 크기를 기반으로 기존 거칠기 모델을 보완하며 이 모델과 함께 사용할 수 있습니다. 표준 거칠기와 마찬가지로 매닝 계수는 구성 요소 또는 하위 구성 요소의 속성이거나 지형 래스터 데이터 세트에서 가져올 수 있습니다.
Mesh generation
하단 및 상단 경계 좌표의 정의만으로 수직 방향의 메시 설정이 단순화되었습니다.
Component transformations
사용자는 이제 여러 하위 구성 요소로 구성된 구성 요소에 회전, 변환 및 스케일링 변환을 적용하여 복잡한 형상 어셈블리 설정 프로세스를 단순화 할 수 있습니다. GMO (General Moving Object) 구성 요소의 경우, 이러한 변환을 구성 요소의 대칭 축과 정렬되도록 신체에 맞는 좌표계에 적용 할 수 있습니다.
Changing the number of threads at runtime
시뮬레이션 중에 솔버가 사용하는 스레드 수를 변경하는 기능이 런타임 옵션 대화 상자에 추가되어 사용 가능한 스레드를 추가하거나 다른 태스크에 자원이 필요한 경우 스레드 수를 줄일 수 있습니다.
Probe-controlled heat sources
활성 시뮬레이션 제어가 형상 구성 요소와 관련된 heat sources로 확장되었습니다. 히스토리 프로브로 열 방출을 제어 할 수 있습니다.
Time-dependent temperature at sources
질량 및 질량 / 운동량 소스의 유체 온도는 이제 테이블 입력을 사용하여 시간의 함수로 정의 할 수 있습니다.
Emissivity coefficients
공극으로의 복사 열 전달을위한 방사율 계수는 이제 사용자가 방사율과 스테판-볼츠만 상수를 지정하도록 요구하지 않고 직접 정의됩니다. 후자는 이제 단위 시스템을 기반으로 솔버에 의해 자동으로 설정됩니다.
Output
등속 필드 솔버 옵션을 사용할 때 유량 속도를 선택한 데이터 로 출력 할 수 있습니다 .
벽 접착력으로 인한 지오메트리 구성 요소의 토크 는 기존 벽 접착력의 출력 외에도 일반 이력 데이터에 별도의 수량으로 출력됩니다.
난류 모델 출력이 요청 될 때 난류 에너지 및 소산과 함께 전단 속도 및 y +가 선택된 데이터로 자동 출력됩니다 .
공기 유입 모델 출력에 몇 가지 수량이 추가되었습니다. 자유 표면을 포함하는 모든 셀에서 혼입 된 공기 및 빠져 나가는 공기의 체적 플럭스가 재시작 및 선택된 데이터로 출력되어 사용자에게 공기가 혼입 및 탈선되는 위치 및 시간에 대한 자세한 정보를 제공합니다. 전체 계산 영역 및 각 샘플링 볼륨 에 대해이 두 수량의 시간 및 공간 통합 등가물 이 일반 히스토리 로 출력됩니다.
솔버의 출력 파일 flsgrf 의 최종 크기 는 시뮬레이션이 끝날 때보 고됩니다.
2 유체 시뮬레이션의 경우, 기존의 출력 수량 유체 체류 시간 및 유체 가 이동 한 거리는 이제 유체 # 1 및 # 2와 유체의 혼합물에 대해 별도로 계산됩니다.
질량 입자의 경우 각 종의 총 부피와 질량이 계산되어 전체 계산 영역, 샘플링 볼륨 및 플럭스 표면에 대한 일반 히스토리 로 출력되어 입자 종 수에 대한 현재 출력을 보완합니다.
예를 들어 사용자가 가스 미순환을 식별하고 연료 탱크의 환기 시스템을 설계하는 데 도움이 되도록 마지막 국부적 가스 압력이 옵션 출력량으로 추가되었습니다. 이 양은 유체가 채워지기 전에 셀의 마지막 간극 압력을 기록하며, 단열 버블 모델과 함께 사용됩니다.
New Customizable Source Routines
사용자 정의 가능한 새로운 소스 루틴이 추가되었으며 사용자의 개발 환경에서 액세스 할 수 있습니다.
소스 루틴 이름
설명
cav_prod_cal
캐비 테이션 생산 및 확산 속도
sldg_uset
슬러지 정착 속도
phchg_mass_flux
증발 및 응축에 의한 질량 흐름
flhtccl
유체#1과#2사이의 열 전달 계수
dsize_cal
2상 유동에서의 동적 낙하 크기 모델의 충돌 및 이탈율
elstc_custom.
점탄성 유체에 대한 응력 방정식의 소스 용어
Brand New User Interface
FLOW-3D의 사용자 인터페이스가 완전히 재설계되어 사용자의 작업 흐름을 획기적으로 간소화하는 최신의 타일 구조를 제공합니다.
Dock widgets 설정
Physics, Fluids, Mesh 및 FAVOR ™를 포함한 모든 설정 작업이 형상 창 주위의 dock widgets으로 변환되어 모델 설정을 단일 탭으로 압축 할 수 있습니다. 이 전환을 통해 이전 버전의 복잡한 트리가 훨씬 깔끔하고 효율적인 메뉴 표시로 바뀌어 모델 설정 탭을 떠나지 않고도 모든 매개 변수에 쉽게 액세스 할 수 있습니다.
New Model Setup icons With our new Model Setup design comes new icons, representing each step of the setup process.New Physics icons Our Physics icons are designed to be easily differentiated from one another at a glance, while providing clear visual representation of each model’s purpose and use.
RSS feed
새 RSS 피드부터 FLOW-3D v12.0 의 시뮬레이션 관리자 탭이 개선되었습니다 . FLOW-3D 를 시작하면 사용자에게 Flow Science의 최신 뉴스, 이벤트 및 블로그 게시물이 표시됩니다.
Configurable simulation monitor
시뮬레이션을 실행할 때 중요한 작업은 모니터링입니다. FLOW-3Dv12.0에서는 사용자가 시뮬레이션을 더 잘 모니터링할 수 있도록 Simulation Manager의 플로팅 기능이 향상되었습니다. 사용자는 시뮬레이션 런타임 그래프를 통해 모니터링할 사용 가능한 모든 일반 기록 데이터 변수를 선택하고 각 그래프에 여러 변수를 추가할 수 있습니다. 이제 런타임에서 사용할 수 있는 일반 기록 데이터는 다음과 같습니다.
최소/최대 유체 온도
프로브 위치의 온도
유동 표면 위치에서의 유량
시뮬레이션 진단(예:시간 단계, 안정성 한계)
Runtime plots of the flow rate at the gates of the large dam / Large dam with flux surfaces at the gates
Conforming mesh visualization
사용자는 이제 새로운 FAVOR ™ 독 위젯을 통해 적합한 메쉬 블록을 시각화 할 수 있습니다 .
Large raster and STL data
데이터를 처리하는 데 걸리는 시간으로 인해 큰 형상 데이터를 처리하는 것은 어려울 수 있습니다. 대형 지오메트리 데이터를 처리하는 데 여전히 상당한 시간이 소요될 수 있지만 FLOW-3D는 이제 이러한 대형 데이터 세트를 백그라운드 작업으로로드하여 사용자가 데이터를 처리하는 동안 완벽하게 응답하고 중단없는 인터페이스에서 계속 작업 할 수 있습니다.
FLOW-3D 의 공기 혼입 모델(air entrainment model)은 자유 표면에서 용해되지 않은 공기 혼입을 시뮬레이션하는 강력한 도구입니다. 제트 및 방수로 충돌시 관찰되는 국부적이고 난류가없는 자유 표면 혼입 기능이 있습니다. 이러한 기능은 엔지니어가 설계시 공기 유입을 예측하고, 공기유입이 안전하게 작동하도록 적절한 수정을 할 수 있게 합니다.
Spillway hydraulics / 여수로 수리장치
여수로 구조는 다양한 작동 조건을 처리 할 수 있도록 설계되어야 합니다. 유동 조건이 설계 범위의 상단에 도달하면 여수로 표면의 불규칙성으로 인해 유동이 분리 될 수 있습니다. 이는 여수로 표면의 압력이 캐비테이션을 일으킬 정도로 낮아지게 합니다. 캐비테이션은 구조물의 강도에 매우 해로우며 치명적인 손상을 초래할 수 있습니다.
공기 유입은 캐비테이션의 가능성을 줄이는 수단입니다. 물이 공기에 존재하면 캐비테이션 영역의 붕괴하는 기포에 감쇠 효과를 추가하여 캐비테이션 손상을 줄입니다. 여수로의 속도가 충분히 높으면 공기를 동반시키고 캐비테이션을 줄이기 위해 폭기 장치를 추가해야합니다.
폭기 흐름의 시뮬레이션과 폭기 장치에서 포획된 공기의 예측.
왼쪽 이미지는 거시적 인 밀도에 의해 착색됩니다. 오른쪽 그래프는 폭기 장치에 유입 된 수분의 일정한 부피와 폭기 장치 이후의 수분 및 공기의 양을 비교 한 것입니다.
아래 동영상은 FLOW-3D에서 공기 유입 과정을 시뮬레이션하는 방법을 보여줍니다. 여기에는 공기혼입 및 드리프트 플럭스 모델의 이론에 대한 세부 정보와 FLOW-3D에서 기본 공기혼입 시뮬레이션을 설정하는 방법에 대한 데모가 포함되어 있습니다.
Fish passage design / 물고기 개체수 유지를 위한 어도 설계
공기가 물로 혼입되면 미생물의 성장을 유지하고 건강한 어류 개체군의 생존을 보장 할 수 있습니다. 그러나 과포화 상태의 용존 기체는 수생 생물에 부정적인 영향을 미치는 수질 문제가 됩니다. 공기 동반 모델의 또 다른 용도는 강의 하류로 방출되는 배수로에서 동반되는 공기의 농도를 결정하기 위해 해양 생물학에서 사용됩니다. 이 모델에 대한 더 자세한 정보는 Air Entrainment 의 Flow Science Report를 다운로드하십시오.
LEED (Leadership in Energy and Environmental Design)는 제 3자가 친환경 건축물 인증을 제공하는 자발적 인증 시스템입니다. LEED 자발적 참여는 리더십, 혁신, 환경 보호 및 사회적 책임을 보여 줍니다. LEED는 건물 소유자와 운영자에게 건물의 성능과 수익에 즉시 영향을 미치는데 필요한 도구를 제공하는 동시에 건물 거주자에게 건강한 실내 공간을 제공합니다.
FLOW-3D는 보고타(콜롬비아)의 사무실 건물에서 “IEQ-Credit2–환기 증가”라는 신뢰를 얻는 데 큰 도움을 주었습니다. 이러한 인정을 받기 위해서는 실외 공기가 ASHRAE의 표준 비율인 30%를 초과한다는 것을 증명해야만 합니다. 이 건물에서 실외 공기는 태양 광선에 의해, 가열되는 지붕 위의 2개의 유리 굴뚝에 의해 발생되는 온도 차이에 의해 발생하는 열 부력의 영향으로 제공됩니다. 이것은 바람이 불지 않는 조건에서 이루어져야 합니다.
The model has been easily created using Google SketchUp® and exported directly to FLOW-3D in STL format. The STL model was imported as STL baffles in FLOW-3D.
외부 에너지 시뮬레이션 소프트웨어를 통해 공간 및 열 하중(또는 손실)의 초기 조건을 파악했습니다. 이 소프트웨어는 태양열 복사, 열 관성, 단열, 내부 부하, 단열, 유리 및 건물의 열 동작을 정의하는 기타 모든 매개 변수를 고려합니다. 열 하중은 FLOW-3D에 포함되어 있어 CFD시뮬레이션이 수행될 때 최종 온도 분포를 제공합니다.
<좌측>건물 내부의 공기 흐름이 간소화됩니다. <우측>건물 내부 온도 분포
왼쪽 그림에서 2개의 열 굴뚝은 바람이 전혀 없는 상태에서도 실외 공기 유입을 유발합니다. 환기 효율, 시간당 공기 변화, 국소 난류, 공기 잔류 시간은 설계 프로세스에서 검사할 수 있는 변수 중 일부입니다. 오른쪽 그림에서 FLOW-3D는 건물 내부의 온도 분포를 제공하여 설계 과정에서 열적 쾌적성 평가를 수행할 수 있도록 합니다. 그릴 크기는 모든 공간에서 원하는 쾌적한 온도를 얻을 수 있도록 조정할 수 있습니다.
외부 공기 온도와 기압은 경계 조건(흡입 및 배출 공기 그릴)에서 설정되었습니다. 나머지는 건물 내부에서 상세한 공기 이동을 분석한 것입니다. 시뮬레이션은 에너지가 정상 상태 조건이 달성될 때까지 수행되었다.
건물 내 공기를 올바르게 분배하고 적절한 쾌적한 온도를 확보하기 위해 건물 구조와 그릴 치수를 조정해야 했습니다. 이 과정을 반복한 후 모델을 검증하여 LEED 인증 조건을 충족한 사례입니다.
하수도 시스템은 액션영화의 도피 루트로 사용되지 않는 한 흥미롭지 않을 것입니다. 폭우로 인해 이산화탄소 수치가 올라갈 때까지 여러분은 그것에 대해 생각조차 하지 않을 것입니다. 불행하게도, 770개 이상의 오래 된 미국 도시들 아래에 있는 하수구 시스템은 심한 폭풍으로 오염 문제를 일으킵니다. 이러한 구형 설계는 하수 및 폭풍 유실을 위한 비용 효율적인 단일 스타일 파이프를 사용했으며 연결된 파이프로 강 및 호수에 하수를 내보냅니다(CSO).
1994년 미국 환경보호청(EPA)은 주로 북동부 및 그레이트 레이크 지역의 관련 지방 자치 단체들에게 CSO관련 문제를 줄이거나 제거하도록 하는 정책을 발표했습니다. (2000년 “Clean Water Act”의 일부로 법률화된 정책). 인디애나 폴리스(Indianapolis)는 가벼운 비 폭풍으로 인해 하수 오물의 백업 및 범람이 발생할 수 있는 도시 중 하나였으므로, 주요 건설 조건에서 2025년까지 문제를 해결하는 것이 필요하였습니다.
인디애나 폴리스는 국제 디자인 회사인 AECOM에 Citizens Energy Group이 건설하고 있는 3개의 깊은 암석 저장 터널 중 첫 번째를 설계할 것을 요청했습니다. 총 25마일인 이 시스템은 대규모 지하 펌프장과 기존의 하수구에서 CSO를 수직으로 떨어뜨리는 연결 구조물을 포함합니다. 첫 번째 터널의 경우, 강우가 가라 앉은 후에 3 개의 커다란 강하 구조물이 CSO를 저장 터널로 전환하여 후속 처리를 수행했습니다.
프로젝트를 해결하기 위해 AECOM은 여러 가능한 낙하 구조물 설계의 동작을 시뮬레이션하기 위해 FLOW-3D를 선택하여, 구축 및 평가 예산이 책정 된 물리적 모델에 대한 재 작업의 필요성을 최소화했습니다. 테스트 결과는 예측 값과 일치하였으므로 재설계가 필요하지 않았습니다. 또한, 이제 AECOM은 유압 설계작업의 첫 번째 단계를 일반적으로 CFD시뮬레이션을 사용합니다.
Large Scale Project on a Tight Delivery Schedule
촉박한 납품 일정에 따른 대규모 프로젝트
20세기에 건설된 하수 처리장은 주거용, 상업용, 환경유출물의 유출로 무엇을 해야 할 것인지에 대한 새로운 인식을 가져다 주었습니다. CSO 방전은 정상적으로 운영되는 동안 처리시설로 직접 이동되며 모든 과정이 양호하게 운영됩니다. 불행하게도, 대규모 폭풍이 발생하는 동안, 발전소들의 초과 용량문제를 피하기 위해 인근 수역으로 과도한 유량을 방출합니다. 이들 배출은 기름과 살충제, 야생동물 배설물에 이르기까지 다양한 오염 물질을 포함합니다. 고무적인 성공의 신호로, 1990 년대에 착공된 새로운 CSO 분리, 저장 및 처리 시설로 오염의 영향에 대해 67 %의 개선을 이루었지만, 여전히 많은 연구가 이루어져야 합니다. 인디애나 폴리스의 경우, 인디애나 폴리스시 공공사업부가 CSO 장기 통제계획을 준비한 2008년에 그러한 노력이 시작되었습니다. 정상적인 처리 공장에서 처리 할 수 있을 때까지 오버플로우가 발생하는 “저장 및 운송”접근법의 핵심은 인디애나 폴리스 터널 저장 시스템 또는 인디애나라고 합니다.
이 시스템의 첫번째 단계는 딥 록 터널 커넥터(DRTC)라고 불리는 1억 8천만달러 가치의 프로젝트입니다. DRTC는 길이 7마일의 18피트 직경의 지하 터널로, 기존의 인디애나 폴리스의 3개의 서버 대 계층 유출 연결의 흐름 경로를 다시 만들 것입니다(그림 1). 목표는 과잉 강우 유출을 기존 하수구와 새 터널 사이의 낙하 구조를 통해 이들 대피소에서 거대한 터널로 안전하게 재배치하고, 폭풍 후 처리를 위해 처리장으로 펌핑 될 수있을 때까지 유지합니다.
Fig. 1. City of Indianapolis Deep Rock Tunnel Connector (DRTC), a “storage and transport” concept being built to handle combined sewage overflow (CSO) during heavy storms. Three vertical drop structures will capture this flow and divert it downwards to 18-foot-diameter storage tunnels running more than 250 feet underground; the tunnels store the CSO until sewage treatment plant capacity becomes available. (Image courtesy Citizens Energy Group)
평균적으로 지표면 아래 250피트 깊이에서, DRTC는 건설과 궁극적인 운영 동안 위의 주변 지역에 대한 혼란을 최소화하도록 설계되었습니다. 그러나 이 프로젝트의 규모와 복잡성은 AECOM의 과제에 긴급성을 더했습니다. 세 장소 각각에 대한 가능한 낙하 구조 설계와 평가, 구조물 설계의 60%를 7개월 이내에 마무리 지었습니다.
이러한 구조물의 목적은 표준 도시 하수 시스템에서 깊은 저장 터널로 하수 흐름을 전달하는 동시에, 효율적 손실( 느린 속도 또는 백업)과 장기적인 도심을 방지하는 것입니다. 각 섹션의 크기와 모양이 유입 흐름의 볼륨 및 속도와 세심하게 일치하지 않을 경우 발생할 수 있는 구조적 손상입니다. AECOM의 수석 기술 전문가인 라이언 에디슨 컨설턴트는 계약의 스케줄링 요구 사항이 유효성 검사를 위해서는, 단 하나의 모델에만 물리적 건물과 테스트 활동을 제한할 것이라는 것을 알게되었습니다. 다른 주요 건설 프로젝트에 15년간 FLOW-3D 시뮬레이션 소프트웨어를 사용해 왔기 때문에, 난류, 과전압 및 에너지 낭비를 예측하는 능력은 충분하지 않고 디자인 프로젝트에 적합하다고 자신했습니다. 또한 여러 검증(what-if) 시나리오를 실행하기 위한 소프트웨어 옵션을 통해 설계 세부 사항을 다시 실행해야 하는 위험을 최소화할 수 있었습니다. 변경 사항이 적용될 경우 상당한 이점은 여러개의 병렬 시공 트랙이 있는 프로젝트에 있습니다. 시간 제약에도 불구하고, 에디슨은 특히 이 도전에 만족했습니다. 왜냐하면 “CFD로 드롭 구조 설계를 만들고 물리학에서 이것들은 너무 큰 구조이기 때문입니다.”라고 그는 말합니다. 그것들은 CFD는 실제로 사용되지 않는데 보통 물리적 모델이나 손으로 계산하는 것으로 이루어집니다.
DRTC 프로젝트를 위해서, 그는 먼저 시뮬레이션된 작동 조건에 대해서 컴퓨터 설계를 테스트할 것입니다. 에디슨은 3차원의 일시적이고 격동적인 흐름 조건을 모델링 할 수 있는 소프트웨어 패키지인 FLOW-3D를 사용했습니다. 각 설계에 대한 계산 메쉬를 변경하지 않고도 여러 설계 지오 메트리를 모델링 할 수 있는 기능이였습니다. 시뮬레이션 데이터로 무장한 에디슨은 그 결과를 아이오와 대학교 II. 시설에서 시험한 1:10 크기의 물리적 모델의 작동 데이터와 비교하였습니다. (후자는 원래 아이오와 유압 연구소라고 불렸지만, 지금은 그룹의 다양한 범위를 반영하여 IIHR-Hydroscience & Engineering으로 알려져 있습니다.)
Zeroing in on the Drop-Structure Challenge
드롭 구조 과제에서 영점 조정
가장 제한적인 DRTC 사이트의 지오 메트리는 CSO 008로 지정된 레귤레이터에서 발생합니다. 기존 CSO 레귤레이터(기울기 약 75피트 아래)를 새 18피트 직경의 수집 터널과 연결하려면, 이 위치에서 150피트 이상의 수직 방향 주행이 필요합니다. 각 낙하 구조에 7백만달러 이상이 소요되는 경우, 프로젝트 관리자들은 물리적 모델이 구축된 후 비용과 시간이 많이 소요되는 재설계가 필요한 가능성을 낮추려고 애썼습니다.
역사적으로 낙하 구조는 이전 프로젝트를 적용하여 설계된 후 축소 모델로 구축되었으며, 테스트만으로도 6개월 이상이 소요될 수 있습니다. 가속화된 이 프로젝트에서, 2009년 가을에 시작한 AECOM의 초기 과제는 두가지 표준 개념 중에서 하나를 선택하는 것이었습니다. 포장-파운드 스타일과 접선 vortex버전, 둘 다 시속 35마일의 폭풍이 몰아치는 물 속에서 속도를 늦추고 통제하기 위해서 직접 계산 및 FLOW-3D에서 결정한 일반 구조 직경 및 구성 요소 크기를 사용한 초기 CFD분석으로, AECOM은 시공 가능성 및 비용 고려 사항을 평가하는 데 사용했습니다. CSO 008의 현장 요구 사항과 비용 효율성을 고려할 때, 시 당국과 AECOM은 접선 소용돌이 낙하 구조를 선택했습니다. 이 설계의 핵심 요소는 흐름을 먼저 환상적인 제트로 유도한 다음, vortex 유도 나선형 흐름을 생성하는 테이퍼(확대) 접근 채널에 의해 공급되는 수직 튜브(드롭 샤프트)입니다. 이 통제 된 하강은 속도가 느려지고 하루 3 억 갤런 (mgd) 이상에 이르는 흐름을 안전하게 처리합니다. 스토리지 터널의 파괴적인 난류를 방지하는 것이 핵심 목표이므로 드롭 샤프트 흐름의 사전 차단이 설계의 핵심입니다.
구조 자체는 6 개의 주요 부분으로 구성됩니다. 1) 접근 채널 (기존의 하수 터널에서 나온 것), 2) 수평 흐름을 넓히고 수직 드롭 샤프트로 수평 흐름을 전달하는 직사각형 전이 테이퍼 채널, 3) 드롭 샤프트 자체 4) 탈 기실 (유량을 수평 방향으로 방향을 바꾸고 공기 유입을 감소시키는), 5) 수직 공기 배출구를 통해 낙하에서 유입 된 공기를 제거하고 적하 유체의 공기 코어가 열려 있고 6) 탈기 챔버와 저장 터널 챔버를 연결하는 파이프 (adit) (그림 2).
Fig. 2. CAD diagram of proposed Indianapolis DRTC combined sewage overflow (CSO) vertical drop structure, showing approach channel, taper channel and vortex dropshaft. Using FLOW-3D CFD analysis software, AECOM simulated the flow behavior, gaining confidence in the system performance prior to physical model testing. (Image courtesy AECOM)
Design and CFD Analysis
설계 및 CFD분석
일반적인 소용돌이 설계는 널리 받아들여지고 있지만, 각 낙하 구조는 최적의 접선 흐름 특성을 보장하기 위해 인디애나 폴리스의 위상에 맞는 적절한 크기를 가져야 했습니다. 특히, 가능한 설계에 대한 AECOM의 계획은 세가지 목표를 가지고 있었습니다. 결합된 접근법과 테이퍼 채널을 짧은 길이로 제한하는 현장, 고유의 제약이 있었는지를 결정합니다. 허용 가능하지만 접근 방식에서 과도한 난류 조건이 발생하지 않았습니다. 테이퍼 채널에 안정적인 흐름 조건이 존재하는지 확인하고 다양한 흐름 조건에서 흐름 안정성을 평가했고, 논리적 기준점은 밀워키 인라인 스토리지 프로젝트라고 불리는 잘 알려지고 문서화된 시스템이었습니다.
Edison은 DRTC 프로젝트 규모에 맞춰 H-4로 지정된 Milwaukee 드롭 구조 설계를 기반으로 초기 설계를 기반으로했습니다. 166 피트의 기본 낙하 길이를 포함하고 체적 유량, 벽, 대칭 및 기타 초기 매개 변수를 지정하는 FLOW-3D 분석을 설정합니다. 그는 우리가 CFD를 통해 발견한 것은 밀워키에서 이 디자인을 사용하면 우리의 어플리케이션에 잘 맞지 않는다는 것이라고 말합니다. FLOW-3D는 이것을 보여 주고 있었기 때문에 CFD를 사용하여 변형을 시도하고 우리의 수정된 디자인을 고안했습니다. 더 넓은 접근 경로, 더 넓은 테이퍼 및/또는 더 깊은 테이퍼 깊이를 사용한 수정은 에디슨은 FLOW-3D에서 각 변동 사항을 설정하는 것이 매우 빠르다고 말합니다. (그림 3,4,5). 개선의 진전은 고무적이었습니다. 시뮬레이션 결과의 높은 수준은 심지어 절삭(침식)을 개선하기 위해 드롭 축의 바닥에 의문스러운 플레이트가 수직 흐름이 수평으로 전환되는 난류 분리 및 감소가되도록 기능을 추가하도록 설득했습니다.
Figs. 3, 4 and 5. Tangential drop structure flow simulated with FLOW-3D. Structure dimensions were optimized through multiple design iterations. (Image courtesy AECOM)
9번째 설계 변동에 대한 FLOW-3D 출력 동작인 V9는 접근 섹션을 확장했으며, 모든 흐름 볼륨 레벨에서 300mg/d까지 양호한 흐름 안정성을 보였으며 유압식 점프는 없었습니다. 그리고 양호한 Froude numners(유체 움직임에 미치는 중력의 영향을 나타내기 위해 사용되는 치수 없는 수량), 2010년 2월부터 AECOM이 물리적 시험과 검증을 위해 선택하였습니다(그림 6). 그 계획은 아이오와 연구소의 시험 결과에 기초하여 CFD와 최적화를 추가하는 것이였습니다.
Fig. 6. Scale model (1:10) of vertical drop structure, tested at University of Iowa IIHR Hydroscience & Engineering facility. (Image courtesy AECOM)
에디슨은 V9에서 결정된 치수 매개 변수에 대해 그 디자인을 아이오와 주에 가져가서 CFD를 이용해 만들었는데 완벽하게 작동했습니다. (II.)직원들은 실제로 무언가를 설치한 것은 이번이 처음이며, 변경하라고 말할 만한 것이 아무것도 없다고 말했습니다. 측정된 데이터는 드롭 샤프트 연결 구조 내의 수면 높이, Adit내 공기 침투의 정량, 벤트 샤프트 위로 공기 흐름을 포함했습니다. 흐름이 증가함에 따라 와류량이 증가함에 따라 축 벽에 부착되어 탈산소까지 원활하게 회전하는 모습이 포착되었습니다(그림 7).
에디슨은 후속 실험을 위해 여러번 시험장을 돌아다녔습니다. 물리적 모델이 처음부터 올바르게 작동했기 때문에 시험 프로그램을 확장할 시간이 있었습니다. “재미 있는 것은 환기구를 움직이는 것과 같이 우리가 궁금했던 것들을 탐구해서 지적으로 그것을 가지고 놀 시간이 있었다는 것입니다.” 에디슨은 예정보다 앞서 있었기 때문에 잔여 프로젝트 시간을 이용해 탈염소와 adit 내의 유압 장치를 조사할 수 있었습니다.
Fig. 7. Operation of scale-model vertical drop structure, showing test run of 300 million gallons per day (mgd). Flow vortex development shows good rotation and attachment to the shaft wall all the way down to the de-aeration chamber. No design modifications were necessary to the simulated design. (Image courtesy AECOM)
Final Results
AECOM은 2010년 7월 DRTC에 대한 전반적인 작업을 마쳤습니다. 2013년 3월부터 18구경 터널을 굴착하기 시작했고, CSO드롭 구조 3개(CFD로 설계된 나머지 2개의 구조물만 있음)는 모두 현재 공사 중입니다.
에디슨의 의견으로는, 토목 공학은 전체적으로 CFD를 채택하는 데 느린 편이었습니다. 이를 입증하기 위해 그는 인천 국제 공항을 처음 방문한 당시 접선 소용돌이 모형의 소위 “묘지”에서 본것을 기술했습니다. 그러나 그는 이들을 다시 처리해야 했다고 말했습니다. 그는 유압 설계를 위한 시뮬레이션 사용으로 판매되는 것을 권장하고 있습니다.
에디슨은 DRTC노력을 요약하면서 “정말 재미 있었습니다. 물리적 모델링이 필요한 위치에 대해 더 자세히 알아보았고, 그렇다면 어떤 경우에는 순수한 RAID기반 설계를 수행할 수 있습니다. 많은 DRTC작업들이 그것의 증거입니다. 물리적 모델은 실제로 필요하지 않았지만 검증을 통해 위험을 줄일 수 있었습니다. 프로젝트에서 이 두가지를 모두 수행할 수 있었다는 것은 믿을 수 없는 일입니다.”라고 말했습니다.
This article first appeared in WaterWorld Magazine.
하수도 시스템은 액션영화의 도피 루트로 사용되지 않는 한 흥미롭지 않을 것입니다. 폭우로 인해 이산화탄소 수치가 올라갈 때까지 여러분은 그것에 대해 생각조차 하지 않을 것입니다. 불행하게도, 770개 이상의 오래 된 미국 도시들 아래에 있는 하수구 시스템은 심한 폭풍으로 오염 문제를 일으킵니다. 이러한 구형 설계는 하수 및 폭풍 유실을 위한 비용 효율적인 단일 스타일 파이프를 사용했으며 연결된 파이프로 강 및 호수에 하수를 내보냅니다(CSO).
1994년 미국 환경보호청(EPA)은 주로 북동부 및 그레이트 레이크 지역의 관련 지방 자치 단체들에게 CSO관련 문제를 줄이거나 제거하도록 하는 정책을 발표했습니다. (2000년 “Clean Water Act”의 일부로 법률화된 정책). 인디애나 폴리스(Indianapolis)는 가벼운 비 폭풍으로 인해 하수 오물의 백업 및 범람이 발생할 수 있는 도시 중 하나였으므로, 주요 건설 조건에서 2025년까지 문제를 해결하는 것이 필요하였습니다.
인디애나 폴리스는 국제 디자인 회사인 AECOM에 Citizens Energy Group이 건설하고 있는 3개의 깊은 암석 저장 터널 중 첫 번째를 설계할 것을 요청했습니다. 총 25마일인 이 시스템은 대규모 지하 펌프장과 기존의 하수구에서 CSO를 수직으로 떨어뜨리는 연결 구조물을 포함합니다. 첫 번째 터널의 경우, 강우가 가라 앉은 후에 3 개의 커다란 강하 구조물이 CSO를 저장 터널로 전환하여 후속 처리를 수행했습니다.
프로젝트를 해결하기 위해 AECOM은 여러 가능한 낙하 구조물 설계의 동작을 시뮬레이션하기 위해 FLOW-3D를 선택하여, 구축 및 평가 예산이 책정 된 물리적 모델에 대한 재 작업의 필요성을 최소화했습니다. 테스트 결과는 예측 값과 일치하였으므로 재설계가 필요하지 않았습니다. 또한, 이제 AECOM은 유압 설계작업의 첫 번째 단계를 일반적으로 CFD시뮬레이션을 사용합니다.
Large Scale Project on a Tight Delivery Schedule
촉박한 납품 일정에 따른 대규모 프로젝트
20세기에 건설된 하수 처리장은 주거용, 상업용, 환경유출물의 유출로 무엇을 해야 할 것인지에 대한 새로운 인식을 가져다 주었습니다. CSO 방전은 정상적으로 운영되는 동안 처리시설로 직접 이동되며 모든 과정이 양호하게 운영됩니다. 불행하게도, 대규모 폭풍이 발생하는 동안, 발전소들의 초과 용량문제를 피하기 위해 인근 수역으로 과도한 유량을 방출합니다. 이들 배출은 기름과 살충제, 야생동물 배설물에 이르기까지 다양한 오염 물질을 포함합니다. 고무적인 성공의 신호로, 1990 년대에 착공된 새로운 CSO 분리, 저장 및 처리 시설로 오염의 영향에 대해 67 %의 개선을 이루었지만, 여전히 많은 연구가 이루어져야 합니다. 인디애나 폴리스의 경우, 인디애나 폴리스시 공공사업부가 CSO 장기 통제계획을 준비한 2008년에 그러한 노력이 시작되었습니다. 정상적인 처리 공장에서 처리 할 수 있을 때까지 오버플로우가 발생하는 “저장 및 운송”접근법의 핵심은 인디애나 폴리스 터널 저장 시스템 또는 인디애나라고 합니다.
이 시스템의 첫번째 단계는 딥 록 터널 커넥터(DRTC)라고 불리는 1억 8천만달러 가치의 프로젝트입니다. DRTC는 길이 7마일의 18피트 직경의 지하 터널로, 기존의 인디애나 폴리스의 3개의 서버 대 계층 유출 연결의 흐름 경로를 다시 만들 것입니다(그림 1). 목표는 과잉 강우 유출을 기존 하수구와 새 터널 사이의 낙하 구조를 통해 이들 대피소에서 거대한 터널로 안전하게 재배치하고, 폭풍 후 처리를 위해 처리장으로 펌핑 될 수있을 때까지 유지합니다.
Fig. 1. City of Indianapolis Deep Rock Tunnel Connector (DRTC), a “storage and transport” concept being built to handle combined sewage overflow (CSO) during heavy storms. Three vertical drop structures will capture this flow and divert it downwards to 18-foot-diameter storage tunnels running more than 250 feet underground; the tunnels store the CSO until sewage treatment plant capacity becomes available. (Image courtesy Citizens Energy Group)
평균적으로 지표면 아래 250피트 깊이에서, DRTC는 건설과 궁극적인 운영 동안 위의 주변 지역에 대한 혼란을 최소화하도록 설계되었습니다. 그러나 이 프로젝트의 규모와 복잡성은 AECOM의 과제에 긴급성을 더했습니다. 세 장소 각각에 대한 가능한 낙하 구조 설계와 평가, 구조물 설계의 60%를 7개월 이내에 마무리 지었습니다.
이러한 구조물의 목적은 표준 도시 하수 시스템에서 깊은 저장 터널로 하수 흐름을 전달하는 동시에, 효율적 손실( 느린 속도 또는 백업)과 장기적인 도심을 방지하는 것입니다. 각 섹션의 크기와 모양이 유입 흐름의 볼륨 및 속도와 세심하게 일치하지 않을 경우 발생할 수 있는 구조적 손상입니다. AECOM의 수석 기술 전문가인 라이언 에디슨 컨설턴트는 계약의 스케줄링 요구 사항이 유효성 검사를 위해서는, 단 하나의 모델에만 물리적 건물과 테스트 활동을 제한할 것이라는 것을 알게되었습니다. 다른 주요 건설 프로젝트에 15년간 FLOW-3D 시뮬레이션 소프트웨어를 사용해 왔기 때문에, 난류, 과전압 및 에너지 낭비를 예측하는 능력은 충분하지 않고 디자인 프로젝트에 적합하다고 자신했습니다. 또한 여러 검증(what-if) 시나리오를 실행하기 위한 소프트웨어 옵션을 통해 설계 세부 사항을 다시 실행해야 하는 위험을 최소화할 수 있었습니다. 변경 사항이 적용될 경우 상당한 이점은 여러개의 병렬 시공 트랙이 있는 프로젝트에 있습니다. 시간 제약에도 불구하고, 에디슨은 특히 이 도전에 만족했습니다. 왜냐하면 “CFD로 드롭 구조 설계를 만들고 물리학에서 이것들은 너무 큰 구조이기 때문입니다.”라고 그는 말합니다. 그것들은 CFD는 실제로 사용되지 않는데 보통 물리적 모델이나 손으로 계산하는 것으로 이루어집니다.
DRTC 프로젝트를 위해서, 그는 먼저 시뮬레이션된 작동 조건에 대해서 컴퓨터 설계를 테스트할 것입니다. 에디슨은 3차원의 일시적이고 격동적인 흐름 조건을 모델링 할 수 있는 소프트웨어 패키지인 FLOW-3D를 사용했습니다. 각 설계에 대한 계산 메쉬를 변경하지 않고도 여러 설계 지오 메트리를 모델링 할 수 있는 기능이였습니다. 시뮬레이션 데이터로 무장한 에디슨은 그 결과를 아이오와 대학교 II. 시설에서 시험한 1:10 크기의 물리적 모델의 작동 데이터와 비교하였습니다. (후자는 원래 아이오와 유압 연구소라고 불렸지만, 지금은 그룹의 다양한 범위를 반영하여 IIHR-Hydroscience & Engineering으로 알려져 있습니다.)
Zeroing in on the Drop-Structure Challenge
드롭 구조 과제에서 영점 조정
가장 제한적인 DRTC 사이트의 지오 메트리는 CSO 008로 지정된 레귤레이터에서 발생합니다. 기존 CSO 레귤레이터(기울기 약 75피트 아래)를 새 18피트 직경의 수집 터널과 연결하려면, 이 위치에서 150피트 이상의 수직 방향 주행이 필요합니다. 각 낙하 구조에 7백만달러 이상이 소요되는 경우, 프로젝트 관리자들은 물리적 모델이 구축된 후 비용과 시간이 많이 소요되는 재설계가 필요한 가능성을 낮추려고 애썼습니다.
역사적으로 낙하 구조는 이전 프로젝트를 적용하여 설계된 후 축소 모델로 구축되었으며, 테스트만으로도 6개월 이상이 소요될 수 있습니다. 가속화된 이 프로젝트에서, 2009년 가을에 시작한 AECOM의 초기 과제는 두가지 표준 개념 중에서 하나를 선택하는 것이었습니다. 포장-파운드 스타일과 접선 vortex버전, 둘 다 시속 35마일의 폭풍이 몰아치는 물 속에서 속도를 늦추고 통제하기 위해서 직접 계산 및 FLOW-3D에서 결정한 일반 구조 직경 및 구성 요소 크기를 사용한 초기 CFD분석으로, AECOM은 시공 가능성 및 비용 고려 사항을 평가하는 데 사용했습니다. CSO 008의 현장 요구 사항과 비용 효율성을 고려할 때, 시 당국과 AECOM은 접선 소용돌이 낙하 구조를 선택했습니다. 이 설계의 핵심 요소는 흐름을 먼저 환상적인 제트로 유도한 다음, vortex 유도 나선형 흐름을 생성하는 테이퍼(확대) 접근 채널에 의해 공급되는 수직 튜브(드롭 샤프트)입니다. 이 통제 된 하강은 속도가 느려지고 하루 3 억 갤런 (mgd) 이상에 이르는 흐름을 안전하게 처리합니다. 스토리지 터널의 파괴적인 난류를 방지하는 것이 핵심 목표이므로 드롭 샤프트 흐름의 사전 차단이 설계의 핵심입니다.
구조 자체는 6 개의 주요 부분으로 구성됩니다. 1) 접근 채널 (기존의 하수 터널에서 나온 것), 2) 수평 흐름을 넓히고 수직 드롭 샤프트로 수평 흐름을 전달하는 직사각형 전이 테이퍼 채널, 3) 드롭 샤프트 자체 4) 탈 기실 (유량을 수평 방향으로 방향을 바꾸고 공기 유입을 감소시키는), 5) 수직 공기 배출구를 통해 낙하에서 유입 된 공기를 제거하고 적하 유체의 공기 코어가 열려 있고 6) 탈기 챔버와 저장 터널 챔버를 연결하는 파이프 (adit) (그림 2).
Fig. 2. CAD diagram of proposed Indianapolis DRTC combined sewage overflow (CSO) vertical drop structure, showing approach channel, taper channel and vortex dropshaft. Using FLOW-3D CFD analysis software, AECOM simulated the flow behavior, gaining confidence in the system performance prior to physical model testing. (Image courtesy AECOM) Prediction of Shrinkage Defects During Investment Casting Process
This article was contributed by Dr. S. Savithri, Senior Principal Scientist at CSIR-NIIST
인베스트먼트 주조공정은 가장 오래된 주조 공정 중 하나로 기원전 4000년 이후에 보편화되었습니다. 이 과정은 용해된 금속을 소모품패턴으로 생성된 세라믹 쉘에 주입하는 과정을 수반합니다. 일찍이 그것은 금, 은, 구리와 청동 합금으로 장신구와 우상을 만드는데 사용되었습니다.
인베스트먼트 주조공정은 1897년 아이오와 주, 위원회 블러프스의 Barabas Frederick Philbrook이 묘사한 대로 치과의사들이 왕관과 인레이를 만들기 위해 그것을 사용하기 시작한 19세기 말 현대 산업공정으로 사용되기 시작했다. 1940년대에는 제2차 세계대전 당시 기존 방법으로는 형성될 수 없거나 지나치게 많은 가공이 필요한 특수 합금의 정밀 순모형 제조 기술에 대한 수요로 인해 투자 주조 공정이 증가하였다.
오늘날 투자 주조 공정은 표면 마감 및 치수 정확도가 우수하여 거의 순 형태에 가까운 철, 비철 및 초합금의 소형 산업용 부품을 생산하는데 주로 사용됩니다.
인베스트먼트 주조 공정은 다음 네 가지 주요 단계로 구성됩니다.
왁스 패턴 생성 후, 패턴 클러스터를 만들기 위해 게이트 시스템으로 청소 및 조립합니다.
나무는 세라믹 쉘을 얻기 위해 미세 모래와 Course한 모래 입자의 슬러리로 번갈아 코팅됩니다.
용기는 건조되고, 왁스를 녹이기 위해 가열되며, 강도를 높이고 주입 준비합니다.
마침내 주조 합금이 용해되어 예열된 쉘에 주입됩니다. 응고 후에 쉘이 파손되어 주조 부품을 얻습니다.
Figure 1. Solid model of the casting geometry
인베스트먼트 주조 공정에서 얻은 부품은 많은 중요한 용도에 사용되므로 내부적인 결함이 없어야 합니다. 투자 주조 공정에서 발생하는 주요 결함은 세라믹 포함, 균열, 변형, 플래시, 주탕불량, 수축, 슬래그 포함, 탕경계등입니다. 얻은 주조물의 품질을 예측하려면 금속-몰드 열 전달계수, 주입 온도 등 다양한 주조 공정 매개 변수의 영향을 연구해야 합니다. 즉, 쉘 두께 및 쉘 열 전달계수가 그것입니다. 현대 컴퓨터 시스템 및 시뮬레이션 소프트웨어의 출현과 함께 금형 충진 및 응고 시뮬레이션은 주조공장에서 결함을 예측하고 설계를 최적화하는데 점점 더 많이 사용되고 있습니다.
이 연구의 주요 목적은 투자 주조 공정에서 주요 요소인 복사 열 전달과 인베스트먼트 주조공정에 고유한 쉘 금형이 FLOW-3D에서 효과적으로 구현될 수 있는지를 조사하는 것입니다. FLOW-3D를 사용하여 간단한 형상을 위한 인베스트먼트 주조공정의 주입 및 응고 시뮬레이션을 수행함으로써 두 구성요소의 서로 다른 효과를 조사합니다. 다양한 위치에서 얻은 온도의 수치는 문헌 [1]에보고 된 실험 결과로 검증됩니다. 복사 열 전달계수, 쉘 몰드 두께, 탕구 및 게이트의 위치에 대한 영향도 조사했습니다.
Figure 2. Shell mold
Methodology
현재 연구에서 사용된 계산 형상은 그림 1에 나와 있습니다. 쉘 몰드는 다음 단계를 사용하여 작성되었습니다.
구성 요소 1로 형상을 FLOW-3D로 가져오고 지정된 셀 크기로 가져온 형상을 중심으로 메쉬 블록을 작성합니다.
“보완”유형의 component1의 첫 번째 하위 구성 요소를 만들어 하위 구성 요소 외부의 모든 항목을 메쉬의 범위까지 확고하게 만듭니다.
솔리드 데이터베이스에서 이 솔리드 블록의 금형 재질 특성을 정의하십시오.
솔리드 특성 GUI의 구성 요소 특성에서 “열 침투 깊이”를 정의하는 옵션이 있습니다. 여기서 쉘 두께 값을 정의 할 수 있습니다.
이제 전처리기를 실행하십시오.
분석 탭> 3D 탭으로 이동 한 다음 이전 단계에서 생성 한 prpgrf 파일을 엽니다. ‘Iso-surface’와 ‘color variable’에서 “열 활성화 구성 요소 볼륨”을 선택하고 “렌더링”을 선택하십시오.
Display에 이제 형상의 셸 부분 만 표시됩니다.
개체 목록 (창의 왼쪽 하단)에서 “구성 요소 1″을 선택하고 “구성 요소 1″을 마우스 오른쪽 단추로 클릭 한 다음 “stl로 내보내기”를 선택하여 이 곡면을 STL 파일로 저장하십시오.
Figure 3. The view of the two mesh blocks for the creation of a void with discretization
쉘 몰드 용 STL 파일을 만든 후 파일을 구성 요소 1로 새 시뮬레이션으로 가져오고 이전에 작성한 주조 형상을 하위 구성 요소로 가져오고 유형을 ‘hole’으로 선택합니다. 쉘 몰드와 함께 주조 형상이 그림 2에 나와 있습니다. 이것은 우리의 계산 영역으로 사용됩니다. 다음은 계산 영역을 cubical/rectangular셀로 분할하기 위한 메쉬를 만드는 것입니다. 메쉬 블록을 작성하여 FLOW-3D에서 메쉬를 생성합니다. 현재의 작업을 위해 우리는 2.5mm의 고정된 셀 크기가 선택된 그림 3에 표시된 균일한 메쉬 옵션을 선택했습니다. 입력 위치 주변에 메시 블록 2가 사용되는 현재 시뮬레이션을 위해 메시 블록 2개가 생성되었습니다. 쉘과 주변 공기 사이의 30°C에서의 열 전달을 고려하여 쉘 주위에 보이드 영역이 정의됩니다. 이 영역은 ‘열 전달 유형 1’이 있는 보이드 영역으로 선택되며 셸과 주변 공기 사이에 열 전달 계수 값이 지정됩니다. 열 전달 유형 1은 방사선을 포함한 종합 열 전달 계수가 됩니다.
쉘 주형에 선택된 재료는 zircon이며 열 특성은 Sabau and Vishwanathan에 의해 수행된 실험에서 얻을 수 있습니다[2]. 표 1은 연구에 사용된 재료에 대해 지정된 값을 보여 줍니다.
MATERIAL
PROPERTY
VALUE
UNIT
Fluid –AluminiumA356
alloy
Density
2437
kg/m³
Thermal conductivity
116.8
W/(mK)
Specific heat
1074
J/(kgK)
Latent heat
433.22
kJ/m³
Liquidus temperature
608
°C
Solidus temperature
552.4
°C
Zircon Mold
Thermal conductivity
1.09
W/(mK)
Specific heat* Density
1.63E+06
J/( m³K)
Initial and boundary conditions used are show in Table 2.
Mold temperature
430°C
Melt pouring temperature
680°C
Filling time
7 s
Interface heat transfer coefficient
850 W/m2K
Heat transfer coefficient between ambient and mold (radiation effect)
30 -100 W/m2K
Table 2. Initial and boundary conditions used for the simulation
탕구저에 들어가는 용융물의 초기 속도와 온도는 메시 블록 2의 상단 경계에서 속도 경계 조건으로 주어집니다. 기본적으로 다른 모든 경계는 대칭 유형으로 설정됩니다.
Results & Discussion
Validation with reported experimental results
충전 및 응고 동안 냉각 곡선을 얻기 위한 실험에서 Sabuet.al[1]에 의해 선택된 네 개의 위치가 검증 목적으로 사용되었습니다. 그들은 C1, C2, S11, S12및 S21로 언급됩니다. C1과 C2지점은 주물의 플레이트의 중심에 있으며 S11, S12및 S21은 모두 쉘에 위치합니다. 이러한 위치에서의 온도 변화는 그림 4와 같습니다.
온도 프로파일의 수치 및 실험결과의 차이가 허용한계 안에 있음을 알 수 있습니다. 프로브 포인트 C1과 C2의 경우, 수치와 실험 결과 사이의 차이는 응고 중에 5%, 응고 후 냉각 시 12% 이내입니다. 쉘의 점에 대한 수치 결과는 실험 결과보다 약 5% 높습니다. 이는 쉘 재료에 열 물리학적 특성을 할당할 때 발생하는 가정과 쉘 열 전달 계수의 값 때문일 수 있습니다.
Fill sequence & solidification pattern for two different sprue locations
두 가지 다른 스프 루 위치의 채우기 순서 및 응고 패턴
2 개의 상이한 탕구 위치에 주물충전 순서는5a 및5b에 나와 있습니다. 최종 탕구가 더 많은 스플라인을 생성하므로 결함으로 이어질 수 있습니다. 탕구가 중간에 놓여지면 흐름은 보다 균일 해지고 두 주조 단면에서 비슷한 온도 분포를 보입니다. 50 % 응고 후의 온도 프로파일의 2D 도면은 두 경우 모두 그림 5c 및 5d에 나와 있습니다. 수축 위치에서 볼 때 두 탕구 위치가 결함을 일으키는 것은 분명합니다.
Figure 5a. Fill sequence at different time intervals when the sprue is located at one end
Figure 5b. Fill sequence at different time intervals when the sprue is located in the middle
Figure 5c. 2D temperature profile after 50% solidification when the sprue is located at one end
Figure 5d. 2D temperature profile after 50% solidification when the sprue is located in the middle
Effect of shell thickness
인베스트먼트 주조에 대한 쉘 두께의 효과를 연구하기 위해 두께가 7.2, 10, 15 및 20 mm인 주물을 선정하였습니다. 그림 6a 및 6b는 주조품의 특정 위치에서 냉각 곡선을 나타내며, 이는 C1으로 나타내고 쉘 몰드 내의 특정 위치에 있으며, 응고 중에 S11로 나타납니다. 세라믹 쉘의 두께가 7.2 mm에서 15 mm로 증가하면 냉각 속도가 감소하여 응고 시간이 길어지는 것을 볼 수 있습니다.
Effect of shell heat transfer coefficient
셸 열 전달 계수는 열이 셸 금형의 외부 벽에서 방사선을 통해 주변 공기로 열을 방출하는 속도를 나타냅니다. 이 효과를 조사하기 위해 열 전달 계수의 값을 20에서 80W/m2K까지 다양하게 했습니다. 7a 및 7b로부터, h의 변화는 주조 재료 및 쉘의 냉각 속도에 중요한 영향을 미친다는 것을 알 수 있습니다. 열 전달 계수가 20에서 80W/m2K로 증가하면 C1에서의 응고 시간이 812 초에서 334 초 (약 44 %)로 감소되었음을 알 수 있습니다. 따라서, h의 값을 변화시키는 것은 주물의 미세 구조에 영향을 미칩니다.
Figure 6a. Temperature profile at location C1 (casting) for the casting geometry where the sprue is located at one end for various shell thickness values
F
Figure 6b. Temperature profile at location S11 (shell) for the casting geometry where the sprue is located at one end for various shell thickness values
Figure 7a. Temperature profile at location C1 (casting) for the casting geometry where the sprue is located at one end for various heat transfer coefficient values between the shell mold & ambient
Figure 7b. Temperature profile at location S11 (shell) for the casting geometry where the sprue is located at one end for various heat transfer coefficient values between the shell mold & ambient
Conclusions
인베스트먼트 주조 공정의 몰드 충진 및 응고 시뮬레이션은 FLOW-3D를 사용하여 수행되었습니다. 주조 공정에 대한 주조 매개변수의 영향을 연구하기 위해 파라메트릭 연구가 수행되었습니다. 본 연구에서 다음과 같은 결론을 도출 할 수 있습니다.
FLOW-3D는 멀티 캐비티 몰드의 주입 및 응고 모델링이 가능합니다. 프로브 위치의 예측 온도 프로파일은 실험 데이터의 허용오차 이내였다.
쉘 두께의 경우, 두 경우 모두 셸의 임계 두께가 있으며, 그 이상으로 열 전달 특성이 역행하는 것으로 확인되었습니다. 셸 두께가 증가함에 따라 응고 시간이 임계 두께까지 증가하여 감소하기 시작했습니다. 원래 형상의 경우 임계 두께는 15~20mm인 반면 수정된 형상의 경우 10mm와 15mm 사이에 있다.
쉘과 대기 사이의 열 전달 계수 h는 열 전달 특성에 가장 큰 영향을 미치는 것으로 나타났습니다. h가 20에서 80W/m2K로 4 배 증가할 때 탕구의 중심에서 응고 시간이 40 % 이상 감소했습니다.
References
Sabau, A.S., Numerical Simulation of the Investment Casting Process, Transactions of the American Foundry Society, vol. 113, Paper No. 05-160, 2005.
Sabau, A.S., and Viswanathan, S., Thermophysical Properties of Zircon and Fused Silica-based Shells used in the Investment Casting Process, Transactions of the American Foundry Society, vol. 112, Paper No. 04-081, 2004.
Design and CFD Analysis
설계 및 CFD분석
일반적인 소용돌이 설계는 널리 받아들여지고 있지만, 각 낙하 구조는 최적의 접선 흐름 특성을 보장하기 위해 인디애나 폴리스의 위상에 맞는 적절한 크기를 가져야 했습니다. 특히, 가능한 설계에 대한 AECOM의 계획은 세가지 목표를 가지고 있었습니다. 결합된 접근법과 테이퍼 채널을 짧은 길이로 제한하는 현장, 고유의 제약이 있었는지를 결정합니다. 허용 가능하지만 접근 방식에서 과도한 난류 조건이 발생하지 않았습니다. 테이퍼 채널에 안정적인 흐름 조건이 존재하는지 확인하고 다양한 흐름 조건에서 흐름 안정성을 평가했고, 논리적 기준점은 밀워키 인라인 스토리지 프로젝트라고 불리는 잘 알려지고 문서화된 시스템이었습니다.
Edison은 DRTC 프로젝트 규모에 맞춰 H-4로 지정된 Milwaukee 드롭 구조 설계를 기반으로 초기 설계를 기반으로했습니다. 166 피트의 기본 낙하 길이를 포함하고 체적 유량, 벽, 대칭 및 기타 초기 매개 변수를 지정하는 FLOW-3D 분석을 설정합니다. 그는 우리가 CFD를 통해 발견한 것은 밀워키에서 이 디자인을 사용하면 우리의 어플리케이션에 잘 맞지 않는다는 것이라고 말합니다. FLOW-3D는 이것을 보여 주고 있었기 때문에 CFD를 사용하여 변형을 시도하고 우리의 수정된 디자인을 고안했습니다. 더 넓은 접근 경로, 더 넓은 테이퍼 및/또는 더 깊은 테이퍼 깊이를 사용한 수정은 에디슨은 FLOW-3D에서 각 변동 사항을 설정하는 것이 매우 빠르다고 말합니다. (그림 3,4,5). 개선의 진전은 고무적이었습니다. 시뮬레이션 결과의 높은 수준은 심지어 절삭(침식)을 개선하기 위해 드롭 축의 바닥에 의문스러운 플레이트가 수직 흐름이 수평으로 전환되는 난류 분리 및 감소가되도록 기능을 추가하도록 설득했습니다.
Figs. 3, 4 and 5. Tangential drop structure flow simulated with FLOW-3D. Structure dimensions were optimized through multiple design iterations. (Image courtesy AECOM)
9번째 설계 변동에 대한 FLOW-3D 출력 동작인 V9는 접근 섹션을 확장했으며, 모든 흐름 볼륨 레벨에서 300mg/d까지 양호한 흐름 안정성을 보였으며 유압식 점프는 없었습니다. 그리고 양호한 Froude numners(유체 움직임에 미치는 중력의 영향을 나타내기 위해 사용되는 치수 없는 수량), 2010년 2월부터 AECOM이 물리적 시험과 검증을 위해 선택하였습니다(그림 6). 그 계획은 아이오와 연구소의 시험 결과에 기초하여 CFD와 최적화를 추가하는 것이였습니다.
Fig. 6. Scale model (1:10) of vertical drop structure, tested at University of Iowa IIHR Hydroscience & Engineering facility. (Image courtesy AECOM)
에디슨은 V9에서 결정된 치수 매개 변수에 대해 그 디자인을 아이오와 주에 가져가서 CFD를 이용해 만들었는데 완벽하게 작동했습니다. (II.)직원들은 실제로 무언가를 설치한 것은 이번이 처음이며, 변경하라고 말할 만한 것이 아무것도 없다고 말했습니다. 측정된 데이터는 드롭 샤프트 연결 구조 내의 수면 높이, Adit내 공기 침투의 정량, 벤트 샤프트 위로 공기 흐름을 포함했습니다. 흐름이 증가함에 따라 와류량이 증가함에 따라 축 벽에 부착되어 탈산소까지 원활하게 회전하는 모습이 포착되었습니다(그림 7).
에디슨은 후속 실험을 위해 여러번 시험장을 돌아다녔습니다. 물리적 모델이 처음부터 올바르게 작동했기 때문에 시험 프로그램을 확장할 시간이 있었습니다. “재미 있는 것은 환기구를 움직이는 것과 같이 우리가 궁금했던 것들을 탐구해서 지적으로 그것을 가지고 놀 시간이 있었다는 것입니다.” 에디슨은 예정보다 앞서 있었기 때문에 잔여 프로젝트 시간을 이용해 탈염소와 adit 내의 유압 장치를 조사할 수 있었습니다.
Fig. 7. Operation of scale-model vertical drop structure, showing test run of 300 million gallons per day (mgd). Flow vortex development shows good rotation and attachment to the shaft wall all the way down to the de-aeration chamber. No design modifications were necessary to the simulated design. (Image courtesy AECOM)
Final Results
AECOM은 2010년 7월 DRTC에 대한 전반적인 작업을 마쳤습니다. 2013년 3월부터 18구경 터널을 굴착하기 시작했고, CSO드롭 구조 3개(CFD로 설계된 나머지 2개의 구조물만 있음)는 모두 현재 공사 중입니다.
에디슨의 의견으로는, 토목 공학은 전체적으로 CFD를 채택하는 데 느린 편이었습니다. 이를 입증하기 위해 그는 인천 국제 공항을 처음 방문한 당시 접선 소용돌이 모형의 소위 “묘지”에서 본것을 기술했습니다. 그러나 그는 이들을 다시 처리해야 했다고 말했습니다. 그는 유압 설계를 위한 시뮬레이션 사용으로 판매되는 것을 권장하고 있습니다.
에디슨은 DRTC노력을 요약하면서 “정말 재미 있었습니다. 물리적 모델링이 필요한 위치에 대해 더 자세히 알아보았고, 그렇다면 어떤 경우에는 순수한 RAID기반 설계를 수행할 수 있습니다. 많은 DRTC작업들이 그것의 증거입니다. 물리적 모델은 실제로 필요하지 않았지만 검증을 통해 위험을 줄일 수 있었습니다. 프로젝트에서 이 두가지를 모두 수행할 수 있었다는 것은 믿을 수 없는 일입니다.”라고 말했습니다.
This article first appeared in WaterWorld Magazine.
Integration of CFD Analysis into Die-Cast Process Design
This article was contributed by Alex Reikher, Ph.D., of Shiloh Industries
오늘날의 조직은 오래되고, 잘 구축되었으며, 빠르게 성장하는 새로운 경제로부터 점점 더 많은 압박을 받고 있습니다. 시장의 세계화는 기업들이 그들의 경쟁 우위를 유지하기 위한 방안을 찾고 있습니다. 인터넷 기술의 급격한 발전과 자유로운 정보 교환은 기업이 경쟁 우위를 유지할 수 있는 기간을 단축하는 요인들입니다. 조직이 업계에서 선두 자리를 유지할 수 있는 방법 중 하나는 혁신기술을 시장에 도입하는데 필요한 시간을 줄이는 것입니다. 다이캐스팅 공정 개발 시간 단축이라는 목표를 가지고, FLOW-3D로 모델링하는 것은 Shiloh Industries의 엔지니어링 부서의 핵심 부분이 되었습니다.
우리는 7년이 넘게 다이캐스팅 모델링 도구인 FLOW-3D를 사용하여 예측 결과의 정확성과 신뢰성을 입증할 수 있었습니다. 이러한 결과는 실제 주조 결함, 온도 분포 및 흐름 패턴과 좋은 상관 관계를 가지고 있습니다.
Shiloh Industries의 새로운 프로젝트는 게이트와 러너의 컨셉개발에서 시작되며 대략적인 slow shot profile, shot 실린더 직경, 최소 환기 영역 및 프로세스압력 요구 조건사항이 필요합니다. 유동분석은 최상의 유동패턴을 개발하고 공기 유입을 최소화하기 위해 수행됩니다. 러너설계가 완료된 후 열 분석을 실행하여 waterline 배치를 최적화합니다
Figure 1: The casting part
FLOW-3D의 매력적인 특징은 프로세스의 단계마다 별도의 분석을 실행할 수 있는 기능입니다. 이를 통해 올바른 shot 프로파일, 게이트 디자인 및 워터 라인 위치를 선택할 수 있습니다. 완전히 결합된 흐름과 열 분석은 모든 구성 요소가 잘 작동하는지 확인하기 위해 한번만 수행하면 됩니다. GMO(움직이는 물체)모델을 도입하면 저속 촬영단계에서 샷 슬리브의 최고 플런저 속도를 설정할 수 있습니다. 여기에 설명된 프로젝트에서 부품 설계는 현재 생산 버전에서 크게 변경되었습니다.
부품 형상은 그림 1에 나와 있습니다. 요구되는 주조 품질을 보장하기 위해 충전 및 응고 과정에서 어려움이 있습니다. 예를 들어, 응고 및 후속 냉각 중에 높은 내부 응력이 발생하여 바람직하지 않은 변형력이 발생할 수 있습니다.
설계 프로세스의 초기 단계에서 평가를 위해 21 개의 러너 구성이 제안되었습니다. FLOW-3D는 모든 변형을 평가하는데 사용되었습니다. 그림 2는 고려된 주 디자인 중 일부를 보여줍니다.
Figure 2: Three of the twenty-one runner systems modeled in FLOW-3D
러너 시스템의 초기 평가 기준은 유동패턴이었습니다. 설계 프로세스의 첫 번째 단계가 완료된 후, 추가 평가를 위해 그림 3에 표시된 두 가지의 러너 설계가 추가평가를 위해 채택되었습니다.
Figure 3: Runners selected for further evaluation based on the flow pattern
응고 분석은 두 번째 단계에서 평가하였습니다. 주조물뿐만 아니라 다이의 온도 분포도 분석하였습니다. 그림 4는 최종 러너 시스템 설계를 하여 응고가 끝날 때 부품의 온도 분포를 보여줍니다.
Figure 4: Different views of the final runner system chosen based on temperature distribution in the part at the end of solidification
Conclusion
7년이 넘는 기간 동안 우리는 다이 캐스팅공정 도구인 FLOW-3D를 사용하여 예측 결과의 정확성과 신뢰성을 입증할 수 있었습니다. 이러한 결과는 실제 주조 결함, 온도 분포 및 흐름 패턴과 좋은 상관 관계를 가지고 있습니다.
우리는 다이 캐스팅 공정 시뮬레이션뿐만 아니라 일반적인 CFD 모델링도 FLOW-3D를 사용하고 있습니다. 프로세스 개발 중에 설계 변경을 고객에게 권장해야 하는 경우 FLOW-3D를 사용하면 이러한 변경 사항을 신속하고 안정적으로 평가할 수 있으며 제안된 변경 사항뿐만 아니라 변경사항이 부품 성능에 미칠 영향을 고객에게 제시할 수 있습니다.
Learn more about the versatility and power of modeling metal casting processes with FLOW-3D Cast
Increasing Productivity by Reducing Ejection Times
This article was contributed by Eugene Moore of Hellebusch Tool & Die
시뮬레이션 소프트웨어는 설계자와 엔지니어가 주조 공정의 세부 사항을 이해하고 경쟁사보다 저렴한 비용으로 고품질 부품을 일관성 있게 제작할 수 있게 해주는 유용한 도구입니다. 고압 다이캐스팅에서 시뮬레이션 소프트웨어는 주조 내로 금속을 공급하고 난류로 인한 공기 유입을 방지하기 위해 샷 슬리브 팁의 타이밍을 개선하고 오버 플로우에 대한 가장 효과적인 위치를 식별하는 더 나은 게이팅 시스템을 설계하는 데 사용됩니다. 이 기사에서는 프로세스 시간을 단축하기 위해 부품을 다이에서 배출하기 전에 시간을 줄이는 방법을 살펴 보겠습니다.
비스킷은 주조 과정에서 고형화 된 마지막 장소이기 때문에 우리의 노력에 집중할 수 있는 자연스러운 곳이며 따라서 부품을 언제 꺼낼 수 있는지를 결정합니다. 따라서 비스킷의 응고 시간을 줄일 수 있다면 전반적인 공정 시간을 줄일 수 있습니다. 이를 수행하는 한 가지 방법은 유체와 접촉하는 영역의 양을 늘려 샷팁을 통해 금속에서 더 많은 열을 제거하는 것입니다. 이 경우 정확하게 적용할 수는 없지만, 아래에 표시된 정상상태 대류 방정식을 사용하면 이 접근법의 기초가 가장 쉽게 표시됩니다.
이 방정식에서 열 흐름은 대류 열 전달 계수이고, 금속 팁과 샷 팁 온도의 차이이며, 금속과 접촉하는 샷 팁의 표면적입니다. 그림 1에서 볼 수 있듯이 오늘날 시장에서 볼 수 있는 다양한 형태의 플런저 팁이 금속과 접촉하는 표면적을 증가시키도록 설계되었습니다. Figure 1: Plunger tips varying in size and surface area [1]
비스킷에서 제거된 열을 증가시키는 또 다른 방법은 비스켓에서 샷 팁과 금속 사이의 온도 차이를 조절하는 것입니다. 이는 그림 2 에서처럼 냉각 선을 팁에 추가하여 수행됩니다. 이 접근법의 단점은 피스톤 어셈블리에 상당한 복잡성을 추가한다는 것입니다. Figure 2: Cooling within plunger tip [2]
New design
이 기사에서 FLOW-3D Cast를 사용하여 새로운 플런저 팁 디자인을 분석하고 수정되지 않은 표준 원통형 팁과 비교했습니다. 그림 3에서와 같이 끝 부분에 별 모양의 컷 아웃이 있는 원통형 팁으로 구성된 수정된 팁은 수정되지 않은 샷 팁보다 20 % 더 많은 표면적을 갖습니다. 팁은 분석을 위해 물로 냉각되지 않습니다. Figure 3: Shape of the modified tip to give 20% increase in area
Analysis
각 샷 팁 디자인에 대해 충진 (샷 팁 모션 포함) 및 응고 (플로우 미포함) 시뮬레이션을 실행했습니다. 모든 다른 매개 변수는 사례간에 동일합니다. 주요 관심사는 두 가지입니다. 즉, 충전 중 흐름 패턴과 전반적인 응고 시간입니다. 샷 팁 디자인이 파동 및 공기 유입을 유발하는 경우 팁 또는 샷 슬리브 프로파일을 다시 설계해야 하기 때문에 충진 중 흐름 패턴이 중요합니다.
첫 번째 비교는 그림 4에 표시된 샷 슬리브의 흐름 패턴입니다. 이 그림은 수정 된 팁이 있거나 없는 샷 슬리브 중 유체의 이미지를 보여 주며 팁의 모양이 샷 슬리브에 영향을 주지 않는 것으로 나타났습니다. 흐름 패턴. 샷 프로파일에 거의 영향을 주지 않기 때문에 응고에 집중할 수 있습니다. Figure 4: Flow patterns in the shot sleeves from both tips.
두 번째 비교는 응고 시간입니다. 그림 5는 시간의 함수로서 팁의 평균 온도, 시간의 함수로서 금속으로부터 팁으로의 열 유속 및 추출시의 액체 금속의 온도 프로파일을 비교합니다 Figure 5: The above time plots show the average temperature in tip on the upper left hand corner and the heat flux from the metal to the tip in the upper right hand corner. The images below this show the metal temperature within the biscuit of the two castings.
그림 5에서 볼 수 있듯이, 그래프는 금속에서 더 많은 열을 추출했기 때문에 수정 된 팁의 평균 온도가 더 높음을 보여줍니다. 이것은 또한 열 유속 플롯에 표시됩니다. 그래프 아래의 이미지는 비스킷과 탄환의 경계면에서 액체금속을 보여줍니다. 데이터는 수정된 팁을 사용하여 열 제거가 12.7 % 증가한 것을 보여줍니다.
Conclusions
샷 팁 디자인은 주조 부품의 응고 시간에 눈에 띄는 영향을 미칩니다. 시뮬레이션 소프트웨어는 효과를 분석하고 이 지식을 사용하여 프로세스 매개변수를 최적화하는 방법을 제공합니다.
FLOW-3D CAST는 엔지니어가 주조 공정 설계 개념을 빠르고 효율적으로 탐구하고 최적화 할 수 있게 합니다. FLOW-3D CAST를 사용하여 위해 귀금속 주조를 포함한 다양한 주물 공정은 보석을 만들기 위해 모델링하고 최적화 할 수 있습니다. 완벽성이 중요한 보석 산업에서 FLOW-3D CAST의 매우 정확한 주입 및 응고 시뮬레이션은 표면 산화물, 공기 침투, 수축 및 기공을 포함한 중요한 보석 주조 결함을 예측합니다. Flow Science의 고도로 숙련 된 기술 엔지니어 팀은 포괄적인 교육은 물론, 전화, 전자 메일 및 금속 주조 산업에 대한 전문가의 지원을 제공합니다.
Centrifugal Casting of a Jewelry Ring Tree
FLOW-3D CAST simulation of centrifugal casting with non-inertial reference frame motion on the top and stationary motion on the bottom left.
Gravity Pour Ring Casting
Gravity casting of a jewelry tree containing 2 rings. Fluid is injected into the tree in an asymmetrical manner – the center of inflow source is offset the center of the crucible.
Single Ring Centrifugal Casting
하나의 백금 링을 FLOW-3D CAST를 사용하여 원심력으로 주조했습니다. 주조 공정에서 중요한 단계인 충진 과정은 주조 엔지니어가 주조 링의 표면 결함, 공기 유입 및 주입 온도를 자세하게 볼 수 있게 시뮬레이션 됩니다.
중력 주조는 큰 부품(일반적으로 철, 청동, 황동 또는 알루미늄)을 만드는 데 사용됩니다. 사형 주조 및 영구 금형을 포함한 대부분의 주조 공장 주조 공정은 FLOW-3D CAST를 사용하여 모델링 할 수 있습니다. 주입 프로세스는 고압 다이 캐스팅에 비해 덜하지만 과도한 공기 주입으로 인한 공기 유입으로 인해 품질이 저하될 수 있습니다. 주입하는 동안 잠재적 결함의 위치와 온도의 변화 뿐만 아니라, 용탕 표면의 움직임도 정확하게 예측됩니다. 충진이 완료된 후 용탕의 응고 및 수축을 모델링 할 수도 있습니다.
Accurate Filling Simulations
주조 공정에서 주입 작업은 결함들이 라이저로 이동하는지, 또는 부품에 갇힌 채로 남아 있는지 여부와 같은 주입 패턴 및 관련 결함을 분석하는 작업으로 이루어집니다. 시뮬레이션 분석을 사용하면 설계의 효율성을 검증하고 비용을 절감하면서 생산에 들어가기 전에 설계를 테스트할 수 있습니다. 주입의 정확성은 산화물의 결함과 갇힌 공기의 위치를 추적하는 데 중요할 뿐만 아니라, 응고 결과의 핵심입니다. 올바른 주입 패턴은 주입 마지막의 올바른 열 분포를 의미합니다. 이 열 분포는 응고 분석의 기초가 됩니다.
BMW Injector Casting Process – Innovative ingate system for gravity casting. Courtesy BMW and Flow Science Deutschland.
Solidification of Castings for Foundry Applications
편석, 열응력, 마이크로 및 매크로 기공 등 응고와 관련된 다양한 결함들이 있습니다. 정확한 응고 결과를 얻기 위한 중요한 첫번째 단계는, 정확한 주입입니다. 정확한 주입은 응고 모델링의 초기 조건인 올바른 열 프로필을 캡처하는데, FLOW-3D CAST는 주조 부품을 보다 신속하게 설계하고 폐기율을 낮출 수 있는 많은 응고 관련 결함을 감지할 수 있습니다.
FLOW-3D는 미국 뉴멕시코주(New Mexico) 로스알라모스(Los Alamos)에 있는 Flow Scicence, Inc에서 개발한 범용 전산유체역학(Computational Fluid Dynamics) 프로그램입니다. 로스알라모스 국립연구소의 수치유체역학 연구실에서 F.Harlow, B. Nichols 및 T.Hirt 등에 의해 개발된 MAC(Marker and Cell) 방법과 SOLA-VOF 방식을 기초로 하여, Hirt 박사가 1980년에 Flow Science, Inc사를 설립하여 계속 프로그램을 발전시켰으며 1985년부터 FLOW-3D를 전세계에 배포하였습니다.
유체의 3차원 거동 해석을 수행하는데 사용되는 CFD모형은 몇몇 있으나, 유동해석에 적용할 물리모델 선정은 해석의 정밀도와 밀접한 관계가 있으므로, 해석하고자 하는 대상의 유동 특성을 분석하여 신중하게 결정하여야 합니다.
FLOW-3D는 자유표면(Free Surface) 해석에 있어서 매우 정확한 해석 결과를 제공합니다. 해석방법은 자유표면을 포함한 비정상 유동 상태를 기본으로 하며, 연속방정식, 3차원 운동량 보전방정식(Navier-Stokes eq.) 및 에너지 보존방정식 등을 적용할 수 있습니다.
FLOW-3D는 유한차분법을 사용하고 있으며, 유한요소법(FEM, Finite Element Method), 경계요소법(Boundary Element Method)등을 포함하여 자유표면을 포함하는 유동장 해석(Fluid Flow Analysis)에서 공기와 액체의 경계면을 정밀하게 표현 가능합니다.
유체의 난류 해석에 대해서는 혼합길이 모형, 난류 에너지 모형, RNG(Renormalized Group Theory) k-ε 모형, k-ω 모형, LES 모형 등 6개 모형을 적용할 수 있으며, 자유표면 해석을 위하여 VOF(Volume of Fluid) 방정식을 사용하고, 격자 생성시 사용자가 가장 쉽게 만들 수 있는 직각형상격자는 형상을 더욱 정확하게 표현하기 위해 FAVOR(Fractional Area Volume Obstacle Representation) 기법을 각 방정식에 적용하고 있습니다.
FLOW-3D는 비압축성(Incompressible Fluid Flow), 압축성 유체(Compressible Fluid Flow)의 유동현상 뿐만 아니라 고체와의 열전달 현상을 해석할 수 있으며, 비정상 상태의 해석을 기본으로 합니다.
FLOW-3D v12.0은 모델 설정을 간소화하고 사용자 워크 플로우를 개선하는 GUI(그래픽 사용자 인터페이스)의 설계 및 기능에 있어 중요한 변화를 가져왔습니다. 최첨단 Immersed Boundary Method는 FLOW-3Dv12.0솔루션의 정확도를 높여 줍니다. 다른 특징적인 주요 개발에는 슬러지 안착 모델, 2-유체 2-온도 모델, 사용자가 자유 표면 흐름을 훨씬 더 빠르게 모델링 할 수 있는 Steady State Accelerator등이 있습니다.
물리 및 수치 모델
Immersed Boundary Method
힘과 에너지 손실에 대한 정확한 예측은 솔리드 바디 주변의 흐름과 관련된 많은 엔지니어링 문제를 모델링하는 데 중요합니다. FLOW-3D v12.0의 릴리스에는 이러한 문제 해결을 위해 설계된 새로운 고스트 셀 기반 Immersed Boundary Method (IBM)가 포함되어 있습니다. IBM은 내부 및 외부 흐름을 위해 벽 근처 해석을 위해 보다 정확한 솔루션을 제공하여 드래그 앤 리프트 힘의 계산을 개선합니다.
Two-field temperature for the two-fluid model
2유체 열 전달 모델은 각 유체에 대한 에너지 전달 공식을 분리하도록 확장되었습니다. 이제 각 유체에는 고유한 온도 변수가 있어 인터페이스 근처의 열 및 물질 전달 솔루션의 정확도를 향상시킵니다. 인터페이스에서의 열 전달은 시간의 표 함수가 될 수 있는 사용자 정의 열 전달 계수에 의해 제어됩니다.
슬러지 침전 모델 / Sludge settling model
중요 추가 기능인 새로운 슬러지 침전 모델은 도시 수처리 시설물 응용 분야에 사용하면 수처리 탱크 및 정화기의 고형 폐기물 역학을 모델링 할 수 있습니다. 침전 속도가 확산된 위상의 방울 크기에 대한 함수인 드리프트-플럭스 모델과 달리, 침전 속도는 슬러지 농도의 함수이며 기능적인 형태와 표 형태로 모두 입력 할 수 있습니다.
Steady-state accelerator for free surface flows
이름이 암시하듯이, 정상 상태 가속기는 안정된 상태의 솔루션에 대한 접근을 가속화합니다. 이는 작은 진폭의 중력과 모세관 현상을 감쇠하여 이루어지며 자유 표면 흐름에만 적용됩니다.
Void particles
보이드 입자가 버블 및 위상 변경 모델에 추가되었습니다. 보이드 입자는 항력과 압력 힘을 통해 유체와 상호 작용하는 작은 기포의 역할을 하는 붕괴된 보이드 영역을 나타냅니다. 주변 유체 압력에 따라 크기가 변경되고 시뮬레이션이 끝난 후 최종 위치는 공기 침투 가능성을 나타냅니다.
Sediment scour model
침전물의 정확성과 안정성을 향상시키기 위해 침전물의 운반과 침식 모델을 정밀 조사하였다. 특히, 침전물 종에 대한 질량 보존이 크게 개선되었습니다.
Outflow pressure boundary condition
고정 압력 경계 조건에는 이제 압력 및 유체 비율을 제외한 모든 유량이 해당 경계의 상류에 있는 흐름 조건을 반영하는 ‘유출’ 옵션이 포함됩니다. 유출 압력 경계 조건은 고정 압력 및 연속성 경계 조건의 혼합입니다.
Moving particle sources
시뮬레이션 중에 입자 소스는 이동할 수 있습니다. 시간에 따른 변환 및 회전 속도는 표 형식으로 정의됩니다. 입자 소스의 운동은 소스에서 방출 된 입자의 초기 속도에 추가됩니다.
Variable center of gravity
중력 및 비 관성 기준 프레임 모델에서 시간 함수로서의 무게 중심의 위치는 외부 파일의 표로 정의할 수 있습니다. 이 기능은 연료를 소모하는 로켓을 모델링하고 단계를 분리할 때 유용합니다.
공기 유입 모델
가장 간단한 부피 기반 공기 유입 모델 옵션이 기존 질량 기반 모델로 대체되었습니다. 질량 기반 모델은 부피와 달리 주변 유체 압력에 따라 부피가 변화하는 동안 흡입된 공기량이 보존되기 때문에 물리학적 모델입니다.
Tracer diffusion / 트레이서 확산
유동 표면에서 생성된 추적 물질은 분자 및 난류 확산 과정에 의해 확산될 수 있으며, 예를 들어 실제 오염 물질의 거동을 모방합니다.
모델 설정
시뮬레이션 단위
이제 온도를 포함하여 단위계 시스템을 완전히 정의해야 합니다. 표준 단위 시스템이 제공됩니다. 또한 사용자는 선택한 옵션에서 질량, 시간 및 길이 단위를 정의하여 편리하며, 사용자 정의된 단위를 사용할 수 있습니다. 사용자는 또한 압력이 게이지 단위로 정의되는지 절대 단위로 정의되는지 여부를 지정해야 합니다. 기본 시뮬레이션 단위는 Preferences(기본 설정)에서 설정할 수 있습니다. 단위를 완벽하게 정의하면 FLOW-3D는 물리적 수량에 대한 기본 값을 정의하고 범용 상수를 설정할 수 있으므로 사용자가 필요로 하는 작업량을 최소화할 수 있습니다.
Shallow water model
천수(shallow water) 모델에서 매닝의 거칠기
Manning의 거칠기 계수는 지형 표면의 전단 응력 평가를 위해 천수(shallow water) 모델에서 구현되었습니다. 표면 결함의 크기를 기반으로 기존 거칠기 모델을 보완하며이 모델과 함께 사용할 수 있습니다. 표준 거칠기와 마찬가지로 매닝 계수는 구성 요소 또는 하위 구성 요소의 속성이거나 지형 래스터 데이터 세트에서 가져올 수 있습니다.
메시 생성
하단 및 상단 경계 좌표의 정의만으로 수직 방향의 메시 설정이 단순화되었습니다.
구성 요소 변환
사용자는 이제 여러 하위 구성 요소로 구성된 구성 요소에 회전, 변환 및 스케일링 변환을 적용하여 복잡한 형상 어셈블리 설정 프로세스를 단순화 할 수 있습니다. GMO (General Moving Object) 구성 요소의 경우, 이러한 변환을 구성 요소의 대칭 축과 정렬되도록 신체에 맞는 좌표계에 적용 할 수 있습니다.
런타임시 스레드 수 변경
시뮬레이션 중에 솔버가 사용하는 스레드 수를 변경하는 기능이 런타임 옵션 대화 상자에 추가되어 사용 가능한 스레드를 추가하거나 다른 태스크에 자원이 필요한 경우 스레드 수를 줄일 수 있습니다.
프로브 제어 열원
활성 시뮬레이션 제어가 형상 구성 요소와 관련된 heat sources로 확장되었습니다. history probes로 열 방출을 제어 할 수 있습니다.
소스에서 시간에 따른 온도
질량 및 질량/모멘트 소스의 유체 온도는 이제 테이블 입력을 사용하여 시간의 함수로 정의 할 수 있습니다.
방사율 계수
공극으로의 복사 열 전달을위한 방사율 계수는 이제 사용자가 방사율과 스테판-볼츠만 상수를 지정하도록 요구하지 않고 직접 정의됩니다. 후자는 이제 단위 시스템을 기반으로 솔버에 의해 자동으로 설정됩니다.
Output
등속 필드 솔버 옵션을 사용할 때 유량 속도를 선택한 데이터로 출력 할 수 있습니다.
벽 접착력으로 인한 지오메트리 구성 요소의 토크는 기존 벽 접착력 출력과 함께 별도의 수량으로 일반 이력 데이터에 출력됩니다.
난류 모델 출력이 요청 될 때 난류 에너지 및 소산과 함께 전단 속도 및 y +가 선택된 데이터로 자동 출력됩니다.
공기 유입 모델 출력에 몇 가지 수량이 추가되었습니다.자유 표면을 포함하는 모든 셀에서 혼입 된 공기 및 빠져 나가는 공기의 체적 플럭스가 재시작 및 선택된 데이터로 출력되어 사용자에게 공기가 혼입 및 탈선되는 위치 및 시간에 대한 자세한 정보를 제공합니다.전체 계산 영역 및 각 샘플링 볼륨 에 대해이 두 수량의 시간 및 공간 통합 등가물이 일반 히스토리 로 출력됩니다.
솔버의 출력 파일 flsgrf 의 최종 크기는 시뮬레이션이 끝날 때 보고됩니다.
2 유체 시뮬레이션의 경우, 기존의 출력 수량 유체 체류 시간 및 유체 가 이동 한 거리는 이제 유체 # 1 및 # 2와 유체의 혼합물에 대해 별도로 계산됩니다.
질량 입자의 경우, 각 종의 총 부피 및 질량이 계산되어 전체 계산 영역, 샘플링 볼륨 및 플럭스 표면에 대한 일반 히스토리 로 출력되어 입자 종 수에 대한 현재 출력을 보완합니다.
최종 로컬 가스 압력 은 사용자가 가스 포획을 식별하고 연료 탱크의 배기 시스템 설계를 지원하는 데 도움이되는 선택적 출력량으로 추가되었습니다.이 양은 유체로 채워지기 전에 셀의 마지막 공극 압력을 기록하며 단열 버블 모델과 함께 사용됩니다.
새로운 맞춤형 소스 루틴
새로운 사용자 정의 가능 소스 루틴이 추가되었으며 사용자의 개발 환경에서 액세스 할 수 있습니다.
소스 루틴 이름
기술
cav_prod_cal
Cavitation 생성과 소산 비율
sldg_uset
슬러지 침전 속도
phchg_mass_flux
증발 및 응축으로 인한 질량 플럭스
flhtccl
유체 # 1과 # 2 사이의 열전달 계수
dsize_cal
2 상 흐름에서 동적 액적 크기 모델의 응집 및 분해 속도
elstc_custom
점탄성 유체에 대한 응력 방정식의 Source Terms
새로운 사용자 인터페이스
FLOW-3D 사용자 인터페이스는 완전히 새롭게 디자인되어 현대적이고 평평한 구조로 사용자의 작업 흐름을 획기적으로 간소화합니다.
Setup dock widgets
Physics, Fluids, Mesh 및 FAVOR ™를 포함한 모든 설정 작업이 지오 메트리 윈도우 주변에서 독 위젯으로 변환되어 모델 설정을 단일 탭으로 요약할 수 있습니다. 이러한 전환으로 인해 이전 버전의 복잡한 접이식 트리가 훨씬 깨끗하고 효율적인 메뉴 프레젠테이션으로 대체되어 사용자는 ModelSetup탭을 떠나지 않고도 모든 매개 변수에 쉽게 액세스 할 수 있습니다.
New Model Setup icons
새로운 모델 설정 디자인에는 설정 프로세스의 각 단계를 나타내는 새로운 아이콘이 있습니다.
New Physics icons
RSS feed
새 RSS 피드부터 FLOW-3D v12.0의 시뮬레이션 관리자 탭이 개선되었습니다. FLOW-3D 를 시작하면 사용자에게 Flow Science의 최신 뉴스, 이벤트 및 블로그 게시물이 표시됩니다.
Configurable simulation monitor
시뮬레이션을 실행할 때 중요한 작업은 모니터링입니다. FLOW-3Dv1.0에서는
사용자가 시뮬레이션을 더 잘 모니터링할 수 있도록 SimulationManager의 플로팅 기능이 향상되었습니다. 사용자는 시뮬레이션 런타임 그래프를 통해 모니터링할 사용 가능한 모든 일반 기록 데이터 변수를 선택하고
각 그래프에 여러 변수를 추가할 수 있습니다. 이제 런타임에서 사용할 수 있는 일반 기록 데이터는
다음과 같습니다.
최소/최대
유체 온도
프로브 위치의 온도
유동 표면 위치에서의 유량
시뮬레이션 진단(예:시간 단계, 안정성 한계)
출입문에 유동 표면이 있는 대형 댐 Runtime plots of the flow rate at the gates of the large dam
Conforming 메쉬 시각화
용자는 이제 새로운 FAVOR ™ 독 위젯을 통해 적합한 메쉬 블록을 시각화 할 수 있습니다.
Large raster and STL data
데이터를 처리하는 데 걸리는 시간 때문에 큰 지오 메트리 데이터를 처리하는 것은 수고스러울 수 있습니다. 대형 지오 메트리 데이터를 처리하는 데는 여전히 상당한 시간이 걸릴 수 있지만, FLOW-3D는 이제 이러한 대규모 데이터 세트를 백그라운드 작업으로 로드하여 사용자가 데이터를 처리하는 동안 완전히 응답하고 중단 없는 인터페이스에서 작업을 계속할 수 있습니다
유체 흐름 문제는 복잡한 기하학적 구조의 자유 표면과 관련되는 경우가 많으며 대부분 매우 일시적입니다. 수력학의 예로는 배수로, 강, 교각 주변, 홍수 범람, 수문, 잠금 장치 및 다수의 기타 구조물의 흐름이 있습니다. 이러한 유형의 흐름을 계산적으로 모델링 하는 능력은 이러한 계산이 정확하고 합리적인 계산 자원으로 수행될 수 있다면 매력적입니다. 유용하게 사용하려면 시뮬레이션은 물리적 모델을 사용하는 것보다 훨씬 빠르고 저렴해야 합니다.
Fluid flow problems often involve free surfaces in complex geometry and in many cases are highly transient. Examples in hydraulics are flows over spillways, in rivers, around bridge pilings, flood overflows, flows in sluices, locks, and a host of other structures. A capability to computationally model these types of flows is attractive if such computations can be done accurately and with reasonable computational resources. To be useful, simulations should be much faster and less expensive than using physical models.
많은 컴퓨터 프로그램은 유체의 역학을 설명하는 편미분 방정식을 풀 수 있습니다. 시뮬레이션에 자유 표면을 포함 할 수있는 프로그램은 많지 않습니다. 그 이유는 Free Surface 경계 문제로 잘 알려진 수학적인 문제입니다. 자유 경계 문제는 다루기 어려운 표면이 이동함에 따라 계산 영역이 변화하는 한편, 그 표면 이동 자체가 계산에 의해 결정된다는 점에 있습니다. 계산 영역의 변화는 그 크기와 모양의 변화뿐만 아니라, 경우에 따라서는 영역의 결합과 분리(즉, 자유 표면의 발생과 소멸)을 포함합니다.
Many computer programs can solve the partial differential equations describing the dynamics of fluids. Not many programs are capable of including free surfaces in their simulations. The difficulty is a classical mathematical one often referred to as the free-boundary problem. A free boundary poses the difficulty that on the one hand the solution region changes when its surface moves, and on the other hand, the motion of the surface is in turn determined by the solution. Changes in the solution region include not only changes in size and shape, but in some cases, may also include the coalescence and break up of regions (i.e., the loss and gain of free surfaces).
이 책에서는 모든 자유 표면을 고려한 유체흐름 현상을 수치 해석용으로 모델링하는 방법에 대해 설명합니다. 이 기술은 VOF (Volume-of-Fluid) 법에 근거한 것으로, 특히 자유 표면 흐름에 적합한 다양한 기능을 제공합니다. 이 책에서는 VOF 법이 자유 표면과 그 발생과 소멸을 해석하는데 가장 자연스럽고 매우 효율적인 방법을 제시합니다.
In this note a computational modeling technique for fluid flows with arbitrary free surfaces is discussed. The technique is based on the Volume-of-Fluid (VOF) technique. This technique has many unique properties that make it especially applicable to flows having free surfaces. The goal of this discussion is to show why the VOF approach offers a natural way to capture free surfaces and their evolution with great efficiency.
VOF 법의 특징을 잘 보여주기 위해 간단하지만 매우 중요한 유동 현상에 관한 문제를 다룹니다. 여기에서는 계단 낙차형상의 낙하류를 예로 들어 있습니다. 개념적으로 간단한 흐름인 동시에 결과의 타당성을 확인하기위한 좋은 실험 데이터도 제공되어 있습니다 (N. Rajaratnam and MR Chamani “Energy Loss at Drops”J. Hydraulic Res. Vol. 33 p.373,1995 참조).
A good recommendation for the VOF method is to demonstrate its capabilities on a simple hydraulic flow problem, one that is far from trivial. The example selected is of flow over a step. This flow has conceptual simplicity and good experimental data available for validation (see N. Rajaratnam and M.R. Chamani, “Energy Loss at Drops,” J. Hydraulic Res. Vol. 33, p.373, 1995).
Prototype Hydraulic Flow with Free Surfaces
그림 1a는 정상 상태에 도달 한 후 흐름의 문제를 보여줍니다. 계단 낙차형상 상부로부터의 월류(액체 또는 스냅 시트)에는 상하 모두의 자유 표면이 있습니다. 월류의 아래쪽에는 월류와 계단 가공면 사이에 웅덩이가 형성되어 있으며, 하류에서는 액체는 평평한 정상 표면에서 오른쪽으로 흐르고 있습니다. 엄밀히 말하면, 웅덩이 영역의 흐름 상태는 정상입니다. 이것은 충돌하는 액체에 의해 풀에 난류 혼합이 발생하고 있기 때문입니다. 그러나 평균적인 구성이 존재하고 그것은 실험에서도 보고됩니다.
Figure 1a shows the flow problem after it has reached a steady-state condition. The overflow (sheet of liquid or nappe) leaving the top of the step has both an upper and lower free surface. At the bottom of the overflow a pool has formed between the overflow and the face of the step, while downstream, liquid is flowing to the right with a flat, steady surface. Strictly speaking, the flow conditions in the pool region are not steady because turbulent mixing is generated in the pool by the impinging fluid. There is, however, an average configuration and that is what is reported in the experiments.
실용적인 목적 유동 흐름은 항상 2 차원입니다. 즉, 그림 1a에서 수직 방향에서는 큰 변화는 없습니다. 현실에서는 웅덩이 위쪽으로 공간을 만들기 위해서는 대기에 여유공간이 필요하고, 그게 없으면 닫힐 것입니다.
For all practical purposes the flow is two-dimensional, that is, it does not have any significant variation in the direction normal to the illustration in Fig. 1a. In actuality, to have an air space above the pool there must be some opening to the atmosphere otherwise it would close up.
계단 낙차형상 상단의 유속은 중요합니다. 즉, 이것은 표면파와 같거나 그 이상의 속도이기 때문에 하류에서의 교란이 영역을 관통하고 상류 흐름 (계단 낙차형상의 왼쪽)에 영향을 줄 수 없습니다. 따라서 이 영역에서의 흐름은 예외적으로 원활하고 정상입니다.
The flow speed at the top of the step is critical, that is, it has a speed equal to or greater than the speed of surface waves, so that no disturbances from downstream can penetrate through this region to affect flow upstream (to the left of the step), which is why the flow is exceptionally smooth and steady in that region.
이 문제는 수치 시뮬레이션과 비교할 수 있는 기하 형상 기능이 많이 있습니다. 예를 들어, 계단 낙차형상의 전후 흐름의 높이, 월류가 바닥에 충돌 할 때의 각도, 월류 아래에 형성되는 웅덩이의 깊이 등입니다. 또한 실용화를 위한 중요한 비교 항목으로는, 계단 낙차형상을 통해 떨어지는 낙하 류에 의해 손실되는 에너지의 양 (운동 에너지와 위치 에너지의 합)가 있습니다.
There are many geometric features in this problem that can be compared with a numerical simulation; such as flow heights before and after the step, the angle of the overflow stream when it strikes the bottom and the depth of the pool formed under the overflow. Additionally, an important comparison for practical applications is the amount of energy (i.e., kinetic plus potential) lost by the flow in passing over the step.
Simulation of Prototype Problem
그림 1a는 시뮬레이션의 결과입니다. 이 예에서는 실험에 사용된 모든 기하 형상 및 물질의 특성이 시뮬레이션에 사용되었습니다. 실험실 테스트에서 사용한 계단 낙차형상의 높이가 62cm에서 액체는 보통의 물 (밀도 = 1.0gm / cc 어떻게 점성 = 0.01dynes / cm)입니다. 계산 영역에 들어가는 물의 깊이는 15.5cm에서 속도가 임계에 가까운 123.0cm/s 였습니다. 물론, 중력은 수직 방향으로 크기는 g = -980cm / s^2입니다.
Figure 1a is from a simulation. For this example all of the geometric and material properties used in the experiments were used in the simulation. The height of the step used in the laboratory test is 62cm and the fluid is ordinary water (density=1.0 gm/cc and dynamic viscosity=0.01dynes/cm). The depth of water entering the computational region was 15.5cm and was given a near critical velocity of 123.0cm/s. Of course, gravity was in the vertical direction with magnitude g=-980cm/s^2.
Figure 1a. Simulation of flow over a step.
Figure 1b. Grid used in simulation.
월류 왼쪽에 있는 웅덩이에 난류가 발생 할 것으로 예상 되었기 때문에, 시뮬레이션에서는 난류 모델 (the Renormalization Group, 즉 RNG 모델)을 사용했습니다. 그 후, 난류 모델을 사용하지 않고 한 시뮬레이션에서도 비슷한 결과를 얻을 수 있었지만, 이것은 그다지 놀라운 일이 아닙니다. 흐름의 중요한 요소의 대부분은 매끄러운 (즉 난류가 아닌) 유입, 유출, 월류 때문입니다.
Because some turbulence was expected to develop in the pool to the left of the overflow, a turbulence model (the Renormalization Group or RNG model) was used in the simulation. Subsequent simulations without a turbulence model produced very similar results, which is not too surprising since most of the important elements of the flow are smooth (i.e., non-turbulent) inflow, overflow and outflow streams.
그림 1b 시뮬레이션 영역은 폭 170cm, 높이 100cm에 가로 80 개, 세로 60 개, 총 4800 개의 셀로 구성되는 같은 크기의 사각형 셀의 격자로 세분화되어 있습니다. 이 격자는 유체 역학의 지배 미분 방정식 (나비에 – 스토크스 방정식)의 유한 차분 근사의 기초로 사용됩니다. 격자 셀의 수와 크기는 흐름 속에서 예측되는 최소의 특성을 파악하는 목적으로 선택되었습니다. 결과를보고 어떤 조정이 필요하다고 생각되는 경우는 숫자를 쉽게 늘리거나 줄일 수 있습니다. 사실, 해상도를 바꾸어 시뮬레이션을 반복하여 계산이 그러한 변화에 영향을 많이 들어 있지 않은지 확인하는 것이 좋습니다.
The simulation region shown in Fig. 1b is 170cm wide and 100cm high and has been subdivided into a grid of equal sized rectangular cells consisting of 80 cells in the horizontal direction and 60 cells in the vertical direction, for a total of 4800 cells. This grid is used as the basis for finite-difference approximations of the governing differential equations of fluid dynamics (the Navier-Stokes equations). The number and size of the grid cells was chosen with the goal of capturing the smallest expected features of the flow. The number can be easily increased or decreased if the results seem to warrant some adjustment. In fact, it is often a good idea to repeat a simulation with a change of resolution to make sure that the solution is not too sensitive to such changes.
왼쪽의 경계는 지정된 속도 경계입니다 (유체의 높이도 지정). 오른쪽의 경계는 유출 경계에서 모든 유량이 경계에 수직 제로 기울기이며, 균일 한 유출이 촉진됩니다. 상하 경계는 단단한 벽으로 세 번째 방향의 경계는 대칭면 (점성 저항 제로의 벽)으로 처리되었습니다. 계단 낙차형상의 표면도 자유-미끄럼(free slip) 경계로 처리되었습니다.
The left boundary was a specified velocity boundary (also with a specified fluid height). The right boundary was an outflow boundary where all flow quantities have a zero gradient normal to the boundary to encourage a uniform outflow. The top and bottom boundaries are rigid walls, while in the third direction the boundaries were treated as planes of symmetry (i.e., walls with zero viscous drag). The surface of the step was also treated as a free-slip boundary.
초기 조건은 예측되는 흐름의 배열을 대략적으로 근사하도록 설정할 수 있었지만, 흐름의 구성은 계산하고 싶은 것 중 하나이기 때문에 유체가 어떻게 분포되는지를 모르는 경우에는 간단한 방법이 필요합니다. 이 예제에서는 비정상 흐름 시뮬레이터를 사용했기 때문에 그림 1a의 계단 낙차형상에 유체의 블록만 있고 왼쪽 경계의 같은 수평 속도와 높이가 할당된 간단한 초기 조건을 정의할 수 있습니다. 시뮬레이션은 이후 정상 흐름으로 발전하고 있지만, 이것은 약 8.0 초 후에 발생합니다. 시뮬레이션은 정상 상태에 도달 한 것을 보장하기 위해, 10.0 초의 시간까지 실행되었습니다. 그림 2는 중간 시간을 두 보여줍니다. 도 2b는 0.2 초, 그림 2c는 0.5 초 시점에서 그림 2d는 마지막 10.0 초 시점을 보여줍니다.
Initial conditions could have been set to roughly approximate the expected flow arrangement, but since the flow configuration is one of the things that one would like to compute, especially for situations where one doesn’t know what the distribution of fluid is likely to be, a simpler approach is needed. Because a transient flow simulator was used for this example a simple initial condition could be defined that consisted of just a block of fluid on top of the step, Fig. 1a with the same horizontal velocity and height assigned to the left boundary. The simulation then followed the development of the steady flow, which occurs after about 8.0s. The simulation was run out to a time of 10.0s to assure that steady conditions had been reached. Figure 2 shows two intermediate times; 2.b at 0.2s and 2.c at 0.5s plus the final time in 2.d at 10.0s.
Figures 2a-2d. Simulation times of 0.0, 0.2, 0.5 and 10.0s.
처음에는 단일 결합하고 있는 자유 표면이었던 것이 액체가 바닥에 충돌한 후 2 개의 독립적인 자유 표면 (상하 스냅 표면)으로 변화하는 것에 주목하십시오. 아래 경계의 충격점의 좌우로 흐름이 분리되도 문제는 없습니다. 이에 대해서는 다음 섹션에서 자세히 설명합니다.
It should be noted that what starts as a single, connected free surface changes to two independent free surfaces (upper and lower nappe surfaces) after the fluid strikes the bottom. No difficulties are experienced with this separation of the flow into portions flowing to the left and right of the impact point on the bottom boundary. This will be discussed at further length in the next section.
실험과 시뮬레이션의 비교는 다음 표와 같으며 매우 잘 일치하고 있습니다.
Comparisons between experiment and simulation are given in the following table and are in excellent agreement.
Comparison Table
Experimental Results
Simulation Results
Outflow Height/Step Height
0.094
0.094
Pool Height/Step Height
0.41
0.41
Angle of Nappe at Bottom
57°
59°
Energy Loss/Initial Energy
0.29
0.296
이러한 결과를 고려하면이 같은 정밀도를 달성하려면 상당한 계산시간이 필요할 것으로 생각될지도 모릅니다. 그러나 실제로는 Pentium 4, 3.20GHz의 데스크톱 컴퓨터의 총 CPU 시간은 단 88 초였습니다. 계산시간이 너무 짧은 것은 설명이 필요하며, 이것은 다음 섹션의 목적입니다.
In view of these results it might be expected that a considerable amount of computational time would be required to achieve such accuracy. In fact, the total cpu time on a desktop Pentium 4, 3.20GHz computer was only 88s. Such a short computational time requires explanation and that is the purpose of the following sections.
Figures 2a-2d. Simulation times of 0.0, 0.2, 0.5 and 10.0s.
Why the VOF Technique Works Well / VOF 법이 적합한 이유
VOF 법의 구조와 그것이 매우 효율적인 방법인 이유를 이해하기 위해 다양한 계산법 중에서도 특히 VOF 법에 대한 몇 가지 기본 개념을 나타냅니다.
There are a few general concepts about computational methods and the VOF technique in particular that can be used to gain an understanding of how and why VOF works so efficiently.
Basic Theory
모든 수치해석 방법에서 흐름의 문제를 단순하게 산술 계산하도록 유한의 수치 세트로 단순화해야합니다. 연속 유체를 이산화된 수치 세트에 근사하기 위해서 일반적으로 사용되는 것이 유체가 차지하는 공간을 격자로 분할하는 방법입니다. 이 격자는 일반적으로 다수의 작은 직사각형의 블록(요소)로 구성됩니다. 이러한 각 요소에 대해 평균화 처리를 실시함으로써 그 요소의 유체의 압력, 밀도, 속도 및 온도의 대표 값을 얻을 수 있습니다.
All numerical methods must use some simplification to reduce a fluid flow problem to a finite set of numerical values that can then be manipulated using elementary arithmetical operations. A typical procedure for approximating a continuous fluid by a discrete set of numerical values is to subdivide the space occupied by the fluid into a grid consisting of a set of small, often rectangular “bricks.” Within each element an averaging process is applied to obtain representative element values for the fluid’s pressure, density, velocity and temperature.
간단한 수식을 사용해, 어느 시간에 걸친 각 요소 값과 인접한 요소의 상호 작용을 근사할 수 있습니다. 예를 들어, 요소의 밀도는 그 요소와 인접 요소 사이에서 (질량 보존에 의한) 질량 유량이 교환된 경우에만 변경됩니다. 요소 사이에서 질량이 교환되는 물질의 속도는 운동량 보존 법칙에 의해 계산되며 일반적으로 나비에-스토크스 방정식으로 표현됩니다. 나비에-스토크스 방정식은 인접한 요소 사이에 작용하는 압력과 점성 응력을 이용하여 요소에서 변화하는 유체 속도를 근사합니다.
Simple equations can be devised to approximate how each element’s values interact with neighboring elements over time. For instance, the density of an element can only change when there is a net flow of mass exchanged between an element and its neighbors (i.e., conservation of mass). The material velocity that carries mass between elements is computed from the conservation of momentum principal, usually expressed in the form of the Navier-Stokes equations, which uses the pressures and viscous stresses acting between neighboring elements to approximate the changing fluid velocities in the elements.
이러한 요소와 인접 요소 사이의 상호 작용에 따른 아이디어는 편미분 방정식 근방의 양의 변화에 의해 생기는 작은 변화의 효과를 평가하는 것과 본질적으로 동일합니다. 공학계의 교과서에서 파생된 작은 컨트롤 볼륨을 사용하여 그 크기를 무한대까지 작게 한 근사치의 극한으로 편미분 방정식이 유도됩니다. 수치 시뮬레이션에서도 같은 방식을 취하고 있지만, 요소 수가 너무 많으면 추적이 어렵게 되어 컨트롤 볼륨의 크기를 최대한 작게 만들 수 없습니다. 실제 시뮬레이션 현상을 해결하는데 충분하고 계산 시간을 최소한으로 억제 할 수 있는 요소수를 설정하는 것이 목표입니다.
This idea of an element interacting with its neighbors is essentially what is meant by a partial differential equation; that is, evaluating the effects of small changes caused by the variation in quantities nearby. Partial differential equations are typically derived in engineering text books as the limit of approximations made with small control volumes whose sizes are then reduced to infinitesimal values. In a numerical simulation the same thing is done except that the control volume sizes cannot be taken to the limit because that would require too many elements to keep track of. In practice, the goal is to use enough elements to resolve the phenomena of interest, and no more, so that computing times are kept to a minimum.
요소에 사용되는 연산은 기본적으로 더하기, 빼기, 곱하기 및 나누기만 포함된 간단한 것입니다. 예를 들어, 요소의 질량의 변화는 일정한 시간 간격에 걸쳐 요소의 측면에서 유입 및 유출된 질량의 가산 및 감산에서 구할 수 있습니다. 그러나 시뮬레이션에서는 이러한 연산을 수천, 때로는 수백만 요소에 대해 매우 짧은 시간 간격에 대해 반복 계산해야합니다. 따라서 이러한 반복 계산의 고속 처리는 컴퓨터가 적합합니다.
Arithmetical operations associated with an element generally involve only simple addition, subtraction, multiplication and division. For instance, the change of mass in an element involves the addition and subtraction of mass entering and leaving through the faces of the element over a fixed interval of time. A simulation requires that these operations be done for thousands or even millions of elements as well as repeated for many small time intervals. Computers are ideal for performing these types of repetitive operations very rapidly.
자유 표면을 수반하는 유체 운동의 시뮬레이션에서는 형상이 변화하는 계산 영역을 다루어야합니다. 이 복잡성에 대응할 수있는 분석 방법이 아래에서 설명하는 VOF 법입니다.
Simulating fluid motion with free surfaces introduces the complexity of having to deal with solution regions whose shapes are changing. A convenient way to deal with this is to use the Volume of Fluid (VOF) technique described next.
The VOF Concept
VOF 법은 각 격자 셀의 체적 중 액체가 차지하는 비율, 즉 체적 점유율을 기록한다는 생각에 근거합니다. 일반적으로 부피 점유율은 F로 표시됩니다. F는 부피 점유율이기 때문에 값이 취할 수있는 범위는 0.0 ~ 1.0입니다.
The VOF technique is based on the idea of recording in each grid cell the fractional portion of the cell volume that is occupied by liquid. Typically the fractional volume is represented by the quantity F. Because it is a fractional volume, F must have a value between 0.0 and 1.0.
액체 내부의 영역에서는 F 값은 1.0이 액체의 외부, 즉 (공기 등) 기체 영역에서 F 값은 0입니다. F 값이 0.0과 1.0 사이에서 변화하는 장소가 자유 표면이 존재하는 위치입니다. 즉 0.0보다 크고 1.0보다 작은 F 값을 가지는 요소는 반드시 표면을 가지고 있습니다.
In interior regions of liquid the value of F would be 1.0, while outside of the liquid, in regions of gas (air for example), the value of F is zero. The location of a free surface is where F changes from 0.0 to 1.0. Thus, any element having an F value lying between 0.0 and 1.0 must contain a surface.
여기서 유의해야 할 것은 VOF 법에서 자유 표면을 직접적으로 정의하는 것이 아니라 벌크 유체의 위치를 정의한다는 점입니다. 이렇게하면 계산상의 어려움을 초래하지 않고 유체 영역을 결합 또는 분할 할 수 있습니다. 자유 표면은 단순히 유체의 체적 점유율이 1.0과 0.0 사이에서 변화하는 장소로 정의됩니다. 이것은 자유 표면을 수반하는 거의 모든 문제에 적용 할 수 VOF 법의 뛰어난 특징이기도합니다.
It is important to emphasize that the VOF technique does not directly define a free surface, but rather defines the location of bulk fluid. It is for this reason that fluid regions can coalesce or break up without causing computational difficulties. Free surfaces are simply a consequence of where the fluid volume fraction passes from 1.0 to 0.0. This is a very desirable feature that makes the VOF technique applicable to just about any kind of free surface problem.
또한 격자의 각 요소에 단일 수치 (F)를 할당하여 유체의 위치를 기록 할 수 있는 점도 VOF 법의 중요한 특징입니다. 이것은 평균값을 기준으로 압력과 속도 등 다른 모든 유체 물성의 기록과 완전히 일치합니다.
Another important feature of the VOF technique is that it records the location of fluid by assigning a single numerical value (F) to each grid element. This is completely consistent with the recording of all other fluid properties in an element such as pressure and velocity components by their average values.
Some Details of the VOF Technique
Figure 3. Surface in 1D column of elements.
정확도를 위해 요소 내에 자유 표면을 배치하는 방법을 갖는 것이 바람직합니다. 인접 요소의 F 값을 고려하면 이를 쉽게 할 수 있습니다. 예를 들어, 열의 일부에 액체가 충전되어있는 1 차원 요소를 상상하십시오 (그림 3). 액체의 표면은 열 중앙 영역의 요소에 있습니다. 이것을 표면 요소라고합니다. 여기에서는 표면 요소를 제외하고 F 값은 0.0 또는 1.0이어야한다고 가정하고 있기 때문에 이를 사용하여 표면의 정확한 위치를 파악할 수 있습니다. 우선, 표면이 표면 또는 바닥을 확인하는 테스트를 실시합니다. 표면요소에 대해 액체가 없을 경우에는 표면으로 간주합니다. 위의 요소에 액체가 들어있는 경우는 물론, 그 표면은 바닥입니다. 윗면에 관해서는 정확한 위치는 표면 요소의 아래쪽에서 위쪽으로 요소의 세로 크기를 F 배 한 거리에있는로 계산합니다. 바닥도 마찬가지로 표면 요소의 상단에서 아래로, 요소의 세로 크기를 F 배 한 거리에 있습니다. 이 방법에 의한 요소의 표면 위치의 특정은 요소 내의 액체의 부피 점유율로 F를 정의한 후에 합니다.
For accuracy purposes it is desirable to have a way to locate a free surface within an element. Considering the F values in neighboring elements can easily do this. For example, imagine a one-dimensional column of elements in which a portion of the column is filled with liquid, Fig. 3. The liquid surface is in an element in the central region of the column, which will be referred to as the surface element. Because we assume the values of F must be either 0.0 or 1.0, except in the surface element, we can use this to locate the exact position of the surface. First a test is made to see if the surface is a top or bottom surface. If the element above the surface element is empty of liquid, the surface must be a top surface. It the element above is full of liquid then, of course, the surface is a bottom surface. For a top surface we compute its exact location as lying above the bottom edge of the surface element by a distance equal to F times the vertical size of the element. A bottom surface is similarly located a distance equal to F times the vertical size of the element below the top edge of the surface element. Locating the surface within an element in this way follows from the definition of F as a fractional volume of liquid in the element.
1 차원 열의 표면 위치 계산은 간단하고 정확하며 계산이 거의 필요없습니다. 그러나 2 차원 및 3 차원의 경우 하나의 표면 셀에 연속적인 표면 방향이 존재할 가능성이 있기 때문에 위치 계산은 조금 복잡해집니다. 그럼에도 불구하고 이를 취급하는 것은 어렵지 않습니다. 그림 4의 이차원의 예는 표면의 위치를 계산할 뿐만 아니라 경사와 곡률도 이해할 수 있는 쉬운 방법을 보여줍니다.
Calculating surface locations in one-dimensional columns is simple, accurate and requires very little arithmetic. In two and three dimensional situations, however, computing a location is a little more complicated because there is a continuous range of surface orientations possible within a surface cell. Nevertheless, dealing with this is not difficult. A two-dimensional example, Fig. 4, will illustrate a simple way to not only compute the location of the surface, but also to get a good idea of its slope and curvature.
Figure 4. Surface in 2D grid of elements.
1 차원의 경우처럼 먼저 인근 요소를 테스트하여 표면의 대략적인 방향을 찾아야합니다. 그림 4는 바깥 쪽의 법선이 상승 방향에 가장 가깝게 됩니다. 이것은 그 방향 밖의 값의 차이가 다른 방향보다 크기 때문입니다. 그럼 거의 수직으로 있는 요소 열에서 표면의 국소적인 높이가 계산됩니다. 그림 4의 2 차원의 경우에는 이러한 높이가 화살표로 표시되어 있습니다. 마지막으로, 표면 요소를 포함하는 컬럼의 높이에 따라 그 요소의 표면의 위치를 확인합니다. 다른 2 개의 높이를 사용하면 국소적인 표면 경사와 표면 곡률을 계산할 수 있습니다.
As in the one-dimensional case, it is first necessary to find the approximate orientation of the surface by testing the neighboring elements. In Fig. 4 the outward normal would be closest to the upward direction because the difference in neighboring values in that direction is larger than in any other direction. Next, local heights of the surface are computed in element columns that lie in the approximate normal direction. For the two-dimensional case in Fig. 4 these heights are indicated by arrows. Finally, the height in the column containing the surface element gives the location of the surface in that element, while the other two heights can be used to compute the local surface slope and surface curvature.
3 차원에서도 동일한 절차를 사용하지만, 표면 요소의 주위에 있는 9개의 열에 대해 열 높이를 요구해야합니다. 필요한 계산은 조금 더 걸리지만, 주된 내용은 열의 간단한 덧셈과 경사와 곡률을 추구하는 열의 높이의 합과 차이가 있습니다. 이 토론을 토대로, 이제 자유 표면을 정의하는 데 필요한 모든 정보를 빠르고 쉽게 평가하기 위해 부분 유체 체적을 사용하는 방법을 알아야합니다.
In three-dimensions the same procedure is used although column heights must be evaluated for nine columns around the surface element. Although a little more computation is needed, it consists primarily of simple summations in the columns and then sums and differences of column heights for evaluating the slope and curvature. Based on this discussion, the reader should now see how the fractional fluid volume can be used to quickly and easily evaluate all the information needed to define free surfaces.
다루어야 할 문제가 앞으로 2 개 남아 있습니다. 하나는 그림 1 및 2와 같은 시뮬레이션은 유체가 존재하는 영역에는 유체 역학만으로 해결합니다. 이것은 VOF 법의 계산 효율이 높은 또 하나의 이유입니다. 계단 형상의 낙하류의 문제로 유체가 차지하는 영역은 계산 격자의 오픈 공간의 절반 이하입니다. 액체를 둘러싼 기체의 흐름을 계산할 필요가 있다면 필요한 계산 시간이 크게 늘어납니다. 그러나 액체만으로 계산을 할 경우 자유 표면 경계 조건을 지정해야합니다. 이 조건은 접선 응력의 소실과 기체의 압력에 동일한 표준 압력을 표면에 추가하는 것입니다.
There are two remaining issues to deal with. One issue is that a simulation like that in Figs. 1 and 2 is only solving for the fluid dynamics in regions where there is fluid. This is another reason for the computational efficiency of the VOF method. The region occupied by fluid in the flow over a step problem is much less than half of the open region in the computational grid. If it were necessary to also solve for the flow of gas surrounding the liquid, then considerably more computational time would be required. In order to perform solutions only in the liquid, however, it is necessary to specify boundary conditions at free surfaces. These conditions are the vanishing of the tangential stress and application of a normal pressure at the surface that equals the pressure of the gas.
두 번째 문제는 자유 표면이 유체와 함께 움직일 때의 움직임과 변형을 유체 점유율 변수 F를 구함으로써 계산해야 한다는 것입니다. 변수 F는 불연속 (주로 0.0 또는 1.0)이기 때문에 계산 격자를 이동할 때 이 불연속성이 유지되도록주의해야합니다. VOF 법은이 목적으로 특수 이류(advection) 알고리즘이 사용되고 있습니다.
A second issue is that movement and deformation of a free surface must be computed by solving for the fraction of fluid variable, F, as it moves with the fluid. Because the variable F is discontinuous (i.e., primarily 0.0 or 1.0) some care must be taken to maintain this discontinuity as it moves through a computational grid. In the VOF method, special advection algorithms are used for this purpose.
Illustration of Free-Surface Tracking by VOF Technique
그림 6a는 이것의 적합 여부를 보여줍니다. 유체의 체적 점유율은 격자 요소마다 균일하게 분류되고 그 요소의 값을 나타냅니다. 자유 표면은 거의 모든 곳에서 선명하게 정의되어 있습니다. 스냅의 가장 낮은 가장 좁은 부분에만 선명한 유체 분포의 손실을 확인할 수 있습니다 (그림 5b). 이것은 예상대로입니다. 이 영역에서는 스냅의 두께는 3 가지 요소보다 작고, 따라서 부분 충전된 표면 요소에 연결된 작은 F 값이 어떤 중심 요소 (값 1.0)에 혼입하기 때문입니다. 계산 목적으로 이 것은 별로 문제가 되지 않습니다. 이 시뮬레이션 방법은 액체 내부의 요소는 순수한 액체 성분과 같은 방식으로 처리되기 때문입니다.
Figure 6a is an illustration of how well this works; the fluid volume fraction is colored uniformly in each grid element to represent its value in that element. The free surface is sharply defined nearly everywhere. Only in the lowest and narrowest part of the nappe is there any noticeable loss of a sharp fluid fraction distribution, Fig. 5b. This was expected because in this region the nappe is less than three elements in thickness and this allows some of the smaller F values associated with partially filled surface elements to mix in with the central element, which should have a value of 1.0. For computational purposes this doesn’t really matter because the simulation method treats elements interior to the liquid as though they are pure liquid elements.
그림 5b에 나타내는 영역에서는 실제 실험에서 난류 및 공기 혼입이 관찰된 것도 지적해 두지 않으면 안됩니다. 따라서 유체 점유율의 값을 1보다 조금 작게 보이는 것이 다소 현실적입니다. 이것은 전혀 의외라는 것은 없습니다. 난류와 공기 유입을 담당하는 풀의 액체 제트의 교점은 난류와 공기 유입의 원인이 되지만, 유체 점유율 값(fluid fraction values )은 액체 내부에 “유입” 원인이 되기 때문에 실수가 아닙니다.
It should also be pointed out that in the region shown in Fig. 5b turbulence and air entrainment are observed in actual experiments. Thus, the appearance of fluid fraction values a little less than unity is somewhat realistic. This is not entirely accidental because the intersection of jet of liquid with a pool, which is responsible for turbulence and air entrainment, is also responsible for the “entrainment” of fluid fraction values into the interior of the liquid.
Figure 5a (left): Fluid fraction values in elements, showing sharpness of surface definition. Figure 5b (right): Close up of fluid fraction values where the overflow hits bottom.
Summary
처음에는 컴퓨터가 단순히 반복적인 산술 연산을 수행하고, 복잡하고 시간에 의존적인 유체 역학 문제에 대해, 현실적인 시뮬레이션을 할 수 있다는 것이 다소 마술처럼 보일 수 있습니다. 이 논의의 목적은 비교적 기본적인 절차로 이를 수행하는 접근법을 설명하는 것입니다.
간단하지만 사소한 유압 흐름 예제를 사용하여 계산된 시뮬레이션이 물리적인 측정 결과와 매우 일치하는 세부 결과를 생성 할 수 있음이 입증되었습니다. VOF (Volume of Fluid) 기술을 기반으로 한 시뮬레이션은 정확하고, 매우 효율적인 것이 추가로 입증되었습니다.
분명하게, 수력 발전소에서 사용되는 것과 같은 복잡한 유압 구조와 관련된 실제 예는 유용한 결과를 얻기 위해서는 이 예에서 사용되는 몇 초 이상의 많은 계산 시간을 소비해야합니다. 그럼에도 불구하고 이러한 결과는 합리적인 시간 (사람과 컴퓨터 모두)에서 수행 될 수 있으며, 실제 실험에서는 거의 불가능한 세부 사항들을 포함합니다. 또한, 지오메트리, 유동 조건 또는 유체 특성의 거의 모든 종류의 변화의 영향을 쉽게 테스트 할 수있는 능력은 시뮬레이션을 사용하는 또 다른 강력한 이유입니다. 기술의 발전에 따라 hydraulic flow 시뮬레이션을 위한 현재 소프트웨어 및 하드웨어는 기존의 물리적 모델링에 비해 상당한 비용 이점을 제공합니다.
At first it may seem somewhat magical that a computer can simply perform repeated arithmetic operations on arrays of numbers and produce a realistic simulation of a complex, time-dependent, fluid dynamics problem. It was the purpose of this discussion to explain an approach that does this with relatively elementary procedures.
Using a simple, but non-trivial, hydraulic flow example it has been demonstrated that computational simulations can produce detailed results in excellent agreement with physical measurements. It has been further demonstrated that the simulation, which was based on the Volume of Fluid (VOF) technique, uses simple approximation methods that are both accurate and efficient.
Clearly, real world examples involving complex hydraulic structures such as those used in hydroelectric power stations, must consume more than the few seconds of computational time used in our example to obtain useful results. Nevertheless, those results can be generated in reasonable times (both man and computer) and contain a richness of detail rarely possible in physical experiments. For examples visit our water and environmental application pages. In addition, the ability to easily test the influence of just about any kind of change in geometry, flow condition or fluid property is another powerful reason to employ simulations. Current software and hardware for hydraulic flow simulations offer a significant cost advantage over traditional physical modeling.
Postscript
The first detailed description of the VOF method was in 1981 by C.W. Hirt and B.D. Nichols, J. Comp. Phys., 39, p.201. All simulations appearing in this article were performed with the commercial software package FLOW-3D developed by Flow Science, Inc. This program uses an enhanced variant of the VOF concept called TruVOF.
본 자료는 국내 사용자들의 편의를 위해 원문 번역을 해서 제공하기 때문에 일부 오역이 있을 수 있어서 원문과 함께 수록합니다. 자료를 이용하실 때 참고하시기 바랍니다.
커튼 코팅에서 액체는 슬롯이나 슬라이드 다이에서 흘러 나와 중력을 받아 수평으로 움직이는 기판에 떨어집니다. 커튼 코팅은 단일 층의 유체 또는 다중 유체 층을 포함 할 수 있으며 사진 필름, 특수 용지 및 포장의 생산에 사용됩니다.
슬롯에서 나오는 유속과 관련된 액체의 물리적 특성과 기질 속도에 따라 액체가 먼저 기질에 접촉하는 접촉선의 안정성은 물론 코팅 두께를 결정합니다. 접촉 라인이 불안정하면 코팅 아래에 웅덩이나 공기 유입이 발생하여 균일하지 않은 코팅 두께 및 기타 결함을 발생시킬 수 있습니다. FLOW-3D는 유량, 유량, 기질 속도 등 다양한 공정 파라미터로 공정을 시뮬레이션하고 점성, 표면 장력, 접착력 등 액체의 물리적 성질을 변화시켜 커튼 코팅의 안정성을 조사하는 데 사용할 수 있습니다.
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
View in article
