• 대한조선학회논문집
• /
• 제31권3호
• /
• pp.80-86
• /
• 1994

자유수면을 항주하는 선박에 의하여 발생되는 비선형 조파현상을 해석하기 위한 수치해석법을 개발하였다. 유동은 비점성, 비압축성으로 가정하고 선체 및 자유수면 형상과 일치하는 좌표계의 생성을 위하여 타원형 편미분방정식을 수치해석하여 물체적합 좌표계를 생성하였으며 변환된 정규격자 물체적합 좌표계에 대한 Euler방정식을 유한차분법(Finite Difference Method)을 이용하여 계산하였다. 수치해석을 위하여 시간에 대한 미분항은 전진차분, 공간에 대한 미분항은 중심차분법으로 이산화하였고 대류항에는 수치계산의 안정을 위해 인위적인 소산(dissipation)항을 첨가하였다. 자유수면의 형상은 매 시간 단계마다 자유수면 경계조건들을 만족하도록 다시 계산되었고 격자점들은 자유수면형상의 변화에 적합하게 다시 생성되도록 하였으며 압력에 대한 Poisson방정식은 반복연산법에 의하여 풀고 그 결과를 이용하여 속도를 외삽하였다. 개발된 수치해석법의 검증을 위해 수식선형인 Wigley 모형에 대한 계산을 Fn=0.250-0.408에 대하여 수행하고, 그 결과를 실험 결과와 비교하여 잘 일치함을 보였다.

• Journal of Ship and Ocean Technology
• /
• 제10권4호
• /
• pp.12-23
• /
• 2006

In this paper, a numerical study is carried out to investigate the turbulent flow around a twin-skeg container ship model with rudders including propeller effects. A commercial CFD code, FLUENT is used with body forces distributed on the propeller disk to simulate the ship stem and wake flows with the propeller in operation. A multi-block, matching, structured grid system has been generated for the container ship hull with twin-skegs in consideration of rudders and body-force propeller disks. The RANS equations for incompressible fluid flows are solved numerically by using a finite volume method. For the turbulence closure, a Reynolds stress model is used in conjunction with a wall function. Computations are carried out for the bare hull as well as the hull with appendages of a twin-skeg container ship model. For the bare hull, the computational results are compared with experimental data and show generally a good agreement. For the hull with appendages, the changes of the stem flow by the rudders and the propellers have been analyzed based on the computed result since there is no experimental data available for comparison. It is found the flow incoming to the rudders has an angle of attack due to the influence of the skegs and thereby the hull surface pressure and the limiting streamlines are changed slightly by the rudders. The axial velocity of the propeller disk is found to be accelerated overall by about 35% due to the propeller operation with the rudders. The area and the magnitude of low pressure on the hull surface enlarge with the flow acceleration caused by the propeller. The propellers are found to have an effect on up to the position where the skeg begins. The propeller slipstream is disturbed strongly by the rudders and the flow is accelerated further and the transverse velocity vectors are weakened due to the flow rectifying effect of the rudder.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제4권2호
• /
• pp.151-161
• /
• 2012

The present study numerically investigates the aerodynamic characteristics of two-dimensional wings in the vicinity of the ground by solving two-dimensional steady incompressible Navier-Stokes equations with the turbulence closure model of the realizable k-${\varepsilon}$ model. Numerical simulations are performed at a wide range of the normalized ground clearance by the chord length ($0.1{\leq}h/C{\leq}1.25$) for the angles of attack ($0^{\circ}{\leq}{\alpha}{\leq}10^{\circ}$) in the prestall regime at a Reynolds number (Re) of $2{\times}10^6$ based on free stream velocity $U_{\infty}$ and the chord length. As the physical model of this study, a cambered airfoil of NACA 4406 has been selected by a performance test for various airfoils. The maximum lift-to-drag ratio is achieved at ${\alpha}=4^{\circ}$ and h / C = 0.1. Under the conditions of ${\alpha}=4^{\circ}$ and h / C = 0.1, the effect of the Reynolds number on the aerodynamic characteristics of NACA 4406 is investigated in the range of $2{\times}10^5{\leq}Re{\leq}2{\times}10^9$. As Re increases, $C_l$ and $C_d$ augments and decreases, respectively, and the lift-to-drag ratio increases linearly.

• 대한조선학회논문집
• /
• 제58권1호
• /
• pp.1-9
• /
• 2021

In the present study, the model-scale Propeller Open Water (POW) tests for the propeller of 176K bulk carrier and 8600TEU container ship were conducted through Computational Fluid Dynamics (CFD) simulation. In order to solve the incompressible viscous flow field, the Reynolds-averaged Navier-Stokes (RaNS) equations were employed as the governing equations. The γ-Reθ(gamma-Re-theta) transition model combined with the SST k-ωturbulence model was introduced to describe the laminar-turbulence transition considering the low Reynolds number of model-scale. Firstly, the flow simulation developing over a flat plate was performed to verify the transition modeling, in which the wall shear stresses were compared with experiments and other numerical results. Then, to investigate the effect of the model, the CFD simulation for the POW test was performed and the simulated propeller performance was validated through comparison with the experiment conducted at Korea Research Institute of Ships & Ocean Engineering (KRISO).

• Ocean Systems Engineering
• /
• 제8권4호
• /
• pp.381-407
• /
• 2018

This paper presents a hybrid numerical approach, which combines a two-phase Navier-Stokes model (NS) and the fully nonlinear potential theory (FNPT), for modelling wave-structure interaction. The former governs the computational domain near the structure, where the viscous and turbulent effects are significant, and is solved by OpenFOAM/InterDyMFoam which utilising the finite volume method (FVM) with a Volume of Fluid (VOF) for the phase identification. The latter covers the rest of the domain, where the fluid may be considered as incompressible, inviscid and irrotational, and solved by using the Quasi Arbitrary Lagrangian-Eulerian finite element method (QALE-FEM). These two models are weakly coupled using a zonal (spatially hierarchical) approach. Considering the inconsistence of the solutions at the boundaries between two different sub-domains governed by two fundamentally different models, a relaxation (transitional) zone is introduced, where the velocity, pressure and surface elevations are taken as the weighted summation of the solutions by two models. In order to tackle the challenges associated and maximise the computational efficiency, further developments of the QALE-FEM have been made. These include the derivation of an arbitrary Lagrangian-Eulerian FNPT and application of a robust gradient calculation scheme for estimating the velocity. The present hybrid model is applied to the numerical simulation of a fixed horizontal cylinder subjected to a unidirectional wave with or without following current. The convergence property, the optimisation of the relaxation zone, the accuracy and the computational efficiency are discussed. Although the idea of the weakly coupling using the zonal approach is not new, the present hybrid model is the first one to couple the QALE-FEM with OpenFOAM solver and/or to be applied to numerical simulate the wave-structure interaction with presence of current.

• Coupled systems mechanics
• /
• 제8권4호
• /
• pp.339-350
• /
• 2019

A numerical simulation of the incompressible multiphase hydraulic jump flow was performed to compare the interface prediction through the use of the three RANS turbulence models: $k-{\varepsilon}$, $RNGk-{\varepsilon}$ and SST $k-{\omega}$. A three dimensional no submerged hydraulic jump and a two dimensional submerged hydraulic jump were modeled. Both the geometry and the mesh were created using the open source Gmsh code. The project's geometry consists of a rectangular channel with length and height differences between the two dimensional and three dimensional simulations. Uniform hexahedral cells were used for the mesh. Three refining meshes were constructed to allow to verify simulation convergence. The Volume of Fluid (abbr. VOF) method was used for treatment of the air-water surface. The turbulence models were evaluated in three distinct set up configurations to provide a greater accuracy in the flow representation. In the two-dimensional analysis of a submerged hydraulic jump simulation, the turbulence model RNG RNG $k-{\varepsilon}$ provided a better interface adjust with the experimental results than the model $k-{\varepsilon}$ and SST $k-{\omega}$. In the three-dimensional simulation of a no-submerged hydraulic jump the k-# showed better results than the SST $k-{\omega}$ and RNG $k-{\varepsilon}$ capturing the height and length of the ledge with a better fit with the experimental results.

• Smart Structures and Systems
• /
• 제30권3호
• /
• pp.221-237
• /
• 2022

This paper is concerned with nonlinear modeling and efficient numerical simulation of electrostrictive materials and structures. Two types of such materials are considered: relaxor ferroelectric ceramics and electrostrictive polymers. For ceramics, a geometrically linear formulation is developed, whereas polymers are studied in a geometrically nonlinear regime. In the paper, we focus on constitutive modeling first. For the reversible constitutive response under consideration, we introduce the augmented Helmholtz free energy, which is composed of a purely elastic part, a dielectric part and an augmentation term. For the elastic part, we involve an additive decomposition of the strain tensor into an elastic strain and an electrostrictive eigenstrain, which depends on the polarization of the material. In the geometrically nonlinear case, a corresponding multiplicative decomposition of the deformation gradient tensor replaces the additive strain decomposition used in the geometrically linear formulation. For the dielectric part, we first introduce the internal energy, to which a Legendre transformation is applied to compute the free energy. The augmentation term accounts for the contribution from vacuum to the energy. In our formulation, the augmented free energy depends not only on the strain and the electric field, but also on the polarization and an internal polarization; the latter two are internal variables. With the constitutive framework established, a Finite Element implementation is briefly discussed. We use high-order elements for the discretization of the independent variables, which include also the internal variables and, in case the material is assumed incompressible, the hydrostatic pressure, which is introduced as a Lagrange multiplier. The elements are implemented in the open source code Netgen/NGSolve. Finally, example problems are solved for both, relaxor ferroelectric ceramics and electrostrictive polymers. We focus on thin plate-type structures to show the efficiency of the numerical scheme and its applicability to thin electrostrictive structures.

• Wind and Structures
• /
• 제34권5호
• /
• pp.407-419
• /
• 2022

Many factors should be considered by architects and designers for designing a tall building. Wind load is one of these important factors that govern the design of tall building structures and can become a serious challenge when buildings tend to be built very tall and slender. On the other hand, through the initial stages of a design process, choosing the design geometry greatly affects the wind-induced forces on a tall building. With this respect, geometric shapes with 3-fold rotational symmetry are one of the applied plan shapes in tall buildings. This study, therefore, aims to investigate the aerodynamic characteristics of 8 different geometrical shapes using Computational Fluid Dynamics (CFD) by measuring the drag and lift forces. A case study approach was conducted in which different building shape models have the same total gross area and the same height of 300 meters. The simulation was an incompressible transient flow that ran 1700 timesteps (85 seconds on the real-time scale). The results show a great difference between wind-induced force performance of buildings with different plan shapes. Generally, it is stated that the shapes with the same area, but with smaller perimeters, are better choices for reducing the drag force on buildings. Applying the lift force, the results show that the buildings with plan shapes that have rounded corners act better in crosswind flow while, those with sharp corners induce larger forces in the same direction. This study delivers more analytical understanding of building shapes and their behavior against the wind force through the parametric modelling.

• Wind and Structures
• /
• 제35권2호
• /
• pp.83-92
• /
• 2022

When the rolling stocks run on the curve, the external rail has to be lifted to a certain level to balance the centrifugal force acting on the train body. Under such a situation, passengers may feel uncomfortable, and the slanted vehicle has the potential overturning risks at high speed. This paper conducted a wind tunnel test in an annular wind tunnel with φ=3.2 m based on a 1/20^th scaled high-speed train (HST) model. The sensitivity of Reynolds effects ranging from Re = 0.37×10⁶ to Re = 1.45×106 was tested based on the incoming wind from U=30 m/s to U=113 m/s. The wind speed covers the range from incompressible to compressible. The impact of roll angle ranging from γ=0° to γ=4° on train aerodynamics was tested. In addition, the boundary layer development was also analyzed under different wind speeds. The results indicate that drag and lift aerodynamic coefficients gradually stabilized and converged over U=70 m/s, which could be regeared as the self-similarity region. Similarly, the thickness of the boundary layer on the floor gradually decreased with the wind speed increase, and little changed over U=80 m/s. The rolling moment of the head and tail cars increased with the roll angle from γ=0° to γ=4°. However, the potential overturning risks of the head car are higher than the tail car with the increase of the roll angle. This study is significant in providing a reference for the overturning assessment of HST.

• Journal of applied mathematics & informatics
• /
• 제41권1호
• /
• pp.51-70
• /
• 2023

The flow of an incompressible Ellis fluid in an inclined asymmetric channel, driven by peristaltic waves was studied under low Reynolds number and long wavelength assumptions. The wave on each side of the channel are assumed to be an infinite train of sinusoidal waves, both having the same constant wave speed and wavelength however, they vary in wave amplitude, channel half width and phase angle. We derived expressions for the axial and transverse velocities, volume flow rate, pressure rise per unit wavelength and streamlines. The effects of varying the wave amplitudes, the phase angle, the channel width, the angle of inclination of the channel as well as the fluid parameters on the flow were analyzed. Trapping conditions were determined and the presence of reflux highlighted using the streamlines for the necessary channel and fluid conditions. By varying the fluid parameters, changes in the fluid that deviated from the Newtonian case resulted in a reduction in the axial velocity in the neighborhood of the center of the channel and a simultaneous increase in the velocity at the periphery of the channel. A nonlinear relation was observed with the pressure rise and the volume flow rate. This nonlinear relation is more pronounced with an increase in the absolute value of the volume flow rate. For Newtonian fluids a linear relation exists between these two variables. The fluid parameters had little effects on the streamlines. However, variations of the wave amplitudes, volume flow, channel width and phase angle had greater effects on the streamlines and hence the trapped region.