• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2012년도 제45회 KOSCO SYMPOSIUM 초록집
• /
• pp.305-307
• /
• 2012

This study was computationally explored how the fuel autoignition reactivity was affected by operating parameters such as fuel, pressure, intake temperatures, engine speed and EGR compositions for HCCI combustion. This is done for DME and CHEMKIN-PRO was used as a solver. At first, influence of the operating parameters and EGR compositions were showed. And then, in order to clarify the mechanism of them on autoignition reactivity, data-sets of kinetic were analyzed to investigate the elementary reaction path for heat release at transient tempeatures by using contribution matrix.

• 대한수학회보
• /
• 제59권2호
• /
• pp.507-528
• /
• 2022

We propose a path integral method to construct a time stepwise local volatility for the stock index market under Dupire's model. Our method is focused on the pricing with the Monte Carlo Method (MCM). We solve the problem of randomness of MCM by applying numerical integration. We reconstruct this task as a matrix equation. Our method provides the analytic Jacobian and Hessian required by the nonlinear optimization solver, resulting in stable and fast calculations.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2001년도 추계학술대회 논문집
• /
• pp.113-116
• /
• 2001

Finite element analysis programs have been for metal forming process design They will become more and more important in understanding forming process For large-scale forging analysis problems, the performance of a linear equation solver is very important for the overall efficiency of the analysis code. With problem size increased, the computation time needs to be reduced, which is spent on setting the system of algebraic equations associated with finite element model Many matrix solvers have been developed and used usefully in finite element program for this purpose.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.115-119
• /
• 2009

We numerically investigated propagation of various waves in the two-phase flows such as sound wave, shock wave, rarefaction wave, and contact discontinuity in terms of pressure, void fraction, velocity and density of the two phases. The waves have been generated by a hydrodynamic shock tube, a pair of symmetric impulsive expansion, impulsive pressure and impulsive void waves. The six compressible two-fluid two-phase conservation laws with interfacial friction terms have been solved in two fractional steps. The first PDE Operator is solved by the HLL scheme and the second Source Operator by the semi-implicit stiff ODE solver. In the HLL scheme, the fastest wave speeds were estimated by the analytic eigenvalues of an approximate Jacobian matrix. We have discussed how the interfacial friction terms affect the wave structures in the numerical solution.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제13권1호
• /
• pp.13-20
• /
• 2009

The IPSAP which is a finite element analysis program has been developed for high parallel performance computing. This program consists of various analysis modules - stress, vibration and thermal analysis module, etc. The M orthogonal block Lanczos algorithm with shiftinvert transformation is used for solving eigenvalue problems in the vibration module. And the multifrontal algorithm which is one of the most efficient direct linear equation solvers is applied to factorization and triangular system solving phases in this block Lanczos iteration routine. In this study, the performance enhancement procedures of the IPSAP are composed of the following stages: 1) communication volume minimization of the factorization phase by modifying parallel matrix subroutines. 2) idling time minimization in triangular system solving phase by partial inverse of the frontal matrix and the LCM (least common multiple) concept.

• International Journal of Control, Automation, and Systems
• /
• 제1권2호
• /
• pp.184-193
• /
• 2003

Biaffine Matrix Inequality (BMI) is known to provide the most general framework in control synthesis, but problems involving BMI's are very difficult to solve because nonconvex optimization should be solved. In the previous paper, we proposed a new solver for problems involving BMI's using Evolutionary Algorithms (EA). In this paper, we solve several control synthesis examples such as Reduced-order control, Simultaneous stabilization, Multi-objective control, $H_{\infty}$ optimal control, Maxed $H_2$ / $H_{\infty}$control design, and Robust $H_{\infty}$ control. Each of these problems is formulated as the standard BMI form, and solved by the proposed algorithm. The performance in each case is compared with those of conventional methods.

• 정보처리학회논문지A
• /
• 제8A권3호
• /
• pp.227-234
• /
• 2001

In this paper we propose a block-type parallel preconditioner for solving large sparse nonsymmetric linear systems, which we expect to be scalable. It is Multi-Color Block SOR preconditioner, combined with direct sparse matrix solver. For the Laplacian matrix the SOR method is known to have a nondeteriorating rate of convergence when used with Multi-Color ordering. Since most of the time is spent on the diagonal inversion, which is done on each processor, we expect it to be a good scalable preconditioner. We compared it with four other preconditioners, which are ILU(0)-wavefront ordering, ILU(0)-Multi-Color ordering, SPAI(SParse Approximate Inverse), and SSOR preconditiner. Experiments were conducted for the Finite Difference discretizations of two problems with various meshsizes varying up to $1025{\times}1024$. CRAY-T3E with 128 nodes was used. MPI library was used for interprocess communications, The results show that Multi-Color Block SOR is scalabl and gives the best performances.

• KIEAE Journal
• /
• 제16권1호
• /
• pp.81-87
• /
• 2016

Purpose: Heat transfers phenomena are described by the second order partial differential equation and its boundary conditions. In a three-dimensional structure of a building, the heat transfer phenomena generally include more than one material, and thus, become complicate. The analytic solutions are useful to understand heat transfer phenomena, but they can hardly be applied in engineering or design problems. Engineers and designers have generally been forced to use numerical methods providing reliable results. Finite volume methods with the unstructured grid system is only the suitable means of the analysis for the complex and arbitrary domains. Method: To obtain an numerical solution, a discretization method, which approximates the differential equations, and the interpolation methods for temperature and heat flux between two or more materials are required. The discretization methods are applied to small domains in space and time, and these numerical solutions form the descretized equations provide approximated solutions in both space and time. The accuracy of numerical solutions is dependent on the quality of discretizations and size of cells used. The higher accuracy, the higher numerical resources are required. The balance between the accuracy and difficulty of the numerical methods is critical for the success of the numerical analysis. A simple and easy interpolation methods among multiple materials are developed. The linear equations are solved with the BiCGSTAB being a effective matrix solver. Result: This study provides an overview of discretization methods, boundary interface, and matrix solver for the 3-dimensional numerical heat transfer including two materials.

• 한국전산유체공학회지
• /
• 제19권2호
• /
• pp.86-92
• /
• 2014

In this paper, prediction of separation trajectory for Two-stage-To-Orbit space launch vehicle has been numerically simulated by using an aerodynamic database based on steady state analysis. Aerodynamic database were obtained for matrix of longitudinal and vertical positions. The steady flow simulations around the launch vehicle have been made by using a 3-D RANS flow solver based on unstructured meshes. For this purpose, a vertex-centered finite-volume method was adopted to discretize inviscid and viscous fluxes. Roe's finite difference splitting was utilized to discretize the inviscid fluxes, and the viscous fluxes were computed based on central differencing. To validate this flow solver, calculations were made for the wind-tunnel experiment model of the LGBB TSTO vehicle configuration on steady state conditions. Aerodynamic database was constructed by using flow simulations based on test matrix from the wind-tunnel experiment. ANN(Artificial Neural Network) was applied to construct interpolation function among aerodynamic variables. Separation trajectory for TSTO launch vehicle was predicted from 6-DOF equation of motion based on the interpolated function. The result of present separation trajectory calculation was compared with the trajectory using experimental database. The predicted results for the separation trajectory shows fair agreement with reference[4] solution.

• 한국항공우주학회지
• /
• 제33권9호
• /
• pp.34-40
• /
• 2005

본 연구에서는 비정상 유동해석을 위한 CFD 코드의 개발을 위해 대각화 ADI 기법을 적용한 정상 해석기법과 내재적 이중시간 전진기법을 도입하였다. 정상상태 Navier-Stokes 방정식의 Jacobian 행렬은 비점성항에 대해서만 적용하였고 여기에 내재적 인공점성 연산자를 첨가하여 블록 5대각 행렬을 유도하였다. 시간단축을 위해 스칼라 5대각 행렬로 대체하였다. 가상시간에 대한 정상상태기법에 실시간에 대한 미분항이 포함된 새로운 잔류항을 정의하였다. 가상시간에 대해 수렴된 해로부터 실시간 해를 구하고 시간에 대해 적분을 수행하는 내재적 이중시간 전진기법을 이용한 비정상 Navier-Stokes 코드를 개발하였다. 이에 대한 검증으로 정지한 유체속에 진동하는 평판문제, 원기둥 후방의 주기적인 Karman 와류생성, 이중원호 익형주위의 충격파 진동문제등을 수치해석하여 이론치, 실험치, 타연구자의 계산결과와 비교, 분석하였다.