747 research outputs found

    Three-Dimensional Smoothed Particle Hydrodynamics Simulation for Liquid Droplet with Surface Tension

    Get PDF
    We provide a basic method of Smoothed Particle Hydrodynamics (SPH) to simulate liquid droplet with surface tension in three dimensions. Liquid droplet is a simple case for surface tension modeling. Surface tension works only on fluid surface. In SPH method, we simply apply the surface tension on the boundary particles of liquid. The particle on the 3D boundary was detected dynamically using Free-Surface Detection algorithm. The normal vector and curvature of the boundary surface were calculated simultaneously with 3D boundary surface reconstruction using Moving Least-Squares (MLS) method. Before the reconstruction, the coordinate system was transformed into a local coordinate system. Afterwards, the surface tension force which depends on curvature of the surface, was calculated and applied on the boundary particles of the droplet. We present the simulation result of droplet motion with gravity force. By using the basic method of SPH for fluid modeling, and a combination of 3D Free-Surface Detection algorithm with MLS method, we can simulate droplet phenomena with expected result.Comment: 9 pages, 4 figure, Selected Paper from the International Symposium on Computational Science 201

    λ‹€μƒμœ λ™ 및 고체 파편물의 수λ ₯학적 거동 해석을 μœ„ν•œ GPU 기반의 SPH-DEM 연계해석 μ½”λ“œ 개발

    Get PDF
    ν•™μœ„λ…Όλ¬Έ (박사) -- μ„œμšΈλŒ€ν•™κ΅ λŒ€ν•™μ› : κ³΅κ³ΌλŒ€ν•™ μ—λ„ˆμ§€μ‹œμŠ€ν…œκ³΅ν•™λΆ€, 2020. 8. κΉ€μ‘μˆ˜.In the late phase severe accident of LWR, the massive corium releases out of the reactor pressure vessel (RPV) and falls to the coolant if the In-Vessel Retention (IVR) strategy fails. The melt jet can be fragmented into debris particles based on the assumption that the ex-vessel pool is sufficiently deep. It is known that there are various three-phase flow issues associated with the fragmented debris particles under the influence of phase change of cavity coolant. In such cases, the vapor phase forms a sharp and dynamic interface with the liquid phase while the transient relocation behavior of debris particles is the main concern. Thus, coupling Lagrangian-based multi-phase CFD techniques and Discrete Element Method (DEM) can be an effective approach in terms of numerical modeling of such behaviors. In this respect, an integrated numerical code for incompressible 3-phase flow has been developed in this study by two-way phase coupling of multi-phase Smoothed Particle Hydrodynamics (SPH) and DEM model. Smoothed Particle Hydrodynamics (SPH) is one of the best-known meshless CFD methods in which the fluid system is represented as the finite number of Lagrangian particles. The SPH code developed in this study proposes a new density estimation model and improves the surface tension model for accurate simulation of incompressible two-phase flow behavior. The demonstration of its applicability has been performed through several V&V simulations including multi-phase dam-break and sloshing simulations. Discrete Element Method (DEM) is a direct simulation method for a rigid body that can analyze the translation, rotation, and collision behavior of solid particles in detail. The soft-sphere collision model with Hertz-Mindlin contact force equations is adopted for developed DEM code in this study. To precisely estimate the wall boundary interactions of bed-formed debris particles, a versatile wall boundary model is newly proposed in this study that also covers the sliding and rolling behavior of solid particles. The inter-particle collision behavior and sliding & rolling motion of particles are well proven in several V&V cases. The numerical code system for incompressible 3-phase flow is newly developed by two-way phase coupling of the above two models (SPH-DEM). The unresolved coupling approach between two methods was adopted for the analysis of the overall behavior of particulate solid debris. The fundamental validation of the phase coupled model was performed for both single-particle behavior and particulate granular flow such as dam-breaking motion of particle-fluid. The SPH-DEM coupled code in this study has been parallelized based on Graphical Process Unit (GPU) in order to overcome the inherent efficiency problem of the Lagrangian-based numerical method. Parallel mapping and reduction are applied for solving discretized summation equations of each SPH particle, solving contact force equations for each DEM particle, and also for solving coupling equations between SPH and DEM particles. The efficiency of code parallelization was evaluated through the scalability analysis based on the benchmark calculation. Finally, the simulation of the vapor-driven leveling behavior of spherical solids was performed as a case study to demonstrate the applicability of the developed code. The time-variant surface shape of solid particles was compared with the benchmark experiments both qualitatively and quantitatively. The effect of gas flow rate on the tendency of leveling behavior also has been analyzed. The developed numerical system in this study is expected to be a good alternative for the simulation of such phenomena that were difficult to handle with traditional numerical methods since the numerical schemes used in the code have a high potential for simulation of complicated physics with highly deformable geometry. In addition, this validated code system can contribute to hydrodynamic modeling studies for severe accident technology by performing numerical experiments on conditions that hard to be conducted on a laboratory scale.κ°€μ••κ²½μˆ˜λ‘œ μ€‘λŒ€μ‚¬κ³  ν›„κΈ° κ³Όμ •μ—μ„œ ν•΅μ—°λ£Œ 용육물 λ…Έλ‚΄ μ–΅λ₯˜ μ „λž΅μ˜ μ‹€νŒ¨λ‘œ 용육 ν•΅μ—°λ£Œκ°€ 용기 λ°”κΉ₯μͺ½ 침수 κ³΅λ™μœΌλ‘œ 재배치되면, κ³΅λ™μ˜ κΉŠμ΄λ‚˜ 제트의 속도 λ“±μ˜ 쑰건에 따라 ν•΅μ—°λ£Œ νŒŒνŽΈν™”κ°€ λ°œμƒν•  수 μžˆλ‹€. μ΄λŸ¬ν•œ ν•΅μ—°λ£Œ νŒŒνŽΈλ¬Όμ—μ„œ λ°œμƒν•˜λŠ” 뢕괴열에 μ˜ν•΄ λƒ‰κ°μž¬μ˜ 상변화가 이루어지면, μ΄μƒμœ λ™ μžμ—°λŒ€λ₯˜μ™€ 고체 파편물이 κ³΅μ‘΄ν•˜λŠ” 3μƒμœ λ™ μ‹œμŠ€ν…œμ΄ λœλ‹€. ν•΅μ—°λ£Œ 용육물의 제트 νŒŒμ‡„, ν•˜κ°• 및 퇴적, 평탄화 λ“± 일련의 κ³Όμ •μ—μ„œ λ‹€μ–‘ν•œ ν˜•νƒœμ˜ 파편물 3상 냉각 거동이 λ°œμƒν•  수 있으며, ν•΅μ—°λ£Œ 파편물의 뢕괴열이 μΆ©λΆ„νžˆ μ œκ±°λ˜μ§€ λͺ»ν•˜λ©΄ 파편물의 응집 및 재용육이 λ°œμƒν•  수 μžˆλ‹€. μ€‘λŒ€μ‚¬κ³  μ™„ν™”μ˜ κ΄€μ μ—μ„œ μ΄λŸ¬ν•œ 3상 냉각 거동에 λŒ€ν•œ μƒμ„Έν•œ 이해가 ν•„μš”ν•˜μ§€λ§Œ, μ΄μƒμœ λ™ μžμ—°λŒ€λ₯˜, μ΄μƒμœ λ™ 열전달, κ³ ν™”, 용육, λΉ„λ“±, μœ λ™-μž…μž μƒν˜Έμž‘μš© 등을 ν¬ν•¨ν•˜λŠ” 거동 자체의 λ³΅μž‘μ„±μœΌλ‘œ 인해 κ΄€λ ¨ ν˜„μƒλ“€μ˜ 예츑 및 ν‰κ°€λŠ” 큰 λΆˆν™•μ‹€μ„±μ„ λ‚΄ν¬ν•˜κ³  μžˆλ‹€. 전톡적인 ν•΅μ—°λ£Œ νŒŒνŽΈλ¬Όμ„ ν¬ν•¨ν•œ 3상 거동 해석은 λ‹€λ₯Έ μ€‘λŒ€μ‚¬κ³  해석과 λ§ˆμ°¬κ°€μ§€λ‘œ κ²½ν—˜μ  방법둠에 μ˜μ‘΄ν•œ 보수적인 μ ‘κ·Ό λ°©λ²•μ΄λ‚˜ 고체 μž…μžμ™€ 이상 유체λ₯Ό λͺ¨λ‘ μ—°μ†μ²΄λ‘œ κ°€μ •ν•˜λŠ” λ‹€μœ μ²΄ λͺ¨λΈμ„ 기반으둜 μ΄λ£¨μ–΄μ‘Œλ‹€. μ΅œκ·Όμ—λŠ” μ΄λŸ¬ν•œ 방법둠듀이 κ°€μ§€λŠ” 본질적인 ν•œκ³„λ₯Ό κ·Ήλ³΅ν•˜κ³ μž 고체 μž…μžλ“€ μ‚¬μ΄μ˜ μΆ©λŒμ΄λ‚˜ νšŒμ „μ„ λ³„κ°œλ‘œ λ‹€λ£¨λŠ” μ΄μ‚°μš”μ†Œλ²•(Discrete Element Method, DEM)κ³Ό 격자 기반의 μ˜€μΌλŸ¬λ¦¬μ•ˆ μ „μ‚°μœ μ²΄ν•΄μ„(CFD) 기법을 μ—°κ³„ν•˜λŠ” ν˜•νƒœλ‘œ μ„œλ‘œ κ°„μ˜ μƒν˜Έμž‘μš©μ„ ν•΄μ„ν•˜λŠ” 연ꡬ듀이 많이 이루어지고 μžˆλ‹€. ν•œνŽΈ, μ΅œκ·Όμ—λŠ” ν•˜λ“œμ›¨μ–΄ 및 μ†Œν”„νŠΈμ›¨μ–΄μ˜ μ„±λŠ₯이 λΉ„μ•½μ μœΌλ‘œ μ’‹μ•„μ§€λ©΄μ„œ 격자(Grid)에 κΈ°λ°˜ν•˜μ§€ μ•Šκ³  질점 ν•˜λ‚˜ν•˜λ‚˜μ˜ μ›€μ§μž„μ„ λ”°λΌκ°€λ©΄μ„œ μœ λ™μ— λŒ€ν•œ 지배방정식을 ν•΄μ„ν•˜λŠ” λΌκ·Έλž‘μ§€μ•ˆ 유체 ν•΄μ„κΈ°λ²•μ˜ μ‘μš©μ΄ λŠ˜μ–΄λ‚˜κ³  μžˆλ‹€. λΌκ·Έλž‘μ§€μ•ˆ 해석 κΈ°λ²•μ—μ„œλŠ” λ‹€μƒμœ λ™ 해석 μ‹œμ— 앑체와 기체 두 상을 μ™„μ „νžˆ λ³„κ°œμ˜ μ˜μ—­μœΌλ‘œ ν•΄μƒν•˜μ—¬ 지배방정식을 ν’€κΈ° λ•Œλ¬Έμ— 계면 마찰λ ₯μ΄λ‚˜ ν•­λ ₯, μ–‘λ ₯ 등에 λŒ€ν•œ λ³„λ„μ˜ 상관식 없이 제 1 원리 기반으둜 μœ λ™μ„ 해석할 수 μžˆμ–΄, μ΄μƒμœ λ™μ— λŒ€ν•œ 보닀 근본적인 해석이 κ°€λŠ₯ν•˜λ‹€. ν•΅μ—°λ£Œ 파편물의 생성 및 3상 냉각 거동과 κ΄€λ ¨λœ ν˜„μƒλ“€μ€ λŒ€λΆ€λΆ„ 기체 μƒμ˜ μƒμ„±μ΄λ‚˜ μ΄μƒμœ λ™ μžμ—°λŒ€λ₯˜μ— 영ν–₯을 λ°›λŠ” ν˜„μƒλ“€λ‘œ 앑체 기체 μ‚¬μ΄μ˜ μΈν„°νŽ˜μ΄μŠ€κ°€ λ³΅μž‘ν•˜κ³  역동적인 κ²½ν–₯이 있기 λ•Œλ¬Έμ—, 앑체-기체 λ‹€μƒμœ λ™μ— 효과적인 λΌκ·Έλž‘μ§€μ•ˆ 기반의 μœ μ²΄ν•΄μ„ 기법과 강체 μ΄μ‚°μš”μ†Œλ²•(DEM)을 μ—°κ³„ν•˜λ©΄ 효과적인 3상 μœ λ™ 해석 체계λ₯Ό ꡬ좕할 수 μžˆλ‹€. ν•˜μ§€λ§Œ, ν•΅μ—°λ£Œ νŒŒνŽΈλ¬Όμ„ ν¬ν•¨ν•œ 3상 냉각 거동과 κ΄€λ ¨ν•˜μ—¬ λΌκ·Έλž‘μ§€μ•ˆ μž…μž 기반 μœ μ²΄ν•΄μ„ 기법을 ν™œμš©ν•œ μ—°κ΅¬λŠ” μ„Έκ³„μ μœΌλ‘œλ„ 아직 μˆ˜ν–‰λœ λ°”κ°€ μ—†λ‹€. μ΄λŸ¬ν•œ ν•„μš”μ„±μ— 따라, λ³Έ μ—°κ΅¬μ—μ„œλŠ” λŒ€ν‘œμ μΈ μž…μž 기반의 μœ μ²΄ν•΄μ„ 방법둠 쀑 ν•˜λ‚˜μΈ μ™„ν™”μž…μžμœ μ²΄λ™μ—­ν•™(Smoothed Particle Hydrodynamics, SPH) 기법과 κ°•μ²΄μ˜ 좩돌, 병진, νšŒμ „ μš΄λ™μ„ μ§μ ‘μ μœΌλ‘œ λ‹€λ£¨λŠ” μ΄μ‚°μš”μ†Œλ²•(DEM)의 연계λ₯Ό 톡해 고체 μž…μžλ₯Ό ν¬ν•¨ν•œ 3상 μœ λ™ 해석을 μœ„ν•œ λΌκ·Έλž‘μ§€μ•ˆ 해석 체계λ₯Ό κ΅¬μΆ•ν•˜μ˜€λ‹€. 고체 μž…μžμ™€ 유체 μ‚¬μ΄μ˜ μ—°κ³„λŠ” mm μ΄ν•˜μ˜ μŠ€μΌ€μΌμ„ κ°€μ§€λŠ” ν•΅μ—°λ£Œ 파편물의 ν˜•μƒμ  νŠΉμ„±μ„ κ³ λ €ν•˜μ—¬ 두 상 μ‚¬μ΄μ˜ 겹침을 ν—ˆμš©ν•˜μ—¬ μš΄λ™λŸ‰ κ΅ν™˜μ„ λͺ¨λΈλ§ν•˜λŠ” 비해상(unresolved) λ°©μ‹μœΌλ‘œ μ΄λ£¨μ–΄μ‘Œλ‹€. λ˜ν•œ, SPH 유체 λͺ¨λΈ, DEM 강체 λͺ¨λΈ, SPH-DEM 연계 λͺ¨λΈ 각각에 λŒ€ν•œ 검증을 λ‹€μ–‘ν•œ μŠ€μΌ€μΌμ—μ„œ λ‹€μ–‘ν•œ μ‹€ν—˜ μ—°κ΅¬λ“€κ³Όμ˜ 비ꡐλ₯Ό 톡해 μˆ˜ν–‰ν•˜μ˜€λ‹€. ν•œνŽΈ, μ™„ν™”μž…μžμœ μ²΄λ™μ—­ν•™(SPH) 기법과 μ΄μ‚°μš”μ†Œλ²•(DEM) λͺ¨λΈμ€ λΌκ·Έλž‘μ§€μ•ˆ 해석 κΈ°λ²•μ˜ νŠΉμ„±μƒ 컴퓨터 μ„±λŠ₯의 비약적인 λ°œμ „μ—λ„ λΆˆκ΅¬ν•˜κ³  μ˜€μΌλŸ¬λ¦¬μ•ˆ 해석 기법에 λΉ„ν•΄ 계산 효율 및 μ‹œκ°„μ— λŒ€ν•œ μƒλŒ€μ μΈ μ œμ•½μ΄ μ‘΄μž¬ν•œλ‹€. 특히 앑체-기체의 μ΄μƒμœ λ™ 해석을 λ‹€λ£° 경우 기체 μƒμ˜ 밀도가 μž‘κΈ° λ•Œλ¬Έμ— λΌκ·Έλž‘μ§€μ•ˆ μœ μ²΄ν•΄μ„ κΈ°λ²•μ—μ„œλŠ” 더 μž‘μ€ μ‹œκ°„ 간격이 μš”κ΅¬λœλ‹€. 이에 λ³Έ μ—°κ΅¬μ—μ„œλŠ” κ·Έλž˜ν”½ 처리 μž₯치 (Graphics Processing Unit, GPU)λ₯Ό ν™œμš©ν•˜μ—¬ SPH 해석, DEM 해석, SPH-DEM 연계 해석이 λͺ¨λ‘ 각 λΌκ·Έλž‘μ§€μ•ˆ μ§ˆμ μ— λŒ€ν•΄ λ™μ‹œμ— μˆ˜ν–‰λ  수 μžˆλ„λ‘ GPU 기반의 μ—°κ³„μ½”λ“œ 병렬화 및 가속화λ₯Ό μˆ˜ν–‰ν•˜μ˜€λ‹€. λ§ˆμ§€λ§‰μœΌλ‘œ κ°œλ°œν•œ λΌκ·Έλž‘μ§€μ•ˆ 기반의 3μƒμœ λ™ 해석 μ²΄κ³„μ˜ μœ μš©μ„± μž…μ¦μ„ μœ„ν•΄ μ›μžλ‘œ μ€‘λŒ€μ‚¬κ³  ν›„κΈ° κ³Όμ •μ—μ„œ λ°œμƒν•  수 μžˆλŠ” ν•΅μ—°λ£Œ 파편물 μΈ΅(debris bed)의 평탄화(self-leveling) 거동에 λŒ€ν•œ 검증 해석을 μˆ˜ν–‰ν•˜μ˜€λ‹€. μ‹œκ°„μ— λ”°λ₯Έ 파편물 μΈ΅ ν‘œλ©΄μ˜ ν˜•μƒ λ³€ν™”λ₯Ό 타 κΈ°κ΄€μ—μ„œ μˆ˜ν–‰λœ κΈ°μ²΄μ£Όμž… μ‹€ν—˜κ³Ό λΉ„κ΅ν•˜λŠ” ν˜•νƒœλ‘œ 검증이 μ΄λ£¨μ–΄μ‘Œλ‹€. 뢄석 κ²°κ³Ό, λ³Έ μ—°κ΅¬μ—μ„œ κ°œλ°œν•œ SPH-DEM 연계해석 μ½”λ“œκ°€ 고체 μž…μž 상을 ν¬ν•¨ν•œ 수λ ₯학적 3상 거동을 μ •λŸ‰μ μœΌλ‘œ, μ •μ„±μ μœΌλ‘œ 잘 ν•΄μ„ν•˜λŠ” 것을 ν™•μΈν•˜μ˜€λ‹€. λ³Έ μ—°κ΅¬μ—μ„œ κ°œλ°œν•œ λΌκ·Έλž‘μ§€μ•ˆ 기반의 SPH-DEM 3μƒμœ λ™ 해석 μ²΄κ³„λŠ” μ›μžλ‘œ μ€‘λŒ€μ‚¬κ³ μ˜ 해석적 κ΄€μ μ—μ„œ 기쑴의 μˆ˜μΉ˜ν•΄μ„ 기법듀이 닀루기 μ–΄λ €μ› λ˜ ν˜„μƒλ“€μ— λŒ€ν•œ λŒ€μ•ˆ λ˜λŠ” μƒν˜Έ λ³΄μ™„μ˜ 역할을 ν•  수 μžˆλ‹€. λ˜ν•œ, λ³Έ μ—°κ΅¬μ—μ„œ κ°œλ°œν•œ μ½”λ“œλŠ” 제 1원리 기반의 물리 법칙을 기반으둜 μœ λ™ 및 κ°•μ²΄μ˜ 거동을 ν•΄μ„ν•˜κΈ° λ•Œλ¬Έμ— μ‹€ν—˜μœΌλ‘œ κ΅¬ν˜„ν•˜κΈ° μ–΄λ €μš΄ μ‘°κ±΄μ΄λ‚˜ μŠ€μΌ€μΌμ— λŒ€ν•œ 수치적 μž¬ν˜„μ΄ κ°€λŠ₯ν•˜κ³ , 이λ₯Ό λ°”νƒ•μœΌλ‘œ 기쑴의 μŠ€μΌ€μΌλ§ 법칙을 κ²€μ¦ν•˜κ±°λ‚˜ μ‹€ν—˜ κ²°κ³Όκ°€ μ—†λŠ” μ˜μ—­μ—μ„œ 수치 데이터λ₯Ό μƒμ„±ν•˜μ—¬ 기쑴의 상관식을 κ°œμ„ ν•˜λŠ”λ° ν™œμš©ν•  수 μžˆλ‹€. μ΄λŸ¬ν•œ μ μ—μ„œ λ³Έ μ—°κ΅¬λŠ” μ›μžλ‘œ μ€‘λŒ€μ‚¬κ³ μ˜ ν•΄μ„μ΄λ‚˜ μ•ˆμ „μ„± 평가와 κ΄€λ ¨ν•˜μ—¬ λΆˆν™•μ‹€μ„±μ„ μ €κ°ν•˜λŠ”λ° κΈ°μ—¬ν•œλ‹€.Chapter 1 Introduction 1 1.1 Background and Motivation 1 1.2 Previous Studies 3 1.2.1 Numerical Studies on Particulate Debris Bed 3 1.2.2 SPH-DEM Phase Coupling 4 1.3 Objectives and Scope 5 Chapter 2 Fluid Phase: Smoothed Particle Hydrodynamics 8 2.1 Smoothed Particle Hydrodynamics (SPH) 9 2.1.1 SPH Particle Approximation 9 2.1.2 SPH Particle Approximation of Derivatives 10 2.1.3 Kernel Function 11 2.1.4 Accuracy of SPH Approximation 12 2.1.5 Governing Equations for Incompressible Flow 14 2.2 Multi-phase SPH Governing Equations 16 2.2.1 Treatment of Multi-Phase Flow 16 2.2.2 Normalized Density Model 18 2.2.3 Continuum Surface Force (CSF) Model 19 2.3 Implementation of SPH Model 21 2.3.1 Algorithm of SPH Code 21 2.3.2 Nearest Neighboring Particles Searching (NNPS) 22 2.3.3 Time Integration 23 2.4 V&V Simulations 24 2.4.1 Rayleigh-Taylor Instability 25 2.4.2 Bubble Terminal Velocity 25 2.4.3 Dam-Break Simulation 25 2.4.4 Centralized Sloshing Simulation 26 Chapter 3 Solid Phase: Discrete Element Method 45 3.1 Discrete Element Method (DEM) 46 3.2 DEM Contact Force 47 3.2.1 Soft-sphere Contact Model 47 3.2.2 Contact Force Model 48 3.2.3 Hertz-Mindlin Contact Force Model 49 3.3 Wall Boundary Conditions 52 3.3.1 Versatile Wall Boundary Model 52 3.3.2 Particle Collision with the Wall 54 3.3.3 Sliding and Rolling on the Wall Boundary 56 3.4 DEM Implementation Algorithm 57 3.4.1 Contact Detection 58 3.4.2 Estimation of Relative Velocity 59 3.4.3 Calculation of Contact Force 60 3.4.4 Wall Boundary Conditions and Time Integration 60 3.5 V&V and Simulations 61 3.5.1 Conservation of Momentum and Angular Momentum 62 3.5.2 Conservation of Energy in Elastic Collision 63 3.5.3 V&V Simulations for Wall Boundary Model 63 3.5.4 Granular Collapse of Spherical Particles 64 Chapter 4 Two-way Phase Coupling of SPH and DEM 76 4.1 Unresolved Coupling of SPH and DEM 76 4.2 Governing Equations 78 4.2.1 SPH Particles: Locally Averaged N-S Equations 78 4.2.2 DEM Particles: Coupling Forces Acting on Solid Particles 80 4.2.3 SPH Particles: Reaction Force from Momentum Exchange 82 4.3 Algorithm of SPH-DEM Coupled Model 83 4.4 V&V Simulations for SPH-DEM Coupled Model 84 4.4.1 Single DEM Particle Behavior 85 4.4.2 Pressure Drop through Packed Bed 87 4.4.3 Granular Flow in Liquid: 3D Dam-Break 89 Chapter 5 GPU Parallelization of Coupled SPH-DEM Code 103 5.1 Parallelization of Governing Equations 104 5.1.1 GPU-based Parallelization 104 5.1.2 Parallelization of SPH-DEM Governing Equations 104 5.2 Parallelization of NNPS and Contact Detection 105 5.3 Results of GPU Parallelization 107 5.3.1 Speedup in Computation Time 107 5.3.2 Parallelization Factors 107 Chapter 6 Code Application to Vapor-Driven Leveling Behavior of Spherical Debris 113 6.1 Self-Leveling Behavior of Debris Bed 114 6.1.1 Self-Leveling Issue in LWR 114 6.1.2 Self-Leveling Behavior in Terms of Debris Coolability 114 6.2 Benchmark Experiment 116 6.3 SPH-DEM Simulation Setup 118 6.3.1 Properties and Simulation Conditions 118 6.3.2 Sequence of SPH-DEM Leveling Simulation 120 6.3.3 Determination of Inclined Angle 121 6.4 Validation Results and Discussions 121 6.4.1 SPH-DEM Simulation Results 121 6.4.2 Validation Result 122 6.4.3 Effect of Gas Flow Rate 122 Chapter 7 Summary 129 7.1 Summary 129 7.2 Recommendations 131 References 134 κ΅­λ¬Έ 초둝 142Docto

    Smoothed Particle Hydrodynamics for Navier-Stokes Fluid Flow Application

    Get PDF
    The aim of this publication is to introduce the particle based computational fluid dynamics (CFD) method smoothed particle hydrodynamics (SPH) and introduce an applicable and valid SPH implementation for practical cases. For this purpose, current research approaches are combined regarding performance and numerical stability.  The principles of the method, the mathematical basics and the discretization of the Navier-Stokes equations are clarified. Furthermore, the implementation of method-specific boundary conditions, wall, inlet and outlet, as well as several correction procedures and a surface tension setup into the present code framework are described. The advantages and validity of the method are shown based on different cases. The free surface fluid behavior of a dam break is compared to experimental data of the time dependent water level of selected positions. A Karman vortex street is validated by its Strouhal number for different Reynolds numbers. The frequency of an oscillating drop is analysed and compared to the analytical solution.  The SPH is utilized for pipe flows influenced by a backward facing step and shows an expected qualitative flow field

    Pairwise Force SPH Model for Real-Time Multi-Interaction Applications

    Get PDF
    In this paper, we present a novel pairwise-force smoothed particle hydrodynamics (PF-SPH) model to allow modeling of various interactions at interfaces in real time. Realistic capture of interactions at interfaces is a challenging problem for SPH-based simulations, especially for scenarios involving multiple interactions at different interfaces. Our PF-SPH model can readily handle multiple kinds of interactions simultaneously in a single simulation; its basis is to use a larger support radius than that used in standard SPH. We adopt a novel anisotropic filtering term to further improve the performance of interaction forces. The proposed model is stable; furthermore, it avoids the particle clustering problem which commonly occurs at the free surface. We show how our model can be used to capture various interactions. We also consider the close connection between droplets and bubbles, and show how to animate bubbles rising in liquid as well as bubbles in air. Our method is versatile, physically plausible and easy-to-implement. Examples are provided to demonstrate the capabilities and effectiveness of our approach

    Enhanced SPH modeling of free-surface flows with large deformations

    Get PDF
    The subject of the present thesis is the development of a numerical solver to study the violent interaction of marine flows with rigid structures. Among the many numerical models available, the Smoothed Particle Hydrodynamics (SPH) has been chosen as it proved appropriate in dealing with violent free-surface flows. Due to its Lagrangian and meshless character it can naturally handle breaking waves and fragmentation that generally are not easily treated by standard methods. On the other hand, some consolidated features of mesh-based methods, such as the solid boundary treatment, still remain unsolved issues in the SPH context. In the present work a great part of the research activity has been devoted to tackle some of the bottlenecks of the method. Firstly, an enhanced SPH model, called delta-SPH, has been proposed. In this model, a proper numerical diffusive term has been added in the continuity equation in order to remove the spurious numerical noise in the pressure field which typically affects the weakly-compressible SPH models. Then, particular attention has been paid to the development of suitable techniques for the enforcement of the boundary conditions. As for the free-surface, a specific algorithm has been designed to detect free-surface particles and to define a related level-set function with two main targets: to allow the imposition of peculiar conditions on the free-surface and to analyse and visualize more easily the simulation outcome (especially in 3D cases). Concerning the solid boundary treatment, much effort has been spent to devise new techniques for handling generic body geometries with an adequate accuracy in both 2D and 3D problems. Two different techniques have been described: in the first one the standard ghost fluid method has been extended in order to treat complex solid geometries. Both free-slip and no-slip boundary conditions have been implemented, the latter being a quite complex matter in the SPH context. The proposed boundary treatment proved to be robust and accurate in evaluating local and global loads, though it is not easy to extend to generic 3D surfaces. The second technique has been adopted for these cases. Such a technique has been developed in the context of Riemann-SPH methods and in the present work is reformulated in the context of the standard SPH scheme. The method proved to be robust in treating complex 3D solid surfaces though less accurate than the former. Finally, an algorithm to correctly initialize the SPH simulation in the case of generic geometries has been described. It forces a resettlement of the fluid particles to achieve a regular and uniform spacing even in complex configurations. This pre-processing procedure avoids the generation of spurious currents due to local defects in the particle distribution at the beginning of the simulation. The delta-SPH model has been validated against several problems concerning fluid-structure interactions. Firstly, the capability of the solver in dealing with water impacts has been tested by simulating a jet impinging on a flat plate and a dam-break flow against a vertical wall. In this cases, the accuracy in the prediction of local loads and of the pressure field have been the main focus. Then, the viscous flow around a cylinder, in both steady and unsteady conditions, has been simulated comparing the results with reference solutions. Finally, the generation and propagation of 2D gravity waves has been simulated. Several regimes of propagation have been tested and the results compared against a potential flow solver. The developed numerical solver has been applied to several cases of free-surface flows striking rigid structures and to the problem of the generation and evolution of ship generated waves. In the former case, the robustness of the solver has been challenged by simulating 2D and 3D water impacts against complex solid surfaces. The numerical outcome have been compared with analytical solutions, experimental data and other numerical results and the limits of the model have been discussed. As for the ship generated waves, the problem has been firstly studied within the 2D+t approximation, focusing on the occurrence and features of the breaking bow waves. Then, a dedicated 3D SPH parallel solver has been developed to tackle the simulation of the entire ship in constant forward motion. This simulation is quite demanding in terms of complexities of the boundary geometry and computational resources required. The wave pattern obtained has been compared against experimental data and results from other numerical methods, showing in both the cases a fair and promising agreement

    Toward the development of a virtual spray test-rig using the Smoothed Particle Hydrodynamics method

    Get PDF
    In this work we present the numerical simulation of air-assisted liquid atomization at high pressure us- ing the Smoothed Particle Hydrodynamics (SPH) method. Different post-processing tools are applied to facilitate the comparison with experimental observations. This allows to quantitatively validate the nu- merical method against the experiment, in terms of (i) frequency of the Kelvin–Helmholtz instability that develops on the jet surface, and (ii) statistical distribution of the jet intact length. The qualitative com- parison also shows a good prediction of the jet global instability and of the fragmented liquid lumps, with regards to length and time scales. In addition, the post-processing tools also give access to the local parameters of the generated spray in the vicinity of the nozzle, which are not easily accessible in a real experiments. Using these tools, 1D profiles and 2D maps of the liquid phase properties such as the vol- ume fraction, the droplet concentration, the Sauter Mean Diameter (SMD) and the droplet sphericity are presented. Because of the Lagrangian nature of the SPH method, it is also possible to monitor the whole atomization cascade as a causal tree, from the primary instabilities to the spray characteristics. This tree contains various information such as the fragmentation spectrum and the breakup activity, which are of great interest for researchers and engineers. Hence, the capability of the Smoothed Particle Hydrody- namics (SPH) method for simulating air-assisted atomization at high ambient pressure is demonstrated as well as its applicability to realistic configurations. This is a first step towards the development of a complete virtual spray test-rig

    Development of GPU-based SPH Framework for Hydrodynamic Interactions With Non-spherical Solid Debris

    Get PDF
    일본의 ν›„μΏ μ‹œλ§ˆ 사고 이후 μ›μžλ‘œ μ€‘λŒ€ 사고에 λŒ€ν•œ μ—°κ΅¬μ˜ ν•„μš”μ„±κ³Ό λŒ€μ²˜ λŠ₯λ ₯ 확보에 λŒ€ν•œ μ€‘μš”μ„±μ΄ 점점 μ¦κ°€ν•˜κ³  μžˆλ‹€. 사고 μ‹œ λ°œμƒν•  수 μžˆλŠ” 노심 용육물 거동에 λŒ€ν•œ ν‰κ°€λŠ” 용육물-콘크리트 μƒν˜Έμž‘μš©(MCCI, Molten Core Concrete Interaction)κ³Ό 증기 ν­λ°œλ‘œλΆ€ν„°μ˜ μ›μžλ‘œ 노심 냉각성 및 건전성에 λ”°λ₯Έ μž¬μž„κ³„ μΈ‘λ©΄μ—μ„œ 맀우 μ€‘μš”ν•˜λ‹€. 특히 OPR 1000의 경우, 사전 좩수 쑰건(Wet cavity condition)을 기본적인 μ›μžλ‘œ μ™Έλ²½ 냉각 λŒ€μ‘ μ „λž΅μœΌλ‘œ μ±„νƒν•¨μœΌλ‘œμ¨ ν•΅μ—°λ£Œ-λƒ‰κ°μž¬ μƒν˜Έμž‘μš©(FCI, Fuel Coolant Interaction) λ°˜μ‘μ΄ ν•„μ—°μ μœΌλ‘œ λ°œμƒν•˜λŠ” κ²ƒμœΌλ‘œ μ•Œλ €μ Έ μžˆλ‹€. [Jin, 2014] FCI ν˜„μƒμ€ μž„μ˜ ν˜•νƒœμ˜ ν•΅μ—°λ£Œ 고체 파편물과 λƒ‰κ°μž¬μ˜ μƒν˜Έμž‘μš©λΏλ§Œ μ•„λ‹ˆλΌ, λƒ‰κ°μž¬ λΉ„λ“± ν˜„μƒ 등도 ν¬ν•¨ν•˜λŠ” λ‹€μœ μ²΄, 닀상 ν˜„μƒμœΌλ‘œ κ·Έ ν˜„μƒμ΄ 맀우 λ³΅μž‘ν•˜λ‹€. 이 κ³Όμ •μ—μ„œ μ›μžλ‘œ 건물 ν•˜λΆ€μ— 고체 파편물이 ν‡΄μ λ˜μ–΄ μž”ν•΄ 측이 ν˜•μ„±λ˜κ³ , κ·Έ 냉각성에 따라 μ‚¬κ³ μ˜ λ‹€μŒ 진행 상황에 영ν–₯을 쀄 수 μžˆλ‹€. μ΄λŸ¬ν•œ λΉ„κ΅¬ν˜• 고체 파편물 거동에 λŒ€ν•œ 이해λ₯Ό μœ„ν•΄ 강체 κ°œλ…μ„ μ μš©ν•œ 고체 해석 μ²΄κ³„λŠ” 쒋은 접근법이 될 수 μžˆλ‹€. λ”°λΌμ„œ λ³Έ μ—°κ΅¬μ—μ„œλŠ” μœ μ²΄μ™€ 고체 κ°„ 수λ ₯학적 μƒν˜Έμž‘μš© 해석을 μœ„ν•΄ μž…μžμœ μ²΄λ™μ—­ν•™(SPH, Smoothed Particle Hydrodynamics) 기법과 강체역학(RBD, Rigid Body Dynamics) 기법을 μ—°κ³„ν•˜μ—¬ λΌκ·Έλž‘μ§€μ•ˆ 해석 체계λ₯Ό κ΅¬μΆ•ν•˜μ˜€λ‹€. μ™„ν™”μž…μžμœ μ²΄λ™μ—­ν•™ 기법은 해석 유체λ₯Ό μœ ν•œκ°œμ˜ μž…μžλ‘œ ν‘œν˜„ν•¨μœΌλ‘œμ¨ μœ λ™μ„ ν•΄μ„ν•˜λŠ” λΌκ·Έλž‘μ§€μ•ˆ 해석 기법 쀑 ν•˜λ‚˜μ΄λ‹€. κ°œλ³„ μž…μžλ“€μ˜ μ›€μ§μž„μœΌλ‘œ μœ λ™μ„ ν•΄μ„ν•˜λ―€λ‘œ λΉ„μ„ ν˜•μ˜ λŒ€λ₯˜ν•­μ— λŒ€ν•œ 계산이 ν•„μš” μ—†μœΌλ©°, μž…μžκ°€ μΆ”κ°€λ˜κ±°λ‚˜ 사라지지 μ•ŠλŠ” ν•œ 해석 κ³„μ˜ 전체 μ§ˆλŸ‰μ€ μžλ™μœΌλ‘œ λ³΄μ‘΄λœλ‹€. μ΄λŸ¬ν•œ λΌκ·Έλž‘μ§€μ•ˆ κΈ°λ²•μ˜ νŠΉμ„±μœΌλ‘œ SPH 방법은 자유 ν‘œλ©΄ μœ λ™, λ‹€μœ μ²΄ μœ λ™, 닀상 μœ λ™, ν˜•νƒœ λ³€ν™”κ°€ 큰 μœ λ™ 등에 λŒ€ν•΄ 해석 μž₯점을 κ°–λŠ”λ‹€. λ³Έ μ—°κ΅¬μ—μ„œλŠ” SPH 기법을 μ μš©ν•œ in-house SOPHIA μ½”λ“œλ₯Ό μ‚¬μš©ν•˜μ—¬ λΉ„μ••μΆ• 닀상 μœ λ™ 해석을 μˆ˜ν–‰ν•˜μ˜€μœΌλ©°, 벀치마크 λ°μ΄ν„°μ™€μ˜ λΉ„κ΅μ—μ„œ 쒋은 검증 해석 κ²°κ³Όλ₯Ό λ³΄μ˜€λ‹€. 강체역학은 μ™Έλ ₯에 μ˜ν•΄ ν˜•νƒœκ°€ λ³€ν•˜μ§€ μ•ŠλŠ” κ°•μ²΄μ˜ κ°œλ…μ„ μ΄μš©ν•˜μ—¬ 고체의 병진 μš΄λ™κ³Ό νšŒμ „ μš΄λ™μ„ ν•΄μ„ν•˜λŠ” 연ꡬ 뢄야이닀. λ³Έ μ—°κ΅¬μ—μ„œλŠ” μ΄μ‚°μš”μ†Œλ²•(DEM, Discrete Element Method) λΆ„μ•Όμ—μ„œ 였랜 μ‹œκ°„ λ™μ•ˆ 널리 μ‚¬μš©λ˜κ³  κ²€μ¦λ˜μ—ˆλ˜ Hertz-Mindlin 좩돌 λͺ¨λΈμ„ μ μš©ν•˜μ—¬ 강체 κ°„ 좩돌 해석을 κ΅¬ν˜„ν•˜μ˜€λ‹€. κ°•μ²΄λŠ” μœ ν•œκ°œμ˜ μž…μžλ“€λ‘œ ν‘œν˜„ν•  수 있으며, 강체 κ°„ μΆ©λŒμ€ 각 강체λ₯Ό κ΅¬μ„±ν•˜κ³  μžˆλŠ” μž…μžμŒμ˜ μž‘μ€ 쀑첩을 기반으둜 κ³„μ‚°λœλ‹€. λ³Έ μ—°κ΅¬μ—μ„œλŠ” μž…μžκΈ°λ°˜μ˜ 강체역학 해석 μ½”λ“œλ₯Ό μ΄μš©ν•˜μ—¬ 단일 강체 및 닀쀑 강체 μΆ©λŒμ— λŒ€ν•΄ 검증 해석을 μˆ˜ν–‰ν•˜μ˜€μœΌλ©°, 해석해 및 벀치마크 데이터 결과와 잘 μΌμΉ˜ν•¨μ„ ν™•μΈν•˜μ˜€λ‹€. μ›μžλ ₯ λΆ„μ•Όμ—μ„œ λ°œμƒν•  수 μžˆλŠ” λΉ„κ΅¬ν˜• 고체와 μœ μ²΄κ°„ μƒν˜Έμž‘μš© 해석을 μœ„ν•΄ μ•žμ„œ μ„€λͺ…ν•œ SPH 기법과 강체역학 연계 해석 μ½”λ“œλ₯Ό κ°œλ°œν•˜μ˜€λ‹€. λ³Έ μ—°κ΅¬μ—μ„œ μ μš©ν•œ μ™„μ „ 해상 방식(Fully resolved approach)은 유체-고체의 상이 λΆ„λ¦¬λ˜μ–΄ 있고, 제 1 원리λ₯Ό λ§Œμ‘±ν•˜λ―€λ‘œ 고체의 ν˜•μƒμ— λ”°λ₯Έ 상관식과 ν‘œλ©΄ 적뢄이 ν•„μš”ν•˜μ§€ μ•Šλ‹€λŠ” μž₯점이 μžˆλ‹€. λ˜ν•œ 고체 κ²½κ³„λ©΄μ—μ„œμ˜ μ •ν™•ν•œ μ••λ ₯ 계산을 μœ„ν•΄ 유체 μž…μž 정보λ₯Ό 기반으둜 λ…Έμ΄λ§Œ μ••λ ₯ 경계 쑰건을 μ μš©ν•˜μ˜€λ‹€. λ³Έ μ—°κ΅¬μ—μ„œλŠ” μ΄λŸ¬ν•œ 해상 λ°©μ‹μ˜ 유체-강체 연계 해석 μ½”λ“œλ₯Ό μ΄μš©ν•˜μ—¬ λΉ„κ΅¬ν˜• 고체와 유체의 수λ ₯학적 μƒν˜Έμž‘μš©μ— λŒ€ν•œ 검증 해석을 μˆ˜ν–‰ν•˜μ˜€μœΌλ©°, μ„ ν–‰λœ μ‹€ν—˜κ³Όμ˜ λΉ„κ΅μ—μ„œ 쒋은 κ²°κ³Όλ₯Ό λ³΄μ˜€λ‹€. μœ λ™ 해석을 μœ„ν•΄ λ³Έ 연ꡬ에 μ μš©ν•œ SPH λ°©λ²•μ—μ„œλŠ” μˆ˜μ‹λ“€μ΄ 맀우 μ„ ν˜•μ μ΄κ³  외연적(Explicit)으둜 계산을 μˆ˜ν–‰ν•˜κΈ° λ•Œλ¬Έμ— 각 μž…μžμ— λŒ€ν•œ 계산이 λ…λ¦½μ μœΌλ‘œ μˆ˜ν–‰λ˜μ–΄λ„ λ¬Έμ œκ°€ μ—†λ‹€. λ”°λΌμ„œ SPH 방법은 계산 병렬화에 μ΅œμ ν™”λœ λ°©λ²•μœΌλ‘œ 잘 μ•Œλ €μ Έ 있으며, λŒ€κ·œλͺ¨ 고해상도 해석을 μœ„ν•΄ μ΄λŠ” ν•„μˆ˜μ μ΄λ‹€. λ˜ν•œ μž…μž 기반의 강체 계산을 μœ„ν•΄μ„œλŠ” 효율적인 계산 μ•Œκ³ λ¦¬μ¦˜μ΄ ν•„μš”ν•˜λ‹€. λ”°λΌμ„œ λ³Έ μ—°κ΅¬μ—μ„œλŠ” λŒ€κ·œλͺ¨ 계산과 높은 μ—°μ‚° νš¨μœ¨μ„±μ„ μœ„ν•΄ κ·Έλž˜ν”½μ²˜λ¦¬μž₯치(GPU, Graphic Processing Unit)λ₯Ό μ΄μš©ν•˜μ—¬ 계산 병렬화λ₯Ό μˆ˜ν–‰ν•˜μ˜€μœΌλ©°, 이λ₯Ό μ΄μš©ν•œ 닀쀑 고체와 유체의 μƒν˜Έμž‘μš© ν•΄μ„μ—μ„œ 쒋은 계산 μ„±λŠ₯을 ν™•μΈν•˜μ˜€λ‹€. λ³Έ μ—°κ΅¬μ—μ„œ μˆ˜ν–‰ν•œ λΉ„κ΅¬ν˜• 고체와 유체의 수λ ₯학적 μƒν˜Έμž‘μš©μ„ μœ„ν•œ GPU 기반의 SPH 해석 μ½”λ“œ 개발 연ꡬλ₯Ό 톡해 μ›μžλ‘œ μ€‘λŒ€μ‚¬κ³  μ‹œ λ°œμƒν•  수 μžˆλŠ” λƒ‰κ°μž¬μ™€ ν•΅μ—°λ£Œ 고체 파편물의 수λ ₯학적 μƒν˜Έμž‘μš© 뿐만 μ•„λ‹ˆλΌ, 고체 파편물 κ°„ 역학적 μƒν˜Έμž‘μš©μ— λŒ€ν•΄ 효율적인 해석 체계λ₯Ό κ°œλ°œν•˜μ˜€λ‹€. 이λ₯Ό 톡해 μŠ΅μ‹ 곡동(wet cavity)μ—μ„œ λ°œμƒν•˜λŠ” ν•΅μ—°λ£Œ 고체 파편물의 퇴적 μž‘μš©, μ“°λ‚˜λ―Έ μ‚¬κ³ λ‘œ μΈν•œ ν•΄μ•ˆ ꡬ쑰물의 수λ ₯학적 μƒν˜Έμž‘μš©, 그리고 침수 사고 μ‹œ μ›μžλ‘œ 건물 λ‚΄ λΆ€μœ λ¬Όμ˜ 거동 λ“± μ›μžλ ₯ λΆ„μ•Όμ—μ„œ λ°œμƒν•  수 μžˆλŠ” λ‹€μ–‘ν•œ 고체-유체의 수λ ₯학적 μƒν˜Έμž‘μš©μ— λŒ€ν•œ 해석적 연ꡬ에 μ μš©ν•˜κ³  κΈ°μ—¬ν•  수 μžˆμ„ κ²ƒμœΌλ‘œ κΈ°λŒ€ν•œλ‹€.Since the Fukushima accident, the necessity for researches on severe accidents and the importance of securing the ability to cope with the accidents have been increasing. The evaluation of the molten core behavior that may occur during the accident is very important in terms of re-criticality according to the coolability and integrity of the reactor core from the MCCI (Molten Core Concrete Interaction) and steam explosion. In the case of OPR 1000, especially, FCI (Fuel Coolant Interaction) is known to occur unconditionally because the wet cavity condition has been adopted as a basic strategy for ex-vessel cooling. [Jin, 2014] FCI is a highly complicated phenomenon, which includes multi-fluid, multi-phase interaction between the arbitrary shape of solid debris and coolant as well as coolant boiling. In this process, the debris bed is formed at the bottom of the containment, and its coolability influences the next phase of the accident. For the understanding on the solid debris behavior, a solid system with a rigid body can be a good approach for the non-spherical solid debris analysis. Therefore, in this study, Smoothed Particle Hydrodynamics (SPH) method and Rigid Body Dynamics (RBD) are coupled in a fully Lagrangian manner for the hydrodynamic interactions between fluid and solid. Smoothed Particle Hydrodynamics (SPH) is one of the Lagrangian-based analysis methods which represents the fluid flow as a finite number of particles. Since the flow is analyzed by the motion of individual particles, there is no need to calculate the nonlinear convective term, and the total mass of the system is automatically conserved as long as particles are not added or removed. Through these Lagrangian nature, it is well known that the SPH method is effective for the free surface flow, multi-fluid and multi-phase flow, and highly deformable flow. In this study, the incompressible multi-phase flow analysis has been performed using the in-house SPH code, SOPHIA code, and V&V simulation results showed good agreement with the benchmark data. Rigid Body Dynamics (RBD) is a research field that analyses the translation and rotation of a solid body by using the concept that a rigid body doesn’t change its shape by external forces. In this study, the collision calculation between rigid bodies is implemented by applying the Hertz-Mindlin contact force model commonly used and verified for a long time in the Discrete Element Method (DEM) field. A rigid body can be expressed as a group of finite particles, and the contact forces between solid bodies are calculated based on the small overlap of the particle pairs. Using the particle-based RBD analysis code implemented in this study, V&V simulations on single- and multi- rigid body collisions have been performed and showed good agreement with the analytical solution and the benchmark data. To analyze the hydrodynamic interactions between non-spherical solids and fluids that can occur in the nuclear field, the integrated code has been developed by coupling RBD with SPH code. Since a fully resolved approach adopted in this study as a phase coupling method satisfies the 1st principle and the fluid-solid phase is entirely separated from each other, there is no need for the surface integral and empirical correlations depending on the solid geometry. In addition, the Neumann pressure boundary condition is implemented for accurate pressure estimation at the solid interface using the fluid particle properties. By applying the resolved SPH-RBD coupled code, V&V simulations were carried out on the hydrodynamic interactions of non-spherical solid-fluid and showed good agreement with the experimental data. In the SPH method, since the numerical expression are highly linear and the calculations are performed explicitly, there is no problem even if the calculations for each particle are performed independently. Therefore, the SPH is well known as an optimized method for parallelization, and it is essential for large scale high-resolution simulations. In addition, an efficient computational algorithm is required for particle-based rigid body calculation. In this study, therefore, the parallelization was performed using a Graphical Processing Unit (GPU) for large-scale calculations and high computational efficiency, and it showed a good performance in analyzing the interactions of a large number of solids and fluids particles. Through the researches on the development of a GPU-based SPH framework for the hydrodynamic interaction of non-spherical solids and fluids in this study, an efficient analysis system has been developed for not only the hydrodynamic interaction of solid corium debris with coolant but also the mechanical interaction between solid debris which can occur at the severe accidents in the nuclear reactor. By using this, it is expected that the integrated code will contribute to analytical researches on various accident scenarios that may occur in the nuclear field such as solid fuel debris sedimentation in the wet cavity, hydrodynamic interactions with coastal structures caused by the Tsunami, and the behavior of floating objects in the reactor building at the flooding accident, etc.Chapter 1 Introduction 1 1.1 Background and Motivation 1 1.2 Previous Studies 3 1.2.1 Numerical Studies on FCI Premixing Jet Breakup 3 1.2.2 Numerical Studies on Fluid-Solid Coupling with RBD 4 1.3 Objectives and Scope 5 Chapter 2 Smoothed Particle Hydrodynamics (SPH) 9 2.1 SPH Overview 9 2.1.1 Basic Concept of SPH 9 2.1.2 SPH Particle Approximation 10 2.1.3 SPH Kernel Function 12 2.1.4 SPH Governing Equations 13 2.2 SPH Multi-phase Models 16 2.2.1 Normalized Density Approach 16 2.2.2 Treatments for Multi-phase Flow 17 2.2.3 Surface Tension Force Model 18 2.3 SPH Code Implementation 20 2.3.1 Nearest Neighbor Particle Search (NNPS) 20 2.3.2 Algorithm of SPH Code 21 2.3.3 Time Integration 21 2.3.4 GPU Parallelization 22 Chapter 3 Rigid Body Dynamics (RBD) 30 3.1 RBD Overview 30 3.2 Collision Models of Rigid Body 31 3.2.1 Monaghan Boundary Force (MBF) Model 31 3.2.2 Ideal Plastic Collision Model 33 3.2.3 Impulse-based Boundary Force (IBF) Model 35 3.2.4 Penalty-based Contact Model 37 3.2.5 Determination of Collision Model 40 3.3 Algorithm of RBD 41 3.3.1 Calculation of Rigid Body Information 41 3.3.2 Contact Detection 42 3.3.3 Contact Normal Calculation 42 3.3.4 Contact Force Calculation 45 3.3.5 Summation of Rigid Body Particles 46 3.3.6 Time Integration 47 3.4 GPU Parallelization 48 3.4.1 Algorithm 1: Atomic Operation 49 3.4.2 Algorithm 2: Sorting 50 3.5 Code V&V Simulations 51 3.5.1 Conservation of Momentum & Angular Momentum 51 3.5.2 Conservation of Kinetic Energy in Elastic Collision 52 3.5.3 Bouncing Block 53 3.5.4 Sliding Block on a Slope 55 3.5.5 Collapse of Stacked Multi-body 57 Chapter 4 Two-way Coupling of SPH-RBD 75 4.1 Resolved Approach 75 4.2 Governing Equations 75 4.2.1 Solid Phase 75 4.2.2 Fluid Phase 78 4.3 Algorithm of SPH-RBD Code 78 4.4 Code V&V Simulations 81 4.4.1 Karman Vortex Problem 81 4.4.2 Water Entry 84 4.4.3 Sinking & Rotating Body 85 4.4.4 Floating & Falling Body 85 4.4.5 Collapse of Stacked Multi-body with Fluid 87 4.4.6 Code Application to Non-spherical Debris Sedimentation 89 Chapter 5 Conclusion 110 5.1 Summary 110 5.2 Recommendations 112 Nomenclature 114 Bibliography 117 κ΅­λ¬Έ 초둝 127λ°•

    An improved SPH scheme for cosmological simulations

    Get PDF
    We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and largely improves the SPH performance in capturing the development of gas-dynamical instabilities. We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations. Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder, thinner and more extended and our results on galaxy clusters show entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard SPH, thus becoming the core of an efficient code for large cosmological simulations.Comment: 21 figures, 2 tables, accepted to MNRA
    • …
    corecore