• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.1
• /
• pp.1-8
• /
• 1996

A potential-based panel method has been developed for numerical computation of wave resistance on a hybrid hydrofoil. Hybrid hydrofoil is composed of a main body, two struts and two hydrofoils. The main body, which is assumed to be an axisymmetric body for the present analysis, is normally used to support displacement of a body with its buoyancy. Normal dipoles and the sources are distributed on the body(main body, struts, hydrofoils) and the sources are distributed on the free surface. Linearized free surface and the radiation conditions are satisfied using the fourth order finite difference operator and the semi-linear pressure Kutta condition is used for the numerical computation of the hydrofoils. Poisson type free surface condition has been used for the numerical computation and hyperboloidal panel method has been used for better numerical accuracy. To verify this numeric method, model tests are performed in circulation water channel. From the comparison of experimental results with numeric ones, the present method can be used as a useful tool for the design of high speed vessels.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2001.10a
• /
• pp.16-20
• /
• 2001

종래에 자동차 차체의 프레스 성형 공정은 설계자의 정성적인 계산이나 다이 (Trial Dic)에 의해 시행 착오 끝에 설계되었다. 그런데, 부정확하게 계산되는 경우가 많아, 아이가 준비된 후 단점이 노출되어 여러 차례 수정 과정을 겪게 되어 시간과 경비를 증대시키는 원인이 되었다. 자동차 회사에서는 이 때문에 금형 설계 단계에서 빠르고 정확하게 금형의 불량을 예측하고 성형성을 정량적으로 평가하기 위한 시스템을 필요로 한다. 그중 비선형 유한요소법(F.E.M)에 의한 시뮬레이션은 정확한 해와 많은 정보를 줄 수 있다는 장점을 가지고 있다. 그러나, 3차원적인 대단히 복잡하고도 치수가 큰 자동차 판넬에 있어서는 수렴성 등의 문제로 적절한 시간 안에 충분한 정밀도를 갖는 해를 구할 수 있는 단계까지는 이르지 못하고 있다.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.12
• /
• pp.35-44
• /
• 2013

In this study, a particle method without using grid was applied for analysing large deformation problems in soil flows instead of using ordinary finite element or finite difference methods. In the particle method, a continuum equation was discretized by various particle interaction models corresponding to differential operators such as gradient, divergence, and Laplacian. Soil behavior changes from solid to liquid state with increasing water content or external load. The Mohr-Coulomb failure criterion was incorporated into the particle method to analyze such three-dimensional soil behavior. The yielding and hardening behavior of soil before failure was analyzed by treating soil as a viscous liquid. First of all, a sand column test without confining pressure and strength was carried out and then a self-standing clay column test with cohesion was carried out. Large deformation from such column tests due to soil yielding or failure was used for verifying the developed particle method. The developed particle method was able to simulate the three-dimensional plastic deformation of soils due to yielding before failure and calculate the variation of normal and shear stresses both in sand and clay columns.

• Geophysics and Geophysical Exploration
• /
• v.18 no.4
• /
• pp.223-231
• /
• 2015

Self-potential (SP) is sensitive to groundwater flow and there are many causes to generate SP. Among many mechanisms of SP, pore-fluid flow in porous media can generate potential without any external current source, which is referred to as electrokinetic potential or streaming potential. When calculating SP responses on the surface due to geothermal fluid within an engineered geothermal system (EGS) reservoir, SP anomaly is usually considered to be generated by fluid injection or production within the reservoir. However, SP anomaly can also result from geothermal water fluid within EGS reservoirs experiencing temperature changes between injection and production wells. For more precise simulation of SP responses, we developed an algorithm being able to take account of SP anomalies produced by not only water injection and production but also the fluid of geothermal water, based on three-dimensional finite-element-method employing tetrahedron elements; the developed algorithm can simulate electrical potential responses by both point source and volume source. After verifying the developed algorithm, we assumed a simple geothermal reservoir model and analyzed SP responses caused by geothermal water injection and production. We are going to further analyze SP responses for geothermal water in the presence of water production and injection, considering temperature distribution and geothermal water flow in the following research.

• The Korean Journal of Rheology
• /
• v.10 no.4
• /
• pp.217-226
• /
• 1998

The anisotropy caused by the fiber orientation, which is inevitably generated by the flow during injection molding of short fiber reinforced polymers, greatly influences dimensional accuracy, mechanical properties, and other quality of the final product. Since the filling stage of the injection molding process plays a vital role in determining fiber orientation, an accurate analysis of flow field for the filling stage is needed. Unbalanced filling occurs when a complex or a multi-cavity mold is used leading to development of regions where the fiber suspension is under compression. It is impossible to make an accurate calculation of the flow field during filling with the analysis assuming incompressible fluid. A mold with four cavities with different filling times was produced to compare the numerical analysis results with the experimental data. There was a good agreement between the experimental and theoretical results when the compressibility of the polymer melt was considered for the numerical simulation. The fiber orientation states for compressible and incompressible fluids were also compared qualitatively as well as quantitatively in this study.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.14 no.2
• /
• pp.151-160
• /
• 2002

A modeling system for three-dimensional flow (FEMOS) has been developed and applied to simulate the tidal currents of Asan Bay. The system can consider tidal flats changing with time and uses a finite element method that can adapt coastline change effectively. The simulation results for Asan Bay with large tidal flats, shallow water depth and high tidal range showed good agreements with the observed currents of long-term variations at the medium layer and short-term variations of vertical profiles. Based on the simulated tidal currents, the horizontal distributions of bottom shear stress were calculated and showed close relation with the change of bottom topography. The system can be used widely to study coastal circulation in the coastal region with complex geography.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.568-576
• /
• 2005

A 3-D finite element model combined with a volume tracking method is presented in this work to simulate the mold filling for casting processes. Especially, the analysis involves an adaptive grid method that is created under a criterion of element categorization of filling states and locations in the total region at each time step. By using an adaptive grid wherein the elements, finer than those in internal and external regions, are distributed at the surface region through refinement and coarsening procedures, a more efficient analysis of transient fluid flow with free surface is achieved. Adaptive grid based on VOF method is developed in tetrahedral element system. Through a 3-D analysis of the benchmark test of the casting process, the efficiency of the proposed adaptive grid method is verified. Developed FE code is applied to a typical industrial part of the casting process such as aluminum road wheel.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.540-540
• /
• 2012

강우 발생시 유역에 집수된 물이 하천에 이르는 경로에 따른 유출은 지표유출, 지표하유출, 지하수유출로 구분된다. 정확한 수문순환 과정의 해석을 위해서는 지표 흐름뿐만 아니라 지표하 및 지하수 흐름의 해석이 중요한 실정이나 일반적으로 실무에서 사용되는 강우-유출해석 모형은 지표유출을 해석하기 위한 모형이 대부분이며, 지표하 유출과 침투량을 산정하는데 어려움이 있다. 일반적인 강우-유출해석 모형은 Horton 방법, NRCS 방법, Green-Ampt 방법에 의해 유효우량을 분리하며, 이 과정은 침투량을 직접적으로 모형화 할 수 없으므로 지표 및 지표하, 지하수 흐름을 복합적으로 해석할 수 있는 모형이 질적이나 양적으로 부족한 실정이다. 이러한 지표하 흐름과 침투량을 산정하기 위하여 FE-FLOW, PM, MS-VMS, GMS, GW-VISTAS, ARGUS 및 MODFLOW와 같은 지하수 모형을 사용하고 있다. 본 연구에서는 지표하유출 산정을 위한 침투량의 비교분석을 위해 현재 가장 범용되는 지하수 유동 모델링 프로그램인 Visual Modlfow 모형과 GMS 모형을 이용하여 침투량 산정을 위한 수치 모의를 진행하였다. 각 모형의 입력자료는 2009년 국립방재연구원에서 수행한 침투실험시설 자료를 이용하여 동일한 조건을 부여하고, 두 모형의 비교를 위해 Visual Modflow에서는 MODFLOW의 기본 해석방법인 유한차분법(FDM)을 이용하고, GMS 모형에서는 3차원 유한요소해석이 가능한 GMS-FEMWATER를 이용하였다. 두 모형의 수치모의 조건으로 2009년 국립방재연구원에서 수행한 침투실험방법과 동일하게 공극률에 따른 투수성 보도블럭의 구분과 50mm/hr, 100mm/hr, 150mm/hr, 200mm/hr의 강우강도별 선행함수조건에 따른 수치모의를 진행하였으며, 수치모의된 침투량의 적정성을 판단하기 위하여 국립방재연구원의 침투실험 결과자료와 비교분석하였다. 침투실험 자료와 각각 수치모의된 침투량을 비교분석한 결과 서로 유사한 경향을 보이고 있으나 초기 침투시 상대오차가 비교적 크게 발생하였다. 이는 수치모형의 경우 수리실험과는 다르게 모의시작과 동시에 해당 강우강도의 침투가 시작되므로 초기 유입 유출량 발생시간의 차이가 종료시간까지 누적 침투량에 미치는 것으로 판단되며, 매개변수에 많은 영향을 받는 것으로 판단된다.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.3
• /
• pp.35-47
• /
• 1996

A surface panel method treating a boundary-value problem of the Dirichlet type is presented to design a three dimensional body with free surface corresponding to a prescribed pressure distribution. An integral equation is derived from Green's theorem, giving a relation between total potential of known strength and the unknown local flux. Upon discretization, a system of linear simultaneous equations is formed including free surface boundary condition and is solved for an assumed geometry. The pseudo local flux, present due to the incorrect positioning of the assumed geometry, plays a role f the geometry corrector, with which the new geometry is computed for the next iteration. Sample designs for submerged spheroids and Wigley hull and carried out to demonstrate the stable convergence, the effectiveness and the robustness of the method. For the calculation of the wave resistance, normal dipoles and Rankine sources are distributed on the body surface and Rankine sources on the free surface. The free surface boundary condition is linearized with respect to the oncoming flow. Four-points upwind finite difference scheme is used to compute the free surface boundary condition. A hyperboloidal panel is adopted to represent the hull surface, which can compensate the defects of the low-order panel method. The design of a 5500TEU container carrier is performed with respect to reduction of the wave resistance. To reduce the wave resistance, calculated pressure on the hull surface is modified to have the lower fluctuation, and is applied as a Dirichlet type dynamic boundary condition on the hull surface. The designed hull form is verified to have the lower wave resistance than the initial one not only by computation but by experiment.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.97-106
• /
• 1989

This paper is on the modeling of two-dimensional groundwater flow, which is the first step of the development of Dynamic System Model for groundwater flow and pollutant transport in subsurface porous media. The particular features of the model are its versatility and flexibility to deal with as many real-world problems as possible. Points as well as distributed sources/sinks are included to represent recharges/pumping and rainfall infiltrations. All sources/sinks can be transient or steady state. Prescribed hydraulic head on the Dirichlet boundaries and fluxes on Neumann or Cauchy boundaries can be time-dependent or constant. Sources/sinks strength over each element and node, hydraulic head at each Dirichlet boundary node and flux at each boundary segment can vary independently of each other. Either completely confined or completely unconfined aquifers, or partially confined and partially unconfined aquifers can be dealt with effectively. Discretization of a compound region with very irregular curved boundaries is made easy by including both quadrilateral and triangular elements in the formulation. Large-field problems can be solved efficiently by including a pointwise iterative solution strategy as an optional alternative to the direct elimination solution methed for the matrix equation approximating the partial differential equation of groundwater flow. The model also includes transient flow through confining leaky aquifers lying above and/or below the aquifer of interest. The model is verified against three simple cases to which analytical solutions are available. The groundwater flow model shall be combined with the model of pollutant transport in subsurface porous media. Then the combined model, with the applications of the Eigenvalue technique and the Dynamic system theory, shall be improved to the Dynamic System Model which can simulate the real groundwater flow and the pollutant transport accurately and effectively for the analyses and predictions.