• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.27-41
• /
• 2008

An unstructured hybrid mesh flow solver has been developed for the simulation of three-dimensional steady and unsteady incompressible flow fields. The incompressible Navier-Stokes equations with an artificial compressibility method were discretized by using a node-based finite-volume method. For the unsteady time-accurate computation, a dual-time stepping method was adopted to satisfy a divergence-free flow field at each physical time step. An implicit time integration method with local time stepping was implemented to accelerate the convergence in the pseudo-time sub-iteration procedure. The one-equation Spalart-Allmaras turbulence model has been adopted to solve high-Reynolds number flow fields. The flow solver was parallelized to minimize the CPU time and to overcome the computational overhead. This method has been applied to calculate steady and unsteady flow fields around submarine configurations and a 3-D infinite cylinder. Validations were made by comparing the predicted results with those of experiments or other numerical results. It was demonstrated that the present method is efficient and robust for the prediction of steady and unsteady incompressible flow fields.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.306-320
• /
• 2010

Modern ultra-large container carriers can be exposed to the unprecedented springing excitation from ocean waves due to their relatively low torsional rigidity. Large deck opening on the deck of container carriers tends to cause warping distortion of hull structure under wave-induced excitation, eventually leading to the higher chance of resonance vibration between its torsional response and incoming waves. To handle this problem, a higher-order B-spline Rankine panel method and Vlasov-beam FE model was directly coupled in the time domain, and the coupled equation was solved by using an implicit iterative method. In order to capture the complicated behavior of thin-walled open section girder, a sophisticated beam-based finite element model was developed, which takes into account warping distortion and shear-on-wall effect. Then, the developed beam model was directly coupled with the time-domain Rankine panel method for hydrodynamic problem by using the fixed-point iteration method. The developed computational scheme was validated through the comparison with the frequency-domain solution on the container carrier model in linear springing regime.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.10a
• /
• pp.323-328
• /
• 2003

A basic design on a flexible riser in a floating type development system for upwelling deep ocean water is presented. In the numerical study, an implicit finite difference algorithm is employed for three-dimensional riser equations. Fluid and geometric non-linearity and bending stiffness are considered and solved by Newton-Raphson iteration. To keep the depth of end point of a flexible and light riser is very important for upwelling deep ocean water in a floating type development system. Weight attached at the end point of the riser in order to keep its intake depth is designed under the strong surface current and the configuration of the riser is predicted. The results of this study can be contributed to the design of the development system in floating type for upwelling deep ocean water.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2002.05a
• /
• pp.161-166
• /
• 2002

Low-tension cables have been increasingly used in recent years due to deep-sea developments and the advent of synthetic cables. In the case of low-tension cables, large displacements may happen due to relatively small restoring forces of tension and thus the effects of fluid and geometric non-linearities become predominant. In this study, three-dimensional (3-D) dynamic behavior of a towed low-tension cable with non-uniform characteristics is numerically analyzed by considering fluid and geometric non-linearities and bending stiffness. A Fortran program is developed by employing a finite difference method. In the algorithm, an implicit time integration and Newton-Raphson iteration are adopted. For the calculation of huge size of matrices, block tri-diagonal matrix method is applied, which is much faster than the well-known Gauss-Jordan method in two point boundary value problems. Some case studies are carried out and the results of numerical simulations are compared with a in-house program of WHOI Cable with good agreements.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.7 s.100
• /
• pp.836-842
• /
• 2005

Nonlinear dynamic characteristics of active piezolaminated plates are investigated with respect to the thermopiezoelastic behaviors. For largely deformed structures with small strain, the incremental total Lagrangian formulation is presented based on the virtual work principles. A multi-field layer-wise finite shell element is proposed for assuring high accuracy and non-linearity of displacement, electric and thermal fields. For dynamic consideration of thermopiezoelastic snap-through phenomena, the implicit Newmark's scheme with the Newton-Raphson iteration is implemented for the transient response of various piezolaminated models with symmetric or eccentric active layers. The bifurcate thermal buckling of symmetric structural models is first investigated and the characteristics of piezoelectric active responses are studied for finding snap-through piezoelectric potentials and the load-path tracking map. The thermoelastic stable and unstable postbuckling, thermopiezoelastic snap-through phenomena with several attractors are proved using the nonlinear time responses for various initial conditions and damping loss factors. Present results show that thermopiezoelastic snap-through phenomena can result in the difficulty of buckling and postbuckling control of intelligent structures.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.12
• /
• pp.1624-1634
• /
• 1997

A numerical study on two-dimensional unsteady transonic cascade flow has been performed by adopting dual time stepping and the k-.omega. turbulence model. An explicit 4 stage Runge-Kutta scheme for the compressible Navier-Stokes equations and an implicit Gauss-Seidel iteration scheme for the k-.omega. turbulence model are proposed for fictitious time stepping. This mixed time stepping scheme ensures the stability of numerical computation and exhibits a good convergence property with less computation time. Typical steady-state convergence accelerating schemes such as local time stepping, residual smoothing and multigrid combined with dual time stepping shows good convergence properties. Numerical results are presented for unsteady laminar flow past a cylinder and turbulent shock buffeting problem for bicircular arc cascade flow is discussed.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.5
• /
• pp.405-412
• /
• 2012

A flow control on airfoil installed a vortex cell for high efficiency wind turbine blade in stationary and dynamic stall conditions have been numerically investigated by solving the compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with Roe's flux-difference splitting and an implicit time-integration method coupled with dual time step sub-iteration. The computed result for the airfoil in the stationary showed that lift-drag ratio increases due to low pressure by the vortex cell. The oscillating airfoil with the vortex cell showed that the magnitude of hysteresis loop is reduced due to the enhanced vortex in the cell.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.11a
• /
• pp.29-34
• /
• 2004

Numerical and experimental studies on pipeline laying for intake Deep Ocean Water are carried out. In the numerical study, an implicit finite difference algorithm is employed for three-dimensional pipe equations. Fluid non-linearity and bending stiffness are considered and solved by Newton-Raphson iteration. Seabed is modeled as elastic foundation with linear spring and damper. Top tension and general configuration of pipeline at a depth are predicted. It is found that control for tension to prevent being large curvature of pipeline is needed on th steep seabed and, it should be considered 23.5 ton of tension at a top of pipe on the process of pipeline laying at 400m of water depth The largest top tension of pipe on condition of the beam sea during pipe laying is shown from the experiment. The results of this study can be contributed to the design of pipeline laying for upwelling deep ocean water.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2002.10a
• /
• pp.191-196
• /
• 2002

The paper presents a numerical and experimental investigation on dynamic behaviors of a towed low tension cable. In the numerical study, an implicit finite difference algorithm is employed for three-dimensional cable equations. Fluid and geometric non-linearity and bending stiffness are considered and solved by Newton-Raphson iteration. Block tri-diagonal matrix method is applied for the fast calculation of the huge size of matrices. In order to verify the numerical results and to see real physical phenomena, an experiment is carried out for a 6m cable in a deep and long towing tank. The cable is towed in two different ways; one is towed at a constant speed and the other is towed at a constant speed with top end horizontal oscillations. Cable tension and shear forces are measured at the top end. Numerical and experimental results are compared with good agreements in most cases but with some differences in a few cases. The differences are due to drag coefficients caused by vortex shedding. In the numerical modeling, non-uniform element length needs to be employed to cope with the sharp variation of tension and shear forces at near top end.

• Smart Structures and Systems
• /
• v.12 no.3_4
• /
• pp.235-250
• /
• 2013

Optimal sensor placement (OSP) technique plays a key role in the structural health monitoring (SHM) of large-scale structures. According to the mathematical background and implicit assumptions made in the triaxial effective independence (EfI) method, this paper presents a novel multi-dimensional OSP method for the Canton Tower focusing on application demands. In contrast to existing methods, the presented method renders the corresponding target mode shape partitions as linearly independent as possible and, at the same time, maintains the stability of the modal matrix in the iteration process. The modal assurance criterion (MAC), determinant of the Fisher Information Matrix (FIM) and condition number of the FIM have been taken as the optimal criteria, respectively, to demonstrate the feasibility and effectiveness of the proposed method. Numerical investigations suggest that the proposed method outperforms the original EfI method in all instances as expected, which is looked forward to be even more pronounced should it be used for other multi-dimensional optimization problems.