• Journal of the Korean Society of Propulsion Engineers
• /
• v.9 no.2
• /
• pp.97-104
• /
• 2005

A simple but efficient grid generation technique by using the modified compressible form of stream function has been formulated. Transformation of a physical plane to a streamline plane, the Von Mises Transformation, has been widely used to solve the differential equations governing flow phenomena, however, limitation arises in low velocity region of boundary layer, mixing layer and wake region where the relatively large grid spacing is inevitable. Modified Von Mises Transformation with simple mathematical adjustment for the stream function is suggested and applied to solve the confined coaxial turbulent jet mixing with simple $\kappa-\epsilon$ turbulence model. Comparison with several experimental data of axial mean velocity, turbulent kinetic energy, and Reynolds shear stress distribution shows quite good agreement in the mixing layer except in the centerline where the turbulent kinetic energy distributions were somewhat under estimated. This formulation is strongly suggested to be utilized specially for free turbulent mixing layers in axisymmetric flow conditions such as the investigation of mixing behavior, jet noise production and reduction for Turbofan engines.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.2 no.1
• /
• pp.1-12
• /
• 1998

Industrial ejector system is a facility to transport, to compress or to pump out a low pressure secondary flow by using a high pressure primary flow. An advantage of the ejector system is in its geometrical simplicity, not having any moving part, compared with other fluid machinery. Most of the previous works have been performed experimentally and analytically. The obtained data. are too insufficient to improve our current understanding on the detailed flow field inside the ejector. In order to provide more comprehensive data on this ejector flow field, two-dimensional computations using Reynolds-averaged Navier-Stokes equations were performed for a very wide range of operating pressure ratio of the supersonic ejector with a secondary throat. The current results showed that the supersonic ejector system has an optimum pressure ratio for the secondary flow total pressure to be minimized. The numerical results clearly revealed the shock system, shock/boundary layer interaction, and secondary flow entrainment inside the supersonic ejector.

• Journal of Ship and Ocean Technology
• /
• v.10 no.4
• /
• pp.1-11
• /
• 2006

The finite volume based multi-block RANS code, WAVIS developed at MOERI is applied to the numerical self-propulsion test. WAVIS uses the cell-centered finite volume method for discretization of the governing equations. The realizable $k-{\epsilon}$ turbulence model with a wall function is employed for the turbulence closure. The free surface is captured with the two-phase level set method and body forces are used to model the effects of a propeller without resolving the detail blade flow. The propeller forces are obtained using an unsteady lifting surface method based on potential flow theory. The numerical procedure followed the self-propulsion model experiment based on the 1978 ITTC performance prediction method. The self-propulsion point is obtained iteratively through balancing the propeller thrust, the ship hull resistance and towing force that is correction for Reynolds number difference between the model and full scale. The unsteady lifting surface code is also iterated until the propeller induced velocity is converged in order to obtain the propeller force. The self-propulsion characteristics such as thrust deduction, wake fraction, propeller efficiency, and hull efficiency are compared with the experimental data of the practical container ship. The present paper shows that hybrid RANS and potential flow based numerical method is promising to predict the self-propulsion parameters of practical ships as a useful tool for the hull form and propeller design.

• Journal of the Korean Solar Energy Society
• /
• v.35 no.1
• /
• pp.35-43
• /
• 2015

A numerical simulation on the wake flow of a wind turbine which is a scaled version of a multi-megawatt wind turbine has been performed. Two different inlet conditions of averaged wind speed including one below and one above the rated wind speed were used in the simulation. Steady-state pitch angles of the blade associated with the two averaged wind speeds were imposed for the simulation. The steady state analysis based on the Reynolds averaged Navier-Stokes equations with the method of frame motion were used for the simulation to find the torque of the rotor and the wake field behind the wind turbine. The simulation results were compared with the results obtained from the wind tunnel testing. From comparisons, it was found that the simulation results on the turbine power are pretty close to the experimental values. Also, the wake results were relatively close to the experimental results but there existed some discrepancy in the shape of velocity deficit. The reason for the discrepancy is considered due to the steady state solution with the frame motion method used in the simulation. However, the method is considered useful for solutions with much reduced calculation time and reasonably good accuracy compared to the transient analysis.

• Tribology and Lubricants
• /
• v.37 no.4
• /
• pp.129-135
• /
• 2021

An air-bearing stage uses externally pressurized air as the lubricant between the stage and the rail. The supporting force generated by the supplied air makes the stage rise and move smoothly with extremely low friction. Mechanical contacts rarely happen, the bearing surfaces do not produce wear particles, and dust is not generated. It also has the advantage of having low energy loss and high precision. Because of its advantages, an air-bearing stage is used in several types of machines that require high precision. In this article, the effect of the pocket depth on the hammering phenomena of the air bearing is studied. An analysis program is developed to calculate the dynamic behavior of the stage by solving the Reynolds equation between the stage and the guideway and the equations of motion on the stage. The acceleration, constant movement, and deceleration are applied to the stage. The stage is modeled as a five-degree-of-freedom system. In the course of the dynamic behavior, the hammering phenomena occur under some special conditions. The deeper the pocket, the more unstable the behavior of the stage, and air hammering occurs when it exceeds a certain depth. In addition, the higher the supply pressure, the more unstable the behavior of the stage. However, hammering occurs even with a shallow pocket depth. Other conditions that affect the hammering phenomena are calculated and discussed.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.6
• /
• pp.474-481
• /
• 2018

This paper is to compare by numerical analysis the flow characteristics and propulsion performance of stern with the shape change of K-duct, a pre-swirl duct developed by Korea Research Institute of Ships & Ocean Engineering (KRISO). First, the characteristics of the propeller and the resistance and self-propulsion before and after the attachment of the K-duct to the ship were verified and the validity of the calculation method was confirmed by comparing this result with the model test results. After that, resistance and self-propulsion calculations were performed by the same numerical method when the K-duct was changed into five different shapes. The efficiency of the other five cases was compared using the delivery horsepower in the model scale and the flow characteristics of the stern were analyzed as the velocity and pressure distributions in the area between the duct end and the propeller plane. For the computation, STAR-CCM +, a general-purpose flow analysis program, was used and the Reynolds Averaged Navier-Stokes (RANS) equations were applied. Rigid Body Motion (RBM) method was used for the propeller rotating motion and SST $k-{\omega}$ turbulence model was applied for the turbulence model. As a result, the tangential velocity of the propeller inflow changed according to the position angle change of the stator, and the pressure of the propeller hub and the cap changes. This regulated the propeller hub vortex. It was confirmed that the vortex of the portion where the fixed blade and the duct meet was reduced by blunt change.

• Journal of Convergence for Information Technology
• /
• v.9 no.8
• /
• pp.155-163
• /
• 2019

Laminar mixed convection of a nanofluid consists of water and $Al_2O_3$ in a horizontal circular tube has been studied numerically. Two-phase mixture model has been used to investigate hydrodynamic and thermal behaviors of the nanofluid with variables physical properties. Three dimensional Navier-Stokes, energy and volume fraction equations have been discretized using the finite volume method. The Brownian motions of nanoparticles have been considered to determine the thermal conductivity and dynamic viscosity of $Al_2O_3$-Water nanofluid, which depend on temperature. The calculated results show good agreement with the previous numerical data. Results show that in a given Reynolds number (Re), increasing solid nanoparticles volume fraction and Richardson number (Ri) increases the convective heat transfer coefficient and wall shear stress.

• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.4
• /
• pp.373-382
• /
• 2019

This paper shows numerical results for the estimation of the propulsor efficiency of Pre-Swirl Duct for 176k bulk carrier as well as its design method. Reynolds averaged Navier-Stokes equations have been solved and the k-epsilon model applied for the turbulent closure. The propeller rotating motion is determined using a sliding mesh technique. The design process is divided into each part of Pre-Swirl Duct, duct and Pre-Swirl Stator. The design of duct was performed first because it is located further upstream than Pre-Swirl Stator. The distribution of velocity through the duct was analyzed and applied for the design of Pre-Swirl Stator. The design variables of duct include duct angle, diameter, and chord length. Diameter, chord length, equivalent angle are considered when designing the Pre-Swirl Stator. Furthermore, a variable pitch angle stator is applied for the final model of Pre-Swirl Duct. The largest reduction rate of the delivered power in model scale is 7.6%. Streamlines, axial and tangential velocities under the condition that the Pre-Swirl Duct is installed were reviewed to verify its performance.

• Advances in aircraft and spacecraft science
• /
• v.6 no.3
• /
• pp.239-256
• /
• 2019

This study develops a dual purpose: i) evaluating the effects of two different Mars atmosphere models (NASA Glenn and GRAM-2001) on aerodynamics of a capsule (Pathfinder) entering the Mars atmosphere, ii) verifying the feasibility of evaluating the ambient density and pressure by means of the methods by McLaughlin and Cassanto, respectively and therefore to re-build the values provided by the models. The method by McLaughlin relies on the evaluation of the capsule drag coefficient, the method by Cassanto relies on the measurement of pressure at a point on the capsule surface in aerodynamic shadow. The study has been carried out computationally by means of: i) a code integrating the equations of dynamics of the capsule for the computation of the entry trajectory, ii) a DSMC code for the solution of the flow field around the capsule in the altitude interval 50-100 km. The models show consistent differences at altitudes higher than about 40 km. It seems that the GRAM-2001 model is more reliable than the NASA Glenn model. In fact, the NASA Glenn model produces, at high altitude, temperatures that seem to be too low compared with those from the GRAM-2001 model and correspondingly very different aerodynamic conditions in terms of Mach, Reynolds and Knudsen numbers. This produces pretty different capsule drag coefficients by the two models as well as pressure on its surface, making not feasible neither the method by McLaughlin nor that by Cassanto, until a single, reliable model of the Mars atmosphere is not established. The present study verified that the implementation of the Cassanto method in Mars atmosphere should rely (such as it is currently) on pressure obtained experimentally in ground facilities.

• Tribology and Lubricants
• /
• v.37 no.5
• /
• pp.172-178
• /
• 2021

An injector is a mechanical device present inside the engine. Its main function is to supply an appropriate volume of fuel into the combustion chamber, which is directly related to the overall engine efficiency of a car. During the operation of an injector, a magnetic force lifts the parts of the injector from closed position to open position which generates a horizontal force on the needle. The horizontal force acts on a different position from that of the center of mass of the needle. Therefore, this causes eccentricity in the needle and the generation of a tilting motion during the lifting operation which can result in wear. However, appropriate selection of clearances for these parts can prevent wear. In this study, lubrication analysis is conducted to determine the optimum clearance of parts with sliding motion inside the injector. The height functions are derived considering the dynamic behavior and relative velocity of the parts. Using the derived height function, the pressure profiles are calculated for the lubricated surfaces from the Reynolds' equation. Subsequently, the fluid reaction forces are calculated. The equations of motions are applied to the fluid reaction forces and external forces are solved to calculate the minimum film thickness between each part with variation in the clearances. Finally, the optimum clearances are determined. The effect of the clearances on the behavior of the moving parts is presented and discussed.