• 한국결정성장학회지
• /
• 제14권4호
• /
• pp.155-159
• /
• 2004

열교환법을 활용한 사파이어 단결정 성장 공정에서 도가니 형상 변화가 결정 온도와 고/액 계면 형태에 미치는 영향에 관해 고찰하기 위해 유한요소법, implicit Euler법, frontal 해석 연산을 활용한 수치해석을 수행하였다. 개발된 컴퓨터 시뮬레이션 기법은 고/액 계면의 형상이 반구 형상에서 평면 형상으로 전환되는 열전달 현상 해석에 효율적이다. 본 연구에서는 고/액 계면의 휨도를 개선하기 위해, 도가니 밑면의 다양한 형상을 고려하였으며, 도가니 형상은 공정 최적화 변수로 고려되어야 한다.

• 대한용접접합학회:학술대회논문집
• /
• 대한용접접합학회 2003년도 추계학술발표대회 개요집
• /
• pp.156-158
• /
• 2003

An implicit finite element implementation for Leblond's transformation plasticity constitutive equations, which are widely used in welded steel structure is proposed in the framework of parallel computing. The implementation is based upon the multiplicative decomposition of deformation gradient and hyper elastic formulation. We examine the efficiency of parallel computation for the finite element analysis of a welded structure using domain-wise multi-frontal solver.

• 대한용접접합학회:학술대회논문집
• /
• 대한용접접합학회 2002년도 춘계학술발표대회 개요집
• /
• pp.161-163
• /
• 2002

An implicit finite element implementation for Leblond's transformation plasticity constitutive equations, which are widely used in welded steel structure is proposed in the framework of parallel computing. The implementation is based upon the updated Lagrangian formulation. We examine the efficiency of parallel compuatation for the finite element analysis of a welded structure using multi-frontal solver.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 춘계 학술대회논문집
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

• 한국항공우주학회지
• /
• 제30권8호
• /
• pp.37-45
• /
• 2002

본 연구에서는 이차원에서 비정상 비점성 유동해석을 위한 비정렬 동적 편자 기법을 개발하였다. 유동해석 기법은 시간에 대해 2차의 정확도를 갖는 내재적인 시간적분법을 사용하였으며, 격자중심의 유한 체적법과 Roe의 풍상차분법을 이용하여 공간에 대한 차분화를 하였다. 시간과 공간에 대한 정확도를 증가시키기 위해서는 해에 따라 원하는 위치에 격자점들을 임의로 추가할 수 있는 비정상 동적 적응격자 기법을 사용하였다. 이를 이용하여 이차원의 2자유도를 갖는 스프링 에어포일 시스템의 와류와의 간섭현상에 따른 공탄성적 변위를 예측하였다.

• Nuclear Engineering and Technology
• /
• 제52권3호
• /
• pp.461-468
• /
• 2020

The improvements of the DeCART/MUSAD code system for uncertainty analysis of HTGR neutronic parameters are presented in this paper. The function for quantifying an uncertainty of critical-spectrumweighted few group cross section was implemented using the generalized adjoint B₁ equation solver. Though the changes between the infinite and critical spectra cause a considerable difference in the contribution by the graphite scattering cross section, it does not significantly affect the total uncertainty. To reduce the number of iterations of the generalized adjoint transport equation solver, the generalized adjoint B₁ solution was used as the initial value for it and the number of iterations decreased to 50%. To reflect the implicit uncertainty, the correction factor was derived with the resonance integral. Moreover, an additional correction factor for the double heterogeneity was derived with the effective cross section of the DH region and it reduces the difference from the complete uncertainty. The code system was examined with the MHTGR-350 Ex.II-2 3D core benchmark. The k_eff uncertainty for Ex.II-2a with only the fresh fuel block was similar to that of the block and the uncertainty for Ex.II-2b with the fresh fuel and the burnt fuel blocks was smaller than that of the fresh fuel block.

• 한국전산유체공학회지
• /
• 제19권2호
• /
• pp.17-23
• /
• 2014

Cavitation causes a great deal of noise, damage to components, vibrations, and a loss of efficiency in devices, such as propellers, pump impellers, nozzles, injectors, torpedoes, etc. Thus, the cavitating flow simulation is of practical importance for many engineering systems. In the present work, a two-phase flow solver based on the homogeneous mixture model has been developed. The solver employs an implicit preconditioning, dual time stepping algorithm in curvilinear coordinates. The flow characteristics around Clark-Y hydrofoil were calculated and then validated by comparing with the experimental data. The lift and drag coefficients with changes of angle of attack and cavitation number were obtained. The results show that cavity length and lift, drag coefficient increase with increasing angle of attack.

• Journal of Mechanical Science and Technology
• /
• 제20권11호
• /
• pp.1961-1971
• /
• 2006

A soft recovery system (SRS) is a device that stops a high speed projectile without damaging the projectile. The SRS is necessary to verify the shock resistant requirements of microelectronics and electro-optic sensors in smart munitions, where the projectiles experience over 20,000 g acceleration inside the barrel. In this study, a computer code for the performance evaluation of a SRS based on ballistic compression decelerator concept has been developed. It consists of a time accurate compressible one-dimensional Euler code with use of deforming grid and a projectile motion analysis code. The Euler code employs Roe's approximate Riemann solver with a total variation diminishing (TVD) method. A fully implicit dual time stepping method is used to advance the solution in time. In addition, the geometric conservation law (GCL) is applied to predict the solutions accurately on the deforming mesh. The equation of motion for the projectile is solved with the four-stage Runge-Kutta time integration method. A small scale SRS to catch a 20 mm bullet fired at 500 m/s within 1,600 g-limit has been designed with the proposed method.

• 한국신재생에너지학회:학술대회논문집
• /
• 한국신재생에너지학회 2009년도 추계학술대회 논문집
• /
• pp.600-605
• /
• 2009

A general purpose viscous flow solver Ansys CFX is used to study a Savonius type wave energy converter in a 3D numerical viscous wave tank. This paper presents the results of a computational fluid dynamics (CFD) analysis of the effect of blade configuration on the performance of 3 bladed Savonius rotors for wave energy extraction. A piston-type wave generator was incorporated in the computational domain to generate the desired incident waves. A complete OWC system with a 3-bladed Savonius rotor was modeled in a three dimensional numerical wave tank and the hydrodynamic conversion efficiency was estimated. The flow over the rotors is assumed to be two-dimensional (2D), viscous, turbulent and unsteady. The CFX code is used with a solver of the coupled conservation equations of mass, momentum and energy, with an implicit time scheme and with the adoption of the hexahedral mesh and the moving mesh techniques in areas of moving surfaces. Turbulence is modeled with the k.e model. Simulations were carried out simultaneously for the rotor angle and the helical twist. The results indicate that the developed models are suitable to analyze the water flows both in the chamber and in the turbine. For the turbine, the numerical results of torque were compared for all the cases.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제17권3호
• /
• pp.197-207
• /
• 2013

The Cahn-Hilliard equation was proposed as a phenomenological model for describing the process of phase separation of a binary alloy. The equation has been applied to many physical applications such as amorphological instability caused by elastic non-equilibrium, image inpainting, two- and three-phase fluid flow, phase separation, flow visualization and the formation of the quantum dots. To solve the Cahn-Hillard equation, many numerical methods have been proposed such as the explicit Euler's, the implicit Euler's, the Crank-Nicolson, the semi-implicit Euler's, the linearly stabilized splitting and the non-linearly stabilized splitting schemes. In this paper, we investigate each scheme in finite-difference schemes by comparing their performances, especially stability and efficiency. Except the explicit Euler's method, we use the fast solver which is called a multigrid method. Our numerical investigation shows that the linearly stabilized stabilized splitting scheme is not unconditionally gradient stable in time unlike the known result. And the Crank-Nicolson scheme is accurate but unstable in time, whereas the non-linearly stabilized splitting scheme has advantage over other schemes on the time step restriction.