• Tunnel and Underground Space
• /
• v.15 no.2 s.55
• /
• pp.145-156
• /
• 2005

MODFLOW, 3-D finite difference code, is widely used to model groundwater flow and has been used to assess the effect of excavations on the groundwater system due to construction of subways and mountain tunnels. The results of numerical analysis depend on boundary conditions, initial conditions, conceptual models and hydrogeological properties. Therefore, its accuracy can only be enhanced using more realistic and field oriented input parameters. In this study, SA(simulated annealing) was used to integrate hydraulic conductivities from a few of injection tests with geophysical reference images. The realizations of hydraulic conductivity random field are obtained and then groundwater flows in each geostatistically equivalent media are analyzed with a numerical simulation. This approach can give probabilistic results of groundwater flow modeling considering the uncertainty of hydrogeological medium. In other words, this approach makes it possible to quantify the propagation of uncertainty of hydraulic conductivities into groundwater flow.

• Tunnel and Underground Space
• /
• v.26 no.3
• /
• pp.181-191
• /
• 2016

In geological storage of $CO_2$, the behavior of $CO_2$ is influenced by pore network of rock. In this study, the drilling cores from Pohang Basin were analyzed quantitatively using three-dimensional images acquired by X-ray micro computed tomography. The porosities of sandstone specimens around 740 m-depth (T1), 780 m-depth (T2) and 810 m-depth (T3) which were target strata were 25.22%, 23.97%, 6.28%, respectively. Equivalent diameter, volume, area, local thickness of pores inside the sandstone specimens were analyzed. As a result, the microstructural properties of T1 and T2 specimens were more suitable for geological storage of $CO_2$ than those of T3 specimens. The result of the study can be used as input data of the site for decision of injection condition, flow simulation and so on.

• Journal of Advanced Marine Engineering and Technology
• /
• v.19 no.1
• /
• pp.60-70
• /
• 1995

In the field of marine transportation the energy saving is one of the most important factors for profit. In order to reduce the fuel oil consumption the ship's propulsion efficiency must be increased as much as possible. The propulsion efficiency depends upon a combination of an engine and a propeller. The propeller has better efficiency as lower rotational speed. This situation led the engine manufacturers to design the engine that has lower speed, longer stroke and a small number of cylinders. Consequently the variation of rotational torque became larger than before because of the longer delay-time in the fuel oil injection process and an increased output per cylinder. As this new trends the conventional mechanical-hydrualic governors for engine speed control have been replaced by digital speed controllers which adopted the PID control or the optimal control algorithm. But these control algorithms have not enough robustness to suppress the variation of the delay-time and the parameter pertubation. In this paper we consider the delay-time and the perturbation of engine parameters as the modeling uncetainties. Next we design the controller which has zero offset in steady state engine speed, based on the two-degree-of-freedom control theory and $\mu$-synthesis. Thd validity of the controller is investigated through the response simulation. We use a personal computer and an analog computer as the digital controller and the engine (plant) part respectively. And, we certify that the designed controller maintains its performance even though the engine parameters may vary.

• Journal of ILASS-Korea
• /
• v.25 no.4
• /
• pp.162-169
• /
• 2020

The worldwide focus on reducing the emissions, fuel and lubricant consumption in T-GDI engines is leading engineers to consider the crankcase ventilation and oil mist separation system as an important means of control. In today's passenger cars, the oil mist separation systems mainly use the inertia effect (e.g. labyrinth, cyclone etc.). Therefore, this study has investigated high efficiency cylinder head-integrated oil-mist separator by using a compact multi-impactor type oil mist separator system to ensure adequate oil mist separation performance. For this purpose, engine dynamometer testing with oil particle efficiency measurement equipment and 3D two-phase flow simulation have been performed for various engine operating conditions. Tests with an actual engine on a dynamometer showed oil aerosol particle size distributions varied depending on operating conditions. For instance, high rpm and load increases bot only blow-by gases but the amount of small size oil droplets. Submicron-sized particles (less than 0.5 ㎛) were also observed. It is also found that the impactor type separator is able to separate nearly no droplets of diameter lower than 3 ㎛. CFD results showed that the complex aerodynamics processes that lead to strong impingement and break-up can strip out large droplets and generate more small size droplets.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.10
• /
• pp.709-716
• /
• 2022

In this paper the Gamma guidance law, which was used for IUS (Inertial Upper Stage), is applied for solid-motor guidance of a upper stage of a satellite launch vehicle. The RCS (Reaction Control System), which activates after burnout of the upper stage, is employed for the convergence of the guidance algorithm and compensation of velocity errors induced by the solid motor. The algorithm is also simplified by replacing the time-consuming numerical integration process to predict final vehicle states with Keplerian trajectories. The performance of the algorithm is evaluated by conducting 3-DOF computer simulations for off-nominal flight conditions. The numerical results show that Gamma guidance can reduce the orbit injection accuracy in comparison with that obtained by applying open-loop commands.

• Journal of Satellite, Information and Communications
• /
• v.12 no.3
• /
• pp.80-85
• /
• 2017

In this paper, a mircrowave breakdown of X-band circular waveguide cavity filter, which occurred during ground test, was introduced, and electro-magnetic field simulation results to identify a root cause, and the analysis of possibility of its occurrence on orbit operation were presented. Filter modeling for simulation was conducted with a commercial tool (FEST3D), and electric fields inside the filter were monitored at the input of 1 W continuous wave. In our observation, strong electric field intensities were monitored on the tuning screws especially at the input of band-edge frequencies. The threshold power levels for the breakdown were also estimated and compared with the input power levels actually injected to the filter. From this estimation, we could figure out that the power exceeding the breakdown threshold was injected to the filter so that strong electric fields were generated and temperature increased high, and this became a root cause of the electrical short. Our further analysis showed that this kind of microwave breakdown is not likely to occur on orbit operation, and multipactor is expected not to occur at the input of band-edge frequencies. As a measure to prevent the microwave breakdown, we suggested to avoid the injection of band-edge frequencies and inject lower power levels to the filter.

• The Korean Journal of Rheology
• /
• v.10 no.4
• /
• pp.202-216
• /
• 1998

Improved version of previous 'Orthotropic' closure approximation, termed 'ORW' has been numerically developed using new homogeneous flow data. Previous 'Orthotropic' closure approximation, i.e., ORF or ORL showed non-physical oscillation for interaction coefficient $C_1$<0.001 at simple shear flow. It also shows non-physcial oscillation and under-prediction compared with 'Distribution Function Calculation' at non-homogeneous flow of center-gated disk. These phenomena are mainly due to the flow data of 'Distribution Function Calculation' which were used for least-square optimization. ORW obtained by fitting flow data of low interaction coefficient does not show non-physical oscillation and results in reasonably good behaviors at non-homogeneous flows as well as homogeneous flows. Fitting function forms have not been found to improve overall behaviors. It has been found that considering all the eigenvalues of orientation tensor (including the third eigenvalues) might end up with a better closure approximation than just considering the first and second eigenvalues. It is, however, very important and yet difficult to select appropriate function forms of eigenvalues. Numerical simulation including coupling and in-plane velocity gradient effects were performed for injection mold filing process with a film-gated strip and a center-gated disk using ORW and various other closure approximations for comparisons. Although ORW is in excellent agreement with 'Distribution Function Calculation', the predicted results seem to have consistent error in comparison with experimental data. The diffusivity term with constant interaction coefficient might have to be further investigated in order to accurately describe orientation states.

• Tunnel and Underground Space
• /
• v.30 no.5
• /
• pp.446-461
• /
• 2020

Numerical simulations of CIEMAT column test in Spain are performed to investigate the coupled thermo-hydro-mechanical (THM) behavior of MX80 bentonite pellets using TOUGH2-FLAC3D. The heater power and injection pressure of water in the numerical simulations are identical to those in the laboratory test. To investigate the applicability of the thermo-hydraulic (TH) model used in TOUGH2 code to prediction of the coupled TH behavior, the simulation results are compared with the observations of temperature and relative humidity with time. The tendencies of the coupled behavior observed in the test are well represented by the numerical models and the simulator in terms of temperature and relative humidity evolutions. Moreover, the performance of the models for the reproduction and prediction of the coupled TH behavior is globally satisfactory compared with the observations. However, the calculated stress change is relatively small and slow due to the limitations of the simple elastic and swelling pressure model used in numerical simulations. It seems that the two models are insufficient to realistically reproduce the complex coupled THM behavior in the bentonite pellets.

• Journal of the Korean Vacuum Society
• /
• v.22 no.2
• /
• pp.66-78
• /
• 2013

In order to develop of 80 kW RF plasma torch system, we achieved three-dimensional simulations for the extraction of more information as temperature in torch and fluid behavior analysis, etc. The position of powder injection tube, the plasma discharge characteristics with various input current and various length of ceramic tube, and the plasma temperature characteristics with process gas flow rate such those was simulated. RF thermal plasma torch designed by simulation was manufactured that was measured to the maximum of 89.3 kW power. The mass production using developed 80 kW RF thermal plasma torch system were investigated by characteristics manufactured of Si nano powder. The mass-production level of Si nano-powder was average of 539 g/hr and high yield rate of 71.6%, respectively. The particle size distribution $D_{99}/D_{50}$ of manufacturing nano-powder was investigated to 1.98 as a good uniform.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.4
• /
• pp.343-351
• /
• 2016

Cryogenic spray characteristics of a nitrogen swirl injector operating in supercritical environment have been numerically investigated. By comparing the equation of states(EOS) used for supercritical condition, SRK EOS was applied to predict the nitrogen thermodynamic property under supercritical environment. A Chung's method was implemented for the calculation of viscosity and conductivity and Takahashi's correlation based on Fuller's Theorem was implemented for the calculation of diffusion coefficient. By injecting the nitrogen with 5 bar differential pressure into 50 bar chamber filled with nitrogen, numerical simulation has been conducted. The dynamic Smagorinsky sub-grid scale (SGS) model has been compared with the algebraic Smagorinsky SGS model using FFT frequency analysis. The instability at the liquid film and gas core inside injector and the propagation of pressure oscillation into the injector has been investigated. The spreading angle of swirl injector obtained by numerical calculation has been validated with experimental result.