• International Journal of Naval Architecture and Ocean Engineering
• /
• 제8권4호
• /
• pp.320-329
• /
• 2016

An implicit method of solving the six degree-of-freedom rigid body motion equations based on the second order Adams-Bashforth-Moulten method was utilised as an improvement over the leapfrog scheme by making modifications to the rigid body motion solver libraries directly. The implementation will depend on predictor-corrector steps still residing within the hybrid Pressure Implicit with Splitting of Operators - Semi-Implicit Method for Pressure Linked Equations (PIMPLE) outer corrector loops to ensure strong coupling between fluid and motion. Aitken's under-relaxation is also introduced in this study to optimise the convergence rate and stability of the coupled solver. The resulting coupled solver ran on a free floating object tutorial test case when converged matches the original solver. It further allows a varying 70%-80% reduction in simulation times compared using a fixed under-relaxation to achieve the required stability.

• 한국정보과학회논문지:시스템및이론
• /
• 제32권9호
• /
• pp.445-456
• /
• 2005

본 논문은 물리적인 힘을 기반으로 유체의 흐름을 실시간으로 시뮬레이션하기 위하여 유체 의 흐름을 지배하는 Wavier-Stokes 방정식에 대한 빠르고 정확한 풀이 기법을 제안한다 본 논문에서는 Navier-Stokes 방정식에 있는 비선형 항의 속도에 대한 초기값을 Stokes 방정식의 해로써 추정한다. 주어진 비선형 미분방정식의 해에 근사하게 초기값을 추정함으로써 정확하고 안정적인 풀이 기법을 만들 수 있었다. 또한 유한차분법(finite difference method)의 암시적(implicit) 방법 중에서 방대한 계산량을 피할 수 있는 ADI(Alternating Direction Implicit) 방법을 사용함으로써 큰 시간 간격(time-step)에 대해서 시스템이 안정적이며 계산속도 또한 빠르다. 실험 결과들은 특히 연기, 구름과 같이 큰 레이놀드 수(Reynolds number)를 가지는 유체에 대해서 탁월한 성능을 보여주었다.

• 한국전산유체공학회지
• /
• 제5권2호
• /
• pp.20-27
• /
• 2000

An unstructured implicit Euler solver is parallelized on a Cray T3E. Spatial discretization is accomplished by a cell-centered finite volume formulation using an upwind flux differencing. Time is advanced by the Gauss-Seidel implicit scheme. Domain decomposition is accomplished by using the k-way n-partitioning method developed by Karypis. In order to analyze the parallel performance of the solver, flows over a 2-D NACA 0012 airfoil and 3-D F-5 wing were investigated.

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

An unstructured implicit Euler solver is parallelized on a Cray T3E. Spatial discretization is accomplished by a cell-centered finite volume formulation using an unpwind flux differencing. Time is advanced by the Gauss-Seidel implicit scheme. Domain decomposition is accomplished by using the k-way N-partitioning method developed by Karypis. In order to analyze the parallel performance of the solver, flows over a 2-D NACA 0012 airfoil and a 3-D F-5 wing were investigated.

• 한국전산유체공학회지
• /
• 제3권2호
• /
• pp.17-26
• /
• 1998

The three-dimensional incompressible Navier-Stokes equations have been solved by a node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method with Jacobi matrix solver is used for the time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetragedra, prisms, pyramids, hexahedra, or mixed-element grid. Inviscid bump flow is solved to check the accuracy of high order convective flux discretisation. And viscous flows around a circular cylinder and a sphere are studied to show the efficiency and accuracy of the solver.

• 한국전산유체공학회지
• /
• 제14권4호
• /
• pp.13-22
• /
• 2009

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second-order scheme.

• 한국전산유체공학회지
• /
• 제21권1호
• /
• pp.10-18
• /
• 2016

Numerical boundary conditions are as important as the governing equations when analyzing the fluid flows numerically. An explicit boundary condition method updates the solutions at the boundaries with extrapolation from the interior of the computational domain, while the implicit boundary condition method in conjunction with an implicit time integration method solves the solutions of the entire computational domain including the boundaries simultaneously. The implicit boundary condition method, therefore, is more robust than the explicit boundary condition method. In this paper, steady compressible 2-Dimensional Navier-Stokes solver is developed. We present the implicit boundary condition method coupled with LU-SGS(Lower Upper Symmetric Gauss Seidel) method. Also, the explicit boundary condition method is implemented for comparison. The preconditioning Navier-Stokes equations are solved on unstructured meshes. The numerical computations for a number of flows show that the implicit boundary condition method can give accurate solutions.

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

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed. A two-fluid three-field model was adopted for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. The hydrodynamic solver is for the 3D component of a nuclear system code and the component-scale analysis tools for transient two-phase flows. The finite volume method and unstructured grid are adopted, which are useful for the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been adapted to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing problems well.

• 한국전산유체공학회지
• /
• 제12권4호
• /
• pp.44-53
• /
• 2007

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed for a 3D component of a nuclear system code and a component-scale analysis tool. A two-fluid three-field model is used for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. An unstructured grid is adopted for realistic simulations of the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been applied to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the modified numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing very well.

• Ocean Systems Engineering
• /
• 제7권3호
• /
• pp.299-318
• /
• 2017

The paper aims at illustrating several key issues and ongoing efforts for development of a reliable fully-Lagrangian particle-based solver for simulation of hydroelastic slamming. Fluid model is founded on the solution of Navier-Stokes along with continuity equations via an enhanced version of a projection-based particle method, namely, Moving Particle Semi-implicit (MPS) method. The fluid model is carefully coupled with a structure model on the basis of conservation of linear and angular momenta for an elastic solid. The developed coupled FSI (Fluid-Structure Interaction) solver is applied to simulations of high velocity impact of an elastic aluminum wedge and hydroelastic slammings of marine panels. Validations are made both qualitatively and quantitatively in terms of reproduced pressure as well as structure deformation. Several remaining challenges as well as important key issues are highlighted. At last, a recently developed multi-scale MPS method is incorporated in the developed FSI solver towards enhancement of its adaptivity.