• 한국강구조학회 논문집
• /
• 제23권5호
• /
• pp.535-546
• /
• 2011

본 연구에서는 다양한 명시적 호장법을 사용하여 공간프레임의 반강접 탄소성 후좌굴 해석을 수행하였다. 이를 위해 이전 연구를 발전시켜 다양한 명시적 알고리즘의 호장법과 명시적, 묵시적 해석법에 동시에 적용 가능한 반강접 탄소성 공간프레임요소를 제안하였다. 다양한 명시적 호장법은 예측단계와 수렴단계에 명시적 해석법인 동적이완법을 적용한 것을 의미한다. 따라서 명시적 호장법에는 명시적(예측단계)-명시적(예측단계) 호장법, 명시적(예측단계)-묵시적(수렴단계) 호장법, 묵시적(예측단계)-명시적(수렴단계) 호장법으로 구분된다. 또한 명시적 호장법에 적용 가능하도록 수정된 반강접 탄소성 공간프레임요소는 오일러리안 유한변형이론에 의해 강체회전변형을 고려하였기 때문에 대변위가 발생하는 기하학적 비선형 문제에 적용될 수 있고, 완전 탄소성 소성힌지 알고리즘에 의한 재료적 비선형성을 고려하였으며, 부재내부에 정적 응축된 회전 및 축방향 성분의 선형 스프링에 의해 접합부 반강접 특성을 반영하였다. 제안된 해석법을 이용하여 검증예제를 수행함으로써 본 연구에서 제안된 다양한 명시적 호장법 및 공간프레임요소의 정확성을 검증한다.

• Wind and Structures
• /
• 제23권4호
• /
• pp.301-311
• /
• 2016

The numerous benefits of Savonius turbine such as simple in structure, has appropriate self-start ability, relatively low operating velocity, water acceptance from any direction and low environmental impact have generated interests among researchers. However, it suffers from a lower efficiency compared to other types of water turbine. To improve its performance, parameters such flow pattern, pressure and velocity in different conditions must be analyzed. For this purpose, a detailed description on the flow field of various types of Savonius rotors is required. This article presents a numerical study on a nonlinear two-dimensional flow over a classic Savonius type rotor and a Benesh type rotor. In this experiment, sliding mesh was used for solving the motion of the bucket. The unsteady Reynolds averaged Navier-Stokes equations were solved for velocity and pressure coupling by using the SIMPLE (Semi-Implicit Method for Pressure linked Equations) algorithm. Other than that, the turbulence model using $k-{\varepsilon}$ standard obtained good results. This simulation demonstrated the method of the flow field characteristics, the behavior of velocity vectors and pressure distribution contours in and around the areas of the bucket.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1996년도 춘계학술발표회논문집(2)
• /
• pp.95-101
• /
• 1996

In this paper, the steady 2-dimensional model for a long horizontal line with different end temperatures undergoing natural convection at very high Rayleigh number is proposed to numerically investigate the heat transfer and flow characteristics. The dimensionless governing equations are solved by using SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm which is developed using control volumes and staggered grids. The numerical results are verified by comparison with the operating PWR test data. The analysis focuses on the effects of variation of the heat transfer rates at the pipe surface, the thermal conductivities of the pipe material and the thickness of the pipe wall on the thermal stratification. The results show that the heat transfer rate at the pipe surface is the controlling parameter. A significant reduction and disappearance of thermal stratification phenomenon is observed at the Biot number of 5.0$\times$10$^{-2}$. The results also show that the increment of the thermal conductivity and thickness of the wall weakens the thermal stratification and somewhat reduces azimuthal temperature gradient in the pipe wall. Those effects are however minor, when compared with those due to the variation of the heat transfer rates at the surface of the pipe wall.

• 정보저장시스템학회논문집
• /
• 제5권1호
• /
• pp.19-35
• /
• 2009

The present work is an experimental and analytical study on a flexible disk rotating close to a rigid rotating disk in open air. In the analytical study, the air flow in the gap between the flexible disk and the rigid disk is modeled using Navier-Stokes and continuity equations while the flexible disk is modeled using the linear plate theory. The flow equations are discretized using the cell centered finite volume method (FVM) and solved numerically with semi-implicit pressure-linked equations (SIMPLE algorithm). The spatial terms in the disk equation are discretized using the finite difference method (FDM) and the time integration is performed using fourth-order Runge-Kutta method. An experimental test-rig is designed to investigate the dynamics of the flexible disk when rotating close to a co-rotating, a counter-rotating and a fixed rigid disk, which works as a stabilizer. The effects of rotational speed, initial gap height and inlet-hole radius on the flexible disk displacement and its vibration amplitude are investigated experimentally for the different types of stabilizer. Finally, the analytical and experimental results are compared.

• Nuclear Engineering and Technology
• /
• 제28권1호
• /
• pp.27-35
• /
• 1996

본 논문에서는, 가압경수로 발전소의 가압기 밀림관내 비정상상태의 열성층현상에 대한 계산 모델을 제안하여 배관내의 온도분포, 유동특성 및 열응력에 대해 연구하였다. 경계면이 시간에 따라 변화가 없거나 정상상태에서 개발된 다른 모델과는 달리 본 모델에서는 고온 및 저온유체 사이의 경계면을 시간의 함수로 가정하였다. 열성층현상에 대한 무차원지배방정식은 SIMPLE 알고리즘을 사용하여 해를 구하였다. 본 수치계산의 결과는 주어진 조건하에서 무차원시간이 약 27.0 일 때 배관의 고온부 및 저온부사이의 최대무차원온도차는 0.78정도이었고, 이때의 열성층 현상에 의한 최대 열응력은 276 MPa로 계산되었다.

• 태양에너지
• /
• 제5권2호
• /
• pp.77-83
• /
• 1985

An analysis is performed to investigate the influence of the buoyancy force and the thickness variation of melting layer in the containment that is filled with phase-change Material surrounding a cylindrical heating tube during melting process. The phase-change material is assumed to be initially solid at its phase-change temperature and the remaining solid at any given time is still at the phase-change temperature and neglecting the effect of heat transfer occuring within the solid. At the start of melting process, the thickness of melting layer is assumed to be a stefan-problem and after the starting process, the change of temperature and velocity is calculated using a two dimensional finite difference method. The governing equations for velocity and temperature are solved by a finite difference method which used SIMPLE (Semi Implicit Method Pressure linked Equations) algorithm. Results are presented for a wide range of Granshof number and in accordance with the time increment and it is founded that two dimensional fluid flow occurred by natural convection decreases the velocity of melting process at the bottom of container. The larger the radius of heating tube, the higher heat transfer is occurred in the melting layer.

• Nuclear Engineering and Technology
• /
• 제30권6호
• /
• pp.581-591
• /
• 1998

Numerical analysis was peformed for the two-dimensional turbulent natural convection for a long horizontal line with different end temperatures. The turbulent model has been applied a standard k-$\varepsilon$ two equation model of turbulence similar to that the proposed by the Launder and Spalding. The dimensionless governing equations are solved by using SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm which is developed using control volumes and staggered grids. The numerical results are verified by comparison with the operating PWR test data. The analysis focuses on the effects of variation of the heat transfer rates at the pipe surface, the thermal conductivities of the pipe material and the thickness of the pipe wall on the thermal stratification. The results show that the heat transfer rate at the pipe surface is the controlling parameter for mitigating of thermal stratification in the long horizontal pipe. A significant reduction and disappearance of the thermal stratification phenomenon is observed at the Biot number of 4.82$\times$10$^{-1}$ . The results also show that the increment of the thermal conductivity and thickness of the wall weakens a little the thermal stratification and somewhat reduces temperature gradient of y-direction in the pipe wall. These effects are however minor, when compared with those due to the variation of the heat transfer rates at the surface of the pipe wall.

• 물과 미래
• /
• 제23권3호
• /
• pp.319-328
• /
• 1990

2차원 하구모델의 매개변수 동정을 BF좌표계로 변환된 ADI-FDM출력 수치모형을 사용하여 수행하였다. 모형에서 기본 방정식을 동수역학식과 수송방정식이 복합된 방정식을 사용하였다. Thompson식이 지배방정식을 사각평면으로 변환 하는데 사용도이ㅓㅆ으며 Thompson의 Elliptic 격자발생기법이 BF 좌표계에서 곡선격자망을 발생하는데 사용되었다. 지배방정식에서 동정될 매개변수는 마찰계수와 분산 정수이다. 출력수치모형은 BF좌표계에서 조석평균 염도 모형이 사용되었다. 최적화 기법으로 영향계수 알고리즘을 적용하였다. 또한 Lumped 모형이 동정에 고려 되었다. 본 연구는 새로운 출력 수치모형이 하구염도 모형의 매개변수 동정에 있어 유용한가를 검토 하는데 중점을 두었다. 제안된 기법은 간단한 가상모형에 실험적인 적용을 통하여 검토되었다. 검토한 결과는 제안된 기법이 하구모형에서 매개변수 동정시 도입될 수 있음을 보여주었다.

• Nuclear Engineering and Technology
• /
• 제51권6호
• /
• pp.1487-1503
• /
• 2019

The methods and performance of a pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT employing the drift-flux model are presented. This code aims at providing an accurate yet fast core thermal-hydraulics solution capability to high-fidelity multiphysics core analysis systems targeting massively parallel computing platforms. The four equation drift-flux model is adopted for two-phase calculations, and numerical solutions are obtained by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models involving turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analyses. ESCOT is parallelized by a domain decomposition scheme that involves both radial and axial decomposition to enable highly parallelized execution. The ESCOT solutions are validated through the applications to various experiments which include CNEN $4{\times}4$, Weiss et al. two assemblies, PNNL $2{\times}6$, RPI $2{\times}2$ air-water, and PSBT covering single/two-phase and unheated/heated conditions. The parameters of interest for validation include various flow characteristics such as turbulent mixing, spacer grid pressure drop, cross-flow, reverse flow, buoyancy effect, void drift, and bubble generation. For all the validation tests, ESCOT shows good agreements with measured data in the extent comparable to those of other subchannel-scale codes: COBRA-TF, MATRA and/or CUPID. The execution performance is examined with a mini-sized whole core consisting of 89 fuel assemblies and for an OPR1000 core. It turns out that it is about 1.5 times faster than a subchannel code based on the two-fluid three field model and the axial domain decomposition scheme works as well as the radial one yielding a steady-state solution for the OPR1000 core within 30 s with 104 processors.