• Tunnel and Underground Space
• /
• v.30 no.5
• /
• pp.473-483
• /
• 2020

This paper examines the vertical stress change concentrated on mine pillar which occurs due to the stress disturbance from opening excavation at room and pillar mine by FLAC3D, a finite difference method (FDM) software. The mesh size combination is decided with a careful consideration of relative error and run-time, then its performance is verified. A series of numerical analyses is conducted and the vertical stress at central pillar was observed for the test cases of 1×1 to 11×11 mine pillars, 40 m to 320 m depth with 40 m difference. The results show that the vertical stress of pillar approaches to the similar value with the value estimated by tributary area theory(TAT) when the development area (N_P) is increased or the height of overburden (H_OB) is decreased, while it is overestimated in the opposite case. Furthermore, it also represents that the vertical stress factor (VSF) converges to a specific value when the depth is increased whille keeping the development area identical.

• Journal of Mechanical Science and Technology
• /
• v.16 no.8
• /
• pp.1102-1108
• /
• 2002

In this study, a topological design sensitivity of the ai. bearing surface (ABS) is suggested by using an adjoint variable method. The discrete form of the generalized lubrication equation based on a control volume formulation is used as a compatible condition. A residual function of the slider is considered as an equality constraint function, which represents the slider in equilibrium. The slider thickness parameters at all grid cells are chosen as design variables since they are the topological parameters determining the ABS shape. Then, a complicated adjoint variable equation is formulated to directly handle the highly nonlinear and asymmetric coefficient matrix and vector in the discrete system equation of air-lubricated slider bearings. An alternating direction implicit (ADI) scheme is utilized for the numerical calculation. This is an efficient iterative solver to solve large-scale problem in special band storage. Then, a computer program is developed and applied to a slider model of a sophisticated shape. The simulation results of design sensitivity analysis (DSA) are directly compared with those of FDM at the randomly selected grid cells to show the effectiveness of the proposed approach. The overall distribution of DSA results are reported, clearly showing the region on the ABS where special attention should be given during the manufacturing process.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2001.05a
• /
• pp.40-45
• /
• 2001

Recently, the use of numerical simulation has been increased rapidly because of the development of high performance computer systems. The present study is aimed to investigate flow characteristics of a two-dimensional sharp plane. Unsteady calculation by FDM(Finite Difference Method) based upon SOLA scheme which was performed at $Re=2{\times}10^4$in viscous incompressible flow within a finite domain on the irregular grid formation. Total numbers of irregular grids are $8{\times}10^4$. The minimum grid size is 1/100 of the plane length L which is the representative length. The inclined angles of every objects are $15^{\circ}, \;30^{\circ}\;and\; 45^{\circ}.$ And, the edge angle of the plane is $30^{\circ}.$ This study discussed the flow characteristics in term of the turbulent intensity, vorticity and frequency analysis. Developed flows show that the periodic Karman vortices occur at the back of the plane.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.17 no.3
• /
• pp.333-339
• /
• 2015

In this paper, the surface settlement resulted from the shallow TBM tunnelling has been numerically simulated. TBM tunnelling is especially used in urban area to avoid serious vibration and noise caused by explosion in NATM. Surface settlement is one of the most important problems encountered in all tunnelling and critical in urban areas. In this study, face vibration of TBM excavation is considered to estimate surface settlement trend according to TBM extrusion. The dynamic excavation forces are calculated by total torque on the TBM cutterhead in mixed-face of soil and weathered rock condition with shallow depth. A 3-dimensional FDM code is employed to simulate TBM tunnelling and mechanical-dynamic coupling analysis is performed. The 3D numerical analysis results showed that dynamic settlement histories and trend of surface settlement successfully. The maximum settlement occurred at the excavation point located at 2.5D behind the face, and the effect of face vibration on the surface settlement was verified in this study.

• Journal of Korean Society of Environmental Engineers
• /
• v.30 no.11
• /
• pp.1146-1153
• /
• 2008

Numerical simulation was carried out to study the trichloroethylene (TCE) degradation by permeable reactive barrier (PRB), and revealed the effect of concentration of TCE, iron medium mass, and concentration of iron-reducing bacteria (IRB). Newly developed model was based on axial dispersion reactor model with chemical and biological reaction terms and was implemented using MATLAB ver R2006A for the numerical solutions of dispersion, convection, and reactions over column length and elapsed time. The reaction terms include reactions of TCE degradation by zero-valent iron (ZVI, Fe$^0$) and ferrous iron (Fe$^{2+}$). TCE concentration in the column inlet was maintained as 10 mg/L. Equation for Fe$^0$ degradation includes only TCE reaction term, while one for Fe$^{2+}$ has chemical and biological reaction terms with TCE and IRB, respectively. Two coupled equations eventually modeled the change of TCE concentration in a column. At Fe$^0$ column, TCE degradation rate was found to be more than 99% from 60 hours to 235 hours, and declined to less than 1% in 1,365 hours. At the Fe$^{2+}$ and IRB mixed column, TCE degradation rate was equilibrated at 85.3% after 210 hours and kept it constant. These results imply that the ferrous iron produced by IRB has lowered the TCE degradation efficiency than ZVI but it can have higher longevity.http://kci.go.kr/kciportal/ci/contents/ciConnReprerSearchPopup.kci#

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.9
• /
• pp.796-805
• /
• 2016

As a part of the first stage in the helicopter flight simulator development, this study numerically evaluates handling qualities of the dynamics model. The flight dynamics model was generated using public information for AS365 N2, the target aircraft of the simulator. The flight simulator is under development as a pilot training and research tool for firefighting missions. The assessment of the model intends to validate general characteristics and suitability before the model is enhanced with flight test data. The evaluation is based on the ADS-33E-PRF(Aeroautical Design Standard Performance Specification Handling Qualities Requirement) criteria, with consideration of category of the aircraft, missions, and environment. The numerical operations follow required or recommended procedures of flight test for compliance demonstration. Evaluation results are evaluated according to the rating specified in maneuverability ADS-33E-PRF. Results have identified to provide a satisfactory platform for flight dynamic model in the general helicopter simulator generated based on the RotorLibFDM, and can be used as a base for basic training and research.

• Journal of the Korea Organic Resources Recycling Association
• /
• v.28 no.4
• /
• pp.43-52
• /
• 2020

As the membrane gas separation technology grows, various models were developed by numerous researchers to describe the separation process. In this work, the counter-current model was compared thoroughly with experimental data. Experimentally, hollow fiber membrane using CA module was prepared for the separation of biogas. The pure gas permeation properties of membrane module for methane, nitrogen, oxygen, and carbon dioxide were measured. The permeance of CO2 and CH4 were 25.82 GPU and 0.65 GPU, respectively. The high CO2/CH4 selectivity of 39.7 was obtained. the separation test for three different simulated mixed gases were carried out after pure gas test, and the gas concentration of the permeate at various stage-cut were measured from CA membrane module. Results showed that the experimental data agreed with the numerical simulation. A mathematical model has implemented in this study for the separation of biogas using a membrane module. The finite difference method (FDM) is applied to calculate the membrane biogas separation behaviors. Futhermore, the counter-current model can be considered as a convenient model for biogas separation process.

• Journal of the Korean GEO-environmental Society
• /
• v.13 no.6
• /
• pp.59-66
• /
• 2012

The purpose of this study is to estimate the behavior and the mechanism of debris flow on the slope, which has specially various gradient plane. The numerical simulation was performed by using the Finite Differential Element method (FDM) based on the equation for the mass conservation and momentum conservation. The mechanism of flow type for debris flow is divided into three flow types which are stony debris flow, immature debris flow, and turbulent water flow, respectively. First, flow discharge, water flow depth, sediment volume concentration was investigated by variable input of flow discharge at the straight slope angle and two step inclined plane. As the input of flow discharge was decrease, flow discharge and water flow depth was increased, after the first coming debris flow only reached at the downstream. As the input of flow discharge was increased, the curve of flow discharge and flow depth was highly fluctuated. As the results of RMS ratio, the flow discharge and flow depth was lower two step slope angle than the straight slope angle. Second, the behavior of debris flow was investigated by the four cases of gradient degree at the downstream of slope angle. The band width of flow discharge and flow depth for $14^{\circ}$ between $16^{\circ}$ was higher than other gradient degree, and fluctuation curve was continuously high after 10 seconds.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.2 no.3
• /
• pp.71-81
• /
• 2000

For the construction of NATM tunnel, it is required a design based on the accurate soil condition from soil investigation. However, in practice, it often designs tunnels without fully understanding the condition. Especially, when soft soil comes up, or ground water breaks out suddenly on the construction, it needs to secure the stability of tunnel by appropriate reinforcing construction according to the results of measurements on field superlatively reflecting the faced situation. This report reviews the mostsuitable stability of tunnel in the construction of soft soil of tunnel by numerical analysis using FDM after re-evaluated the soil properties through back analysis using the results of measurements to simulate abruptly occurred deformation. And applying steel pipe grouting row by row on the wall and the low part of tunnel and also applying the construction method of temporary invert after excavation of the upper part of tunnel, the excavation of soft soil tunnel secured the structural stability of tunnel has been completed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.3
• /
• pp.279-285
• /
• 2011

We perform a numerical study to determine the time of onset of natural convection in a transient hot wire (THW) device for measuring the thermal conductivity of nanofluids. The samples used in this simulation are water-based $Al_2O_3$ nanofluids with volume fractions of 1%, 4%, and 10%, and the properties are calculated by theoretical models and experimental correlations. The THW apparatus using coated wire is modeled by the control-volume-based finite difference method, and the start of natural convection is determined by observing the temperature rise of the wire under a gravity field. The onset time is 11.5 s for water and 41.6 s for water-based $Al_2O_3$ nanofluids predicted by Maxwell thermal conductivity model with a 10% volume fraction. We confirm that the onset time of natural convection of nanofluids in the cylinder increases with the nanoparticle volume fraction. We suggest a correlation for predicting the onset time on the basis of the numerical results. Finally, it is shown that the measurement error due to natural convection is negligible if the measurement using the transient hot wire method is completed before the onset of natural convection in the base fluid.