• Structural Engineering and Mechanics
• /
• v.22 no.1
• /
• pp.1-16
• /
• 2006

This paper presents a virtual boundary element-equivalent collocation method (VBEM) for the plane magnetoelectroelastic solids, which is based on the fundamental solutions of the plane magnetoelectroelastic solids and the basic idea of the virtual boundary element method for elasticity. Besides all the advantages of the conventional boundary element method (BEM) over domain discretization methods, this method avoids the computation of singular integral on the boundary by introducing the virtual boundary. In the end, several numerical examples are performed to demonstrate the performance of this method, and the results show that they agree well with the exact solutions. So the method is one of the efficient numerical methods used to analyze megnatoelectroelastic solids.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2797-2802
• /
• 2007

It has been confirmed that implementation of the no-slip boundary conditions for the lattice-Boltzmann method play an important role in the overall accuracy of the numerical solutions as well as the stability of the solution procedure. We in this paper propose a new algorithm, i.e. the method of the dynamic boundary condition for no-slip boundary condition. The distribution functions on the wall along each of the links across the physical boundary are assumed to be composed of equilibrium and nonequilibrium parts which inherit the idea of Guo's extrapolation method. In the proposed algorithm, we apply a dynamic equation to reflect the computational slip velocity error occurred on the actual wall boundary to the correction; the calculated slip velocity error dynamically corrects the fictitious velocity on the wall nodes which are subsequently employed to the computation of equilibrium distribution functions on the wall nodes. Along with the dynamic selfcorrecting process, the calculation efficiently approaches the steady state. Numerical results show that the dynamic boundary method is featured with high accuracy and simplicity.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.36 no.11
• /
• pp.777-782
• /
• 1987

A new composite method of finite element and boundary integral methods is presented to solve the two dimensional magnetostatic field problems with open boundary. The method can deal with the current source of the boundary integral regin where the boundary integral method is applied, and also Neumann and Dirichlet boundary conditions at the interfacial boundary between the boundary integral region and the finite element region where the finite element method is applied. The new approach has been applied to a simple linear problem to verify the usefulness. It is shown that the proposed algorithm gives more accurate results than the finite element methed under the same elementdiscretization.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.10a
• /
• pp.160-165
• /
• 1997

Structural boundary condition is very important as a part of a structural system because it determines the dynamic characteristics of the structure. It is often to experience that experimental measurements of structural dynamic characteristics are somewhat different from the analytical predictions in which idealized boundary conditions are usually assumed. However, real structural boundary conditions are not so ideal; not perfectly clamped, for instance. Thus this paper introduces a new method to determine the non-ideal structural boundary conditions in the frequency domain. In this method, structural boundary conditions are modeled by both extensional (vertical) and torsional elastic springs. The effective springs are then determined from experimental FRFs (frequency response functions) by using the spectral element method (SEM). For a cantilevered beam experiments are conducted to determine the real boundary conditions in terms of effective springs. Dynamic characteristics (analytically predicted) based on identified boundary conditions are found to be much closer to experimental measurements when compared with those based on ideal boundary conditions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.6
• /
• pp.619-627
• /
• 2004

In most industrial applications, the geometrical complexity is combined with the moving boundaries. These problems considerably increase the computational difficulties since they require, respectively, regeneration and deformation of the grid. As a result, engineering flow simulation is restricted. In order to solve this kind of problems the immersed boundary method was developed. In this study, the immersed boundary method is applied to the numerical simulation of stationary, rotating and oscillating cylinders in the 2-dimensional square cavity. No-slip velocity boundary conditions are given by imposing feedback forcing term to the momentum equation. Besides, this technique is used with a second-order accurate interpolation scheme in order to improve the accuracy of flow near the immersed boundaries. The governing equations for the mass and momentum using the immersed boundary method are discretized on the non-staggered grid by using the finite volume method. The results agree well with previous numerical and experimental results. This study presents the possibility of the immersed boundary method to apply to the complex flow experienced in the industrial applications. The usefulness of this method will be confirmed when we solve the complex geometries and moving bodies.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.77-84
• /
• 2014

The lattice Boltzmann (LB) method has been used to simulate rarefied gas flows in a micro-system as an alternative tool. However, previous results were mainly focused on a simple geometry with flat walls because the LB method is modeled on uniform Cartesian lattices. When previous boundary conditions for the microflows are applied to curved walls, the use of them requires approximation of the curved boundary by a series of stair steps, and introduces additional errors. For macroflows, no-slip curved wall boundary treatments have been developed remarkably in order to overcome these limits. However, the investigations for the slip curved wall boundary have rarely been performed for microflows. In this work, a curved boundary treatment of the LB method for a slip flow has been introduced. The results of the LB method for 2D microchannel and 3D microtube flows are in excellent agreement with the analytical solutions.

• International Journal of CAD/CAM
• /
• v.8 no.1
• /
• pp.11-20
• /
• 2009

The conventional methods of boundary-conformed 2D surfaces generation usually yield some problems. This paper deals with two boundary-conformed 2D surface generation methods, one conventional approach, the linear Coons method, and a new method, boundary-conformed interpolation. In this new method, unidirectional 2D surface has been generated using some of the geometric properties of the given boundary curves. A method of simultaneous displacement of the interpolated curves from the opposite boundaries has been adopted. The geometric properties considered for displacements include weighted combination of angle bisector and linear displacement vectors at all the data-points of the two opposite generating curves. The algorithm has one adjustable parameter that controls the characteristics of transformation of one set of curves from its parents. This unidirectional process has been extended to bi-directional parameterization by superimposing two sets of unidirectional curves generated from both boundary pairs. Case studies show that this algorithm gives reasonably smooth transformation of the boundaries. This algorithm is more robust than the linear Coons method and capable of resolving the 2D boundary-conformed parameterization problems.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.1
• /
• pp.46.1-46.1
• /
• 2019

We present an accurate and efficient method to calculate the gravitational potential of an isolated system in three-dimensional Cartesian and cylindrical coordinates subject to vacuum (open) boundary conditions. Our method consists of two parts: an interior solver and a boundary solver. The interior solver adopts an eigenfunction expansion method together with a tridiagonal matrix solver to solve the Poisson equation subject to the zero boundary condition. The boundary solver employs James's method to calculate the boundary potential due to the screening charges required to keep the zero boundary condition for the interior solver. A full computation of gravitational potential requires running the interior solver twice and the boundary solver once. We develop a method to compute the discrete Green's function in cylindrical coordinates, which is an integral part of the James algorithm to maintain second-order accuracy. We implement our method in the {\tt Athena++} magnetohydrodynamics code, and perform various tests to check that our solver is second-order accurate and exhibits good parallel performance.

• Earthquakes and Structures
• /
• v.9 no.1
• /
• pp.173-194
• /
• 2015

Site response analysis is an important topic in earthquake engineering. A time-domain numerical method called as one-dimensional (1D) finite element artificial boundary method is proposed to simulate the homogeneous plane elastic wave propagation in a layered half space subjected to the obliquely incident plane body wave. In this method, an exact artificial boundary condition combining the absorbing boundary condition with the inputting boundary condition is developed to model the wave absorption and input effects of the truncated half space under layer system. The spatially two-dimensional (2D) problem consisting of the layer system with the artificial boundary condition is transformed equivalently into a 1D one along the vertical direction according to Snell's law. The resulting 1D problem is solved by the finite element method with a new explicit time integration algorithm. The 1D finite element artificial boundary method is verified by analyzing two engineering sites in time domain and by comparing with the frequency-domain transfer matrix method with fast Fourier transform.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.3
• /
• pp.353-361
• /
• 2001

Recently the boundary element method has been developed swiftly. The boundary element method is an efficient and accurate means for analysis of two dimensional elastic crack problems. This paper is concerned with the evaluation and the prediction of the stress intensity factor(SIF) for the crack emanating from the circular hole using boundary element method-neural network. The SIF of the crack emanating from the hole was calculated by using boundary element method. Neural network is used to evaluate and to predict SIF from the results of boundary element method. The organized neural network system (structure of four processing element) was learned with the accuracy 99%. The learned neural network system could be evaluated and predicted with the accuracy of 83.3% and 71.4% (in cases of SIF and virtual SIF). Thus the proposed boundary element method-neural network is very useful to estimate the SIF.