• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.159-173
• /
• 1991

This paper describes a potential-based panel method formulated for the analysis of a super-cavitating two-dimensional hydrofoil. The method employs normal dipoles and sources distributed on the foil and cavity surfaces to represent the potential flow around the cavitating hydrofoil. The kinematic boundary condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the fictitious inner flow region of the foil, and the dynamic boundary condition on the cavity surface is satisfied by requiring thats the potential vary linearly, i.e., the tangential velocity be constant. Green's theorem then results in a potential-based integral equation rather than the usual velocity-based formulation of Hess & Smith type. With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with improved accuracy compared to those of the linearized lilting surface theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. This iteration process is repeated until the results are converged. Characteristics of iteration and discretization of the present numerical method are much faster and more stable than the existing nonlinear theories. The theory shows good correlations with the existing theories and experimental results for the super-cavitating flow. In the region of small angles of attack, the present prediction shows and excellent comparison with the Geurst's linear theory. For the long cavity, the method recovers the trends of the Wu's nonlinear theory. In the intermediate regions of the short super-cavitation, the method compares very well with the experimental results of Parkin and also those of Silberman.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.1
• /
• pp.11-22
• /
• 1987

A plane buoyant jet discharged vertically upward into a crossflow is analyzed by numerical solution of the governing equations of continuity, momentum and constituent transport. The turbulent transport is modelled by the Prandtl's mixing length theory. In the numerical solution procedure, the governing equations are transformed by stream function and vorticity transport, non-dimensionalyzed by discharge velocity, slot width, and parameters representing flow characteristics, and solved by Gauss-Seidel iteration method with successive underrelaxation. The numerical experiments were performed for the region of established flow of buoyant jet in the range of discharge densimetric Froude number of 4 to 32 and in the range of velocity ratio of 8 to 15, which is the ratio of discharge velocity to crossflow velocity. Variations of velocities and temperatures, flow patterns and vorticity patterns of receiving water due to buoyant jet were investigated. Also investigated are the effects of velocity ratio and discharge densimetric Froude number on the trajectories of buoyant jet. Computed are velocities, temperatures and local densimetric Froude numbers along the trajectory of the buoyant jet. Spreading rate and dispersion ratio were analyzed in terms of discharge densimetric Froude number, local densimetric Froude number and distance from the source along the jet trajectory. It was noted that the similarity law holds in both the profiles of velocity and temperatures across the jet trajectory and the integral type analysis of Gaussian distribution is applicable.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.1A
• /
• pp.31-43
• /
• 2009

In this study, the numerical method is presented, which can consider the various train types and can solve the equations of motion for a vehicle-bridge interaction analysis by non-iteration procedure through formulating the coupled equations of motion. The coupled equations of motion for the vehicle-bridge interaction are solved by the Newmark ${\beta}$ of a direct integration method, and by composing the effective stiffness matrix and the effective force vector according to a analysis step, those can be solved with the same manner of the solving procedure of equilibrium equations in static analysis. Also, the effective stiffness matrix is reconstructed by the Skyline method for increasing the analysis effectiveness. The Cholesky's matrix decomposition scheme is applied to the analysis procedure for minimizing the numerical errors that can be generated in directly calculating the inverse matrix. The equations of motion for the conventional trains are derived, and the numerical models of the conventional trains are idealized by a set of linear springs and dashpots with 16 degrees of freedom. The bridge models are simplified by the 3 dimensional space frame element which is based on the Euler-Bernoulli theory. The rail irregularities of vertical and lateral directions are generated by the PSD functions of the Federal Railroad Administration (FRA). The results of the vehicle-bridge interaction analysis are verified by the experimental results for the railway plate girder bridges of a span length with 12 m, 18 m, and the experimental and analytical data are applied to the low pass filtering scheme, and the basis frequency of the filtering is a 2 times of the 1st fundamental frequency of a bridge bending.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.4
• /
• pp.312-326
• /
• 2012

This paper considers a fully coupled 3D BEM-FEM analysis for the ship structural hydroelasticity problem in waves. Fluid flows and structural responses are analyzed by using a 3D Rankine panel method and a 3D finite element method, respectively. The two methods are fully coupled in the time domain using a fixed-point iteration scheme, and a relaxation scheme is applied for improve convergence. In order to validate the developed method, numerical tests are carried out for a barge model. The computed natural frequency, motion responses, and time histories of stress are compared with the results of the beam-based hydroelasticity program, WISH-FLEX, which was thoroughly validated in previous studies. This study extends to a real-ship application, particularly the springing analysis for a 6500 TEU containership. Based on this study, it is found that the present method provides reliable solutions to the ship hydroelasticity problems.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.317-324
• /
• 1994

The circulation of Lake Paldang is analysed numerically as an ultimate goal to develop a vehicle predicting the dispersion and concentration of pollutants and sediment flowed into the lake. In finite difference formulation of 2-D depth averaged governing equations. Abbott's 3-time level scheme is employed and for nonlinear terms time centering iteration technique in time and space is used. Model parameters for shear stresses and eddy diffusivities are determined through measured data in and near the lake. Predicted velocities for steady flow are shown to be close to the measured velocities and further improved by taking into account of wind effect. This indicates that the wind effect is needed for proper circulation analysis and it calls for the inclusion of the wind effect. Simulated results of unsteady flow caused by flood inflows and release through Paldang dam are found to characterize the flow features quite well as expected. This implies that the developed model can be used as a tool to analysing the circulation in the lake.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.14 no.4
• /
• pp.525-535
• /
• 2001

In this paper, a new adaptive analysis scheme for element-free Galerkin method(EFGM) is proposed. The novel point of this scheme is that the triangular cell structure based on the Delaunay triangulation is used in the numerical integration and the node adding/removing process. In adaptive analysis with this scheme, there is no need to divide the integration cell and the memory cell structure. For the adaptive analysis of crack propagation, the reconstruction of cell structure by adding and removing the nodes on integration cells based the estimated error should be carried out at every iteration step by the Delaunay triangulation technique. This feature provides more convenient user interface that is closer to the real mesh-free nature of EFGM. The analysis error is obtained basically by calculating the difference between the values of the projected stresses and the original EFG stresses. To evaluate the performance of proposed adaptive procedure, the crack propagation behavior is investigated for several examples.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.41 no.2
• /
• pp.9-15
• /
• 2004

We present a new unified analytical front surface potential model, which can accurately describe the transitions between the partially-depleted (PD) and the fully-depleted (FD) regimes with an analytical expression for the critical voltage V$_{c}$ delineating the PD and the FD region. It is valid in all regions of operation (from the sub -threshold to the strong inversion) and has the shorter calculation time than the iterative procedure approach. A charge sheet model based on the above explicit surface potential formulation is used to derive a single formula for the drain current valid in all regions of operation. Most of the secondary effects can be easily included in the charge sheet model and the model accurately reproduces various numerical and experimental results. No discontinuity in the derivative of the surface potential is found even though three types of smoothing functions are used. More importantly, the newly introduced parameters used in the smoothing functions do not strongly depend on the process parameter.

• Journal of the Korean earth science society
• /
• v.37 no.5
• /
• pp.277-285
• /
• 2016

We introduce a depth scaling strategy to improve the accuracy of frequency-domain elastic full waveform inversion (FWI) using the new pseudo-Hessian matrix for seismic data without low-frequency components. The depth scaling strategy is based on the fact that the damping factor in the Levenberg-Marquardt method controls the energy concentration in the gradient. In other words, a large damping factor makes the Levenberg-Marquardt method similar to the steepest-descent method, by which shallow structures are mainly recovered. With a small damping factor, the Levenberg-Marquardt method becomes similar to the Gauss-Newton methods by which we can resolve deep structures as well as shallow structures. In our depth scaling strategy, a large damping factor is used in the early stage and then decreases automatically with the trend of error as the iteration goes on. With the depth scaling strategy, we can gradually move the parameter-searching region from shallow to deep parts. This flexible damping factor plays a role in retarding the model parameter update for shallow parts and mainly inverting deeper parts in the later stage of inversion. By doing so, we can improve deep parts in inversion results. The depth scaling strategy is applied to synthetic data without lowfrequency components for a modified version of the SEG/EAGE overthrust model. Numerical examples show that the flexible damping factor yields better results than the constant damping factor when reliable low-frequency components are missing.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.69-77
• /
• 1991

A collision-free trajectory planning algorithm using an iterative learning concept is proposed for dual robot arms in a 3-D common workspace to accurately follow their specified paths with constant velocities. Specifically, a collision-free trajectory minimizing the trajectory error is obtained first by employing the linear programming technique. Then the total operating time is iteratively adjusted based on the maximum trajectory error of the previous iteration so that the collision-free trajectory has no deviation from the specified path and also that the operating time is near-minimal. To show the validity of the proposed algorithm, a numerical example is presented based on two planar robots.

• The KIPS Transactions:PartB
• /
• v.16B no.4
• /
• pp.289-298
• /
• 2009

Active Contour Model, that is, Snake algorithm is effective for detection and tracking the objects. However, this algorithm has some drawbacks; numerous parameters must be designed(weighting factors, iteration steps, etc.), a reasonable initialization must be available and moreover suffers from numerical instability. Therefore we propose a novel Energy Corrected Snake(ECS) algorithm which improved on external energy of Snake algorithm for detection and tracking the moving object more effectively. The proposed algorithm uses the difference image, getting when the object is moving. It copies four direction images from the difference image and performs the accumulating compute to erasing image noise, so that it gets external energy steadily. Then external energy united with contour that is computed by internal energy. Consequently we can detect and track the moving object more speedily and easily. To show the effectiveness of the proposed algorithm, we experiment on 3 situations. The experimental results showed that the proposed algorithm outperformed by 6$\sim$9% of detection rate and 6$\sim$11% of tracker detection rate compared with the Snake algorithm.