• The Pure and Applied Mathematics
• /
• v.14 no.4
• /
• pp.283-287
• /
• 2007

An extension of the contraction mapping theorem is provided in a Banach space setting to approximate fixed points of operator equations. Our approach is justified by numerical examples where our results apply whereas the classical contraction mapping principle cannot.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.135-140
• /
• 2000

Numerical calculations are carried out in order to evaluate the performance of low-Re Reynolds stress model based on SSG model for a swirling turbulent flow in a pipe. The results are compared with those of $\kappa-\epsilon$ model and GL model, and the experimental data. The finite volume method is used for the discretization, and the power-law scheme is employed as a numerical scheme. The SIMPLE algorithm is used for velocity-Pressure correction in the governing equations.

• Journal of applied mathematics & informatics
• /
• v.6 no.2
• /
• pp.455-466
• /
• 1999

In this note we present a numerical scheme for finding an approxximate solution of an equation which can be viewed as a model for spatial diffusion of age-depednent biological populations. Discretization of the model yields a linear system with a block tridi-agonal matrix. Our main concern will be discussion of stability for this scheme by examining the eigenvalues of the block tridiagonal matrix. Numerical results are presented.

• Korean Journal of Mathematics
• /
• v.7 no.2
• /
• pp.245-269
• /
• 1999

Multigrid method and acceleration by conjugate gradient method for first-order system least squares (FOSLS) using bilinear finite elements are developed for various boundary value problems of planar linear elasticity. They are two-stage algorithms that first solve for the displacement flux variable, then for the displacement itself. This paper focuses on solving for the displacement flux variable only. Numerical results show that the convergence is uniform even as the material becomes nearly incompressible. Computations for convergence factors and discretization errors are included. Heuristic arguments to improve the convergences are discussed as well.

• Journal of applied mathematics & informatics
• /
• v.28 no.5_6
• /
• pp.1331-1345
• /
• 2010

This paper describes a numerical solution of Navier-Stokes equations. It includes algorithms for discretization by finite element methods and a posteriori error estimation of the computed solutions. In order to evaluate the performance of the method, the numerical results are compared with some previously published works or with others coming from commercial code like ADINA system.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.7 no.2
• /
• pp.65-73
• /
• 2003

In this work, we use a numerical method to solve the nonlinear integral equation of the second kind when the kernel of the integral equation in the logarithmic function form or in Carleman function form. The solution has a computing time requirement of $0(N^2)$, where (2N +1) is the number of discretization points used. Also, the error estimate is computed.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.5-9
• /
• 2010

A numerical wind tunnel simulation is performed in order to predict wind loads acting on a building. The aim of the present study is to suggest a guideline for the numerical wind tunnel analysis, which could provide more detail wind load distributions compared to the wind code and expensive wind tunnel experiments. To validate the present numerical simulation, wind-induced loads on a 6 m cube model is predicted. Atmospheric boundary layer is used as a inlet boundary condition. Various effect of numerical methods are investigated such as size of computational domain, grid density, turbulence model and discretization scheme. The appropriate procedure for the numerical wind tunnel analysis is suggested through the present study.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.1049-1056
• /
• 2005

A conservative pressure-based finite-volume numerical method has been developed for computing flow and heat transfer by using an unstructured grid system. The method admits arbitrary convex polyhedra. Care is taken in the discretization and solution procedures to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are found by a novel second-order accurate spatial discretization. Momentum interpolation is used to prevent pressure checkerboarding and the SIMPLE algorithm is used for pressure-velocity coupling. The resulting set of coupled nonlinear algebraic equations is solved by employing a segregated approach, leading to a decoupled set of linear algebraic equations fer each dependent variable, with a sparse diagonally dominant coefficient matrix. These equations are solved by an iterative preconditioned conjugate gradient solver which retains the sparsity of the coefficient matrix, thus achieving a very efficient use of computer resources.

• International Journal of Aerospace System Engineering
• /
• v.2 no.2
• /
• pp.67-72
• /
• 2015

The main concern of this study is to integrate the EARSM into an industrial RANS solver in conjunction with the $k-{\omega}$ model, as proposed by Hellsten (EARSMKO2005). In order to improve the robustness, particular limiters are introduced to turbulent conservative variables, and a suitable full-approximation storage (FAS) multi-grid (MG) strategy is designed to incorporate turbulence model equations. The present limiters and MG strategy improve both robustness and efficiency significantly but without degenerating accuracy. Two discretization approachs for velocity gradient on cell interfaces are implemented and compared with each other. Numerical results of a three-dimensional supersonic square duct flow show that the proper discretization of velocity gradient improves the accuracy essentially. To assess the capability of the resulting EARSM $k-{\omega}$ model to predict complex engineering flow, the case of Common Research Model (CRM, Wing-Body) is performed. All the numerical results demonstrate that the resulting model performs well and is comparable to the standard two-equation models such as SST $k-{\omega}$ model in terms of computational effort, thus it is suitable for industrial applications.

• Geomechanics and Engineering
• /
• v.18 no.3
• /
• pp.235-245
• /
• 2019

This paper presents a novel methodology for face stability assessment of rock tunnels under water table by combining the kinematical approach of limit analysis and numerical simulation. The tunnels considered in this paper are excavated in fractured rock masses characterized by the Hoek-Brown failure criterion. In terms of natural rock deposition, a more convincing case of depth-dependent mi, GSI, D and ${\sigma}_c$ is taken into account by proposing the horizontally layered discretization technique, which enables us to generate the failure surface of tunnel face point by point. The vertical distance between any two adjacent points is fixed, which is beneficial to deal with stability problems involving depth-dependent rock parameters. The pore water pressure is numerically computed by means of 3D steady-state flow analyses. Accordingly, the pore water pressure for each discretized point on the failure surface is obtained by interpolation. The parametric analysis is performed to show the influence of depth-dependent parameters of $m_i$, GSI, D, ${\sigma}_c$ and the variation of water table elevation on tunnel face stability. Finally, several design charts for an undisturbed tunnel are presented for quick calculations of critical support pressures against face failure.