• 대한조선학회논문집
• /
• 제39권4호
• /
• pp.1-10
• /
• 2002

초기쇄파의 수치모사에는 지금까지 경계적분법이 주로 쓰여왔고, 이 방법은 과도한 계산시간의 문제를 제외하고는 어느 정도 성공적이라고 할 수 있다. 본 논문에서는 쇄파실험을 수치모사하기 위한 새로운 수치기법을 보였다. 이 수치기법은 고차 스펙트럴/경계요소법과 경계적분법을 순차적으로 사용하며, 계산시간을 현저히 줄여준다. 조파 및 파 에너지 집중과정은 고차 스펙트럴/경계요소법에 의해 효율적으로 수치모사되고, 파의 전복과정만이 경계적분법에 의해 계산된다. 계산예에서 높은 입자속도와 가속도 등 쇄파의 두드러진 특성이 보여졌다.

• 한국정밀공학회지
• /
• 제20권8호
• /
• pp.175-184
• /
• 2003

A method of the shape design sensitivity analysis based on the boundary integral equation formulation is presented for two-dimensional inhomogeneous thermal conducting solids with multiple domains. Shape variation of the external and interface boundary is considered. A sensitivity formula of a general performance functional is derived by taking the material derivative to the boundary integral identity and by introducing an adjoint system. In numerical analysis, state variables of the primal and adjoint systems are solved by the boundary element method using quadratic elements. Two numerical examples of a compound cylinder and a thermal diffuser are taken to show implementation of the shape design sensitivity analysis. Accuracy of the present method is verified by comparing analyzed sensitivities with those by the finite difference. As application to the shape optimization, an optimal shape of the thermal diffuser is found by incorporating the sensitivity analysis algorithm in an optimization program.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1994년도 봄 학술발표회 논문집
• /
• pp.128-136
• /
• 1994

In this study, dynamic boundary element method using mass matrix is derived, using fundamental solutions for the semi-infinite domain. In constituting boundary integral equations for the dynamic equilibrium condition, inertia term in the form of domain integral is transformed into boundary integral form. Corresponding system equations are derived, and a boundary element program is developed. In addition, equations for free vibration is formulated, and eigenvalue analysis is performed. The results from the dynamic boundary element analysis for a tunnel problem are compared with those from the finite element analysis. According to the comparison, boundary element method using mass matrix is consistent with the results of finite element method. Consequently, in tunnel vibration problems, it results in reasonable solution compared with other methods where relatively higher degree of freedoms are employed.

• Kyungpook Mathematical Journal
• /
• 제62권3호
• /
• pp.557-581
• /
• 2022

The anisotropic-diffusion convection equation with exponentially variable coefficients is discussed in this paper. Numerical solutions are found using a combined Laplace transform and boundary element method. The variable coefficients equation is usually used to model problems of functionally graded media. First the variable coefficients equation is transformed to a constant coefficients equation. The constant coefficients equation is then Laplace-transformed so that the time variable vanishes. The Laplace-transformed equation is consequently written as a boundary integral equation which involves a time-free fundamental solution. The boundary integral equation is therefore employed to find numerical solutions using a standard boundary element method. Finally the results obtained are inversely transformed numerically using the Stehfest formula to get solutions in the time variable. The combined Laplace transform and boundary element method are easy to implement and accurate for solving unsteady problems of anisotropic exponentially graded media governed by the diffusion convection equation.

• 대한수학회논문집
• /
• 제34권3호
• /
• pp.1015-1027
• /
• 2019

A class of systems of singularly perturbed convection diffusion type equations with integral boundary conditions is considered. A numerical method based on a finite difference scheme on a Shishkin mesh is presented. The suggested method is of almost first order convergence. An error estimate is derived in the discrete maximum norm. Numerical examples are presented to validate the theoretical estimates.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집A
• /
• pp.341-348
• /
• 2001

A Mixed Volume and Boundary Integral Equation Method is applied for the effective analysis of plane elastostatic problems in unbounded solids containing general anisotropic inclusions and voids or isotropic inclusions. 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. This new method can also be applied to general two-dimensional elastodynamic and elastostatic problems with arbitrary shapes and number of anisotropic inclusions and voids or isotropic inclusions. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with orthotropic inclusions and voids or isotropic inclusions, it will be established that this new method is very accurate and effective for solving plane elastic problems in unbounded solids containing general anisotropic inclusions and voids or isotropic inclusions.

• Structural Engineering and Mechanics
• /
• 제15권5호
• /
• pp.535-550
• /
• 2003

In this paper, a boundary-type meshfree method, the boundary radial point interpolation method (BRPIM), is presented for solving boundary value problems of two-dimensional solid mechanics. In the BRPIM, the boundary of a problem domain is represented by a set of properly scattered nodes. A technique is proposed to construct shape functions using radial functions as basis functions. The shape functions so formulated are proven to possess both delta function property and partitions of unity property. Boundary conditions can be easily implemented as in the conventional Boundary Element Method (BEM). The Boundary Integral Equation (BIE) for 2-D elastostatics is discretized using the radial basis point interpolation. Some important parameters on the performance of the BRPIM are investigated thoroughly. Validity and efficiency of the present BRPIM are demonstrated through a number of numerical examples.

• 한국해양공학회지
• /
• 제32권6호
• /
• pp.447-457
• /
• 2018

In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

• 충청수학회지
• /
• 제7권1호
• /
• pp.141-152
• /
• 1994

We consider a numerical scheme for solving a nonlinear boundary integral equation (BIE) obtained by reformulation the nonlinear boundary value problem (BVP). We give a simple alternative to the standard collocation method for the nonlinear BIE. This method consists of one conventional linear system and another coupled linear system resulting from an auxiliary BIE which is obtained by differentiating both side of the nonlinear interior integral equations. We obtain an analogue BIE through the perturbation of the fundamental solution of Laplace's equation. We procure the super-convergence of approximate solutions.

• 소음진동
• /
• 제7권6호
• /
• pp.975-984
• /
• 1997

A direct boundary element method (DBEM) is developed for thin bodies whose surfaces are rigid or compliant. The Helmholtz integral equation and its normal derivative integral equation are adoped simultaneously to calculate the pressure on both sides of the thin body, instead of the jump values across it, to account for the different surface conditions of each side. Unlike the usual assumption, the normal velocity is assumed to be discontinuous across the thin body. In this approach, only the neutral surface of the thin body has to be discretized. The method is validated by comparison with analytic and/or numerical results for acoustic scattering and radiation from several surface conditions of the thin body; the surfaces are rigid when stationary or vibrating, and part of the interior surface is lined with a sound-absoring material.