• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• 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.

• Journal of Wetlands Research
• /
• v.23 no.2
• /
• pp.107-115
• /
• 2021

In order to analyze wave propagation, tidal current, and sediment transport in the vicinity of the Gaeya open channel, it was classified into before(CASE1W) and after(CASE2W) installation of various artificial structures, and the calculation results for CASE1W and CASE2W were compared. For wave propagation, the results of incident and reflected waves were derived using the SWAN numerical model, and the tidal current velocity results were derived using the FLOW2DH numerical model for tidal current. The results of the SWAN numerical model and the FLOW2DH numerical model became the input conditions for the SEDTRAN numerical model that predicts sediment transport, and the maximum bed shear stress and suspended sediment concentration distribution near the Gaeya open channel were calculated through the SEDTRAN numerical model. As a result of the calculation of the SWAN numerical model, the wave height of CASE2W was increased by 40~50 % compared to CASE1W because the incident wave was diffracted and superimposed and the reflected wave was generated by about 7 km long northen jetty. As a result of the calculation of the FLOW2DH numerical model, According to the northen breakwater, the northen jetty and Geumrando, CASE2W was calculated 10~30 % faster than CASE1W in the tidal current of the Gaeya open channel. As a result of the calculation of the SEDTRAN numerical model, the section where the maximum bed shear stress is 1.0 N/m² or more and the suspended concentration is 80mg/L or more was widely distributed in the Gaeya open channel from the marine environment by the complex wave field(incident wave, reflected wave and tidal wave) and the installation of various artificial structures. it is believed that a sedimentation phenomenon occurred in the Gaeya open channel.

• KIEAE Journal
• /
• v.16 no.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.

• Progress in Superconductivity and Cryogenics
• /
• v.17 no.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.

• Explosives and Blasting
• /
• v.36 no.3
• /
• pp.10-18
• /
• 2018

The dynamic behaviour of the rock mass under the dynamic load is different from the static application of the maximum load of the same size. An experimental approach to investigating rock behavior under dynamic loads is more difficult than that under static conditions in control of dynamic loads, measurement and analysis of the results. Numerical methods are less constrained by performing the experiments numerically, rather than experimental ones, so they can be very powerful analytical tool at the design stage. However, even if the algorithms of the analysis method are appropriate, careful analysis is required because the calculation results may vary largely depending on input data and boundary conditions. In this paper, when investigating the behavior of rock structures under dynamic load numerically, the effects of boundary conditions, dynamic load and calculation time step, and dynamic load characteristics on the calculation results were reviewed to provide guidance on setting up boundary conditions and calculation time step related to dynamic analysis.

• Structural Engineering and Mechanics
• /
• v.64 no.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.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.27 no.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.

• Journal of the Korean Society for Precision Engineering
• /
• v.21 no.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.

• Solar Energy
• /
• v.9 no.3
• /
• pp.64-72
• /
• 1989

In this study, the computer simulation for the heat transfer in pulse combustion water heater is performed. The attention is focused to the effects of the installation of corebuster in the flue tube on heat transfer. The energy equations are established for both wall and gas side in the combustion chamber, flue way, exhaust chamber and muffler, and the numerical calculation is executed. Zone method takes longer computer calculation time compared with semi-zone method. Semi-zone method is chosen for numerical calculation. As a result of this study, it is found that the installation of the core buster in flue tube increases total heat transfer. It is also found that the total heat transfer is increased with the increasing of the ratio of the cross section area of corebuster to that of the flue tube. However, the heat transfer effect is negligible for the area ratio above 0.5.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.8
• /
• pp.1206-1212
• /
• 2001

The kinematic analysis of Stewart platform manipulator(SPM) is carried out in order to reduce the calculation time for its forward kinematic solution when the iterative numerical method is employed. The kinematic equations for three substructures of the 6-3 SPM are newly derived by introducing Denavit-Hartenberg link parameters and using kinematic constraints associated with the SPM and substructure kinematics. It is shown that the forward kinematics can be easily solved from three nonlinear equations with three unknown variables only, leading to a great reduction in calculation time.