• Advances in Energy Research
• /
• 제5권2호
• /
• pp.91-105
• /
• 2017

Thermal hydraulic (TH) analysis of nuclear power reactors is utmost important. In this way, the numerical codes that preparing TH data in reactor core are essential. In this paper, a subchannel analysis of a Russian pressurized water reactor (WWER1000) core with enhanced numerical code is carried out. For this, in fluid domain, the mass, axial and lateral momentum and energy conservation equations for desired control volume are solved, numerically. In the solid domain, the cylindrical heat transfer equation for calculation of radial temperature profile in fuel, gap and clad with finite difference and finite element solvers are considered. The dependence of material properties to fuel burnup with Calza-Bini fuel-gap model is implemented. This model is coupled with Isotope Generation and Depletion Code (ORIGEN2.1). The possibility of central hole consideration in fuel pellet is another advantage of this work. In addition, subchannel to subchannel and subchannel to rod connection data in hexagonal fuel assembly geometry could be prepared, automatically. For a demonstration of code capability, the steady state TH analysis of a the WWER1000 core is compromised with Thermal-hydraulic analysis code (COBRA-EN). By thermal hydraulic parameters averaging Fuel Assembly-to-Fuel Assembly method, the one sixth (symmetry) of the Boushehr Nuclear Power Plant (BNPP) core with regular subchannels are modeled. Comparison between the results of the work and COBRA-EN demonstrates some advantages of the presented code. Using the code the thermal modeling of the fuel rods with considering the fission gas generation would be possible. In addition, this code is compatible with neutronic codes for coupling. This method is faster and more accurate for symmetrical simulation of the core with acceptable results.

• 한국공간구조학회:학술대회논문집
• /
• 한국공간구조학회 2008년도 춘계 학술발표회 논문집
• /
• pp.107-112
• /
• 2008

본 연구에서는 이산계열 대공간구조물의 크기최적화에 따른 후 좌굴거동의 변화에 대하여 조사하고 그 결과를 기술하였다. 본 연구에서는 35개의 부재를 가지는 얕은 아치형 2차원 공간트러스 구조물을 이용하여 연구를 수행하였다. 먼저 트러스의 최적의 부재 크기 패턴을 조사하기 위하여 수학적 프로그래밍 기법을 도입하였다. 이때 최소화해야하는 트러스의 전체 부재의 중량을 목적함수로 이용하고 하중이 가해지는 절점에서 발생하는 변위 값과 각 부재에서 발생하는 응력 값을 허용치 이하로 제한하는 제약조건으로 사용하였다. 크기최적화에서 도출된 트러스 구조물의 후좌굴 거동은 통합 비선형해석기법을 이용하여 조사 분석하였다. 본 연구에서 제시한 결과는 이산계열 대공간구조물의 설계에 기본 데이터로 유용하게 사용될 것으로 판단된다.

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

In this paper, prediction of separation trajectory for Two-stage-To-Orbit space launch vehicle has been numerically simulated by using an aerodynamic database based on steady state analysis. Aerodynamic database were obtained for matrix of longitudinal and vertical positions. The steady flow simulations around the launch vehicle have been made by using a 3-D RANS flow solver based on unstructured meshes. For this purpose, a vertex-centered finite-volume method was adopted to discretize inviscid and viscous fluxes. Roe's finite difference splitting was utilized to discretize the inviscid fluxes, and the viscous fluxes were computed based on central differencing. To validate this flow solver, calculations were made for the wind-tunnel experiment model of the LGBB TSTO vehicle configuration on steady state conditions. Aerodynamic database was constructed by using flow simulations based on test matrix from the wind-tunnel experiment. ANN(Artificial Neural Network) was applied to construct interpolation function among aerodynamic variables. Separation trajectory for TSTO launch vehicle was predicted from 6-DOF equation of motion based on the interpolated function. The result of present separation trajectory calculation was compared with the trajectory using experimental database. The predicted results for the separation trajectory shows fair agreement with reference[4] solution.

• 한국가스학회지
• /
• 제12권2호
• /
• pp.57-62
• /
• 2008

배관 엘보우의 내호면(intrados)의 과 내부에 국부적으로 두께 감육이 발생한 경우, 내압과 엘보우를 닫는 방향으로의 굽힘하중을 부가하여 파손 모드를 연구하였다. 탄소성해석 시 반력-변위 곡선이 세 그룹으로 나뉘므로 각 그룹의 한 경우씩을 해석하여 소성붕괴에 의한 파손모드의 차이를 확인하였다. 이를 위해 주요 부위에서 하중-국부적응력 곡선이 어떻게 변화하는지 결정하여, 이로부터 관찰된 파손모드와 비교하여 설명하였다. 감육폭이 $90^{\circ}$인 경우 배관은 엘보우 측면부터 소성붕괴가 시작되었으며, $360^{\circ}$인 경우 내호면으로부터 소성변형이 시작되어 서로 다른 파손모드를 보여주었다. 배관의 감육측정에 의한 건전성 평가 시 이와 같은 파손 모드의 차이점을 고려하여 평가를 실시하여야 한다.

• 한국전산유체공학회지
• /
• 제16권4호
• /
• pp.7-13
• /
• 2011

A comparative analysis of thermal models in the lattice Boltzmann method(LBM) for the simulation of laminar natural convection in a square cavity is presented. A HYBRID method, in which the thermal equation is solved by the Navier-Stokes equation method while the mass and momentum conservation are resolved by the lattice Boltzmann method, is introduced and its merits are explained. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with a deferred correction method to ensure stability of the solutions. The HYBRID method and the double-population method are applied to the simulation of natural convection in a square cavity and the predicted results are compared with the benchmark solutions given in the literatures. The predicted results are also compared with those by the conventional Navier-Stokes equation method. In general, the present HYBRID method is as accurate as the Navier-Stokes equation method and the double-population method. The HYBRID method shows better convergence and stability than the double-population method. These observations indicate that this HYBRID method is an efficient and economic method for the simulation of incompressible fluid flow and heat transfer problem with the LBM.

• 전자공학회논문지CI
• /
• 제47권4호
• /
• pp.82-87
• /
• 2010

CPU를 능가하는 GPU의 연산능력 향상으로 범용 계산에 그래픽 프로세서를 사용하는 GP-GPU연구가 활발히 전개되고 있으며, 그 응용분야가 확대되고 있다. 본 논문에서는 전자기학 관련 분야에서 널리 사용되는 FDTD 알고리즘을 nVIDIA에서 제공하는 소프트웨어 플랫폼인 CUDA를 사용하여 구현한다. FDTD 알고리즘의 주요 연산과정을 병렬화하고, 그래픽 카드 내각기 다른 메모리의 사용에 따라 최적화하며, 단일 프로세서에서 FDTD 알고리즘을 실행시킨 경우와 비교하여 그 성능 향상 정도를 측정한다. 실험결과 단일 프로세서로 구현하였을 때에 비해 실행시간이 45배까지 향상됨을 확인할 수 있었다.

• 대한전기학회논문지:전기물성ㆍ응용부문C
• /
• 제51권11호
• /
• pp.551-558
• /
• 2002

A time-domain combined field integral equation (CFIE) is presented to obtain the transient scattering response from arbitrarily shaped three-dimensional conducting bodies. This formulation is based on a linear combination of the time-domain electric field integral equation (EFIE) with the magnetic field integral equation (MFIE). The time derivative of the magnetic vector potential in EFIE is approximated using a central finite difference approximation and the scalar potential is averaged over time. The time-domain CFIE approach produces results that are accurate and stable when solving for transient scattering responses from conducting objects. The incident spectrum of the field may contain frequency components, which correspond to the internal resonance of the structure. For the numerical solution, we consider both the explicit and implicit scheme and use two different kinds of Gaussian pulses, which may contain frequencies corresponding to the internal resonance. Numerical results for the EFIE, MFIE, and CFIE are presented and compared with those obtained from the inverse discrete Fourier transform (IDFT) of the frequency-domain CFIE solution.

• 한국전산유체공학회지
• /
• 제17권1호
• /
• pp.94-100
• /
• 2012

In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

• International Journal of Aeronautical and Space Sciences
• /
• 제12권4호
• /
• pp.318-330
• /
• 2011

This paper investigates the convergence characteristics of the modified artificial compressibility method proposed by Turkel. In particular, a focus is mode on the convergence characteristics due to variation of the preconditioning factor (${\alpha}_u$) and the artificial compressibility (${\beta}$) in conjunction with an upwind method. For the investigations, a code using the modified artificial compressibility is developed. The code solves the axisymmetric incompressible Reynolds averaged Navier-Stokes equations. The cell-centered finite volume method is used in conjunction with Roe's approximate Riemann solver for the inviscid flux, and the central difference discretization is used for the viscous flux. Time marching is accomplished by the approximated factorization-alternate direction implicit method. In addition, Menter's k-${\omega}$ shear stress transport turbulence model is adopted for analysis of turbulent flows. Inviscid, laminar, and turbulent flows are solved to investigate the accuracy of solutions and convergence behavior in the modified artificial compressibility method. The possible reason for loss of robustness of the modified artificial compressibility method with ${\alpha}_u$ >1.0 is given.

• 대한조선학회지
• /
• 제21권3호
• /
• pp.43-50
• /
• 1984

A body form, which experiences minimum vertical wave-exciting forces in the vicinity of a prescribed wave frequency in water of finite depth, is obtained by an approximate method. Its configuration has the symmetry with respect to the vertical axis, expressed in terms of exponential functions. By distributing three-dimensional pulsating sources and dipoles on the immersed surface of the body, a velocity potential is determined and subsequently hydrodynamic forces including the 2nd-order time-mean drift forces are calculated. The dynamic behavior of the body moored in irregular waves is investigated numerically by using central difference method. Hereby irregular wave trains are simulated with examining its repeatability by comparing the resulting spectrum with original one. Numerical results indicated that the body form obtained from the present analysis possesses in general a favorable hydrodynamic characteristics in comparison with a spherical buoy and that the maximum excursion of the body can be significantly reduced by setting pre-tension of an appropriate amount in the mooring cable.