• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.5
• /
• pp.441-449
• /
• 2015

In this paper we perform a linearized buckling analysis using the Kirchhoff plate theory and the von Karman nonlinear strain-displacement relation. Design sensitivity analysis(DSA) expressions for plane elasticity and buckling problems are derived with respect to Young's modulus and thickness. Using the design sensitivity, we can formulate the topology optimization method for minimizing the compliance and maximizing eigenvalues. We develop a topology optimization method applicable to plate buckling problems using the prestress for buckling analysis. Since the prestress is needed to assemble the stress matrix for buckling problem using the von Karman nonlinear strain, we introduced out-of-plane motion. The design variables are parameterized into normalized bulk material densities. The objective functions are the minimum compliance and the maximum eigenvalues and the constraint is the allowable volume. Through several numerical examples, the developed DSA method is verified to yield very accurate sensitivity results compared with the finite difference ones and the topology optimization yields physically meaningful results.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.11-23
• /
• 1997

This paper contains a new approach to blade section design method for marine propellers. The hydrodynamic characteristics of 2-D section are highly influenced by its geometrical parameters i.e., thickness and camber distributions and leading edge radius etc. To consider fully turbulent flow field near 2-D section. the finite volume method with k-${\varepsilon}$ turbulent model which solve Reynolds time averaged Navier-Stokes(RANS) equation is applied. In this study, O-type grid system that can provide many calculation points on blade surface is used. The results were compared with those of the experiment of NACA0012 to confirm the accuracy of the developed codes. The goal of this study is the development of a blade section with high efficiency and low drag. To achieve this, we carried out the tests of lift, drag and cavitation characteristics in cavitation tunnel. The results of experiment were compared with numerical results in order to validate the proposed blades design method. By comparing the numerical results with the experiments, we found that the new blade section, KH28 allows superior performance in efficiency and cavitation avoidance characteristics. We further investigated the blade section design method and an application study of this section, KH28 to apply to the marine propeller. In order to improve the accuracy of numerical results on prediction of lift and drag, we conclude here that the 2-layer boundary model must be used.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.3
• /
• pp.252-264
• /
• 2018

The calculation of steady-state cavitating flows around Supercavitating Underwater Bodies (SUB's), which consist of a circular disk head (cavitator), a conical fore-body, a cylindrical middle-body and either a boat-tail or a flare-tail, are carried out. To calculate the axisymmetric cavitating flow, used is a commercial computational fluid dynamics code based on the finite volume method, Fluent. From the analysis of numerical results, the cavity and drag, affected by the fore-body and tail of the SUB's, are investigated. Firstly, the effect of the fore-body shape is investigated with the same disk cavitator and a cylindrical rear-body of fixed diameter. Then with the same cavitator and a fixed fore-body, the effect of the rear-body shape is investigated. Before the cavity generated by the cavitator covers the slant of fore-bodies sufficiently, the larger the cone angle of the fore-body(i.e., the shorter the slant length), the larger the drag and the slower the development of cavity. After the cavity covers the fore-body completely so that the pressure drag component of the body is vanished, the characteristics of drag-velocity curves are identical. Also, as the tail angle is bigger, the cavity generated by the cavitator is suppressed further and the drag becomes larger. The peak of the drag appears for the flare-tail, i.e., when the tail angle is positive(+). On the contrary, the trough of the drag appears for the boat-tail, i.e., when the tail angle is negative(-). When the tail angle is 5 degrees, the peak of the drag appears at the body speed of 80m/s and the value of the drag is 43% larger than that at the design speed of 100m/s. When the tail angle is -5 degrees, the trough of the total drag appears at 75m/s and that drag is 30% smaller than that of the cavitator, which means the rest of the body has a negative drag.

• Journal of the Korean Vacuum Society
• /
• v.19 no.4
• /
• pp.249-255
• /
• 2010

Methods of the characteristics evaluation of turbo-molecular pumps (TMP) are well-defined in the international measurement standards such as ISO, PNEUROP, DIN, JIS, and AVS. The Vacuum Center in the Korea Research Institute of Standards and Science (KRISS) has recently designed, constructed, and established the integrated characteristics evaluation system of TMPs based on the international documents by continuously pursuing and acquiring the reliable international credibility through measurement perfection. The measurement of TMP pumping speed is normally performed with the throughput and orifice methods dependent on the mass flow regions. However, in the UHV range of the molecular flow region, the high uncertainties of the gauges, mass flow rates, and conductance are too critical to precisely accumulate reliable data. In order to solve the uncertainty problems of pumping speeds in the UHV range, we introduced a SRG with 1% accuracy and a constant volume flow meter (CVFM) to measure the finite mass flow rates down to $10^{-1}$ Pa-L/s with 3% uncertainty for the throughput method. In this way we have performed the measurement of pumping speed down to $10^{-4}$ Pa with an uncertainty of less than 6% for a 1000 L/s TMP. In this article we suggest that the CVFM has an ability to measure the conductance of the orifice experimentally with flowing the known mass through the orifice chambers, so that we may overcome the discontinuity problem encountering during introducing two measurement methods in one pumping speed evaluation sequence.

• Journal of Korea Water Resources Association
• /
• v.45 no.1
• /
• pp.39-52
• /
• 2012

Surface-subsurface interactions are an intrinsic component of the hydrologic response within a watershed. In general, these interactions are considered to be one of the most difficult areas of the discipline, particularly for the modeler who intends simulate the dynamic relations between these two major domains of the hydrological cycle. In essence, one major complexity is the spatial and temporal variations in the dynamically interacting system behavior. The proper simulation of these variations requires the need for providing an appropriate coupling mechanism between the surface and subsurface components of the system. In this study, an approach for modelling surface-subsurface flow and transport in a fully intergrated way is presented. The model uses the 2-dimensional diffusion wave equation for sheet surface water flow, and the Boussinesq equation with the Darcy's law and Dupuit-Forchheimer's assumption for variably saturated subsurface water flow. The coupled system of equations governing surface and subsurface flows is discretized using the finite volume method with central differencing in space and the Crank-Nicolson method in time. The interactions between surface and subsurface flows are considered mass balance based on the continuity conditions of pressure head and exchange flux. The major module consists of four sub-module (SUBFA, SFA, IA and NS module) is developed.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.12
• /
• pp.5-15
• /
• 2012

In this study, a series of CFD (Computational Fluid Dynamics) numerical analyses were performed in order to evaluate the thermal performance of six full-scale closed-loop vertical ground heat exchangers constructed in a test bed located in Wonju. The circulation HDPE pipe, borehole and surrounding ground formation were modeled using FLUENT, a finite-volume method (FVM) program, for analyzing the heat transfer process of the system. Two user-defined functions (UDFs) accounting for the difference in the temperatures of the circulating inflow and outflow fluid and the variation of the surrounding ground temperature with depth were adopted in the FLUENT model. The relevant thermal properties of materials measured in laboratory were used in the numerical analyses to compare the thermal efficiency of various types of the heat exchangers installed in the test bed. The simulation results provide a verification for the in-situ thermal response test (TRT) data. The CFD numerical back-analysis with the ground thermal conductivity of 4 W/mK yielded better agreement with the in-situ thermal response tests than with the ground thermal conductivity of 3 W/mK.

• Journal of the Korea Society for Simulation
• /
• v.27 no.4
• /
• pp.19-26
• /
• 2018

The aerial launching UAV(Unmanned Aerial Vehicle) mainly uses a set of folding tandem wings to maximize flight performance and minimize the space required for mounting in a mothership. This folding tandem wing has a unique aerodynamic problem that is different from the general type of fixed wing aircraft, such as the rear wing interference problem caused by the wing of the front wing wake and vortex, and the imbalance of the pivot moment applied to the front and rear wings when the wing is deployed. In this paper, we have modeled and simulated various cases through computational fluid dynamics based on the finite volume method and analyzed various aerodynamic phenomena of the tandem wing type aircraft. We find that the front wing shall be installed higher than the rear for minimizing the wake influence and the rear wing can be deployed faster than the front because of the pivot moment due to aerodynamic forces. Also, considering the pivot moment due to aerodynamic force, the rear wing can be deployed much faster than the front wing. Therefore, it is necessary to consider it when developing the wing deploy mechanism.

• Journal of Korea Society of Industrial Information Systems
• /
• v.24 no.2
• /
• pp.33-40
• /
• 2019

Many companies have tried to develop wind power generators with a larger capacity, smaller size and lighter weight. High temperature superconducting (HTS) generators are more suitable for wind power systems because they can reduce volume and weight compared with conventional generators. However, the HTS generator has problems such as huge vacuum vessel and the difficulty of repairing the HTS field coils. These problems can be overcome through the modularization of the HTS field coil. The HTS module coil require a current leads (CLs) for deliver DC current, which causes a large heat transfer load. Therefore, CLs should be designed optimally for reducing the conduction and Joule heat loads. This paper deals with a structural design and thermal analysis of a module coil for a 750 kW-class HTS generator. The conduction and radiation heat loads of the module coils were analysed using a 3D finite element method program. As a result, the total thermal load was less than the cooling capacity of the cryo-cooler. The design results can be effectively utilized to develop a superconducting generator for wind power generation systems.

• Tunnel and Underground Space
• /
• v.30 no.6
• /
• pp.509-518
• /
• 2020

This study constructed a new simulation platform, including mesh generation process, numerical simulation, and post-processing for results analysis based on exploration data to perform real-scale numerical analysis considering the actual geological structure efficiently. To build the simulation platform, we applied for open-source programs. The source code is open to be available for code modification according to the researcher's needs and compatibility with various numerical simulation programs. First, a three-dimensional model(3D) is acquired based on the exploration data obtained using a drone. Then, the domain's mesh density was adjusted to an interpretable level using Blender, the free and open-source 3D creation suite. The next step is to create a 3D numerical model by creating a tetrahedral volume mesh inside the domain using Gmsh, a finite element mesh generation program. To use the mesh information obtained through Gmsh in a numerical simulation program, a converting process to conform to the program's mesh creation protocol is required. We applied a Python code for the procedure. After we completed the stability analysis, we have created various visualization of the study using ParaView, another open-source visualization and data analysis program. We successfully performed a preliminary stability analysis on the full-scale Dokdo model based on drone-acquired data to confirm the usefulness of the proposed platform. The proposed simulation platform in this study can be of various analysis processes in future research.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.26 no.7
• /
• pp.922-930
• /
• 2020

Marine accidents involving fishing boats, caused by a loss of stability, have been increasing over the last decade. One of the main reasons for these accidents is a sudden wind attacks. In this regard, the wind loads acting on the ship hull need to be estimated accurately for safety assessments of the motion and maneuverability of the ship. Therefore, this study aims to develop a computational model for the inlet boundary condition and to numerically estimate the wind load acting on a fishing boat. In particular, wind loads acting on a fishing boat at the wind speed profile boundary condition were compared with the numerical results obtained under uniform wind speed. The wind loads were estimated at intervals of 15° over the range of 0° to 180°, and i.e., a total of 13 cases. Furthermore, a numerical mesh model was developed based on the results of the mesh dependency test. The numerical analysis was performed using the RANS-based commercial solver STAR-CCM+ (ver. 13.06) with the k-ω turbulent model in the steady state. The wind loads for surge, sway, and heave motions were reduced by 39.5 %, 41.6 %, and 46.1 % and roll, pitch, and yaw motions were 48.2 %, 50.6 %, and 36.5 %, respectively, as compared with the values under uniform wind speed. It was confirmed that the developed inlet boundary condition describing the wind speed gradient with respect to height features higher accuracy than the boundary condition of uniform wind speed. The insights obtained in this study can be useful for the development of a numerical computation method for ships.