• International Journal of Naval Architecture and Ocean Engineering
• /
• v.2 no.3
• /
• pp.139-145
• /
• 2010

In the present study, we conducted resistance test, propeller open water test and self-propulsion test for a ship's resistance and propulsion performance, using computational fluid dynamics techniques, where a Reynolds-averaged Navier-Stokes equations solver was employed. For convenience of mesh generation, unstructured meshes were used in the bow and stern region of a ship, where the hull shape is formed of delicate curved surfaces. On the other hand, structured meshes were generated for the middle part of the hull and the rest of the domain, i.e., the region of relatively simple geometry. To facilitate the rotating propeller for propeller open water test and self-propulsion test, a sliding mesh technique was adopted. Free-surface effects were included by employing the volume of fluid method for multi-phase flows. The computational results were validated by comparing with the existing experimental data.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.1
• /
• pp.22-26
• /
• 2011

The resistance performances of a small leisure yacht with two different keels, center and double, respectively, are investigated using a model test at circulating water channel and CFD analysis. Flow patterns around the keels are observed using a tuft test to make clear the relation between the resistance performance and the flow characteristics around the keels. The results show that the keel does not affect free surface flows and that the double keel yacht has better performance compared to a single keel yacht in oblique condition from the resistance point of view.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.12
• /
• pp.3386-3394
• /
• 1994

An experimental study has been performed on the flow over a rotating disk, where the diameter of the disk is 500 mm and the maximum vertical deviation of the upper surface is $50 \mu{m}$ for the whole range of the angular velocity up to 3400 rpm. The flow visualization experiment for the wall jet flow induced by impinging circular jet is carried out using schlieren system and measurements are made by 3-hole and 5-hole pitot tubes. Schlieren photographs show that as the rotating speed increases the wall jet flow becomes more stable and the size of the largest eddies becomes smaller. Measurements for impinging jet flow on the stationary disk verify the accuracy of the present experiment, and those for free rotating disk flow display the existence of transition region from laminar to turbulent flows. Measurements for impinging jet flow on the rotating disk exhibit the interaction between the wall jet and the viscous pumping effect, which explains the decay in size of turbulent eddies illustrated by the schlieren photographs.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.11
• /
• pp.1538-1551
• /
• 1997

Kinetic energy conservation for fixed and moving grids is examined in time-accurate finite element computation of fully unsteady inviscid flows. As numerical algorithms, fractional step method (FSM) and modified SIMPLE are used. To simulate the flow in moving grid system, arbitrary Lagrangian-Eulerian (ALE) method is adopted. In the present study, the energy conserving time integration rule for finite element algorithm is proposed and discussed schematically. It is shown that the discretization by Crank-Nicolson in time and Galerkin (central difference) in space must be used to ensure energy conservation. The developed code has been tested for a standing vortex in fixed or moving grid system, sloshing in a tank and propagation of a solitary wave, and has been shown to be a completely energy conserving algorithm.

• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.85-92
• /
• 1996

Numerical simulations are made for the three-dimensional flow around a wing in ground effect craft haying the complex geometry. A numerical tool is developed for the primary design of hull and wing shape of practical Wing-In-Ground effect(WIG) stop. The finite-difference method is utilized to descretize the governing equations and pressure field is obtained by using Marker-And-Cell(MAC) method. The air and water flows are simultaneously simulated in the time-marching solution procedure for the Navier-Stokes equation. The porosity technique and the density function are devised for the implementation of the three-dimensional body-boundary and the free-surface conditions, respectively. In this paper, a craft is modeled simply by three blocks containing a wing mounted on a main body horizontally, with the endplate. The numerical calculations of a WIG advancing in a calm water are performed and the WIG-generated wave profiles are also obtained. In the final paper, details of the numerical methods employed for the present study and calculated results are discussed.

• Journal of computational fluids engineering
• /
• v.5 no.1
• /
• pp.33-42
• /
• 2000

The LiBr-H₂O absorption process on a horizontal tube has been analyzed numerically. The flow field, which was calculated in the authors' previous study by solving the fully elliptic Navier-Stokes equations with accurate free-surface-tracking method, is used to solve the temperature and concentration distributions in the absorption film. With the assumption that the absorbent is linear, calculations have been made for various inlet temperature and flow-rate conditions. For low inlet temperature, the absorption rate is large in the upstream region but the mean temperature also increases and as a result the absorption decreases as the film flows to downstream while high-inlet-temperature case does the opposite. The difference in the absorption rate due to the inlet temperature change becomes smaller in the downstream than that in the upstream. For large flow rate, the heat transfer to the wall becomes poor due to the thick film and so does the absorption rate. The analyses have also been carried out for multiple tube arrangement and the results show that the absorption rate converges after a few tube rows.

• Ocean and Polar Research
• /
• v.25 no.3
• /
• pp.287-302
• /
• 2003

The seasonal variations in the circulation of the water mass in the East Sea/Japan Sea have been simulated using a free surface primitive ocean model, RIAMOM (RIAM Ocean Model), comparing the results from GFDL-MOM1 (Geophysical Fluid Dynamics Laboratory Modular Ocean Model, version 1.1, hereafter MOM) with the GDEM (Generalized Digital Environmental Model) data. Both models appear to successfully reproduce the distinct features of circulation in the East Sea/Japan Sea, such as the NB (Nearshore Branch) flowing along the Japanese coast, the EKWC (East Korean Warm Current) flowing northward along the Korean coast, and the NKCC/LCC (North Korean Cold Current/Liman Cold Current) flowing southwestward along Korean/Russian coast. RIAMOM has shown better performance, compared to MOM, in terms of the realistic simulation of the flow field in the East Sea/Japan Sea; RIAMOM has produced more rectified flows on the coastal region, for example, the narrower and stronger NKCC/LCC than MOM has. There is however obvious differences between the model results and the GDEM data in terms of the calculation of the water mass; both models have shown a tendency to overpredict temperature and underpredict salinity below 50m; more diffusive forms of thermocline and halocline have been simulated than noted in GDEM data.

• Journal of the Korean Society of Visualization
• /
• v.8 no.1
• /
• pp.31-38
• /
• 2010

An experimental study to evaluate dynamic structures of flow and turbulence characteristics in bubble-driven liquid flow in a rectangular tank with a varying flow rate of compressed air is conducted. Liquid flow fields are measured by time-resolved particle image velocimetry (PIV) with fluorescent tracer particles to eliminate diffused reflections, and by an image intensifier to acquire enhanced clean particle images. Instantaneous vector fields are investigated by using the two frame cross-correlation function and bad vectors are eliminated by magnitude difference technique. By proper orthogonal decomposition (POD) analysis, the energy distributions of spatial and temporal modes are acquired. When Reynolds number increases, bubble-induced turbulent motion becomes dominant rather than the recirculating flow near the side wall. The total kinetic energy transferred to the liquid from the rising bubbles shows a nonlinear relation regarding the energy input because of the interaction between bubbles and free surface.

• Ocean Systems Engineering
• /
• v.1 no.1
• /
• pp.31-57
• /
• 2011

Impact pressure due to sloshing is of great concern for the ship owners, designers and builders of the LNG carriers regarding the safety of LNG containment system and hull structure. Sloshing of LNG in partially filled tank has been an active area of research with numerous experimental and numerical investigations over the past decade. In order to accurately predict the sloshing impact load, a new numerical method was developed for accurate resolution of violent sloshing flow inside a three-dimensional LNG tank including wave breaking, jet formation, gas entrapping and liquid-gas interaction. The sloshing flow inside a membrane-type LNG tank is simulated numerically using the Finite-Analytic Navier-Stokes (FANS) method. The governing equations for two-phase air and water flows are formulated in curvilinear coordinate system and discretized using the finite-analytic method on a non-staggered grid. Simulations were performed for LNG tank in transverse and longitudinal motions including horizontal, vertical, and rotational motions. The predicted impact pressures were compared with the corresponding experimental data. The validation results clearly illustrate the capability of the present two-phase FANS method for accurate prediction of impact pressure in sloshing LNG tank including violent free surface motion, three-dimensional instability and air trapping effects.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.4
• /
• pp.1182-1196
• /
• 2014

In problems of global hydroelastic ship response in severe seas including the whipping problem, we need to know the hydrodynamic forces acting on the ship hull during almost arbitrary ship motions. In terms of ship sections, some of them can enter water but others exit from water. Computations of nonlinear free surface flows, pressure distributions and hydrodynamic forces in parallel with the computations of the ship motions including elastic vibrations of the ship hull are time consuming and are suitable only for research purposes but not for practical calculations. In this paper, it is shown that the slamming forces can be decomposed in two components within three semi-analytical models of water entry. Only heave motion is considered. The first component is proportional to the entry speed squared and the second one to the body acceleration. The coefficients in these two components are functions of the penetration depth only and can be precomputed for given shape of the body. During the exit stage the hydrodynamic force is proportional to the acceleration of the body and independent of the body shape for bodies with small deadrise angles.