• 대한수학회지
• /
• 제50권5호
• /
• pp.1129-1163
• /
• 2013

Based on finite element discretization, two linearization approaches to the defect-correction method for the steady incompressible Navier-Stokes equations are discussed and investigated. By applying $m$ times of Newton and Picard iterations to solve an artificial viscosity stabilized nonlinear Navier-Stokes problem, respectively, and then correcting the solution by solving a linear problem, two linearized defect-correction algorithms are proposed and analyzed. Error estimates with respect to the mesh size $h$, the kinematic viscosity ${\nu}$, the stability factor ${\alpha}$ and the number of nonlinear iterations $m$ for the discrete solution are derived for the linearized one-step defect-correction algorithms. Efficient stopping criteria for the nonlinear iterations are derived. The influence of the linearizations on the accuracy of the approximate solutions are also investigated. Finally, numerical experiments on a problem with known analytical solution, the lid-driven cavity flow, and the flow over a backward-facing step are performed to verify the theoretical results and demonstrate the effectiveness of the proposed defect-correction algorithms.

• 한국항만학회지
• /
• 제13권1호
• /
• pp.167-174
• /
• 1999

Numerical analysis of two-fluid flows for both water and air is carried out. Free-Surface flows with an arbitrary deformation have been simulated around two dimensional submerged hydrofoil. The computation is performed using a finite volume method with unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell-wise local mesh refinement. the integration in space is of second order based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels The linear equation systems are solved by conjugate gradient type solvers and the non-linearity of equations is accounted for through picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations the continuity equation the conservation equation of one species and the equations or two turbulence quantities.

• 대한수학회보
• /
• 제48권2호
• /
• pp.377-386
• /
• 2011

For a smooth algebraic curve C of genus g $\geq$ 4, let $SU_C$(r, d) be the moduli space of semistable bundles of rank r $\geq$ 2 over C with fixed determinant of degree d. When (r,d) = 1, it is known that $SU_C$(r, d) is a smooth Fano variety of Picard number 1, whose rational curves passing through a general point have degree $\geq$ r with respect to the ampl generator of Pic($SU_C$(r, d)). In this paper, we study the locus swept out by the rational curves on $SU_C$(r, d) of degree < r. As a by-product, we present another proof of Torelli theorem on $SU_C$(r, d).

• Journal of Mechanical Science and Technology
• /
• 제14권7호
• /
• pp.789-795
• /
• 2000

Free-surface flows with an arbitrary deformation, induced by a submerged hydrofoil, are simulated numerically, considering two-fluid flows of both water and air. The computation is performed by a finite volume method using unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell wise local mesh refinement. The integration in space is of second order, based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels. The linear equations are solved by conjugate gradient type solvers, and the non-linearity of equations is accounted for through Picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations, the continuity equation, the conservation equation of one species, and the equations for two turbulence quantities. Finally, a comparison is quantitatively made at the same speed between the computation and experiment in which the grid sensitivity is numerically checked.

• 대한수학회지
• /
• 제58권4호
• /
• pp.835-847
• /
• 2021

We study behavior of numerical solutions for a nonlinear eigenvalue problem on ℝⁿ that is reduced from a dispersion managed nonlinear Schrödinger equation. The solution operator of the free Schrödinger equation in the eigenvalue problem is implemented via the finite difference scheme, and the primary nonlinear eigenvalue problem is numerically solved via Picard iteration. Through numerical simulations, the results known only theoretically, for example the number of eigenpairs for one dimensional problem, are verified. Furthermore several new characteristics of the eigenpairs, including the existence of eigenpairs inherent in zero average dispersion two dimensional problem, are observed and analyzed.

• 대한지하수환경학회지
• /
• 제6권1호
• /
• pp.23-32
• /
• 1999

Richards 방정식(RE) 해를 구하는 새로운 수치해석적인 방법으로 해적응격자(SAG)법을 개발하였다. SAG 법은 격자생성법을 이용하여 해의 구배가 큰 영역에 더 많은 수의 격자가 밀집되도록 일정한 수의 격자를 자동으로 재분배한다. 이 방법은 좌표변환기법을 이용하여 지배방정식인 RE를 새로운 좌표에서의 RE로 변환하고 유한차분법을 적용하여 방정식의 해를 구한다. 이때 격자점들의 이동은 변환된 RE에 수식으로 반영되기 때문에 고정된 격자점 을 갖는 변환된 영역에서는 해를 구하는 과정에서 내삽법이 불필요하게 된다. 따라서 SAG법은 불포화대에서의 지하수 침투과정을 모사할 때 습윤전선의 이동과 관련하여 발생하는 수치해석적 난제를 크게 개선할 수 있는 방법이다. SAG법과 고정격자를 이용하는 기존의 수정 Picard법을 비교하기 위하여 1차원 침투문제에 대한 수치실험을 실시하였다. 41개의 격자점을 이용한 SAG법은 201점의 고정격자법과 비교하였을 때, 해의 정확도에서는 비슷한 값을 보였으며 계산시간은 반으로 줄어들었다. SAG해의 질량평형과 수렴도는 시간간격 ($\Delta$t)과 해적응 격자생성에 사용된 가중모수 (${\gamma}$)에 크게 영향을 받는 것으로 나타났다. 따라서 SAG법을 이용하여 질량보존적이며 동시에 수렴하는 해를 구하기 위해서는 $\Delta$t와 ${\gamma}$를 자동으로 재조정하는 과정이 요구되며, 이러한 과정은 계산시간을 증가시키는 요인으로 작용할수도 있을 것이다. 본 연구에서 제시된 방법은 특별히 시간에 따른 전선의 이동을 다루는 불포화 유동 및 오염물 거동 문제에서 유용하게 사용될 수 있을 것으로 기대된다.