• Journal of computational fluids engineering
• /
• v.20 no.1
• /
• pp.32-38
• /
• 2015

A 3-dimensional compressible turbulent boundary layer solver has been developed. A time marching method is used to integrate the turbulent boundary layer equations. While the direct integration of the boundary layer equations is performed for unseparated flow regions, the inverse integration is performed for separated flow regions. The program is verified for flows that have analytical solutions or other numerical results. The solver will be merged with an Euler solver for viscous-inviscid interaction.

• Proceedings of the KIEE Conference
• /
• 1999.11b
• /
• pp.101-103
• /
• 1999

This paper presents an analysis of 2-dimensional(2-D) multi-regions problem using boundary integral equation method(BIEM). When compared with finite element method(FEM), there are only a few unknown variables in BIEM because it implements numerical analysis only for the surface or boundary of a model. As a result, a lot of computational memory and time can be saved. Procedure to analyze 2-D multi-regions problem using potentials and its derivatives in a boundary as unknown variables, first, numerical analysis is performed for each of subregions. And then interface continuity condition is applied to the interface between them and Gauss Quadrature Formula are adopted to solve singular integral in a boundary in this paper.

• Korean Journal of Computational Design and Engineering
• /
• v.22 no.1
• /
• pp.28-36
• /
• 2017

Sharing of CAD data plays a key role in the PLM and a lightweight model is widely used for visualizing and sharing a large data. The lightweight model is mainly composed of triangular elements to minimize file size. There is no problem at all to visually confirm the shape based on these triangular elements but there is a limit to numerically calculate the exact position on the curve or surface. In this paper, a boundary curve generation method using triangular elements is proposed to increase the utilization of lightweight models. After matching connectivity of triangular elements, boundary element edges are extracted. Boundary curves are generated by connecting of these boundary element edges. This proposed method has been tested on several models to demonstrate the feasibility.

• Structural Engineering and Mechanics
• /
• v.39 no.5
• /
• pp.669-682
• /
• 2011

A Chebyshev spectral method (CSM) for the dynamic analysis of non-uniform Timoshenko beams under various boundary conditions and concentrated masses at their ends is proposed. The matrix-based Chebyshev spectral approach was used to construct the spectral differentiation matrix of the governing differential operator and its boundary conditions. A matrix condensation approach is crucially presented to impose boundary conditions involving the homogeneous Cauchy conditions and boundary conditions containing eigenvalues. By taking advantage of the standard powerful algorithms for solving matrix eigenvalue and generalized eigenvalue problems that are embodied in the MATLAB commands, chebfun and eigs, the modal parameters of non-uniform Timoshenko beams under various boundary conditions can be obtained from the eigensolutions of the corresponding linear differential operators. Some numerical examples are presented to compare the results herein with those obtained elsewhere, and to illustrate the accuracy and effectiveness of this method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.7
• /
• pp.1120-1131
• /
• 2003

A recently developed numerical method based on a mixed volume and boundary integral equation method is applied to calculate the accurate stress intensity factors at the crack tips in unbounded isotropic solids in the presence of multiple anisotropic inclusions and cracks subject to external loads. Firstly, it should be noted that this newly developed numerical method does not require the Green's function for anisotropic inclusions to solve this class of problems since only Green's function for the unbounded isotropic matrix is involved in their formulation for the analysis. Secondly, this method takes full advantage of the capabilities developed in FEM and BIEM. In this paper, a detailed analysis of the stress intensity factors are carried out for an unbounded isotropic matrix containing an orthotropic cylindrical inclusion and a crack. The accuracy and effectiveness of the new method are examined through comparison with results obtained from analytical method and volume integral equation method. It is demonstrated that this new method is very accurate and effective for solving plane elastostatic problems in unbounded solids containing anisotropic inclusions and cracks.

• Journal of Mechanical Science and Technology
• /
• v.15 no.5
• /
• pp.671-680
• /
• 2001

Boundary condition is one of the major factors to influence the numerical stability and solution accuracy in numerical analysis. One of the most important physical boundary conditions in the flowfield analysis is the wall boundary condition imposed on the body surface. To solve a two-dimensional Euler equation, totally four numerical wall boundary conditions should be prescribed. Two of them are supplied by the flow tangency condition. The other two conditions, therefore, should be prepared additionally in a suitable way. In this paper, four different sets of wall boundary conditions are proposed and then applied to solve high-speed flowfields around a quarter circle geometry. A two-dimensional compressible Euler solver is prepared based on the finite volume method. This solver hires three different upwind schemes; Steger-Warmings flux vector splitting, Roes flux difference splitting, and Lious advection upstream splitting method. It is found that the way to specify the additional numerical wall boundary conditions strongly affects the overall stability and accuracy of the upwind schemes in high-speed flow calculation. The optimal wall boundary conditions should be also chosen very carefully depending on the numerical schemes used to solve the problem.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.7 s.262
• /
• pp.620-627
• /
• 2007

In this paper, we adapt a modified internal energy non-equilibrium first-order extrapolation thermal boundary condition to the thermal lattice Boltzmann model (TLBM). This model is the double populations approach to simulate hydrodynamic and thermal fields. The bounce-back boundary condition which is a traditional boundary condition of lattice Boltzmann method has only a first order in numerical accuracy at the boundary and numerical instability. A non-equilibrium first-order extrapolation boundary condition has been verified to be of better numerical stability than the bounce-back boundary condition and this boundary condition is proved to be of second-order accuracy for the flat boundaries. The two-dimensional natural convection flow in a square cavity with Pr=0.71 and various Rayleigh numbers are simulated. The results are found to be in good agreement with those of previous studies.

• Nuclear Engineering and Technology
• /
• v.47 no.5
• /
• pp.546-558
• /
• 2015

Chord length sampling method in Monte Carlo simulations is a method used to model spherical particles with random sampling technique in a stochastic media. It has received attention due to the high calculation efficiency as well as user convenience; however, a technical issue regarding boundary effect has been noted. In this study, after analyzing the distribution characteristics of spherical particles using an explicit method, an alternative chord length sampling method is proposed. In addition, for modeling in finite media, a correction method of the boundary effect is proposed. Using the proposed method, sample probability distributions and relative errors were estimated and compared with those calculated by the explicit method. The results show that the reconstruction ability and modeling accuracy of the particle probability distribution with the proposed method were considerably high. Also, from the local packing fraction results, the proposed method can successfully solve the boundary effect problem. It is expected that the proposed method can contribute to the increasing of the modeling accuracy in stochastic media.

• Advances in Computational Design
• /
• v.9 no.3
• /
• pp.213-227
• /
• 2024

In this study, a framework for coupling of the convolution quadrature time-domain boundary element method (CQBEM) and image-based finite element method (IMFEM) is presented for 2-D elastic wave propagation. This coupling method has three advantages: 1) the finite element modeling for heterogeneous areas can be performed without difficulties by using digital data for the analysis model, 2) wave propagation in an infinite domain can be calculated with high accuracy by using the CQBEM, and 3) a small time-step size can be used. In general, a small time-step size cannot be used in the classical time-domain boundary element method. However, the CQBEM used in this analysis can address a small time-step size. This makes it possible to couple the CQBEM and image-based FEM which require a small-time step size. In this study, the formulation and validation of the pro-posed method are described and confirmed by solving fundamental elastic wave scattering problems. As a numerical example, elastic wave scattering in inhomogeneous media is demonstrated using the proposed coupling method.

• Korean Journal of Remote Sensing
• /
• v.25 no.5
• /
• pp.455-464
• /
• 2009

Lee(2009) proposed the boundary-adaptive despeckling method using a Bayesian model which is based on the lognormal distribution for image intensity and a Markov random field(MRF) for image texture. This method employs the Point-Jacobian iteration to obtain a maximum a posteriori(MAP) estimate of despeckled imagery. The boundary-adaptive algorithm is designed to use less information from more distant neighbors as the pixel is closer to boundary. It can reduce the possibility to involve the pixel values of adjacent region with different characteristics. The boundary-adaptive scheme was comprehensively evaluated using simulation data and the effectiveness of boundary adaption was proved in Lee(2009). This study, as an extension of Lee(2009), has suggested a modified iteration algorithm of MAP estimation to enhance computational efficiency and to combine classification. The experiment of simulation data shows that the boundary-adaption results in yielding clear boundary as well as reducing error in classification. The boundary-adaptive scheme has also been applied to high resolution Terra-SAR data acquired from the west coast of Youngjong-do, and the results imply that it can improve analytical accuracy in SAR application.