• Nuclear Engineering and Technology
• /
• v.55 no.2
• /
• pp.452-459
• /
• 2023

The core power control is an important issue for the study of dynamic characteristics in China initiative accelerator driven subcritical system (CiADS), which has direct impact on the control strategy and safety analysis process. The CiADS is an experimental facility that is only controlled by the proton beam intensity without considering the control rods in the current engineering design stage. In order to get the optimized operation scheme with the stable and reliable features, the variation of beam intensity using the continuous and periodic control approaches has been adopted, and the change of collimator and the adjusting of duty ratio have been proposed in the power control process. Considering the neutronics and the thermal-hydraulics characteristics in CiADS, the physical model for the core power control has been established by means of the point reactor kinetics method and the lumped parameter method. Moreover, the multi-inputs single-output (MISO) logical structure for the power control process has been constructed using proportional integral derivative (PID) controller, and the meta-heuristic algorithm has been employed to obtain the global optimized parameters for the stable running mode without producing large perturbations. Finally, the verification and validation of the control method have been tested based on the reference scenarios in considering the disturbances of spallation neutron source and inlet temperature respectively, where all the numerical results reveal that the optimization method has satisfactory performance in the CiADS core power control scenarios.

• Steel and Composite Structures
• /
• v.50 no.5
• /
• pp.563-581
• /
• 2024

This paper presents a novel finite element model for the free vibration analysis of variable-thickness functionally graded porous (FGP) microplates resting on Pasternak's medium in the hygro-thermal environment. The governing equations are established according to refined higher-order shear deformation plate theory (RPT) in construction with the modified couple stress theory. For the first time, three-node triangular elements with twelve degrees of freedom for each node are developed based on Hermitian interpolation functions to describe the in-plane displacements and transverse displacements of microplates. Two laws of variable thickness of FGP microplates, including the linear law and the nonlinear law in the x-direction are investigated. Effects of thermal and moisture changes on microplates are assumed to vary continuously from the bottom surface to the top surface and only cause tension loads in the plane, which does not change the material's mechanical properties. The numerical results of this work are compared with those of published data to verify the accuracy and reliability of the proposed method. In addition, the parameter study is conducted to explore the effects of geometrical and material properties such as the changing law of the thickness, length-scale parameter, and the parameters of the porosity, temperature, and humidity on the free vibration response of variable thickness FGP microplates. These results can be applied to design of microelectromechanical structures in practice.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.5 no.1
• /
• pp.21-30
• /
• 1985

The tension leg platform (TLP) is a kind of compliant structures, and is also a type of moored stable platform with a buoyancy exceeding the weight because of having tensioned vertical anchor cables. In this paper, among the various kinds of tension leg structures, Deep Oil Technology (DOT) TLP was analyzed because it has large-displacement portions of the immersed surface such as vertical corner pontoons and small-diameter elongated members such as cross-bracing. It also has results of hydraulic model tests, comparable with theorectical analysis. Because of the vertical axes of symmetry in the three vertical buoyant legs and because there are no larger horizontal buoyant members between these three vertical members, it was decided to develop a numerical algorithm which would predict the dynamic response of the DOT TLP using the previously developed numerical algorithm Floating Vessel Response Simulation (FVRS) for vertically axisymmetric bodies of revolution. In addition, a linearized hydroelastic Morison equation subroutine would be developed to account for the hydrodynamic pressure forces on the small member cross bracing. Interaction between the large buoyant members or small member cross bracings is considered to be negligible and is not included in the analysis. The dynamic response of the DOT TLP in the surge mode is compared with the results of the TLP algorithm for various combinations of diffraction and Morison forces and moments. The results which include the Morison equation are better than the results for diffraction only. This is because the vertically axisymmetric buoyant members are only marginally large enough to consider diffractions effects. The prototype TLP results are expected to be more inertially dominated.

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