• Bulletin of the Society of Naval Architects of Korea
• /
• v.26 no.4
• /
• pp.27-34
• /
• 1989

A potential-based panel method is formulated for the analysis of a partially cavitating 2-dimensional hydrofoil. The method employs dipoles and sources distributed on the foil surface to represent the lifting and cavity problems, respectively. The kinematic boundry condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the inner flow region of the foil. The dynamic boundary condition on the cavity surface is satisfied by requiring that the potential vary linearly, i.e., the velocity be constant. Green's theorem then results in a potential-based boundary value problem rather than a usual velocity-based formulation. With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with more improved accuracy than the zero-thickness hydrofoil 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. It was found that five iterations are necessary to obtain converged values, while only two iterations are sufficient for engineering purpose.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.454-458
• /
• 2007

Nonpoint source (NPS) pollution is a serious problem causing the degradation of soil and water quality. Concentrated overland flow is the primary transport mechanism for a large amount of NPS pollutants from hillslope areas to downslope areas in a watershed. In this study, a soil erosion model, nLS model, to identify transitional overland flow regions (i.e., ephemeral gully head areas) was developed using the kinematic wave overland flow theory. Spatial data, including digital elevation models (DEMs), soil, and landcover, were used in the GIS-based model algorithm. The model was calibrated and validated using gully head locations in a large agricultural watershed, which were identified using 1-m aerial photography. The model performance was better than two previous approaches; the overall accuracy of the nLS model was 72 % to 87 % in one calibration subwatershed and the mean overall accuracy was 75 to 89 % in four validation subwatersheds, showing that the model well predicted potential transitional erosion areas at different watersheds. However, the user accuracy in calibration and validation was still low. To improve the user accuracy and study the effects of DEM resolution, finer resolution DEMs may be preferred because DEM grid is strongly sensitive to estimating model parameters. Information gained from this study can improve assessing soil erosion process due to concentrated overland flow as well as analyze the effect of microtopographic landscapes, such as riparian buffer areas, in NPS control.

• Ocean Systems Engineering
• /
• v.3 no.3
• /
• pp.219-236
• /
• 2013

The accurate prediction of motion in waves of a marine vehicle is essential to assess the maximum sea state vs. operational requirements. This is particularly true for small crafts, such as Autonomous Surface Vessels (ASV). Two different numerical methods to predict motions of a SWATH-ASV are considered: an inviscid strip theory initially developed at MIT for catamarans and then adapted for SWATHs and new a hybrid strip theory, based on the numerical solution of the radiation forces by an unsteady viscous, non-linear free surface flow solver. Motion predictions obtained by the viscous flow method are critically discussed against those obtained by potential flow strip theory. Effects of viscosity are analyzed by comparison of sectional added mass and damping calculated at different frequencies and for different sections, RAOs and motions response in irregular waves at zero speed. Some relevant conclusions can be drawn from this study: influence of viscosity is definitely non negligible for SWATH vessels like the one presented: amplitude of the pitch and heave motions predicted at the resonance frequency differ of 20% respectively and 50%; in this respect, the hybrid method with fully non-linear, viscous free surface calculation of the radiation forces turns out to be a very valuable tool to improve the accuracy of traditional strip theories, without the burden of long computational times requested by fully viscous time domain three dimensional simulations.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.5
• /
• pp.545-565
• /
• 2018

A numerical code is developed based on potential flow theory to investigate nonlinear sloshing in rectangular Liquefied Natural Gas (LNG) tanks under forced excitation. Using this code, internal free-surface elevation and sloshing loads on liquid tanks can be obtained both in time domain and frequency domain. In the mathematical model, acceleration potential is solved in the calculation of pressure on tanks and the artificial damping model is adopted to account for energy dissipation during sloshing. The Boundary Element Method (BEM) is used to solve boundary value problems of both velocity potential and acceleration potential. Numerical calculation results are compared with published results to determine the efficiency and accuracy of the numerical code. Sloshing properties in partially filled rectangular and membrane tank under translational and rotational excitations are investigated. It is found that sloshing under horizontal and rotational excitations share similar properties. The first resonant mode and excitation frequency are the dominant response frequencies. Resonant sloshing will be excited when vertical excitation lies in the instability region. For liquid tank under rotational excitation, sloshing responses including amplitude and phase are sensitive to the location of the center of rotation. Moreover, experimental tests were conducted to analyze viscous effects on sloshing and to validate the feasibility of artificial damping models. The results show that the artificial damping model with modifying wall boundary conditions has better applicability in simulating sloshing under different fill levels and excitations.

• Transactions of Materials Processing
• /
• v.7 no.3
• /
• pp.233-245
• /
• 1998

A new approach has been proposed for the incremental analysis of the nonsteady state large deformation of planar anisotropic elastic-plastic sheet forming. A mathematical brief review of a constitutive law for the incremental deformation theory has been presented from flow theory using the minimum plastic work path for elastic-plastic material. Since the material embedded coordinate system(Lagrangian quantity) is used in the proposed theory the stress integration procedure is completely objective. A new return mapping algorithm has been also developed from the general midpoint rule so as to achieve numerically large strain increment by successive control of yield function residuals. Some numerical tests for the return mapping algorithm were performed using Barlat's six component anisotropic stress potential. Performance of the proposed algorithm was shown to be good and stable for a large strain increment, For planar anisotropic sheet forming updating algorithm of planar anisotropic axes has been newly proposed. In order to show the effectiveness and validity of the present formulation earing simulation for a cylindrical cup drawing and front fender stamping analysis are performed. From the results it has been shown that the present formulation can provide a good basis for analysis for analysis of elastic-plastic sheet metal forming processes.

• Wind and Structures
• /
• v.20 no.1
• /
• pp.15-36
• /
• 2015

Characterization of wind flows over a complex terrain, especially mountain-gorge terrain (referred to as the very complex terrain with rolling mountains and deep narrow gorges), is an important issue for design and operation of long-span bridges constructed in this area. In both wind tunnel testing and numerical simulation, a transition section is often used to connect the wind tunnel floor or computational domain bottom and the boundary top of the terrain model in order to generate a smooth flow transition over the edge of the terrain model. Although the transition section plays an important role in simulation of wind field over complex terrain, an appropriate shape needs investigation. In this study, two principles for selecting an appropriate shape of boundary transition section were proposed, and a theoretical curve serving for the mountain-gorge terrain model was derived based on potential flow theory around a circular cylinder. Then a two-dimensional (2-D) simulation was used to compare the flow transition performance between the proposed curved transition section and the traditional ramp transition section in a wind tunnel. Furthermore, the wind velocity field induced by the curved transition section with an equivalent slope of $30^{\circ}$ was investigated in detail, and a parameter called the 'velocity stability factor' was defined; an analytical model for predicting the velocity stability factor was also proposed. The results show that the proposed curved transition section has a better flow transition performance compared with the traditional ramp transition section. The proposed analytical model can also adequately predict the velocity stability factor of the wind field.

• Journal of Ship and Ocean Technology
• /
• v.10 no.4
• /
• pp.1-11
• /
• 2006

The finite volume based multi-block RANS code, WAVIS developed at MOERI is applied to the numerical self-propulsion test. WAVIS uses the cell-centered finite volume method for discretization of the governing equations. The realizable $k-{\epsilon}$ turbulence model with a wall function is employed for the turbulence closure. The free surface is captured with the two-phase level set method and body forces are used to model the effects of a propeller without resolving the detail blade flow. The propeller forces are obtained using an unsteady lifting surface method based on potential flow theory. The numerical procedure followed the self-propulsion model experiment based on the 1978 ITTC performance prediction method. The self-propulsion point is obtained iteratively through balancing the propeller thrust, the ship hull resistance and towing force that is correction for Reynolds number difference between the model and full scale. The unsteady lifting surface code is also iterated until the propeller induced velocity is converged in order to obtain the propeller force. The self-propulsion characteristics such as thrust deduction, wake fraction, propeller efficiency, and hull efficiency are compared with the experimental data of the practical container ship. The present paper shows that hybrid RANS and potential flow based numerical method is promising to predict the self-propulsion parameters of practical ships as a useful tool for the hull form and propeller design.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.41 no.12
• /
• pp.931-939
• /
• 2013

In this paper, we try to find the lifting-line method which is applicable to the conceptual design of aircraft wings, and analyze the accuracy and coverage of the method. Two methods that are extended from the lifting-line theory of Prandtl are selected. One of the methods is Weissinger's method which imposes the velocity boundary condition at the control points located at the quarter chord, and the other is Phillips's method which combines the three-dimensional vortex lifting law. Calculations are performed for an elliptic wing, a swept back wing, and a tapered unswept wing with dihedral angle and geometric twist. The aerodynamic data of the potential flow such as spanwise distributions of circulation and downwash, lift and induced drag are obtained through calculations, and these data are compared with theoretical results and wind tunnel test data. As a result, Weissinger's method showed good accuracy and reliability regardless of wing shapes, but Phillips's method revealed inaccurate results for a swept back wing.

• Journal of Ocean Engineering and Technology
• /
• v.10 no.3
• /
• pp.96-104
• /
• 1996

Hamilton's principle is used to derive Euler-Lagrange equations for free surface flow problems of incompressible ideal fluid. The velocity field is chosen to satisfy the continuity equation a priori. This approach results in a hierarchial set of governing equations consist of two evolution equations with respect to two canonical variables and corresponding boundary value problems. The free surface elevation and the Lagrange's multiplier are the canonical variables in Hamilton's sense. This Lagrange's multiplier is a velocity potential defined on the free surface. Energy is conserved as a consequence of the Hamiltonian structure. These equations can be applied to waves in water of finite depth including generalization of Hamilton's equations given by Miles and Salmon.

• Geotechnical Engineering
• /
• v.12 no.2
• /
• pp.95-106
• /
• 1996

In this study, a constitutive model which can describe the anisotropic and plastic behaviors of natural cohesive soils, was developed based on anisotropic bounding surface theory. The model was fomulated by the concepts of the improved anisotropic bounding surface function, nonassociated flow rule with new plastic potential function, anisotropic hardening rule, and new mapping rule governing the plastic behavior inside bounding sutraface. Comparing with the results of Ku consolidation and triaxial shearing tests, the predictions by the proposed model agree quite well with real soil responses.