• Journal of Advanced Navigation Technology
• /
• v.26 no.6
• /
• pp.504-509
• /
• 2022

In this paper, we deal with the stability condition of linear time-varying interval discrete systems with time-varying delays and unstructured uncertainty. For the time-varying interval discrete system which has interval matrix as its system matrices, time-varying delay time within some interval value and unstructured uncertainty which can include non-linearity and be expressed by only its magnitude, the stability condition is proposed. Compared with the previous result derived by using a upper bound solution of the Lyapunov equation, the new result is derived by the form of simple inequality based on Lyapunov stability condition and has the advantage of being more effective in checking stability. Furthermore, the proposed condition is very comprehensive, powerful and inclusive the previously published conditions of various linear discrete systems, and can be expressed by the terms of magnitudes of the time-varying delay time and uncertainty, and bounds of interval matrices. The superiority of the new condition is shown in the derivation, and the usefulness and advantage of the proposed condition are examined through numerical example.

• Journal of Advanced Navigation Technology
• /
• v.25 no.6
• /
• pp.551-556
• /
• 2021

In this paper, we deal with the stability condition of linear interval discrete systems with time-varying delays and unstructured uncertainty. For the interval discrete system which has interval matrix as its system matrices, time-varying delay time within some interval value and unstructured uncertainty which can include non-linearity and be expressed by only its magnitude, the stability condition is proposed. Compared with the previous result derived by using a upper bound solution of the Lyapunov equation, the new results are derived by the form of simple inequality based on Lyapunov stability condition and have the advantage of being more effective in stability application. Furthermore, the proposed stable conditions are very comprehensive and powerful, including the previously published stable conditions of various linear discrete systems. The superiority of the new condition is proven in the derivation process, and the utility and superiority of the proposed condition are examined through numerical example.

• Journal of Advanced Navigation Technology
• /
• v.27 no.6
• /
• pp.871-876
• /
• 2023

In this paper, we deal with the stable conditions when two uncertainties exist simultaneously in a linear discrete time-varying interval system with time-varying delay time. The interval system is a system in which system matrices are given in the form of an interval matrix, and this paper targets the system in which the delay time of these interval system matrices and state variables is time-varying. We propose the system stability condition when there is simultaneous unstructured uncertainty that includes nonlinearity and only its magnitude and uncertainty in the system matrix of delayed state variables. The stable bounds for two types of uncertainty are derived as an analytical equation. The proposed stability condition and bounds can include previous stability condition for various linear discrete systems, and the values such as time-varying delay time variation size, uncertainty size, and range of interval matrix are all included in the conditional equation. The new bounds of stability are compared with previous results through numerical example, and its effectiveness and excellence are verified.

• Journal of Korean Association for Spatial Structures
• /
• v.24 no.3
• /
• pp.43-51
• /
• 2024

Dynamic responses of nuclear power plant structure subjected to earthquake loads should be carefully investigated for safety. Because nuclear power plant structure are usually constructed by material of reinforced concrete, the aging deterioration of R.C. have no small effect on structural behavior of nuclear power plant structure. Therefore, aging deterioration of R.C. nuclear power plant structure should be considered for exact prediction of seismic responses of the structure. In this study, a machine learning model for seismic response prediction of nuclear power plant structure was developed by considering aging deterioration. The OPR-1000 was selected as an example structure for numerical simulation. The OPR-1000 was originally designated as the Korean Standard Nuclear Power Plant (KSNP), and was re-designated as the OPR-1000 in 2005 for foreign sales. 500 artificial ground motions were generated based on site characteristics of Korea. Elastic modulus, damping ratio, poisson's ratio and density were selected to consider material property variation due to aging deterioration. Six machine learning algorithms such as, Decision Tree (DT), Random Forest (RF), Support Vector Machine (SVM), K-Nearest Neighbor (KNN), Artificial Neural Networks (ANN), eXtreme Gradient Boosting (XGBoost), were used t o construct seispic response prediction model. 13 intensity measures and 4 material properties were used input parameters of the training database. Performance evaluation was performed using metrics like root mean square error, mean square error, mean absolute error, and coefficient of determination. The optimization of hyperparameters was achieved through k-fold cross-validation and grid search techniques. The analysis results show that neural networks present good prediction performance considering aging deterioration.

• Journal of the Society of Naval Architects of Korea
• /
• v.61 no.5
• /
• pp.289-296
• /
• 2024

Recently, with the acceleration of global warming, the importance of carbon neutrality is being emphonasized. In response to this, various technologies are being developed to achieve carbon neutrality, with Carbon Capture, Utilization, and Storage (CCUS) being a prominent example. Research has been conducted on the injection technology to environmentally discharge carbon captured by using one of the carbon capture technologies, HAK-CRS. In this study, numerical simulations were performed using commercial software, ANSYS Fluent, to understand the dispersed multiphase flow between a water jet and CaCO₃ particles. The analysis focused on the differences in particle behavior when injecting CaCO₃ with different flow rates. It was observed that as the mass flow rate of CaCO₃ increased, there was a tendency for the particles to deviate from the flow of the jet and rapidly fall, influenced significantly by gravity. These results indicate that, for CaCO₃ particles to disperse widely, the flow rate should not be excessively high. Given the potential adverse impacts on marine ecosystems due to the high density of CaCO₃, research on injection technology should also be conducted.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.

• Journal of the Society of Disaster Information
• /
• v.9 no.4
• /
• pp.449-461
• /
• 2013

As the hydrofluoric acid leak in Gumi-si, Gyeongsangbuk-do or hydrochloric acid leak in Ulsan, Gyeongsangnam-do demonstrated, chemical related accidents are mostly caused by large amounts of volatile toxic substances leaking due to the damages of storage tank or pipe lines of transporter. Safety assessment is the most important concern because such toxic material accidents cause human and material damages to the environment and atmosphere of the surrounding area. Therefore, in this study, a hydrofluoric acid leaked from a storage tank was selected as the study example to simulate the leaked substance diffusing into the atmosphere and result analysis was performed through the numerical Analysis and diffusion simulation of ALOHA(Areal Location of Hazardous Atmospheres). the results of a qualitative evaluation of HAZOP (Hazard Operability)was looked at to find that the flange leak, operation delay due to leakage of the valve and the hose, and toxic gas leak were danger factors. Possibility of fire from temperature, pressure and corrosion, nitrogen supply overpressure and toxic leak from internal corrosion of tank or pipe joints were also found to be high. ALOHA resulting effects were a little different depending on the input data of Dense Gas Model, however, the wind direction and speed, rather than atmospheric stability, played bigger role. Higher wind speed affected the diffusion of contaminant. In term of the diffusion concentration, both liquid and gas leaks resulted in almost the same $LC_{50}$ and ALOHA AEGL-3(Acute Exposure Guidline Level) values. Each scenarios showed almost identical results in ALOHA model. Therefore, a buffer distance of toxic gas can be determined by comparing the numerical analysis and the diffusion concentration to the IDLH(Immediately Dangerous to Life and Health). Such study will help perform the risk assessment of toxic leak more efficiently and be utilized in establishing community emergency response system properly.

• Journal of Korea Water Resources Association
• /
• v.48 no.12
• /
• pp.1065-1075
• /
• 2015

In the estuary where the structure such as river-mouth weir has been installed, the flow is developed very complicatedly due to river water from upstream, tide of the sea and floodgate operation. Especially, if basin outlets more than one exists in one estuary, the boundary conditions will be significantly more complex form. Saemangeum(SMG) project area in Korea is the most typical example. There are Mankyung river and Dongjin river in upstream. The water of them inflows into SMG project area. In the downstream, river flow was drained from inland to sea over the SMG sea dike through the sluice. The connecting channel was located between Mankyung and Dongjin basins. It functions not only as transportation by ship in ordinary period but also as flood sharing by sending flood flow to each other in flood period. Therefore, in order to secure the safety against flood, it is very important to understand the flood sharing capacity for connecting channel. In this study, the flood control effect was analyzed using numerical simulation. Delft3D was used to numerical simulation and simulated period was set up with neap tide, in which the maximum flood stage occurred due to poor drainage. Actually, three connecting channels were designed in land use plan of the SMG Master Plan, but they were simplified to a single channel for conciseness of analysis in this study. According to the results of numerical analysis, the water level difference between two basins was increased and the maximum flood stage at dike sluice was also upraised depending on decrease of conveyance. And the velocity induced by same water level difference was decreased when the conveyance became smaller. In certain conveyance above, there was almost no flood control effect. Therefore, if the results of this study are considered for design of connecting channel, it will be expected to draw the optimal conveyance for minimizing dredging construction cost while maximizing the flood control effect.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.11 no.4
• /
• pp.327-351
• /
• 2009

Laboratory tests are performed in this paper to investigate the brittle failure characteristics of over-stressed rocks taken in deep depth. Also, numerical simulation performed using that the so-called CWFS(Cohesion Weakening Frictional Strengthening) model is known to predict brittle failure phenomenon reasonably well. The most typical rock types of Korean peninsula - granite and gneiss - were used for testing. Results of uniaxial compression tests showed that the crack initiation stress was about 41 % to 42% of the uniaxial compressive strength regardless of rock types, where as, the crack damage stress of granite was about 75%, and that of gneiss was about 97%. Through the damage-controlled test, strength parameters of each rock were obtained as a function of damage degree. After the peak, the crack damage stress and the maximum stress were decreased, The cohesion was decreased and the friction angle was increased with increase of rock damage. Before reaching the peak, the elastic modulus was slightly increased, while decreased after the peak. Poisson's ratio was increased as the damage of rock proceeds. Comparison of uniaxial compression tests and damage-controlled tests shows the crack initiation stress estimated from the damage-controlled test fluctuated within the range of crack initiation stress obtained from the uniaxial compression test; the crack damage stress was less than that estimated from the uniaxial compression test. In order to predict the critical depth that brittle failure occurs, numerical simulations using the CWFS model were performed for an example site. Material parameters obtained from the laboratory tests mentioned above were used for CWFS simulation. Comparison between the critical depth predicted from the numerical simulation using the CWFS model and that predicted by using the damage index proposed by Martin et al.(l999), showed that critical depth cannot be reasonably predicted by the currently used damage index except for circular tunnels. A modified damage index was proposed by the author which takes the shape of tunnels other than circular into account.

• The Korean Journal of Rheology
• /
• v.11 no.2
• /
• pp.82-90
• /
• 1999

Electrorheological(ER) fluid is a material that shows the dramatic change of rheological properties under an electric field and responds reversibly in a few milliseconds. ER fluid's response to an electric field along with its fast switching capability allows ER devices to be precisely controlled. The real application with ER fluid, however, has many limitations to be overcome; temperature fluctuation, moisture, dust, aggregation, precipitation, and low yield stress, for example. The magnitude and the characteristics of yield stress of ER fluid plays an important role in practical applications. In this research, a dynamic simulation on the squeezing flow of the ER fluid was carried out. Numerical simulation on isolated chains was performed to find out the effect of hydrodynamic and electrostatic force depending on the chain location, the squeezing rate, and the chain structure. Suspension model that is composed of a large number of particles was also investigated. The increase of normal stresses as well as the existence of a yield stress at an earlier stage could be observed, and the effective control of the normal stresses could be achieved at an optimal condition of the hydrodynamic force and the electrostatic force.