• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.4
• /
• pp.1-10
• /
• 2002

For the time-domain simulation of incipient breaking waves, usually the boundary integral method has been used so far, and it seems to be successful except a problem of too much computation time. The present paper shows a new computation technique for the simulation of breaking wave experiment. This technique uses the high-order spectral/boundary element method and the boundary integral method in sequence, and reduces the computation time remarkably. The wave generation and energy focusing process is efficiently simulated by the high-order spectral/boundary element method. Only the wave over-turning process is simulated by the boundary integral method. In the example calculation result, salient features of breaking waves such as high particle velocities and accelerations are shown.

• Journal of the Korean Society for Precision Engineering
• /
• v.20 no.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.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.04a
• /
• 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
• /
• v.62 no.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.

• Communications of the Korean Mathematical Society
• /
• v.34 no.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.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• 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
• /
• v.15 no.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.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.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.

• Journal of the Chungcheong Mathematical Society
• /
• v.7 no.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.

• Journal of KSNVE
• /
• v.7 no.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.