FLOW-3D 소프트웨어 제품군의 모든 제품은 2023R1에서 IT 관련 개선 사항을 받았습니다. FLOW-3D 2023R1은 이제 Windows 11 및 RHEL 8을 지원합니다. 누락된 종속성을 보고하도록 Linux 설치 프로그램이 개선되었으며 더 이상 루트 수준 권한이 필요하지 않으므로 설치가 더 쉽고 안전해집니다. 또한 워크플로를 자동화한 사용자를 위해 입력 파일 변환기에 명령줄 인터페이스를 추가하여 스크립트 환경에서도 워크플로가 업데이트된 입력 파일로 작동하는지 확인할 수 있습니다.
확장된 PQ 2 분석
제조에 사용되는 유압 시스템은 PQ 2 곡선을 사용하여 모델링할 수 있습니다. 장치의 세부 사항을 건너뛰고 흐름에 미치는 영향을 포함하기 위해 질량-운동량 소스 또는 속도 경계 조건을 사용하여 유압 시스템을 근사화하는 것이 편리한 단순화인 경우가 많습니다. 기존 PQ 2 분석 모델을 확장하여 이러한 유형의 기하학적 단순화를 허용하면서도 여전히 현실적인 결과를 제공합니다. 이것은 시뮬레이션 시간과 모델 복잡성의 감소로 해석됩니다.
FLOW-3D 2022R2 의 새로운 기능
FLOW-3D 2022R2 제품군 의 출시와 함께 Flow Science는 워크스테이션과 FLOW-3D 의 HPC 버전 을 통합하여 단일 노드 CPU 구성에서 다중 구성에 이르기까지 모든 유형의 하드웨어 아키텍처를 활용할 수 있는 단일 솔버 엔진을 제공합니다. 노드 병렬 고성능 컴퓨팅 실행. 추가 개발에는 점탄성 흐름을 위한 새로운 로그 구조 텐서 방법, 지속적인 솔버 속도 성능 개선, 고급 냉각 채널 및 팬텀 구성 요소 제어, 향상된 연행 공기 기능이 포함됩니다.
통합 솔버
FLOW-3D 제품을 단일 통합 솔버로 마이그레이션하여 로컬 워크스테이션 또는 고성능 컴퓨팅 하드웨어 환경에서 원활하게 실행했습니다.
많은 사용자가 노트북이나 로컬 워크스테이션에서 모델을 실행하지만 고성능 컴퓨팅 클러스터에서 더 큰 모델을 실행합니다. 2022R2 릴리스에서는 통합 솔버를 통해 사용자가 HPC 솔루션에서 OpenMP/MPI 하이브리드 병렬화의 동일한 이점을 활용하여 워크스테이션 및 노트북에서 실행할 수 있습니다.
솔버 성능 개선
멀티 소켓 워크스테이션
멀티 소켓 워크스테이션은 이제 매우 일반적이며 대규모 시뮬레이션을 실행할 수 있습니다. 새로운 통합 솔버를 통해 이러한 유형의 하드웨어를 사용하는 사용자는 일반적으로 HPC 클러스터 구성에서만 사용할 수 있었던 OpenMP/MPI 하이브리드 병렬화를 활용하여 모델을 실행할 수 있는 성능 이점을 볼 수 있습니다.
낮은 수준의 루틴으로 벡터화 및 메모리 액세스 개선
대부분의 테스트 사례에서 10%에서 20% 정도의 성능 향상이 관찰되었으며 일부 사례에서는 20%를 초과하는 런타임 이점이 있었습니다.
정제된 체적 대류 안정성 한계
시간 단계 안정성 한계는 모델 런타임의 주요 동인입니다. 2022R2에서는 새로운 시간 단계 안정성 한계인 3D 대류 안정성 한계를 숫자 위젯에서 사용할 수 있습니다. 실행 중이고 대류가 제한된(cx, cy 또는 cz 제한) 모델의 경우 새 옵션은 30% 정도의 일반적인 속도 향상을 보여주었습니다.
압력 솔버 프리 컨디셔너
경우에 따라 까다로운 흐름 구성의 경우 과도한 압력 솔버 반복으로 인해 실행 시간이 길어질 수 있습니다. 어려운 경우 2022R2에서는 모델이 너무 많이 반복될 때 FLOW-3D가 자동으로 새로운 프리 컨디셔너를 활성화하여 압력 수렴을 돕습니다. 테스트의 런타임이 1.9배에서 335배까지 빨라졌습니다!
점탄성 유체에 대한 로그 형태 텐서 방법
점탄성 유체에 대한 새로운 솔버 옵션을 사용자가 사용할 수 있으며 특히 높은 Weissenberg 수치에 효과적입니다.
활성 시뮬레이션 제어 확장
능동 시뮬레이션 제어 기능은 연속 주조 및 적층 제조 응용 프로그램과 주조 및 기타 여러 열 관리 응용 프로그램에 사용되는 냉각 채널에 일반적으로 사용되는 팬텀 개체를 포함하도록 확장되었습니다.
연행 공기 기능 개선
디퓨저 및 유사한 산업용 기포 흐름 응용 분야의 경우 이제 대량 공급원을 사용하여 물 기둥에 공기를 도입할 수 있습니다. 또한 혼입 공기 및 용존 산소의 난류 확산에 대한 기본값이 업데이트되었으며 매우 낮은 공기 농도에 대한 모델 정확도가 향상되었습니다.
Investigation of the Turbulent Schmidt Number Effects On Numerical Modelling Of Vortex-Type Stormwater Retention Ponds
S. M. Yamini1; H. Shamloo2; S. H. Ghafari3 1M.Eng., Dep. of Civil Engineering K.N. Toosi University of Technology, Valiasr St., Tehran, Iran. smyamini@alumni.kntu.ac.ir 2Associate Professor, Dep. of Civil Engineering K.N. Toosi University of Technology, Valiasr St., Tehran, Iran. hshamloo@kntu.ac.ir 3Ph.D., Dep. of Civil Engineering Univ. of Tehran, Enqelab St., Tehran, Iran. sarvenazghafari@ut.ac.ir
Abstract
정확하고 신뢰할 수 있는 CFD 모델링 결과를 얻는 것은 이러한 시뮬레이션에서 입력의 중요성 때문에 종종 정밀 조사의 대상입니다.
난류 모델링이 RANS(Reynolds-Averaged Navier-Stokes) 방정식을 기반으로 하는 경우 난류 스칼라 전송을 추정하려면 난류 흐름에서 질량 1에 대한 운동량 확산의 비율로 정의되는 난류 슈미트 수(Sct)의 정의가 필요합니다.
그러나 이 매개변수는 난류 흐름의 속성이므로 보편적인 값이 허용되지 않았습니다. 우수 저류지의 수치 연구에서 적절한 Sct를 설정하는 실제 역할은 수력 효율의 평가가 추적자 테스트의 출력 질량 농도를 기반으로 하기 때문에 가장 중요합니다.
본 연구에서는 FLOW-3D를 사용하여 와류형 우수 저류지의 여러 수치 시뮬레이션을 체계적으로 수행했습니다. 다양한 난류 슈미트 수의 범위는 메쉬 감도를 조사하기 위해 다른 수의 계산 셀에 의해 수행된 수치 시뮬레이션에 도입되었습니다.
또한 사용자 정의 또는 자동 계산 값으로 최대 난류 혼합 길이의 영향을 평가했습니다. 이 연구의 결과는 실험 결과와 밀접한 일치를 제공하는 Sct= 0.625와 함께 수리학적 직경의 7%와 동일한 최대 난류 혼합 길이의 일정한 값을 갖는 확립된 수치 모델입니다.
특히 수치적 무차원 RDT 곡선의 피크 값은 극적으로 감소하여 실험 결과와 거의 일치했습니다. 이것은 FLOW-3D가 난류 유동의 와류형 물리학에서 질량 확산도를 적절하게 예측하는 상당한 능력을 가지고 있다는 결론을 내립니다.
– Achieving accurate and reliable CFD modelling results often is the subject of scrutiny because of the importance of the inputs in those simulations. If turbulence modelling is based on Reynolds-Averaged Navier-Stokes (RANS) equations, estimating the turbulent scalar transport requires the definition of the turbulent Schmidt number (Sct), defined as the ratio of momentum diffusivity to mass one in a turbulent flow. However, no universal value has been accepted for this parameter as it is a property of turbulent flows.
The practical role of establishing a suitable Sct in numerical studies of stormwater retention ponds is of the utmost importance because the assessment of the hydraulic efficiency of them is based on output mass concentration of tracer tests. In this study, several numerical simulations of a vortex-type stormwater retention pond were systematically carried out using FLOW-3D. A range of various turbulent Schmidt numbers were introduced in numerical simulations performed by different number of computational cells to investigate mesh sensitivity.
Moreover, the effects of maximum turbulent mixing length as a user-defined or automatically computed value were assessed. The outcome of this study is an established numerical model with a constant value of maximum turbulent mixing length equal to 7% of the hydraulic diameter along with Sct= 0.625 which provides a close agreement with experimental results.
Noticeably, the peak values of numerical dimensionless RDT curves are dramatically decreased, resulted in a close match with experimental results. This concludes that FLOW-3D has a considerable ability to appropriately predict mass diffusivity in vortex-type physics of turbulent flows.
[1] C. Gualtieri, A. Angeloudis, F. Bombardelli, S. Jha, and T. Stoesser, “On the Values for the Turbulent Schmidt Number in Environmental Flows,” Fluids, vol. 2, p. 17, 2017. [2] Å. Adamsson, L. Bergdahl, and S. Lyngfelt, “Measurement and three-dimensional simulation of flow in a rectangular detention tank,” Urban Water Journal, vol. 2, no. 4, pp. 277-287, 2005/12/01 2005, doi: 10.1080/15730620500386545. [3] C. Gualtieri, “Numerical simulation of flow and tracer transport in a disinfection contact tank,” 2006. [4] S. Khan, B. Melville, and A. Shamseldin, Modeling the Layouts of Stormwater Retention Ponds using Residence Time. 2009, pp. 77-83. [5] F. Martínez-Solano, P. L. I. Rey, C. Gualtieri, and P. López-Jiménez, “Modelling flow and concentration field in rectangular water tanks,” 2010. [6] W. B. Rauen, A. Angeloudis, and R. A. Falconer, “Appraisal of chlorine contact tank modelling practices,” Water Research, vol. 46, no. 18, pp. 5834-5847, 2012/11/15/ 2012, doi: https://doi.org/10.1016/j.watres.2012.08.013.
[7] J. Zhang, A. Tejada-Martínez, and Q. Zhang, “Evaluation of LES and RANS for Determining Hydraulic Performance of Disinfection Systems for Water Treatment,” Journal of Fluids Engineering, vol. 136, 05/15 2014, doi: 10.1115/1.4027652. [8] J. Zhang, A. E. Tejada-Martínez, and Q. Zhang, “Developments in computational fluid dynamics-based modeling for disinfection technologies over the last two decades: A review,” Environmental Modelling & Software, vol. 58, pp. 71- 85, 2014/08/01/ 2014, doi: https://doi.org/10.1016/j.envsoft.2014.04.003. [9] C. Gualtieri and F. Salzano, “DIscussion on “The effect of baffle spacing on hydrodynamics and solute transport in serpentine contact tanks”,” Journal of Hydraulic Research, vol. 52, pp. 152-154, 02/28 2014, doi: 10.1080/00221686.2013.877528. [10] A. Angeloudis, T. Stoesser, R. A. Falconer, and D. Kim, “Flow, transport and disinfection performance in small- and full-scale contact tanks,” Journal of Hydro-environment Research, vol. 9, no. 1, pp. 15-27, 2015/03/01/ 2015, doi: https://doi.org/10.1016/j.jher.2014.07.001. [11] A. Angeloudis, T. Stoesser, C. Gualtieri, and R. A. Falconer, “Contact Tank Design Impact on Process Performance,” Environmental Modeling & Assessment, vol. 21, no. 5, pp. 563-576, 2016/10/01 2016, doi: 10.1007/s10666-016-9502- x. [12] D. Valero and D. B. Bung, “Sensitivity of turbulent Schmidt number and turbulence model to simulations of jets in crossflow,” Environmental Modelling & Software, vol. 82, pp. 218-228, 2016/08/01/ 2016, doi: https://doi.org/10.1016/j.envsoft.2016.04.030. [13] F. Sonnenwald, I. Guymer, and V. Stovin, “Computational fluid dynamics modelling of residence times in vegetated stormwater ponds,” Proceedings of the Institution of Civil Engineers – Water Management, vol. 171, pp. 1-11, 11/07 2017, doi: 10.1680/jwama.16.00117. [14] F. Sonnenwald, I. Guymer, and V. Stovin, “A CFD-Based Mixing Model for Vegetated Flows,” Water Resources Research, vol. 55, no. 3, pp. 2322-2347, 2019, doi: https://doi.org/10.1029/2018WR023628. [15] S. B. Pope, Turbulent Flows. Cambridge, UK: Cambridge University Press, 2000. [16] R. Rossi and G. Iaccarino, “Numerical simulation of scalar dispersion downstream of a square obstacle using gradienttransport type models,” Atmospheric Environment, vol. 43, no. 16, pp. 2518-2531, 2009/05/01/ 2009, doi: https://doi.org/10.1016/j.atmosenv.2009.02.044. [17] R. Chowdhury, M. Ahadi, K. A. Mazurek, G. Putz, D. Bergstrom, and C. Albers, “Physical Scale and Computational Modeling in the Development of a Vortex-Type Stormwater Retention Pond,” in World Environmental and Water Resources Congress 2016, 2016, pp. 388-397. [18] V. Yakhot and L. M. Smith, “The renormalization group, the ɛ-expansion and derivation of turbulence models,” Journal of Scientific Computing, vol. 7, no. 1, pp. 35-61, 1992/03/01 1992, doi: 10.1007/BF01060210. [19] Flow Science, Inc., FLOW-3D User manual. Santa Fe, NM, USA. (2015). [20] M. M. Bishop, J. M. Morgan, B. Cornwell, and D. K. Jamison, “Improving the Disinfection Detention Time of a Water Plant Clearwell,” Journal AWWA, vol. 85, no. 3, pp. 68-75, 1993, doi: https://doi.org/10.1002/j.1551- 8833.1993.tb05958.x. [21] F. L. Hart, “Improved Hydraulic Performance of Chlorine Contact Chambers.,” Jounal of Water Pollution Control Federation, vol. 51(12), pp. 2868–2875, 1979.
Won-Ik Cho, Peer Woizeschke Bremer Institut für angewandte Strahltechnik GmbH, Klagenfurter Straße 5, Bremen 28359, Germany
Received 30 July 2020, Revised 3 October 2020, Accepted 18 October 2020, Available online 1 November 2020.
Abstract
Molten pool flow and heat transfer in a laser welding process using beam oscillation and filler wire feeding were calculated using computational fluid dynamics (CFD). There are various indirect methods used to analyze the molten pool dynamics in fusion welding. In this work, based on the simulation results, the surface fluctuation was directly measured to enable a more intuitive analysis, and then the signal was analyzed using the Fourier transform and wavelet transform in terms of the beam oscillation frequency and buttonhole formation. The 1st frequency (2 x beam oscillation frequency, the so-called chopping frequency), 2nd frequency (4 x beam oscillation frequency), and beam oscillation frequency components were the main components found. The 1st and 2nd frequency components were caused by the effect of the chopping process and lumped line energy. The beam oscillation frequency component was related to rapid, unstable molten pool behavior. The wavelet transform effectively analyzed the rapid behaviors based on the change of the frequency components over time.
Korea Abstract
빔 진동 및 필러 와이어 공급을 사용하는 레이저 용접 공정에서 용융 풀 흐름 및 열 전달은 CFD (전산 유체 역학)를 사용하여 계산되었습니다. 용융 용접에서 용융 풀 역학을 분석하는 데 사용되는 다양한 간접 방법이 있습니다.
본 연구에서는 시뮬레이션 결과를 바탕으로 보다 직관적 인 분석이 가능하도록 표면 변동을 직접 측정 한 후 빔 발진 주파수 및 버튼 홀 형성 측면에서 푸리에 변환 및 웨이블릿 변환을 사용하여 신호를 분석했습니다.
1 차 주파수 (2 x 빔 발진 주파수, 이른바 초핑 주파수), 2 차 주파수 (4 x 빔 발진 주파수) 및 빔 발진 주파수 성분이 발견 된 주요 구성 요소였습니다. 1 차 및 2 차 주파수 성분은 쵸핑 공정과 집중 라인 에너지의 영향으로 인해 발생했습니다.
빔 진동 주파수 성분은 빠르고 불안정한 용융 풀 동작과 관련이 있습니다. 웨이블릿 변환은 시간 경과에 따른 주파수 구성 요소의 변화를 기반으로 빠른 동작을 효과적으로 분석했습니다.
1 . 소개
융합 용접에서 용융 풀 역학은 용접 결함과 시각적 이음새 품질에 직접적인 영향을 미칩니다. 이러한 역학을 연구하기 위해 고속 카메라를 사용하는 직접 방법과 광학 또는 음향 신호를 사용하는 간접 방법과 같은 다양한 측정 방법을 사용하여 여러 실험 방법을 고려했습니다. 시간 도메인의 원래 신호는 특별히 주파수 도메인에서 변환 된 신호로 변환되어 용융 풀 동작에 영향을 미치는 주파수 성분을 분석합니다. Kotecki et al. (1972)는 고속 카메라를 사용하여 가스 텅스텐 아크 용접에서 용융 풀을 관찰했습니다. [1]. 그들은 120Hz 리플 DC 출력을 가진 용접 전원을 사용할 때 용융 풀 진동 주파수가 120Hz임을 보여주었습니다. 전원을 끈 후 진동 주파수는 용융 풀의 고유 주파수를 나타내는 용융 풀 크기와 관련이 있습니다. 진동은 응고 중에 용접 표면 스케일링을 생성했습니다. Zacksenhouse and Hardt (1983)는 레이저 섀도 잉 동작 측정 기술을 사용하여 가스 텅스텐 아크 용접에서 완전히 관통 된 용융 풀의 동작을 측정했습니다 [2] . 그들은 2.5mm 두께의 강판에서 6mm 풀 반경 (고정 용접)에 대해 용융 풀의 고유 주파수가 18.9Hz라는 것을 발견했습니다. Semak et al. (1995) 고속 카메라를 사용하여 레이저 스폿 용접에서 용융 풀 및 키홀 역학 조사 [3]. 그들은 깊이가 약 3mm이고 반경이 약 3mm 인 용융 풀에서 200Hz의 낮은 체적 진동 주파수를 관찰했습니다. 0.45mm Aendenroomer와 den Ouden (1998)은 강철의 펄스 가스 텅스텐 아크 용접에서 용융 풀 진동을보고했습니다 [4] . 그들은 침투 깊이에 따라 진동 모드 변화를 보였고 주파수는 50Hz에서 150Hz 사이에서 변화했습니다. 주파수는 완전히 침투 된 용융 풀에서 더 낮았습니다. Hermans와 den Ouden (1999)은 단락 가스 금속 아크 용접에서 용융 풀 진동을 분석했습니다. [5]. 그들은 용융 풀의 단락 주파수와 고유 주파수가 같을 때 부분적으로 침투 된 용융 풀의 경우 공정 안정성이 향상되었음을 보여주었습니다. Yudodibroto et al. (2004)는 가스 텅스텐 아크 용접에서 용융 풀 진동에 대한 필러 와이어의 영향을 조사했습니다 [6] . 그들은 금속 전달이 특히 부분적으로 침투 된 용융 풀에서 진동 거동을 방해한다는 것을 보여주었습니다. Geiger et al. (2009) 레이저 키홀 용접에서 발광 분석 [7]. 신호의 주파수 분석을 사용하여 용융 풀 (1.5kHz 미만)과 키홀 (약 3kHz)에 해당하는 진동 주파수 범위를 찾았습니다. Kägeler와 Schmidt (2010)는 레이저 용접에서 용융 풀 크기의 변화를 관찰하기 위해 고속 카메라를 사용했습니다 [8] . 그들은 용융 풀에서 지배적 인 저주파 진동 성분 (100Hz 미만)을 발견했습니다. Shi et al. (2015) 고속 카메라를 사용하여 펄스 가스 텅스텐 아크 용접에서 용융 풀 진동 주파수 분석 [9]. 그들은 용접 침투 깊이가 작을수록 용융 풀의 진동 빈도가 더 높다는 것을 보여주었습니다. 추출 된 진동 주파수는 완전 용입 용접의 경우 85Hz 미만 이었지만 부분 용입 용접의 경우 110Hz에서 125Hz 사이였습니다. Volpp와 Vollertsen (2016)은 레이저 키홀 역학을 분석하기 위해 광학 신호를 사용했습니다 [10] . 그들은 공간 레이저 강도 분포로 인해 0.8에서 154 kHz 사이의 고주파 범위에서 피크를 발견했습니다. 위에서 언급 한 실험적 접근법은 공정 조건, 측정 방법 및 측정 된 위치에 따라 수십 Hz에서 수십 kHz까지 광범위한 용융 풀 역학에 대한 결과를 보여 주었다는 점에 유의해야합니다.
융합 용접에서 용융 풀 역학을 연구하기 위해 분석 접근 방식도 사용되었습니다. Zacksenhouse와 Hardt (1983)는 2.5mm 두께의 강판에서 대칭형 완전 관통 용융 풀의 고유 진동수를 계산했습니다 [2] . 매스 스프링 해석 모델을 사용하여 용융 풀 반경 6mm (고정 용접)에 대해 20.4Hz (실험에서 18.9Hz)의 고유 진동수와 3mm 풀 반경 (연속 용접)에 대해 40Hz의 고유 진동수를 예측했습니다. ). Postacioglu et al. (1989)는 원통형 용융 풀과 키홀을 가정하여 레이저 용접의 용융 풀에서 키홀 진동의 고유 진동수를 계산했습니다 .. 특정 열쇠 구멍 모양의 경우 약 900Hz의 기본 주파수가 계산되었습니다. Postacioglu et al. (1991)은 또한 레이저 용접에서 용접 속도를 고려하기 위해 타원형 용융 풀의 고유 진동수를 계산했습니다 [12] . 그들은 타원형 용융 풀의 모양이 고유 진동수에 영향을 미친다는 것을 보여주었습니다. 고유 진동수는 축의 길이 비율이 낮았으며, 즉 타원의 반장 축과 반 단축의 비율이 낮았습니다. Kroos et al. (1993)은 축 대칭 용융 풀과 키홀을 가정하여 레이저 키홀 용접의 동적 거동에 대한 이론적 모델을 개발했습니다 .. 키홀 폐쇄 시간은 0.1ms였으며 안정성 분석은 약 500Hz의 주파수에서 공진과 같은 진동을 예측했습니다. Maruo와 Hirata (1993)는 완전 관통 아크 용접에서 용융 풀을 모델링했습니다 [14] . 그들은 녹은 웅덩이가 정적 타원 모양을 가지고 있다고 가정했습니다. 그들은 고유 진동수와 진동 모드 사이의 관계를 조사하고 용융 풀 크기가 감소함에 따라 고유 진동수가 증가한다는 것을 보여주었습니다. Klein et al. (1994)는 원통형 키홀 모양을 사용하여 완전 침투 레이저 용접에서 키홀 진동을 연구했습니다 [15] . 그들은 점성 감쇠로 인해 키홀 진동이 낮은 kHz 범위로 제한된다는 것을 보여주었습니다. Klein et al. (1996)은 또한 레이저 출력의 작은 변동이 강한 키홀 진동으로 이어질 수 있음을 보여주었습니다[16] . 그들은 키홀 진동의 주요 공진 주파수 범위가 500 ~ 3500Hz라는 것을 발견했습니다. Andersen et al. (1997)은 고정 가스 텅스텐 아크 용접 [17] 에서 고정 된 원통형 모양을 가정하여 용융 풀의 고유 진동수를 예측 했으며 완전 용입 용접에서 용융 풀 폭이 증가함에 따라 감소하는 것으로 나타났습니다. 3.175mm 두께의 강판의 경우 주파수는 20Hz ~ 100Hz 범위였습니다. 위에 표시된 분석 방법은 일반적으로 단순한 용융 풀 모양을 가정하고 고유 진동수를 계산했습니다. 이것은 단순한 용융 풀 모양으로 고정 용접 공정을 분석하는 데 충분하지만 대부분의 용접 사례를 설명하는 과도 용접 공정에서 용융 풀 역학 분석에는 적합하지 않습니다.
반면에 수치 접근 방식은 고온 및 강한 빛과 같은 실험적 제한없이 자세한 정보를 제공하기 때문에 용융 풀 역학을 분석하는 이점이 있습니다. 전산 유체 역학 (CFD)의 수치 시뮬레이션 기술이 발전함에 따라 용융 풀 역학 분석에 대한 많은 연구가 수행되었습니다. 실제 용융 표면 변화는 VOF (체적 부피) 방법을 사용하여 계산할 수 있습니다. Cho et al. (2010) CO 2 레이저-아크 하이브리드 용접 공정을 위한 수학적 모델 개발 [18], 구형 방울이 생성 된 금속 와이어의 용융 과정이 와이어 공급 속도와 일치한다고 가정합니다. 그들은 필러 와이어가 희석되는 용융 풀 동작을 보여주었습니다. Cho et al. (2012)는 높은 빔 품질과 높은 금속 흡수율로 인해 업계에서 널리 사용되는 디스크 레이저 키홀 용접으로 수학적 모델을 확장했습니다 [19] . 그들은 열쇠 구멍에서 레이저 광선 번들의 다중 반사를 고려하고 용융 풀에서 keyholing과 같은 빠른 표면 변화를 자세히보고했습니다. 최근 CFD 시뮬레이션은 험핑 (Otto et al., 2016 [20] ) 및 기공 (Lin et al., 2017 [21] )과 같은보다 구체적인 현상을 분석하는데도 사용되었습니다 .) 레이저 용접에서. 그러나 용융 풀 역학과 관련된 연구는 거의 수행되지 않았습니다. Ko et al. (2000)은 수치 시뮬레이션을 사용하여 가스 텅스텐 아크 용접 풀의 동적 거동을 조사했습니다 [22] . 그들은 완전히 침투 된 용융 풀이 부분적으로 침투 된 풀보다 낮은 주파수에서 진동한다는 것을 보여주었습니다. 진동은 수십 분의 1 초 내에 무시할 수있는 크기로 감쇠되었습니다. Geiger et al. (2009)는 또한 수치 시뮬레이션을 사용하여 레이저 용접에서 용융 풀 거동을 보여주었습니다 [7]. 그들은 계산 된 증발 속도를 주파수 분석에 사용하여 공정에서 나오는 빛의 실험 결과와 비교했습니다. 판금 레이저 용접에서 중요한 공간 빔 진동 및 추가 필러 재료가있는 공정에 대한 용융 풀 역학에 대한 연구도 불충분합니다. Hu et al. (2018)은 금속 전달 메커니즘을 밝히기 위해 전자빔 3D 프린팅에서 와이어 공급 모델링을 수행했습니다. 그들은 주로 열 입력에 의해 결정되는 액체 브리지 전이, 액적 전이 및 중간 전이의 세 가지 유형의 금속 전달 모드를 보여주었습니다 .. Meng et al. (2020)은 레이저 빔 용접에서 용융 풀에 필러 와이어에 의해 추가 된 추가 요소의 전자기 교반 효과를 모델링했습니다. 용가재의 연속적인 액체 브릿지 이동이 가정되었고, 그 결과 전자기 교반의 영향이 키홀 깊이에 미미한 반면 필러 와이어 혼합을 향상 시켰습니다 [24] . Cho et al. (2017) 용접 방향에 수직 인 1 차원 빔 진동과 용접 라인을 따라 공급되는 필러 와이어를 사용하여 레이저 용접을위한 시뮬레이션 모델 개발 [25]. 그들은 시뮬레이션을 사용하여 특정 용접 현상, 즉 용융 풀의 단추 구멍 형성을 보여주었습니다. Cho et al. (2018)은 다중 반사 수와 전력 흡수량의 푸리에 변환을 사용하여 주파수 영역에서 소위 쵸핑 주파수 (2 x 빔 발진 주파수) 성분을 발견했습니다 [26] . 그러나 그들은 용융 풀 역학을 분석하기 위해 간접 신호를 사용했습니다. 따라서보다 직관적 인 분석을 위해서는 표면의 변동을 직접 측정해야합니다.
이 연구는 이전 연구에서 개발 된 레이저 용접 모델을 사용하여 3 차원 과도 CFD 시뮬레이션을 수행하여 빔 진동 및 필러 와이어 공급을 포함한 레이저 용접 공정에서 용융 풀 역학을 조사합니다. 용융 된 풀 표면의 시간적 변화는 시뮬레이션 결과에서 추출되었습니다. 추출 된 데이터는 주파수 영역뿐만 아니라 시간-주파수 영역에서도 분석되었습니다. 신호 처리를 통해 도출 된 결과는 특징적인 용융 풀 역학을 나타내며 빔 진동 주파수 및 단추 구멍 형성 측면에서 레이저 용접의 역학을 줄일 수있는 잠재력을 제공합니다.
2 . 방법론
그림 1도 1은 용접 방향에 수직 인 1 차원 빔 진동과 용접 라인을 따라 공급되는 필러 와이어를 사용하는 레이저 용접 프로세스의 개략적 설명을 보여줍니다. 1mm 두께의 알루미늄 합금 (AlSi1MgMn) 시트는 시트 표면에 초점을 맞춘 멀티 kW 파이버 레이저 (YLR-8000S, IPG Photonics, USA)를 사용하여 용접되었습니다. 시트는 에어 갭이있는 맞대기 이음으로 정렬되었습니다. 1 차원 스캐너 (ILV DC-Scanner, Ingenieurbüro für Lasertechnik + Verschleiss-Schutz (ILV), 독일)를 사용하여 레이저 빔의 1 차원 정현파 진동을 실현했습니다. 이 스캔 시스템에서 최대 진동 폭은 250Hz의 진동 주파수에서 1.4mm입니다. 오정렬에 대한 공차를 개선하기 위해 동일한 최대 너비 값이 사용되었습니다. 와이어 공급 시스템은 1을 공급했습니다. 2mm 직경의 알루미늄 합금 (AlSi5) 필러 와이어를 일정한 공급 속도로 에어 갭을 채 웁니다. 1mm 에어 갭의 경우 와이어 이송 속도는 용접 속도의 1.5 배 값으로 설정되었으며 참조 실험 조건은 문헌에서 얻었습니다 (Schultz, 2015 참조).[27] ).
CFD 시뮬레이션은 레이저 용접에서 열 전달 및 용융 풀 동작을 계산하기 위해 수행되었습니다. 그림 2 는 CFD 시뮬레이션을위한 계산 영역을 보여줍니다. 실온에서 1.2mm 직경의 필러 와이어가 공급되고 레이저 빔이 진동했습니다. 1mm 두께의 공작물이 용접 속도로 왼쪽에서 오른쪽으로 이동했습니다. 0.1mm의 최소 메쉬 크기가 도메인에서 생성되었습니다. 침투 깊이가 더 깊은 이전 연구의 메쉬 테스트 결과는 0.2mm 이하의 메쉬 크기로 시뮬레이션 정확도가 확보 된 것으로 나타 났으므로 [28] 본 연구에서 사용 된 메쉬 크기가 적절할 수 있습니다. 도메인을 구성하는 세포의 수는 약 120 만 개였습니다. 1 번 테이블사용 된 레이저 용접 매개 변수를 보여줍니다. 용융 풀 역학 측면에서 다양한 진동 주파수와 에어 갭 크기가 고려되었으며 12 개의 용접 사례가 표 2 에 나와 있습니다. 표 3 은 시뮬레이션에 사용 된 알루미늄 합금과 순수 알루미늄 (Cho et al., 2018 [26] )의 표면 장력 계수를 제외하고 온도와 무관 한 열-물리적 재료 특성을 보여줍니다 . 여기서 표면 장력 계수는 액체 온도에서 온도와 표면 장력 계수 사이의 선형 관계를 가진 유일한 온도 의존적 특성이었습니다.
시뮬레이션을 위해 단상 뉴턴 유체와 비압축성 층류가 가정되었습니다. 질량, 운동량 및 에너지 보존의 지배 방정식을 해결하여 계산 영역에서 속도, 압력 및 온도 분포를 얻었습니다. VOF 방법은 자유 표면 경계를 찾는 데 사용되었습니다. 스칼라 보존 방정식을 추가로 도입하여 용융 풀에서 충전재의 부피 분율을 계산했습니다. 시뮬레이션에 사용 된 레이저 용접의 수학적 모델은 다음과 같습니다. 레이저 빔은 가우스와 같은 전력 밀도 분포를 기반으로 697 개의 광선 에너지 번들로 나뉩니다. 광선 추적 방법을 사용하여 다중 반사를 고려했습니다. 재료에 대한 레이저 빔의 반사 (또는 흡수) 에너지는 프레 넬 반사 모델을 사용하여 계산되었습니다. 온도에 따른 흡수율의 변화를 고려 하였다. 혼합물의 흡수율은베이스 및 충전제 물질 분획의 가중 평균을 사용하여 계산되었습니다. 반동 압력과 부력도 고려되었습니다. 경계 조건으로 에너지와 압력의 균형은 VOF 방법으로 계산 된 자유 표면에서 고려되었습니다. 레이저 용접 모델과 지배 방정식은 FLOW-3D v.11.2 (2017), Flow Science, Inc.에서 유한 차분 방법과 유한 체적 방법을 사용하여 이산화되고 해결되었습니다. 경계 조건으로 에너지와 압력의 균형은 VOF 방법으로 계산 된 자유 표면에서 고려되었습니다. 레이저 용접 모델과 지배 방정식은 FLOW-3D v.11.2 (2017), Flow Science, Inc.에서 유한 차분 방법과 유한 체적 방법을 사용하여 이산화되고 해결되었습니다. 경계 조건으로 에너지와 압력의 균형은 VOF 방법으로 계산 된 자유 표면에서 고려되었습니다. 레이저 용접 모델과 지배 방정식은 FLOW-3D v.11.2 (2017), Flow Science, Inc.에서 유한 차분 방법과 유한 체적 방법을 사용하여 이산화되고 해결되었습니다.[29] . 계산에는 48GB RAM이 장착 된 Intel® Xeon® 프로세서 E5649로 구성된 워크 스테이션이 사용되었습니다. 계산 시스템을 사용하여 0.2 초 레이저 용접을 시뮬레이션하는 데 약 18 시간이 걸렸습니다. 지배 방정식 (Cho and Woizeschke, 2020 [30] ) 및 레이저 용접 모델 (Cho et al., 2018 [26] )에 대한 자세한 설명은 부록 A 에서 확인할 수 있습니다 .
그림 3 은 용융 풀 변동의 직접 측정에 대한 개략적 설명을 보여줍니다. 용융 풀의 역학을 분석하기 위해 시뮬레이션 중에 용융 풀 표면의 시간적 변동 운동을 측정했습니다. 상단 및 하단 표면 모두에서 10kHz의 샘플링 주파수로 변동을 측정 한 반면, 측정 위치는 X 축의 레이저 빔 위치에서 2mm 떨어진 용접 중심선에있었습니다. 그림 4시간 신호를 분석하는 데 사용되는 푸리에 변환 및 웨이블릿 변환의 개략적 설명을 보여줍니다. 측정 된 시간 신호는 고속 푸리에 변환 (FFT) 방법을 사용하여 주파수 영역으로 변환되었습니다. 결과는 측정 기간 동안 평균화 된 주파수 성분의 크기를 보여줍니다. 웨이블릿 변환 방법은 시간-주파수 영역에서 국부적 인 특성을 찾는 데 사용되었습니다. 결과는 주파수 구성 요소의 크기뿐만 아니라 시간 변화도 보여줍니다.
3 . 결과
이 연구 에서는 표 2에 표시된 12 가지 용접 사례 를 시뮬레이션했습니다. 그림 5 는 3 차원 시뮬레이션 결과를 평면도 와 바닥면으로 보여줍니다. 결과는 용융 된 풀의 거동에 따라 분류 할 수 있습니다 : 단추 구멍 형성 없음 (녹색), 안정 또는 불안정 단추 구멍 있음 (파란색), 불안정한 단추 구멍으로 인한 구멍 결함 (빨간색). 일반적인 열쇠 구멍보다 훨씬 큰 직경을 가진 단추 구멍은 레이저 용접의 특정 진동 조건에서 나타날 수 있습니다 (Vollertsen, 2016 [31]). 진동 주파수가 증가함에 따라 용접 이음 부 코스 및 스케일링 측면에서 시각적 이음새 품질이 향상되었습니다. 고주파에서 스케일링은 무시할 수있을 정도 였고 코스는 균질했습니다. 언더컷 결함의 발생도 감소했습니다. 그러나 관통 결함 부족 (case 7, case 10)이 나타났다. 에어 갭은 단추 구멍 형성에 중요했습니다. 에어 갭 크기가 증가함에 따라 단추 구멍이 더 쉽게 형성되었지만 구멍 결함으로 더 쉽게 남아 있습니다. 안정적인 단추 구멍 형성은 고려 된 공극 조건의 좁은 영역에서만 나타납니다.
그림 6 은 시뮬레이션과 실험에서 융합 영역의 모양을 보여줍니다. 버튼 홀이없는 경우 1, 불안정한 버튼 홀 형성이있는 경우 8, 안정적인 버튼 홀 형성이있는 경우 11의 3 가지 경우에 대해 시뮬레이션 결과와 실험 결과를 비교하여 유사성을 나타냈다. 본 연구에서 고려한 용접 조건의 경우 표면 품질 결과는 Fig. 5 와 같이 큰 차이를 보였으 나 단면 융착 영역 [26] 과 형상은 큰 차이를 보이지 않았다.
무화과. 7 과 8 은 각각 100Hz와 250Hz의 진동 주파수에서 시뮬레이션 결과를 기반으로 분석 된 용융 풀 역학과 시뮬레이션 및 실험 결과를 보여줍니다. 이전 연구에서 볼 수 있듯이 레이저 빔의 진동 주파수는 단추 구멍 형성과 밀접한 관련이 있습니다 (Cho et al., 2018 [26] 참조 ). 그림 7 (a) 및 (b)는 각각 시뮬레이션 및 실험을 기반으로 한 진동 주파수 100Hz에서 대표적인 용융 풀 동작을 보여줍니다. 완전히 관통 된 키홀 및 버튼 홀 형성은 관찰되지 않았으며 응고 후 거친 비드 표면이 남았습니다. 그림 7(c)와 (d)는 각각 윗면과 바닥면의 표면 변동에 대한 시뮬레이션 결과를 기반으로 한 용융 풀 역학 분석을 보여줍니다. 샘플링 데이터는 상단 표면이 공작물의 상단 표면 위치에서 평균적으로 변동하는 반면 하단 표면은 공작물의 하단 표면 위치에서 평균적으로 변동하는 것으로 나타났습니다. 표면 변동의 푸리에 변환 및 웨이블릿 변환 결과는 명확한 1 차 주파수 (2 x 빔 발진 주파수, 이른바 초핑 주파수, Cho et al., 2018 [26] 참조 ) 및 2 차 주파수 (4 x 빔 발진)를 보여줍니다. 주파수) 두 표면의 구성 요소, 그러나 바닥 표면과 첫 번째에 대한 결과주파수 성분이 더 강합니다. 반면 그림 8 (a)와 (b)에서 보는 바와 같이 250Hz의 진동 주파수에서 시뮬레이션과 실험 결과는 안정된 버튼 홀 형성과 응고 후 매끄러운 비드 표면을 나타냈다. 그림 8 의 샘플링 신호의 진폭은 그림 7 의 진폭 보다 작으며 푸리에 변환 및 웨이블릿 변환의 결과에서 중요한 주파수 성분이 발견되지 않았습니다.
Fig. 9 는 진동 주파수 200Hz에서 시뮬레이션 결과를 바탕으로 분석 된 용융 풀 역학과 시뮬레이션 및 실험 결과를 보여준다. 이 주파수에서 Fig. 9 (a)와 (b) 에서 보는 바와 같이 , 시뮬레이션과 실험 모두에서 불안정한 buttonhole 거동이 관찰되었다. 바닥면에서 샘플링 데이터의 푸리에 변환 및 웨이블릿 변환의 결과 빔 발진 주파수 성분이 발견되었습니다.
4 . 토론
시뮬레이션 및 실험 결과는 비드 표면 품질이 향상되고 빔 진동 주파수가 증가함에 따라 버튼 홀이 형성되는 것으로 나타났습니다. 표면의 변동 데이터에 대한 푸리에 변환 및 웨이블릿 변환의 결과에 따라 다음과 같은 주요 주파수 구성 요소가 발견되었습니다. 1 차 및 2 차버튼 홀 형성이없는 주파수, 불안정한 용융 풀 거동이있는 빔 진동 주파수, 안정적인 버튼 홀 형성이있는 중요한 주파수 성분이 없습니다. 이들 중 불안정한 용융 풀 동작과 관련된 빔 진동 주파수 성분은 완전히 관통 된 키홀과 반복적으로 생성 및 붕괴되는 불안정한 버튼 홀의 특성으로 인해 웨이블릿 변환 결과에서 명확한 실선 형태로 나타나지 않았습니다. 분석 결과는 윗면보다 바닥면에서 더 분명했습니다. 이는 필러 와이어 공급 및 키홀 링 공정에서 강한 하향 흐름으로 인해 용융 풀 역학이 바닥 표면 영역에서 더 강했기 때문입니다. 진동 주파수가 증가함에 따라 용융 풀 역학과 상단 표면과 하단 표면 간의 차이가 감소했습니다.
첫 번째 주파수 (2 x 빔 진동 주파수)는이 연구에서 관찰 된 가장 분명한 구성 요소였습니다. Schultz et al. (2018)은 또한 실험을 통해 동일한 성분을 발견했습니다 [32] , 용융 풀 표면 운동에 대한 푸리에 분석을 수행했습니다. 첫 번째 주파수 성분은 빔 발진주기 당 두 개의 주요 이벤트가 있음을 의미합니다. 이것은 레이저 빔이 빔 진동주기 당 두 번 와이어를 절단하거나 절단하는 프로세스와 일치합니다. 용융 된 와이어 팁은 낮은 진동 주파수에서 고르지 않고 날카로운 모서리를 갖는 것으로 나타났습니다 (Cho et al., 2018 [26] ). 이것은 첫 번째 원인이 될 수 있습니다.용융 된 풀에서 지배적이되는 주파수 성분. 진동 주파수가 증가하면 용융 된 와이어 팁이 더 균일 해 지므로 효과가 감소합니다. 용접 방향으로의 정현파 횡 방향 빔 진동을 통한 에너지 집중도 빔 진동주기 당 두 번 발생합니다. 그림 10 은 발진 주파수에 따른 레이저 빔의 라인 에너지 (단위 길이 당 에너지)의 변화를 보여줍니다. 그림 10 b) 의 라인 에너지 는 레이저 출력을 공정 속도로 나누어 계산했습니다. 여기서 처리 속도는(w이자형엘디나는엔지에스피이자형이자형디)2+(디(에스나는엔유에스영형나는디ㅏ엘wㅏV이자형나는엔에프나는지.10ㅏ))디티)2. 낮은 발진 주파수에서 라인 에너지는 발진 폭의 양쪽 끝에 과도하게 집중됩니다. 이러한 집중된 에너지는 과도한 키홀 링 프로세스를 초래하므로 언더컷 결함이 나타날 수있는 높은 흐름 역학이 발생합니다. 진동 주파수가 증가함에 따라 집중 에너지는 더 작은 조각으로 나뉩니다. 따라서 높은 진동 주파수에서 과도한 키홀 링 및 수반되는 언더컷 결함의 발생이 감소되었습니다. 위에서 언급 한 두 가지 현상 (불균일 한 와이어 팁과 집중된 라인 에너지)은 빔 발진주기 당 두 번 발생하며 발진 주파수가 증가하면 그 효과가 감소합니다. 따라서 저주파 에서 2 차 주파수 성분 (4 x 빔 발진 주파수)이 나타나는 것은이 두 현상의 동시 작용입니다.
두 가지 현상 중 첫 번째 주파수 에 대한 주된 효과 는 집중된 라인 에너지입니다. Cho et al. (2018)은 전력 흡수 데이터를 푸리에 변환을 사용하여 분석했을 때 1 차 주파수 성분이 더 우세 해졌고, 2 차 주파수 성분은 발진 주파수가 증가함에 따라 상대적으로 약화 되었음을 보여주었습니다 [26] . 용융 된 와이어 팁은 또한 빈도가 증가함에 따라 더욱 균일 해졌습니다. 결과는 진동 주파수의 증가가 용융 풀에 대한 와이어의 영향을 제거하는 것으로 나타났습니다. 따라서 발진 주파수가 증가함에 따라 라인 에너지 집중의 영향 만 남을 수 있습니다. 그림 10 과 같이, 집중 선 에너지가 작은 조각으로 분할되기 때문에 효과도 감소하지만 최대 값이 변경되지 않았기 때문에 여전히 효과적입니다.
빔 진동 주파수 성분은 불안정한 단추 구멍 및 열쇠 구멍 붕괴를 수반하는 불안정한 용융 풀 동작과 관련이 있습니다. 언더컷 결함이있는 케이스 8 (발진 주파수 200Hz)에서 발진 주파수 성분이 관찰되었습니다. 이것은 특히 완전히 관통 된 열쇠 구멍과 불안정한 단추 구멍에서 불안정한 용융 풀 동작을 보여주었습니다. 경우 10 (진동 주파수 250Hz)의 경우 상대적으로 건강한 비드가 형성 되었으나, 도 11 (a) 와 같이 웨이블릿 변환 결과에서 t1의 시간 간격으로 진동 주파수 성분이 관찰되었다 . 이 시간 간격 t1의 용융 풀 거동은 그림 11에 나와 있습니다.(비). 완전히 관통 된 열쇠 구멍이 즉시 무너지는 것이 분명하게 관찰되었습니다. 이것은 진동 주파수 성분이 불안정한 용융 풀 거동과 밀접한 관련이 있음을 보여줍니다. 발견 된 주파수 성분으로부터 완전히 관통 된 열쇠 구멍과 같은 불안정한 용융 풀 거동을 예측할 수 있습니다. 완전히 관통 된 키홀이 반복적으로 붕괴되기 때문에 빔 진동 주파수 성분은 그림 9 (d) 와 같이 웨이블릿 변환 결과에서 명확한 실선 형태로 보이지 않습니다 .
Cho and Woizeschke (2020)에 따르면 단추 구멍 형성은 자체 지속 가능한 카테 노이드처럼 작용하기 때문에 용융 풀 역학을 감소시킬 수 있습니다 [30] . 그림 12 는 버튼 홀 형성 측면에서 t2의 시간 간격에서 용융 풀 거동의 변화를 보여줍니다. 단추 구멍은 t2의 간헐적 인 부분에만 형성되었습니다. 1st 이후이 시간 동안 웨이블릿 변환의 결과로 주파수 성분이 사라졌고, 버튼 홀 형성은 용융 풀 역학을 줄이는 데 효과적이었습니다. 따라서, 웨이블릿 변환의 결과로 주파수 성분이 지워지는 것을 관찰함으로써 버튼 홀 형성을 예측할 수있다. 이와 관련하여 웨이블릿 변환 기술은 시간에 따른 용융 풀 변화를 나타낼 수 있습니다. 이 기술은 향후 용융 풀 동작을 모니터링하는 데 사용될 수 있습니다.
5 . 결론
CFD 시뮬레이션 결과를 사용하여 빔 진동 및 필러 와이어 공급을 통한 레이저 용접에서 용융 풀 역학을 분석 할 수있었습니다. 용융 풀 표면의 변동 데이터의 푸리에 변환 및 웨이블릿 변환은 여기서 용융 풀 역학을 분석하는 데 사용되었습니다. 결과는 다음과 같은 결론으로 이어집니다.1.
1 차 주파수 (2 x 빔 발진 주파수, 이른바 초핑 주파수), 2 차 주파수 (4 x 빔 발진 주파수) 및 빔 발진 주파수 성분은 푸리에 변환 및 웨이블릿 변환 분석에서 발견 된 주요 성분이었습니다.2.
1 차 주파수와 2 차 주파수 성분 의 출현은 두 가지 사건, 즉 레이저 빔에 의한 필러 와이어의 절단 공정과 집중된 레이저 라인 에너지의 효과의 결과였습니다. 이는 빔 진동주기 당 두 번 발생했습니다. 따라서 두 번째 주파수 성분은 동시 작용으로 인해 발생했습니다. 빔 진동 주파수 성분은 불안정한 용융 풀 동작과 관련이 있습니다. 구성 요소는 열쇠 구멍과 단추 구멍의 붕괴와 함께 나타났습니다.삼.
낮은 발진 주파수에서는 1 차 주파수와 2 차 주파수 성분이 함께 나타 났지만 발진 주파수가 증가함에 따라 그 크기가 함께 감소했습니다. 집중 선 에너지는 주파수가 증가함에 따라 최대 값이 변하지 않는 반면, 잘게 잘린 선단이 평평 해져 그 효과가 사라졌기 때문에 쵸핑 프로세스보다 더 큰 영향을 미쳤습니다.4.
용융 풀 거동의 빠른 시간적 변화는 웨이블릿 변환 방법을 사용하여 분석되었습니다. 따라서이 방법은 열쇠 구멍 및 단추 구멍의 형성 및 붕괴와 같은 일시적인 용융 풀 변화를 해석하는 데 사용할 수 있습니다.
CRediT 저자 기여 성명
조원익 : 개념화, 방법론, 소프트웨어, 검증, 형식 분석, 조사, 데이터 큐 레이션, 글쓰기-원고, 글쓰기-검토 및 편집. Peer Woizeschke : 감독, 프로젝트 관리, 작문-검토 및 편집.
경쟁 관심의 선언
저자는이 논문에보고 된 작업에 영향을 미칠 수있는 경쟁적인 재정적 이해 관계 나 개인적 관계가 없다고 선언합니다.
감사의 말
이 작업은 알루미늄 합금 용접 역량 센터 (Centr-Al)에서 수행되었습니다. Deutsche Forschungsgemeinschaft (DFG, 프로젝트 번호 290705638 , “용접 풀 캐비티를 생성하여 레이저 깊은 용입 용접에서 매끄러운 이음매 표면”) 의 자금은 감사하게도 인정됩니다.
-대량 보존 방정식,(A1)∇·V→=미디엄˙에스ρ어디, V→속도 벡터입니다. ρ밀도이고 미디엄˙에스필러 와이어를 공급하여 질량 소스의 비율입니다. 단위미디엄에스단위 부피당 질량입니다. WFS (와이어 공급 속도) 및 필러 와이어의 직경과 같은 매스 소스 및 필러 와이어 조건,디w계산 영역에서 다음과 같은 관계가 있습니다.(A2)미디엄=∫미디엄에스디V=미디엄0+씨×ρ×W에프에스×π디w24×티어디, 미디엄총 질량, 미디엄0초기 총 질량, V볼륨입니다.씨단위 변환 계수입니다. 티시간입니다.
-운동량 보존 방정식,(A3)∂V→∂티+V→·∇V→=−1ρ∇피+ν∇2V→−케이V→+미디엄˙에스ρ(V에스→−V→)+지어디, 피압력입니다. ν동적 점도입니다. 케이뭉툭한 영역의 다공성 매체 모델에 대한 항력 계수, V에스→질량 소스에 대한 속도 벡터입니다. 지신체 힘으로 인한 신체 가속도입니다.
-에너지 절약 방정식,(A4)∂h∂티+V→·∇h=1ρ∇·(케이∇티)+h˙에스어디, h특정 엔탈피입니다. 케이열전도율, 티온도이고 h˙에스특정 엔탈피 소스로, Eq 의 질량 소스와 연관됩니다 . (A1) . 계산 영역의 총 에너지,이자형다음과 같이 계산됩니다.(A5)이자형=∫미디엄에스h에스디V=∫미디엄에스씨Vw티w디V어디, 씨Vw질량 원의 비열, 티w질량 소스의 온도입니다.
또한, 엔탈피 기반 연속체 모델을 사용하여 고체-액체 상 전이를 고려했습니다.
-VOF 방정식,(A6)∂에프∂티+∇·(V→에프)=에프˙에스어디, 에프유체가 차지하는 부피 분율이며 0과 1 사이의 값을 가지며 에프˙에스질량의 소스와 연결된 유체의 체적 분율의 비율 식. (A1) . 질량 공급원에 해당하는 부피 분율은 다음에 할당됩니다.에프에스.
-스칼라 보존 방정식,(A7)∂Φ∂티+∇·(V→Φ)=Φ˙에스어디, Φ필러 와이어의 스칼라 값입니다. 셀의 유체가 전적으로 필러 와이어로 구성된 경우Φ1이고 유체에 대한 필러 와이어의 부피 분율에 따라 0과 1 사이에서 변경됩니다. Φ˙에스Eq 에서 질량 소스에 연결된 스칼라 소스의 비율입니다 . (A1) . 스칼라 소스는 전적으로 필러 와이어이기 때문에 1에 할당됩니다. 확산 효과는 고려되지 않았습니다.
흡수율을 계산하기 위해 프레 넬 반사 모델을 사용했습니다. ㅏ=1−ρ씨재료의 표면 상에 도시 된 바와 같이 수학 식. (A8) 원 편광 빔의 경우.(A8)ㅏ=1−ρ씨=1−12(ρ에스+ρ피)어디,ρ에스=(엔1씨영형에스θ−피)2+큐2(엔1씨영형에스θ+피)2+큐2,ρ에스=(피−엔1에스나는엔θ티ㅏ엔θ)2+큐2(피+엔1에스나는엔θ티ㅏ엔θ)2+큐2,피2=12{[엔22−케이22−(엔1에스나는엔θ)2]2+2엔22케이22+[엔22−케이22−(엔1에스나는엔θ)2]},큐2=12{[엔22−케이22−(엔1에스나는엔θ)2]2+2엔22케이22−[엔22−케이22−(엔1에스나는엔θ)2]}.어디, 복잡한 인덱스 엔1과 케이1반사 지수와 공기의 흡수 지수이며 엔2과 케이2공작물을위한 것입니다. θ입사각입니다. 도시 된 바와 같이 수학 식. (A9)에서 , 혼합물의 흡수율은 식에서 얻은 모재 및 필러 와이어 분획의 가중 평균이됩니다 . (A7) .(A9)ㅏ미디엄나는엑스티유아르 자형이자형=Φㅏw나는아르 자형이자형+(1−Φ)ㅏ비ㅏ에스이자형어디, ㅏ비ㅏ에스이자형과 ㅏw나는아르 자형이자형각각 비금속과 필러 와이어의 흡수율입니다.
자유 표면 경계에서의 반동 압력 에이 싱은 Eq. (A10) .(A10)피아르 자형(티)≅0.54피에스ㅏ티(티)=0.54피0이자형엑스피(엘V티−티비아르 자형¯티티비)어디, 피에스ㅏ티포화 압력, 피0대기압입니다. 엘V기화의 잠열, 티비끓는 온도이고 아르 자형¯보편적 인 기체 상수입니다.
참고 문헌
D.J. Kotecki, D.L. Cheever, D.G. Howden Mechanism of ripple formation during weld solidification Weld. J., 51 (8) (1972), pp. 386s-391s Google Scholar [2] M. Zacksenhouse, D.E. Hardt Weld pool impedance identification for size measurement and control J. Dyn. Syst. Meas. Control, 105 (3) (1983), pp. 179-184 CrossRefView Record in ScopusGoogle Scholar [3] V.V. Semak, J.A. Hopkins, M.H. McCay, T.D. McCay Melt pool dynamics during laser welding J. Phys. D, 28 (1995), pp. 2443-2450 CrossRefView Record in ScopusGoogle Scholar [4] A.J.R. Aendenroomer, G. den Ouden Weld pool oscillation as a tool for penetration sensing during pulsed GTA welding Weld. J., 77 (5) (1998), pp. 181s-187s Google Scholar [5] M.J.M. Hermans, G. den Ouden Process behavior and stability in short circuit gas metal arc welding Weld. J., 78 (4) (1999), pp. 137-141 View Record in ScopusGoogle Scholar [6] B.Y.B. Yudodibroto, M.J.M. Hermans, Y. Hirata, G. den Ouden Influence of filler wire addition on weld pool oscillation during gas tungsten arc welding Sci. Technol. Weld. Join., 9 (2) (2004), pp. 163-168 View Record in ScopusGoogle Scholar [7] M. Geiger, K.-H. Leitz, H. Koch, A. Otto A 3D transient model of keyhole and melt pool dynamics in laser beam welding applied to the joining of zinc coated sheets Prod. Eng. Res. Dev., 3 (2009), pp. 127-136 CrossRefView Record in ScopusGoogle Scholar [8] C. Kägeler, M. Schmidt Frequency-based analysis of weld pool dynamics and keyhole oscillations at laser beam welding of galvanized steel sheets Phys. Procedia, 5 (2010), pp. 447-453 ArticleDownload PDFView Record in ScopusGoogle Scholar [9] Y. Shi, G. Zhang, X.J. Ma, Y.F. Gu, J.K. Huang, D. Fan Laser-vision-based measurement and analysis of weld pool oscillation frequency in GTAW-P Weld. J., 94 (2015), pp. 176s-187s Google Scholar [10] J. Volpp, F. Vollertsen Keyhole stability during laser welding—part I: modelling and evaluation Prod. Eng.-Res. Dev., 10 (2016), pp. 443-457 CrossRefView Record in ScopusGoogle Scholar [11] N. Postacioglu, P. Kapadia, J. Dowden Capillary waves on the weld pool in penetration welding with a laser J. Phys. D, 22 (1989), pp. 1050-1061 CrossRefView Record in ScopusGoogle Scholar [12] N. Postacioglu, P. Kapadia, J. Dowden Theory of the oscillations of an ellipsoidal weld pool in laser welding J. Phys. D, 24 (1991), pp. 1288-1292 CrossRefView Record in ScopusGoogle Scholar [13] J. Kroos, U. Gratzke, M. Vicanek, G. Simon Dynamic behaviour of the keyhole in laser welding J. Phys. D, 26 (1993), pp. 481-486 View Record in ScopusGoogle Scholar [14] H. Maruo, Y. Hirata Natural frequency and oscillation modes of weld pools. 1st Report: weld pool oscillation in full penetration welding of thin plate Weld. Int., 7 (8) (1993), pp. 614-619 CrossRefView Record in ScopusGoogle Scholar [15] T. Klein, M. Vicanek, J. Kroos, I. Decker, G. Simon Oscillations of the keyhole in penetration laser beam welding J. Phys. D, 27 (1994), pp. 2023-2030 CrossRefView Record in ScopusGoogle Scholar [16] T. Klein, M. Vicanek, G. Simon Forced oscillations of the keyhole in penetration laser beam welding J. Phys. D, 29 (1996), pp. 322-332 View Record in ScopusGoogle Scholar [17] K. Andersen, G.E. Cook, R.J. Barnett, A.M. Strauss Synchronous weld pool oscillation for monitoring and control IEEE Trans. Ind. Appl., 33 (2) (1997), pp. 464-471 View Record in ScopusGoogle Scholar [18] W.-I. Cho, S.-J. Na, M.-H. Cho, J.-S. Lee Numerical study of alloying element distribution in CO2 laser-GMA hybrid welding Comput. Mater. Sci., 49 (2010), pp. 792-800 ArticleDownload PDFView Record in ScopusGoogle Scholar [19] W.-I. Cho, S.-J. Na, C. Thomy, F. Vollertsen Numerical simulation of molten pool dynamics in high power disk laser welding J. Mater. Process. Technol., 212 (2012), pp. 262-275 ArticleDownload PDFView Record in ScopusGoogle Scholar [20] A. Otto, A. Patschger, M. Seiler Numerical and experimental investigations of humping phenomena in laser micro welding Phys. Procedia, 83 (2016), pp. 1415-1423 ArticleDownload PDFView Record in ScopusGoogle Scholar [21] R. Lin, H.-P. Wang, F. Lu, J. Solomon, B.E. Carlson Numerical study of keyhole dynamics and keyhole-induced porosity formation in remote laser welding of Al alloys Int. J. Heat Mass Trans., 108 (2017), pp. 244-256 ArticleDownload PDFView Record in ScopusGoogle Scholar [22] S.H. Ko, C.D. Yoo, D.F. Farson, S.K. Choi Mathematical modeling of the dynamic behavior of gas tungsten arc weld pools Metall. Mater. Trans. B., 31B (2000), pp. 1465-1473 CrossRefView Record in ScopusGoogle Scholar [23] R. Hu, X. Chen, G. Yang, S. Gong, S. Pang Metal transfer in wire feeding-based electron beam 3D printing: modes, dynamics, and transition criterion Int. J. Heat Mass Transf., 126 (2018), pp. 877-887 ArticleDownload PDFView Record in ScopusGoogle Scholar [24] X. Meng, A. Artinov, M. Bachmann, M. Rethmeier Theoretical study of influence of electromagnetic stirring on transport phenomena in wire feed laser beam welding J. Laser Appl., 32 (2020), Article 022026 CrossRefGoogle Scholar [25] W.-I. Cho, V. Schultz, F. Vollertsen Simulation of the buttonhole formation during laser welding with wire feeding and beam oscillation L. Overmeyer, U. Reisgen, A. Ostendorf, M. Schmidt (Eds.), Proceedings of the Lasers in Manufacturing, German Scientific Laser Society, Munich, Germany (2017) Google Scholar [26] W.-I. Cho, V. Schultz, P. Woizeschke Numerical study of the effect of the oscillation frequency in buttonhole welding J. Mater. Process. Technol., 261 (2018), pp. 202-212 ArticleDownload PDFView Record in ScopusGoogle Scholar [27] V. Schultz, T. Seefeld, F. Vollertsen Bridging Large Air Gaps by Laser Welding with Beam Oscillation International Conference on Application of Lasers in Manufacturing, New Delhi, India (2015), pp. 31-32 CrossRefGoogle Scholar [28] W.-I. Cho, S.-J. Na Impact of wavelengths of CO2, disk, and green lasers on fusion zone shape in laser welding of steel J. Weld. Join., 38 (3) (2020), pp. 235-240 CrossRefView Record in ScopusGoogle Scholar [29] FLOW-3D User Manual. 2017. Version 11.2.1.06, Flow Science Inc. Google Scholar [30] W.-I. Cho, P. Woizeschke Analysis of molten pool behavior with buttonhole formation in laser keyhole welding of sheet metal Int. J. Heat Mass Transf., 152 (2020), Article 119528 ArticleDownload PDFView Record in ScopusGoogle Scholar [31] F. Vollertsen Loopless production: definition and examples from joining 69th IIW Annual Assembly and International Conference, Melbourne, Australia (2016) Google Scholar [32] V. Schultz, W.-I. Cho, A. Merkel, P. Woizeschke Deep penetration laser welding with high seam surface quality due to buttonhole welding Proc. of the IIW Annual Assembly, Com. IV, Bali, Indonesia (2018) IIW-Doc. IV-1390-18
CFD (전산 유체 역학) 전문가가 필요하고 때로는 실행하는데 몇 주가 걸리는 믹싱 시뮬레이션의 시대는 오래 전입니다. 컴퓨팅 및 관련 기술의 엄청난 도약에 힘 입어 Ansys, Comsol 및 Flow Science와 같은 회사는 엔지니어의 데스크톱에 사용하기 쉬운 믹싱 시뮬레이션을 제공하고 있습니다.
“병렬화 및 고성능 컴퓨팅의 발전과 템플릿화는 비전문 화학 엔지니어에게 정확한 CFD 시뮬레이션을 제공했습니다.”라고 펜실베이니아 주 피츠버그에있는 Ansys Inc.의 수석 제품 마케팅 관리자인 Bill Kulp는 말합니다 .
흐름 개선을위한 실용적인 지침이 필요하십니까? 다운로드 화학 처리의 eHandbook을 지금 흐름 도전 싸우는 방법!
예를 들어, 회사는 휴스턴에있는 Nalco Champion과 함께 프로젝트를 시작했습니다. 이 프로젝트는 시뮬레이션 전문가가 아닌 화학 엔지니어에게 Ansys Fluent 및 ACT (분석 제어 기술) 템플릿 기반 시뮬레이션 앱에 대한 액세스 권한을 부여합니다. 새로운 화학 물질을위한 프로세스를 빠르고 효율적으로 확장합니다.
Giving Mixing Its Due
“화학 산업은 CFD와 같은 계산 도구를 사용하여 많은 것을 얻을 수 있지만 혼합 프로세스는 단순하다고 가정하기 때문에 간과되는 경우가 있습니다. 그러나 최신 수치 기법을 사용하여 우수한 성능을 달성하는 흥미로운 방법이 많이 있습니다.”라고 Flow Science Inc. , Santa Fe, NM의 CFD 엔지니어인 Ioannis Karampelas는 말합니다 .
이러한 많은 기술이 회사의 Flow-3D Multiphysics 모델링 소프트웨어 패키지와 전용 포스트 프로세서 시각화 도구 인 FlowSight에 포함되어 있습니다.
“모든 상업용 CFD 패키지는 어떤 형태의 시각화 도구와 번들로 제공되지만 FlowSight는 매우 강력하고 사용하기 쉽고 이해하기 쉽게 설계되었습니다. 예를 들어, 프로세스를 재 설계하려는 엔지니어는 다양한 설계 변경의 효과를 평가하기 위해 매우 직관적인 시각화 도구가 필요합니다.”라고 그는 설명합니다.
이 접근 방식은 실험 측정을 얻기 어려운 공정 (예 : 쉽게 측정 할 수없는 매개 변수 및 독성 물질의 존재로 인해 본질적으로 위험한 공정)을 더 잘 이해하고 최적화하는데 특히 효과적입니다.
동일한 접근 방식은 또한 믹서 관련 장비 공급 업체가 고객 요구에 맞게 제품을보다 정확하게 개발하고 맞춤화하는 데 도움이되었습니다. “이는 불필요한 프로토 타이핑 비용이나 잠재적 인 과도한 엔지니어링을 방지합니다. 두 가지 모두 일부 공급 업체의 문제였습니다.”라고 Karampelas는 말합니다.
CFD 기술 자체는 계속해서 발전하고 있습니다. 예를 들어, 수치 알고리즘의 관점에서 볼 때 구형 입자의 상호 작용이 열 전달을 적절하게 모델링하는 데 중요한 다양한 문제에 대해 이산 요소 모델링을 쉽게 적용 할 수있는 반면, LES 난류 모델은 난류 흐름 패턴을 정확하게 시뮬레이션하는 데 이상적입니다.
컴퓨팅 리소스에 대한 비용과 수요에도 불구하고 Karampelas는 난류 모델의 전체 제품군을 제공 할 수있는 것이 중요하다고 생각합니다. 특히 LES는 이미 대부분의 학계와 일부 산업 (예 : 전력 공학)에서 선택하는 방법이기 때문입니다. .
그럼에도 불구하고 CFD의 사용이 제한적이거나 비실용적 일 수있는 경우는 확실히 있습니다. 여기에는 나노 입자에서 벌크 유체 증발을 모델링하는 것과 같이 관심의 규모가 다른 규모에 따라 달라질 수있는 문제와 중요한 물리적 현상이 아직 알려지지 않았거나 제대로 이해되지 않았거나 아마도 매우 복잡한 문제 (예 : 모델링)가 포함됩니다. 음 펨바 효과”라고 Karampelas는 경고합니다.
반면에 더욱 강력한 하드웨어와 업데이트 된 수치 알고리즘의 출현은 CFD 소프트웨어를 사용하여 과다한 설계 및 최적화 문제를 해결하기위한 최적의 접근 방식이 될 것이라고 그는 믿습니다.
“복잡한 열교환 시스템 및 새로운 혼합 기술과 같이 점점 더 복잡한 공정을 모델링 할 수있는 능력은 가까운 장래에 가능할 수있는 일을 간단히 보여줍니다. 수치적 방법 사용의 주요 이점은 설계자가 상상력에 의해서만 제한되어 소규모 믹서에서 대규모 반응기 및 증류 컬럼에 이르기까지 다양한 화학 플랜트 공정을 최적화 할 수있는 길을 열어 준다는 것입니다. 실험적 또는 경험적 접근 방식은 항상 관련성이 있지만 CFD가 미래의 엔지니어를위한 선택 도구가 될 것이라고 확신합니다.”라고 그는 결론을 내립니다.
본질적으로 Lagrangian 입자는 복잡한 흐름에서 물리량을 추적하는 독특한 방법을 가지고 있습니다. 이들의 속성은 메시 해상도에 의해 덜 제한되며, 동시에 질량, 운동량 및 열 전달을 통해 유체 및 고체와 함께 매우 세부적이고 사실적으로 상호 작용할 수 있습니다. 후 처리(Post Processing) 측면에서 입자는 시각화를 향상 시킬 수 있습니다.
FLOW-3D의 Lagrangian 입자 모델
FLOW-3D의 입자 모델은 전기장 효과 및 유체 흐름과의 양방향 커플 링을 포함하여 마커에서 크기와 밀도가 다른 질량 입자로 진화했습니다. 이 모델은 공기 중의 오염 물질, 금속 함유물 및 분리기에서 포착되는 파편을 추적하는데 성공적으로 적용되었습니다. 최근에는 FLOW-3D의 입자 모델이 기능을 확장하기 위한 큰 변화가 있었습니다. 현재 모델에서 입자는 기본 기능에 따라 클래스로 그룹화됩니다.
마커 입자 는 단순한 질량이 없는 마커로 유체 흐름을 추적하는 데 가장 적합합니다.
질량 입자 는 모래 알갱이 또는 내포물과 같은 고체 물체를 나타냅니다.
액체 입자 는 유체로 만들어지며 모든 유체 속성을 상속합니다.
가스 입자 는 주변 유체의 온도 및 압력 부하에 따라 크기가 변하는 기포를 나타냅니다.
보이드 입자 는 가스 입자와 유사하지만 그 특정 기능은 붕괴된 기포를 표시하고 추적하는 것입니다. 이는 다른 응용 분야에서 주조시 금형 충전 중에 생성되는 잠재적 다공성 결함을 예측하는 데 유용합니다.
프로브 입자 는 해당 위치에서 변수 값을 기록하고 보고하는 진단 장치로 사용됩니다. 다른 클래스의 입자로 만들 수 있습니다.
사용자 입자 는 소스 코드에서 사용자 정의 함수를 통해 사용자 정의를 할 수 있습니다.
각 입자 클래스에는 드래그 계수 및 각 숫자 입자가 물리적 입자의 구름을 나타낼 수 있는 매크로 입자 계수와 같이 클래스의 모든 입자에 적용되는 속성이 있습니다. 사용자 클래스의 입자에는 사용자가 사용자 정의 할 수 있는 세 가지 추가 속성이 있습니다.
다양한 크기와 밀도의 입자를 나타내는 재료 입자 클래스 내에서 여러 종을 정의 할 수 있습니다. 주변 유체와의 열 전달은 모든 재료 입자, 즉 질량, 액체, 가스, 보이드 및 사용자 입자에 적용되는 또 다른 기능입니다.
가스 입자의 압력은 상태 방정식과 온도 변화에 따른 변화를 사용하여 계산됩니다. 기체 입자가 유체가 없는 표면을 벗어나면 기체 영역에 부피를 추가합니다.
액체 입자의 유체는 응고 뿐만 아니라 증발 및 응축으로 인해 상 변화를 겪을 수 있습니다. 응고된 입자는 질량 입자와 유사한 고체 물체로 작동하지만 일단 들어가서 다시 녹으면 유체로 변환됩니다. 또한 2 유체 상 변화 모델이 활성화되면 액체 입자가 기체 내에서 이동하면서 증발 및 응축될 수 있으므로 스프레이 냉각 모델링에 유용합니다.
각 파티클 클래스는 FLOW-3D POST 에서 별도의 개체로 시각화 할 수 있습니다. 속도, 온도, 입자 수명 또는 고유 ID와 같은 개별 입자 속성을 색상에 사용할 수 있습니다. 표시된 입자 크기는 클래스 내에서의 변화를 반영합니다.
Lagrangian 입자를 직접 금속 증착에 적용
직접 금속 증착은 동일한 금속의 분말 스트림이 주입되는 고체 금속 기판에 용융 풀을 형성하기 위해 레이저를 사용하는 적층 제조 공정의 한 유형입니다. 분말 입자가 풀 내부에서 녹고, 풀이 다시 응고되면 일반적으로 두께가 0.2-0.8mm이고 너비가 1-2mm 인 고형화된 금속 층이 형성됩니다.
laser/powder gun 어셈블리가 기판 표면을 계속 스캔하므로 복잡한 모양을 층별로 만들 수 있습니다. 레이저 출력, 속도 및 분말 공급 사이의 적절한 균형은 공정의 성공과 효율성을 위해 중요합니다. 엔지니어의 주요 관심 사항은 다음과 같습니다.
용융 풀의 크기와 모양
금속 흐름 및 그 내부의 냉각 속도
응고된 층의 형상
이 섹션에서 설명하는 시뮬레이션은 이러한 특성을 정확하게 예측합니다. 레이저와 기판의 움직임은 좌표계를 레이저에 부착함으로써 반전됩니다. Inconel 718 합금의 기판은 10mm/s의 일정 속도로 움직입니다. 레이저는 1.8kW의 출력으로 반경 1mm의 원형 열원으로 모델링됩니다. 3 개의 파우더 건은 0.684 g/s의 속도로 레이저 충돌 점에서 고체 금속 입자를 전달합니다. 각 건은 크기가 2 x 2 mm이고 초당 입자 비율은 105 입니다.
입자는 액체 입자 클래스를 사용하여 모델링됩니다. 모든 입자의 직경은 40 μm입니다. 매크로 입자 배율 10은 시뮬레이션에서 입자 수를 줄이는데 사용됩니다. 3백만 개의 물리적 입자를 나타내는 매 초당 시뮬레이션에서 3 x 105 개의 숫자 입자가 생성됩니다. 입자의 초기 온도는 480°C입니다. 즉, 풀에 충돌하기 전에 고체 상태입니다.
시뮬레이션은 분말을 첨가하기 전에 용융 풀이 형성 될 수 있도록, 시작한 후 2초 후에 입자 소스를 활성화하여 10초 동안 실행했습니다. 일단 풀에 들어가면 입자가 녹아 금속으로 전환되어 금속의 부피가 증가하여 궁극적으로 레이저에서 하류의 재응고 금속 층을 형성합니다. 용융 풀 모양은 대칭 평면에 표시됩니다.
새로운 Lagrangian 입자 모델은 FLOW-3D의 현재 기능을 크게 확장 할 뿐만 아니라 금속의 핵심 가스 버블 추적과 같은 향후 확장을 위한 강력한 개발 플랫폼을 만듭니다.
때로는 유체 입자가 이동한 거리가 중요한 경우도 있습니다. FLOW-3D에서 사용자는 모델 설정 ‣ 출력 위젯에서 유체가 이동한 거리에 대한 출력을 요청할 수 있습니다. 이 기능은 유체가 흐름 영역(경계 또는 질량 소스를 통해)에 들어간 시간 또는 유체가 도메인을 통해 이동한 거리를 계산합니다. 이 기능은 모든 시뮬레이션에도 사용할 수 있으며, 특별한 모델을 사용할 필요가 없으며, 흐름에도 영향을 미치지 않습니다. 이 모델을 사용하려면 출력 위젯으로 이동하고 추가 출력 섹션에서 “Distance traveled by fluid” 옆의 체크상자를 선택하십시오.
유체 도착 시간을 아는 것은 종종 유용합니다. 예를 들어 주조 시뮬레이션에서 주입 시간을 결정하는 데 사용할 수 있습니다. 제어 볼륨은 충전 프로세스 동안 여러 번 채워지고 비워지기 때문에 계산 셀이 채워지는 처음과 마지막 시간 모두 기록되고, 후 처리를 위해 저장될 수 있습니다. 이 작업은 출력 위젯과 추가 출력 섹션 내에서 유체 도착 시간 확인란을 선택하여 수행됩니다.
때로는 유체가 계산 영역 내에서 보내는 시간인 체류시간을 아는 것이 유용합니다. 이는 출력 ‣ Output ‣ Additional Output ‣ Fluid residence time 확인란을 선택하여 수행합니다. 여기서 S로 지정된 이 변수에 대한 전송 방정식은 단위 소스 항과 함께 Solve됩니다.
여기에서 t는 시간이며 u는 유체 속도입니다.
S의 단위는 시간이다. 계산 도메인에 들어가는 모든 유체에 대한 S의 초기값은 0입니다.
의 값은 항상 second order체계를 가진 데이터로부터 근사치를 구합니다.
이 출력 옵션은 1 유체 및 2 유체 유량 모두에 사용할 수 있습니다.
노트
경계 조건 또는 소스에서 도메인으로 유입되는 유체가 이미 도메인에 있는 유체와 혼합될 때 체류가 감소하는 것처럼 보일 수 있습니다.
벽면 접촉 시간 출력은 (1)개별 유체 요소가 특정 구성 요소와 접촉하는 시간 및 (2)특정 구성 요소가 유체와 접촉하는 시간을 추적합니다. 이 모델은 액체 금속이 모래 오염물과 접촉했을 때 오염과 상관 관계가 있는 proxy 변수를 제공하기 위한 것입니다. 이 출력은 최종 주조물에서 오염된 유체가 어디에 있는지 확인하는 데 사용될 수 있습니다. 접촉 시간 모델의 또 다른 해석은, 예를 들어, 용해를 통해 다소 일정한 비율로 화학물질을 방출하는 물에 잠긴 물체에 의한 강의 물의 오염입니다.
모델은 Model Setup ‣ Output ‣ Wall contact time 박스를 확인하여 활성화됩니다. 또한 Model Setup ‣ Output ‣ Geometry Data section의 각 구성요소에 대해 해당 구성요소를 계산에 포함하기 위해 반드시 설정해야 하는 Contact time flag가 있습니다.
추가 정보
Wall Contact Time with Fluid and Component Properties: Contact Time with Fluid for more information on the input variables를 참조하십시오.
노트
이 모델은 실제 구성 요소, 즉 고체, 다공성 매체, 코어 가스 및 충전 퇴적물 구성 요소로 제한됩니다. 접촉 시간은 유체 # 1과 관련해서만 계산됩니다.
Model Setup ‣ Output ‣ Geometry Data ‣ Forces 옵션을 활성화하면 부품에 대한 압력, 전단력, 탄성 및 벽 접착력을 History Data에 출력할 수 있습니다.
압력을 가지지 않은 셀(즉, 도메인 외부에 있거나 다른 구성 요소 안에 있는 셀)이 구성 요소 주변의 각 셀에 대한 압력 영역 제품을 합산하는 동안 어떻게 처리되는지를 제어하는 압력 계산에 대한 몇 가지 추가 옵션이 있습니다. 기본 동작은 이러한 셀에서 사용자 정의 기준 압력을 사용하는 것입니다. 지정되지 않은 경우 기준 압력은 초기 무효 압력인 PVOID로 기본 설정됩니다. 또는, 코드는 Reference pressure is code calculated 옵션을 선택하여 구성요소의 노출된 표면에 대한 평균 압력을 사용할 수 있습니다.
마지막으로, 일반 이동 물체의 경우, 규정된/제약을 받는 대로 물체를 이동시키는 힘을 나타내는 잔류 힘의 추가 출력이 있습니다.
Torques
Model Setup ‣ Output ‣ Force 옵션이 활성화되면 구성 요소의 토크가 계산되고 History Data에 출력됩니다. 토크는 힘-모멘트에 대한 기준점 X, 힘-모멘트에 대한 기준점 Y, 정지 구성 요소에 대한 힘-모멘트 입력에 대한 기준점 Z에 의해 지정된 지점에 대해 보고됩니다. 참조점의 기본 위치는 원점입니다.
General Moving Objects에는 몇 가지 추가 참고 사항이 있습니다. 첫째, 토크는 (1) 6-DOF 동작의 질량 위치 중심 또는 (2)고정축 및 고정점 회전의 회전 축/점에 대해 보고됩니다. 힘에서 행해지는 것과 마찬가지로, 규정된/제한된 바와 같이 물체를 이동시키는 토크를 나타내는 잔류 토크의 출력도 있습니다.
노트
힘 및 토크 출력은 각 지오메트리 구성 요소의 일반 히스토리 데이터에 기록됩니다. 출력은 개별 힘/토크 기여 (예: 압력, 전단, 탄성, 벽 접착) 및 개별 기여도의 합으로 계산된 총 결합력/토크로 제공됩니다.
일반 이동 객체의 부력과 안정성에 대한 정보는 각 구성 요소에 대해 모델 설정 Setup 출력 ‣ 기하학적 데이터 ‣ 부력 중심 및 도량형 높이 옵션을 활성화하여 History Data에서 출력할 수 있습니다. 이렇게 하면 구성 요소의 중심 위치와 중심 높이가 출력됩니다.
깊이 기준 유압 데이터를 요청하려면 출력 ‣ 고급으로 이동한 후 유압 데이터 옆의 확인란을 선택하십시오(심층 평균 값과 중력을 -Z 방향으로 가정).
이 옵션은 FLOW-3D가 유압 시뮬레이션에 유용할 수 있는 추가 깊이 평균 데이터를 출력하도록 합니다.
Flow depth
Maximum flow depth
Free surface elevation
Velocity
Offset velocity
Froude number
Specific hydraulic head
Total hydraulic head
이 수량 각각에 대해 하나의 값 이 메쉬의 모든 (x, y) 위치에서 계산되고 수직 열의 모든 셀에 저장됩니다 (이 수량이 깊이 평균이기 때문에 z 방향으로 데이터의 변화가 없습니다). 변수는 정확도를 보장하기 위해주기마다 계산됩니다. 모든 경우에, 깊이 평균 속도, z- 방향 의 중력 가속도, 유체 깊이, 및 컬럼 내 유체의 최소 z- 좌표입니다.
자유 표면 고도는 수직 기둥의 맨 위 유체 요소에 있는 자유 표면의 z-좌표로 계산됩니다.
The Froude number 은
식으로 계산됩니다.
유체 깊이는 깊이 평균 메쉬 열의 모든 유체의 합으로 계산됩니다.
특정 유압 헤드
및 총 유압 헤드
변수는 다음에서 계산됩니다.
노트
깊이 기준 유압 출력 옵션은 예리한 인터페이스가 있고 중력이 음의 z 방향으로 향할 때에만 유체 1에 유효합니다.
유압 헤드 계산은 스트림 라인이 평행하다고 가정한다는 점을 유념해야 합니다. 예를 들어 플럭스 표면이 재순환 흐름 영역에 배치되는 경우 이 문제가 발생할 수 있습니다. 이 경우, 유량 표면에서 보고된 유량 평균 유압 헤드는 헤드의 계산에서 흐름 방향이 무시되기 때문에 예상보다 클 수 있습니다.
Total Hydraulic Head 3D(총 유압 헤드 3D)
또한 총 유압 헤드 3D 옵션을 확인하여 국부적(3D) 속도 필드, 플럭스 표면에서의 유압 에너지(배플 참조) 및 플럭스 기반 유압 헤드를 사용하여 유체 1의 총 헤드를 계산할 수 있다. 3D 계산은 국부 압력을 사용하여 수행되며(즉, 압력이 유체 깊이와 관련이 있다고 가정하지 않음) 원통 좌표와 호환됩니다.
노트
유압 헤드 계산은 스트림 라인이 평행하다고 가정한다는 점을 유념해야 한다. 예를 들어 플럭스 표면이 재순환 흐름 영역에 배치되는 경우 문제가 발생할 수 있습니다. 이 경우, 플럭스 표면에서 보고된 유량 평균 유압 헤드는 헤드의 계산 시 흐름 방향이 무시되기 때문에 예상보다 클 수 있습니다.
3D 유압 헤드 계산은 입력 파일에 중력이 정의되지 않은 경우 중력 벡터의 크기를 1로 가정합니다.
Flux-averaged hydraulic head
특정 위치 (즉, 배플)의 플럭스 평균 유압 헤드는 다음과 같이 계산됩니다.
유압 헤드 계산에서는 유선이 평행하다고 가정합니다. 예를 들어 플럭스 표면이 재순환 흐름 영역에 배치된 경우 (예: 아래에 표시된 것과 같이) 문제가 될 수 있습니다.
이 경우 플럭스 표면에 보고된 플럭스 평균 유압 헤드는 헤드 계산 시 흐름 방향이 무시되므로 예상보다 클 수 있습니다.
FLOW-3D에는 History Probes, Flux surface, Sampling Volumes의 세 가지 주요 측정 장치가 있습니다. 이러한 장치를 시뮬레이션에 추가하는 방법은 모델 설정 섹션에 설명되어 있습니다(측정 장치 참조). 이들의 출력은 기록 데이터 편집 시간 간격으로 flsgrf 파일의 일반 기록 데이터 카탈로그에 저장됩니다. 이러한 결과는 Analyze ‣ Probe 탭에서 Probe Plots을 생성하여 액세스할 수 있습니다.
히스토리 프로브를 생성하는 단계는 모델 설정 섹션에 설명되어 있습니다(기록 프로브 참조). 시뮬레이션에 사용된 물리 모델에 따라 각각의 History Probe에서 서로 다른 출력을 사용할 수 있습니다. 프로브를 FSI/TSE로 지정하면 유한 요소 메시 안에 들어가야 하는 위치에서 응력/스트레인 데이터만 제공한다. 유체 프로브가 솔리드 형상 구성 요소에 의해 차단된 영역 내에 위치하는 경우, 기하학적 구조와 관련된 수량(예: 벽 온도)만 계산된다. 일반적으로 프로브 좌표에 의해 정의된 위치에서 이러한 양을 계산하려면 보간이 필요하다.
Flux-averaged hydraulic head when 3D Hydraulic Head is requested from additional output options
Hydraulic energy flow when hydraulic data output is requested
Total number of particles of each defined species in each particle class crossing flux surface when the particle model is active
Flow rate for all active and passive scalars this includes scalar quantities associated with active physical models (eg. suspended sediment, air entrainment, ect.)
노트
유속과 입자수의 기호는 유동 표면을 설명하는 함수의 기호에 의해 정의된 대로 흐름이나 입자가 플럭스 표면의 음에서 양으로 교차할 때 양의 부호가 됩니다.
플럭스 표면은 각 표면의 유량과 입자 수가 정확하도록 그들 사이에 적어도 두 개의 메쉬 셀이 있어야 합니다.
유압 데이터 및 총 유압 헤드 3D 옵션을 사용할 때는 유압 헤드 계산이 스트림 라인이 평행하다고 가정한다는 점을 유념해야 한다. 예를 들어 플럭스 표면이 재순환 흐름 영역에 배치되는 경우 이 문제가 발생할 수 있습니다. 이 경우, 유량 표면에서 보고된 유량 평균 유압 헤드는 헤드의 계산에서 흐름 방향이 무시되기 때문에 예상보다 클 수 있습니다.
샘플링 볼륨은 해당 범위 내에서 볼륨을 측정하는 3 차원 데이터 수집 영역입니다. 샘플링 볼륨을 만드는 단계는 모델 설정 섹션에 설명되어 있습니다(샘플링 볼륨 참조). 각 샘플링 볼륨의 계산 수량은 다음과 같습니다.
시료채취량 내에서 #1 유체 총량
시료채취량 내 #1 유체질량 중심
샘플링 용적 가장자리에 위치한 솔리드 표면을 포함하여 샘플링 용적 내의 모든 벽 경계에 작용하는 좌표계의 원점에 상대적인 유압력 및 모멘트.
샘플링 용적 내 총 스칼라 종량: 이것은 부피 적분으로 계산되므로 스칼라 양이 질량 농도를 나타내면 샘플링 용적 내의 총 질량이 계산된다. 거주 시간과 같은 일부 종의 경우, 평균 값이 대신 계산됩니다.
샘플링 볼륨 내의 입자 수: 각 샘플링 볼륨 내에 있는 각 입자 등급의 정의된 각 종별 입자 수(입자 모델이 활성화된 경우)
운동 에너지, 난류 에너지, 난류 소실율 및 와류에 대한 질량 평균
표본 체적의 6개 경계 각각에서 열 유속: 유체 대류, 유체 및 고체 성분의 전도 및 유체/구성 요소 열 전달이 포함됩니다. 각 플럭스의 기호는 좌표 방향에 의해 결정되는데, 예를 들어, 양방향의 열 플럭스도 양수입니다. 출력에서 확장 또는 최대 디버그 수준을 선택하지 않는 한 이러한 디버그 수준은 fsplt에 자동으로 표시되지 않습니다.
증발하는 빗방울에서 남은 잔류의 물은 새로 씻은 자동차에서 좋지 못할 수 있습니다. 그러나, 동일한 증발 공정은, 예를 들어, 드롭 잔류 물이 인쇄 된 이미지 또는 텍스트의 일부가되는 잉크젯 인쇄에서 유리할 수있다. 그러나 동일한 증발 과정이 어떤 경우엔 도움이 될 수 있습니다 예를 들면, 잉크 찌꺼기가 인쇄 된 이미지나 텍스트의 일부가 되는 잉크젯 인쇄가 그렇습니다.
액체 방울의 증발로 인한 잔류의 물이 예상치 못한 방식으로 나타날 수 있습니다. 커피 링 얼룩이 잘 알려진 예이며, 커피의 잔류의 물이 물방울의 바깥 쪽 가장자리에 모여 얇은 원형 링 얼룩이 남습니다. 이 현상은 흥미로운 유체역학적인 과정의 결과입니다. 커피 링 얼룩이 형성 되려면 액체가 증착 된 고체 표면에 고정 된 접촉선이 있어야합니다. 고정 된 접촉선은 액체 방울이 고체 기판과 교차하는 액체 방울의 외부의 가장자리가 방울이 증발함에 따라 정지 상태를 유지함을 의미합니다. 증발은 기판의 열에 의해 발생하며 방울의 얇은 외부의 가장자리에서 가장 크게 생깁니다. 표면 장력은 액체가 증발하면서 손실 된 액체를 대체하기 위해 가장자리를 향해 발생하게 됩니다. 이는 결국 더 많은 용질을 가장자리로 운반하며 모든 액체가 증발 한 후, 결과적으로 커피 링 얼룩을 형성하게하는 더 높은 농도의 용질 잔류 물을 생성합니다.
모델링 접근법
FLOW-3D v12.0의 최신 업데이트로 인해 ‘접촉선의 고정’ 모델이 개발되었으며, 소프트웨어의 기능이 표면 장력 중심의 애플리케이션으로도 광범위하게 확장되었습니다. 표면 접촉의 고정 및 비고정 특성은 잉크젯 인쇄, 코팅 및 스프레이 냉각에서 중요한 역할을 합니다. 습윤 특성에 대한 표면 공법은 미세 유체 장치에서 액체 샘플의 이동을 제어하는 데 사용될 수 있습니다. 모델의 주요 특징은 방울의 가장자리를 고정 위치에 고정하는 수단을 제공하는 것입니다. 형상 구성 요소 및 하위 구성 요소중에 표면에 ‘고정’ 속성을 지정할 수 있습니다. 유체의 접촉선은 처음 표면과 접촉하는 곳에 고정됩니다. 전방 속도를 0으로 유지하면 고정이 적용됩니다. 유체는 접촉선과 표면을 따라 이동하는 것이 아니라 롤오버하여 접촉점을 지나야만 이동할 수 있습니다.
커피 링 얼룩 검증
그림 1은 평평한 수평 표면에 놓인 원형 물방울의 결과를 보여줍니다. 표면은 30 ℃의 일정한 온도로 유지됩니다. 초기 유체 온도는 20 ℃이고 주변 공극의 온도는 일정한 20 ℃입니다. 유체는 밀도 0.967 g/cm3, 점도 0.02022 poise, 비열 1.645e+07 cm2/s/K, 열전도도 1.2964e+4 g*cm/s3/K, 표면 장력 계수 33.15 g/cm2의 일반적인 잉크를 나타냅니다.
액적 표면의 초기 곡률 반경은 7.5e-03 cm이고, 차지하는 공간은 반경 4.5e-03 cm의 원이며, 겉보기의 초기 접촉각은 37.87 도입니다. 그림 1-a를 참조하시기 바랍니다. 지정된 정적 접촉각은 0 도입니다.
정압에 의한 상변화 모델이 활성화됩니다. 공극 내의 증기 분압은 0이고 상변화 수용 계수는 Rsize = 0.01 입니다.
잉크가 건조될 때 기판 상에 고체가 잔류하는 물이 형성되는 것을 포착하기 위해 잔류 물 모델도 켜집니다. 유체에 용해 된 안료의 농도는 초기 농도 0.01 g/cm3 이고 최대 농도 rmax = 1.1625 g/cm3 에서 운반이 가능한 스칼라로 표시됩니다. 용해 된 안료는 질량 평균을 기준으로 안료의 단위질량당 0.05 poise의 속도로 유체의 순 점도를 향상시킵니다.
이 공정은 3.0 도의 방위 방향으로 하나의 셀에 걸쳐있는 축 대칭 원통형 메쉬로 모델링됩니다. (x 간격 = 6e-05 cm, z 간격 = 4e-05 cm.)
그림 1은 유체가 증발함에 따라 접촉선이 고정 된 상태를 유지하고 있음을 보여줍니다. 0 도의 정적 접촉각 조건은 액적의 중심을 향한 압력 구배를 가져오고, 이는 접촉선 방향으로의 유동을 생성합니다. 용해 된 안료의 농도는 증발로 인해 자유 표면 근처에서 증가하며, 흐름을 따라 농도는 접촉선을 향해 더욱 재분배합니다. (그림 2). 액체가 계속 증발함에 따라, 남아있는 액체의 안료 농도는 증가합니다. 농도가 최대 rmax에 도달하면, 과잉된 안료는 고체가 잔류하는 물로 전환됩니다.
접촉선 근처의 유체가 먼저 건조되어 고체가 잔류하는 물이 남습니다. 해당 영역의 유체에 안료 농도가 높기 때문에 고체가 잔류하는 물의 특징인 ‘커피 링’ 패턴이 기판 표면에 생성됩니다. (그림 3 및 4). 안료의 총 질량(용해 + 건조 잔류 물)은 초기 질량의 0.025 % 이내로 보존됩니다.
물방울 벽의 검증
접촉선 고정 응용의 두 번째 예는 수직의 벽에 고정 된 한 방울의 액체 알루미늄의 거동입니다. 유체 밀도는 2.7 g / cm3, 표면 장력 계수 200 g / cm2 및 점도 0.27 poise입니다. 정적 접촉각은 0 도입니다.
초기의 겉보기의 접촉각이 90도가 되도록 반경 0.5cm의 물방울을 수직 벽에 놓습니다 (그림 5). 7e+06 cm/s2의 중력 크기는 표면 장력의 복원 작용을 없애고 액적이 눈에 띄도록 변형시키기 위하여 인위적으로 향상되었습니다. 결과들은 비슷한 크기의 물방울에 대한 실험 결과와의 질적 비교를 포함하여 그림 5에서 보여줍니다.
요약
FLOW-3D의 접촉선 고정 모델은 표면 장력 및 벽의 접착 기능을 확장하여 표면 공법에서 복잡한 상호 작용을 모델링합니다. 접촉선 고정이 실제로 응용되는 분야에 관하여 더 많은 예시와 추가적인 참조를 찾으신다면여기에서 찾을 수 있습니다.
FLOW–3D의 침전물 제거 및 드리프트 플럭스 모델은 침전 탱크, 이차 클라리퍼 및 기타 유사한 구조물을 따라 순환하는 입자의 정착 과정을 모델링 하는 데 매우 유용한 도구입니다. 이러한 모델을 사용하면 사용자가 입자 지름, 드래그 계수 등과 같은 개별 특성을 사용하여 모델에 여러개의 구성을 입력할 수 있습니다.
표준 출력물에는 유체 내 침전물 농도와 포장된 침전물 축적 및 가까운 기하학의 단단한 표면에 누적되는 비율이 포함됩니다. 전체적으로 시스템의 추가적인 유압 특성(예:전체적인 유압 효율과 잔류 시간)은 스칼라 추적 모델의 동시 사용을 통해 쉽게 얻을 수 있다. 이동 물체(밸브나 게이트 개방의 작용 등)를 포함하여 단순하고 효율적으로 솔리드 객체를 표현하는 FLOW-3DFAVOR™방법을 사용하면 재메싱 작업에 소요되는 시간이 최소화된 기하학적 매개 변수 공간을 매우 효율적으로 연구할 수 있습니다.
FLOW-3D에서 벽, 다공성 벽 또는 플럭스 측정 표면의 역할을 할 수 있는 완화된 객체는 이러한 유형의 시뮬레이션에 사용할 수 있는 모델링 접근 방식에 추가적인 유연성을 제공합니다.
이 튜토리얼 매뉴얼은 FLOW-3D 를 처음 사용하는 사용자에게 그래픽 사용자 인터페이스(GUI)의 주요 구성 요소를 쉽게 익히도록 하고, 다양한 시뮬레이션의 설정 및 실행 방법을 안내하기 위한 것입니다.
이 매뉴얼에 있는 실습과정은FLOW-3D의 기본 사항을 다루기 위한 것입니다. 이 매뉴얼에서 제시하는 문제는 다양한 주제를 설명하고, 발생할 수 있는 많은 질문을 해결하기 위해 선정되었습니다. 이 매뉴얼의 실습과정은FLOW-3D를 실행하는 컴퓨터에 앉아 사용하는 것이 가장 좋습니다.
CFD 사용 철학에 대한 간단한 섹션 다음에는 중요 파일과 시뮬레이션 파일을 실행하는 방법이 소개되어 있습니다. 이 소개 섹션 다음에는 모델 설정, 시뮬레이션 실행 및 포스트 프로세스, Simulation Manager 탐색 방법에 대한 설명이 있습니다. 이러한 각 단계에 대한 자세한 내용은 모델 설정, 컴퓨팅 결과 및 후처리 장에서 확인할 수 있습니다.
1.CFD 사용에 대한 철학
CFD (Computational Fluid Dynamics)는 유체 흐름(질량, 운동량 및 에너지 보존)에 대한 지배 방정식의 컴퓨터 솔루션입니다. 지정된 지배방정식은 이론 장에 설명된 Numerical방법을 사용하여 이산화되고 계산됩니다.
CFD 소프트웨어를 사용하는 것은 여러 면에서 실험을 설정하는 것과 유사합니다. 실제 상황을 시뮬레이션하기 위해 실험을 올바르게 설정하지 않으면, 그 결과는 실제 상황을 반영하지 않습니다. 같은 방법으로 수치 모델이 실제 상황을 정확하게 나타내지 않으면, 그 결과는 실제 상황을 반영하지 않습니다. 사용자는 어떤 것이 중요한지, 어떻게 표현해야 하는지를 결정해야 합니다. 시작하기 전에 다음과 같은 질문을 하는 것이 중요합니다.
CFD 계산에서 무엇을 알고 싶습니까?
중요한 현상을 포착하기 위해 규모와 Mesh는 어떻게 설계되어야 하는가?
실제 물리적 상황을 가장 잘 나타내는 경계 조건은 무엇입니까?
어떤 종류의 유체를 사용해야합니까?
이 문제에 어떤 유체 특성이 중요합니까?
다른 어떤 물리적 현상이 중요합니까?
초기 유체 상태는 어떻게 됩니까?
어떤 단위 시스템을 사용해야합니까?
모델링 되는 문제가 실제 상황을 가능한 한 유사하게 나타내는지 확인하는 것이 중요합니다. 사용자는 복잡한 시뮬레이션 작업을 해결 가능한 부분으로 나누는 것이 좋습니다.
복잡한 물리 효과를 추가하기 전에, 간단하고 쉽게 이해할 수 있는 근사값으로 점차적으로 시작하여 프로세스 진행하십시오. 간단한 손 계산(베르누이 방정식, 에너지 균형, 파동 전파, 경계층 성장 등)은 물리 및 매개 변수를 선택하는데 도움이 되고, 결과와 비교할 수 있는 점검항목을 제공합니다.
CFD의 장단점을 이해하면 분석을 진행하는데 도움이 될 수 있습니다. CFD는 다음과 같은 경우 탁월한 분석 옵션입니다.
기하 구조, 물리학 또는 필요한 상세 수준으로 인해 표준 엔지니어링 계산이 유용하지 않은 경우가 많습니다.
실제 실험은 비용이 많이 소요됩니다.
실험에서 수집할 수 있는 것보다 유체흐름에 대한 자세한 정보가 필요한 경우 유용합니다.
위험하거나 적대적인 조건, 확장이 잘되지 않는 프로세스 등으로 인해 정확한 실험 측정을 하기가 어려운 경우
복잡한 흐름 정보에 대한 커뮤니케이션
CFD는 다음과 같은 경우에 덜 효과적입니다.
솔루션이 계산 리소스가 매우 많이 소요되거나, 도메인 크기를 줄이기 위한 가정 또는 해결되지 않은 물리적 현상을 설명하기 위한 반 임계 모델이 필요한 경우
CFD 시뮬레이션에 대한 입력이 되는 중요한 물리적 현상이 알려지지 않은 경우
물리적 현상이 잘 이해되지 않거나 매우 복잡한 경우
CFD를 사용할 때 명심해야 할 몇 가지 중요한 참고 사항이 있습니다.
CFD는 규정된 초기 및 경계 조건에 따라 지정된 지배 방정식의 수치해석 솔루션입니다. 따라서 모델 설정, 즉 어떤 방정식을 풀어야 하는지, 재료 특성, 초기 조건 및 경계 조건이, 가능한 한 물리적 상황과 최대한 일치해야 합니다.
방정식의 수치 해는 일반적으로 어떤 종류의 근사치를 필요로 합니다. 물리적 모델에 대한 가정과 해결방법을 검토한 후 사용하는 것이 좋습니다.
디지털 컴퓨터는 숫자가 유한 정밀도로 이진수로 표시되는 방식으로 인해 반올림 오류가 발생합니다. 이는 문제를 악화시키기 때문에 매우 근소한 숫자의 차이를 계산해야 하는 상황을 피하십시오. 이러한 상황의 예는 시뮬레이션 도메인이 원점에서 멀리 떨어져 있을 때입니다.
2.중요한 파일
FLOW-3D시뮬레이션과 관련된 많은 파일이 있습니다. 가장 중요한 것들이 아래에 설명되어 있습니다. 모든 prepin.*파일의 명칭에서 prepin는 파일 형식을 의미하며, 별표시*위치는 시뮬레이션 이름을 의미합니다. (예 : prepin.example_simulation.)
·prepin.*: 시뮬레이션용 입력 파일입니다. 시뮬레이션 설정을 설명하는 모든 입력 변수가 포함되어 있습니다.
·prpgrf.*: 이것은 전 처리기 출력 파일입니다. 여기에는 계산된 초기 조건이 포함되며 시뮬레이션을 실행하기 전에 설정을 확인하는 데 사용될 수 있습니다.
·flsgrf.*: 솔버 출력 파일입니다. 시뮬레이션의 최종 결과가 포함됩니다.
·prperr.*, report.*, prpout.*: 이 파일들은 Preprocessor Diagnostic Files.
·hd3err.*, hd3msg.*, hd3out.*: 이 파일들은 Solver Diagnostic Files.
모든 시뮬레이션 파일은 단일 폴더에 함께 유지하므로, 설명이 될 수 있는 시뮬레이션 이름을 사용하는 것이 좋습니다. 그러나 매우 긴 파일 이름은 운영 체제에 따라 문제가 될 수 있습니다.
노트
시뮬레이션 이름이 inp(즉, 입력 파일이 있다면 prepin.inp) 출력 및 진단 파일은 모두 .dat이름을 갖습니다. 예: flsgrf.dat.
모든 입력 파일은 네트워크 위치의 컴퓨터 대신 로컬 디렉토리에 저장하는 것이 좋습니다. 이것은 솔버가 더 빠르게 실행되고 GUI의 응답 속도가 빨라지며 실행중인 시뮬레이션을 방해하는 네트워크 문제 가능성을 제거합니다.
3.시뮬레이션 관리자
FLOW-3D시뮬레이션 관리자의 탭은 주로 시뮬레이션을 실행할 수 있도록 시뮬레이션 환경을 구성하고 실행 시뮬레이션에 대한 상태 정보를 표시하는데 사용됩니다.
작업 공간 (Workspaces)
작업 공간(Workspaces)은 Simulation Manager의 필수 부분이며 파일을FLOW-3D에서 처리하는 방식입니다. 기본적으로 시뮬레이션을 포함하고 구성하는 폴더입니다. 몇 가지 예를 들면 시뮬레이션과 또 다른 작업 공간인 검증 사례를 포함하도록 할 수 있습니다:
새로운 작업 공간 만들기
이 튜토리얼에서는 작성하려는 시뮬레이션을 포함할 작업 공간(Workspaces)을 작성하십시오.
1.File -> New workspace로이동
2.작업 공간 이름으로 Tutorial를 입력하십시오.
3.기본 위치는 현재 사용자의 홈 디렉토리에 있습니다. 다른 곳에서 찾을 수 있지만 기본 위치가 우리의 목적에 적합합니다.
4.하위 디렉토리를 사용하여 작업 공간 이름 만들기확인란을 선택합니다. 이렇게 하면 파일 시스템에서 작업 공간에 대한 새로운 하위 디렉토리가 만들어져 시뮬레이션 파일을 훨씬 쉽게 구성할 수 있습니다.
5.확인을 눌러 새 작업 공간을 작성하십시오. 이제 포트폴리오에표시됩니다.
작업 공간 닫기
포트폴리오를 정리하고 탐색하기 쉽도록 필요 없는 작업공간을 닫는 것이 편리합니다. 작업 공간을 닫으면 포트폴리오에서 해당 작업 공간만 제거됩니다. 그러나, 컴퓨터에서 작업 공간을 삭제하지는 않습니다.
작업 공간을 닫으려면
1.기존 작업 공간을 마우스 오른쪽 버튼으로 클릭하고 작업 Close Workspace를선택하십시오. 또는 포트폴리오에서 작업 공간을 선택 (왼쪽 클릭) 하고 Delete 키를 누를 수 있습니다.
2.작업 공간을 닫을 것인지 묻는 메세지가 표시됩니다. 예를선택하십시오.
3.포트폴리오는더 이상 닫힌 작업 공간을 포함하지 않습니다.
기존 작업 공간 열기
오래된 작업 공간을 열어야 할 때가 있을 것입니다. 예를 들어, 새 프로젝트에 유사한 시뮬레이션을 작성하기 전에 기존 시뮬레이션의 설정을 검토할 수 있습니다. 기존 작업 공간을 열려면
1.File -> Open Workspace…를 선택하십시오
2.작업 공간 파일이 있는 디렉토리를 찾으십시오. Tutorial.FLOW-3D_Workspace.
3.작업 공간을 로드 하려면 OK을 누르십시오.
작업 공간에서 시뮬레이션 작업
작업 공간을 사용하는 방법을 알았으니, 여기에 시뮬레이션을 추가해 봅시다.
Example를 추가하십시오
작업 공간에 작업 시뮬레이션을 추가하는 가장 간단한 방법은 포함된 예제 시뮬레이션 중 하나를 추가하는 것입니다. FLOW-3D의 다양한 기능을 사용하는 방법을 보여주기 위해 설계된 간단하고 빠른 시뮬레이션입니다. 기존 작업 공간에 예제를 추가하려면 다음을 수행하십시오.
1.포트폴리오에서 원하는 작업 공간을 강조 표시하십시오
2.File -> Add example…을선택하십시오. 또는 작업공간을 마우스 오른쪽 버튼으로 클릭하고 예제 추가…를 선택할 수 있습니다.
3.예제대화 상자에서 예제를 선택하고 열기를누르십시오. 자연 대류(Natural Convection)예제를 선택했습니다.
4.새 시뮬레이션대화 상자가 열립니다.
5.디렉토리가 작업 공간 위치에 있는지 확인하는 것이 좋으므로 기본 시뮬레이션 이름과 위치를 잘 확인하는 것이 좋습니다. FLOW-3D는 모든 시뮬레이션 파일을 이 작업 공간 디렉토리의 별도 하위 디렉토리에 배치하여 파일 구성을 쉽게 만들어 줍니다.
6.시뮬레이션을 위한 단위 시스템을 선택하십시오. 표준 단위 시스템이 권장되지만 각 단위를 독립적으로 선택하기 위해 사용자 지정 단위 시스템을 선택할 수 있습니다.
7.확인을 눌러 새 시뮬레이션을 작업 공간에 추가하십시오.
작업 공간에서 시뮬레이션 제거
작업 공간에서 시뮬레이션을 제거해야 하는 경우가 있습니다 (이는 작업 공간에서 시뮬레이션을 제거만 하며, 컴퓨터에서 시뮬레이션을 삭제하지는 않습니다). 작업 공간에서 시뮬레이션을 제거하려면 다음을 수행하십시오.
1.작업 공간에서 기존 시뮬레이션을 마우스 오른쪽 버튼으로 클릭하고 (이 경우 이전 섹션에서 추가 한 예제 사용) 시뮬레이션 제거를선택하십시오. 또는 작업 공간에서 시뮬레이션을 선택 (왼쪽 클릭)하고 Delete 키를 누를 수 있습니다.
2.작업 공간에는 더 이상 시뮬레이션이 포함되지 않습니다.
모든 작업 공간 및 디스크에서 시뮬레이션 삭제
작업 공간에서 시뮬레이션을 제거하는 것 외에도 디스크에서 모든 시뮬레이션 파일을 삭제해야 할 수도 있습니다. 작업 공간에서 시뮬레이션을 제거하고 디스크에서 시뮬레이션 파일을 삭제하려면 다음을 수행하십시오.
1.작업 공간에서 기존 시뮬레이션을 마우스 오른쪽 단추로 클릭하고 (이 경우 이전 섹션에서 추가 한 예제 사용) 모든 작업 공간 및 디스크에서 시뮬레이션 삭제를선택하십시오.
2.시뮬레이션 디렉토리에서 삭제할 파일을 선택할 수 있는 창이 나타납니다. 삭제할 파일을 선택한 다음 확인을 눌러 해당 파일을 삭제하거나 취소를 눌러 작업을 중단하십시오.
3.OK를 선택한 경우 선택한 작업 공간은 더 이상 시뮬레이션을 포함하지 않습니다. 선택한 작업 공간의 모든 시뮬레이션 파일은 디렉토리에서 삭제됩니다.
경고
이 작업은 취소할 수 없으므로 계속하기 확인 후 파일을 삭제해야 합니다.
작업 공간에 기존 시뮬레이션 추가
기존 시뮬레이션을 작업 공간에 추가하려면 다음을 수행하십시오.
1.열린 작업 공간을 마우스 오른쪽 버튼으로 클릭하고 기존 시뮬레이션 추가…를선택합니다. 작업 공간을 선택한 다음 File->Add Existing Simulation…을 선택할 수도있습니다.
2.prepin.*파일 위치로 이동하여 열기를선택하십시오.
3.시뮬레이션이 이제 작업 공간에 나타납니다.
작업 공간에 새로운 시뮬레이션 추가
대부분의 경우 기존 시뮬레이션을 사용하는 대신 새 시뮬레이션을 작성하게 됩니다. 작업 공간에 새로운 시뮬레이션을 추가하려면:
1.기존 작업 공간을 마우스 오른쪽 버튼으로 클릭하고 새 시뮬레이션 추가…를선택하십시오.
2.시뮬레이션 이름을 입력하라는 message가 표시됩니다. 이 예제에서는 heat transfer example불러오십시오.
3.그런 다음 드롭다운 목록을 사용하여 시뮬레이션을 위한 단위 시스템을 결정합니다. 사용 가능한 옵션은 질량, 길이, 시간, 전기요금 각각 g, cm, s, coul기준의 Kg, m, s, CGS입니다. 또한 엔지니어링 단위도 사용할 수 있으며, slug, ft, s의 기초 단위가 있지만, 전기 충전을 위한 단위는 없습니다. 이러한 옵션 중 어느 것도 해당되지 않는 경우, 질량, 길이, 시간 및 전기요금에 대한 기준 등을 사용자 정의하여 사용자 지정 단위 시스템을 사용할 수 있습니다.
4.온도 단위는 드롭다운 목록을 사용하여 지정해야 합니다. 사용 가능한 옵션은 SI 및 CGS 단위의 경우 Celsius 및 Kelvin, 엔지니어링 단위의 경우 Fahrenheit 및 Rankine입니다. Custom units(사용자 정의 단위) 옵션을 선택한 경우, 사용 가능한 온도 단위는 질량 및 길이에 대해 선택한 기본 단위에 따라 변경됩니다.
노트
새 시뮬레이션의 시뮬레이션 단위는 신중하게 선택하십시오. 일단 설정하면 단위를 변경할 수 없습니다.
5.이 시뮬레이션에 사용된 템플릿이 기본 템플릿이 됩니다. 템플릿은 포함된 설정을 새 시뮬레이션에 적용하는 저장된 값 세트입니다. 다른 템플릿을 사용해야하는 경우 찾아보기 아이콘 ( )을 클릭하여 사용 가능한 템플릿 목록에서 선택하십시오.
6.기본 시뮬레이션 이름과 위치는 디렉토리가 작업 공간 위치에 있는지 확인하는 것이 좋습니다. FLOW-3D는 모든 시뮬레이션 파일을 이 작업 공간 디렉토리의 별도 하위 디렉토리에 배치하여 파일 구성을 훨씬 쉽게 만듭니다. 시뮬레이션을 다른 위치에 저장하려면 찾아보기 아이콘 ( )을 사용하여 원하는 위치로 이동하십시오.
7.확인을 클릭하여 작업 공간에 새 시뮬레이션을 추가하십시오.
heat transfer example
다른 옵션
우리는 지금 이러한 옵션을 사용하지 않는 동안, 이 시뮬레이션을 마우스 오른쪽 버튼으로 클릭하여 추가 옵션에 대한 액세스를 제공합니다.
일반적으로 사용되는 Add Simulation Copy…그리고 Add Restart Simulation…을 추가합니다. 첫 번째 옵션은 기존 시뮬레이션의 사본을 작성하고, 두 번째 옵션은 기존 시뮬레이션을 복사하고 원래 시뮬레이션의 결과를 다시 시작 시뮬레이션의 초기 조건으로 사용하도록 다시 시작 옵션을 구성합니다.
추가 정보
재시작 시뮬레이션에 대한 자세한 내용은 도움말에서 모델 설정 장의 재시작 섹션을 참조하십시오.
전처리 및 시뮬레이션 실행
시뮬레이션 전처리
시뮬레이션 전처리는 초기 조건을 계산하고 입력 파일에서 일부 진단 테스트를 실행합니다. 문제가 올바르게 구성되었는지 확인하거나 전 처리기의 진단 정보가 필요한 경우에 유용합니다. 시뮬레이션을 실행하기 전에 전처리할 필요가 없습니다. 시뮬레이션을 전처리 하려면
1.작업 공간에서 시뮬레이션을 마우스 오른쪽 버튼으로 클릭하고 Preprocess Simulation->Local을선택합니다. 이 경우 입력 파일 heat transfer example이 아직 완전히 정의되지 않았으므로 작업 공간에서 예제 문제를 선택하십시오.
2.전처리 프로세스가 시작되고 Simulation Manager 하단의 텍스트 창에 일부 정보가 인쇄된 후 성공적으로 완료됩니다. 포트폴리오에서 시뮬레이션 이름 옆의 아이콘도 시뮬레이션이 성공적으로 처리되었음을 나타내도록 변경됩니다.
추가 정보
자세한 내용은 도움말의 컴퓨팅 결과 장의 전처리 섹션을 참조하십시오.
시뮬레이션 실행
시뮬레이션을 실행하면 입력 파일에 정의된 문제에 대한 지배 방정식(물리적 모델, 형상, 초기 조건, 경계 조건 등)이 해석됩니다. 시뮬레이션을 실행하려면
1.작업 공간에서 시뮬레이션을 마우스 오른쪽 버튼으로 클릭하고 Run Simulation->Local을 선택하십시오. 이 경우 입력 파일 heat transfer example이 아직 완전히 정의되지 않았으므로 작업 공간에서 예제 문제를 선택하십시오.
2.솔버가 시작되고 시뮬레이션 관리자 하단의 텍스트 창에 일부 정보가 인쇄되고 플롯이 업데이트 된 후 성공적으로 완료됩니다. 포트폴리오에서 시뮬레이션 이름 옆의 아이콘도 시뮬레이션이 성공적으로 실행되었음을 나타내도록 변경됩니다. 또한 솔버가 실행되는 동안 큐에 시뮬레이션이 나타나는 것을 볼 수 있으며, 완료되면 사라집니다.
추가 정보
시뮬레이션 실행 및 진단 읽기에 대한 자세한 내용은 도움말의 컴퓨팅 결과 장에서 솔버 실행 섹션을 참조하십시오.
작업 공간에서 모든 시뮬레이션 실행
작업 공간을 마우스 오른쪽 버튼으로 클릭하고 Simulate Workspace->Local을 선택하여 작업 공간에서 모든 시뮬레이션을 실행할 수도 있습니다.
추가 정보
자세한 내용은 컴퓨팅 결과 장에서 솔버 실행 섹션을 참조하십시오.
대기열
사전 처리 또는 실행에 작업이 제출되면 큐의 맨 아래에 시뮬레이션이 자동으로 추가됩니다. 그런 다음 솔버에 사용 가능한 라이센스 및 계산 리소스가 있으면 시뮬레이션이 사전 처리되거나 실행됩니다. 대기열에 있지만 아직 전처리 또는 실행되지 않은 시뮬레이션은 대기열 맨 아래의 컨트롤을 사용하여 대기열에서 다시 정렬하거나 대기열에서 제거할 수 있습니다.
추가 정보
자세한 내용은 컴퓨팅 결과 장을 참조하십시오.
파일 시스템에서 파일 찾기
어떤 이유로 구조물 파일에 액세스해야 하는 경우 (아마 *.STL 폴더에 파일을 배치해야 함) 표시된 파일 경로를 시뮬레이션 입력 파일로 클릭하여 파일 시스템의 해당 위치로 이동할 수 있습니다.
4.모델 설정
Model Setup(모델 설정) 탭은 시뮬레이션 관리자에서 현재 선택한 시뮬레이션에 대한 입력 매개 변수를 정의하는 곳입니다. 여기에는 전역설정, 물리학 모델, 유체, 기하학, 메싱, 구성요소 특성, 초기 조건, 경계 조건, 출력 옵션 및 숫자가 포함된다.
이 섹션은 물에 잠긴 모래(위; 파랑)의 바닥에서 가열된 구리 블록(위; 빨간색)에 의해 발생하는 열 기둥(아래)을 보여주는 간단한 시뮬레이션 설정 방법을 안내합니다.
이 튜토리얼은 방법이나 모델이 어떻게 작동하는지, 옵션을 선택한 이유 등에 대한 포괄적인 논의를 의도한 것이 아니며, 이 특정 시뮬레이션을 설정하기 위해 수행해야 할 사항에 대한 간략한 개요일 뿐입니다. 여기서 행해지는 것에 대한 방법/모델과 추론의 세부사항은 사용 설명서의 다른 장에서 확인할 수 있습니다.
시작하려면 새 작업 공간을 작성하고 새 시뮬레이션을 추가하십시오. 이를 수행하는 방법에 대한 지침은 새 작업 공간 작성 및 작업 공간에 새 시뮬레이션 추가를 참조하십시오.
탐색
모델 설정은 주로 빨간색으로 표시된 처음 9 개의 아이콘의 탐색을 통해 수행됩니다. 각 아이콘은 시뮬레이션의 특정 측면을 구성하기 위한 위젯을 엽니다. Global에서 시작하여Numerics로 끝나는 다음 섹션은 각 위젯의 목적을 보여줍니다.
통제 수단
다음은 FLOW-3D 사용자 인터페이스의 그래픽 디스플레이 영역에서 사용되는 마우스 컨트롤입니다.
행동
버튼/키
동작
기술
회전
왼쪽
길게 클릭
마우스 왼쪽 버튼을 클릭 한 채로Meshing & Geometry창에서 마우스를 움직입니다. 그에 따라 모델이 회전합니다.
줌
중간 버튼/스크롤 휠
스크롤/클릭 한 상태
마우스를 앞뒤로 움직여 확대/축소하려면 가운데 휠을 굴리거나 마우스 가운데 버튼을 클릭 한 상태로 유지하십시오.
팬
우측
길게 클릭
마우스 오른쪽 버튼을 클릭 한 채로 창에서 마우스를 움직입니다. 모델이 마우스와 함께 움직입니다.
객체에 초점 설정
해당 없음
객체 위에 커서를 놓기
커서를 개체 위로 가져 가면 마우스 오른쪽 버튼 클릭 메뉴를 통해 추가 조작을 위해 개체가 활성화됩니다. 개체가 활성화되면 강조 표시됩니다. Meshing & Geometry탭에서Tools->Mouse Hover Selection환경 설정 이 활성화된 경우에만 수행됩니다.
선택
왼쪽
더블 클릭
객체를 두 번 클릭하면 마우스 오른쪽 버튼 메뉴를 통해 추가 조작을 위해 객체를 선택하고 활성화합니다. Meshing & Geometry탭에서Tools->Mouse Hover Selection환경 설정 이 비활성화 된 경우에만 활성화됩니다.
액세스 객체 속성
우측
딸깍 하는 소리
강조 표시된 객체를 마우스 오른쪽 버튼으로 클릭하면 객체 식별, 표시/숨기기, 활성화/비활성화, 투명도 조정 등의 옵션 목록이 표시됩니다.
커서 좌표 반환 (프로브)
왼쪽
Shift + 클릭
Shift 키를 누르면 커서가 대상으로 바뀝니다. Shift 키를 누른 상태에서 클릭하면 화면의 왼쪽 하단에 표시된 표면의 좌표가 표시됩니다.
피벗 점 배치
왼쪽
cntrl + 클릭
Ctrl 키를 누르고 있으면 커서가 피벗 아이콘으로 바뀝니다. Ctrl 키를 누른 상태에서 클릭하여 피벗 점을 설정하십시오. 뷰가 피벗 점을 중심으로 회전합니다. 토글사용자 정의 피벗피벗 점을 끕니다. 보기 창 위의 버튼을 누릅니다.
도움이 되는 툴바 옵션도 있습니다. 옵션의 목적을 찾으려면 아이콘 위로 마우스를 가져갑니다.
글로벌
이 매뉴얼에 대한 시뮬레이션을 만들려면 원하는 작업 공간을 마우스 오른쪽 단추로 클릭하고 새 시뮬레이션 추가…를 선택하십시오. 매뉴얼 섹션의 새 시뮬레이션 추가 작업 공간에 설명된 대로 이름을 ‘heat transfer example’로 지정하고 작업 공간에 추가하십시오. SI와 Kelvin을 각각 단위 시스템과 온도로 선택합니다. 일단 설정되면 시뮬레이션을 위한 단위는 변경할 수 없다는 점을 기억하십시오.
글로벌 아이콘을 클릭하여 글로벌 위젯을 여십시오. 여기에서 정의된 단위가 표시되고 시뮬레이션 완료 시간이 설정됩니다. 이 시뮬레이션의 경우 완료 시간을 200 초로 설정하십시오. 시뮬레이션에 대한 중요한 세부 정보는 여기 노트 필드에도 추가할 수 있습니다.
추가 정보
자세한 내용은 모델 설정 장의 전역 섹션을 참조하십시오.
물리
물리아이콘을 클릭하여 물리 위젯을 엽니다.
이 문제의 경우, 하나의 유체, 자유 표면, 경계 및 비압축/제한 압축의 기본 설정이 모두 정확합니다.
관련 물리 메커니즘(즉, 추가 지배 방정식 또는 지배 방정식 용어)은 물리 위젯에서 정의됩니다. 모델을 활성화하려면 해당 모델의 아이콘을 마우스 왼쪽 버튼으로 클릭하고 ‘활성화‘를 선택하십시오. 이 시뮬레이션을 위해서는 다음 모델을 활성화해야 합니다.
·Density evaluation(밀도 평가): 이 모델은 열 기둥을 생성하는 밀도 변화를 설명합니다. 다른 양(예: 온도 또는 스칼라)의 함수로 평가된 밀도를 선택하고 Include volumetric thermal expansion 상자를 선택하십시오.
·Gravity and non-inertial reference frame(중력 및 비 관성 기준 프레임): 중력을 나타내는 힘이 추가되므로 Z 중력 성분에 -9.81을 입력하십시오.
· Heat transfer(열 전달): 이 모델은 유체와 고체 물체 사이의 열 전달을 설명합니다. 이 시뮬레이션의 경우 First order for the Fluid internal Energy advection를 선택하고 Fluid to solid heat transfer를 활성화하려면 확인란을 선택하십시오. 나머지 옵션은 기본값으로 두어야합니다.
· Viscosity and turbulence(점성 및 난류): 이 모델은 유체의 점성 응력을 설명합니다. Viscous flow 옵션을 선택하고 나머지 옵션은 기본값으로 두십시오.
추가 정보
자세한 내용은 모델 설정 장의 물리 섹션을 참조하십시오.
유체
유체의 속성은 모델 설정 탭의 유체 위젯에 정의되어 있습니다. 유체 위젯은 수직 도구 모음에서 Fluids f3d_fluids_icon 아이콘을 클릭하여 액세스할 수 있습니다. 먼저 유체 옵션 1 이 속성 옵션으로 선택되어 있는지 확인하십시오. 유체 1의 속성은 수동으로 입력할 수 있지만 일반적인 유체의 속성을 설정하는 빠른 방법은 재료 속성로드 버튼을 클릭하여 재료 데이터베이스에서 유체를 로드하는 것입니다. 다음으로, 원하는 재료를 탐색하십시오. 이 경우 Fluids->Liquids->Water_at_20_C를 선택하고 Load를 클릭하십시오.
이 시뮬레이션에는 데이터베이스에 없는 특성인 체적 열 팽창 계수가 필요합니다. 밀도 하위 탭에서 을 입력하십시오. 최종 속성 세트는 다음과 같아야 합니다.
추가 정보
자세한 내용은 모델 설정 장의 유체 섹션을 참조하십시오.
Geometry(기하)
기하형상아이콘을 클릭하여 물리 위젯을 엽니다.
이 시뮬레이션을 위해 생성해야 하는 두 가지 형상은 구리 블록과 모래층이 있습니다. 둘 다 프리미티브를 사용하여 작성합니다. 보다 현실적인 시뮬레이션은 Primitives, Stereolithography(STL) Geometry File (s)/또는 Raster File (s)을 사용하여 지오메트리를 정의할 수 있습니다.
구리 블록을 만들려면 먼저 지정된 형상 아이콘을 클릭하여 작성합니다. 구리 블록을 x 및 y 방향 원점에서 +/- 2cm 연장하고 z 방향으로 0-4cm 연장합니다. 나머지 옵션은 그대로 두고 블럭을 솔리드로 만들고 새 구성 요소에 추가합니다.
하위 구성 요소 정의를 마치고 구성 요소 정의로 이동하려면 확인을 선택하십시오. 자동으로 열린 구성요소 추가 대화상자에서 Type as General(솔리드)을 그대로 두고 Name(이름) 필드에 Copper block을 입력한 다음 OK(확인)를 선택하여 구성요소 정의를 완료하십시오.
아이콘을 다시 클릭하여 베드 하위 구성 요소를 작성하십시오. 아래 표시된 범위를 사용하고 컴포넌트에 추가 선택 사항을 새 컴포넌트(2)로 설정하십시오.
하위 구성 요소 정의를 마치고 구성 요소 정의로 이동하려면 확인을 선택하십시오. 대화 형으로 이름 필드에서Bed를 입력한 후 구성요소 정의를 마칩니다. 최종 형상은 다음과 같이 표시됩니다.
새 구성 요소를 추가하면 가로 및 세로 방향으로 그래픽 표시 창에 길이 스케일이 자동으로 생성됩니다. 이 눈금자 도구를 사용하여 생성된 기하학적 객체의 범위를 빠르게 측정할 수 있습니다.
노트
표시 영역에는 지오메트리 모양 정의만 표시되므로 객체가 솔리드인지 구멍인지에 대한 정보는 표시되지 않습니다. 옵션을 사용하여 Mesh 후에 나중에 수행할 수 있습니다.
추가 정보
자세한 내용은 도움말 모델 설정 장의 형상 섹션을 참조하십시오.
구성 요소 속성
열전달 모델은 고체 구성 요소의 전도 방정식을 해결하기 위해 재료 특성이 필요합니다. 이러한 속성은 이 아이콘을 클릭하여 구성 요소 속성 위젯에서 설정합니다.
각 구성 요소에는 솔리드 특성 및 표면 특성이 정의 되어 있어야합니다. 구리 블록에 대해 이를 설정하려면 먼저 형상 위젯에서 구성 요소 1: copper block 요소를 선택하십시오. 그런 다음 컴포넌트 특성 위젯에서 솔리드 특성을 선택하고 다음과 같이 특성을 정의하십시오.
여기에서 두 번째 구성 요소(베드)에 대해 설명된 구성 요소 특성 정의를 위한 대체 방법을 사용할 수 있습니다. 이 방법에서는 구성 요소 2: 베드 구성 요소를 클릭하고 재료 필드 옆에 있는 재료 특성로드아이콘을 선택하여 시작합니다. 다음으로 재료를 탐색합니다. 이 경우 Solids->Sands->Sand_Quartz를 선택하고 Load를 선택하십시오.
추가 정보
l 자세한 내용은 모델 설정 장의 유체 섹션을 참조하십시오.
l 주어진 물리적 모델에 필요한 속성에 대한 자세한 내용은 모델 참조 장을 참조하십시오.
Meshing(메싱)
Mesh는 Mesh 위젯에서 생성 및 정의되며, 위젯을 통해 액세스 할 수 있습니다. 이 아이콘을 눌러새 Mesh를 추가합니다. Mesh의 범위를 형상에 빠르게 적용하려면 형상에 맞추기 라디오 버튼을 선택하고 오프셋 라디오 버튼을 백분율로 유지합니다. 블록 속성에서 셀 크기를 0.004로 설정하십시오.
Mesh 상단은 z 방향으로 위쪽으로 확장해야 합니다. Z-Direciton탭을 선택하고 Mesh Plane 2에 0.2를 입력합니다.
이 시뮬레이션은 2D가 될 것입니다. 동일한 프로세스에 따라 Y 방향 범위를 -0.005 및 0.005 로 설정하십시오. 그리고 합계 셀을 1로 설정하십시오.
최종 Mesh는 그래픽 디스플레이 창 바로 위의 Mesh->Flow Mesh->View 모드 드롭 다운 메뉴에서 옵션을 변경하여 다른 방식으로 볼 수 있습니다. 그리드 라인 마다 그리드 선을 표시합니다 옵션은 Mesh Plane의 옵션만 표시됩니다 Plane Mesh 및 개요 옵션은 Mesh의 범위를 보여줍니다.
또한 솔버가Mesh의 최종 지오메트리를 인식하는 방법은 FAVOR TM 알고리즘을 사용하여 형상 정의를 면적 분수 및 부피 분수로 변환합니다. 이렇게 하려면아이콘을 클릭한 다음 생성을 선택하십시오.
잠시 후 회색 영역이 고체 물질을 나타내는 아래와 같은 형상을 표시해야 합니다.
추가 정보
l Mesh에 대한 자세한 내용은 모델 설정 장의 Mesh 섹션을 참조하십시오.
l FAVORTM및FAVORize 옵션에 대한 자세한 내용은 모델 설정 장의 Reviewing the FAVORized Geometry and Mesh 섹션을 참조하십시오.
경계 조건
FLOW-3D는 구성 요소 유형 및 활성 물리적 모델에 기초한 구성 요소에 적절한 경계 조건을 자동으로 적용합니다. 그러나 경계 조건 위젯에서 Mesh 블록면의 경계 조건은 각 Mesh 블록에 대해 수동으로 설정해야 합니다( ).
이 매뉴얼의 경우 경계 조건 중 3 가지가 경계조건( X Min , X Max, Z Max 경계)을 기본 대칭 조건조건부터 변경해야 합니다.
·X Min :
o경계 조건 위젯의 경계 섹션 아래에 있는 X Min 목록을 클릭하십시오. Type에서 경계 유형을 Velocity로 설정하고 X 속도에 대해 0.001을 입력하십시오.
o경계 조건 위젯의 경계 섹션 아래에 있는 Z 최대 목록을 클릭하십시오. 경계 유형을 압력으로 설정하고 압력에 대해 0을 입력하십시오.
o다음으로 유체 분율을 0.0으로 설정하십시오.
o마지막으로 온도를 298K로 맞춘다.
추가 정보
자세한 내용은 모델 설정 장의 Mesh 경계 조건 섹션을 참조하십시오.
초기 조건
도메인 내부의 솔리드 객체(구성 요소)와 유체 모두에 대해 초기 조건을 설정해야 합니다.
· 구성 요소 :이 시뮬레이션에서 솔리드 객체에 필요한 유일한 초기 조건은 초기 온도입니다. 이것은 각 구성 요소에 대한 위젯에 설정되어 있는 구성 요소 속성에 대해 수행한 것과 유사한 방식으로 구성 요소를 등록합니다. 구성 요소 속성을 설정할 때 이전과 동일한 방법으로 구성 요소 1의 초기 온도를 350K로 설정하고 구성 요소 2의 초기 온도를 298K로 설정하십시오.
유체: 유체의 초기 조건을 설정하기 위해 조금 더 설정해야 합니다. 이 경우 유체 구성, 온도, 속도 및 압력 분포를 모두 설정해야 합니다. 유체 초기 조건은 초기 위젯을 설정하고 초기 f3d_initial_icon를 클릭하면 열립니다.
아이콘을 선택한 후 유체 목록에서 압력을 선택하고 온도를 298K로 설정합니다. x, y, z 속도를 0.0으로 설정하십시오.
다음으로, 높이/볼륨 목록과 유체 높이 사용 드롭다운 버튼을 선택합니다. 유체 높이를 0.15로 설정하십시오.
추가 정보
자세한 내용은 모델 설정 장의 초기 조건 섹션을 참조하십시오.
출력
FLOW-3D 옵션에는 결과 파일에 기록될 데이터와 출력 위젯에서 발견된 빈도를 제어하는 7가지 데이터 유형이 있습니다. 출력아이콘을 클릭합니다.
다른 데이터 유형은 다음과 같습니다.
·Restart: 모든 흐름 변수. 기본 출력 주기는 시뮬레이션 시간의 1/10입니다.
·Selected: 사용자가 선택한 흐름 변수 만. 기본 출력 주기는 시뮬레이션 시간의 1/100입니다.
·History: 하나의 변수와 시간의 변화를 보여주는 데이터. 예는 시간 단계 크기, 평균 운동 에너지, 배플에서의 유속 등을 포함합니다. 기본 출력 주기 = 시뮬레이션 시간의 1/100.
·Short print: hd3msg.*파일에 텍스트 진단 데이터가 기록 됩니다. 기본 출력 주기는 시뮬레이션 시간의 1/100입니다.
·Long print : hd3out.*파일에 텍스트 진단 데이터가 기록 됩니다. 기본 출력 주기는 시뮬레이션 시간의 1/10입니다.
·Solidification: 응고 모델이 활성화 된 경우에만 사용 가능합니다.
·FSI TSE: 변형 가능한 솔리드에 대한 추가 출력 옵션.
일반적으로 이 시뮬레이션에는 기본 출력 속도가 적합합니다. 그러나 Selected Data의 일부 추가 구성은 유용합니다. Selected data interval을 0.5로 설정한 다음 Fluid 온도, Fluid velocity, Macroscopic density 및 Wall 온도 옆에 있는 상자를 선택합니다. 그러면 이러한 값이 0.5초마다 출력됩니다.
추가 정보
자세한 내용은 모델 설정 장의 출력 섹션을 참조하십시오.
Numerics
기본 Numerics 옵션은 대부분의 시뮬레이션에서 잘 작동하므로 기본 옵션에서 벗어나야 하는 충분한 이유가 없는 경우에는 현재 그대로 두는 것이 가장 좋습니다.
이것으로 모델 설정 섹션에서 시작된 예제 문제의 설정을 마칩니다. 이제 실행할 준비가 되었으므로 전처리 및 시뮬레이션 실행의 단계에 따라 시뮬레이션을 실행하십시오.
추가 정보
자세한 내용은 모델 설정 장의Numerics옵션 섹션을 참조하십시오.
일반 시뮬레이션 설정 점검 목록
시뮬레이션을 설정하는 데 필요한 단계에 대한 개략적인 개요가 아래에 나와 있습니다. 이 목록은 포괄적인 목록이 아닙니다. 일반적인 단계, 고려해야 할 몇 가지 중요한 사항 및 권장되는 설정 순서를 간단히 설명하는 안내서일 뿐입니다.
시작하기 전에
1.물리적 문제의 다이어그램을 그리기 및 주석 달기 : 이 다이어그램에는 기하학적 치수, 유체의 위치, 관련 힘, 움직이는 물체의 속도, 관련 열 전달 메커니즘 등이 포함되어야 합니다. 완성된 다이어그램은 문제에 대한 모든 관련 엔지니어링 정보로 인한 물리적 문제에 대한 이미지여야 합니다.
2.모델링 접근법 결정: 주석이 달린 다이어그램을 가이드로 사용하여 문제점에 접근하는 방법을 결정 : 문제가 되는 유체의 수, 혼화 가능한 경우, 하나 이상의 유체에서 방정식을 풀어야하는 경우 및 압축성이 중요한지 파악하여 시작하십시오. 그런 다음 어떤 물리적 메커니즘이 중요한지 결정하십시오. 이러한 각 옵션 (예: 유체 유형, 열 전달 메커니즘 등)에 대한 관련 엔지니어링 정보를 다이어그램에 추가하십시오. 물리적 메커니즘이 포함되거나 무시된 이유를 정당화하려고 합니다. 이를 통해 시뮬레이션 프로세스 초기에 오류를 수정하는 데 시간이 거의 걸리지 않는 초기에 실수를 잡을 수 있습니다.
3.다이어그램에 계산 영역을 그리고, 계산 영역의 가장자리에 있는 물리적 상황 설명 : 경계의 물리적 상황을 가장 잘 나타내는 경계 조건 유형을 기록합니다. 사용 가능한 경계 조건 유형이 경계의 물리적 상황에 대한 합리적인 근사치가 아닌 경우 이 경계를 다른 곳으로 이동해야 합니다.
모델 설정 : 일반
1.문제, 시뮬레이션의 목적, 사례 번호 등을 설명하는 메모를 추가하십시오. 메모는 향후 사용자 또는 나중에 참조할 수 있도록 설정을 설명하고 정당화하는 데 도움이 됩니다. 시뮬레이션의 목적, 분석 방법 등을 논의해야합니다.
2.사용할 솔버와 프로세서 수를 선택하십시오.
3.단위 시스템 선택: 소규모 문제를 모델링 할 때는 작은 단위 (예 : mm-gm-msec)를 사용하고 규모가 큰 문제는 큰 단위 (예 : SI)를 사용하십시오. 이를 통해 기계 정밀도로 인한 반올림 오류를 방지할 수 있습니다.
4.유체 수, 인터페이스 추적 옵션 및 유량 모드를 선택하십시오. 주석이 달린 다이어그램을 이 단계의 지침으로 사용하십시오. 유체의 수는 질량, 운동량 및 에너지 보존을 관장하는 방정식이 유체 분율 1을 나타내는) 또는 유체 분획 (유체 1 및 유체 2)이 있는 영역에서 해결되는지 여부를 나타냅니다. 인터페이스 추적 옵션은 유체 분율의 변화가 급격한지 또는 확산되어야 하는지 여부를 정의하는 반면, 흐름 모드는두 유체 문제에서 처리되는 영역을 정의합니다.
5.마감 조건 정의: 시뮬레이션 종료 시점을 선택합니다. 시간, 채우기 비율 또는 기타 정상 상태 측정을 기반으로 할 수 있습니다.
6.기존 결과에서 시뮬레이션을 다시 시작하는 방법 정의 (선택 사항): 기존 결과 파일에서 시뮬레이션을 다시 시작할 때 다시 시작 옵션이 적용됩니다. 재시작 옵션은 재시작 소스 파일에서 가져온 정보와 시뮬레이션의 초기 조건을 사용하여 재설정되는 정보를 정의합니다.
모델 설정 : 물리
1.주석이 달린 다이어그램을 기반으로 관련 실제 모델 활성화
모델 설정 : 유체
1.유체의 속성 정의 1: 주석이 달린 다이어그램을 가이드로 사용하여 활성 물리적 모델에 대한 적절한 물리적 속성을 정의하십시오.
2.유체 2의 속성 정의 (사용하는 경우): 주석이 달린 다이어그램을 가이드로 사용하여 활성 물리적 모델에 적절한 물리적 속성을 정의하십시오.
3.인터페이스의 속성 정의: f = 1과 f = 0의 영역 사이의 인터페이스 속성을 정의하십시오. 여기에는 표면 장력, 상 변화 및 확산에 대한 특성이 포함됩니다.
모델 설정 : Mesh 및 형상
1.모든 STL 파일의 오류 점검: ADmesh, netfabb Studio 또는 유사한 프로그램을 사용하여 모든 STL 파일의 오류를 점검하십시오. 이는 모델 설정에 시간을 소비하기 전에 형상 정의와 관련된 문제를 파악하는 데 도움이 됩니다.
2.모든 하위 구성 요소 및 구성 요소 가져 오기 및 정의 : 주석이 달린 다이어그램에 설명 된 대로 실제 사례와 일치하도록 3D 솔리드 형상을 정의합니다. 최종 결과는 물리적 형상의 정확한 복제본이어야 합니다. 각 부분에 설명적인 이름을 사용하고 대량 소스가 될 구성 요소를 포함하십시오.
3.모든 구성 요소의 속성 정의: 주석이 달린 다이어그램에 그려진 내용을 기반으로 각 구성 요소의 모든 재료 속성, 표면 속성, 모션 속성 등을 정의합니다. 경계 조건이 정의될 때까지 질량 소스 특성을 정의하기를 기다리십시오.
4.스프링과 로프 및 각각에 대한 관련 속성을 정의합니다.
5.주석이 달린 다이어그램에 설명된 시뮬레이션 도메인과 일치하도록 Mesh를 정의하십시오. 도메인의 모서리가 다이어그램에서 식별된 위치에 있는지 확인하십시오. 또한 인터페이스 (셀이있는 셀과 셀이다른 셀 이 있는 셀)를 식별하려면 세 개의 셀이 필요합니다. ). 최소 5 개의 셀이 예상되는 가장 얇은 연속 영역에 맞도록 충분히 작은 셀을 사용하십시오. 과 .
6.지오메트리를 정의하는 모든 배플 정의
7.경계 조건, 질량 소스, 질량 모멘텀 소스, 밸브 및 벤트 정의: 경계 조건 (질량 소스, 질량 모멘텀 소스, 밸브 및 벤트 포함)은 모든 방정식을 풀기 위해 주어진 위치에서 솔루션을 규정합니다. 주석이 달린 다이어그램을 사용하여 각 경계 (또는 소스 등)에 지정된 내용이 유동 솔루션, 열 전달 솔루션, 전위 등에 대한 현실과 일치하는지 확인하십시오.
8.유체 및 구성 요소의 초기 조건을 정의합니다. 초기 조건은 모든 방정식 (유량 솔루션, 열 전달 솔루션, 전위 등)에 대해 모든 영역에서 솔루션을 규정합니다. .주석이 달린 다이어그램을 사용하여 초기 조건에 지정된 내용이 현재 현실에 대한 근사치인지 확인하십시오. 유체 영역뿐만 아니라 구성 요소의 초기 조건을 설정해야 합니다.
9.모든 측정 장치 정의 (샘플링 볼륨, 플럭스 표면 및 히스토리 프로브)
모델 설정 : 출력
1.출력 기준 (시간, 채우기 비율 또는 응고된 비율)을 선택하십시오.
2.재시작 데이터에 추가할 출력을 선택하십시오.
3.선택한 데이터에 기록할 정보를 선택하십시오.
4.재시작, 선택, 히스토리, 짧은 인쇄 및 긴 인쇄 데이터의 출력 속도 정의 : 기본 속도는 재시작 및 긴 인쇄 데이터의 경우 (10개 출력)/(시뮬레이션 종료 시간) 및 선택한 기록, 짧은 인쇄 데이터의 경우 (100개 출력)/(시뮬레이션 종료 시간)입니다.
모델 설정 : 숫자
1.기본값이 아닌 필수 숫자 옵션을 선택 FLOW-3D의 숫자 옵션은 고급 사용자를 대상으로 하며, 지배 방정식을 해결하는 데 사용되는 숫자 근사치 및 방법을 상당히 제어할 수 있습니다. 이러한 옵션 중 일부를 잘못 사용하면 솔루션에 문제가 발생할 수 있으므로 일반적으로 이 옵션의 기능을 먼저 이해하고 조정의 정당성을 갖추지 않고는 이러한 설정을 조정하지 않습니다.
5.FLOW-3D에서 후 처리
이 섹션에서는 FLOW-3D에 통합된 포스트 프로세서를 사용하는 방법에 대해 설명합니다. 보다 강력한 외부 포스트프로세서FlowSight에 대한 튜토리얼은FlowSight설명서를 참조하십시오. 또한 이 섹션에서는 Flow Over A Weir 예제 문제를 실행하여 생성된 결과 파일을 사용합니다. 이 예제 문제를 실행하는 방법에 대한 지침은 예제 추가 및 시뮬레이션 사전 처리 및 실행을 참조합니다.
FlowSight 사용에 대한 기본 참조는FlowSight의 Help->helpLocal Help 메뉴에서 액세스하는FlowSight사용자 설명서입니다.
추가 정보
기존 플롯
기존 플롯은 솔버가 자동으로 생성하는 사전 정의된 플롯입니다. 사용자 정의 플롯은 아래의 사용자 정의 플롯 섹션에 설명되어 있습니다.
1.분석 탭을 클릭하십시오. FLOW-3D 결과 대화 상자가 표시됩니다; 메세지가 나타나지 않으면 (분석 탭이 열림) 결과 파일 열기를 선택하여 동일한 대화 상자를 엽니다.
2.기존 라디오 버튼을 선택하십시오. 데이터 파일 경로 상자에 두 가지 유형의 파일이 표시됩니다 (있는 경우). 이름이 prpplt.*있는 파일 에는 전처리flsplt.*기에 의해 자동으로 작성된 플롯이 포함되고 이름이 있는 파일에는 입력 파일에 사전 지정된 플롯 뿐만 아니라 후 처리기에 의해 자동으로 작성된 플롯이 포함됩니다.
3. 확인을 선택flsplt.Flow_Over_A_Weir하고 클릭하십시오. 그러면 디스플레이 탭이 자동으로 열립니다.
4.사용 가능한 플롯 목록이 오른쪽에 나타납니다. 목록에서 해당 플롯의 이름을 클릭하면 특정 플롯을 볼 수 있습니다. 플롯 26 이 아래에 나와 있습니다.
커스텀 플롯
1.분석탭으로 돌아갑니다. 대화 상자를 열려면 결과 파일 열기를선택하십시오.
2.전체 출력 파일을 보려면 사용자 정의단일 선택 단추를 선택하십시오. 전체 출력 파일에는 prpgrf.*파일과 파일이 포함됩니다 flsgrf.*. 시뮬레이션이 실행되었으므로 전 처리기 출력 파일이 삭제되어 flsgrf파일에 통합되었습니다.
3.flsgrf.Flow_Over_A_Weir대화 상자 에서 파일을 선택하고 확인을클릭하십시오.
이제 분석 탭이 표시됩니다. 시뮬레이션 결과를 시각화 하는 방법에는 여러 가지가 있습니다. 사용 가능한 플롯 유형은 다음과 같습니다.
·Custom : 이 매뉴얼 의 FLSINP 파일을사용하여플롯합니다. 사용자정의섹션의 출력 코드를 사용하여 출력 플롯을 수동으로 수정하는 데 사용할 수 있습니다. 이것은 고급 옵션입니다.
·프로브 :개별 셀, 경계, 구성 요소 및 도메인 전체(전역) 변수 대 시간에 대한 그래픽 및 텍스트 출력을 표시합니다. 자세한 내용은 프로브플롯및 프로브 : 특정시점의데이터와시간을 참조하십시오.
·1-D :셀 데이터는 X, Y 또는 Z 방향의 셀 라인을 따라 볼 수 있습니다. 플롯 제한은 공간 및 시간에 모두 적용할 수 있습니다. 자세한 내용은 1-D 플롯및 1-D : 라인을따른데이터대시간을 참조하십시오.
·2-D :셀 데이터는 XY, YZ 또는 XZ 평면에서 볼 수 있습니다. 플롯 제한은 공간 및 시간에 모두 적용할 수 있습니다. 속도 벡터 및 입자를 추가할 수 있습니다. 자세한 내용은 2 차원플롯및 2 차원 : 평면의데이터와시간의데이터를 참조하십시오.
·3-D :유체와 고체의 표면 플롯을 생성하고 셀 데이터로 채색 할 수 있습니다. 속도 벡터, 입자 (있는 경우) 및 유선과 같은 추가 정보를 추가할 수 있습니다. 플롯 제한은 공간 및 시간에 모두 적용할 수 있습니다. 자세한 내용은 3D 플롯및 3D : 표면의데이터대시간을 참조하십시오.
·텍스트 출력 :cell-by-cell 재시작, 선택 및 응고 데이터를 텍스트 파일에 쓸 수 있습니다. 자세한 내용은 텍스트출력및 텍스트 : ASCII 형식의공간데이터출력대시간을 참조하십시오.
·중립 파일 :재시작 및 선택된 데이터는 별도의 텍스트 파일에 정의 된 지정된 지점(보간 또는 셀 중심)에서 출력 될 수 있습니다. 자세한 내용은 중립파일 : 사용자정의좌표에서의공간데이터출력대시간을 참조하십시오.
·FSI TSE :유한 요소 유체 / 고체 상호 작용 및 열 응력 진화 물리학 패키지에서 출력됩니다. 자세한 내용은 FSI / TSE : 표면의구조데이터와시간을 참조하십시오.
3 차원 도표
1.Analyze -> 3-D탭을 선택하십시오.
2.Iso-surface = Fraction of fluid를선택하십시오. 이것은 표면을 그리는 데 사용되는 변수입니다. 선택한 등면변수에 대한 등고선 값기준을 충족하는 모든 셀을 통해 표면이 그려집니다. 유체의 분율이기본값이며 유체 표면이 표시됩니다.
3.색상 변수 = 압력을선택하십시오. 이 선택은 등위면의 색을 지정하는 데 사용되는 변수를 결정합니다 (이 경우 유체 표면은 압력에 의해 색이 그려집니다).
4.Component iso-surface overlay = Solid volume을선택하십시오. 솔리드 볼륨은 유체와 함께 솔리드 구성 요소를 표시합니다. 이전 단계에서는 체적 분수의 보완을등위면으로 선택하여 이 작업을 수행했지만 이 옵션을 사용하면 유체와 고체 표면을 동시에 플롯 할 수 있습니다.
5.이동 시간 프레임의최소 및 최대 위치들 (0 내지 1.25 초)에 슬라이더 위치.
6.렌더버튼을 클릭하여 디스플레이탭으로 전환하고 t = 0.0에서 1.25 초 사이에 일련의 11 플롯을 생성하여 압력에 의해 채색된 유체 표면과 위어 구조를 보여줍니다. 데이터 다시 시작이 선택되었으므로 11 개의 플롯이 있습니다.
7.사용 가능한 플롯이 사용 가능한 시간 프레임목록에 나열됩니다. 다음을 클릭하여 시간 프레임 사이를 이동하거나 시간 프레임을 두 번 클릭하여 표시하십시오. 첫 번째 및 마지막 시간 프레임은 다음과 같아야 합니다.
8.Analyze -> 3-D 탭으로돌아가서 Data Source그룹에서 Selected data라디오 버튼을 선택하십시오.
9.시간 프레임선택기의 두 슬라이더가 모두 오른쪽에 있으므로 마지막 시간 프레임 만 생성됩니다. 사용 가능한 시간 프레임이 많고 렌더링하는데 시간이 오래 걸리므로 선택한 데이터를 선택하면 인터페이스에서 자동으로 수행됩니다. 사용 가능한 모든 시간 프레임을 렌더링 하려면 왼쪽 슬라이더를 Time Frame Min = 0 으로이동하십시오.
10. 렌더링버튼을 클릭하십시오. 몇 초 안에 뷰가 디스플레이창으로 전환되고 101 개의 플롯이 사용 가능한 시간 프레임목록에 나열됩니다. 시간 프레임 사이를 이동하려면 다음을반복해서 클릭하십시오.
대칭 흐름 표시
위어 중심 아래로 대칭 평면을 사용하여 시뮬레이션을 설정했으므로 위어 구조의 절반만 시뮬레이션되고 표시됩니다. 프리젠테이션 목적으로 대칭 모델의 두 반쪽을 모두 표시할 수 있습니다.
1.아래와 같이 Analyze -> 3-D탭으로 돌아가서 Open Symmetry Boundaries확인란을 선택하십시오.
2.렌더링을클릭하십시오. 유체 표면이 디스플레이탭의 대칭 경계에서 열린 상태로 나타납니다.
3.화면 위의 도구 모음 메뉴에서 도구 -> 대칭을 선택하십시오.
4.대화 상자에서 Y 방향확인란을 선택하여 Y = 0 평면에서 결과를 미러링합니다.
5.적용및 닫기를선택하십시오.
6.마지막 시간 프레임을 두 번 클릭하십시오. 디스플레이는 아래와 같이 전체 위어 구조를 보여줍니다.
3 차원 애니메이션 만들기
다음 단계는 3 차원 유체 표면의 애니메이션을 만드는 것입니다. 애니메이션은 사용 가능한 시간 프레임 목록의 프레임에서 만든 동영상입니다. 애니메이션의 시각적 효과를 향상시키려면 모든 프레임에 공통 색상 스케일을 적용하는 것이 좋습니다.
1.분석 -> 3-D탭으로 돌아갑니다.
2.윤곽 제한그룹 상자에서 전역라디오 버튼을 모두 선택하십시오.
3.렌더를클릭 하여 다시 그리고 디스플레이탭으로 돌아갑니다.
4.도구 -> 대칭 -> Y 방향 -> 적용선택을 반복하여 Y = 0 평면에서 결과를 반영합니다.
5.선택 도구 -> 애니메이션 ->러버 밴드 캡처를다음과 같이 선택 확인 Mesh지가 나타납니다 그것을 읽은 후.
6.마우스 왼쪽 버튼을 클릭 한 상태에서 드래그하여 애니메이션을 적용할 화면 부분을 선택하십시오. 선택한 영역 주위에 선택 상자가 나타납니다.
7.디스플레이 창 위에서 빨간색 캡처버튼을 선택하십시오. 애니메이션을 시작하는 대화 상자가 나타납니다.
8.애니메이션의 기본 이름은 out.avi입니다. 아래에 표시된 것처럼 보다 구체적인 이름이 권장됩니다.
9.기본 프레임 속도는 초당 10 프레임입니다. 이 시뮬레이션의 마감 시간은 1.25 초이고, 일정한 시간 간격으로 100 개의 플롯이 있으므로 ‘실제‘속도는 초당 80 프레임입니다. 너무 빠를 수 있으므로 대신 5를입력 하고 확인을누르십시오.
10. 각 시간 프레임이 표시창에 렌더링 되고 비트 맵 파일이 시뮬레이션 디렉토리에 작성됩니다. 이 프로세스가 완료되면 다음 대화 상자가 나타납니다.
프로세스의 다음 단계를 시작하려면 확인버튼을 클릭하십시오. 새로운 프로세스 (BMP2VAI.exe)가 시작되고 압축 방법을 선택할 수 있는 새로운 비디오 압축창이 나타납니다. 다른 창 뒤에 숨겨져 있으면 앞으로 가져옵니다.
애니메이션의 기본 압축은 압축되지 않습니다. 파일 크기가 너무 커서 뷰어에 로드 할 수 없으므로 대부분의 애니메이션에는 권장되지 않습니다. Windows를 사용하는 경우 Microsoft Video 1을, Linux를 사용하는 경우 Cinepak을선택하십시오. 여기에서 선택하는 것은 컴퓨터에서 사용할 수 있는 비디오 코덱과 비디오를 표시하는 데 사용하는 기계에서 사용할 수 있는 것입니다.
애니메이션 속도가 데이터 속도에 의해 제한되지 않도록 데이터 속도확인란을 선택 취소하십시오.
압축 프로세스를 시작하려면 확인을 클릭하십시오. 압축이 완료되면 다음 대화 상자가 나타납니다.
확인을클릭하십시오. 애니메이션 프로세스가 완료되었습니다.
Windows 탐색기에서 .avi파일을 찾는 가장 빠른 방법 은 시뮬레이션 관리자탭으로 이동하여 시뮬레이션 입력 파일링크를 클릭하는 것 입니다.
.avi파일 을 두 번 클릭하여 애니메이션을 재생 하십시오. 이전에 선택한 압축 형식을 읽을 수 있는 올바른 코덱이 설치되어 있지 않으면 오픈 소스 다중 코덱 비디오 플레이어 설치를 고려하십시오.
2 차원 도표
1.Analyze -> 2-D탭을 선택하십시오. 이 시뮬레이션의 결과를 보는 데 가장 유용한 평면은 평면 Y = 0.0에있는 위어 중심선의 XZ 평면입니다.
2.XZ 평면라디오 버튼을 선택하십시오.
3.Y 제한 슬라이더를 모두 Y = 0.25 (Y = 0.0에 가장 가까운 셀 중심 y 좌표)로 드래그 합니다. 또한 동일한 위치가 J = 2로 식별되어 해당 셀이 도메인에서 두 번째임을 나타냅니다. 첫 번째 셀 (J = 1)은 Mesh 외부에 있으며 경계 조건 속성을 계산하는 데 사용됩니다. 기본 윤곽 변수는 압력이며 기본 속도 벡터는 기본적으로 선택됩니다. 솔리드 형상은 모든 2D 플롯과 함께 자동으로 표시되므로 3D 플롯과 같이 활성화 할 필요가 없습니다.
4.벡터 옵션을클릭하고 X = 2및 Z = 2를입력하십시오. 벡터는 이제 다른 모든 셀에 플롯 됩니다. 벡터 옵션을 적용하려면 확인을 선택하십시오.
5.Y = 0 평면에서 2 차원 압력 플롯의 시간 시퀀스를 생성하려면 렌더링을 클릭하십시오. T = 0.0 초 (왼쪽) 인 다음과 유사한 그래픽이 나타납니다. T = 0.125 초 (중간); 그리고 T = 1.25 초 (오른쪽).
6.디스플레이화면의 오른쪽 상단에 있는 형식버튼을 선택하십시오.
7.선 색상, 벡터 길이 및 화살촉 크기 변경과 같은 다양한 옵션을 시험해보십시오. 변경 사항을 보려면 적용을 선택하십시오.완료되면 재설정및 확인을 선택하여 기본 설정으로 돌아가서 대화 상자를 닫습니다. 모든 플롯에 대해 선호하는 옵션 세트가 있는 경우 저장버튼을 선택하여 저장할 수 있습니다.
1 차원 도표
분석 -> 1-D탭을 선택하십시오. 이 탭에서는 하나 이상의 플롯 시간에서 셀 행을 따라 압력, 유체 깊이, 유체 상승 및 속도와 같은 셀별 출력 변수의 꺾은 선형 차트 플롯을 사용할 수 있습니다.
데이터 소스로 선택을 선택합니다. 사용 가능한 변수는 이제 더 빈번한 플로팅을 위해 선택된 변수 만 표시합니다.
자유 변수 표고를 데이터 변수로 선택하십시오. 유압 데이터는 출력탭에서 선택되었으므로 사용할 수 있습니다.
이 시뮬레이션의 흐름 방향은 주로 x 축과 평행하므로 X 방향을선택하십시오.
Y 방향 슬라이더를 0.25(J = 2)로 이동하여 Y 방향에서 흐름 중심선에 가장 가까운 셀이 표시됩니다.
기본적으로 전체 X 범위가 표시됩니다. 플롯의 범위를 제한하려는 경우 X 방향슬라이더를 이동할 수 있습니다. Z 방향슬라이더의 위치는 주어진 x, y 위치에서 z 셀의 각 열에 대해 하나의 자유 표면 높이만 기록되므로 중요하지 않습니다.시간 프레임 슬라이더는 0초와 1.25초여야 합니다.
렌더링을클릭하십시오. t = 0.0에서 t = 1.25s까지의 시리즈 플롯이 디스플레이탭의 플롯 목록에 나열됩니다. 이러한 플롯을 볼 수 있는 여러 가지 모드가 있습니다. 기본 모드는 단일 모드이며 형식버튼 아래의 드롭 다운 상자에 표시됩니다.
다양한 시간에 유체 표면 높이의 플롯을 비교하려면 드롭 다운 상자에서 오버레이 모드를 선택하십시오.
오른쪽 창에서 플롯 1, 13및 101을 선택하려면 클릭하십시오. 플롯 이름에는 또한 기록된 시간이 표시됩니다 (t = 0.0, 0.15s 및 1.25 초). 출력은 아래와 같이 나타납니다.
이 플롯을 비트 맵 또는 포스트 스크립트 파일에 저장하려면 출력버튼을 선택하십시오.
확인 화면에 플롯오버레이 플롯을 캡처하는 확인란을 (그리고 단 하나의 출력 파일을).
쓰기버튼을 선택하여 이미지 파일을 만듭니다.
결과 이미지 파일은 시뮬레이션 디렉토리에 있으며 (시뮬레이션 관리자탭 에서이 파일을 찾는 방법을 기억하십시오) 이름이 지정한 plots_on_screen.bmp됩니다.
프로브 플롯
1. 분석 -> 프로브탭을 선택하십시오. 시간 기록 플롯은이 탭에서 변수 대 시간의 라인 그래프 또는 텍스트 출력으로 생성됩니다. FLOW-3D에는 데이터 소스그룹에서 선택되는 세 가지 유형의 시간 종속 데이터가 있습니다.
·공간 데이터 :재시작및 선택된 데이터소스. 단일 x, y, z 셀 중심 좌표의 시간 종속 값이 표시됩니다. 값은 시간과 관련하여 통합되거나 시간과 관련하여 차별화되거나 이동 평균 (시간)으로 통합될 수 있습니다.
·일반 history 데이터 :.글로벌 수량은 시간에 따라 다릅니다. 일반적인 양은 평균 운동 에너지, 시간 단계 및 대류 볼륨 오류입니다. 또한 이 데이터 유형에는 모델 설정 -> 메싱 및 지오메트리 탭에서 이러한 옵션을 선택한 경우 지정된 측정 위치(배플, 샘플링 볼륨, 히스토리 프로브)의 모든 데이터와 이동 또는 정지 상태의 솔리드 및 스프링/로프를 위한 통합 출력이 포함됩니다.
·Mesh-dependent data :메쉬 경계에서 시간에 따른 수량(계산 또는 사용자 지정)입니다. 일반적인 수량은 경계에서의 유량 및 경계에서의 지정된 유체 높이입니다.
2.데이터 원본에서 일반 기록 라디오 버튼을 선택합니다. X, Y 및 Z 데이터 점 슬라이더가 회색으로 바뀝니다. 이는 일반 기록 데이터가 특정 셀과 연결되어 있지 않기 때문입니다.
3.목록에서 질량 평균 유체 평균 운동 에너지를 선택하십시오.
4. 단위를 선택하여 플로팅 단위 대화 상자를 엽니다.
5. 플롯에 단위 표시를 선택하십시오.
6. SI, CGS, slugs/feet/seconds 또는 pounds/inches/seconds를 선택하여 원하는 단위 시스템으로 결과를 변환하고 출력합니다. 장치를 표시하고 변환하려면 모델 설정 -> 일반 탭에서 장치 시스템을 선택해야 합니다. 이전 단계에서 이 항목을 확인했으며, 지오메트리 및 유체 특성은 centimeters/grams/seconds 시스템에서 지정되었습니다.
7.Plotting Units 대화 상자를 닫으려면 OK를 선택하십시오.
8.데이터의 그래픽 출력을 생성하려면 렌더를 선택하십시오. 출력은 시간에 따른 영역의 모든 유체에 대한 질량 평균 평균 운동 에너지를 보여줍니다. 이전 단계에서 선택한 사항에 따라 단위 레이블과 함께 그림이 나타납니다. 플롯은 총 운동 에너지가 일부 평균값 주위에서 진동하고 있음을 나타냅니다. 진동이 작아짐에 따라 시뮬레이션은 정상 상태 흐름에 접근합니다.
9.분석 -> 프로브 탭으로 돌아갑니다.
10. 출력 양식 그룹에서 텍스트를 선택하여 그래프를 텍스트 데이터로 출력한 다음 렌더링을 다시 선택하십시오.
11. 나타나는 텍스트 대화 상자에서 다른 이름으로 저장 버튼을 선택하여 출력을 텍스트 파일로 저장할 수 있습니다.
12. 출력 창을 닫으려면 계속을 선택하십시오.
텍스트 출력
1.Analyze -> Text Output 탭을 선택하십시오.
2.텍스트 출력 은 셀별 데이터 ( 다시 시작 또는 선택됨 ) 만 출력 할 수 있고 (구성 요소, 측정 스테이션 또는 글로벌 데이터 없음) 둘 이상의 셀을 선택할 수 있다는 점을 제외하고 프로브 탭 과 동일한 방식으로 작동합니다. 각 플롯 시간에 대한 출력 데이터. 셀은 슬라이더를 사용하여 3D 블록에서 선택됩니다. 기본 공간 범위는 전체 도메인으로 설정됩니다.
FLOW-3D의 침전물 이송 모델을 사용하여 세굴 및 침전물을 평가할 수 있으며, 여기서 3차원 유량 구성 요소가 세굴 프로세스를 주도하고 있습니다. Flow-3D의 유체역학 모델은 유체물리학을 설명하는 정전기적이지 않은 레이놀즈-평균화된 Navier-Stokes 방정식을 완벽하게 해결합니다. 유체역학적 솔버는 침전물 운반 모듈과 완전히 결합되어 있어 침전물 운반 및 비접착 토양의 부유식 침식, 인포테인먼트 및 침식을 시뮬레이션합니다(Wei et al., 2014). 베드로드, 인포테인먼트 및 정착 프로세스에 사용되는 모든 경험적 관계는 완전히 사용자 정의 가능하며, 최대 10개의 침전물 종(곡물 크기, 질량 밀도, 임계 전단 응력 등 서로 다른 특성을 가진)을 정의할 수 있습니다. FLOW-3D는 짧은 경과 시간 척도에 대한 국부적 스쿠어를 시뮬레이션하는 데 이상적입니다.
FLOW-3D‘s Sediment Transport model can be used to evaluate scour and deposition, where three-dimensional flow components are driving the scouring process. FLOW-3D’s hydrodynamic model solves the full unsteady non-hydrostatic Reynolds-averaged Navier-Stokes equations that describe the flow physics. The hydrodynamic solver is fully coupled with a sediment transport module that simulates bedload and suspended sediment transport, entrainment and erosion for non-cohesive soils (Wei et al., 2014). All empirical relationships used in bedload, entrainment and settling processes are fully customizable, and up to 10 different sediment species (with different properties such as grain size, mass density and critical shear stress) can be defined. FLOW-3Dis ideal for simulating local scour over short episodic time scales.
Modeling Capabilities – Unsteady 3D mobile bed modeling – Bedload and suspended sediment transport – Non-cohesive sediment – 10 individual grain size fractions – Suspended sediment settling and entrainment – Critical angle of repose
Applications – River and coastal morphodynamics – Bridge pier and abutment scour – Local scour at hydraulic structures – Sedimentation basins – Reservoir flushing
Sediment Transport Model
Sentral Transport 모델은 8.0 버전(Brethour, 2009년)에서 처음 도입되었으며, 11.1 버전(Wei et al., 2014년), 가장 최근에는 12.0 버전(Flow Science, 2019년)에서 광범위한 개정을 거쳤습니다. 숫자 모델에서 시뮬레이션된 물리적 프로세스의 개략도가 아래에 나와 있습니다.
The Sediment Transport model was first introduced in version 8.0 (Brethour, 2009), and has gone through extensive revisions in version 11.1 (Wei et al., 2014), and most recently in version 12.0 (Flow Science, 2019). A schematic of the physical processes simulated in the numerical model is illustrated below.
수치 모델에서 침전물은 포장된 Bed로서 일시 중단된 상태로 존재할 수 있습니다. 포장된 Bed는 PRIPT™ 기법을 사용하여 복잡한 솔리드 경계(Hirt 및 Sicilian, 1985)에 표현된 지울 수 없는 솔리드 객체입니다. 이것은 유체역학 용해기의 고체 물체를 나타내는 데 사용되는 방법과 동일합니다. 포장된 Bed의 형태학적 변화는 침전물 질량의 보존에 의해 좌우됩니다.
In the numerical model, sediment can exist as packed bed and in a suspended state. A packed bed is an erodible solid object that is represented using the FAVOR™ technique for complex solid boundaries (Hirt and Sicilian, 1985). This is the same method used to represent solid objects in the hydrodynamic solver. The morphological change in the packed bed is governed by the conservation of sediment mass.
형태학적 변경은 모형에 숫자로 표시되는 여러 가지 물리적 프로세스에 의해 제어됩니다. 이러한 프로세스에는 베드로드 운송, 인포테인먼트 및 증착이 포함됩니다. 베드로드 이송은 침전물이 서스펜션에 전달되지 않고 채널을 따라 횡방향으로 이동하는 물리적 과정입니다. 인포테인먼트란 난류 에디가 패킹 베드 상단의 곡물을 제거하고 일시 중단된 상태로 전환하는 과정입니다. 포장이란 곡물이 현수막에서 안착되어 포장된 침대에 퇴적하는 과정입니다. 수치 모델에서 이것은 일시 중단된 상태에서 포장된 베드 상태로의 전환입니다.
The morphological changes are governed by several different physical processes that are represented numerically in the model. These processes include bedload transport, entrainment and deposition. Bedload transport is the physical process of sediment moving laterally along the channel without being carried into suspension. Entrainment is the process by which turbulent eddies remove the grains from the top of the packed bed and transition to the suspended state. Packing is the process of grains settling out of suspension and depositing onto the packed bed. In the numerical model, this is the transition from the suspended to the packed bed state.
인포테인먼트 및 패킹의 상대적 비율은 포장된 베드와 부유 상태 사이의 침전물 질량 교환을 제어합니다. 이 모델은 Meyer-Peter Müler(1948), Nielsen(1992) 또는 Van Rijn(1984)의 방정식을 사용하여 베드 인터페이스가 포함된 각 메시 셀에서 베드로드 전송을 계산합니다. 메쉬 셀에서 이웃의 각 메쉬 셀로 이동하는 곡물의 양을 결정하기 위해 하위 메쉬 방법이 사용됩니다. 인포테인먼트에서 곡물의 리프팅 속도는 Winterwerp 등(1992)의 방정식을 사용하여 계산됩니다. 안착 속도는 Soulsby(1997년)를 사용하여 계산됩니다. 베드 인터페이스가 포함된 메시 셀에서 인터페이스의 위치, 방향 및 면적을 계산하여 베드 전단 응력, 무차원 전단 응력, 베드로드 전송 속도 및 인포테인먼트 속도를 결정합니다. 3D 난류 흐름의 베드 전단 응력은 표준 벽 함수를 사용하여 중간 곡물 크기에 비례하는 베드 표면 거칠기를 고려하여 평가됩니다.
The relative rates of entrainment and packing control the exchange of sediment mass between the packed bed and suspended states. The model calculates bedload transport in each mesh cell containing the bed interface using the equation of Meyer-Peter Müller (1948), Nielsen (1992) or Van Rijn (1984). A sub-mesh method is employed to determine the amount of grains moving from the mesh cell into each mesh cell in its neighbor. The lifting velocity of grains in entrainment is calculated using the equation of Winterwerp et al. (1992). The settling velocity is calculated using Soulsby (1997). In the mesh cells containing the bed interface, location, orientation and area of the interface are calculated to determine the bed shear stress, dimensionless shear stress, bedload transport rates and entrainment rates. Bed shear stress in 3D turbulent flows is evaluated using the standard wall function with consideration of bed surface roughness that is proportional to the median grain size.
부유된 침전물은 유체의 스칼라 질량 농도로 표시됩니다. 농도는 주어진 셀에서 균일한 것으로 가정되며 유체 셀 밀도 및 점도와 결합됩니다. 각 종에 대해, 부유 침전물 농도는 수송 방정식을 풀어서 계산됩니다.
The suspended sediment is represented as a scalar mass concentration in the fluid. The concentration is assumed to be uniform in a given cell and is coupled with the fluid cell density and viscosity. For each species, the suspended sediment concentration is calculated by solving a transport equation.
Validations
다음 5가지 검증 사례는 실험 데이터와 FLOW-3D의 침전물 이송 모델의 시뮬레이션 결과를 비교합니다.
마오(1986년) Mao는 수중 수평 파이프라인 아래 침대의 무서운 프로파일을 얻기 위해 실험 작업을 수행했습니다. 아래 그림은 FLOW-3D를 사용하여 얻은 결과와 실험 결과를 비교합니다.
Chatterjee et al. (1994)
수평 제트 침수로 인해 국부적인 스쿠어 프로파일을 얻기 위한 실험 작업이 수행되었습니다. 아래 그림은 scour구멍 깊이와 둔부 높이에 대한 실험 대 FLOW-3D의 숫자 결과를 시간의 함수로 비교합니다. 이 애니메이션은 scour구멍과 둔부 높이가 최대 1시간 내에 안정된 상태에 도달한다는 것을 보여줍니다.
Gladstone et al. (1998)
In these experiments the propagation and deposition patterns of particle-laden flows were studied. The plot below compares experimental versus FLOW-3D simulation results from three different setups, labeled case A (100% 0.025mm size particles), case D (50% 0.069mm and 50% 0.025mm size particles), and case G (100% 0.069mm size particles).
Faruque et al. (2006)
이 논문에서, 저자들은 실험을 통해 3차원 벽면 제트기를 물에 잠기게 함으로써 국부적인 악취를 연구했습니다. 아래 표는 세 가지 서로 다른 테일워터 비율에 대한 scour 구멍의 3D 형태학적 변화에 대한 실험과 FLOW-3D 수치 결과를 비교합니다.
In this paper, the authors studied local scour by submerged three-dimensional wall jets via experiments. The table below compares the experimental versus FLOW-3D numerical results for 3D morphological changes in the scour hole for three different tailwater ratios.
References
Brethour, J.M., Hirt, C.W., 2009, Drift Model for Two-Component Flows, FSI-14-TN-83, Flow Science, Inc.
Chatterjee, S.S., Ghosh, S.N., and Chatterjee M., 1994, Local scour due to submerged horizontal jet, Journal of Hydraulic Engineering, 120(8), pp. 973-992.
Faruque, M.A.A., Sarathi, P., and Balachandar R., 2006, Clear Water Local Scour by Submerged Three-Dimensional Wall Jets : Effect of Tailwater Depth, Journal of Hydraulic Engineering, 132(6), pp. 575-580.
Flow Science, 2019, FLOW-3D Version 12.0 User Manual, Santa Fe, NM: Flow Science, Inc. https://www.flow3d.com
Gladstone, C., Phillips, J.C., and Sparks R.S.J., 1998, Experiments on bidisperse, constant-volume gravity currents: propagation and sediment deposition, Sedimentology 45, pp. 833-843.
Hirt, C.W. and Sicilian, J.M., 1985, A porosity technique for the definition of obstacles in rectangular cell meshes, 4th International Conference on Numerical Ship Hydrodynamics, Washington, D.C.
Khosronejad, A., Kang, S., & Sotiropoulos, F., 2012. Experimental and computational investigation of local scour around bridge piers, Advances in Water Resources, 37, pp. 73-85.
Mao, Y., 1986. The interaction between a pipeline and an erodible bed, PhD thesis, Institute of Hydrodynamics and Hydraulic Engineering, Technical University of Denmark, Lyngby, Denmark.
Meyer-Peter, E. and Müller, R., 1948, Formulas for bed-load transport, Proceedings of the 2nd Meeting of the International Association for Hydraulic Structures Research. pp. 39–64.
Nielsen, P., 1992, Coastal bottom boundary layers and sediment transport (Vol. 4). World scientific.
Soulsby, R., 1997, Dynamics of Marine Sands, Thomas Telford Publications, London.
Van Rijn, L. C., 1984, Sediment Transport, Part I: Bed load transport, Journal of Hydraulic Engineering 110(10), pp. 1431-1456.
Wei, G., Brethour, J.M., Grüenzner M., and Burnham, J., 2014, The Sediment Scour Model in FLOW-3D, Technical Note FSI-14-TN-99, Flow Science, Inc.
Winterwerp, J.C., Bakker, W.T., Mastbergen, D.R. and Van Rossum, H., 1992, Hyperconcentrated sand-water mixture flows over erodible bed, Journal of Hydraulic Engineering, 118(11), pp. 1508–1525.
FLOW-3D 는 완전 3 차원 유동, 2 차원 깊이 평균 (천수(shallow water)) 흐름 및 3 차원/2 차원 깊이 혼합 평균 흐름을 위한 포괄적인 난류 모델링 제품군을 제공합니다.
그림1. 부두 하류에서의 와류 검출을위한 Q- 기준의 3D지도
FLOW-3D 에는 8 가지 난기류 옵션이 있습니다.
Prandtl 혼합 길이 모델은 3 차원 난류 효과를 설명하기 위한 가장 초기의 시도 중 하나입니다. 가장 복잡한 모델이 아니며 더 이상 널리 사용되지 않습니다. FLOW-3D 는 주로 학술 연구에서의 유용성을 포함합니다.
소위 1 방정식 모델은 난기류를 나타내는 초기 노력이기도 합니다. 시간 평균 난류 운동 에너지 k를 계산하고 모든 위치에서 알려진 난류 혼합 길이 LT 가 필요합니다. LT 는 일반적으로 미리 알려지지 않기 때문에, one-equation 모델은 복잡한 유량을 모델링하는 데 적합하지 않습니다.
표준 k-ε 모델 (Harlow & Nakayama 1967)은 난류 운동 에너지 k와 소산 속도 ε를 계산하고 난류 혼합 길이 LT를 동적으로 찾는 2 방정식 모델입니다. 이것은 업계 표준이며 광범위한 흐름을 표현하는데 유용하다는 것이 발견되었습니다 (Rodi 1980).
재 정규화 그룹 (RNG) k-ε 모델 (Yakhot & Orszag 1986, Yakhot & Smith 1992)은 2 방정식 k-ε 모델의 보다 견고한 버전이며, 대부분의 산업에서의 문제에 권장됩니다. 표준 k-ε 모델의 기능을 확장하여 과도기 난류, 곡선 흐름, 벽 열 전달 및 물질 전달의 더 나은 적용 범위를 제공합니다.
k-ω 2 방정식 모델 (Wilcox 1988, 1998, 2008)은 두 번째 변수를 난류 소산 ε이 아니라 ω ≡ ε / k로 정의한다 (Kolmogorov 1942). Wilcox는 1988 년부터 k-ω 2 방정식 모델을 개선했으며 1998 년에는 자유 전단 유동에 대한 모델의 정확성을 크게 개선 한 새로운 계수를 도입했습니다. FLOW-3D 의 k-ω 2 방정식 모델은 제트, 후류 및 플럼을 퍼 뜨리는 것과 같은 유선형 압력 구배를 갖는 자유 전단 흐름을 모델링하는 데 적합합니다. LES 모델은 평균 난류 운동 에너지를 나타내기 위해 스칼라를 사용하지 않고 대부분의 난류 변동을 직접 해결합니다. 그것은 2 방정식 모델보다 훨씬 더 미세한 메쉬 해상도를 필요로하며 난기류에 대한 보다 광범위한 통계를 제공합니다.
2-D 심도 평균 얕은물 난류 모델은 대수적인 완전 난류 속도를 가정합니다. 첫 번째 옵션은 일정한 항력 계수 CD를 가정하며, 이는 공간적으로 변화 할 수 있습니다.
두 번째 2-D 심도 평균 천수(shallow water) 난류 모델은 항력 계수 CD 를 유체 깊이와 공간적으로 가변되는 표면 거칠기의 동적 함수로 만듭니다.
완전한 3-D 테스트를 통해 LES 모델 출력을 시간 평균화하면 2 방정식 Reynolds REN (Reynolds Averaged Navier-Stokes) 모델 (표준 k-ε, RNG k-ε 및 k- ω).
아래의 물고기 통로 비디오에 나와 있습니다.
난류 시뮬레이션 – 모델 비교 / Turbulence Simulations – A Model Comparison
첫 번째 비디오에서는 FLOW-3D 의 LES (Large-Eddy Simulation) 난류 모델을 사용하여 어류 통과를 시뮬레이션하여 속도 변동의 크기를 분석합니다. 두 번째 비디오는 동일한 시뮬레이션의 시간 평균 결과를 보여줍니다. 여기에서 간단한 사용자 정의는 시간 평균 LES가 레이놀즈 평균 Navier-Stokes (RANS) 난류 모델 결과와 매우 유사하다는 것을 보여줍니다. 세 번째 비디오는 RNG (Renormalized Group) k-ε 난류 모델을 사용하여 시연하기 위해 동일한 시뮬레이션을 사용합니다. RNG k-ε 모델은 대부분의 레이놀즈 평균 Navier-Stokes (RANS) 난류 모델과 마찬가지로 속도 변동을 등방성 스칼라 값으로 처리하여 시간에 따른 속도 변동을 감쇠시킵니다. 결과는 두 번째 비디오에서 볼 수 있듯이 직접 LES 결과를 시간 평균하여 찾은 결과와 유사합니다.
참고 문헌
Driver, DM and Seegmiller, HL, 1985, AIAA Journal (23), 163-171의 다양한 채널 유동에서 재 부착하는 난류 전단 층의 특징 .
Harlow, FH and Nakayama, PI, 1967, 난류 수송 방정식 , 유체 역학 (10), 2323-2332.
Harlow, FH 및 Nakayama, PI, 1968, 난기류 에너지 감쇠율의 전송 , Los Alamos Scienti fi c 실험실 보고서 LA-3854.
Kolmogorov, AN, 1942, 비압축성 유체에서의 난류 운동 방정식 , Izvestia Academy of Sciences, 소련; Physics (6), 56-58.
Pope, S. B, 2000, Turbulent Flows , Cambridge University Press.
Rodi, W., 1980, 난류 모델과 유압 장치의 적용 : 최첨단 검토 , 국제 유압 연구 협회 (IAHR), 델프트, 네덜란드.
Saffman, PG, 1970, Inhomogeneous Turbulent Flow의 모델 , Royal Society London A (317), 417-433의 절차.
Speziale, CG, Abid, R., and Anderson, EC, 1992, AIAA 저널 (30), 324-331, 벽 근처 난류에 대한 2 방정식 모델의 중요성 평가.
Wilcox, DC, 1988, AIAA Journal (26), 1299-1310의 진보 된 난류 모델에 대한 스케일 결정 방정식의 재평가.
FLOW-3D를 이용하여 수처리분야의 엔지니어가 정확하게 접촉조 모델링을 할 수 있습니다. 스칼라 전송 모델, 파티클 모델과 화학 모델은 차세대 시각화 도구인 FlowSight를 이용하여, 엔지니어가 정확하게 혼합 유압 효율을 결정하는 데드 존과 재순환 영역의 위치를 시각화할 수 있습니다.
FLOW-3D는 FAVOR™ 방법을 사용하므로, 형상 및 구조 요소가 수정된 형상을 re-meshing하는 어려움이 없습니다. 새로운 자동 등격자 생성 기능으로, 자동 격자 생성으로 일반 직각 격자 방식의 사용 편리성을 가지면서도 형상을 정확하게 구현할 수 있습니다.
이것은 수처리 시설물의 소독조 유동 해석 결과입니다. 염소 농도는 체적 측정으로 시간에 따른 농도 C (t)와 함께, 측정 위치에서의 단위 시간당 E (t)농도의 변화를 보여줍니다.