• Journal of the Korean Society of Marine Environment & Safety
• /
• v.24 no.6
• /
• pp.796-802
• /
• 2018

The effects of inner liquid sloshing on vessel motions are a well-known factor. It was investigated experimentally and numerically. In this regard, the study of many efforts to reduce natural phenomena of vessel motions by adopting special devices especially for roll motions. Among many devices, inserting baffles in the inner liquid tank is very common. In this study, one investigated the vessel motions with inner sloshing tanks with baffles inside. For the numerical simulation, one employed a dynamically coupled program between boundary-element-method-based vessel motion analysis program and a particle-based computational fluid dynamics program. Comparing corresponding experimental results validated the dynamically coupled program. The validated coupled program was used to simulate vessel motions, including sloshing effects with various lengths of inner baffles. The simulation results show that not only the filling ratio of inner liquid, but also the length of clearance due to baffles influenced the vessel motions. The significant point of this study was that the natural frequency of vessel motions can be maintained irrespective of the amount of filling ratio through adjustment of the clearance. In a future study, the effects of various numbers of baffles with various clearances would be conducted to percuss the possibility of vessel motion control with inner liquid sloshing effects.

• Journal of the Institute of Electronics Engineers of Korea CI
• /
• v.48 no.1
• /
• pp.90-100
• /
• 2011

We present a parallel bi-conjugate gradient (Bi-CG) matrix solver for large scale Bio-FET simulations based on recent graphics processing units (GPUs) which can realize a large-scale parallel processing with very low cost. The proposed method is focused on solving the Poisson equation in a parallel way, which requires massive computational resources in not only semiconductor simulation, but also other various fields including computational fluid dynamics and heat transfer simulations. As a result, our solver is around 30 times faster than those with traditional methods based on single core CPU systems in solving the Possion equation in a 3D FDM (Finite Difference Method) scheme. The proposed method is implemented and tested based on NVIDIA's CUDA (Compute Unified Device Architecture) environment which enables general purpose parallel processing in GPUs. Unlike other similar GPU-based approaches which apply usually 32-bit single-precision floating point arithmetics, we use 64-bit double-precision operations for better convergence. Applications on the CUDA platform are rather easy to implement but very hard to get optimized performances. In this regard, we also discuss the optimization strategy of the proposed method.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.9
• /
• pp.69-76
• /
• 2018

The present study analyzed the effect of film hole position of 45 degree ribbed cooling channel on film cooling performance of gas turbine blades. We also investigated the influence of the ribs under the fixed blowing ratio. Three-dimensional numerical model was constructed and extensive simulation was conducted using the commercial code (Fluent ver. 17.0) under steady-state condition. Base on the simulation results, We investigated the cooling effectiveness, flow velocity, streamline, and pressure coefficient. Moreover, We analyzed the effect of cooling hole position on ejection of the secondary flow caused by the rib structure. From the results, It was found that internal flow of the cooling channel forms a vortex pair in the counterclockwise from the top side, and clockwise from the bottom side. For the channels with ribs, the vortex flow generated by the ribs caused a higher pressure difference near the hole outlet, resulting in at least 12% higher cooling effectiveness than the channel without ribs. Additionally, when the hole is located on the left side of the ribbed channel (Rib-Left), it can be found that the secondary flow generated by the ribs hits against wall surface near the hole to form a flow in the direction of the hole inclination angle. Therefore, It is considered that the region where the cooling gas discharged to the blade surface stays in the main flow boundary layer is wider than the other cases. In this case, The largest pressure coefficient difference was observed near the outlet of the hole, and as a result, the discharge of the cooling gas was accelerated and the cooling efficiency was slightly increased.

• Journal of the Korean Vacuum Society
• /
• v.22 no.2
• /
• pp.66-78
• /
• 2013

In order to develop of 80 kW RF plasma torch system, we achieved three-dimensional simulations for the extraction of more information as temperature in torch and fluid behavior analysis, etc. The position of powder injection tube, the plasma discharge characteristics with various input current and various length of ceramic tube, and the plasma temperature characteristics with process gas flow rate such those was simulated. RF thermal plasma torch designed by simulation was manufactured that was measured to the maximum of 89.3 kW power. The mass production using developed 80 kW RF thermal plasma torch system were investigated by characteristics manufactured of Si nano powder. The mass-production level of Si nano-powder was average of 539 g/hr and high yield rate of 71.6%, respectively. The particle size distribution $D_{99}/D_{50}$ of manufacturing nano-powder was investigated to 1.98 as a good uniform.

• Journal of Bio-Environment Control
• /
• v.29 no.3
• /
• pp.285-292
• /
• 2020

Because the inner environment of greenhouse has a direct impact on crop production, many studies have been performed to develop technologies for controlling the environment in the greenhouse. However, it is difficult to apply the technology developed to all greenhouses because those studies were conducted through empirical experiments in specific greenhouses. It takes a lot of time and cost to develop the models that can be applicable to all greenhouse in real situation. Therefore studies are underway to solve this problem using computer-based simulation techniques. In this study, a model was developed to predict the inner environment of glass greenhouse using CFD simulation method. The developed model was validated using primary and secondary heating experiment and daytime greenhouse inner temperature data. As a result of comparing the measured and predicted value, the mean temperature and uniformity were 2.62℃ and 2.92%p higher in the predicted value, respectively. R² was 0.9628, confirming that the measured and the predicted values showed similar tendency. In the future, the model needs to improve by applying the shape of the greenhouse and the position of the inner heat exchanger for efficient thermal energy management of the greenhouse.

• Transactions of the KSME C: Technology and Education
• /
• v.3 no.1
• /
• pp.23-28
• /
• 2015

In this paper, we performed the design optimization dual-shell and tube heat exchanger on exhaust waste heat recovery for gas heat pump using CFD and RSM. CFD analysis is useful to design the complex structure such as double shell and tube heat exchanger. By computer simulation, engineers can assess the feasibility of the given design factors and change them to get a better design. But if one wishes to perform complex analysis on the simulation, such dual-shell and tube heat exchanger for GHP, the computational time can become overwhelming. CFD is powerful but it takes a lot of time for complex structure. Therefore, the CFD analysis is minimized by the optimization using the RSM method. As a result, the number of baffle and tube are optimized by 6 baffles and 25 tubes for heat transfer and flow friction. And then pressure drop and heat transfer is improved about 12.2%. We confirm the design optimization using CFD and RSM is useful on complex structure of heat exchanger.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.216-227
• /
• 2004

A new cavitating model by using bubble size distribution based on bubbles-mass has been proposed. Both liquid and vapor phases are treated with Eulerian framework as a mixture containing minute cavitating bubbles. In addition vapor phase consists of various sizes of vapor bubbles, which are distributed to classes based on their mass. The bubble number-density for each class was solved by considering the change of the bubble-mass due to phase change as well as generation of new bubbles due to heterogeneous nucleation. In this method, the bubble-mass is treated as an independent variable, and the other dependent variables are solved in spatial coordinates and bubble-mass coordinate. Firstly, we employed this method to calculate bubble nucleation and growth in stationary super-heated liquid nitrogen, and bubble collapse in stationary sub-cooled one. In the case of bubble growth in super-heated liquid, bubble number-density of the smallest class based on its mass is increased due to the nucleation. These new bubbles grow with time, and the bubbles shift to larger class. Therefore void fraction of each class is increased due to the growth in the whole class. On the other hand, in the case of bubble collapse in sub-cooled liquid, the existing bubbles are contracted, and then they shift to smaller class. It finally becomes extinct at the smallest one. Secondly, the present method is applied to a cavitating flow around NACA00l5 foil. Liquid nitrogen and liquid oxygen are employed as working fluids. Cavitation number, $\sigma$, is fixed at 0.15, inlet velocities are changed at 5, 10, 20 and 50m/s. Inlet temperatures are 90K in case of liquid nitrogen, and 90K and 1l0K in case of liquid oxygen. 110K of oxygen is corresponding to the 90K of nitrogen because of the same relative temperature to the critical one, $T_{r}$=$T/T_c^{+}$. Cavitating flow around the NACA0015 foils was properly analyzed by using bubble size distribution. Finally, the method is applied to a cavitating flow in an inducer of the LE-7A hydrogen turbo-pump. This inducer has 3 spiral foils. However, for simplicity, 2D calculation was carried out in an unrolled channel at 0.9R cross-section. The channel moves against the fluid at a peripheral velocity corresponding to the inducer revolutions. Total inlet pressure, $Pt_{in}$, is set at l00KPa, because cavitation is not generated at a design point, $Pt_{in}$=260KPa. The bubbles occur upstream of the foils and collapse between them. Cavitating flow in the inducer was successfully predicted by using the bubble size distribution.

• Plant Journal
• /
• v.9 no.3
• /
• pp.38-42
• /
• 2013

Submarine gas fields have focused because of the increasing fuel cost, the environmental regulations, and the safety & NIMBY problems. LNG-FPSO which is available for acid gas removal, recovery of the condensate & LPG and Liquefaction in topside process is one of high technology offshore structures. On the other hands, it is necessary to verify the pre-treatment efficiency by the ship motion and to apply to the design for LNG-FPSO. This study is to develop the pre-treatment technology for LNG-FPSO as taking account to the process efficiency by ship motion effects and the area optimization. Based on the simulation results, it founds that hybrid process shows the low circulate rate, the low heat duty and the small size of column dimensions compared to typical amine process. It will be verified the process efficiency in the various conditions by sea states as performing the 6-DOF motion test and CFD simulation.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.6
• /
• pp.532-539
• /
• 2008

Kwon et al.(2008) carried out a hydraulic experiment in order to evaluate Manning coefficient, which implicates flow resistance due to bottom friction as well as drag caused by the squared piers higher than water depth and arranged with equal intervals, under the flow condition with a constant drag coefficient, $Re>10^4$. And, based on the equation of equilibrium, they proposed a formula for the equivalent resistant coefficient including empirical drag interaction coefficient obtained by using the experimental results. In this study, the hydraulic experiment was simulated using FLOW-3D, a 3-dimensional computational fluid dynamic code. The computations were compared with the experiment results as well as the semi-theoretical formula, and the comparisons show a good agreement. From the agreement, it was confirmed that when flow resistance bodies were higher than water depth, Manning n value increases with 2/3 power of water depth as shown in the theoretical formula and that drag interaction coefficient was dominated by their intervals.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.2
• /
• pp.92-108
• /
• 2017

Present work shows a numerical method to analysis of interaction analysis between solitary wave and onshore bridge. Numerical simulation is carried out by TWOPM-3D (three-dimensional one-field model for immiscible two-phase flows), which is based on Navier-Stokes solver. To do this, the solitary wave is generated numerically in numerical wave channel, and numerical results and experimental results were compared and analyzed in order to verify the applicability of force acting on an onshore bridge. From this, we discussed precisely the characteristics of horizontal and vertical forces (uplift and downward forces) changes including water level and velocity changes due to the variation of solitary wave height, water depth, onshore bridge's location and type, and number of girder. Furthermore, It is revealed that the maximum horizontal and vertical forces acting on the girder bridge show different varying properties according to the number of girder, although each maximum force acting on the girder bridge is proportional to the increasement of incident solitary wave height, and the entrained air in the fluid flow affects the vertical force highly.