• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.3
• /
• pp.624-635
• /
• 1990

Analytical study based on linear fracture mechanics was conducted on propagation behavior of inclined surface crack in semi-infinite elastic body. The analytical model was assumed to be inclined surface crack under plane strain condition upon which Hertzian stress was superimposed. Supposing continuous distribution of dislocation and applying Erdogan-Gupta's method to this crack problem, the stress intensity factors $K_{I}$ and $K_{II}$) at the crack-tip were obtained for various Hertzian contact positions. Analytic results have shown that driving force for crack growth is $K_{I}$ for non-lubricated condition and $K_{II}$ for fluid and boundary lubricated condition. The coefficient of friction at the hertzian contact and crack surfaces plays an important role in predicting the direction of crack propagation. It is also found that the maximum effective stress intensity factor exists at cracks of a certain specific length depending on lubricated condition.ion.n.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.4
• /
• pp.10-17
• /
• 2002

A study on the dynamic ground effect on three-dimensional wings is done using an indirect boundary element method(unsteady panel method). An integral equation is obtained by applying Green's theorem on all surfaces of the fluid domain. Constant strength dipole and source panels arc distributed on a wing's surface. The wake sheet is represented by constant strength dipoles. At each time step, a row of wake panels is assumed to be convected from the trailing edge of the wing. The tip vortex behind wings in dynamic ground effect moves outward. The amplitudes of the aerodynamic coefficients for the wings in dynamic ground effect are augmented much more comparing to the case in static ground effect.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.3
• /
• pp.258-268
• /
• 2008

The radiation potential of a ship advancing in waves is studied using the 3D time-domain forward-speed free-surface Green function and the Green integral equation. Numerical solutions are obtained by making use of the 2nd order BEM(Boundary Element Method) which make it possible to take account of the line integral along the waterline in a rigorous manner. The 6 degree of freedom motion memory functions of a hemisphere and the Wigley seakeeping model obtained by direct integration of the time-domain 3D potentials over the wetted surface are presented for various Froude numbers.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.13 no.1
• /
• pp.47-52
• /
• 2010

The wave energy absorption efficiency and the first-order and the time-mean second-order wave loads of a three-dimensional bottom-mounted oscillating water column (OWC) chamber structure are studied. The potential problem is solved by making use of a hybrid Green integral equation associated with the finite-waterdepth free-surface Green function outside a twin chamber and the Rankine Green function inside taking account of the fluctuating air pressure inside the chamber. Numerical results of the primary wave energy converting efficiency and the oscillating and steady wave loads of a three-dimensional bottom-mounted OWC pilot plant have been presented.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.6
• /
• pp.1621-1628
• /
• 1990

Numerical solutions are presented for the heat transfer from radiating and convecting fins. Consideration is given to thin, annular fins attached to a tube surface for which the temperature is constant. Fin to fin, fin to base, and fin to environment radiative interactions are considered. It is assumed that the radiating surface is diffuse-gray, the environment is black, and the surrounding fluid is transparent. The radiation terms are formulated by using Poljak's net-radiation methoad. The mathematical description of the simultaneously heat transport by conduction, convection, and radiation leads to a nonlinear integro-differential equation. This has been solved for a wide range of the pertinent physical parameters by using finite difference method and iteration method based on the Newton-Raphson technique. The temperature distributions, heat transfer rates, fin efficiencies, and fin effectivenesses are presented in dimensionless form. The results definitely indicate that the use of fins leads to a significant increase in heat transfer compared with the unfinned tube.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.31 no.2
• /
• pp.202-209
• /
• 1998

A three-dimensional numerical method based on the Green's integral equation is developed to predict the motion response and drift force in waves for the ocean monitoring facilities. In this method, we use source and doublet distribution, and triangular and rectangular eliments. To eliminate the irregular frequency phenomenon, the method of improved integral equation is applied and the time-mean drift force is calculated by the method of direct pressure integration over the body surface. To conform the validity of the present numerical method, some calculations for the floating sphere are performed and it is shown that the present method provides sufficiently reliable results. As a calculation example for the real facilities, the motion response and the drift force of the vertical cylinder type ocean monitoring buoy with 2.6 m diameter and 3,77 m draft are calculated and discussed. The obtained results of motion response can be used to determine the shape and dimension of the buoy to reduce the motion response, and other data such as the effect of motion reduction due to a damper can be predictable through these motion calculations. Also, the calculation results of drift force can be used in the design procedure of mooring system to predict the maximum wave load acting on the mooring system. The present method has, in principle, no restriction in the application to the arbitrary shape facilities. So, this method can be a robust tool for the design, installation, and operation of various kinds of the floating-type ocean monitoring facilities.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1145-1158
• /
• 1994

A three-dimensional, finite difference, numerical model with free surface was developed on ${\sigma}$-coordinate. A semi-implicit numerical scheme in time has been adopted for computational efficiency. The scheme is essentially independent of the stringent stability criteria (CFL condition) for explicit schemes of external surface gravity wave. Implicit algorithm was applied for vertical shear stress, Coriolis force and pressure gradient terms. The reliability of the model with vertically variable grid system was checked by the comparison of simulation results with analytic solution of wind-driven currents in a one-dimensional channel. Sensitivity analysis of differencing parameters was carried out by applying the model to the calculation of wind-driven currents in a square lake.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.6
• /
• pp.583-590
• /
• 2016

This paper presents a dynamic crack propagation algorithm with Rayleigh damping effect based on the MLS(Moving Least Squares) Difference Method. Dynamic equilibrium equation and constitutive equation are derived by considering Rayliegh damping and governing equations are discretized by the MLS derivative approximation; the proportional damping, which has not been properly treated in the conventional strong formulations, was implemented in both the equilibrium equation and constitutive equation. Dynamic equilibrium equation including time relevant terms is integrated by the Central Difference Method and the discrete equations are simplified by lagging the velocity one step behind. A geometrical feature of crack is modeled by imposing the traction-free condition onto the nodes placed at crack surfaces and the effect of movement and addition of the nodes at every time step due to crack growth is appropriately reflected on the construction of total system. The robustness of the proposed numerical algorithm was proved by simulating single and multiple crack growth problems and the effect of proportional damping on the dynamic crack propagation analysis was effectively demonstrated.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.14 no.5
• /
• pp.458-465
• /
• 2003

In this paper, we analyze the transient electromagnetic response from three-dimensional(3-D) dielectric bodies using a time-domain electric field integral equation formulation. The solution method in this paper is based on the Galerkin＇s method that involves separate spatial and temporal testing procedures. Triangular patch basis functions are used for spatial expansion and testing functions for arbitrarily shaped 3-D dielectric structures. The time-domain unknown coefficients of the equivalent electric and magnetic currents are approximated as an orthonormal basis function set that is derived from the Laguerre functions. These basis functions are also used as the temporal testing. Numerical results involving equivalent currents and far fields computed by the proposed method are presented.

• Journal of Ocean Engineering and Technology
• /
• v.16 no.3
• /
• pp.19-27
• /
• 2002

The motion of a floating OWC chamber in waves is studied taking account of fluctuating air pressure in the air chamber. An atmospheric pressure drop occurs across the upper opening of the chamber which causes not only hydrodynamic but also pneumatic added mass and damping forces to the floating chamber. A velocity potential in the water due to the free surface oscillating pressure patch is added to the conventional radiation-diffraction potential problem. the potential problem inside the chamber is formulated by making use of the Green integral equation associated with the Rankine Green function wile the outer problem with the Kelvin Green function. The two integral equations are solved simultaneously by making use of a matching boundary condition at the lower opening of the chamber to the outer water region. The chamber motion in the frequency domain is calculated for various values of parameters related to the atmospheric pressure drop. The present methods can also be sued for the analysis of air-cushion vehicle motion as well as for the design of a floating OWC wave energy absorber.