• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.33 no.3
• /
• pp.171-178
• /
• 2020

In this study, we develop MPI-OpenMP hybrid parallelization for multibody peridynamic simulations. Peridynamics is suitable for analyzing complicated dynamic fractures and various discontinuities. However, compared with a conventional finite element method, nonlocal interactions in peridynamics cost more time and memory. In multibody peridynamic analysis, the costs increase due to the additional interactions that occur when computing the nonlocal contact and ghost interlayer models between adjacent bodies. The costs become excessive when further refinement and smaller time steps are required in cases of high-velocity impact fracturing or similar instances. Thus, high computational efficiency and performance can be achieved by parallelization and optimization of multibody peridynamic simulations. The analytical code is developed using an Intel Fortran MPI compiler and OpenMP in NURION of the KISTI HPC center and parallelized through MPI-OpenMP hybrid parallelization. Further parallelization is conducted by hybridizing with OpenMP threads in each MPI process. We also try to minimize communication operations by model-based decomposition of MPI processes. The numerical results for the impact fracturing of multiple bodies show that the computing performance improves significantly with MPI-OpenMP hybrid parallelization.

• Journal of Biomedical Engineering Research
• /
• v.25 no.4
• /
• pp.295-300
• /
• 2004

Coil embolization technique has been used recently to treat cerebral aneurysms. When a giant or a multilobular aneurysm are treated by roils, filling an aneurysm sac completely with coils is difficult and partial blocking of an aneurysm sac is inevitable. Blood flow characteristics, which nay affect the embolization process of an aneurysm sac, are changed by the locations of coils for the Partially blocked aneurysms. Blood flow fields are also influenced by the geometry of a parent vessel. In order to suggest the coil locations effective for aneurysm embolization, the blood flow fields of lateral aneurysm models were analyzed for the different coil locations and parent vessel geometries. Three dimensional pulsatile flow fields are analyzed by numerical methods considering non-Newtonian viscosity characteristics of blood. Flow rate into the aneurysm sac (inflow rate) and wall shear stress, which are suspected as flow dynamic factors influencing aneurysm embolization, are also calculated. Inflow rates were smaller and the low wall shear stress zones were larger in the neck blocked models compared to the dome blocked models. Smaller inflow and larger low wall shear stress zones in the distal neck blocked model imply that the distal neck should be the effective coil locations for aneurysm embolization.

• Journal of Korean Society on Water Environment
• /
• v.32 no.3
• /
• pp.243-252
• /
• 2016

In cases of water pollution accidents, accurate prediction for arrival time and concentration of contaminants in a river is essential to take proper measures and minimize their impact on downstream water intake facilities. It is critical to fully understand the behavior characteristics of contaminants on river surface, especially in case of oil spill accidents. Therefore, in this study, the effects of main parameters of advection and diffusion of contaminants were analyzed and validated by comparing the results of Lagrangian particle tracking (LPT) simulation of Environmental Fluid Dynamic Code (EFDC) model with those of Global Position System (GPS)-equipped drifter experiment. Prevention scenario modeling was accomplished by taking cases of movable weir operation into account. The simulated water level and flow velocity fluctuations agreed well with observations. There was no significant difference in the speed of surface particle movement between 5 and 10 layer modeling. Therefore, 5 layer modeling could be chosen to reduce computational time. It was found that full three dimensional modeling simulated wind effects on surface particle movements more sensitively than depth-averaged two dimensional modeling. The diffusion range of particles was linearly proportional to horizontal diffusivity by sensitivity analysis. Horizontal diffusivity estimated from the results of GPS-equipped drifter experiment was 0.096 m²/sec, which was considered to be valid for applying the LPT module in this area. Finally, the scenario analysis results showed that particle movements could be stagnant when discharge from the upstream weir was reduced, implying the possibility of securing time for mitigation actions such as oil boom installation and wiping oil contaminants. The outcomes of this study can help improve the prediction accuracy of particle tracking simulation to establish the most suitable mitigation plan considering the combination of movable weir operation.

• The Korean Journal of Nuclear Medicine
• /
• v.37 no.5
• /
• pp.288-300
• /
• 2003

Purpose: To obtain regional blood flow and tissue-blood partition coefficient with time-activity curves from ${H_2}^{15}O$ PET, fitting of some parameters in the Kety model is conventionally accomplished by nonlinear least squares (NLS) analysis. However, NLS requires considerable compuation time then is impractical for pixel-by-pixel analysis to generate parametric images of these parameters. In this study, we investigated several fast parameter estimation methods for the parametric image generation and compared their statistical reliability and computational efficiency. Materials and Methods: These methods included linear least squres (LLS), linear weighted least squares (LWLS), linear generalized least squares (GLS), linear generalized weighted least squares (GWLS), weighted Integration (WI), and model-based clustering method (CAKS). ${H_2}^{15}O$ dynamic brain PET with Poisson noise component was simulated using numerical Zubal brain phantom. Error and bias in the estimation of rCBF and partition coefficient, and computation time in various noise environments was estimated and compared. In audition, parametric images from ${H_2}^{15}O$ dynamic brain PET data peformed on 16 healthy volunteers under various physiological conditions was compared to examine the utility of these methods for real human data. Results: These fast algorithms produced parametric images with similar image qualify and statistical reliability. When CAKS and LLS methods were used combinedly, computation time was significantly reduced and less than 30 seconds for $128{\times}128{\times}46$ images on Pentium III processor. Conclusion: Parametric images of rCBF and partition coefficient with good statistical properties can be generated with short computation time which is acceptable in clinical situation.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.23 no.4
• /
• pp.393-399
• /
• 2017

In a reciprocating compressor, the piston and connecting rod are important parts. Excess mechanical stress on these parts may cause damage, and broken parts are expensive and difficult to replace. Therefore, it is necessary to analyze the mechanical stress affecting durability and longevity. The main purpose of this study was to identify locations of maximum stress on pistons and connecting rods. Based on dynamic calculation of the working process of a specific air compressor, an analysis of piston and connecting rod performance has been completed. A three-dimensional model for the air compressor's pistons and connecting rods was built separately, and FEM analysis of these components was carried out using a numerical method. The pistons were loaded by pressure which was changed according to crankshaft angle without thermal boundary conditions. The simulation results were used to predict and estimate stress concentration as well as the value of this stress on pistons and connecting rods. The maximum equivalent stress calculated are over 190 MPa on pistons and 123 MPa on connecting rods at crank angle $135^{\circ}$ and $225^{\circ}$ but these are under tensile yield strength. Besides, the calculated safety factors of connecting rods and pistons is higher than 1. Moreover, the results obtained can be used to provide manufacturers with references to optimize the design of pistons and connecting rods for reciprocating compressors.

• Ocean Systems Engineering
• /
• v.1 no.4
• /
• pp.285-296
• /
• 2011

Spar platforms have several advantages for deploying wind turbines in offshore for depth beyond 120 m. The merit of spar platform is large range of topside payloads, favourable motions compared to other floating structures and minimum hull/deck interface. The main objective of this paper is to present the response analysis of the slack moored spar platform supporting 5MW wind turbine with bottom keel plates in regular and random waves, studied experimentally and numerically. A 1:100 scale model of the spar with sparD, sparCD and sparSD configuration was studied in the wave basin ($30{\times}30{\times}3m$) in Ocean engineering department in IIT Madras. In present study the effect of wind loading, blade dynamics and control, and tower elasticity are not considered. This paper presents the details of the studies carried out on a 16 m diameter and 100 m long spar buoy supporting a 90 m tall 5 MW wind turbine with 3600 kN weight of Nacelle and Rotor and 3500 kN weight of tower. The weight of the ballast and the draft of the spar are adjusted in such a way to keep the centre of gravity below the centre of buoyancy. The mooring lines are divided into four groups, each of which has four lines. The studies were carried out in regular and random waves. The operational significant wave height of 2.5 m and 10 s wave period and survival significant wave height of 6 m and 18 s wave period in 300 m water depth are considered. The wind speed corresponding to the operational wave height is about 22 knots and this wind speed is considered to be operating wind speed for turbines. The heave and surge accelerations at the top of spar platform were measured and are used for calculating the response. The geometric modeling of spar was carried out using Multisurf and this was directly exported to WAMIT for subsequent hydrodynamic and mooring system analysis. The numerical results were compared with experimental results and the comparison was found to be good. Parametric study was carried out to find out the effect of shape, size and spacing of keel plate and from the results obtained from present work ,it is recommended to use circular keel plate instead of square plate.

• Journal of the Korean Geotechnical Society
• /
• v.19 no.4
• /
• pp.167-178
• /
• 2003

Damage and overbreak of the remaining rock induced by blasting can not be avoided during tunnel construction which may result in either short-term or long-term tunnel instability. Therefore, in this paper, a methodology to take into account the effect of blast-induced damage in tunnel stability assessment is proposed. Dynamic numerical analysis was executed to evaluate damage and overbreak of the remaining rock for the most common blasting pattern in road tunnel. Rock damage was quantified by utilizing the damage variable factor which is adopted proposed in continuum damage mechanics. The damaged rock stiffness and the damaged failure criteria are used to consider the effect of rock damage in tunnel stability analysis. The damaged geological strength index of the damaged rock was newly proposed from the relationship between deformation modulus and geological strength index. Also the Hoek-Brown failure criteria of the damaged rock was obtained using the damaged geological strength index. Analysing the tunnel stability with the consideration of the blast-induced damage of remaining rock, it was found that the extend of plastic zone and deformation increased compared to the undamaged rock. Therefore the short-term or long-term tunnel stability will be threatened when the rock damage from blasting is ignored in the tunnel stability analysis.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.7 no.1
• /
• pp.75-90
• /
• 1995

A mathematical model for the dynamic analysis of a riser system with the inclusion of internal flow and nonlinear effects due to large structural displacements is developed to investigate the effect of internal flow on marine riser dynamics. The riser system accounts fir the nonlinear boundary conditions and includes a steady flow inside the pipe which is modeled as an extensible or inextensible. tubular beam subject to nonlinear three dimensional hydrodynamic loads such as current or wave excitation. Galerkin's finite element approximation and time incremental operator are implemented to derive the matrix equation of equilibrium for the finite element system and the extensibility or inextensibility condition is used to reduce degree of freedom of the system and the required computational time in the case of a nonlinear model. The algorithm is implemented to develop computer programs used in several numerical applications. The investigations of the effect of infernal flow on riser vibration due to current or wave loading are performed according to the change of various parameters such as top tension, internal flow velocity, current velocity, wave period, and so on. It is found that the effect of internal flow can be controlled by the increase of top tension. However, careful consideration has to be given in the design point particularly for the long riser under the harmonic loading such as waves. And it is also found that the consideration of nonlinear effects due to large structural displacements increases the effect of internal flow on riser dynamics.

• Journal of the Korea Organic Resources Recycling Association
• /
• v.18 no.3
• /
• pp.77-86
• /
• 2010

In this study, the model of the indirect wind suction waste sorting machine for characteristics of the screening of waste was studied using computational fluid dynamics and the drag coefficient for the model and the suction wind speed were obtained. The wind separator are developing by installing a cyclone air outlet to the suction blower impeller waste is selective in a way that does not pass the features and characteristics of the inlet pipe of the pressure loss and separation efficiency can have a significant impact on. Using Wind separator for selection of waste in the waste prior research on the aerodynamic properties are essential. For plastic cases, it is reasonable to take the drag coefficient between 0.8 and 1.0, and for cans, compression depending on whether the cans, the drag coefficient is in the range from 0.2 to 0.7. The separation efficiency of waste as change suction speed was the highest efficiency when the suction speed was 25~26 m/s. Shape of the inlet, depending on how the transfer pipe of the duct pressure loss occurs because the inlet velocity changes through the appropriate design standards to allow for continued research is needed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.5
• /
• pp.1907-1914
• /
• 2013

Textile-type heat exchangers installed on the tunnel walls for facilitating ground source heat pump systems, so called "energy textile", was installed in an abandoned railroad tunnel around Seocheon, South Korea. To evaluate thermal performance of the energy textile, a series of long-term monitoring was performed by artificially applying daily intermittent cooling and heating loads on the energy textile. In the course of the experimental measurement, the inlet and outlet fluid temperatures of the energy textile, pumping rate, temperature distribution in the ground, and air temperature inside the tunnel were continuously measured. From the long-term monitoring, the heat exchange rate was recorded as in the range of 57.6~143.5 W per one unit of the energy textile during heating operation and 362.3~558.4 W per one unit during cooling operation. In addition, the heat exchange rate of energy textile was highly sensitive to a change in air temperature inside the tunnel. The field measurements were verified by a 3D computational fluid dynamics analysis (FLUENT) with the consideration of air temperature variation inside the tunnel. The verified numerical model was used to evaluate parametrically the effect of drainage layer in the energy textile.