• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.31-38
• /
• 2016

The FDS-HCIB method is expanded to simulate water-oil-air flows around oil booms under relative motion, which is intended to increase the thickness of contained oil. The FDS scheme captures discontinuity in the density field and abrupt change of the tangential velocity across an interface without smearing. The HCIB method handles relative motions of thin oil booms with ease. To validate the developed FDS-HCIB code for water-oil-air flow around a moving body, the computed results are compared with the reported experimental results on the shape, length, and thickness of the oil slicks under towing. It is observed that the increase in pressure field between two barriers lifts the oil slick and the interfacial wave propagates and reflects as one barrier gets closer to the other barrier.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.3
• /
• pp.27-36
• /
• 2002

Numerical stability is studied based on numerics and mathematics. It is frequently observed in the region where velocity is zero. In that region, the Euler equation have numerous solutions and, thus, it is impossible to determine a unique solution with only governing equations. However, a unique solution can be determined by additional outer flow conditions or outer numerical discontinuity calculation since the information or a unique solution under undisturbed conditions is lost by disturbances. In this reason, the numerical scheme comsistent with Euler equations cannot remove shock instability completely.

• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.8-15
• /
• 2012

A code developed using the flux-difference splitting scheme on the hybrid Cartesian/immersed boundary method is applied to simulate three-dimensional internal waves. The material interface is regarded as a moving contact discontinuity and is captured on the basis of mass conservation without any additional treatment across the interface. Inviscid fluxes are estimated using the flux-difference splitting scheme for incompressible fluids of different density. The hybrid Cartesian/immersed boundary method is used to enforce the boundary condition for a moving three-dimensional body. Immersed boundary nodes are identified within an instantaneous fluid domain on the basis of edges crossing a boundary. The dependent variables are reconstructed at the immersed boundary nodes along local normal lines to provide the boundary condition for a discretized flow problem. The internal waves are simulated, which are generated by an pitching ellipsoid near an material interface. The effects of density ratio and location of the ellipsoid on internal waves are compared.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.33 no.2
• /
• pp.187-192
• /
• 2013

Fatigue crack was detected from a temperature change around surface crack using the thermographic technique. Thermal gradient across the crack decreased very much due to thermal resistance of contact surface in the crack. Heat diffusion flow passing through the discontinuity was visualized in temperature by infrared camera to find and locate the crack. A fatigue crack specimen(SM-45C), which was prepared according to KS specification and notched in its center to initiate fatigue crack from the notch tip, was heated by halogen lamp at the end of one side to generate a heat diffusion flow in lateral direction. A abrupt jump in temperature across the fatigue crack was observed in thermographic image, by which the crack could be located and sized from temperature distribution.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.9
• /
• pp.12-26
• /
• 2004

A new treatment of cell-interface flux in AUSM-type methods is introduced to reduce the numerical dissipation. Through analysis of TVD limiters, a criterion for the more accurate prediction of cell-interface state is found out and M-AUSMPW+ is developed by determining the transferred property newly and appropriately within the criterion. The superiority of M-AUSMPW+ is clearly revealed in multi-dimensional flow problems. It can eliminate numerical dissipation effectively in a non-flow aligned grid system. As a result, M-AUSMPW+ is shown to be much more accurate and effective than other previous schemes in multi-dimensional problems. Through a stationary contact discontinuity, a vortex flow, a shock wave/boundary layer interactions and viscous shock tube problems, it is verified that accuracy of M-AUSMPW+ is improved.

• Journal of the Korean Geotechnical Society
• /
• v.38 no.9
• /
• pp.45-52
• /
• 2022

Because discontinuity in the rock mass and contact of soil-structure interaction exhibits coupled thermal-hydromechanical (THM) behavior, it is necessary to develop an interface element based on the full governing equations. In this study, we derive force equilibrium, fluid continuity, and energy equilibrium equations for the interface element. Additionally, we present a stiffness matrix of the elastoplastic mechanical model for the interface element. The developed interface element uses six nodes for displacement and four nodes for water pressure and temperature in a two-dimensional analysis. The fully coupled THM analysis for fluid injection into a fault can model the complicated evolution of injection pressure due to decreasing effective stress in the fault and thermal contraction of the surrounding rock mass. However, the result of hydromechanical analysis ignoring thermal phenomena overestimates hydromechanical variables.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.2A
• /
• pp.207-214
• /
• 2008

This study presents a way to validate the quality level of the proposed connection type and verify the experimental test, and performs a 3D nonlinear analysis corresponding to the experimental test. Two mixed-structure beams were cast and tested under a four-point static loading. Force-displacement relation, force-strain relation, force-opening width, and failure mode were observed from comparing the numerical results of the adopted FE model. Nonlinear analysis of mixed structures was carried out by utilizing the contact elements of a general purpose structural analysis computer program (ABAQUS). The results of numerical and experimental simulation show that the proposed L-shaped connection has greater stiffness under flexural loading and better structural performance with regard to the connection.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.1 no.1 s.1
• /
• pp.139-155
• /
• 2001

We investigate the vertical velocity models beneath the newly installed broadband seismic network of KMA (Korea Meteorological Administration) by using receiver function inversion technique. The seismic phases are primarily P-to-S conversions and reverberations generated at the two highest impedance interfaces like the Moho (crust-mantle boundary) and the sediment-basement contact. We obtained the teleseismic P-wave receiver functions, which were derived from teleseismic records of Seoul (SEO), Inchon (INCN), Tejeon (TEJ) , Sosan (SOS/SES), Kangnung (KAN), Ulchin (ULC/ULJ), Taegu (TAG), Pusan (PUS), and Ullung-do (ULL) stations. For Kwangju (KWA/KWJ) and Chunchon (CHU) stations, the Moho conversion Ps arrivals and waveforms of radial receiver functions are azimuthally inconsistent and unclear. From the receiver function inversion result, we found that crustal thickness is 29 km at INCN, SEO, and SOS (SES) stations, 28 km at KAN station in the Kyonggi Massif, 32 km at TEJ station in Okchon Folded Belt, 34 km at TAG, 33 km at PUS station in the Kyongsang Basin, 32 km at KWJ station (readjusted station by prior KWA station) included in the Youngdong-Kwangju Depression Zone, 28 km at ULC station in the eastern margin of the Ryongnam Massif, and 17 km at ULL station in the Ullung Island of the East Sea, respectively. The Moho configuration of INCN, SOS, KWJ, and KAN stations show a laminated smooth transition zone with a 3-5 km thick. The upper crusts(${\sim}5km$) of KAN, ULC, and PUS stations show complex structures with a high velocity. The unusually thick crusts are found at the TAG and PUS stations in the Kyongsang Basin compared to the thin (29-32 km) crust of the western part (INCN, SEO, SOS, TEJ, and KWA stations) The crustal thickness beneath Ullung Island (ULL station) shows the suboceanic crust with about 17 km thickness and complex with a high velocity layer of the upper crust, and the amplitudes of Incoming Ps waves from the western direction are relatively large compared to those from othor directions.

• Journal of the Korean earth science society
• /
• v.23 no.4
• /
• pp.357-368
• /
• 2002

In this study we have tried to explain the barrel-shaped morphology for young supernova remnants considering the dynamical effects of the ejecta. We consider the magnetic field amplification resulting from the Rayleigh-Taylor instability near the contact discontinuity. We can generate the synthetic radio image assuming the cosmic-ray pressure and calculate the azimuthal intensity ratio (A) to enable a quantitative comparison with observations. The postshock magnetic field are amplified by shearing, stretching, and compressing at the R-T finger boundary. The evolution of the instability strongly depends on the deceleration of the ejecta and the evolutionary stage of the remnant. the strength of the magnetic field increases in the initial phase and decreases after the reverse shock passes the constant density region of the ejecta. However, some memory of the earlier phases of amplification is retained in the interior even when the outer regions turn into a blast wave. The ratio of the averaged magnetic field strength at the equator to the one at the pole in the turbulent region can amount to 7.5 at the peak. The magnetic field amplification can make the large azimuthal intensity ratio (A=15). The magnitude of the amplification is sensitive to numerical resolution. This mens the magnetic field amplification can explain the barrel-shaped morphology of young supernova remnant without the dependence of the efficiency of the cosmic-ray acceleration on the magnetic field configuration. In order for this mechanism to be effective, the surrounding magnetic field must be well-ordered. The small number of barrel-shaped remnants may indicate that this condition rarely occurs.