• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2401-2405
• /
• 2008

In the present study, a three-dimensional least square/level set based two-phase flow code was developed for the simulation of three-dimensional sloshing problems using finite element discretization. The present method can be utilized for the analysis of a free surface flow problem in a complex geometry due to the feature of FEM. Since the finite element method is employed for the spatial discretization of governing equations, an unstructured mesh can be naturally adopted for the level set simulation of a free surface flow without an additional load for the code development except that solution methods of the hyperbolic type redistancing and advection equations of the level set function should be devised in order to give a bounded solution on the unstructured mesh. From the numerical experiments of the present study, it is shown that the proposed method is both robust and accurate for the simulation of three-dimensional sloshing problems.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2558-2561
• /
• 2008

In the present study, a parallel Taylor-Galerkin/level set based two-phase flow code was developed using finite element discretization and domain decomposition method based on MPI (Message Passing Interface). The proposed method can be utilized for the analysis of a large scale free surface problem in a complex geometry due to the feature of FEM and domain decomposition method. Four-step fractional step method was used for the solution of the incompressible Navier-Stokes equations and Taylor-Galerkin method was adopted for the discretization of hyperbolic type redistancing and advection equations. A Parallel ILU(0) type preconditioner was chosen to accelerate the convergence of a conjugate gradient type iterative solvers. From the present parallel numerical experiments, it has been shown that the proposed method is applicable to the simulation of large scale free surface flows.

• Journal of applied mathematics & informatics
• /
• v.39 no.5_6
• /
• pp.623-641
• /
• 2021

In this paper, numerical treatment of singularly perturbed parabolic delay differential equations is considered. The considered problem have small delay on the spatial variable of the reaction term. To treat the delay term, Taylor series approximation is applied. The resulting singularly perturbed parabolic PDEs is solved using Crank Nicolson method in temporal direction with non-standard finite difference method in spatial direction. A detail stability and convergence analysis of the scheme is given. We proved the uniform convergence of the scheme with order of convergence O(N^-1 + (∆t)²), where N is the number of mesh points in spatial discretization and ∆t is mesh length in temporal discretization. Two test examples are used to validate the theoretical results of the scheme.

• Journal of applied mathematics & informatics
• /
• v.39 no.5_6
• /
• pp.655-676
• /
• 2021

A parameter uniform numerical scheme is proposed for solving singularly perturbed parabolic partial differential-difference convection-diffusion equations with a small delay and advance parameters in reaction terms and spatial variable. Taylor's series expansion is applied to approximate problems with the delay and advance terms. The resulting singularly perturbed parabolic convection-diffusion equation is solved by utilizing the implicit Euler method for the temporal discretization and finite difference method for the spatial discretization on a uniform mesh. The proposed numerical scheme is shown to be an ε-uniformly convergent accurate of the first order in time and second-order in space directions. The efficiency of the scheme is proved by some numerical experiments and by comparing the results with other results. It has been found that the proposed numerical scheme gives a more accurate approximate solution than some available numerical methods in the literature.

• Nuclear Engineering and Technology
• /
• v.53 no.12
• /
• pp.3861-3878
• /
• 2021

The steady-state simplified spherical harmonics equations (SP_N equations) are a higher order approximation to the neutron transport equations than the neutron diffusion equation that also have reasonable computational demands. This work extends these results for the analysis of transients by comparing of two formulations of time-dependent SP_N equations considering different treatments for the time derivatives of the field moments. The first is the full system of equations and the second is a diffusive approximation of these equations that neglects the time derivatives of the odd moments. The spatial discretization of these methodologies is made by using a high order finite element method. For the time discretization, a semi-implicit Euler method is used. Numerical results show that the diffusive formulation for the time-dependent simplified spherical harmonics equations does not present a relevant loss of accuracy while being more computationally efficient than the full system.

• Nuclear Engineering and Technology
• /
• v.16 no.2
• /
• pp.47-57
• /
• 1984

Mathematically-rigorous time-volume averaged conservation equations were simplified to established the differential equations of THERMIT-6S, which is a two-fluid 3-D code. The difference equations of THERMIT-6S were obtained by discretizing the proceeding set of differential equations. The spatial discretization is characterized by a first-order spatial scheme, donor cell method, and staggered mesh layout. For time discretization, a first order semi-implicit scheme treats implictly sonic terms and terms relating to local transport phenomena and explicitly convective terms. The results were linearized by the Newton-Raphson method. In order to construct the reduced pressure equation, the linearized equations were manipulated so that all variables are coupled between mesh cells through only the pressure variable. By simulating numerically the OPERA-15 experiment, it was found that THERMIT-6S is a very powerful code in predicting reactor behavior after sodium boiling including flow coastdown, reversal flow and flow oscillation.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.13 no.4
• /
• pp.391-396
• /
• 2003

This paper proposes a new approach that converts continuous-valued attributes to categorical-valued ones considering the distribution of target attributes(classes). In this approach, It can be possible to get optimal interval boundaries by considering the distribution of data itself without any requirements of parameters. For each attributes, the distribution of target attributes is projected to one-dimensional space. And this space is clustered according to the criteria like as the density value of each target attributes and the amount of overlapped areas among each density values of target attributes. Clusters which are made in this ways are based on the probabilities that can predict a target attribute of instances. Therefore it has an interval boundaries that minimize a loss of information of original data. An improved performance of proposed discretization method can be validated using C4.5 algorithm and UCI Machine Learning Data Repository data sets.

• The Journal of the Korea Contents Association
• /
• v.10 no.9
• /
• pp.97-106
• /
• 2010

Various time series representation methods have been suggested in order to process time series data efficiently and effectively. SAX is the representative time series representation method combining segmentation and discretization techniques, which has been successfully applied to the time series classification task. But SAX requires a large number of segments in order to represent the meaningful dynamic patterns of time series accurately, since it loss the dynamic property of time series in the course of smoothing the movement of time series. Therefore, this paper suggests a new time series representation method that combines PIPs detection and Persist discretization techniques. The suggested method represents the dynamic movement of high-diemensional time series in a lower dimensional space by detecting PIPs indicating the important inflection points of time series. And it determines the optimal discretizaton ranges by applying self-transition and marginal probabilities distributions to KL divergence measure. It minimizes the information loss in process of the dimensionality reduction. The suggested method enhances the performance of time series classification task by minimizing the information loss in the course of dimensionality reduction.

• KIEAE Journal
• /
• v.16 no.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.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.45 no.3
• /
• pp.233-240
• /
• 2017

The flowfield around transonic wing-body configuration was simulated using in-house CFD code and compared with the experimental data to understand the influence of several features of CFD(Computational Fluid Dynamics) ; grid dependency, turbulence models, spatial discretization, and viscosity. The wing-body configuration consists of a simple planform RAE Wing 'A' with an RAE 101 airfoil section and an axisymmetric body. The in-house CFD code is a compressible Euler/Navier-Stokes solver based on unstructured grid. For the turbulence model, the $k-{\omega}$ model, the Spalart-Allmaras model, and the $k-{\omega}$ SST model were applied. For the spatial discretization method, the central differencing scheme with Jameson's artificial viscosity and Roe's upwind differencing scheme were applied. The results calculated were generally in good agreement with experimental data. However, it was shown that the pressure distribution and shock-wave position were slightly affected by the turbulence models and the spatial discretization methods. It was known that the turbulent viscous effect should be considered in order to predict the accurate shock wave position.