• Ocean Systems Engineering
• /
• v.4 no.2
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.

• Journal of Welding and Joining
• /
• v.29 no.2
• /
• pp.64-71
• /
• 2011

In this paper, crack propagation analyses in the inner diameter (ID) repair weld of the dissimilar metal weldment of a nozzle were performed using a finite element alternating method (FEAM). To calculate the theoretical solution for the crack tip stress intensity factor, a weak type singular integral equation consisted of crack surface traction and dislocation density function was constructed and solved in conjunction with the FEAM. A two-dimensional axisymmetric finite element nozzle model was prepared and ID repair welding was simulated. An initial crack, 10% depth of weld thickness, was assumed and crack propagation trajectory from the initial crack to the 75% depth of thickness was calculated using the FEAM. Crack growth versus time curve was also calculated and compared with the curves obtained from ASME code method. With the method constructed in this paper, crack propagation trajectory and crack growth time were calculated automatically and effectively.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.16 no.12 s.103
• /
• pp.1246-1254
• /
• 2005

Theoretical investigation is made on the electromagnetic fields generated by an impulsive point current source, fur the VV, HV, and VH problems at the interface between an isotropic upper half-space medium and a normally uniaxial lower half-space medium. The electric fields of these problems are associated only with the extraordinary-wave components in the Fourier-Laplace domain. Applying the Cagniard-de Hoop method to each problem, the time-domain solutions of the wave fields are obtained. The fields of the VV case can be expressed in explicit(integral-free) forms. The fields of the HV and VH cases are not integral-free, but the major singularities implicit in the integral solutions can be analytically extracted. The interfacial far fields in the frequency domain are determined by the singularities in the time domain.

• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.3 no.3
• /
• pp.185-190
• /
• 1974

Floor panel heating system is popular in Korea as dwelling house heating system. There are two methods for keeping floor surface warm. One method is delivering warm air under the floor such as Roman Hypocaust and Korean traditional Ondol. The other method is imbedding hot water pipes into the concrete floor slab. This paper gives basic equations for steady and transient solutions of heat conduction from hurried pipes. For steady-state solution, fin Efficiency Method and Sink and Source Method were introduced. Sink and Source Method is applied to transient state and basic solution is given in the form of Exponential Integral Function. Numerical solutions can be solved easily by digital computer from these equations.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.20 no.1
• /
• pp.18-24
• /
• 2017

This paper is a part of developing a computer code that can be used to generate synthetic IR images by calculating the outgoing infrared signal from objects. To predict the reflected component that is a part of the outgoing IR signal, such as those components reflected from the target surface by the solar and sky irradiations, it is necessary to calculate the complicated BRDF values for considering the directional surface reflection characteristics. Since the calculation of reflectance using the BRDF requires a large amount of computation time due to the hemispherical integral term, it is frequently restricted in applying for a real-time prediction of IR signal. In this research, the simplified method for calculating IR reflected component has been proposed by replacing the integral terms into two parts, a directionally uniform component and a step function representing the specular component, to reduce computation time. The proposed method is proved to result in very fast calculation of the BRDF (up to 600 times faster calculations) for most of the surfaces with minimal loss of the accuracy.

• Computational Structural Engineering
• /
• v.9 no.4
• /
• pp.189-198
• /
• 1996

The fracture behaviour of concrete gravity dam mainly due to uplift pressure acting at the crack face was studied. Triangular type and parabolic type distribution of the uplift pressure including uniform type were first considered in case of calculating stress intensity factor(SIF) by the surface integral method. Second, the directions of crack propagation according to the uplift pressure distribution were pursued by FRANC(FRacture ANalysis Code). Third, critical crack lengths according to the uplift pressure distribution under the overflow depth were calculated. The SIF values obtained from the surface integral method were compared with those by FRANC and relatively good agrements could be obtained between both of them. And it could be seen that the direction of crack propagation in case of triangular pressure distribution was a little benter to the dam base than the one by the uniform type. Maximum critical crack lengths under the overflow depth were obtained at about 2/5-1/2 of the dam height.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.36 no.2
• /
• pp.70-79
• /
• 2024

In this study, we introduce a linear spectral method capable of simulating wave generation and transformation caused by a moving bottom in a 3-dimensional space. The governing equations are linear dynamic free-surface boundary conditions and linear kinematic free-surface boundary conditions, which are solved in Fourier space. Solved velocity potential and free-surface displacement should satisfy continuity equation and kinematic bottom boundary condition. For numerical analysis, a 4th order Runge-Kutta method was utilized to analyze the time integral. The results obtained in Fourier space can be converted into velocity potential and free-surface displacement in a real space using inverse Fourier transform. Regular waves generated by various types of moving bottoms were simulated with the linear spectral method. Additionally, obliquely generated regular waves using specified bottom movements were simulated. The results obtained from the spectral method were compared to analytical solutions, showing good agreement between the two.

• Journal of Mechanical Science and Technology
• /
• v.18 no.8
• /
• pp.1435-1450
• /
• 2004

Experimental data are presented which describe the effects of a combustor-level high free-stream turbulence on the near-wall flow structure and heat/mass transfer on the endwall of a linear high-turning turbine rotor cascade. The end wall flow structure is visualized by employing the partial- and total-coverage oil-film technique, and heat/mass transfer rate is measured by the naphthalene sublimation method. A turbulence generator is designed to provide a highly-turbulent flow which has free-stream turbulence intensity and integral length scale of 14.7% and 80mm, respectively, at the cascade entrance. The surface flow visualizations show that the high free-stream turbulence has little effect on the attachment line, but alters the separation line noticeably. Under high free-stream turbulence, the incoming near-wall flow upstream of the adjacent separation lines collides more obliquely with the suction surface. A weaker lift-up force arising from this more oblique collision results in the narrower suction-side corner vortex area in the high turbulence case. The high free-stream turbulence enhances the heat/mass transfer in the central area of the turbine passage, but only a slight augmentation is found in the end wall regions adjacent to the leading and trailing edges. Therefore, the high free-stream turbulence makes the end wall heat load more uniform. It is also observed that the heat/mass transfers along the locus of the pressure-side leg of the leading-edge horseshoe vortex and along the suction-side corner are influenced most strongly by the high free-stream turbulence. In this study, the end wall surface is classified into seven different regions based on the local heat/mass transfer distribution, and the effects of the high free-stream turbulence on the local heat/mass transfer in each region are discussed in detail.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1478-1485
• /
• 1993

A shape optimization problem is formulated to determine the optimal position of heating lines in a compression molding die. The objective of the problem is that the cavity surface would be maintained by a prescribed uniform temperature. A boundary integral equation for the sensitivity of the temperature in terms of hole position is derived using the method of shape design sensitivity analysis. The boundary element method is employed to analyze the temperature and sensitivity field of the die. The sensitivity calculation algorithm is incorporated in an optimization routine. To demonstrate a numerical implementation, an example problem arising in thermal design of a compression molding die is dealt with, showing that the number of heating lines chosen for the die strongly affects the ultimate uniformity of the cavity surface temperature.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.2
• /
• pp.1-13
• /
• 1986

The nonlinear hydrodynamic forces acting on a two-dimensional circular cylinder, oscillating with large amplitude in the free surface, are calculated by using the Semi-Lagrangian Time-Step-ping Method used by O.M. Faltinsen. In present calculation the position and the potential value of free surface are calculated using the exact kinematic and dynamic free surface boundary condition. At each time step an integral equation is solved to obtain the value of potential and normal velocity along the boundaries, consisting of both the body surface and the free surface. Some effort was devoted to the elimination of instability arising in the range of high frequency. Numerical simulations were performed up to the 3rd or 4th period which seems to be enough for the transient effect to die out. Each harmonic component and time-mean force are obtained by the Fourier transform of forces in time domain. The results are compared with others' experimental and theoretical results. Particularly, the calculation shows the tendency that the acceleration-phase 1st-harmonic component(added mass) increases as the motion amplitude increases and a reverse tendency in the velocity-phase 1st-harmonic component(damping coefficient). The Yamashita's experimental result also shows the same tendency. In general, the present result show relatively good agreement with the Yamashita's experimental result except for the time-mean force.