• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제2권2호
• /
• pp.81-95
• /
• 1998

The superconvergence of the Sloan iterate obtained from a Galerkin method for the approximate solution of the singular integral equation based on the use of two sets of orthogonal polynomials is investigated. The discrete Sloan iterate using Gaussian quadrature to evaluate the integrals in the equation becomes the Nystr$\ddot{o}$m approximation obtained by the same rules. Consequently, it is impossible to expect the faster convergence of the Sloan iterate than the discrete Galerkin approximation in practice.

• Journal of applied mathematics & informatics
• /
• 제6권2호
• /
• pp.503-512
• /
• 1999

The purpose of this paper is to introduce the (Toeplitz) quadrature method for solving fredholm integral equations of the second kind with mildly singular kernels. We are presented some numerical examples for the computation of the error estimate using the MathCad package.

• East Asian mathematical journal
• /
• 제17권1호
• /
• pp.33-45
• /
• 2001

In this paper an improved version of a fixed point theorem of the present author [3] in Banach algebras is obtained under the weaker conditions with a different method and using measure of non-compactness. The newly developed fixed point theorem is further-applied to certain nonlinear integral equations of mixed type for proving the existence theorems and stability of the solution in Banach algebras.

• 한국전산유체공학회지
• /
• 제3권1호
• /
• pp.11-21
• /
• 1998

A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition. The velocity field is obtained by using the Biot-Savart integral. The Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well established for potential flow analysis. The present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).

• 대한전기학회논문지:전기물성ㆍ응용부문C
• /
• 제51권11호
• /
• pp.551-558
• /
• 2002

A time-domain combined field integral equation (CFIE) is presented to obtain the transient scattering response from arbitrarily shaped three-dimensional conducting bodies. This formulation is based on a linear combination of the time-domain electric field integral equation (EFIE) with the magnetic field integral equation (MFIE). The time derivative of the magnetic vector potential in EFIE is approximated using a central finite difference approximation and the scalar potential is averaged over time. The time-domain CFIE approach produces results that are accurate and stable when solving for transient scattering responses from conducting objects. The incident spectrum of the field may contain frequency components, which correspond to the internal resonance of the structure. For the numerical solution, we consider both the explicit and implicit scheme and use two different kinds of Gaussian pulses, which may contain frequencies corresponding to the internal resonance. Numerical results for the EFIE, MFIE, and CFIE are presented and compared with those obtained from the inverse discrete Fourier transform (IDFT) of the frequency-domain CFIE solution.

• 한국정밀공학회지
• /
• 제16권4호통권97호
• /
• pp.197-204
• /
• 1999

The boundary element method was used for studying singularities of an interface crack with contact zones. The iterative procedure is applied to estimate the contact zone size. Because the contact zone size was extremely small in a tension field, a large number of Gaussian points were used for numerical integration of the Kernels. Stress extrapolation method and J-integral were used ofr determining stress intensity factors. When the interface crack was assumed to have opened tips, oscillatory singularities appear near the tips of the interface crack. But the interface crack with contact zone which Comninou suggested had no oscillatory behavior. The contact zone size under shear loading was much larger than that under tensile. The stress intensity factors computed by stress extrapolation method were close to those of Comninou's solution. And the stress intensity factor evaluated by J-integral was similar to that by stress extrapolation method.

• 대한기계학회논문집
• /
• 제19권4호
• /
• pp.981-989
• /
• 1995

A new weight function approach to determine SIF(stress intensity factor) using single-layer potential has been presented. The crack surface displacement field was represented by one boundary integral term whose kernel was modified from Kelvin's fundamental solution. The proposed method enables the calculation of SIF using only one SIF solution without any modification for the crack geometries symmetric in two-dimensional plane such as a center crack in a plate with or without an internal hole, double edge cracks, circumferential crack or radial cracks in a pipe. The application procedure to those crack problems is very simple and straightforward with only one SIF solution. The necessary information in the analysis is two reference SIFs. The analysis results using present closed-form solution were in good agreement with those of the literature.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1994년도 가을 학술발표회 논문집
• /
• pp.38-45
• /
• 1994

P-version finite element model for the computation of stress intensity factors in two dimensional cracked panels by J-integral method is presented. The proposed model is based on high order theory and hierarchical shape function. The displacements fields are defined by integrals of Legendre polynomials which can be classified into three part such as basic mode, side mode, integral mode. The stress intensity factors are computed by J-integral method. The example models for validating the proposed p-version model are centrally cracked panel, single and double edged crack in a rectangular panel under pure Mode I. And the analysis results are compared with those by the h-version of FEM and empirical solutions in literatures. Very good agreement with the existing solution are shown.

• 대한전기학회논문지:전기물성ㆍ응용부문C
• /
• 제51권11호
• /
• pp.559-566
• /
• 2002

In this paper, we present an accurate and stable method for the solution of the transient electromagnetic response from the conducting wire structures using the time domain integral equation. By using an implicit scheme with the central finite difference approximation for the time domain electric field integral equation, we obtain the transient response from a wire scatterer illuminated by a plane wave and a conducting wire antenna with an impressed voltage source. Also, we consider a wire above a 3-dimensional conducting object. Numerical results are presented, which show the validity of the presented methodology, and compared with a conventional method using backward finite difference approximation.

• 대한기계학회논문집
• /
• 제18권1호
• /
• pp.121-126
• /
• 1994

A mutual integral, which has the conservation property, is applied to the problem of a crack in an isotropic elastic material. The stress intensity factors $K_{I}, K_{II}, K_{III}$ and T-stress for the problem in an infinite medium are easily obtained by using the mutual integral without solving the boundary value problem. The auxiliary solutions necessary in the proposed method are taken from the known asymptotic solutions. This method is amenable to numerical evaluation of the stress intensity factors and T-stress if the crack in a finite medium is considered.