• International Journal of Naval Architecture and Ocean Engineering
• /
• 제5권4호
• /
• pp.546-558
• /
• 2013

The research performed in this paper was carried out to investigate the numerical analysis of the sheet cavitation on marine propeller. The method is boundary element method (BEM). Using the Green's theorem, the velocity potential is expressed as an integral equation on the surface of the propeller by hyperboloid-shaped elements. Employing the boundary conditions, the potential is determined via solving the resulting system of equations. For the case study, a DTMB4119 propeller is analyzed with and without cavitating conditions. The pressure distribution and hydrodynamic performance curves of the propellers as well as cavity thickness obtained by numerical method are calculated and compared by the experimental results. Specifically in this article cavitation changes are investigate in both the radial and chord direction. Thus, cross flow variation has been studied in the formation and growth of sheet cavitation. According to the data obtained it can be seen that there is a better agreement and less error between the numerical results gained from the present method and Fluent results than Hong Sun method. This confirms the accurate estimation of the detachment point and the cavity change in radial direction.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제18권1호
• /
• pp.27-41
• /
• 2014

In this paper, we present an efficient numerical method for multiphase image segmentation using a multiphase-field model. The method combines the vector-valued Allen-Cahn phase-field equation with initial data fitting terms containing prescribed interface width and fidelity constants. An efficient numerical solution is achieved using the recently developed hybrid operator splitting method for the vector-valued Allen-Cahn phase-field equation. We split the modified vector-valued Allen-Cahn equation into a nonlinear equation and a linear diffusion equation with a source term. The linear diffusion equation is discretized using an implicit scheme and the resulting implicit discrete system of equations is solved by a multigrid method. The nonlinear equation is solved semi-analytically using a closed-form solution. And by treating the source term of the linear diffusion equation explicitly, we solve the modified vector-valued Allen-Cahn equation in a decoupled way. By decoupling the governing equation, we can speed up the segmentation process with multiple phases. We perform some characteristic numerical experiments for multiphase image segmentation.

• 대한조선학회논문집
• /
• 제53권1호
• /
• pp.1-9
• /
• 2016

It is important to research and understand the physical phenomenon around a semi-submersible offshore structure on waves and currents because the wave run-up and load occurs owing to the waves and currents. In this study, the numerical simulations are performed about flow around a fixed semi-submersible offshore structure. The Modified Marker-density method is adopted in the present computation procedure, this method is one of the various methods to define the free-surface. The present computation results are compared with existing experimental and numerical simulation(VOF method) results. And, the computation results are relatively coincident with the existing results of model test and numerical simulation by VOF method.

• 대한조선학회논문집
• /
• 제46권3호
• /
• pp.232-248
• /
• 2009

Numerical methodology to solve ship springing problem, which is basically fluid-structure interaction problem, was explored in this study. Solution of this hydroelasticity problem was sought by coupling higher order B-spline Rankine panel method and finite element method in time domain, each of which is introduced for fluid and structure domain respectively. Even though varieties of different combinations in terms of numerical scheme are possible and have been tried by many researchers to solve the problem, no systematic study regarding the characteristics of each scheme has been done so far. Here, extensive case studies have been done on the numerical schemes especially focusing on the iteration method, FE analysis of beam-like structure, handling of forward speed problem and so on. Two different iteration scheme, Newton style one and fixed point iteration, were tried in this study and results were compared between the two. For the solution of the FE-based equation of motion, direct integration and modal superposition method were compared with each other from the viewpoint of its efficiency and accuracy. Finally, calculation of second derivative of basis potential, which is difficult to obtain with accuracy within grid-based method like BEM was discussed.

• 한국전산구조공학회논문집
• /
• 제23권4호
• /
• pp.369-378
• /
• 2010

포화된 다공질매체의 수치해석에서는 일반적으로 고체영역과 유체영역을 동시에 고려한 혼합유한요소해석(Mixed Finite Element Analysis)이 쓰인다. 여기서 고체영역과 유체영역에서의 변수를 계산하기 위해서는 직접법(Direct Method) 또는 반복법(Iterative method)을 사용할 수 있으나, 각 구성물질의 상이한 물리적 특성 때문에 수치안정성을 확보하기 위해서는 대부분 스태거드 방법(Staggered method)이 제안된다. 본 논문에서는 수치안정성을 높인 스태거드 방법에서 영역 분할기법 중 하나인 FETI(Finite Element Tearing and Interconnecting)기법을 고체영역에 접목시켜 수치효율성을 증대시키는 방법이 제안되었다. 고체영역에서 라그랑지 승수와 Conjugated Gradient Method를 이용해 영역 분할이 진행되고 MPI(Message Passing Interface) 라이브러리를 사용하여 수치 효율성을 검증하였다.

• Geomechanics and Engineering
• /
• 제14권3호
• /
• pp.271-281
• /
• 2018

A jointed rock slope stability evaluation was simulated by a discontinuous deformation analysis numerical method to investigate the process and safety factors for different crack distributions and different overloading situations. An optimized method using Discontinuous Deformation Analysis for Rock Failure (DDARF) is presented to perform numerical investigations on the jointed rock slope stability evaluation of the Dagangshan hydropower station. During the pre-processing of establishing the numerical model, an integrated software system including AutoCAD, Screen Capture, and Excel is adopted to facilitate the implementation of the numerical model with random joint network. These optimizations during the pre-processing stage of DDARF can remarkably improve the simulation efficiency, making it possible for complex model calculation. In the numerical investigations on the jointed rock slope stability evaluations using the optimized DDARF, three calculation schemes have been taken into account in the numerical model: (I) no joint; (II) two sets of regular parallel joints; and (III) multiple sets of random joints. This model is capable of replicating the entire processes including crack initiation, propagation, formation of shear zones, and local failures, and thus is able to provide constructive suggestions to supporting schemes for the slope. Meanwhile, the overloading numerical simulations under the same three schemes have also been performed. Overloading safety factors of the three schemes are 5.68, 2.42 and 1.39, respectively, which are obtained by analyzing the displacement evolutions of key monitoring points during overloading.

• 한국해양공학회지
• /
• 제26권6호
• /
• pp.7-13
• /
• 2012

This paper presents a numerical study of the wave loads acting on offshore structures using a Cartesian-grid-based flow simulation method. Finite volume discretization with a volume-of-fluid (VOF) method is adopted to solve two-phase Navier-Stokes equations. Among the many variations of the VOF method, the CICSAM scheme is applied. The body boundary conditions are satisfied using a porosity function, and wave generation is carried out by using transient (wave or damping) zone approaches. In order to validate the present numerical method, three different basic offshore structures, including a sphere, Pinkster barge, and Wigley model, are numerically investigated. First, diffraction and radiation problems are solved using the present numerical method. The wave exciting and drift forces from the diffraction problems are compared with potential-based solutions. The added mass and wave damping forces from the radiation problems are also compared with the potential results. Next, the wave-induced motion responses of the structures are calculated and compared with the existing experimental data. The comparison results are fairly good, showing the validity of the present numerical method.

• 한국지반환경공학회 논문집
• /
• 제23권12호
• /
• pp.17-23
• /
• 2022

본 논문은 Top-down 방식과 Bottom-up 방식으로 시공된 프리캐스트 판넬식 옹벽의 안전율 변화를 PLAXIS 2D 소프트웨어를 사용하여 수치해석하였다. 수치해석 조건은 시공방식, 네일 및 판넬의 고정방법, 뒤채움재 다짐상태, 강우조건 등이 있다. Top-down 방식과 Bottom-up 방식의 차이점은 뒤채움재의 유무와 지반굴착의 방식으로, Top-down 방식은 지반을 수직으로 굴착한 뒤 시공하나 Bottom-up 방식은 지반을 경사를 주어 굴착한 뒤 모르타르 및 뒤채움재를 이용하여 판넬을 부착한다. 수치해석 결과 Top-down 방식은 Bottom-up 방식과 비교하였을 때, 우기를 제외한 모든 경우에 더 높은 안전률을 확보하였다. Top-down 방식의 우기 시 안전율 저하는 Top-down 방식으로 시공 시 프리캐스트 판넬을 부착하기 위해 사용한 모르타르가 주된 원인이 된 것으로 판단하였다.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2004년도 봄 학술발표회 논문집
• /
• pp.333-340
• /
• 2004

This paper deals with geometric nonlinear analyses using a new meshfree technique which improves the numerical integration accuracy. The new method called the Petrov-Galerkin natural element method (PGNEM) is based on the Voronoi diagram and the Delaunay triangulation which is based on the same concept used for conventional natural element method called the Bubnov-Galerkin natural element method (BGNEM). But, unlike BGNEM, the test shape function is differently chosen from the trial shape function. In the linear static analysis, it is ensured that the numerical integration error of the PGNEM is remarkably reduced. In this paper, the PGNEM is applied to large deformation problems, and the accuracy of the proposed numerical technique is verified through the several examples.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 추계학술대회
• /
• pp.474-479
• /
• 2004

According to our previous study, it is confirmed that the Petrov-Galerkin natural element method (PGNEM) completely resolves the numerical integration inaccuracy in the conventional Bubnov-Galerkin natural element method (BG-NEM). This paper is an extension of PG-NEM to two-dimensional nonlinear dynamic problem. For the analysis, a constant average acceleration method and a linearized total Lagrangian formulation is introduced with the PG-NEM. At every time step, the grid points are updated and the shape functions are reproduced from the relocated nodal distribution. This process enables the PG-NEM to provide more accurate and robust approximations. The representative numerical experiments performed by the test Fortran program, and the numerical results confirmed that the PG-NEM effectively and accurately approximates the nonlinear dynamic problem.