• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.3
• /
• pp.151-157
• /
• 2021

This paper is devoted to a convergence study of the nonlocal integral operator in peridynamics. The implicit formulation can be an efficient approach to obtain the static/quasi-static solution of crack propagation problems. Implicit methods require constly large-matrix operations. Therefore, convergence is important for improving computational efficiency. When the radial influence function is utilized in the nonlocal integral equation, the fractional Laplacian integral equation is obtained. It has been mathematically proved that the condition number of the system matrix is affected by the order of the radial influence function and nonlocal horizon size. We formulate the static crack problem with peridynamics and utilize Newton-Raphson methods with a preconditioned conjugate gradient scheme to solve this nonlinear stationary system. The convergence behavior and the computational time for solving the implicit algebraic system have been studied with respect to the order of the radial influence function and nonlocal horizon size.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.19 no.4
• /
• pp.427-435
• /
• 2008

In this paper, a stable marching-on in time(MOT) method with 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 derivatives in the EFIE and MFIE are approximated using a central finite difference scheme and other terms are averaged over time. This time domain CFIE approach produces results that are accurate and stable when solving for transient scattering responses from conducting objects. Numerical results with the proposed MOT scheme are presented and compared with those obtained from the conventional method and the inverse discrete Fourier transform(IDFT) of the frequency domain CFIE solution.

• Journal of KSNVE
• /
• v.6 no.5
• /
• pp.607-616
• /
• 1996

A combined computational fluid dynamics(CFD)-Kirchhoff method is presented for predicting high-speed impulsive noise generated by a hovering blade. Two types of Kirchhoff integral formula are used; one for the classical linear Kirchhoff formulation and the other for the nonlinear Kirchhoff formulation. An Euler finite difference solver is solved first to obtain the flow field close to the blade, and then this flow field is used as an input to a Kirchhoff formulation to predict the acoustic far-field. These formulas are used at Mach numbers of 0.90 and 0.95 to investigate the effectiveness of the linear and nonlinear Kirchhoff formulas for delocalized flow. During these calculiations, the retarded time equation is also carefully examined, in particular, for the cases of the control surface located outside of the sonic cylinder, where multiple roots are obtained. Predicted results of acoustic far-field pressure with the linear Kirchhoff formulation agree well with experimental data when the control surface is at the certain location(R=1.46), but the correlation is getting worse before or after this specific location of the control surface due to the delocalized nonlinear aerodynamic flow field. Calculations based on the nonlinear Kirchhoff equation using a linear sonic cylinder as a control surface show a reasonable agreement with experimental data in negative amplitudes for both tip Mach numbers of 0.90 and 0.95, except some computational integration problems over a shock. This concliudes that a nonlinear formulation is necessary if the control surface is close to the blade and the flow is delocalized.

• Structural Engineering and Mechanics
• /
• v.48 no.2
• /
• pp.241-255
• /
• 2013

The receding contact problem for two elastic layers whose elastic constants and heights are different supported by two elastic quarter planes is considered. The lower layer is supported by two elastic quarter planes and the upper elastic layer is subjected to symmetrical distributed load whose heights are 2a on its top surface. It is assumed that the contact between all surfaces is frictionless and the effect of gravity force is neglected. The problem is formulated and solved by using Theory of Elasticity and Integral Transform Technique. The problem is reduced to a system of singular integral equations in which contact pressures are the unknown functions by using integral transform technique and boundary conditions of the problem. Stresses and displacements are expressed depending on the contact pressures using Fourier and Mellin formula technique. The singular integral equation is solved numerically by using Gauss-Jacobi integration formulation. Numerical results are obtained for various dimensionless quantities for the contact pressures and the contact areas are presented in graphics and tables.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.5
• /
• pp.1458-1467
• /
• 1996

A direct differentiationmethod is presented for the shape design sensitivity analysis of axisymmeetric elastic solids. Based on the exisymmetric boundary integralequaiton formulation, a new boundary ntegral equatio for sensitivity analysis is derived by taking meterial derivative to the same integral identity that was used in the adjoint variable melthod. Numerical implementation is performed to show the applicaiton of the theoretical formulation. For a simple example with analytic solution, the sensitivities by present method are compared with analytic sensitivities. As an application to the shape optimization, an optimal shape of a gas turbine disc toinimize the weight under stress constraints is found by incorporating the sensitivity analysis algorithm in an optimizatio program.

• Journal of electromagnetic engineering and science
• /
• v.12 no.4
• /
• pp.260-270
• /
• 2012

This work aims to present a theoretical analysis of the electric and magnetic surface current densities of a circular microstrip antenna (CMSA) as a body of revolution. The rigorous analysis of these problems begins with the application of the equivalence principle, which introduces an unknown electric current density on the conducting surface and both unknown equivalent electric and magnetic surface current densities on the dielectric surface. These current densities satisfy the integral equations (IEs) obtained by canceling the tangential components of the electric field on the conducting surface and enforcing the continuity of the tangential components of the fields across the dielectric surface. The formulation of the radiation problems is based on the combined field integral equation. This formulation is coupled with the method of moments (MoMs) as a numerical solution for this equation. The numerical results of the electric and magnetic surface current densities on the outside boundary of a CMSA excited by $TM_{11^-}$ and $TM_{21^-}$ modes are presented. The radiation pattern is calculated numerically in the two principle planes for a CMSA and gives a good results compared with measured results published by other research workers.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.15 no.8
• /
• pp.765-774
• /
• 2004

In this paper, a novel approximate solution for scattering by a very thin planar homogeneous dielectric scatterer with an arbitrary shape is formulated. This solution is based on a volumetric integral equation and is expressed in terms of Fourier transform. It is shown that the obtained solution is reduced to an exact solution for an infinite dielectric slab. For 2D, or 3D scatterers, the formulation is verified numerically. Especially fur edge-on TM polarized wave incidence a closed-form solution of backscattering from a thin dielectric half-plane is formulated, which is very accurate for wide range of normalized surface impedance except very low impedances(│η│〈0.5).

• Nuclear Engineering and Technology
• /
• v.24 no.3
• /
• pp.319-328
• /
• 1992

The Galerkin formulation of the finite element method is applied to the integral law of the first-order form of the one-group neutron transport equation in one-dimensional spherical geometry. Piecewise linear or quadratic Lagrange polynomials are utilized in the integral law for the angular flux to establish a set of linear algebraic equations. Numerical analyses are performed for the scalar flux distribution in a heterogeneous sphere as well as for the criticality problem in a uniform sphere. For the criticality problems in the uniform sphere, the results of the finite element method, with the use of continuous finite elements in space and angle, are compared with the exact solutions. In the heterogeneous problem, the scalar flux distribution obtained by using discontinuous angular and spatical finite elements is in good agreement with that from the ANISN code calculation.

• Journal of electromagnetic engineering and science
• /
• v.19 no.1
• /
• pp.6-12
• /
• 2019

An IE-FFT algorithm is implemented and applied to the electromagnetic (EM) solution of perfect electric conducting (PEC) scattering problems. The solution of the method of moments (MoM), based on the magnetic field integral equation (MFIE), is obtained for PEC objects with closed surfaces. The IE-FFT algorithm uses a uniform Cartesian grid to apply a global fast Fourier transform (FFT), which leads to significantly reduce memory requirement and speed up CPU with an iterative solver. The IE-FFT algorithm utilizes two discretizations, one for the unknown induced surface current on the planar triangular patches of 3D arbitrary geometries and the other on a uniform Cartesian grid for interpolating the free-space Green's function. The uniform interpolation of the Green's functions allows for a global FFT for far-field interaction terms, and the near-field interaction terms should be adequately corrected. A 3D block-Toeplitz structure for the Lagrangian interpolation of the Green's function is proposed. The MFIE formulation with the IE-FFT algorithm, without the help of a preconditioner, is converged in certain iterations with a generalized minimal residual (GMRES) method. The complexity of the IE-FFT is found to be approximately $O(N^{1.5})$and $O(N^{1.5}logN)$ for memory requirements and CPU time, respectively.

• Journal of the Institute of Electronics Engineers of Korea TC
• /
• v.39 no.4
• /
• pp.201-208
• /
• 2002

In this paper, we present a new formulation using time-domain electric field integral equation (TD-EFIE) to obtain transient scattering response from arbitrarily shaped three-dimensional conducting bodies. The time derivative of the magnetic vector potential is approximated with a central finite difference and the scalar potential is time averaged by dividing it into two terms. This approach with an implicit method using central difference results in accurate and more stable transient scattering responses from conducting objects. Detailed mathematical steps are included and several numerical results are presented and compared with the inverse discrete Fourier transform (IDFT) of the frequency-domain solution.