• Bulletin of the Korean Chemical Society
• /
• v.33 no.11
• /
• pp.3805-3809
• /
• 2012

Results for molecular dynamics simulation method of small liquid drops of argon (N = 1200-14400 molecules) at 94.4 K through a Lennard-Jones intermolecular potential are presented in this paper as a preliminary study of drop systems. We have calculated the density profiles ${\rho}(r)$, and from which the liquid and gas densities ${\rho}_l$ and ${\rho}_g$, the position of the Gibbs' dividing surface $R_o$, the thickness of the interface d, and the radius of equimolar surface $R_e$ can be obtained. Next we have calculated the normal and transverse pressure tensor ${\rho}_N(r)$ and ${\rho}_T(r)$ using Irving-Kirkwood method, and from which the liquid and gas pressures ${\rho}_l$ and ${\rho}_g$, the surface tension ${\gamma}_s$, the surface of tension $R_s$, and Tolman's length ${\delta}$ can be obtained. The variation of these properties with N is applied for the validity of Laplace's equation for the pressure change and Tolman's equation for the effect of curvature on surface tension through two routes, thermodynamic and mechanical.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.3 no.4
• /
• pp.263-273
• /
• 2011

The diffuse interface model of Saurel et al. (2008) is used for the computation of compressible cavitating flows around underwater missiles. Such systems use gas injection and natural cavitation to reduce drag effects. Consequently material interfaces appear separating liquid and gas. These interfaces may have a really complex dynamics such that only a few formulations are able to predict their evolution. Contrarily to front tracking or interface reconstruction method the interfaces are computed as diffused numerical zones, that are captured in a routinely manner, as is done usually with gas dynamics solvers for shocks and contact discontinuity. With the present approach, a single set of partial differential equations is solved everywhere, with a single numerical scheme. This leads to very efficient solvers. The algorithm derived in Saurel et al. (2009) is used to compute cavitation pockets around solid bodies. It is first validated against experiments done in cavitation tunnel at CNU. Then it is used to compute flows around high speed underwater systems (Shkval-like missile). Performance data are then computed showing method ability to predict forces acting on the system.

• Nuclear Engineering and Technology
• /
• v.54 no.8
• /
• pp.3117-3129
• /
• 2022

Deeply understanding the phase change of thin liquid sodium film inside wick pore is very important for further studying high-temperature sodium heat pipe's heat transfer. For the first time, the evaporation and condensation of thin liquid sodium film are investigated by the Lennard-Jones potential of molecular dynamics. Based on the startup and normal operation of the sodium heat pipe, three different cases are simulated. First, the equilibrium is achieved and the Mass Accommodation Coefficients of the three cases are 0.3886, 0.2119, 0.2615 respectively. Secondly, the non-equilibrium is built. The change of liquid film thickness, the number of gas atoms, the net evaporation flux (J_net), the heat transfer coefficient (h) at the liquid-gas interface are acquired. Results indicate that the magnitude of the J_net and the h increase with the basic equilibrium temperature. In 520-600 K (the startup of the heat pipe), the h has approached 5-6 W m^-2 K^-1 while liquid film thickness is in 11-13 nm. The fact shows that during the initial startup of the sodium heat pipe, the thermal resistance at the liquid-gas interface can't be negligible. This work is the complement and extension for macroscopic investigation of heat transfer inside sodium heat pipe. It can provide a reference for further numerical simulation and optimal design of the sodium heat pipe in the future.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.11
• /
• pp.1132-1138
• /
• 2008

A program that can simulate gas generator after burning liquid rocket engines was developed along with presenting the characteristics of the engines. The program was written in Matlab and used GUI interface so that many users can use it without any difficulties. The results of the program was compared with the real engine which was developed by the LRE advanced country. Most of the parameters concurred within 1% error expect for the pressure at the turbopump. The reasons of the large differences were supposed that pressure decreases at the schematics were smaller than that of the real engines.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.658-659
• /
• 2008

Numerical simulation is performed for a microdroplet deposition on the pre-patterned micro-structure. The liquid-air interface is tracked by level set method improved by incorporating the ghost fluid approach based on a sharp-interface representation. The method is further extended to treat the contact angle condition at an immersed solid surface. The present computation of a patterning process using microdroplet ejection demonstrates that the multiphase characteristics between the liquid-gas-solid phases can be used to overcome the patterning error.

• Journal of Korea Foundry Society
• /
• v.33 no.4
• /
• pp.157-161
• /
• 2013

A ribbon-type polycrystalline silicon wafer was directly fabricated from liquid silicon via a novel technique for both a fast growth rate and large grain size by exploiting gas pressure. Effects of processing parameters such as moving speed of a dummy bar and the length of the solidification zone on continuous casting of the silicon wafer were investigated. Silicon melt extruded from the growth region in the case of a solidification zone with a length of 1cm due to incomplete solidification. In case of a solidification zone wieh a length of 2 cm, on the other hand, continuous casting of the wafer was impossible due to the volume expansion of silicon derived from the liquid-solid transformation in solidification zone. Consequently, the optimal length of the solidification zone was 1.5 cm for maintaining the position of the solid-liquid interface in the solidification zone. The silicon wafer could be continuously casted when the moving speed of the dummy bar was 6 cm/min, but liquid silicon extruded from the growth region without solidification when the moving speed of the dummy bar was ${\geq}$ 9 cm/min. This was due to a shift of the position of the solid-liquid interface from the solidification zone to the moving area. The present study reports experimental findings on a new direct growth system for obtaining silicon wafers with both high quality and productivity, as a candidate for an alternate route for the fabrication of ribbon-type silicon wafers.

• Nuclear Engineering and Technology
• /
• v.53 no.6
• /
• pp.1821-1833
• /
• 2021

The gas-liquid counter-current flow limitation (CCFL) is closely related to efficient and safety operation of many equipment in industrial cycle. Air-water countercurrent flow experiments were performed in a tube with diameter of 25 mm to understand the triggering mechanism of CCFL. A parallel electrode probe was utilized to measure film thickness whereby the time domain and frequency domain characteristics of liquid film was obtained. The amplitude of the interface wave is small at low liquid flow rate while it becomes large at high liquid flow rate after being disturbed by the airflow. The spectral characteristic curve shows a peak-shaped distribution. The crest exists between 0 and 10 Hz and the amplitude decreases with the frequency increase. The analysis of visual observation and characteristic of film thickness indicate that two flooding mechanisms were identified at low and high liquid flow rate, respectively. At low liquid flow rate, the interfacial waves upward propagation is responsible for the formation of CCFL onset. While flooding at high liquid flow rate takes place as a direct consequence of the liquid bridging in tube due to the turbulent flow pattern. Moreover, it is believed that there is a transition region between the low and high liquid flow rate.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.12
• /
• pp.3360-3371
• /
• 1995

Hydrodynamic characteristics of a stratified flow at a Y-junction were examined analytically. Gas is supplied through the main horizontal channel and liquid is introduced into the gas stream from the Y-shaped bottom branch. Analysis was performed with irrational flow and inviscid fluid assumptions. The Stokes' inverse transformation technique was adopted to convert the real x-y plane into the x-.psi. plane. The potential flow equation was solved numerically in the transformed (x-.psi.) plane and the interface profile, pressure distribution and the streamlines were obtained. The effects of the inlet conditions, injection angle and the gravity on the flow characteristics were also examined. To check the validity of the present method, the previous resultant the two-dimensional obtuse wedge flow was compared. The inverse transformation technique turned out to be also very useful to predict the hydrodynamic characteristics of a stratified flow with the pressure variation at a Y-shaped mixing junction.

• Nuclear Engineering and Technology
• /
• v.53 no.7
• /
• pp.2221-2228
• /
• 2021

One-dimensional oxygen transport relation is indispensable to study the oxygen distribution in the LBE-cooled system with an oxygen control device. In this paper, a numerical research is carried out to study the oxygen transport characteristics in a gas-phase oxygen control device, including the static case and dynamic case. The model of static oxygen control is based on the two-phase VOF model and the results agree well with the theoretical expectation. The model of dynamic oxygen control is simplified and the gas-liquid interface is treated as a free surface boundary with a constant oxygen concentration. The influences of the inlet and interface oxygen concentration, mass flow rate, temperature, and the inlet pipe location on the mass transfer characteristics are discussed. Based on the results, an oxygen mass transport relation considering the temperature dependence and velocity dependence separately is obtained. The relation can be used in a one-dimensional system analysis code to predict the oxygen provided by the oxygen control device, which is an important part of the integral oxygen mass transfer models.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.2123-2128
• /
• 2004

In this study, analysis is made for the effects of groove shape on the thermal performance of a axial groove heat pipe. The mathematical models of two-phase flow in grooved heat pipe are presented for the capillary limitation in steady state. Generally, the heat pipe performance depends on the capillary pressure and liquid flow. The friction force of liquid flow through the groove increases with the groove width decreased, and then the capillary pressure is improved in the gas-liquid interface of groove. Therefore, the optimal groove width shaper exists for the maximum thermal performance of heat pipe. In this paper, the optimal groove shape and scale are presented by considering both capillary pressure and liquid flow.