• 대한조선학회논문집
• /
• 제36권1호
• /
• pp.70-81
• /
• 1999

선체주위의 점성유동을 계산하기 위해서는 수치계산을 위한 3차원 공간 격자계가 필요하다. 본 논문에서는 타원형 미분 방정식인 Poisson 방정식의 해를 이용하여 3차원 공간 격자계를 구성하는 방법을 개발하였다. 2차원에서 사용되던 Sorenson방법을 3차원으로 확장하여 격자계의 분포 및 교차 각도를 지정하는 함수를 정의하게 하였다. 미분방정식의 해는 weighting function scheme과 modified strongly implicit procedure를 사용하여 구하였고, 3차원 내부 격자계와 경계면과의 매끄러운 연결을 위해 trans-finite interpolation을 이용하였다. 적용예로서 컨테이너 운반선과 대형 유조선 주위의 난류유동 계산을 위한 공간 격자계를 보였다.

• Nuclear Engineering and Technology
• /
• 제45권6호
• /
• pp.707-720
• /
• 2013

We analyze piecewise homogenization with flux-weighted cross sections and preservation of averaged currents at the boundary of the homogenized domain. Introduction of a set of flux discontinuity ratios (FDR) that preserve reference interface currents leads to preservation of averaged region reaction rates and fluxes. We consider the class of numerical discretizations with one degree of freedom per volume and per surface and prove that when the homogenization and computing meshes are equal there is a unique solution for the FDRs which exactly preserve interface currents. For diffusion submeshing we introduce a Jacobian-Free Newton-Krylov method and for all cases considered obtain an 'exact' numerical solution (eight digits for the interface currents). The homogenization is completed by extending the familiar full assembly homogenization via flux discontinuity factors to the sides of regions laying on the boundary of the piecewise homogenized domain. Finally, for the familiar nodal discretization we numerically find that the FDRs obtained with no submesh (nearly at no cost) can be effectively used for whole-core diffusion calculations with submesh. This is not the case, however, for cell-centered finite differences.

• 대한기계학회논문집B
• /
• 제28권5호
• /
• pp.580-586
• /
• 2004

The ray effects of finite volume method (FVM) or discrete ordinate method (DOM) are known to show a non-physical oscillation in solution of radiative heat transfer on a boundary. This wiggling behavior is caused by the finite discretization of the continuous control angle. This article proposes a combined procedure of the Monte-Carlo and finite-volume method (CMCFVM) for solving radiative heat transfer in absorbing, emitting, and an-isotropically scattering medium with an isolated boundary heat source. To tackle the problem, which is especially pronounced in a medium with an isolated heat source, the CMCFVM is suggested here and successfully applied to a two-dimensional circular geometry.

• International Journal of Ocean System Engineering
• /
• 제1권3호
• /
• pp.136-143
• /
• 2011

In this paper, three-dimensional free-surface flows are simulated by using two different numerical methods, the constrained interpolation profile (CIP)-based and finite volume (FV)-based methods. In the CIP-based method, the governing equations are solved on stationary staggered Cartesian grids by a finite difference method, and an immersed boundary technique is applied to deal with wave-body interactions. In the FV-based method, the governing equations are solved by applying collocated finite volume discretization, and body-fitted meshes are used. A free-surface boundary is considered as the interface of the multi-phase flow with air and water, and a volumeof-fluid (VOF) approach is applied to trace the free surface. Among many variations of the VOF-type method, the tangent of hyperbola for interface capturing (THINC) and the compressive interface capturing scheme for arbitrary meshes (CICSAM) techniques are used in the CIP-based method and FV-based method, respectively. Numerical simulations have been carried out for dam-breaking and wave-body interaction problems. The computational results of the two methods are compared with experimental data and their differences are observed.

• Ocean Systems Engineering
• /
• 제5권4호
• /
• pp.301-318
• /
• 2015

The present investigations introduce the shell-finite element discretization for the dynamics of slender marine pipelines. A long catenary pipeline, corresponding to a particular Steel Catenary Riser (SCR), is investigated under long-standing cyclic loading. The long structure is divided into smaller tubular parts which are discretized with 8-node planar shell elements. The transient analysis of each part is carried out by the implicit time integration scheme, within a Finite Elements (FE) solver. The time varying external loads and boundary conditions on each part are the results of a prior solution of an integrated line-dynamics model. The celebrated FE approximation can produce a more detailed stress distribution along the structural surface than the simplistic "line-dynamics" approach.

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

An efficient numerical method to solve the unsteady incompressible Navier-Stokes equations is developed. A fully implicit time advancement is employed to avoid the CFL(Courant-Friedrichs-Lewy) restriction, where the Crank-Nicholson discretization is used for both the diffusion and convection terms. Based on a block LU decomposition, velocity-pressure decoupling is achieved in conjunction with the approximate factorization. Main emphasis is placed on the additional decoupling of the intermediate velocity components with only n th time step velocity The temporal second-order accuracy is Preserved with the approximate factorization without any modification of boundary conditions. Since the decoupled momentum equations are solved without iteration, the computational time is reduced significantly. The present decoupling method is validated by solving the turbulent minimal channel flow unit.

• 대한수학회논문집
• /
• 제20권3호
• /
• pp.563-578
• /
• 2005

In [Z. Cai and B. C. Shin, SIAM J. Numer. Anal. 40 (2002), 307-318], we developed the discrete first-order system least squares method for the second-order elliptic boundary value problem by directly approximating $H(div){\cap}H(curl)-type$ space based on the Helmholtz decomposition. Under general assumptions, error estimates were established in the $L^2\;and\;H^1$ norms for the vector and scalar variables, respectively. Such error estimates are optimal with respect to the required regularity of the solution. In this paper, we study solution methods for solving the system of linear equations arising from the discretization of variational formulation which possesses discrete biharmonic term and focus on numerical results including the performances of multigrid preconditioners and the finite element accuracy.

• 한국소음진동공학회논문집
• /
• 제22권8호
• /
• pp.791-799
• /
• 2012

In this paper free vibration analysis of symmetric and cross-ply elastic laminated shells based on FSDT was performed through discretization of equations of motion and boundary condition. Structural model of laminated composite cylindrical shells subjected to a combination of magnetic and thermal fields is developed via Hamilton's variational principle. These coupled equations of motion are based on the electromagnetic equations(Faraday, Ampere, Ohm, and Lorenz equations) and thermal equations which are involved in constitutive equations. Variations of dynamic characteristics of composite shells with applied magnetic field, temperature gradient, and stacking sequence are investigated and pertinent conclusions are derived.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 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.

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

Three-Dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. 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. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for the various contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreements.