• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.27-37
• /
• 1991

In this paper the nonlinear free surface flows for an axisymmetric submerged body oscillating beneath the free surface are solved and the forces acting on the body are calculated. A boundary integral method is applied to solve the axisymmetric boundary value problem and the Runge-Kutta 4-th order method is used for the time stepping of the free surface location. The nonlinear forces acting on the axisymmetric body are computed and compared with published results.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.8
• /
• pp.2524-2539
• /
• 1996

A design process and an axisymmetric boundary integral equation method for precise structural analysis of the aircraft gas turbine rotor disk were developed. This axisymmetric boundary integral equation method for stress and steady-state thermal analysis was improved in solution accuracy by appling an implicit technique for Cauchy principal value evaluation, a subelement technique for weak singular integral evaluation and a double exponentical integral technoque for internal point solution near boundary surfaces. Stresses, temperatures, low cycle fatigue lifes and critical speeds for the turbine rotor disk of the thrust 1421 N class turbojet engine were analysed in a pratical calculation model problem.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.25-34
• /
• 1993

The long-term behavior, such as in excavation problems of weak medium, can be dealt with by the elasto-viscoplasticity models. In this paper, a combined formulation of elasto-viscoplasticity using boundary elements and finite elements without using internal cells is presented. The domain integral introduced due to the viscoplastic stresses is transformed into a boundary integral applying direct integration in cylindrical coordinates. The results of the developed boundary element analysis are compared with those from the explicit solution and from the finite element analysis. It is observed that the boundary element analysis without internal cells results in some error because of its deficiency in handling the nonlinearity in local stress concentration. Therefore, a coupled analysis of boundary elements and finite elements, in which finite elements are used in the area of stress concentration, is developed. The coupled method is applied to a time dependent inelastic problem with semi-infinite boundaries. It results in reasonable solution compared with other methods where relatively higher degree of freedoms are employed. Thus, it is concluded that the combined analysis may be used for such problems in the effective manner.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.3 no.3
• /
• pp.170-175
• /
• 1991

In order to calculate reflection coefficients of surface waves due to underwater obstacles, the hybrid boundary integral equation method is used. Numerical results based on the linear element are verified to the exact solutions for a flat bottom. Reflection coefficient and transmission coefficients for a step are compared to Young's results and the results by EFEM (Kirby et al., 1987), in which reflection coefficients decrease to zero as the periods of incident waves decrease. Reflection coefficients for a sinusiodal hump located on a constant depth increase due to the interaction between humps.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.9 no.1
• /
• pp.9-16
• /
• 2005

Jet vane is an useful component which is installed at the end of a nozzle for the purpose of the posture control and the secure controlling stability during the initial launching of a rocket. During several seconds from its initial launching moment, the JV driving part is heated due to the direct contact of the vane with the combusted gas and the vane is ablated mechanically or chemically. In this study, as the fundamental study for the thermal analysis of jet vane, the heat transfer into a jet vane which is located in the uniform supersonic flow field is calculated. For this, boundary layer integral method and finite difference method are used simultaneously. Based on the thermal boundary conditions derived from the analysis, the transient heat conduction in the vane is also calculated.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.13 no.3
• /
• pp.21-28
• /
• 2009

The precise analysis of soil-pile-structure interaction requires a proper description of soil layer, pile, and structure. In commonly used finite element simulations, mesh boundaries should match the material discontinuity line. However, in practice, the geometry of soil profiles and piles may be so complex that mesh alignment becomes a wasteful and difficult task. To overcome these difficulties, a different integration method is adopted in this paper, which enables easy integration over a regular element with material discontinuity regardless of the location of the discontinuity line. By applying this integration method, the mesh can be generated rapidly and in a highly structured manner, leading to a very regular stiffness matrix. The influence of the shape of the soil profile and piles on the response is examined, and the validity of the proposed soil-pile structure interaction analysis method is demonstrated through several examples. It is seen that the proposed analysis method can be easily used on soil-pile-structure interaction problems with complex interfaces between materials to produce reliable results regardless of the material discontinuity line.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.4
• /
• pp.775-786
• /
• 2002

A Mixed Volume and Boundary Integral Equation Method is applied for the effective analysis of elastic wave scattering problems and 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. In the formulation of this method, the continuity condition at each interface is automatically satisfied, and in contrast to finite element methods, where the full domain needs to be discretized, this method requires discretization of the inclusions only. Finally, this method takes full advantage of the pre- and post-processing capabilities developed in FEM and BIEM. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with orthotropic inclusions and voids or isotropic inclusions, and the analysis of plane wave scattering problems in unbounded isotropic matrix with isotropic inclusions and voids, it will be established that this new method is very accurate and effective for solving plane wave scattering problems and plane elastic problems in unbounded solids containing general anisotropic inclusions and voids/cracks or isotropic inclusions.

• Composites Research
• /
• v.23 no.4
• /
• pp.7-13
• /
• 2010

A volume integral equation method (VIEM) is introduced for the solution of elastostatic problems in an unbounded isotropic elastic solids containing interacting multiple anisotropic inclusions subject to remote uniaxial tension. The method is applied to two-dimensional problems involving long parallel cylindrical inclusions. A detailed analysis of stress field at the interface between the matrix and the central inclusion is carried out for square and hexagonal packing of the inclusions. Effects of the number of anisotropic inclusions and various fiber volume fractions on the stress field at the interface between the matrix and the central inclusion are also investigated in detail. The accuracy of the method is validated by solving the single inclusion problem for which solutions are available in the literature.

• Geophysics and Geophysical Exploration
• /
• v.27 no.3
• /
• pp.181-195
• /
• 2024

This tutorial explains that the static effect in the magnetotelluric (MT) survey is a physical phenomenon caused by charges accumulated on the boundaries of subsurface inhomogeneities. To facilitate understanding of the physical phenomenon, differences between static induction and charge accumulation on the boundary are explained and analyzed with help of schematic illustrations. Subsequently, from the electromagnetic (EM) integral equation formulation, it is clearly shown that the secondary electric field due to charges accumulated on the interface in the presence of the primary field appears as the static effect. Therefore, except in the cases of the layered earth or a two-dimensional earth with transverse magnetic (TM) mode excitation, the static effect always exists in MT responses and further, it is not 'static' but rather frequency dependent. Despite the fact that the static effect is a secondary electric field due to inhomogeneity, inevitable under-sampling in the frequency and spatial domains prevent the effect from being handled properly in numerical inversion. Therefore, considering the practical aspects of the MT survey, which cannot be a continuous measurement covering the entire survey area over a wide frequency band, a three-dimensional (3-D) inversion incorporating the static shift as a constraint with the Gaussian distribution is introduced. To enhance understanding of the integral equation EM modeling, the formulation of the 3-D integral equation and mathematical analyses of the Green tensor and scattering current are described in detail in the Appendix.

• Geophysics and Geophysical Exploration
• /
• v.27 no.1
• /
• pp.51-56
• /
• 2024

In this paper, the vector gravity and gravity gradient tensor of an elliptical disk are derived. The vector gravity of an elliptical disk is defined by differentiating the gravitational potential due to the elliptical disk expressed by a double integral with respect to each axial direction. The vector gravity defined by the double integral is then transformed into a line integral of a closed curve along the elliptical disk boundary using the complex Green's theorem. Finally, vector gravity due to the elliptical disk is derived by 1D parametric numerical integration along the elliptical disk boundary. The xz, yz, zz components of the gravity gradient tensor due to the elliptical disk are obtained by differentiating the vector gravity with respect to vertical direction. The xx, yy, xy components are derived by differentiating the horizontal components of the vector gravity in the form of a double integral with respect to horizontal directions and then using the complex Green's theorem.