• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.2
• /
• pp.65-72
• /
• 2001

This paper demonstrates 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 an 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.

• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.167-169
• /
• 2002

The extended Born, or localized nonlinear approximation of integral equation (IE) solution has been applied to inverting single-hole electromagnetic (EM) data using a cylindrically symmetric model. The extended Born approximation is less accurate than a full solution but much superior to the simple Born approximation. When applied to the cylindrically symmetric model with a vertical magnetic dipole source, however, the accuracy of the extended Born approximation is greatly improved because the electric field is scalar and continuous everywhere. One of the most important steps in the inversion is the selection of a proper regularization parameter for stability. Occam's inversion (Constable et al., 1987) is an excellent method for obtaining a stable inverse solution. It is extremely slow when combined with a differential equation method because many forward simulations are needed but suitable for the extended Born solution because the Green's functions, the most time consuming part in IE methods, are repeatedly re-usable throughout the inversion. In addition, the If formulation also readily contains a sensitivity matrix, which can be revised at each iteration at little expense. The inversion algorithm developed in this study is quite stable and fast even if the optimum regularization parameter Is sought at each iteration step. Tn this paper we show inversion results using synthetic data obtained from a finite-element method and field data as well.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method 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 derived from the mathematical vector identity. 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 calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.27 no.6
• /
• pp.782-786
• /
• 1994

Integrated finite difference method (IFDM) is used to solve one dimensional free oscillation problem in the coastal area. To evaluate the solution accuracy of IFDM in free oscillation analysis, two finite difference equations based on area discretization method and point discretization method are derived from the governing equations of free oscillation, respectively. The difference equations are transformed into a generalized eigenvalue problem, respectively. A numerical example is presented, for which the analytical solution is available, for comparing IFDM to conventional finite difference equation (CFDM), qualitatively. The eigenvalue matrices are solved by sub-space iteration method. The numerical results of the two methods are in good agreement with analytical ones, however, IFDM yields better solution than CFDM in lower modes because IFDM only includes first order differential operator in finite difference equation by Green's theorem. From these results, it is concluded that IFDM is useful for the free oscillation analysis in the coastal area.

• Journal of Energy Engineering
• /
• v.8 no.4
• /
• pp.540-545
• /
• 1999

The classical Fourier heat conduction equation is invalid at temperatures near absolute zero or at very early times in highly transient heat transfer processes. In such situations, a hyperbolic equation model for heat conduction based on the modified Fourier law is introduced because the wave nature of heat propagation becomes dominant. The Fourier model and the hyperbolic model for heat conduction are analyzed by using the Green's function technique together with the integral transform. Analytical expressions for the heat flux and temperature distributions in a finite slab subjected to a periodic surface heating at one of its surfaces are presented and the results obtained from each model are compared with each other. The thermal wave implied b the hyperbolic model is shown to travel through a medium and to reflect back toward the origin at the other insulated surface. On the other hand, the heat by the Fourier model propagates at an infinite speed instantaneously after a thermal disturbance is felt throughout the medium.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.10a
• /
• pp.187-191
• /
• 2003

The motion and time-mean drift force of a 2-D floating BBDB in waves are studied with and without taking account of fluctuating air pressure in the air chamber. It has been found numerically that the drift for a of the BBDB is in the reverse direction of propagation of the incident waves over specific frequency ranges as found by McCormick through his experiment work. The drift force is calculated by Pinkster's near-field method. Since Maruo's formula method for the drift force is always positive, Maruo's formula is only approximate and should be replaced by the correct near-field method.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.5 no.4
• /
• pp.546-558
• /
• 2013

The research performed in this paper was carried out to investigate the numerical analysis of the sheet cavitation on marine propeller. The method is boundary element method (BEM). Using the Green's theorem, the velocity potential is expressed as an integral equation on the surface of the propeller by hyperboloid-shaped elements. Employing the boundary conditions, the potential is determined via solving the resulting system of equations. For the case study, a DTMB4119 propeller is analyzed with and without cavitating conditions. The pressure distribution and hydrodynamic performance curves of the propellers as well as cavity thickness obtained by numerical method are calculated and compared by the experimental results. Specifically in this article cavitation changes are investigate in both the radial and chord direction. Thus, cross flow variation has been studied in the formation and growth of sheet cavitation. According to the data obtained it can be seen that there is a better agreement and less error between the numerical results gained from the present method and Fluent results than Hong Sun method. This confirms the accurate estimation of the detachment point and the cavity change in radial direction.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.18 no.11
• /
• pp.1669-1676
• /
• 1993

In this paper, a thin-straight monopole antenna attached to a conducting box is analyzed using a moment method with the pocklington integral expression for the exact Green's function. A modeling of a conducting box is based on a Wire-Grid method for an efficient calculation procedure. A system of linear equation is obtained using the piecewise sinusoidal basis function. And a Junction basis function is enforcing to represent the physical current on the edge of the conducting box. Especially, the junction basis functions are overlapped for the current continuity preserving on the vertexes. Numerical results are obtained for the current distribution on the unified antenna/conducting box system, input impedance and radiation pattern. The results are compared with the known data.

• Journal of the Korea Safety Management & Science
• /
• v.11 no.2
• /
• pp.95-101
• /
• 2009

The present study is to estimate the effect of wave height affecting at the front face of breakwater systems due to specification of submarine trench such as distance from breakwater to dredged area and width of dredge. The wave diffraction field, which is important hydraulic factor in the ocean, is considered to be two dimensional(2D) plane and the configuration of the submarine dredge on the sea bed designated by single horizontal long-rectangular pit system according to the various specific conditions of dredged locations. The numerical simulation is performed by using Green function based on the boundary integral equation and meshed at moving boundary conditions. The results of present numerical simulations are illustrated by applying the normal incidence. It is shown that the ratios of wave height at the front face of breakwater was varied by dependance of distant from breakwater to dredged area and width of dredge. It means that, when the navigation channel or pit breakwater is dredged on seabed, engineers have to consider the specification of dredge. This study can effectively be utilized for safety assessment to various breakwater systems in the ocean field and provided for safety construction of offshore structure.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.06a
• /
• pp.1752-1761
• /
• 2000

This paper describes the far-field estimation using the near-field measurement data. Measurement in far-field region gives us the acoustical characteristics of the source but in general measurement is made in near-field such as acoustic water tank or anechoic chamber, so far-field acoustical characteristics of the source should be predicted from near-field data. In this case, the number of measurement points in the near field which relates to the accuracy of the predicted field and the amount of data processing, should be optimized. Existing papers say that measurement points is proportional to kL and depends on geometry and directivity of the source. But they do not give us any definite criterion for the required number of measurement points. Boundary Collocation Method which is one of the far-field prediction methods, is analyzed based on Helmholtz integral equation and Green function and it has been found that the number of measurement points is optimized as 0.54kL which is about one half of the existing results.