• 한국전산유체공학회지
• /
• 제19권4호
• /
• pp.20-28
• /
• 2014

An efficient solution algorithm for simulating free surface problem is presented. Navier-Stokes equations for variable density incompressible flow are employed as the governing equation on Cartesian meshes. In order to describe the free surface motion efficiently, VOF(Volume Of Fluid) method utilizing THINC(Tangent of Hyperbola for Interface Capturing) scheme is employed. The most time-consuming part of the current free surface flow simulations is the solution step of the linear system, derived by the pressure Poisson equation. To solve a pressure Poisson equation efficiently, the PCG(Preconditioned Conjugate Gradient) method is utilized. This study showed that the proper application of the preconditioner is the key for the efficient solution of the free surface flow when its pressure Poisson equation is solved by the CG method. To demonstrate the efficiency of the current approach, we compared the convergence histories of different algorithms for solving the pressure Poisson equation.

• 대한조선학회논문집
• /
• 제32권3호
• /
• pp.62-71
• /
• 1995

비압축성 가정하에 Navier-Stokes 방정식을 이용하여 비정상 점성유동을 수치해석하기 위해서는 매시간 단계에서 타원형 압력 Poisson 방정식의 해를 구해야 하며, 이에 많은 계산시간이 소요된다. 본 논문에서는 직접법을 이용하여 압력 Poisson 방정식을 수치해석하였으며, 분할수 증가에 따른 소요시간 문제를 다루었다. Green 정리를 압력 Poisson 방정식에 적용하면 주어진 문제는 경계치문제로 변환되고, convolution 형의 영역적분은 F.F.T.를 이용하여 계산시간을 단축할 수 있어, 직접법 이용시 소요시간은 경계치문제의 해를 구하는 데에 좌우된다. 직접법의 검증을 위하여 해석해를 알고 있는 경우에 대하여 수치해석하였고, 물체경계조건과 정합문제에 관하여 수치해석 하였는데. 분할수가 (n.n) 시 O($n^{3}$) 미만의 계산시간으로 수치해석할 수 있었다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 추계 학술대회논문집
• /
• pp.96-101
• /
• 1998

Various discretiation methods of Laplacian operator in the Pressure-Poisson equation are investigated for the solution of incompressible Navier-Stokes equations using the non-staggered grid. Laplacian operators previously proposed by other researchers are applied to a Driven-Cavity problem. The computational results are compared with those of Ghia. The results show the characteristics of the discrete Laplacian operators.

• 대한조선학회논문집
• /
• 제33권4호
• /
• pp.48-59
• /
• 1996

비정상 점성유동 수치해석단계는 연속방정식을 만족시키는 공간해석단계와 시간전진단계로 구분할 수 있다. 본 연구에서는 공간해석단계의 압력 Poisson 방정식의 해를 구하는데 스팩트럴법을 이용하였다. 압력 Poisson 방정식의 최고차미분항을 Fourier 급수로 전개하면 압력 및 압력의 1차미분항의 Fourier 급수의 적분으로 표현되므로 Gibb's 현상을 제거할 수 있어, 비주기성인 경우에도 스팩트럴법을 적용할 수 있다. 수치해법의 검증을 위하여 2차원 원주상체 및 날개주위 비정상 점성유동을 수치해석하였고, 그 결과를 비교하여 보았다.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제9권4호
• /
• pp.390-403
• /
• 2017

Smoothed Particle Hydrodynamics (SPH) method has a good adaptability for the simulation of free surface flow problems. There are two forms of SPH. One is weak compressible SPH and the other one is incompressible SPH (ISPH). Compared with the former one, ISPH method performs better in many cases. ISPH based on Rankine source solution can perform better than traditional ISPH, as it can use larger stepping length by avoiding the second order derivative in pressure Poisson equation. However, ISPH_R method needs to solve the sparse linear matrix for pressure Poisson equation, which is one of the most expensive parts during one time stepping calculation. Iterative methods are normally used for solving Poisson equation with large particle numbers. However, there are many iterative methods available and the question for using which one is still open. In this paper, three iterative methods, CGS, Bi-CGstab and GMRES are compared, which are suitable and typical for large unsymmetrical sparse matrix solutions. According to the numerical tests on different cases, still water test, dam breaking, violent tank sloshing, solitary wave slamming, the GMRES method is more efficient than CGS and Bi-CGstab for ISPH method.

• Ocean Systems Engineering
• /
• 제1권3호
• /
• pp.207-222
• /
• 2011

A stabilized incompressible smoothed particles hydrodynamics (ISPH) method is utilized to simulate free falling rigid body into water domain. Both of rigid body and fluid domain are modeled by SPH formulation. The proposed source term in the pressure Poisson equation contains two terms; divergence of velocity and density invariance. The density invariance term is multiplied by a relaxed parameter for stabilization. In addition, large eddy simulation with Smagorinsky model has been introduced to include the eddy viscosity effect. The improved method is applied to simulate both of free falling vessels with different materials and water entry-exit of horizontal circular cylinder. The applicability and efficiency of improved method is tested by the comparisons with reference experimental results.

• 대한기계학회논문집B
• /
• 제27권4호
• /
• pp.458-465
• /
• 2003

The main objective of the research is to develop a research code solving transient incompressible Navier-Stokes equation. In this research code, Adams-Bashforth method was applied to the convective terms of the navier stokes equation and the splitted equations were discretized spatially by finite element methods to solve the complex geometry problems easily. To reduce the divergence on the boundaries of pressure poisson equation due to the unsuitable pressure boundary conditions, multi step approximation pressure boundary conditions derived from the boundary linear momentum equations were used. Simulations of Lid Driven Flow and Flow over Cylinder were conducted to prove the accuracy by means of the comparison with results of the previous workers.

• 대한조선학회논문집
• /
• 제31권3호
• /
• pp.87-99
• /
• 1994

Reynolds Averaged Navier-Stokes 방정식을 수치해석하여 Wigley 선형 주위의 난류유동을 계산하였다. 정규격자상에서 공간의 이산화는 2차 정도의 유한차분법을, 시간의 적분에는 4단계 Runge-kutta법을 이용하였다. 시간 증분을 크게 하기 위하여 Jameson의 잔류항 평균 기법을 사용하였다. 압력 Poisson 방정식으로부터 압력장을 구하였고 난류닫음 조건을 만족시키기 위하여 Baldwin-Lomax의 난류 모형을 사용하였다. 수치계산은 레이놀드수가 $4.5{\times}10^6$에서 수행하였고 계산된 속도와 압력분포는 실험 결과와 비교적 잘 일치하는 것을 확인하였다.

• Interaction and multiscale mechanics
• /
• 제3권2호
• /
• pp.157-171
• /
• 2010

In this paper, a theoretical method has been developed for the electric double layer interaction under condition of the variable dielectric permittivity of water. Using Poisson-Boltzmann equation (PBE), for one plate and two plates having similar or dissimilar constant charge or constant potential, we have investigated the electric double layer potential, its gradient and the disjoining pressure as well as the effect of variation of dielectric permittivity on these parameters. It has been assumed that plates are separated by a specific distance and contain a liquid solution in between. It is shown that reduction of the dielectric permittivity near the interfaces results in compression of electric double layers and affects the potential and its gradient which leads to a decreased electrostatic repulsion. In addition, it is shown that variation of dielectric permittivity in the case of higher electrolyte concentration, leads to a greater change in potential distribution between two plates.

• 대한조선학회논문집
• /
• 제58권5호
• /
• pp.271-280
• /
• 2021

The code developed using a pressure-based method for unified conservation laws of incompressible/compressible fluids is expanded to handle moving or deforming body boundaries using the hybrid Cartesian/immersed boundary method. An instantaneous pressure field is calculated from a pressure Poisson equation for the whole fluid domain, including the compressible gas region. The polytropic gas is assumed for the compressible fluid so that the energy equation is decoupled. Immersed boundary nodes are identified based on edges crossing body boundaries. The velocity vector is reconstructed at the immersed boundary node using an interpolation along the assigned local normal line. The developed code is validated by comparing the time histories of pressure and wave elevation for sloshing in a rectangular and a membrane-type tank. The validated code is applied to simulate air cushion effects in a rectangular tank under sway motion. Time variations of pressure fields are analyzed in detail as the air pocket pulsates. It is shown that the contraction and expansion of the air pocket dominate the pressure loads on the wall of the tank. The present results are in good agreement with other experimental and computational results for the amplitude and the decay of the pressure oscillations measured at the pressure gauges.