• Geomechanics and Engineering
• /
• v.22 no.1
• /
• pp.97-108
• /
• 2020

The time-consuming and less objectivity are the main problems of conventional micromechanical parameters calibration method of Particle Flow Code simulations. Thus this study aims to address these two limitation of the conventional "trial-and-error" method. A new calibration method for the linear parallel bond model (CM-LPBM) is proposed. First, numerical simulations are conducted based on the results of the uniaxial compression tests on limestone. The macroscopic response of the numerical model agrees well with the results of the uniaxial compression tests. To reduce the number of the independent micromechanical parameters, numerical simulations are then carried out. Based on the results of the orthogonal experiments and the multi-factor variance analysis, main micromechanical parameters affecting the macro parameters of rocks are proposed. The macro-micro parameter functions are ultimately established using multiple linear regression, and the iteration correction formulas of the micromechanical parameters are obtained. To further verify the validity of the proposed method, a case study is carried out. The error between the macro mechanical response and the numerical results is less than 5%. Hence the calibration method, i.e., the CM-LPBM, is reliable for obtaining the micromechanical parameters quickly and accurately, providing reference for the calibration of micromechanical parameters.

• Journal of Advanced Marine Engineering and Technology
• /
• v.28 no.5
• /
• pp.818-826
• /
• 2004

The spray-induced mixing characteristics and thermal decomposition of aqueous urea solution into ammonia have been studied to design optimum sizes and geometries of the mixing chamber in SCR(Selective Catalytic Reduction) system. The cold flow tests about the urea-injection nozzle were performed to clarify the parameters of spray mixing characteristics such as mean diameter and velocity of drops and spray width determined from the interactions between incoming air and injected drops. Discrete particle model in Fluent code was adopted to simulate spray-induced mixing process and the experimental results on the spray characteristics were used as input data of numerical calculations. The simulation results on the spray-induced mixing were verified by comparing the spray width extracted from the digital images with the simulated Particle tracks of injected drops. The single kinetic model was adopted to predict thermal decomposition of urea solution into ammonia and solved simultaneously along with the verified spray model. The hot air generator was designed to match the flow rate and temperature of the exhaust gas of the real engines The measured ammonia productions in the hot air generator were compared with the numerical predictions and the comparison results showed good agreements. Finally, we concluded that the design capabilities for sizing optimum mixing chamber were established.

• Journal of Korean Society for Atmospheric Environment
• /
• v.27 no.6
• /
• pp.681-691
• /
• 2011

The size distributions of smoke particles from fire are prerequisite for the studies on fire detection and adverse health effects. Above the flame of the fire, coagulation dominates and the smoke particles grow from 1 to 50 nm up to 100 to 3,000 nm, sizes ranging from the free-molecular regime to the continuum regime. The characteristics of the agglomeration of the smoke particles are well known, independently for each of the free-molecular and continuum regimes. However, there are not many systematic studies in the entire regime by the complexity of the mechanisms. The purpose of this work is to find the characteristics of the development of the size distribution of smoke particles by agglomeration in the entire size range covering the free-molecular regime, via transition regime, to the near-continuum and continuum regime for each variation of parameters such as fractal dimension, primary particle size and dimensionless coagulation time. In this work, the dynamic equation for the discrete-size spectrum of the particles was solved using a nodal method based on the modification of a sectional method. In the calculation, the collision frequency function for the entire regime, which is derived by using the concept of collision volume and general enhancement function, was applied. The self-preserving size distribution for the entire regime is compared with the ones for the free-molecular or continuum regimes for each variation of the parameters.

• Journal of the Korean Geotechnical Society
• /
• v.33 no.7
• /
• pp.43-53
• /
• 2017

Road subsidence, frequently occurring in urban areas, is caused by collapsing of surface layer due to underground cavities followed by a loss of soils. To better understand this phenomenon, the mechanism of cavity formation should be identified firstly. Two kinds of possible subsidence mechanisms were established through previous case studies and the numerical analyses based on Distinct Element Method were conducted for each of these mechanisms. It was confirmed that particle loss and surface settlement can develop differently depending on slit size, void ratio, and particle shape among the various factors influencing the road subsidence. The result demonstrated that the effects of varying cavity diameter and depth could be quantified as a damage chart.

• Wind and Structures
• /
• v.27 no.4
• /
• pp.213-221
• /
• 2018

In the desert and Gobi regions with strong wind and large sediment discharge, sand transporting engineering is more effective than sand blocking and sand fixing measures in sand prevention. This study uses the discrete phase model of 3D numerical simulation to study the motion trail, motion state and distribution rule of sand particles with different grain diameters when the included angle between the main shaft of the feather-row lateral transportation sand barrier and the wind direction changes, and conducts a comparison in combination with the wind tunnel test and the flow field rule of common sand barrier. According to the comparison, when wind-sand incoming flow passes through a feather-row sand barrier, sand particles slow down and deposit within the deceleration area under the resistance of the feather-row sand barrier, move along the transportation area formed by the transportation force, and accumulate as a ridge at the tail of the engineering. With increasing wind speed, the eolian erosion of the sand particles to the ground and the feather-row sand barrier is enhanced, and the sand transporting quantity and throw-over quantity of the feather-row sand barrier are both increased. When sand particles with different grain diameters bypass the feather-row sand barrier, the particle size of the infiltrating sands will increase with the included angle between the main shaft of the feather-row sand barrier and the wind direction. The obtained result demonstrates that, at a constant wind speed, the flow field formed is most suitable for the lateral transportation of the wind-drift flow when the included angle between the main shaft of the feather-row sand barrier lateral transportation engineering and the wind speed is less than or equal to $30^{\circ}$.

• Advances in concrete construction
• /
• v.13 no.1
• /
• pp.83-99
• /
• 2022

The two, NBD and SCB tests using gypsum circular discs each containing a single notch have been experimentally accomplished in a rock mechanics laboratory. These specimens have also been numerically modelled by a two-dimensional particle flow which is based on Discrete Element Method (DEM). Each testing specimen had a thickness of 5 cm with 10 cm in diameter. The specimens' lengths varied as 2, 3, and 4 cm; and the specimens' notch angles varied as 0°, 45° and 90°. Similar semi-circular gypsum specimens were also prepared each contained one edge notch with angles 0° or 45°. The uniaxial testing machine was used to perform the experimental tests for both NBD and SCB gypsum specimens. At the same time, the numerical simulation of these tests were performed by PFC2D. The experimental results showed that the failure mechanism of rocks is mainly affected by the orientations of joints with respect to the loading directions. The failure mechanism and fracturing patterns of the gypsum specimens are directly related to the final failure loading. It has been shown that the number of induced tensile cracks showing the specimens' tensile behavior, and increases by decreasing the length and angle of joints. It should be noted that the fracture toughness of rocks' specimens obtained by NBD tests was higher than that of the SCB tests. The fracture toughness of rocks usually increases with the increasing of joints' angles but increasing the joints' lengths do not change the fracture toughness. The numerical solutions and the experimental results for both NDB and SCB tests give nearly similar fracture patterns during the loading process.

• Bulletin of the Korean Chemical Society
• /
• v.24 no.6
• /
• pp.771-776
• /
• 2003

The structural and dynamic properties of small n-alkane clusters embedded in a mesoscopic solvent are investigated. The solvent interactions are taken into account through a multi-particle collision operator that conserves mass, momentum and energy and the solvent dynamics is updated at discrete time intervals. The cluster molecules interact among themselves and with the solvent molecules through intermolecular forces. The properties of n-heptane and n-decane clusters interacting with the mesoscopic solvent molecules through repulsive Lennard-Jones interactions are studied as a function of the number of the mesoscopic solvent molecules. Modifications of both the cluster and solvent structure as a result of cluster-solvent interactions are considered. The cluster-solvent interactions also affect the dynamics of the small n-alkane clusters.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.106-112
• /
• 2004

A kinetic theory analysis is made of low-speed gas flows around a micro-plate. The Boltzmann equation simplified by a collision model is solved by means of a finite difference approximation with the discrete ordinate method. The method does not suffer from statistical noise which is common in particle based methods and requires much less amount of computational effort. Calculations are made for flows around a micro-scale flat plate with a finite length of 20 microns. The method is assessed by comparing the results with those from several different methods and available experimental data.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.1-7
• /
• 2004

A kinetic theory analysis is made of low-speed rarefied gas flows around a flat plate. The Boltzmann equation simplified by a collision model is solved by means of a finite difference approximation with the discrete ordinate method. The method does not suffer from statistical noise which is common in particle based methods and requires much less amount of computational effort. Calculations are made for flows around a micro-scale flat plate with a finite length of 20 microns. The method is assessed by comparing the results with those from several different methods and available experimental data.

• ETRI Journal
• /
• v.45 no.3
• /
• pp.433-447
• /
• 2023

Critical issues such as connection congestion, long transmission delay, and packet loss become even worse during epidemic, disaster, and so on. In this study, a link load balancing method is proposed to address these issues on the data plane, a plane of the software-defined network (SDN) architecture. These problems are NP-complete, so a meta-heuristic approach, discrete particle swarm optimization, is used with a novel hybrid cost function. The superiority of the proposed method over existing methods in the literature is that it provides link and switch load balancing simultaneously. The goal is to choose a path that minimizes the connection load between the source and destination in multipath SDNs. Furthermore, the proposed work is dynamic, so selected paths are regularly updated. Simulation results prove that with the proposed method, streams reach the target with minimum time, no loss, low power consumption, and low memory usage.