• 대한기계학회논문집B
• /
• 제20권3호
• /
• pp.1015-1027
• /
• 1996

In this paper, a numerical procedure for the calculation of the incompressible Navier-Stokes equations in a generalized nonorthogonal body fitted coordinate system is proposed and is validated through three test problems. Present numerical procedure derives the pressure equation by using the pressure substitution method on the regular grid system, and discretized momentum equations are based on the covariant velocity components. Cavity flow, backward facing step flow, and two dimensional channel flow with a sinusoidal wavy wall are chosen as three test problems. Numerical solutions obtained by present procedure shows a good agreement with previous numerical and/or experimental results. Convergence rate is also satisfactory.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.22-24
• /
• 2008

This paper presents a numerical model of internal flows in a lubricant supplying and scavenging flow path of an aero-engine lubrication system. The numerical model was built in the General Analysis Software of Aero-engine Lubrication System, GASLS, developed by Northwestern Polytechnical University. The lubricant flow flux, pressure and temperature distribution at steady state were calculated. GASLS is a general purpose computer program employed a 1-D steady state network algorithm for analyzing flowrates, pressures and temperatures in a complex flow network. All kinds of aero-engine lubrication systems can be divided into finite correlative typical elements and nodes from which the calculation network be developed in GASLS. Special emphasis is on how to use combined elements which is a type of typical elements to replace some complex components like bearing bores, accessory gearboxes or heat exchangers. This method can reduce network complexity and improve calculation efficiency. Final computational results show good agreement with experimental data.

• 한국습지학회지
• /
• 제23권2호
• /
• pp.107-115
• /
• 2021

서해 개야수로 인근의 파랑전파, 해수유동, 퇴적물이동을 분석하기 위해서 각종 인공구조물 설치 전(CASE1W)과 후(CASE2W)로 분류하고, CASE1W와 CASE2W에 대한 계산결과를 비교하였다. 파랑전파에 대해서는 SWAN 수치모형을 이용하여 입사파와 반사파의 결과를 도출하였고, 해수유동에 대해서는 FLOW2DH 수치모형을 이용하여 해수유동에 따른 유속 결과를 도출하였다. SWAN 수치모형과 FLOW2DH 수치모형의 결과는 퇴적물이동을 예측하는 SEDTRAN 수치모형의 입력조건이 되어 개야수로 인근의 최대 저면전단응력과 부유사 농도분포를 계산하였다. SWAN 수치모형 계산결과, CASE2W의 경우 약 7 km 길이의 북측 도류제에 의해서 입사파가 회절 및 중첩되고, 반사파가 생성되어 개야수로 인근의 파고를 CASE1W에 비해 40~50 % 증가시켰다. FLOW2DH 수치모형 계산결과, 북방파제, 북측 도류제 및 금란도에 의해서 개야수로의 유속이 CASE1W과 대비하여 CASE2W가 10~30 % 빠르게 계산되었다. SEDTRAN 수침모형의 계산결과, 복합 파랑장(입사파, 반사파, 조석)에 따른 해양환경과 각종 인공구조물의 설치에 의해서 개야수로의 최대 저면전단응력이 1.0 N/m² 이상인 구간과 부유사농도가 80 mg/L 이상인 구간이 넓게 분포되었다는 것은 개야수로에 퇴적현상이 발생한 것이라고 판단된다.

• KIEAE Journal
• /
• 제16권6호
• /
• pp.159-166
• /
• 2016

Purpose: The architectural structures in the engineering field include more than one material, and the heat transfer through these multiple materials becomes complicated. More or less, the analytic solutions obtained by the hand calculation can provide the limited information of heat transfer phenomena. However, the engineers have generally been forced to obtain reliable results than those of the hand calculation. The numerical calculation such as a finite volume methods with the unstructured grid system is only the suitable means of the analysis for the complex and arbitrary domains that consists of multiple materials. In this study, a new numerical code is developed to provide temperature distributions in the multiple material domains, and the results of this code are compared with the validation cases in ISO10211. Method: Finite volume methods with an unstructured grid is employed. In terms of numerical methods, the heat transfer conduction coefficient is not defined on the surface of the cell between different material cells. The heat transfer coefficient is properly defined to accurately mimic the heat transfer through the surface. The boundary conditions of heat flux considering radiation or heat convection are also developed. Result: The comparison between numerical results and ISO 10211 cases. We are confirmed that the numerical method provides the proper temperature distributions, and the heat transfer equation and its boundary conditions are developed properly.

• 한국초전도ㆍ저온공학회논문지
• /
• 제17권2호
• /
• pp.45-49
• /
• 2015

It has been well known that a toroid is the inevitable shape for a high temperature superconducting (HTS) coil as a component of a large scale superconducting magnetic energy storage system (SMES) because it is the best option to minimize a magnetic field intensity applied perpendicularly to the HTS wires. Even though a perfect toroid coil does not have a perpendicular magnetic field, for a practical toroid coil composed of many HTS pancake coils, some type of perpendicular magnetic field cannot be avoided, which is a major cause of degradation of the HTS wires. In order to suggest an optimum design solution for an HTS SMES system, we need an accurate, fast, and effective calculation for the magnetic field, mechanical stresses, and stored energy. As a calculation method for these criteria, a numerical calculation such as an finite element method (FEM) has usually been adopted. However, a 3-dimensional FEM can involve complicated calculation and can be relatively time consuming, which leads to very inefficient iterations for an optimal design process. In this paper, we suggested an intuitive and effective way to determine the maximum magnetic field intensity in the HTS coil by using an analytic and statistical calculation method. We were able to achieve a remarkable reduction of the calculation time by using this method. The calculation results using this method for sample model coils were compared with those obtained by conventional numerical method to verify the accuracy and availability of this proposed method. After the successful substitution of this calculation method for the proposed design program, a similar method of determining the maximum mechanical stress in the HTS coil will also be studied as a future work.

• 화약ㆍ발파
• /
• 제36권3호
• /
• pp.10-18
• /
• 2018

동적 하중을 받는 암반의 역학적 거동은, 같은 크기의 최대 하중이라도 정적으로 가해지는 경우와는 다른 특성을 보인다. 동적 하중 하에서의 거동 특성을 규명하기 위한 실험적 접근 방법은 동적 하중의 제어와 계측 및 해석에 있어서 정적 조건에서의 실험 방법보다 더 많은 어려움이 있다. 수치해석적 방법은 물리적 실험이 아니라 수치해석적으로 실험을 실시함으로써 물리적 제약을 덜 받으므로 설계 단계에서 매우 유력한 해석 도구가 될 수 있다. 그러나 수치해석방법은 해석방법의 알고리즘이 적정하더라도, 입력 자료와 경계조건의 설정에 따라 계산 결과가 많이 달라질 수 있으므로 해석 시 세심한 주의가 필요하다. 본 논문에서는 동적 하중을 받는 암반 구조물의 거동을 수치해석적으로 검토할 때, 경계 조건, 동적 하중과 계산 시간 간격, 동적 하중 특성이 계산 결과에 미치는 영향을 검토하여, 동적 해석 시 경계조건과 계산 시간 간격의 설정 지침을 제공하고자 하였다.

• Structural Engineering and Mechanics
• /
• 제64권1호
• /
• pp.81-92
• /
• 2017

Many materials in engineering exhibit different modulus in tension and compression, which are known as bi-modulus materials. Based on the bi-modulus elastic theory, a modified semi-analytical model, by introducing a stress function, is established in this paper to study the mechanical response of a bi-modulus cylinder placed in an axisymmetric temperature field. Meanwhile, a numerical procedure to calculate the temperature stresses in bi-modulus structures is developed. It is proved that the bi-modulus solution can be degenerated to the classical same modulus solution, and is in great accordance with the solutions calculated by the semi-analytical model proposed by Kamiya (1977) and the numerical solutions calculated both by the procedure complied in this paper and by the finite element software ABAQUS, which demonstrates that the semi-analytical model and the numerical procedure are accurate and reliable. The result shows that the modified semi-analytical model simplifies the calculation process and improves the speed of computation. And the numerical procedure simplifies the modeling process and can be extended to study the stress field of bi-modulus structures with complex geometry and boundary conditions. Besides, the necessity to introduce the bi-modulus theory is discussed and some suggestions for the qualitative analysis and the quantitative calculation of such structure are proposed.

• Nuclear Engineering and Technology
• /
• 제56권3호
• /
• pp.1002-1012
• /
• 2024

To remove insoluble fission products, which could possibly cause reactor instability and significantly reduce heat transfer efficiency from primary system of molten salt reactor, a helium bubbling method is employed into a passive molten salt fast reactor. In this regard, two-phase flow behavior of molten salt and helium bubbles was investigated experimentally because the helium bubbles highly affect the circulation performance of working fluid owing to an additional drag force. As the helium flow rate is controlled, the change of key thermal-hydraulic parameters was analyzed through a two-phase experiment. Simultaneously, to assess the applicability of numerical model for the analysis of two-phase flow behavior, the numerical calculation was performed using the OpenFOAM 9.0 code. The accuracy of the numerical analysis code was evaluated by comparing it with the experimental data. Generally, numerical results showed a good agreement with the experiment. However, at the high helium injection rates, the prediction capability for void fraction of helium bubbles was relatively low. This study suggests that the multiphaseEulerFoam solver in OpenFOAM code is effective for predicting the helium bubbling but there exists a room for further improvement by incorporating the appropriate drag flux model and the population balance equation.

• 조명전기설비학회논문지
• /
• 제27권1호
• /
• pp.46-52
• /
• 2013

The current distribution passing through grounding electrode is required for calculating an impedance of grounding electrode using the electromagnetic field model. In this paper the numerical calculation for currents passing through a grounding electrode as a function of frequency was given. The proposed approach is based on the wire antenna model(AM) in the frequency domain. The Pocklington's equation driven from the wire antenna theory was numerically calculated by the Galerkin's method. The triangle function was applied to both the basis function and the weighting function. The current distribution of a horizontal ground electrode was simulated in MATLAB. Also these results were compared with the data obtained from the CDEGS HIFREQ calculation.

• 한국정밀공학회지
• /
• 제21권1호
• /
• pp.189-196
• /
• 2004

A modified generalized particle algorithm, MGPA, was suggested to improve the calculation efficiency of standard SPH Method in numerical analysis of high speed impact behavior. MGPA had a new weight function to reduce computation time. The efficiency of this method was proven through calculation for the sample problems of one dimensional rod impact problem and two dimensional plate impact problem. The MGPA method reduced the calculation error and stress oscillation near the boundaries. The validity of this approach was shown by the comparison with ABAQUS results in two dimensional plate impact problem.