• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 2019.05a
• /
• pp.185-188
• /
• 2019

In this paper, the towing stability of the LNG bunker barge is estimated. Currently, LNG bunkering barge is being developed as an infrastructure for the bunkering of LNG (Liquefied Natural Gas), an eco-friendly energy source. Since the LNG bunker barge are in the form of towed ship connected to the tow line, the towing stability of the LNG bunker barge is very important for the safety of not only the LNG bunker barge but also the surrounding sailing vessels. The numerical code for towing simulation was developed to estimate the towing stability of the LNG bunker barge at the initial design stage. The MMG(Manoeuvring Mathematical Group) model was applied to the equations of motion and the empirical formula was applied to the maneuvering coefficients so that they could be used in the initial design stage. To validity of the developed numerical code, it was compared with published calculation and model test results. Towing simulations were carried out according to with and without stern skeg of the LNG bunker barge using the developed numerical code. Through the results of the simulations, the appropriateness of the stern skeg area designed was confirmed.

• International Journal of Air-Conditioning and Refrigeration
• /
• v.10 no.3
• /
• pp.117-128
• /
• 2002

This paper describes the smoke filling process of a fire field model based on a self-deve-loped SMEP (Smoke Movement Estimating Program) code to the simulation of fire induced flows in the two types of atrium space containing a ceiling heat flux. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k-$\varepsilon$ turbulence model with buoyancy production term. Also it solves the radiation equation using the discrete ordinates method. Compressibility is assumed and the perfect gas law is used. Comparison of the calculated upper-layer average tempera-ture and smoke layer clear height with the zone models has shown reasonable agreement. The zone models used are the CFAST and the NBTC one-room. For atrium fires with ceiling glass the ceiling heat flux by solar heat causes a high smoke temperature near the ceiling. However, it has no effect on the smoke movement such as the smoke layer clear heights that are important in fire safety. In conclusion, the smoke layer clear heights that are important in evacuation activity except the early of a fire were not as sensitive as the smoke layer tem-perature to the nature of ceiling heat flux condition. Thus, a fire sensor in atrium with ceiling glass has to consider these phenomena.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.3247-3252
• /
• 2007

Endwall losses contribute significantly to the overall losses in modern turbomachinery, especially when aerodynamic airfoil load and pressure ratio are increased. Hence, reducing the extend and intensity of the secondary flow structures helps to enhance overall efficiency. From the large range of viable approaches, a promising combination positioning and height of endwall contouring was chosen. The objective of this study is to document the three-dimensional flow in a turbine cascade in terms of streamwise vorticity, total pressure loss distribution and static pressure distribution on the endwall and blade surface and to propose an appropriate positioning and height of the endwall contouring which show best secondary, overall loss reduction among the simulated endwall. The flow through the gas turbine were numerically analyzed using three dimensional Navier-Stroke equations with a commercial CFD code ANSYS CFX-10. The result shows that the overall loss is reduced near the flat endwall rather than contoured endwall, and the case of contoured endwall installed at 30% from leading edge with height of 25% for span showed best performance.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.5
• /
• pp.403-413
• /
• 2010

Transition location of boundary layers over airfoils is predicted by using PSE(Parabolized Stability Equations) and $e^N$-method. Growth rates of disturbances are obtained from the PSE analysis and the N-factor curves are calculated by integrating the growth rates. The computational code developed in the present study is validated by comparing the computed results with the well known data for the cases of flat plate boundary layers and airfoils. Predictions of transition location are made for the boundary layers over NACA0012, NLF(1)-0414F, and NLF(1)-0416 airfoil. Predicted transition locations are found to be in good agreement with the experimental data.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.17-23
• /
• 2011

The turbulent flow and vortex shedding phenomena around pantograph panhead of high speed train were investigated and compared with available experimental data and other simulations. The pantograph head was simplified to be a square-cross-section pillar and assumed to be no interference with other bodies. The Reynolds number (Re) was 22,000. The LES(large eddy simulation) of FDS code was applied to solve the momentum equations and the Wener-Wengle wall model was employed to solve the near wall turbulent flow. Smagorinsky model($C_s$=0.2) was used as SGS(subgrid scale) model. The total grid numbers were about 9 millions and the analyzed domain was divided into 12 multi blocks which were communicated with each other by MPI. The time-averaged mainstream flows were calculated and well compared with experimental data. The phased-averaged quantities had also a good agreement with experimental data. The near-wall turbulence should be carefully treated by wall function or direct resolution to get successful application of LES methods.

• Geomechanics and Engineering
• /
• v.10 no.4
• /
• pp.423-436
• /
• 2016

Currently, there is a great interest in the coupling between multiphase fluid flow and geomechanical effects in hydrocarbon reservoirs and surrounding rocks. The ideal solution for this coupled problem is to introduce the geomechanical effects through the stress analysis solution and implement an algorithm, which assures that the equations governing the flow and stress analyses are obeyed in each time step. This paper deals with the implementation of a program (FORTRAN90 interface code), which was developed to couple conventional reservoir (ECLIPSE) and geomechanical (ABAQUS) simulators, using a partial coupling algorithm. The explicit coupled hydro-mechanical behavior of Iranian field during depletion and $CO_2$ injection is studied using the soils consolidation procedure available in ABAQUS. Time dependent reservoir pressure fields obtained from three dimensional compositional reservoir models were transferred into finite element reservoir geomechanical models in ABAQUS as multi-phase flow in deforming reservoirs cannot be performed within ABAQUS. The FEM analysis of the reservoir showed no sign of plastic strain under production and $CO_2$ injection scenarios in any part of the reservoir and the stress paths do not show a critical behavior.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.297-306
• /
• 2014

The paper presents the numerical analysis of wave run-up over rubble-mound breakwaters covered by antifer units using a technique integrating Computer-Aided Design (CAD) and Computational Fluid Dynamics (CFD) software. Direct application of Navier-Stokes equations within armour blocks, is used to provide a more reliable approach to simulate wave run-up over breakwaters. A well-tested Reynolds-averaged Navier-Stokes (RANS) Volume of Fluid (VOF) code (Flow-3D) was adopted for CFD computations. The computed results were compared with experimental data to check the validity of the model. Numerical results showed that the direct three dimensional (3D) simulation method can deliver accurate results for wave run-up over rubble mound breakwaters. The results showed that the placement pattern of antifer units had a great impact on values of wave run-up so that by changing the placement pattern from regular to double pyramid can reduce the wave run-up by approximately 30%. Analysis was done to investigate the influences of surface roughness, energy dissipation in the pores of the armour layer and reduced wave run-up due to inflow into the armour and stone layer.

• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.81-87
• /
• 2015

For Underwater vehicles, Unwanted roll excursions are inevitable as they are caused by induced propeller torque, disturbances, and banking motion during turns. To estimate the manoeuvring performance of underwater vehicle, it is necessary to obtain the roll coefficient of body. This paper was covered estimation of roll coefficient of underwater vehicle using STAR-CCM+, commercial CFD(Computational Fluid Dynamics) code. The RANS equations for incompressible fluid flows was solved numerically by using a finite volume method. An MRF(Moving Reference Frame) Method was Also adopted for rotations of body. For the validation, the flow around a DARPA SUBOFF bare hull model was simulated and good agreement with experiments was obtained. And Pure roll coefficients were calculated and campared with the experimental data which were presented by Seoul National University. Finally, an underwater vehicle model with propeller was simulated and analyzed for estimation of roll coefficient variation caused by induced propeller torque.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.22 no.5
• /
• pp.852-858
• /
• 2013

In numerical control (NC) machining, a machining error in equipment generally occurs for a variety of reasons. If there is a change in direction in the NC code, the characteristics of the automatic acceleration or deceleration function cause an overlap of each axis of the acceleration and deceleration zones, which in turn causes a shift in the actual processing path. Many studies have been conducted for error calibration of the edge as caused by automatic acceleration or deceleration in NC machining. This paper describes a geometric interpretation of the shape and processing characteristics of the operating NC device. The paper then describes a way to determine a feedrate that achieves the desired tolerance by using linear and parabolic profiles. Experiments were conducted by the validate equations using a three-axis NC machine. The results show that the machining errors were smaller than the machine resolution. The results also clearly demonstrate that the NC machine with the developed system can successfully predict machining errors induced with a change in direction.

• Computers and Concrete
• /
• v.19 no.6
• /
• pp.631-639
• /
• 2017

This study experimentally investigated the flexural capacity of a concrete beam reinforced with a newly developed GFRP bar that overcomes the lower modulus of elasticity and bond strength compared to a steel bar. The GFRP bar was fabricated by thermosetting a braided pultrusion process to form the outer fiber ribs. The mechanical properties of the modulus of elasticity and bond strength were enhanced compared with those of commercial GFRP bars. In the four-point bending test results, all specimens failed according to the intended failure mode due to flexural design in compliance with ACI 440.1R-15. The effects of the reinforcement ratio and concrete compressive strength were investigated. Equations from the code were used to predict the deflection, and they overestimated the deflection compared with the experimental results. A modified model using two coefficients was developed to provide much better predictive ability, even when the effective moment of inertia was less than the theoretical $I_{cr}$. The deformability of the test beams satisfied the specified value of 4.0 in compliance with CSA S6-10. A modified effective moment of inertia with two correction factors was proposed and it could provide much better predictability in prediction even at the effective moment of inertia less than that of theoretical cracked moment of inertia.