• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.12
• /
• pp.1738-1746
• /
• 2002

In the FSI (Fluid-Structure Interaction) problems, two different governing equations are to be solved together. One is fur the fluid and the other for the structure. Furthermore, a kinematic constraint should be imposed along the boundary between the fluid and the structure. We use the combined formulation, which incorporates both the fluid and structure equations of motion into a single coupled variational equation so that it is not necessary to calculate the fluid force on the surface of structure explicitly when solving the equations of motion of the structure. A two-dimensional channel flow divided by a Bernoulli-Euler beam is considered and the dynamic response of the beam under the influence of channel flow is studied. The Navier-Stokes equations are solved using a P2P1 Galerkin finite element method with ALE (Arbitrary Lagrangian-Eulerian) algorithm. The internal structural damping effect is not considered in this study and numerical results are compared with a previous work fer steady case. In addition to the Reynolds number, two non-dimensional parameters, which govern this fluid-structure system, are proposed. It is found that the larger the dynamic viscosity and density of the fluid are, the larger the damping of the beam is. Also, the added mass is found to be linearly proportional to the density of the fluid.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.5
• /
• pp.476-483
• /
• 2007

One of the main function of Ieodo Ocean Research Station is to service the information about the weather and fishing grounds condition which are collected through calibrating convection flow and ocean current around the station. However, due to the influence of the station's structure below sea level, it is difficult to obtain the exact flow data. Therefore, it is required to research on the effect of the structure and the method to evaluate and revise the observed data. In this paper, as a basic study, it deals with the algorithm that simulate the interaction between ocean current and the station structure, followed by discussions about the way to applicate the algorithm. Through 3-dimensional computational fluid dynamics analyses (using Navier-Stokes equations with K-turbulence model), the influence of the station and submerged rocks are quantitatively evaluated, and we would suggest methods how to obtain accurate flow information from the measured rough data.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.7
• /
• pp.1-8
• /
• 2005

For study on the unsteady wall interference effect, flows around a circular cylinder in closed test-section wind tunnels have been numerically investigated by solving compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with the Roe's flux-difference splitting and an implicit time-integration method coupled with dual time-step sub-iteration. The computed results showed that the unsteady pressure gradient over the cylinder is enhanced by the wall interference, and as a result the fluctuations of lift and drag are augmented. The drag is further increased because of the lower base pressure. The vortex shedding frequency is also increased by the wall interference. The pressure on the test section wall shows the harmonics having the shedding frequency contained in the wall effect.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.2
• /
• pp.875-882
• /
• 2019

The effects of Triangle Groove Strips (TGS) on Vortex-induced Vibration (VIV) suppression of marine riser are numerically investigated using Computational Fluid Dynamics (CFD) method. The range of Reynolds number in simulations is 4.0 × 10⁴ < Re < 1.2 × 10⁵. The two-dimensional unsteady Reynolds-Averaged Navier-Stokes (RANS) equations and Shear Stress Transport (SST) k-ω turbulence model are used to calculate the flow around marine riser. The Newmark-β method is employed for evaluating the structure dynamics of marine riser. The effect of the height ratio (ε) of TGS on VIV suppression is evaluated. The amplitude responses, frequency responses, vortex patterns and the flow around the structures are discussed in detail. With the increase of the height ratio of TGS, the suppression effect of TGS on VIV suppression is improved firstly and then weakened. When ε=0.04, the suppression effect of TGS is the best. Compared with the VIV responses of smooth marine riser, the amplitude ratio is reduced by 38.9%, the peak of the lift coefficient is reduced by 69% and the peak of the drag coefficient is reduced by 40% when Re=6.0 × 10⁴. With the increase of Reynolds number, the suppression effect of TGS on VIV suppression is improved firstly and then weakened. When the Reynolds number is 7.0 × 10⁴, the amplitude ratio can be reduced by 40.1%. As to the large-amplitude vibration cases, the TGS show nice suppression effect on VIV.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.28-40
• /
• 2002

Turbulent flows around two-dimensional wing sections in ground effect are analysed by incompressible RANS equations and a finite difference method. The Baldwin-Lomax algebraic turbulence model is used to simulate high Reynolds number flows. The main purpose of this study is to clarify the two-dimensional ground effect and its flow characteristics due to different ground boundary conditions, i.e., moving and fixed bottom boundary. As a first step, to validate the present numerical code, the computational result of Clark-Y(t/C 11.7%) is compared with published numerical results and experimental data. Then, NACA4412 section in ground effect is calculated for various ground clearances with two bottom boundary conditions. According to the computational results, the difference in the lift and moment simulated with the two bottom boundary conditions is negligible, but the drag force simulated by the fixed bottom is to some extent smaller than that by the moving bottom. Therefore, it can be concluded that the drag force measured in a wind tunnel with the fixed bottom could be smaller than that with the moving bottom.

• International Journal of Fluid Machinery and Systems
• /
• v.9 no.4
• /
• pp.370-381
• /
• 2016

As a single-channel pump is used for wastewater treatment, this particular pump type can prevent performance reduction or damage caused by foreign substances. However, the design methods for single-channel pumps are different and more difficult than those for general pumps. In this study, a design optimization method to improve the hydrodynamic performance of a single-channel pump impeller is implemented. Numerical analysis was carried out by solving three-dimensional steady-state incompressible Reynolds-averaged Navier-Stokes equations using the shear stress transport turbulence model. As a state-of-the-art impeller design method, two design variables related to controlling the internal cross-sectional flow area of a single-channel pump impeller were selected for optimization. Efficiency was used as the objective function and was numerically assessed at twelve design points selected by Latin hypercube sampling in the design space. An optimization process based on a radial basis neural network model was conducted systematically, and the performance of the optimum model was finally evaluated through an experimental test. Consequently, the optimum model showed improved performance compared with the base model, and the unstable flow components previously observed in the base model were suppressed remarkably well.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.4
• /
• pp.274-280
• /
• 2017

Generally, global performance analysis in offshore platforms is performed using potential-based numerical tools, which neglect hydrodynamic viscous effects. In comparison with the potential theory, computational fluid dynamics (CFD) methods can take into account the viscous effects by solving the Navier-Stokes equation using the finite-volume method. The open-source field operation and manipulation (OpenFOAM) C++ libraries are employed for a finite volume method (FVM) numerical analysis. In this study, in order to apply CFD to the global performance analysis of a hull-mooring coupled system, we developed a numerical wave basin to analyze the global performance problem of a floating body with a catenary mooring system under regular wave conditions. The mooring system was modeled using a catenary equation and solved in a quasi-static condition, which excluded the dynamics of the mooring lines such as the inertia and drag effects. To demonstrate the capability of the numerical basin, the global performance of a barge with four mooring lines was simulated under regular wave conditions. The simulation results were compared to the analysis results from a commercial mooring analysis program, Orcaflex. The comparison included the motion of the barge, catenary shape, and tension in the mooring lines. The study found good agreement between the results from the developed CFD-based numerical calculation and commercial software.

• Smart Structures and Systems
• /
• v.20 no.6
• /
• pp.709-728
• /
• 2017

This disquisition proposes a nonlocal strain gradient beam theory for thermo-mechanical dynamic characteristics of embedded smart shear deformable curved piezoelectric nanobeams made of porous electro-elastic functionally graded materials by using an analytical method. Electro-elastic properties of embedded curved porous FG nanobeam are assumed to be temperature-dependent and vary through the thickness direction of beam according to the power-law which is modified to approximate material properties for even distributions of porosities. It is perceived that during manufacturing of functionally graded materials (FGMs) porosities and micro-voids can be occurred inside the material. Since variation of pores along the thickness direction influences the mechanical and physical properties, so in this study thermo-mechanical vibration analysis of curve FG piezoelectric nanobeam by considering the effect of these imperfections is performed. Nonlocal strain gradient elasticity theory is utilized to consider the size effects in which the stress for not only the nonlocal stress field but also the strain gradients stress field. The governing equations and related boundary condition of embedded smart curved porous FG nanobeam subjected to thermal and electric field are derived via the energy method based on Timoshenko beam theory. An analytical Navier solution procedure is utilized to achieve the natural frequencies of porous FG curved piezoelectric nanobeam resting on Winkler and Pasternak foundation. The results for simpler states are confirmed with known data in the literature. The effects of various parameters such as nonlocality parameter, electric voltage, coefficient of porosity, elastic foundation parameters, thermal effect, gradient index, strain gradient, elastic opening angle and slenderness ratio on the natural frequency of embedded curved FG porous piezoelectric nanobeam are successfully discussed. It is concluded that these parameters play important roles on the dynamic behavior of porous FG curved nanobeam. Presented numerical results can serve as benchmarks for future analyses of curve FG nanobeam with porosity phases.

• Journal of the Korean Society of Environmental Restoration Technology
• /
• v.12 no.3
• /
• pp.98-105
• /
• 2009

The objective of this study is to develop numerical simulation model for analysing the wind speed reduction effect by coastal forest belts. The horizontally homogeneous turbulent flow equations, which are derived from the Reynolds-averaged Navier-Stokes method, both above the tree canopy and within the canopy were first formulated, and a first-order closure scheme with the capability of accounting the bulk momentum transport term within the canopy was employed. The averaged equations were solved numerically by finite difference method, FTCS (forward time centered space) scheme. The proposed model was also used to numerically investigate the effects of structural characteristic of forest belt on the wind speed. The effects of maximum leaf area density were evaluated, with the leaf area density of $1.0m^2/m^3$, $2.0m^2/m^3$, $3.0m^2/m^3$, and $4.0m^2/m^3$. Vertical distributions of leaf area, both uniform and varied distribution with a height, were also considered. A comparison of wind profile indicated that there was in good agreements between simulated and measured wind speed. Also, the results showed horizontal wind speed decreased under a height of the tree with increasing maximum leaf area density. In conclusion, in applications where computational efficiency and simplicity are desirable, the proposed numerical model has of great capability to determine the vertical turbulent momentum transport and wind profile in the costal forest belt.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.1
• /
• pp.1-13
• /
• 2008

In the present work, wave transformation by vertical columns and its wave forces acting on them are discussed using a direct 3-D numerical model based on the VOF (Volume Of Fluid) method. The numerical results for wave transformations and wave forces are critically compared to an advanced experimental data, and provide the verification of the numerical model used in the present study. Overall model-data comparisons are good. After verification of the numerical model, it is used to simulate wave fields around dual vertical columns with arbitrary cross section, and the characteristics of nonlinear wave forces and wave transformations according to the variations of different cross section types of vertical columns, an interval of vertical columns and incident wave angle are discussed.