• Journal of applied mathematics & informatics
• /
• 제9권1호
• /
• pp.1-13
• /
• 2002

In this paper we present numerical methods fur computations of nonequilibrium hypersonic flow of air around bodies including chemical reaction effects and present numerical result of the flow over concave corners. We developed implicit finite difference method to overcome numerical difficulties with the lack of resolution behind the shock and near the body. Using our method we were able to find details of the flow properties near the shock and body and were able to continue the computation of the flow for a long distance from the corner of the body.

• 한국기계연구소 소보
• /
• 통권17호
• /
• pp.167-176
• /
• 1987

It is well known that, about 70-90% of the total drag of a ship is due to the viscous effects. Large amount of theoretical and/or experimental studies on the topic have already been performed. More studies are needed, however, to develop efficient numerical methods which are useful for practical ship design. The present study deals with the theories and numerical methods essential for understanding of real fluid characteristics. Actual ships are not considered because enough computer capacities were not available at the time. Numerical methods, however, are developed to describe complicated ship geometries, transition processes and turbulent boundary layers. The present study can serve as a good start for estimations of viscous ship resistance if an high speed computer is available.

• Earthquakes and Structures
• /
• 제7권2호
• /
• pp.159-178
• /
• 2014

Although the family methods with unconditional stability and numerical dissipation have been developed for structural dynamics they all are implicit methods and thus an iterative procedure is generally involved for each time step. In this work, a new family method is proposed. It involves no nonlinear iterations in addition to unconditional stability and favorable numerical dissipation, which can be continuously controlled. In particular, it can have a zero damping ratio. The most important improvement of this family method is that it involves no nonlinear iterations for each time step and thus it can save many computationally efforts when compared to the currently available dissipative implicit integration methods.

• 대한수학회보
• /
• 제52권5호
• /
• pp.1495-1515
• /
• 2015

This paper deals with the numerical methods for the reconstruction of the source term in a linear parabolic equation from final overdetermination. We assume that the source term has the form f(x)h(t) and h(t) is given, which guarantees the uniqueness of the inverse problem of determining the source term f(x) from final overdetermination. We present the regularization methods for reconstruction of the source term in the whole real line and with Neumann boundary conditions. Moreover, we show the connection of the solutions between the problem with Neumann boundary conditions and the problem with no boundary conditions (on the whole real line) by using the extension method. Numerical experiments are done for the inverse problem with the boundary conditions.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.410-415
• /
• 2001

An analytical and numerical examination of second-order fractional-step methods and boundary condition for the incompressible Navier-Stokes equations is presented. In this study, the compatibility condition for pressure Poisson equation and its boundary conditions, stability, and numerical accuracy of canonical fractional-step methods has been investigated. It has been found that satisfaction of compatibility condition depends on tentative velocity and pressure boundary condition, and that the compatible boundary conditions for type D method and approximately compatible boundary conditions for type P method are proper for divergence-free velocity for type D and approximately divergence-free for type P method. Instability of canonical fractional-step methods is induced by approximation of implicit viscous term with explicit terms, and the stability criteria have been founded with simple model problems and numerical experiments of cavity flow and Taylor vortex flow. The numerical accuracy of canonical fractional-step methods with its consistent boundary conditions shows second-order accuracy except $D_{MM}$ condition, which make approximately first-order accuracy due to weak coupling of boundary conditions.

• 천문학논총
• /
• 제25권3호
• /
• pp.71-76
• /
• 2010

The current status of numerical simulations of galaxy formation is reviewed. After a description of the main numerical simulation techniques, I will present several applications in order to illustrate how numerical simulations have improved our understanding of the galaxy formation process.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제24권3호
• /
• pp.243-291
• /
• 2020

In this paper, multi-block generalized backward differentiation methods for numerical solutions of ordinary differential and differential algebraic equations are introduced. This class of linear multi-block methods is implemented as multi-block boundary value methods (MB₂ VMs). The root distribution of the stability polynomial of the new class of methods are determined using the Wiener-Hopf factorization of a matrix polynomial for the purpose of their correct implementation. Numerical tests, showing the potential of such methods for output of multi-block of solutions of the ordinary differential equations in the new approach are also reported herein. The methods which output multi-block of solutions of the ordinary differential equations on application, are unlike the conventional linear multistep methods which output a solution at a point or the conventional boundary value methods and multi-block methods which output only a block of solutions per step. The MB₂ VMs introduced herein is a novel approach at developing very large scale integration methods (VLSIM) in the numerical solution of differential equations.

• Nuclear Engineering and Technology
• /
• 제16권1호
• /
• pp.21-28
• /
• 1984

수치확산을 포함한 truncation오차의 줄임은 수치해석의 중요한 과제가 되어왔다. Stream line방법이 교차수치 확산과 비확산형의 truncation 오차를 제거하기 위하여 고안되었다. 또한, stream line방법과 유한 차분법이 합쳐진 2단계 stream line방법이 비압축성 난류유동의 지배 방정식을 풀기 위하여 고안되었다. 이 방법은 유한 차분법과 비교되었으며, 두 방법 모두 실험자료와 비교되었다. 그리고, 두 방법의 truncation 오차를 비교하기 위하여 truncation 오차 분석이 행해졌다

• Korean Journal of Mathematics
• /
• 제13권2호
• /
• pp.127-137
• /
• 2005

The objective of numerical analysis is to devise and analyze efficient algorithms or numerical methods for equations arising in mathematical modeling for science and engineering. In this article, we present some recent topics in computational mathematics, specially in the finite element method and overview the development of the mixed finite element method in the context of second order elliptic and parabolic problems. Multiscale methods such as MsFEM, HMM, and VMsM are included.

• Structural Engineering and Mechanics
• /
• 제65권1호
• /
• pp.121-128
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.