• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2002.05a
• /
• pp.181-183
• /
• 2002

The accurate permeability for preform is critical to model and design the impregnation of fluid resin in the composite manufacturing process. In this study, the in-plane and transverse permeability for a woven fabric are predicted numerically through the coupled flow model which combines microscopic with macroscopic flow. The microscopic and macroscopic flow which are flows within the micro-unit and macro-unit cell, respectively, are calculated by using 3-D CVFEM(control volume finite element method). To avoid checker-board pressure field and improve the efficiency on numerical computation, A new interpolation function for velocity is proposed on the basis of analytic solutions. The permeability of plain woven fabric is measured through unidirectional flow experiment and compared with the permeability calculated numerically. Based on the good agreement of the results, the relationships between the permeability and the structures of preform such as the fiber volume fraction and stacking effect can be understood. The reverse and the simple stacking are taken in account. Unlike past literatures, this study is based on more realistic unit cell and the improved prediction of permeability can be achieved. It is observed that in-plane flow is more dominant than transverse flow in the real flow through preform and the stacking effect of multi-layered preform is negligible. Consequently, the proposed coupled flow model can be applied to modeling of real composite materials processing.

• Journal of Mechanical Science and Technology
• /
• v.16 no.7
• /
• pp.966-974
• /
• 2002

An analytic approach has been employed to study condensate film thickness distribution inside cave-shaped cavity of a flat plate heat pipe. The results indicate that the condensate film thickness largely depends on mass flow rate and local velocity of condensate. The increasing rate of condensate film for circular region reveals about 50% higher value than that of vertical region. The physical properties of working fluid affect significantly the condensate film thickness, such as the condensate film thickness for the case of FC-40 are 5 times larger than that of water. In comparison with condensation on a vertical wall, the average heat transfer coefficient in the cave-shaped cavity presented 10∼15% lower values due to the fact that the average film thickness formed inside the cave-shaped cavity was larger than that of the vertical wall with an equivalent flow length. A correlation formula which is based on the condensate film analysis for the cave-shaped cavity to predict average heat transfer coefficient is presented. Also, the critical minimum fill charge ratio of working fluid based on condensate film analysis has been predicted, and the minimum fill charge ratios for FC-40 and water are about Ψ$\_$crit/=3∼7%, Ψ$\_$crit/=0.5∼1.3% respectively, in the range of heat flux q"=5∼90kW/㎡.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.1
• /
• pp.56-65
• /
• 2018

Pressure oscillation caused by the compressibility of entrapped air in dam break flow is analyzed using an open source code, which is a two-phase compressible code for non-isothermal immiscible fluids. Since compressible flows are computed based on a pressure-based method, the code can handle the equation of state of barotropic fluid, which is virtually incompressible. The computed time variation of pressure is compared with other experimental and computational results. The present result shows good agreements with other results until the air is entrapped. As the entrapped air bubbles pulsate, pressure oscillations are predicted and the pressure oscillations damp out quickly. Although the compressibility parameter of water has been varied for a wide range, it has no effects on the computed results, because the present equation of state for water is so close to that of incompressible fluid. Grid independency test for computed time variation of pressure shows that all results predict similar period of pressure oscillation and quick damping out of the oscillation, even though the amplitude of pressure oscillation is sensitive to the velocity field at the moment of the entrapping. It is observed that as pressure inside the entrapped air changes quickly, the pressure field in the neighboring water adjusts instantly, because the sound of speed is much higher in water. It is confirmed that the period of pressure oscillation is dominated by the added mass of neighboring water. It is found that the temperature oscillation of the entrapped air is critical to the quick damping out of the oscillations, due to the fact that the time averaged temperature inside the entrapped air is higher than that of surrounding water, which is almost constant.

• Tribology and Lubricants
• /
• v.37 no.6
• /
• pp.253-260
• /
• 2021

Metal pipes are used in a wide range of applications, from plumbing systems of large construction sites to small devices such as medical tools. When a liquid is enforced to flow through a metal pipe, a higher flow rate is beneficial for higher efficiency. Using high pressures can enhance the flow rate yet can be harmful for medical applications. Thus, we consider an optimal geometrical design to increase the flow rate in medical devices. In this study, we focus on cannulas, which are widely used small metal pipes for surgical procedures, such as liposuction. We characterize the internal flow rate driven by a negative pressure and explore its dependence on the key design parameters. We quantitatively analyze the suction characteristics for each design variable by conducting computational fluid dynamics simulations. In addition, we build a suction performance measurement system which enables the translational motion of cannulas with pre-programmed velocity for experimental validation. The inner diameter, section geometry, and hole configuration are the design factors to be evaluated. The effect of the inner diameter dominates over that of section geometry and hole configuration. In addition, the circular tube shape provides the maximum flow rate among the elliptical geometries. Once the flow rate exceeds a critical value, the rate becomes independent of the number and width of the suction holes. Finally, we introduce a slippery liquid-infused nanoporous surface (SLIPS) coating using nanoparticles and hydrophobic lubricants that effectively improves the flow rate and antifouling property of cannulas without altering the geometrical design parameter.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.187-190
• /
• 2006

In the present study, effects of tree-stream turbulence and surface trip wire on the flow past a sphere at $Re\;=\;0.4\;{\times}\;10^5\;{\sim}\;2.8\;{\times}\;10^5$ are investigated through wind tunnel experiments. Various types of grids are installed upstream of the sphere in order to change the tree-stream turbulence intensity. In the case of surface trip wire, 0.5mm and 2mm trip wires are attached from $20^{\circ}\;{\sim}\;90^{\circ}$ at $10^{\circ}$ interval along the streamwise direction. To investigate the flow around a sphere, drag measurement using a load cell, surface-pressure measurement, surface visualization using oil-flow pattern and near-wall velocity measurement using an I-type hot-wire probe are conducted. In the variation of free-stream turbulence, the critical Reynolds number decreases and drag crisis occurs earlier with increasing turbulence intensity. With increasing Reynolds number, the laminar separation point moves downstream, but the reattachment point after laminar separation and the main separation point are fixed, resulting in constant drag coefficient at each free-stream turbulence intensity. At the supercritical regime, as Reynolds number is further increased, the separation bubble is regressed but the reattachment and the main separation points are fixed. In the case of surface trip wire directly disturbing the boundary layer flow, the critical Reynolds number decreases further with trip wire located more downstream. However, the drag coefficient after drag crisis remains constant irrespective of the trip location.

• The KIPS Transactions:PartA
• /
• v.14A no.4
• /
• pp.191-196
• /
• 2007

In the field of computer graphics, Navier-Stokes equations would be used for realistic simulations of smokes and currents. However, implementations derived from these equations are hard to achieve for real-time simulations, mainly due to its massive and complex calculations. Thus, there have been various attempts to approximate these equations for real-time simulation of smokes and others. When the advection terms of the equations are approximated by the Semi-Lagrange methods, the fluid density can be rapidly reduced and small-scale vorticity phenomena are easy to be missed, mainly due to the numerical losses over time. In this paper, we propose an improved numerical method to approximately calculate the advection terms, and thus eliminate these problems. To calculate the advection terms, our method starts to set critical regions around the target grid points. Then, among the grid points in a specific critical region, we search for a grid point which will be advected to the target grid point, and use the velocity of this grid point as its advection vector. This method would reduce the numerical losses in the calculation of densities and vorticity phenomena, and finally can implement more realistic smoke simulations. We also improve the overall efficiency of vector calculations and related operations through GPU-based implementation techniques, and thus finally achieve the real-time simulation.

• Nuclear Engineering and Technology
• /
• v.42 no.4
• /
• pp.382-393
• /
• 2010

Thermal mixing by steam jets in a pool is dominantly influenced by a turbulent water jet generated by the condensing steam jets, and the proper prediction of this turbulent jet behavior is critical for the pool mixing analysis. A turbulent jet flow induced by a steam jet discharged through a vertical upward single hole into a subcooled water pool was subjected to computational fluid dynamics (CFD) analysis. Based on the small-scale test data derived under a horizontal steam discharging condition, this analysis was performed to validate a CFD method of analysis previously developed for condensing jet-induced pool mixing phenomena. In previous validation work, the CFD results and the test data for a limited range of radial and axial directions were compared in terms of profiles of the turbulent jet velocity and temperature. Furthermore, the behavior of the turbulent jet induced by the steam jet through a horizontal single hole in a subcooled water pool failed to show the exact axisymmetric flow pattern with regards to an overall pool mixing, whereas the CFD analysis was done with an axisymmetric grid model. Therefore, another new small-scale test was conducted under a vertical upward steam discharging condition. The purpose of this test was to generate the velocity and temperature profiles of the turbulent jet by expanding the measurement ranges from the jet center to a location at about 5% of $U_m$ and 10 cm to 30 cm from the exit of the discharge nozzle. The results of the new CFD analysis show that the recommended CFD model of the high turbulent intensity of 40% for the turbulent jet and the fine mesh grid model can accurately predict the test results within an error rate of about 10%. In this work, the turbulent jet model, which is used to simply predict the temperature and velocity profiles along the axial and radial directions by means of the empirical correlations and Tollmien's theory was improved on the basis of the new test data. The results validate the CFD model of analysis. Furthermore, the turbulent jet model developed in this study can be used to analyze pool thermal mixing when an ellipsoidal steam jet is discharged under a high steam mass flux in a subcooled water pool.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.11
• /
• pp.895-903
• /
• 2007

CMBT(Curved Moving Boundary Treatment) is a newly developed scheme for the treatment of a no slip condition on the curved solid wall of moving obstacle in a flow field. In our research CMBT was used to perform LBM simulation of a flow over a moving circular cylinder to determine the flow feature and aerodynamics characteristic of the cylinder. To ascertain the applicability of CMBT on the complex shape of the obstacle, it was first simulated for the case of the flow over a fixed circular cylinder in a channel and the results were compared against the solution of Navier-Stokes equation with deforming mesh technique. The simulations were performed in a moderate range of reynolds number at each moving cylinder to identify the flow feature and aerodynamic characteristics of circular cylinder in a channel. The drag coefficients of the cylinder were calculated from the simulation results. We have numerically confirmed that the critical reynolds number for vortex shedding is ar Re=250 and the result is the same as the case of fixed cylinder. As the cylinder approaching to one wall, the 2nd vortex is developed by interacting with the wall boundary-layer vorticity. As the velocity ratio increase the third vortex are generated by interacting with the 2nd vortexes developed on the upper and lower wall boundary layer. The resultant $C_d$ decrease as reynolds number increasing and the Cd approached to a value when Re>1000.

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

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