• 한국산학기술학회논문지
• /
• 제15권6호
• /
• pp.3950-3960
• /
• 2014

본 연구에서는 수치모의를 통해서 월류 흐름이 존재하는 수몰 사각 실린더의 항력 특성에 대하여 분석하였다. 모의의 신뢰성을 검토하기 위하여 실험자료와 비교하였으며 실험에서 측정하기 어려운 실린더 접촉면의 압력에 대한 분석을 통해서 상대 수심에 따른 항력의 특성을 분석하였다. 3차원 동수역학 모형을 이용한 수몰 사각 실린더의 항력 계산 결과는 실험자료의 상대 수심의 변화에 따른 항력계수의 변화를 유사하게 모의하고 있음을 확인할 수 있었다. 수치모의 결과 분석에 의하면 수몰 사각 실린더에 작용하는 항력은 대부분 압력이며 상대 수심이 증가함에 따라 전단력의 비중은 감소하였다. 실린더 접촉면의 압력계수 분석 결과에 의하면 상대 수심이 낮은 경우에는 전면부에 높은 압력계수가 형성되고 후면부에 낮은 압력계수가 형성되어 결과적으로 높은 항력계수가 나타남을 확인하였다. 상대수심이 증가하면 전면부의 압력계수는 감소하고 후면부의 압력계수는 증가하여 2차원 흐름 내의 사각 실린더와 유사한 양상을 나타낸다. 정수압 영향 분석에 의하면 전면부와 후면부의 수위 차에 의한 정수압은 항력에 미치는 영향이 제한적이며 사각 실린더에 의해 형성되는 국부적인 수위와 함께 3차원적인 흐름에 의해 형성되는 동수압의 영향이 크다는 것을 확인하였다.

• 한국진공학회:학술대회논문집
• /
• 한국진공학회 2014년도 제46회 동계 정기학술대회 초록집
• /
• pp.221-221
• /
• 2014

The present work proposes an improved numerical simulator for design and modification of large area capacitively coupled plasma (CCP) processing chamber. CCP, as notoriously well-known, demands the tremendously huge computational cost for carrying out transient analyses in realistic multi-dimensional models, because electron dissociations take place in a much smaller time scale (${\Delta}t{\approx}10-8{\sim}10-10$) than time scale of those happened between neutrals (${\Delta}t{\approx}10-1{\sim}10-3$), due to the rf drive frequencies of external electric field. And also, for spatial discretization of electron flux (Je), exponential scheme such as Scharfetter-Gummel method needs to be used in order to alleviate the numerical stiffness and resolve exponential change of spatial distribution of electron temperature (Te) and electron number density (Ne) in the vicinity of electrodes. Due to such computational intractability, it is prohibited to simulate CCP deposition in a three-dimension within acceptable calculation runtimes (<24 h). Under the situation where process conditions require thickness non-uniformity below 5%, however, detailed flow features of reactive gases induced from three-dimensional geometric effects such as gas distribution through the perforated plates (showerhead) should be considered. Without considering plasma chemistry, we therefore simulated flow, temperature and species fields in three-dimensional geometry first, and then, based on that data, boundary conditions of two-dimensional plasma discharge model are set. In the particular case of SiH4-NH3-N2-He CCP discharge to produce deposition of SiNxHy thin film, a cylindrical showerhead electrode reactor was studied by numerical modeling of mass, momentum and energy transports for charged particles in an axi-symmetric geometry. By solving transport equations of electron and radicals simultaneously, we observed that the way how source gases are consumed in the non-isothermal flow field and such consequences on active species production were outlined as playing the leading parts in the processes. As an example of application of the model for the prediction of the deposited thickness uniformity in a 300 mm wafer plasma processing chamber, the results were compared with the experimentally measured deposition profiles along the radius of the wafer varying inter-electrode gap. The simulation results were in good agreement with experimental data.

• Computers and Concrete
• /
• 제2권4호
• /
• pp.249-266
• /
• 2005

In the present paper a transient three-dimensional thermo-mechanical model for concrete is presented. For given boundary conditions, temperature distribution is calculated by employing a three-dimensional transient thermal finite element analysis. Thermal properties of concrete are assumed to be constant and independent of the stress-strain distribution. In the thermo-mechanical model for concrete the total strain tensor is decomposed into pure mechanical strain, free thermal strain and load induced thermal strain. The mechanical strain is calculated by using temperature dependent microplane model for concrete (O$\check{z}$bolt, et al. 2001). The dependency of the macroscopic concrete properties (Young's modulus, tensile and compressive strengths and fracture energy) on temperature is based on the available experimental database. The stress independent free thermal strain is calculated according to the proposal of Nielsen, et al. (2001). The load induced thermal strain is obtained by employing the biparabolic model, which was recently proposed by Nielsen, et al. (2004). It is assumed that the total load induced thermal strain is irrecoverable, i.e., creep component is neglected. The model is implemented into a three-dimensional FE code. The performance of headed stud anchors exposed to fire was studied. Three-dimensional transient thermal FE analysis was carried out for three embedment depths and for four thermal loading histories. The results of the analysis show that the resistance of anchors can be significantly reduced if they are exposed to fire. The largest reduction of the load capacity was obtained for anchors with relatively small embedment depths. The numerical results agree well with the available experimental evidence.

• 지질공학
• /
• 제33권1호
• /
• pp.151-167
• /
• 2023

본 연구는 LNG저장탱크가 설치될 느슨한 포화사질지반을 대상으로 개정된 액상화 평가법과 UBC3D-PLM 모델을 이용한 1차원 유효응력해석에 의한 액상화 평가법을 비교한 것이다. 이를 위해 여러가지 실내 및 현장시험을 실시하여 필요한 Parameter를 산정하였다. 검토결과, 지진응답해석결과와 SPT N 값을 이용하는 개정 액상화 평가법은 액상화 발생 가능성을 높게 평가하였지만, 다양한 액상화 저항인자를 고려할 수 있는 유효응력해석법은 액상화에 다소 안정한 것으로 분석되었다. UBC3D-PLM 모델을 이용한 1차원 유한요소해석을 할 경우 보다 간편하게 액상화 안정성 검토가 가능하였고, 액상화 보강 영역을 최소화 할 수 있었다. 또한, LNG저장탱크의 기초를 고려한 2·3차원 수치해석 시에는 액상화 발생 시 내진설계 및 거동특성을 규명하는 것에 활용될 수 있을 것으로 기대된다.

• 대한기계학회논문집A
• /
• 제26권3호
• /
• pp.576-585
• /
• 2002

A three-dimensional Lagrangian explicit time-integration finite element code for analyzing the dynamic impact phenomena was developed. It uses four node tetrahedral elements. In order to consider the effects of strain rate hardening, strain hardening and thermal softening, which are frequently observed in high-velocity deformation phenomena, Johnson-Cook model is used as constitutive model. For more accurate and robust contact force computation, the defense node contact algorithm was adopted and implemented. In order to evaluate the performance of the newly developed three-dimensional hydrocode NET3D, numerical simulations of the oblique impact of mild steel plate by mild steel sphere were carried out. Ballistic limit about various oblique angle between 0 degree and 80 degree was estimated through a series of simulations with different initial velocities of sphere. Element eroding by equivalent plastic strain was applied to mild steel spheres and targets. Ballistic limits and fracture characteristics obtained from simulation were compared with experimental results conducted by Finnegan et al. From numerical studies, the following conclusions were reached. (1) Simulations could successfully reproduce the key features observed in experiment such as tensile failure termed "disking"at normal impacts and outwards bending of partially formed plus segments termed "hinge-mode"at oblique impacts. (2) Simulation results fur 60 degrees oblique impact at 0.70 km/s and 0.91 km/s were compared with experimental results and Eulerian hydrocode CTH simulation results. The Lagrangian code NET3D is superior to Eulerian code CTH in the computational accuracy. Agreement with the experimentally obtained final deformed cross-sections of the projectile is excellent. (3) Agreement with the experimental ballistic limit data, particularly at the high-obliquity impacts, is reasonably good. (4) The simulation result is not very sensitive to eroding condition but slightly influenced by friction coefficient.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.31-40
• /
• 2011

Unsteady sheet cavitation on a three-dimensional twisted hydrofoil was studied using an unsteady Reynolds-averaged Navier-Stokes equations solver based on a cell-centered finite volume method. As a verification test of the computational method. non-cavitating and cavitating flow over a modified NACA66 foil section was simulated and validated against existing experimental data. The numerical uncertainties of forces and pressure were evaluated for three levels of mesh resolution. The computed pressure on the foil and the cavity shedding behavior were validated by comparing with existing experimental data. The cavity shedding dynamics by re-entrant jets from the end and sides of the cavity were investigated.

• International Journal of Fluid Machinery and Systems
• /
• 제4권1호
• /
• pp.133-139
• /
• 2011

In this paper, Based on Two-dimensional Flow Theory, adopting quasi-orthogonal method and point-by-point integration method to design the impeller of the low specific speed centrifugal pump by code, and using RANS (Reynolds Averaged N-S) Equation with a standard k-${\varepsilon}$ two-equation turbulence model and log-law wall function to solve 3D turbulent flow field in the impeller of the low specific speed pump. An analysis of the influences of the blade profile on velocity distributions, pressure distributions and pump performance and the investigation of the flow regulation pattern in the impeller of the centrifugal pump are presented. And the result shows that this method can be used as a new way in low speed centrifugal pump impeller design.

• 한국광학회지
• /
• 제31권2호
• /
• pp.59-70
• /
• 2020

3차원 플래시 라이다 시스템(3D flash LIDAR system)에서의 대기 산란을 해석하기 위해 몬테 카를로 복사 전달(Monte Carlo radiative transfer, MCRT) 방법을 바탕으로 수정된 수치 해석 모델인 MCRT 행렬 방법을 논의한다. MCRT 방법을 바탕으로 라이다 신호의 복사 전달 함수를 행렬 형태로 구성하며, 이는 특성 응답에 해당한다. 근축 근사에 기반하여 본 특성 응답 행렬의 중첩 및 합성곱 연산을 활용함으로써 확장된 전반적인 플래시 라이다의 전산 모사 모델을 개발한다. MCRT 행렬 방법은 기존의 몬테 카를로 기반 방법들에서 과도하게 증가할 수 있는 개별 라이다 신호의 추적을 대폭 경감시킨다. 그 결과 본 방법은 다양한 산란 조건 및 라이다 시스템 구성 환경에서도 그 신호 응답을 빠르게 획득할 수 있는 특징을 지닌다. 본 논문에서는 MCRT 행렬 방법에 기반한 전산 모델을 이용하여 상이한 대기 환경 조건에서 동작하는 3차원 플래시 라이다 시스템을 그 산란 조건, 즉, 그 가시거리에 따른 산란 계수를 달리하며 모사하고, 플래시 라이다 신호의 신호대잡음비의 악화, 신호 오류, 시공간적 확산 및 시간 지연 등 시스템상에서의 산란 효과에 의해 나타나는 다양한 현상을 수치적으로 분석한다. MCRT 행렬 방법은 자율 주행을 위한 플래시 라이다 시스템을 포함해 다양한 라이다 시스템을 분석하는데 매우 효과적으로 사용될 수 있을 것으로 기대된다.

• 한국해안·해양공학회논문집
• /
• 제20권1호
• /
• pp.1-13
• /
• 2008

본 연구에서는 복수의 대형 연직 주상구조물을 중심으로 구조물에 작용하는 비선형파력과 구조물에 의한 비선형파랑변형을 검토하며, 수치해석의 이론으로 쇄파나 강비선형 파랑을 고정도로 해석할 수 있는 Navier-Stokes방정식에 근거한 3차원 VOF(Volume Of Fluid)법을 적용하였다. 본 수치해석의 타당성을 검증하기 위해 파랑변형과 파력에 대한 기존의 수리실험결과와 본 수치해석결과를 비교 검토하였으며, 이로부터 2기로 구성된 임의형상의 연직 주상구조물에 대해 구조물 형상, 구조물 배치간격, 파의 입사각도 변화에 따른 파력과 파랑변형의 특성을 논하였다.

• 한국콘크리트학회:학술대회논문집
• /
• 한국콘크리트학회 2001년도 봄 학술발표회 논문집
• /
• pp.431-436
• /
• 2001

Lateral confining effect due to the existence of the shear reinforcements in R.C. beam is investigate in a numerical way. For the purpose, a three dimensional constitutive model of concrete is developed based on the elasto-plasticity using non-associated plastic flow rule to control the excessive inelastic dilatancy. The plastic flow direction is determined based on the associated plastic flow direction in a way to adjust the directional angle between the two normal vector components along the hydrostatic and deviatoric axis in a meridian plane in which the loading function prescribed. The current formulation is combined with the four parameter elasto-plastic triaxial concrete model recently developed. The resulting elasto-plastic triaxial concrete model predicts the fundamental behaviors of concrete under different confining levels and the 4-points flexural test of a beam with shear reinforcements, compares with the experimental results.