• 대한조선학회지
• /
• 제24권3호
• /
• pp.9-16
• /
• 1987

In this paper, the fluid flow in the two-dimensional tank is analyzed by the Finite Difference Method. The Navier-Stokes equation is modified for the tank fixed coordinate system. For the treatment of the free surface, the Volume of Fluid Method by Hirt and Nichols is adopted. The continuity equation and the Poisson equation which is derived from the Navier-Stokes equation to find the pressure are solved by the Successive-Line-Overrelaxation Method. The comparison of the calculated results with experimental data show a favorable agreement. The fluid flow in the two-dimensional tank can be predicted reasonably before the free surface reaches breaking by this numerical method.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 2006년 제4회 한국유체공학학술대회 논문집
• /
• pp.349-352
• /
• 2006

The conventional segregated finite element formulation produces a small and simple matrix at each step than in an integrated formulation. And the memory and cost requirements of computations are significantly reduced because the pressure equation for the mass conservation of the Navier-Stokes equations is constructed only once if the mesh is fixed. However, segregated finite element formulation solves Poisson equation of elliptic type so that it always needs a pressure boundary condition along a boundary even when physical information on pressure is not provided. On the other hand, the conventional integrated finite element formulation in which the governing equations are simultaneously treated has an advantage over a segregated formulation in the sense that it can give a more robust convergence behavior because all variables are implicitly combined. Further it needs a very small number of iterations to achieve convergence. However, the saddle-paint-type matrix (SPTM) in the integrated formulation is assembled and preconditioned every time step, so that it needs a large memory and computing time. Therefore, we newly proposed the P2PI semi-segregation formulation. In order to utilize the fact that the pressure equation is assembled and preconditioned only once in the segregated finite element formulation, a fixed symmetric SPTM has been obtained for the continuity constraint of the present semi-segregation finite element formulation. The momentum equation in the semi-segregation finite element formulation will be separated from the continuity equation so that the saddle-point-type matrix is assembled and preconditioned only once during the whole computation as long as the mesh does not change. For a comparison of the CPU time, accuracy and condition number between the two methods, they have been applied to the well-known benchmark problem. It is shown that the newly proposed semi-segregation finite element formulation performs better than the conventional integrated finite element formulation in terms of the computation time.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제10권3호
• /
• pp.329-347
• /
• 2018

The Smoothed Particle Hydrodynamics (SPH) method has proved to have great potentials in dealing with the wave-structure interactions. Compared with the Weakly Compressible SPH (WCSPH) method, the ISPH approach solves the pressure by using the pressure Poisson equation rather than the equation of state. This could provide a more stable and accurate pressure field that is important in the study of wave-structure interactions. This paper improves the solid boundary treatment of ISPH by using a high accuracy Simplified Finite Difference Interpolation (SFDI) scheme for the 2D wave-structure coupling problems, especially for free-moving structure. The proposed method is referred as the ISPH_BS. The model improvement is demonstrated by the documented benchmark tests and laboratory experiment covering various wave-structure interaction applications.

• Advances in Computational Design
• /
• 제9권1호
• /
• pp.39-52
• /
• 2024

If the governing differential equation arising from engineering problems is treated as an analytic, continuous and derivable function, it can be expanded by one point as a series of finite numbers. For the function to be zero for each value of its domain, the coefficients of each term of the same power must be zero. This results in a recursive relationship which, after applying the natural conditions or the boundary conditions, makes it possible to obtain the values of the derivatives of the function with acceptable accuracy. The elastoplastic analysis of an inhomogeneous thick sphere of metallic materials with linear variation of the modulus of elasticity, yield stress and Poisson's ratio as a function of radius subjected to internal pressure is presented. The Beltrami-Michell equation is established by combining equilibrium, compatibility and constitutive equations. Assuming axisymmetric conditions, the spherical coordinate parameters can be used as principal stress axes. Since there is no analytical solution, the natural boundary conditions are applied and the governing equations are solved using a proposed new method. The maximum effective stress of the von Mises yield criterion occurs at the inner surface; therefore, the negative sign of the linear yield stress gradation parameter should be considered to calculate the optimal yield pressure. The numerical examples are performed and the plots of the numerical results are presented. The validation of the numerical results is observed by modeling the elastoplastic heterogeneous thick sphere as a pressurized multilayer composite reservoir in Abaqus software. The subroutine USDFLD was additionally written to model the continuous gradation of the material.

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

In this paper, incompressible two-dimensional Navier-Stokes equations are numerically solved for the study of steady laminar flow around a body with the wall effect. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. To investigate the wall effect, numerical computations are carried out for the NACA 0012 section at the various blockage ratios. The pressure and skin friction on the foil surface, velocity pronto in its wake and drag coefficient are investigated as functions of the blockage ratio.

• 대한조선학회논문집
• /
• 제33권2호
• /
• pp.13-29
• /
• 1996

이 논문은 자유표면을 포함하는 시리즈 60($C_B=0.6$)선형 주위 유동장에 대한 계산결과를 보여준다. 지배 방정식으로는 3차원 Navier-Stokes 방정식을 사용하고, 높은 레이놀즈수에서의 난류 유동장을 계산하기 위하여 Baldwin-Lomax난류모형을 채용하였으며 계산시간을 줄이기 위해 물체 표면에서는 벽법칙을 채용하였다. 지배 방정식은 유한 차분법에 의해 차분화 되었으며, 음해법[1]에 의해, 압력 Poisson방정식은 완화법(successive-over-relaxation method)에 의해 프로그램을 구성하였다. 자유표면 유동을 정확히 계산하기 위해서는, 동역학적 자유표면 경계조건식의 수치해법이 매우 중요하다. 이 논문에서는 세 가지의 수치해법을 채용하여 그 결과를 실험결과와 비교하였다. 결론적으로, 계산된 저항계수($C_F,\;C_P$와 파형은 실험 값과 잘 일치하고 있다.

• 대한조선학회논문집
• /
• 제51권4호
• /
• pp.328-333
• /
• 2014

A hybrid particle-mesh method based on the vorticity-velocity formulation for solving the incompressible Navier-Stokes equations is a combination of the Vortex-In-Cell(VIC) method for convection and the penalization method for diffusion. The key feature of the numerical methods is to determine velocity and vorticity fields around a solid body on a temporary grid, and then the time evolution of the flow is computed by tracing the convection of each vortex element using the Lagrangian approach. Assuming that the vorticity and velocity fields are to be computed in time domain analysis, pressure fields are estimated through a complete set of solutions at present time step. It is possible to obtain vorticity and velocity fields prior to any pressure calculation since the pressure term is eliminated in the vorticity-velocity formulation. Therefore, pressure field is explicitly treated by solving a suitable Poisson equation. In this paper, we propose a simple way to numerically implement the vorticity-velocity-pressure formulation including a penalty term. For validation of the proposed numerical scheme, we illustrate the early development of viscous flows around an impulsive started circular cylinder for Reynolds number of 9500.

• 한국가시화정보학회지
• /
• 제21권2호
• /
• pp.14-23
• /
• 2023

The pressure difference across stenotic blood vessels is a commonly used clinical metric for diagnosing many cardiovascular diseases. At present, most clinical pressure measurements rely solely on invasive catheterization. In this study, we propose a novel method for non-invasive pressure estimation using the incompressible Navier-Stokes equations and a 3D multi-path integration approach. We verify spatio-temporal convergence on an in-silico dataset of a cylindrical straight pipe phantom with steady and pulsatile flow fields. We then evaluate the proposed method on an in vitro dataset of reconstructed control, pre-operative, and post-operative carotid artery cases acquired from 4D flow MRI. The performance of our method is compared to existing approaches based on the pressure Poisson equation and work-energy relative pressure. The results demonstrate the proposed method's high accuracy, robustness to spatio-temporal subsampling, and reduced sensitivity to noise, highlighting its great potential for non-invasive pressure estimation.

• 대한조선학회논문집
• /
• 제35권4호
• /
• pp.1-10
• /
• 1998

본 논문에서는 비압축성 Newtonian 점성유동에서 초기에 순간 출발하는 2차원 실린더 주위의 유동을 해석하기 위해서, 와도를 기저로 한 수치해석기법을 제안하고 있다. Helmholtz 분리 형태로 표현된 Navier-Stokes방정식에서 유도되는 와도전달방정식과 압력방정식, 그리고 벡터등식에서 유도되는 속도-와도 관계식을 이 문제의 지배방정식으로 택하고, 경계조건으로는 물체표면에서 와도와 압력의 연성관계와 힘의 평형을 고려한 동적와도경계조건과 동적압력조건이 제시된다. 이 지배방정식과 경계조건을 수치적으로 처리하기 위하여, 와도와 압력이 연성되어 있는 경계조건은 Wu등(1994)이 제안한 대로, 연성관계를 유지한 채로 식을 분리하는 방법을 이용하였고, 와도전달 방정식은 유한체적법으로 계산하였으며, 그 식에 포함된 대류항을 처리하는 방법으로 TVD 방법을 이용하였다. 속도는 Biot-Savart적분항이 포함된 벡터등식에서 panel방법으로 구하고, 압력방정식은 형태가 Poisson방정식이므로 역시 panel방법을 이용하였다. 계산에 사용된 격자로 정규격자를 이용하고, 결과를 다른 수치적, 해석적 결과와 비교하여 그 타당성을 검증하였다.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 2003년도 추계학술대회 논문집 전기물성,응용부문
• /
• pp.220-222
• /
• 2003

In this paper, we analyze on the discharge characteristics through 1-D simulations in an at plasma display panel discharge cell. The model is based on a Poisson' equation, continuity and drift-diffusion equation. Results are presented in a 95% neon, 5% xenon gas mixture, for a gap length of 100us and a gas pressure of 400Torr at ambient temperature. Results for other gap length are also discussed. As a result, an increase of the gap cause increase of luminous efficiency but with larger sustaining voltage.