• Tribology and Lubricants
• /
• v.38 no.3
• /
• pp.73-83
• /
• 2022

This paper is concerned with an analysis of a surface edge crack emanated from a sharp contact edge. For a geometrical model, a square wedge is in contact with a half plane whose materials are identical, and a surface perpendicular crack initiated from the contact edge exists in the half plane. To analyze this crack problem, it is necessary to evaluate the stress field on the crack line which are induced by the contact tractions and pseudo-dislocations that simulate the crack, using the Bueckner principle. In this Part I, the stress filed in the half plane due to the contact is re-summarized using an asymptotic analysis method, which has been published before by the author. Further focus is given to the stress field in the half plane due to a pseudo-edge dislocation, which will provide a stress solution due to a crack (i.e. a continuous distribution of edge dislocations) later, using the Burgers vector. Essential result of the present work is the corrective functions which modify the stress field of an infinite domain to apply for the present one which has free surfaces, and thus the infiniteness is no longer preserved. Numerical methods and coordinate normalization are used, which was developed for an edge crack problem, using the Gauss-Jacobi integration formula. The convergence of the corrective functions are investigated here. Features of the corrective functions and their application to a crack problem will be given in Part II.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.27 no.2
• /
• pp.77-82
• /
• 2023

The structural analysis module is an essential part of any integrated structural system. Diverse integrated systems today require, from the analysis module, efficient real-time responses to real-time input such as earthquake signals, extreme weather-related forces, and man-made accidents. An integrated system may also be for the entire life span of a civil structure conceived during the initial conception, developed throughout various design stages, effectively used in construction, and utilized during usage and maintenance. All these integrated systems' essential part is the structural analysis module, which must be automated and computationally efficient so that responses may be almost immediate. The finite element method is often used for structural analysis, and for automation, many effective finite element meshes must be automatically generated for a given analysis. A computationally efficient finite element mesh generation scheme based on the r-h method of mesh refinement using strain deviations from the values at the Gauss points as error estimates from the previous mesh is described. Shape factors are used to sort out overly distorted elements. A standard cantilever beam analyzed by four-node plane stress elements is used as an example to show the effectiveness of the automated algorithm for a time-domain dynamic analysis. Although recent developments in computer hardware and software have made many new applications in integrated structural systems possible, structural analysis still needs to be executed efficiently in real-time. The algorithm applies to diverse integrated systems, including nonlinear analyses and general dynamic problems in earthquake engineering.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.301-308
• /
• 2013

In this paper, an accuracy analysis of parallel method based on non-overlapping domain decomposition method is carried out. In this approach, proposed by Tak et al.(2013), the decomposed subdomains do not overlap each other and the connection between adjacent subdomains is determined via simple connective finite element named interfacial element. This approach has two main advantages. The first is that a direct method such as gauss elimination is available even in a singular problem because the singular stiffness matrix from floating domain can be converted to invertible matrix by assembling the interfacial element. The second is that computational time and storage can be reduced in comparison with the traditional finite element tearing and interconnect(FETI) method. The accuracy of analysis using proposed method, on the other hand, is inclined to decrease at cross points on which more than three subdomains are interconnected. Thus, in this paper, an accuracy analysis for a novel non-overlapping domain decomposition method with a variety of subdomain numbers which are interconnected at cross point is carried out. The cause of accuracy degradation is also analyze and establishment of countermeasure is discussed.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.13 no.7
• /
• pp.667-678
• /
• 2002

In this paper, a new radiating structure with a multi-layered two-dimensional metallic disk array was proposed for shaping the flat-topped element pattern. It is an infinite periodic planar array structure with metallic disks finitely stacked above the radiating circular waveguide apertures. The theoretical analysis was in detail performed using rigid full-wave analysis, and was based on modal representations for the fields in the partial regions of the array structure and for the currents on the metallic disks. The final system of linear algebraic equations was derived using the orthogonal property of vector wave functions, mode-matching method, boundary conditions and Galerkin's method, and also their unknown modal coefficients needed for calculation of the array characteristics were determined by Gauss elimination method. The application of the algorithm was demonstrated in an array design for shaping the flat-topped element patterns of $\pm$20$^{\circ}$ beam width in Ka-band. The optimal design parameters normalized by a wavelength for general applications are presented, which are obtained through optimization process on the basis of simulation and design experience. A Ka-band experimental breadboard with symmetric nineteen elements was fabricated to compare simulation results with experimental results. The metallic disks array structure stacked above the radiating circular waveguide apertures was realized using ion-beam deposition method on thin polymer films. It was shown that the calculated and measured element patterns of the breadboard were in very close agreement within the beam scanning range. The result analysis for side lobe and grating lobe was done, and also a blindness phenomenon was discussed, which may cause by multi-layered metallic disk structure at the broadside. Input VSWR of the breadboard was less than 1.14, and its gains measured at 29.0 GHz. 29.5 GHz and 30 GHz were 10.2 dB, 10.0 dB and 10.7 dB, respectively. The experimental and simulation results showed that the proposed multi-layered metallic disk array structure could shape the efficient flat-topped element pattern.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.4 s.74
• /
• pp.429-440
• /
• 2006

This paper comprises two specific objectives. The first is to examine the applicability of ordinary kriging interpolation(OK) to the p-adaptivity of the finite element method that is based on variogram modeling. The second objective Is to present the adaptive procedure by the hierarchical p-refinement in conjunction with a posteriori error estimator using the modified S.P.R. (superconvergent patch recovery) method. The ordinary kriging method that is one of weighted interpolation techniques is applied to obtain the estimated exact solution from the stress data at the Gauss points. The weight factor is determined by experimental and theoretical variograms for interpolation of stress data apart from the conventional interpolation methods that use an equal weight factor. In the p-refinement, the analytical domain has to be refined automatically to obtain an acceptable level of accuracy by increasing the p-level non-uniformly or selectively. To verify the performance of the modified S.P.R. method, the new error estimator based on limit value has been proposed. The validity of the proposed approach has been tested with the help of some benchmark problems of linear elastic fracture mechanics such as a centrally cracked panel, a single edged crack, and a double edged crack.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.146-153
• /
• 2004

The Numerical study has been carried out to investigate the effects of chemical reaction and thermal radiation on the rocket plume flow-field at various altitudes. The theoretical formulation is based on the Navier-Stokes equations for compressible flows along with the infinitely fast chemistry and thermal radiation. The governing equations were solved by a finite volume fully-implicit TVD(Total Variation Diminishing) code which uses Roe's approximate Riemann solver and MUSCL(Monotone Upstream-centered Schemes for Conservation Laws) scheme. LU-SGS (Lower Upper Symmetric Gauss Seidel) method is used for the implicit solution strategy. An equilibrium chemistry module for hydrocarbon mixture with detailed thermo-chemical properties and a thermal radiation module for optically thin media were incorporated with the fluid dynamics code. In this study, kerosene-fueled rocket was assumed operating at O/F ratio of 2.34 with a nozzle expansion ratio of 6.14. Flight conditions considered were Mach number zero at ground level, Mach number 1.16 at altitude 5.06km and Mach number 2.9 at altitude 17.34km. Numerical results gave the understandings on the detailed plume structures at different altitude conditions. The diffusive effect of the thermal radiation on temperature field and the effect of chemical recombination during the expansion process could be also understood. By comparing the results from frozen flow and infinitely fast chemistry assumptions, the excess temperature of the exhaust gas resulting from the chemical recombination seems to be significant and cannot be neglected in the view point of performance, thermal protection and flow physics.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.66-71
• /
• 2013

In this paper, the supersonic flows around space launch vehicles have been numerically simulated by using a 3-D RANS flow solver. The focus of the study was made for investigating plume-induced flow separation(PIFS). For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras model was employed for the closure of turbulence. The Gauss-Seidel iteration was used for time integration. To validate the flow solver, calculation was made for the 0.04 scale model of the Saturn-5 launch vehicle at the supersonic flow condition without exhaust plume, and the predicted results were compared with the experimental data. Good agreements were obtained between the present results and the experiment for the surface pressure coefficient and the Mach number distribution inside the boundary layer. Additional calculations were made for the real scale of the Saturn-5 configuration with exhaust plume. The flow characteristics were analyzed, and the PIFS distances were validated by comparing with the flight data. The KSLV-1 is also simulated at the several altitude conditions. In case of the KSLV-1, PIFS was not observed at all conditions, and it is expected that PIFS is affected by the nozzle position.

• Proceedings of the KIEE Conference
• /
• 1996.07c
• /
• pp.1446-1449
• /
• 1996

Thin films of diamond-like carbon(DLC) have been deposited using a magnetron plasma-enhanced chemical vapor deposition(PECVD) method with an rf(13.56 MHz) plasma of $C_{3}H_{8}$. From the Langmuir probe I-V characteristics, it can be observed that increasing the magnetic field yields an increase of the temperature($T_e$) and density($N_e$) of electron. At a magnetic field of 82 Gauss, the estimated values of $T_e$ and $N_e$ are approximately $1.5\;{\times}\;10^5$ K(13.5 eV) and $1.3\;{\times}\;10^{11}\;cm^{-3}$, respectively. Such a highly dense plasma can be attributed to the enhanced ionization caused by the cyclotron motion of electrons in the presence of a magnetic field. On the other hand, the negative dc self-bias voltage($-V_{sb}$) decreases with an increasing magnetic field, which is irrespective of gas pressure in the range of $1{\sim}7$ mTorr. This result is well explained by a theoretical model considering the variation of $T_e$. Deposition rates of DLC films increases with a magnetic field. This may be due to the increased mean free path of electrons in the magnetron plasma. Structures of DLC films are examined by using various techniques such as FTIR and Raman spectroscopy. Most of hydrocarbon bonds in DLC films prepared consist of $sp^3$ tetrahedral bonds. Increasing the rf power leads to an enhancement of cross-linking of carbon atoms in DLC films. At approximately 140 W, the maximum film density obtained is about 2.4 $g/cm^3$.

• Journal of Korean Society of Industrial and Systems Engineering
• /
• v.40 no.4
• /
• pp.137-146
• /
• 2017

For the competitive business environment under purchase dependence, this paper proposes a new approximate calculation of order fill rate which is a probability of satisfying a customer order immediately using the existing inventory. Purchase dependence is different to demand dependence. Purchase dependence treats the purchase behavior of customers, while demand dependence considers demand correlation between items, between regions, or over time. Purchase dependence can be observed in such areas as marketing, manufacturing systems, and distribution systems. Traditional computational methods have a difficulty of the curse of dimensionality for the large cases, when deriving the stationary joint distribution which is utilized to calculate the order fill rate. In order to escape the curse of dimensionality and protect the solution from diverging for the large cases, we develop a greedy iterative search algorithm based on the Gauss-Seidel method. We show that the greedy iterative search algorithm is a dependable algorithm to derive the stationary joint distribution of on-hand inventories in the retailer system by conducting a comparison analysis of a greedy iterative search algorithm with the simulation. In addition, we present some managerial insights such as : (1) The upper bound of order fill rate can be calculated by the one-item pure system, while the lower bound can be provided by the pure system that consists of all items; (2) As the degree of purchase dependence declines while other conditions remain same, it is observed that the difference between the lower and upper bounds reduces, the order fill rate increases, and the order fill rate gets closer to the upper bound.

• Investigative Magnetic Resonance Imaging
• /
• v.20 no.4
• /
• pp.215-223
• /
• 2016

Purpose: The management of metal-induced field inhomogeneities is one of the major concerns of distortion-free magnetic resonance images near metallic implants. The recently proposed method called "Slice Encoding for Metal Artifact Correction (SEMAC)" is an effective spin echo pulse sequence of magnetic resonance imaging (MRI) near metallic implants. However, as SEMAC uses the noisy resolved data elements, SEMAC images can have a major problem for improving the signal-to-noise ratio (SNR) without compromising the correction of metal artifacts. To address that issue, this paper presents a novel reconstruction technique for providing an improvement of the SNR in SEMAC images without sacrificing the correction of metal artifacts. Materials and Methods: Low-rank approximation in each coil image is first performed to suppress the noise in the slice direction, because the signal is highly correlated between SEMAC-encoded slices. Secondly, SEMAC images are reconstructed by the best linear unbiased estimator (BLUE), also known as Gauss-Markov or weighted least squares. Noise levels and correlation in the receiver channels are considered for the sake of SNR optimization. To this end, since distorted excitation profiles are sparse, $l_1$ minimization performs well in recovering the sparse distorted excitation profiles and the sparse modeling of our approach offers excellent correction of metal-induced distortions. Results: Three images reconstructed using SEMAC, SEMAC with the conventional two-step noise reduction, and the proposed image denoising for metal MRI exploiting sparsity and low rank approximation algorithm were compared. The proposed algorithm outperformed two methods and produced 119% SNR better than SEMAC and 89% SNR better than SEMAC with the conventional two-step noise reduction. Conclusion: We successfully demonstrated that the proposed, novel algorithm for SEMAC, if compared with conventional de-noising methods, substantially improves SNR and reduces artifacts.