• The KSFM Journal of Fluid Machinery
• /
• v.15 no.1
• /
• pp.27-35
• /
• 2012

A numerical simulation for particle collection efficiency in a wire-plate electrostatic precipitator (ESP) has been performed. Method of characteristics and finite differencing method (MOC-FDM) were employed to obtain electric field and space charge density, and lattice boltzmann method (LBM) was used to predict the Electrohydrodynamic (EHD) flow according to the ion convection. Large eddy simulation (LES) was considered for turbulent flow and particle simulation was performed by discrete element method (DEM) which considered field charging, electric force, drag force and wall-collision. One way coupling from FDM to LBM was used with small and low density particle assumption. When the charged particle collided with the collecting plate, particle-wall collision was calculated for re-entertainment effect and the effect of gravity force was considered.

• 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.

• Korean Journal of Agricultural Science
• /
• v.50 no.3
• /
• pp.457-468
• /
• 2023

The purpose of this study is to develop a high-speed combine harvester. The performance was evaluated by composing a dynamic simulation model of a threshing cylinder and analyzing the amount of threshed rice grain during threshing operations. The rotational speed of the threshing cylinder was set at 10 rpm intervals from 500 rpm until 540 rpm, based on the rated rotational speed of 507 rpm. The rice stem model was developed using the EDEM software using measured rice stem properties. Multibody dynamics software was utilized to model the threshing cylinder and tank comprising five sections below the threshing cylinder, and the threshing performance was evaluated by weighing the grain collected in the threshing tank during threshing simulations. The simulation results showed that section 1 and 2 threshed more grains compared to section 3 and 4. It was also found that when the threshing speed was higher, the larger number of grains were threshed. Only simulation was conducted in this study. Therefore, the validation of the simulation model is required. A comparative analysis to validate the simulation model by field experiment will be conducted in the future.

• Proceedings of the Korean Society of Computational and Applied Mathematics Conference
• /
• 2003.09a
• /
• pp.1-1
• /
• 2003

In this talk, a reduced-order modeling methodology based on centroidal Voronoi tessellations (CVT's)is introduced. CVT's are special Voronoi tessellations for which the generators of the Voronoi diagram are also the centers of mass (means) of the corresponding Voronoi cells. The discrete data sets, CVT's are closely related to the h-means clustering techniques. Even with the use of good mesh generators, discretization schemes, and solution algorithms, the computational simulation of complex, turbulent, or chaotic systems still remains a formidable endeavor. For example, typical finite element codes may require many thousands of degrees of freedom for the accurate simulation of fluid flows. The situation is even worse for optimization problems for which multiple solutions of the complex state system are usually required or in feedback control problems for which real-time solutions of the complex state system are needed. There hava been many studies devoted to the development, testing, and use of reduced-order models for complex systems such as unsteady fluid flows. The types of reduced-ordered models that we study are those attempt to determine accurate approximate solutions of a complex system using very few degrees of freedom. To do so, such models have to use basis functions that are in some way intimately connected to the problem being approximated. Once a very low-dimensional reduced basis has been determined, one can employ it to solve the complex system by applying, e.g., a Galerkin method. In general, reduced bases are globally supported so that the discrete systems are dense; however, if the reduced basis is of very low dimension, one does not care about the lack of sparsity in the discrete system. A discussion of reduced-ordering modeling for complex systems such as fluid flows is given to provide a context for the application of reduced-order bases. Then, detailed descriptions of CVT-based reduced-order bases and how they can be constructed of complex systems are given. Subsequently, some concrete incompressible flow examples are used to illustrate the construction and use of CVT-based reduced-order bases. The CVT-based reduced-order modeling methodology is shown to be effective for these examples and is also shown to be inexpensive to apply compared to other reduced-order methods.

• Journal of Power System Engineering
• /
• v.5 no.4
• /
• pp.53-60
• /
• 2001

This paper considers a modeling for the MIMO magnetic bearing system. The rotor is flexible and has a complex shape. To obtain the nominal plant transfer functions, we perform a numerical analysis by using the finite element method(F.E.M.) for the rotor's dynamics, and make a nominal model by reducing the modes from the results. And, we have experimented on the frequency response by a closed-loop identification method, and compared it with the simulation's result on the closed-loop control system.

• Transactions of Materials Processing
• /
• v.14 no.6 s.78
• /
• pp.547-552
• /
• 2005

In general, the sheet metal and die are described by finite elements for the simulation of the metal forming processes. Because the characteristics as continuum of the sheet metal are represented with triangles and rectangles, the errors occur inevitably in finite element analysis. Many contact schemes to describe the deformation modes exactly have been introduced in order to decrease these errors. In this study, a scheme for continuous contact treatment is proposed in order to consider the realistic behavior of contact phenomena during the forming process. The discrete mesh causes stepwise propagation of contact nodes of the sheet even though the contact region of the real forming process is altered very smoothly. It gives rise to convergence problem in case that the process, for example bending process, is sensitive to the contact between the sheet and the tools. The analysis of the unconstrained cylindrical bending process without blank holder is also presented in order to investigate the effect of the proposed algorithm.

• Elastomers and Composites
• /
• v.53 no.4
• /
• pp.207-212
• /
• 2018

Taylor rod impact tests have been the subject of many theoretical and experimental investigations. This paper discusses the numerical methods for simulating the Taylor impact test, which is widely used to obtain constitutive equations and failure conditions under high-velocity collisions of materials. With this in mind, a particle-based MPM (material point method) for linear viscoelastic solid materials was implemented, and MPM simulations for viscoelastic deformation behavior were numerically verified and confirmed by comparing the MPM and FEM results. In addition, this modeling and numerical approach could be extended to more complex viscoelastic models for basic understanding and to analyze the deformation and fracture behavior of more complicated viscoelastic material systems.

• Journal of the Korean Institute of Gas
• /
• v.7 no.4 s.21
• /
• 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 Disaster and Security
• /
• v.12 no.2
• /
• pp.57-63
• /
• 2019

The road sinks in the sewer line or subway section are affected by the ground characteristics. Therefore, it is necessary to accurately identify the relationship between the soil properties and the ground motion in the area where cavities occurred in order to establish a countermeasure against the road sink. In this paper, simulation was performed by using EDEM program, which is one of the Discrete Element Method programs, for sandy soil and clayey soil, which are most common in alluvial deposits, with different locations and sizes of cavities in the underground. As a result, it was found that the sink size occurred more in the sandy soil than in the cohesive soil. Deeper and larger cavity is more likely to occur the road sink In the sand soil model while road sink in the clay model is easy to occur when the cavity is more shallower.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.3
• /
• pp.177-183
• /
• 2018

This study simulated ice load and the motion response of a moored semi-submersible rig in pack-ice conditions using a finite element method. Ice flows of random size and shape were modeled, and interactions for ice-sea, ice-structure, ice-ice were simulated using a simplified method. Parameters for the simplified method such as drag force coefficient and the pressure-penetration relation were obtained based on the result of detailed analysis using the coupled Eulerian-Lagrangian method. The mooring lines were modeled by spring elements based on their stiffness. As a result of the simulation over 1,400 seconds, the force and motion response of the rig were obtained and validated using discrete elements and compared with the results found by the Krylov State Research Centre.