• Journal of the Korean Institute of Gas
• /
• 제7권4호
• /
• pp.44-52
• /
• 2003

In this work dynamic heat transfer in a CPFS (cable penetration fire stop) system built in the firewall of nuclear power plants is three-dimensionally investigated to develop a test-simulator that can be used to verify effectiveness of the sealant. Dynamic heat transfer in the fire stop system is formulated in a parabolic PDE (partial differential equation) subjected to a set of initial and boundary conditions. First, the PDE model is divided into two parts; one corresponding to heat transfer in the axial direction and the other corresponding to heat transfer on the vertical planes. The first PDE is converted to a series of ODEs (ordinary differential equations) at finite discrete axial points for applying the numerical method of SOR (successive over-relaxation) to the problem. The ODEs are solved by using an ODE solver In such manner, the axial heat flux can be calculated at least at the finite discrete points. After that, all the planes are separated into finite elements, where the time and spatial functions are assumed to be of orthogonal collocation state at each element. The initial condition of each finite element can be obtained from the above solution. The heat fluxes on the vertical planes are calculated by the Galerkin FEM (finite element method). The CPFS system was modeled, simulated, and analyzed here. The simulation results were illustrated in three-dimensional graphics. Through simulation, it was shown clearly that the temperature distribution was influenced very much by the number, position, and temperature of the cable stream, and that dynamic heat transfer through the cable stream was one of the most dominant factors, and that the feature of heat conduction could be understood as an unsteady-state process.

• Journal of Korean Society of Transportation
• /
• 제23권2호
• /
• pp.161-172
• /
• 2005

This study is to develop Bus Transit Route Analysis and Evaluation Model that can product the quantitative performance measures for Bus Transit Route Network Design. So far, in Korea, there are no so many models that evaluate a variety of other performance measures or service quality that are of concern to both the transit users and operator because of lower-level bus database system and transit route network analysis algorithm's limit. The BTRAEM in this research differ from the previous approach in that the BTRAEM employs a multiple path transit trip assignment model that explicitly considers the transfer and different travel time after boarding. And we develop input-output data structure and quantitative performance measure for the BTRAEM. In the numerical experimental applying BTRAEM to Mandl transit network, We got the meaningful results on performance measure of bus transit route network. In the future, we expect BTRAEM to give a good solution in real transit network.

• Proceedings of the Korea Water Resources Association Conference
• /
• 한국수자원학회 2019년도 학술발표회
• /
• pp.258-258
• /
• 2019

계산 격자에 기반하여 천수 흐름을 모의할 때, 그 격자에 담긴 물의 양을 정확하게 파악할 필요가 있다. 예를 들어, 초기조건으로 수위가 부여된다면 계산격자의 기하 특성에 맞추어 흐름 변수인 수심이나 흐름 단면적으로 바꾸어야하기 때문이다. 필요에 따라서는 모의 결과를 수위로 보이거나 격자 속 수심을 계산에 사용할 수도 있으므로 그 역변환도 고려되어야 한다. 2차원의 삼각형 계산격자에 대해서는 물의 부피와 수위 관계(volume/free-surface relationship)가 이미 정확(exact)하게 구명되어 있다(Hwang, 2017, J. KWRA). 그런데 1차원 문제의 횡단면에서 흐름 단면적과 수위의 관계(area/free-surface relationship)는 수위로부터 면적 환산에 대해서는 정확하나 그 역변환은 그렇지 않다. 매 시간 단계에서 갱신된 흐름 단면적으로부터 수위를 환산하기 위해 미리 작성된 면적-수위 자료를 이용한 선형 보간이 적용된다(Goodell, 2011, The RAS Solution). 이때, 환산 정확도는 자료의 해상도에 의존된다. 다행히 하천 횡단면 대부분을 채워 흐르는 홍수모의에서는 이 문제가 그리 심각하지는 않다. 심지어 수위가 복단면 저수로 턱에 걸쳐있어 흐름단면적이 급변하는 경우에도 환산 수위의 정확도는 크게 훼손되지 않는다. 그러나 미미한 환산 오차일지라도 그로 인해 수위가 저수로 턱을 넘거나 그보다 작을 수 있다. 이 경우, 홍수터의 잠김여부에 따라 수면폭(top width)이 실제 계산 결과에 비해 크게 달라질 수밖에 없다. 수면폭 오차는 그것을 이용하여 결정되는 수리 수심(hydraulic depth)이나 평균 하상고(mean bed level)의 산정에도 전파된다. 이 연구에서는 하천 횡단면에서 수위와 흐름 단면적 사이의 환산 정확도를 크게 높일 수 있는 기법을 제시하였다. 먼저 하천 횡단면에서 주어진 수위에 대해 흐름 단면적을 산정할 수 있는 알고리듬을 보였다. 또한, 횡단면에서 수위와 흐름 단면적의 관계가 단조 증가 함수(monotonically increasing function)임에 착안하여 그 역변환에 대해 해 찾기(root finding) 방법의 하나인 Brent 기법을 적용하였다. 이 기법은 주어진 구간에서 도함수가 알려져 있지 않은 경우에 대해서도 효과적으로 해를 찾을 수 있는 것으로 알려져 있다(Press et al., 2002, Numerical Recipes in C, 2nd Ed.). 내성천 하류 수계의 333개 단면에서 수면폭에 대한 상대 오차를 살펴보면, 선형 보간에 의한 기존 방법으로는 면적-수위 자료의 수가 1,000개가 되어도 그 최대치가 1% 이내에 들지 않은 반면, 이 연구에서 제시한 기법으로 면적-수위 자료 없이도 1% 이내로 줄어드는 것을 확인하였다. 다만, 반복 계산에 의한 계산 시간의 증대를 피할 수 없다. 미리 작성된 면적-수위 자료를 이용하면 계산 비용을 줄일 수 있으며, 약 35개의 구간으로 나누었을 때 비용 대비오차가 적절하였다. 이 연구는 한국건설기술연구원(주요사업 과제번호: 20190116-001)의 지원에 의한 것이다.

• Journal of Korean Tunnelling and Underground Space Association
• /
• 제24권6호
• /
• pp.745-756
• /
• 2022

The road tunnel is a semi-closed space that is blocked on all sides except the entrance and exit, and in the event of a fire, the smoke of the fire spreads longitudinally due to heat buoyancy caused by the fire and air currents that always exist in the tunnel. To solve this problem, smoke removal facilities are installed in road tunnels to secure a safe evacuation environment by controlling the direction of movement of smoke or directly smoking at fire points. In urban areas, the service level of urban roads decreases due to the increase in traffic due to the increase in population, and as a solution, the construction of underground roads in urban areas is increasing. When a fire occurs during hydrogen leakage through TPRD of a hydrogen fuel cell vehicle (FCEV), the fire intensity depends on the amount of leakage, and the maximum fire intensity depends on the orifice diameter of the TPRD. Considering the TPRD orifice diameter of 1.8 mm, this study analyzed the diffusion distance of fire smoke according to the wind speed of the roadway and the opening interval of the large exhaust port when the maximum fire intensity was 15 MW. As a result, it was analyzed that air flow in the tunnel could be controlled if the wind speed of the road in the tunnel was less than 1.25 m/s, and smoke could be controlled within 200 m from the fire if the damper interval was 50 m and 100 m.

• Geomechanics and Engineering
• /
• 제17권4호
• /
• pp.343-354
• /
• 2019

The present study demonstrates the application of seismic petrophysics and amplitude versus angle (AVA) forward modeling to identify the reservoir fluids, discriminate their saturation levels and natural gas composition. Two case studies of the Lumshiwal Formation (mainly sandstone) of the Lower Cretaceous age have been studied from the Kohat Sub-basin and the Middle Indus Basin of Pakistan. The conventional angle-dependent reflection amplitudes such as P converted P ($R_{PP}$) and S ($R_{PS}$), S converted S ($R_{SS}$) and P ($R_{SP}$) and newly developed AVA attributes (${\Delta}R_{PP}$, ${\Delta}R_{PS}$, ${\Delta}R_{SS}$ and ${\Delta}R_{SP}$) are analyzed at different gas saturation levels in the reservoir rock. These attributes are generated by taking the differences between the water wet reflection coefficient and the reflection coefficient at unknown gas saturation. Intercept (A) and gradient (B) attributes are also computed and cross-plotted at different gas compositions and gas/water scenarios to define the AVO class of reservoir sands. The numerical simulation reveals that ${\Delta}R_{PP}$, ${\Delta}R_{PS}$, ${\Delta}R_{SS}$ and ${\Delta}R_{SP}$ are good indicators and able to distinguish low and high gas saturation with a high level of confidence as compared to conventional reflection amplitudes such as P-P, P-S, S-S and S-P. In A-B cross-plots, the gas lines move towards the fluid (wet) lines as the proportion of heavier gases increase in the Lumshiwal Sands. Because of the upper contacts with different sedimentary rocks (Shale/Limestone) in both wells, the same reservoir sand exhibits different response similar to AVO classes like class I and class IV. This study will help to analyze gas sands by using amplitude based attributes as direct gas indicators in further gas drilling wells in clastic successions.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• 제32권2호
• /
• pp.103-108
• /
• 2019

Non-equilibrium molecular dynamics simulation on the thermal boundary resistance(TBR) of an aluminum(Al)/silicon(Si) interface was performed in the present study. The constant heat flux across the Si/Al interface was simulated by adding the kinetic energy in hot Si region and removing the same amount of the energy from the cold Al region. The TBR estimated from the sharp temperature drop at the interface was independent of heat flux and equal to $5.13{\pm}0.17K{\cdot}m^2/GW$ at 300K. The simulation result was experimentally confirmed by the time-domain thermoreflectance technique. A 90nm thick Al film was deposited on a Si(100) wafer using an e-beam evaporator and the TBR on the film/substrate interface was measured using the time-domain thermoreflectance technique based on a femtosecond laser system. A numerical solution of the transient heat conduction equation was obtained using the finite difference method to estimate the TBR value. Experimental results were compared to the prediction and discussions on the nanoscale thermal transport phenomena were made.

• Journal of the Society of Disaster Information
• /
• 제17권1호
• /
• pp.120-127
• /
• 2021

Purpose: In this study, the effect of changes in the variables related to water droplets on the spray density on the floor in the analysis of the water mist fire protection system using FDS was analyzed. Method: When the spray of the water mist nozzle was analyzed in FDS, Particles Per Seconds, Particle Velocity, Size Distribution, and Spray Pattern Shape that can be set in relation to droplets were input to review the analyzed results. Result: In the analysis results, when the number of particles per second was set above a certain value, the spray density of the floor was similar. In the case of Particle Velocity, as the velocity decreases, the spray density of the central portion increases but decreases at a distance of 0.15m or more. From the analysis of the change in the size distribution function, it was found that an increase in the 𝛾 value increases the spray density of the central part, but the value at a remote location decreases. Compared to the result of applying the Gaussian distribution, the median value decreases dramatically when the uniform distribution is applied, but the value at the adjacent position increases. Conclusion: Variables related to the droplet properties of the FDS affect the spray density of the floor. Therefore, in order to increase the reliability of results before performing analyses such as fire suppression or cooling, a sufficient review of input variables is required.

• Structural Engineering and Mechanics
• /
• 제76권6호
• /
• pp.687-708
• /
• 2020

Openings in steel plate shear walls (SPSWs) are usually used for decorative designs, crossing locations of multiple utilities and/or structural objectives. However, earlier studies showed that generating an opening in an SPSW has a negative effect on the cyclic performance of the SPSW. Therefore, this study proposes tripling or doubling the steel-sheet-plate (SSP) layer and stiffening the opening of the SPSW to provide a solution to undesirable opening effects, improve the SPSW performance and provide the infill option of potential strengthening measures after the construction stage. The study aims to investigate the impact of SSP doubling with a stiffened opening on the cyclic behaviour, expand the essential data required by structural designers and quantify the SPSW performance factors. Validated numerical models were adopted to identify the influence of the chosen parameters on the cyclic capacity, energy dissipation, ductility, seismic performance factors (SPF) and stiffness of the suggested method. A finite Element (FE) analysis was performed via Abaqus/CAE software on half-scale single-story models of SPSWs exposed to cyclic loading. The key parameters included the number of SSP layers, the opening size ratios corresponding to the net width of the SSP, and the opening shape. The findings showed that the proposed assembly method found a negligible influence in the shear capacity with opening sizes of 10, 15, 20%. However, a deterioration in the wall strength was observed for openings with sizes of 25% and 30%. The circular opening is preferable compared with the square opening. Moreover, for all the models, the average value of the obtained ductility did not show substantial changes and the ultimate shear resistance was achieved after reaching a drift ratio of 4.36%. Additionally, the equivalent sectional area of the SSP in the twin and triple configuration of the SPSWs demonstrated approximately similar results. Compared with the single SSP layer, the proposed configuration of the twin SSP layer with a stiffened opening suggest to more sufficiency create SSP openings in the SPSW compared to that of other configurations. Finally, a tabular SPF quantification is exhibited for SPSWs with openings.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• 제22권1호
• /
• pp.9-17
• /
• 2021

Curved beams are increasingly used in buildings, vehicles, ships, and aircraft, which has resulted in considerable effort being directed toward developing an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic circular arches has been the subject of a large number of investigations. One of the efficient procedures for the solution of ordinary differential equations or partial differential equations is the differential quadrature method DQM. This method has been applied to a large number of cases to overcome the difficulties of the complex computer algorithms, as well as excessive use of storage due to conditions of non-linear geometries, loadings, or material properties. This study uses DQM to analyze the in-plane vibration of the circular arches considering the effects of midsurface extension and rotatory inertia. Fundamental frequency parameters are calculated for the member with various parameter ratios, boundary conditions, and opening angles. The solutions from DQM are compared with exact solutions or other numerical solutions for cases in which they are available and given to analyze the effects of midsurface extension and rotatory inertia on the frequency parameters of the circular arches.

• Journal of the Korea Society for Simulation
• /
• 제30권1호
• /
• pp.11-20
• /
• 2021

As unmanned underwater systems have recently emerged, a high-speed underwater channel modeling technique, which is one of the most important techniques in the system, has received a lot of attention. In this paper, we proposed a high-speed sound propagation model and verified the applicability through quantitative performance analyses. We used a high-order finite difference method (FDM) for wave propagation modeling in the water, and a domain decomposition method was adopted using multiple general-purpose graphics processing units (GPUs) to increase the calculation efficiency. We compared the results of the model we proposed with the analytic solution in the half-infinite media and results of the Virtual Timeseries Experiment (VirTEX) model, which is based on the ray method. Finally, we analyzed the performance of the model quantitatively using numerical examples. Through quantitative analyses of the improvement in computational performance, we confirmed that the computational speed increases linearly as the number of GPUs increases. The computation times are increased by 2 times and 8 times, respectively, when the domain size of computation and the maximum frequency are doubled. We expect that the proposed high-speed underwater channel modeling technique is able to contribute to the enhancement of national defense as an underwater communication channel model and analysis tool to develop the underwater communication technique for the unmanned underwater system.