• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.8
• /
• pp.757-762
• /
• 2011

We investigate the transient temperature response in biological tissue whose surface is exposed to alternately varying sinusoidal oscillation. Based on the Pennes bio-heat equation, we apply numerical analysis using a finite element method to find the effects of the physical properties of the skin layers. Three layers of tissue-epidermis, dermis, and subcutaneous-are considered as the solution region. We investigate the effects of different properties of the skin layers on the temperature profile. We also investigate the effects of the perfusion rate for the dermis, which is the most sensitive layer. The results show that the temperature profile of tissue depth has a discontinuous point when different physical properties are used.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.14 no.1
• /
• pp.76-85
• /
• 2002

The floating breakwater generally has an excellent water exchanging capability, but with some lowering harbor tranquility due to the wave transmission underneath floating body. In the initial stage of design, it is thus required to investigate several alternatives of breakwater layout for their performance of harbor tranquility. The present study aims to formulate a sort of reflection-transmission boundary condition of floating breakwater so that the existing numerical method using time dependent mild slope equation can still be applied to the case of floating breakwaters. The two and three dimensional tests were each performed to demonstrate the performance of the boundary condition. It was found that the reflection and transmission characteristics around the breakwater were well reproduced by the boundary condition. Finally, the reflection-transmission boundary condition were applied to a floating breakwater installed in an imaginary harbor with an irregular shape and bottom topography. The results surely showed that the present numerical method can effectively used in practical works related to the real sea construction of floating breakwaters.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.35 no.9
• /
• pp.983-989
• /
• 2011

A numerical algorithm for predicting fretting wear was developed using the boundary element method (BEM). A contact analysis was performed numerically using the relation between the elastic displacement and uniformly distributed loading of a rectangular patch on a semi-infinite solid. Geometrical updating based on nodal wear depths was performed. The wear depths were computed using the Archard's equation for sliding wear. In order to investigate the efficiency of BEM for predicting fretting wear, a problem involving a two-dimensional cylinder on a flat contact was analyzed, comparing it with the simulation model proposed by McColl et al. that was based on the finite element method. The developed method was then applied to the analysis of a spherical contact and it was shown that the developed simulation technique could efficiently predict fretting wear. Moreover, the effect of a step cycle on the solution obtained by the developed method was investigated.

• Transactions of the Korean hydrogen and new energy society
• /
• v.32 no.5
• /
• pp.365-373
• /
• 2021

In this study, evaporation of LN2 non-spreading pool on concrete plate was dealt with experimentally. The thermophysical properties of concrete, which is a composite material, were obtained by minimizing the difference between the numerical analysis results obtained from the assumed properties and the results from experiments. The thermal energy required for evaporation of the liquid pool is supplied from the concrete plate and the wall of the container. As a result of the measurement, the thermal energy flowing in from the wall was negligible compared to the one supplied from the concrete plate. It was found that the measured evaporation rate of the liquid pool by the heat energy supplied through the concrete plate agrees well with the PTC model except for the initial section of the experiment. The validity of the semi-infinite assumption and the one-dimensional assumption, which are the main conditions of the PTC model, was also verified through experiments. The evaporation rate model in the non-spreading pool discussed in this study can provide a basic frame for the one in the spreading pool, which is a meaningful result considering that the spreading pool is very realistic compared to the non-spreading pool.

• Tribology and Lubricants
• /
• v.38 no.3
• /
• pp.84-92
• /
• 2022

In Part I, developed was a method to obtain the stress field due to an edge dislocation that locates in an elastic half plane beneath the contact edge of an elastically similar square wedge. Essential result was the corrective functions which incorporate a traction free condition of the free surfaces. In the sequel to Part I, features of the corrective functions, F_kij,(k = x, y;i,j = x,y) are investigated in this Part II at first. It is found that F_xxx(ŷ) = F_xyx(ŷ) where ŷ = y/η and η being the location of an edge dislocation on the y axis. When compared with the corrective functions derived for the case of an edge dislocation at x = ξ, analogy is found when the indices of y and x are exchanged with each other as can be readily expected. The corrective functions are curve fitted by using the scatter data generated using a numerical technique. The algebraic form for the curve fitting is designed as F_kij(ŷ) = $\frac{1}{\hat{y}^{1-{\lambda}}I+yp}$$\sum_{q=0}^{m}{\left}$$\left[A_q\left(\frac{\hat{y}}{1+\hat{y}} \right)^q \right]$ where λ_I=0.5445, the eigenvalue of the adhesive complete contact problem introduced in Part I. To investigate the exponent of F_kij, i.e.(1 - λ_I) and p, Log|F_kij|(ŷ)-Log|(ŷ)| is plotted and investigated. All the coefficients and powers in the algebraic form of the corrective functions are obtained using Mathematica. Method of analyzing a surface perpendicular crack emanated from the complete contact edge is explained as an application of the curve-fitted corrective functions.

• Geomechanics and Engineering
• /
• v.31 no.5
• /
• pp.477-488
• /
• 2022

Pile composite foundation (PCF) has been commonly applied in practice. Existing research has focused primarily on semi-infinite media having equal pile lengths with little attention given to the effects of inclined bedrock and dissimilar pile lengths. This investigation considers the effects of inclined bedrock on vertical loaded PCF with dissimilar pile lengths. The pile-soil system is decomposed into fictitious piles and extended soil. The Fredholm integral equation about the axial force along fictitious piles is then established based on the compatibility of axial strain between fictitious piles and extended soil. Then, an iterative procedure is induced to calculate the PCF characteristics with a rigid cap. The results agree well with two field load tests of a single pile and numerical simulation case. The settlement and load transfer behaviors of dissimilar 3-pile PCFs and the effects of inclined bedrock are analyzed, which shows that the embedded depth of the inclined bedrock significantly affects the pile-soil load sharing ratios, non-dimensional vertical stiffness N₀/wdE_s, and differential settlement for different length-diameter ratios of the pile l/d and pile-soil stiffness ratio k conditions. The differential settlement and pile-soil load sharing ratios are also influenced by the inclined angle of the bedrock for different k and l/d. The developed model helps better understand the PCF characteristics over inclined bedrock under vertical loading.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.6
• /
• pp.549-558
• /
• 2012

A topology optimization method for phononic crystals is developed for the design of sound barriers, using the level set approach. Given a frequency and an incident wave to the phononic crystals, an optimal shape of periodic inclusions is found by minimizing the norm of transmittance. In a sound field including scattering bodies, an acoustic wave can be refracted on the obstacle boundaries, which enables to control acoustic performance by taking the shape of inclusions as the design variables. In this research, we consider a layered structure which is composed of inclusions arranged periodically in horizontal direction while finite inclusions are distributed in vertical direction. Due to the periodicity of inclusions, a unit cell can be considered to analyze the wave propagation together with proper boundary conditions which are imposed on the left and right edges of the unit cell using the Bloch theorem. The boundary conditions for the lower and the upper boundaries of unit cell are described by impedance matrices, which represent the transmission of waves between the layered structure and the semi-infinite external media. A level set method is employed to describe the topology and the shape of inclusions. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. Through several numerical examples, the applicability of the proposed method is demonstrated.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.14 no.3
• /
• pp.11-22
• /
• 2010

An analysis method is proposed for the transient linear or nonlinear analysis of dynamic interactions between a flexible dam body and reservoir impounding compressible water under earthquake loadings. The coupled dam-reservoir system consists of three substructures: (1) a dam body with linear or nonlinear behavior; (2) a semi-infinite fluid region with constant depth; and (3) an irregular fluid region between the dam body and far field. The dam body is modeled with linear and/or nonlinear finite elements. The far field is formulated as a displacement-based transmitting boundary in the frequency domain that can radiate energy into infinity. Then the transmitting boundary is transformed for the direct coupling in the time domain. The near field region is modeled as a compressible fluid contained between two substructures. The developed method is verified and applied to various earthquake response analyses of dam-reservoir systems. Also, the method is applied to a nonlinear analysis of a concrete gravity dam. The results show the location and severity of damage demonstrating the applicability to the seismic evaluation of existing and new dams.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.31 no.4
• /
• pp.423-434
• /
• 1995

It was used boundary element method(BEM) and analysed axisymmetric problem to solve hydraulic cylinder for large vessel acting uniform internal pressure(25N/m super(2)) within elastic limit. This paper was utilized the carbon steel tubes for machine structural purposed model, inner radius was 150mm and outer radius was 250mm, axial length was semi-infinite and the isoparametric element was used. The important results obtained in this study were summarized as follows. Radial, tangential and shearing stress occured the maximum stresses(48, -20 and 34MPa) at the inner radius and the minimum stresses(32, -4 and 18MPa) at the outer radius of the hydraulic cylinder for large vessel. But negative signs have meaning compressive stress and stress diminution ratio was about 0.15MPa/mm. The use of isoparametric element raised accuracy and the increment of input data lessened the error in internal point but computer run-time was increased. The double node was improved the internal solutions to settle discontinuity at corner and the double exponential formula lessened error of stress value at boundary neighborhood. And then coincidence between the analytical and exact results is found to be fairly good, showing that the proposed analytical by BEM is reliable.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2013.02a
• /
• pp.334-334
• /
• 2013

In the past decade, the infrared detectors based on intersubband transition in quantum dots (QDs) have attracted much attention due to lower dark currents and increased lifetimes, which are in turn due a three-dimensional confinement and a reduction of scattering, respectively. In parallel, focal plane array development for infrared imaging has proceeded from the first to third generations (linear arrays, 2D arrays for staring systems, and large format with enhanced capabilities, respectively). For a step further towards the next generation of FPAs, it is envisioned that a two-dimensional metal hole array (2D-MHA) structures will improve the FPA structure by enhancing the coupling to photodetectors via local field engineering, and will enable wavelength filtering. In regard to the improved performance at certain wavelengths, it is worth pointing out the structural difference between previous 2D-MHA integrated front-illuminated single pixel devices and back-illuminated devices. Apart from the pixel linear dimension, it is a distinct difference that there is a metal cladding (composed of a number of metals for ohmic contact and the read-out integrated circuit hybridization) in the FPA between the heavily doped gallium arsenide used as the contact layer and the ROIC; on the contrary, the front-illuminated single pixel device consists of two heavily doped contact layers separated by the QD-absorber on a semi-infinite GaAs substrate. This paper is focused on analyzing the impact of a two dimensional metal hole array structure integrated to the back-illuminated quantum dots-in-a-well (DWELL) infrared photodetectors. The metal hole array consisting of subwavelength-circular holes penetrating gold layer (2DAu-CHA) provides the enhanced responsivity of DWELL infrared photodetector at certain wavelengths. The performance of 2D-Au-CHA is investigated by calculating the absorption of active layer in the DWELL structure using a finite integration technique. Simulation results show the enhanced electric fields (thereby increasing the absorption in the active layer) resulting from a surface plasmon, a guided mode, and Fabry-Perot resonances. Simulation method accomplished in this paper provides a generalized approach to optimize the design of any type of couplers integrated to infrared photodetectors.