• Title/Summary/Keyword: Pressure Poisson Equation

Search Result 51, Processing Time 0.017 seconds

AN EFFICIENT ALGORITHM FOR INCOMPRESSIBLE FREE SURFACE FLOW ON CARTESIAN MESHES (직교격자상에서 효율적인 비압축성 자유표면유동 해법)

  • Go, G.S.;Ahn, H.T.
    • Journal of computational fluids engineering
    • /
    • v.19 no.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.

A Numerical Analysis on the solution of Poisson Equation by Direct Method (직접법을 이용한 Poisson 방정식 수치해법에 관하여)

  • Y.S. Shin;K.P. Rhee
    • Journal of the Society of Naval Architects of Korea
    • /
    • v.32 no.3
    • /
    • pp.62-71
    • /
    • 1995
  • In the numerical analysis of incompressible unsteady Navier-stokes equation, large time is required for solving the pressure Poisson equation of the elliptic type at each time step. In this paper, a numerical analysis by the direct method is carried out to solve the pressure Poisson equation and the computing time is analyzed as mesh size increases. The pressure Poisson equation can be transformed to the boundary value problem by the Green theorem. The computing time for the convolution type of the domain integral can be reduced by using F.F.T. and the computing time in the direct method depends entirely on obtaining the solution of the boundary value problem. The numerical analysis on the known solutions is carried out and compared for the verification of the direct method. And the numerical analysis on the body boundary and domain decomposition problem are carried out with the computing time less than O($n^{3}$) in the (n.n) mesh.

  • PDF

Discretization of Pressure-Poisson Equation for Solving Incompressible Navier-Stokes Equations Using Non-Staggered Grid (정규격자를 사용한 비압축성 Navier-Stokes 방정식의 수치해석을 위한 압력 Poisson 방정식의 이산화)

  • Kim Y. G.;Kim H. T.;Kim J. J.
    • 한국전산유체공학회:학술대회논문집
    • /
    • 1998.11a
    • /
    • 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.

  • PDF

Application of Spectral Method to Two-Dimensional Unsteady Viscous Flow Analysis (스펙트럴법을 적용한 2차원 비정상 점성유동해석)

  • Shin, Y.S.
    • Journal of the Society of Naval Architects of Korea
    • /
    • v.33 no.4
    • /
    • pp.48-59
    • /
    • 1996
  • The numerical step in the unsteady viscous flow analysis can be divided in the space analysis step satisfying continuity equation and the time marching step. In this study the spectral method is applied to solve the pressure Poisson equation in the space analysis step. If the highest order differential term of the pressure Poisson equation is transformed by Fourier series, pressure arid its first derivatives can be expressed by the integral form of Fourier series. So Gibb's phenomena can be eliminated and the spectral method can be applied to non-periodic problems. The numerical analysis of unsteady viscous flow around 2-dimensional circular cylinder and wing is carried out and compared for verification.

  • PDF

Comparison of different iterative schemes for ISPH based on Rankine source solution

  • Zheng, Xing;Ma, Qing-wei;Duan, Wen-yang
    • International Journal of Naval Architecture and Ocean Engineering
    • /
    • v.9 no.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.

Simulation of free falling rigid body into water by a stabilized incompressible SPH method

  • Aly, Abdelraheem M.;Asai, Mitsuteru;Sonoda, Yoshimi
    • Ocean Systems Engineering
    • /
    • v.1 no.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.

Finite Element Analysis of Incompressible Transient Navier-Stokes Equation using Fractional-Step Methods (Fractional-Step법을 이용한 비압축성 비정상 Navier-Stokes 방정식의 유한 요소해석)

  • Kim, Hyung-Min;Lee, Shin-Pyo
    • Transactions of the Korean Society of Mechanical Engineers B
    • /
    • v.27 no.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.

Computation of Turbulent Flow around Wigley Hull Using 4-Stage Runge-Kutta Scheme on Nonstaggered Grid (정규격자계와 4단계 Range-Kutta법을 사용한 Wigley선형 주위의 난류유동계산)

  • Suak-Hp Van;Hyoung-Tae Kim
    • Journal of the Society of Naval Architects of Korea
    • /
    • v.31 no.3
    • /
    • pp.87-99
    • /
    • 1994
  • Reynolds Averaged Navier-Stokes equations are solved numerically for the computation of turbulent flow around a Wigley double model. A second order finite difference method is applied for the spatial discretization on the nonstaggered grid system and 4-stage Runge-Kutta scheme for the numerical integration in time. In order to increase the time step, residual averaging scheme of Jameson is adopted. Pressure field is obtained by solving the pressure-Poisson equation with the appropriate Neumann boundary condition. For the turbulence closure, 0-equation turbulence model of Baldwin-Lomax is used. Numerical computation is carried out for the Reynolds number of 4.5 million. Comparisons of the computed results with the available experimental data show good agreements for the velocity and pressure distributions.

  • PDF

Electric double layers interactions under condition of variable dielectric permittivity

  • Payam, Amir Farrokh;Fathipour, Morteza
    • Interaction and multiscale mechanics
    • /
    • v.3 no.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.

Simulation of a Pulsating Air Pocket in a Sloshing Tank Using Unified Conservation Laws and HCIB Method (통합보존식 해석과 HCIB 법을 이용한 슬로싱 탱크 내부 갇힌 공기에 의한 압력 진동 모사)

  • Shin, Sangmook
    • Journal of the Society of Naval Architects of Korea
    • /
    • v.58 no.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.