• Journal of Advanced Marine Engineering and Technology
• /
• v.12 no.2
• /
• pp.21-34
• /
• 1988

Machinery enclosures are widely adopted to reduce the noise emission in various fields of application. Emitted noise, which is due to the vibration of enclosure's outer surface, is composed of two kinds of sound with different path of propagation. One is the "structure-borne sound", while the other is "air-borne sound". In order to get a most efficient machinery enclouser a prudent consideration upon the above structure-borne and air-borne sound is required, as the guiding principle of contermeasure for each noise is quite different. The controlling of input vibration and its isolation are major subjects for the structure-borne sound, and the specifications of absorbing members and damping panels are the major related matters for the air-borne sound. Hence, it seems very efficient to separate the total sounds into two categories with a great accuracy when one think of further reduction of noise from the existing enclosure, although its separating methods have not been made clear for many years. Author proposes an application method of experimental modal analysis to extract the structure-borne sound from the measured total radiation sound, as the air-borne sound is deduced by the vectorial difference between the measured total radiation sound and the calculated structure-borne sound. In order to calculate the correct structure-borne sound by the excitation at an arbitrary point on the enclosure structure, it is important to decide 1) how to estimate the enclosure's surface vibration velocity and 2) how to compute the radiation sound which is considered as the effect of vibration modes of enclosure surface. The former can be solved with total frequency response function calculated by the application of experimental modal analysis. The latter is to be solved by the author's new approaches for radiation sound computation by means of the Rayleigh's integral equation and the boundary-element method applied complex surface vibration velocity. As a first step, structure-borne sound by the excitation at an arbitry point on the rectangular plate with fixed edges, has been calculated to verified the reliability of the developed computation methods. The results of calculation show good agreements with those of the actual measurements.actual measurements.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.26 no.1
• /
• pp.94-100
• /
• 2015

Coronal mass ejection (CME) released due to solar flare explosion cause geomagnetic disturbance. The induced current by massive geomagnetic disturbance can cause damage to the transformer. The calculated geoelectric field is a major parameter of the geomagnetically induced current (GIC). The method applying a Fourier transform has a high accuracy but it needs all data measured for 24 hours. And the other method applying a integral equation can calculate in real time but it requires to check an accuracy. To reduce the gap between the calculated results of two methods, it adjusts the integration section. As a result, the correlation between two calculated geoelectric fields is high, and the event time and direction of the calculated current is the same as that of the measured current, and it's accuracy rate is above 92 percent.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.4
• /
• pp.15-26
• /
• 1990

This paper describes a potential-baoed panel method formulated for the analysis cf a supercavitating two-dimensional symmetri strut. The method employs normal dipoles and sources distributed on the foil and cavity surfaces to represent the potential flow around the cavitating hydrofoil. The kinematic boundary condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the fictitious inner flow region of the foil, and the dynamic boundary condition on the cavity surface is satisfied by requiring that the potential vary linearly, i.e., the tangential velocity be constant. Green's theorem then results in a potential-based integral equation rather than the usual velocity-based formulation of Hess & Smith type, With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with improved accuracy compared to those of the linearized lifting surface theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. This iteration process is repeated until the results are converged.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.18 no.9
• /
• pp.1023-1029
• /
• 2007

In this paper, we propose the method to solve the BLT equations using Laguerre polynomials in time domain. The solution of BLT equations is obtained by recursive, differential and integral properties of Laguerre polynomials. The verification of the proposed method is tested by applying it to the two-wired transmission line with resistors and capacitors, which is illuminated by the electromagnetic plane wave pulse. And the result is compared with the corresponding transient responses obtained from inverse fast Fourier transform(IFFT) of the frequency domain solutions of BLT equations.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.154-160
• /
• 2007

This paper compares the characteristics of three different algorithms for three-dimensional (3D) magnetotelluric (MT) modeling. These methods are developed by Mackie et al. (1994), Sasaki (1999) and Nam et al. (2007). The first and second methods are based on the finite difference method (FDM), while the last one the finite-element method (FEM). MT responses, apparent resistivities and phases, for a COMMEMI 3D-2 model show a good agreement with integral equation solutions and only minor discrepancies are found over the anomalous bodies in the 3D model. The computation time of the two methods based on FDM is short and the static divergence correction contributes to speed up. The FEM modeling scheme is accurate but slow.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.123-130
• /
• 2002

A general formulation for shape design sensitivity analysis over three dimensional beam structure is developed based on a variational formulation of the beam in linear elasticity. Sensitivity formula is derived based on variational equations in cartesian coordinates using the material derivative concept and adjoint variable method for the displacement and Von-Mises stress functionals. Shape variation is considered for the beam shape in general 3-dimensional direction as well as for the orientation angle of the beam cross section. In the sensitivity expression, the end points evaluation at each beam segment is added to the integral formula, which are summed over the entire structure. The sensitivity formula can be evaluated with generality and ease even by employing piecewise linear design velocity field despite the bending model is fourth order differential equation. For the numerical implementation, commercial software ANSYS is used as analysis tool for the primal and adjoint analysis. Once the design variable set is defined using ANSYS language, shape and orientation variation vector at each node is generated by making finite difference to the shape with respect to each design parameter, and is used for the computation of sensitivity formula. Several numerical examples are taken to show the advantage of the method, in which the accuracy of the sensitivity is evaluated. The results are found excellent even by employing a simple linear function for the design velocity evaluation. Shape optimization is carried out for the geometric design of an archgrid and tilted bridge, which is to minimize maximum stress over the structure while maintaining constant weight. In conclusion, the proposed formulation is a useful and easy tool in finding optimum shape in a variety of the spatial frame structures.

• Journal of Korea Water Resources Association
• /
• v.41 no.8
• /
• pp.761-772
• /
• 2008

The object of this study is to develop the model that solves the numerically difficult problems in hydraulic engineering and to demonstrate the applicability of this model by means of various test examples, such as, verification in the gradually varied unsteady condition, three steady flow problems with the change of bottom slope with exact solution, and frictional bed with analytical solution. The governing equation of this model is the integral form of the Saint-Venant equation satisfying the conservation laws, and finite volume method with the Riemann solver is used. The evaluation of the mass and momentum flux with the HLL Riemann approximate solver is executed. MUSCL-Hancock scheme is used to achieve the second order accuracy in space and time. This study introduce the new and simple technique to discretize the source terms of gravity and hydrostatic pressure force due to longitudinal width variation for the balance of quantity between nonlinear flux and source terms. The results show that the developed model's implementation is accurate, robust and highly stable in various flow conditions with source terms, and this model is reliable for one-dimensional applications in hydraulic engineering.

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.24 no.2
• /
• pp.83-93
• /
• 1988

Recent Experimental results show that the J integral can be effectively used to obtain a valid parameter for predicting plane strain and plane stress fracture. However, only a few research results have been reported for the effect of thickness where the plane strain state can not be assumed. A purpose of this study is to find the behavior of fracture touhness and tearing modulus varing the specimen thickness. The type of specimen in the present study is compact tension (CT). The thicknesses of the low carbon steel specimens that are used in the experiments are 5, 10, 15, 20 and 25mm. The measurement of crack length is taken by optical measurements method. From the study, the followings are found; 1) The fracture toughness and the tearing modulus which are obtained by using Yoon's and Simpson's formula show more conservative than that by using Rice's and Merkel's. 2) The fracture toughness is increase in specimen thickness which is reached 15mm. Beyond this thickness the fracture toughness is decreased in specimen thickness. 3) In the case of CT specimen with the thickness ranging from 5 to 25mm, the tearing modulus which is applied the same J integral equation is almost constant. 4) By using Yoon's formula, the correlation of the plane slress fracture toughness J sub(C) with specimen thickness B is expressed as the following formula. J sub(C)/J sub(IC)=1.7-15.1(B/W)+112.9(B/W) super(2) -301.3(B/W) super(3) +260.6(B/W) super(4)

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.17 no.7 s.110
• /
• pp.623-628
• /
• 2006

This paper presents the characteristics of an antenna factor of two kinds of EMI dipole antennas with a coaxial cable balun used in the frequency range between 30 and 1,000 MHz. The integral equation for unknown current distribution is solved by the Galerkin's method of moments with piecewise sinusoidal functions. An antenna factor for EMI dipole antennas with the coaxial cable balun is derived by using the power loss concepts. We can realize two kinds of EMI dipole antennas with appropriate antenna factors in the frequency range from 30 to 1,000 MHz: 150-cm dipole length($30{\sim}300 MHz$) and 30cm dipole length($300{\sim}1,000 MHz$). To check th ε validity of the theoretical analysis, the complex antenna factor was measured using by reference antenna methods. It is shown that the calculated complex antenna factor is good agreement with experimental results.