Numerical Simulation of Two-Phase Sloshing Flow in LNG Tank Using Finite-Analytic Level-Set Method

Author(s):  
Kai Yu ◽  
Hamn-Ching Chen ◽  
Jang Whan Kim ◽  
Young-Bum Lee

Impact pressure due to sloshing is of great concern for the ship owners, designers and builders of the LNG carriers regarding the safety of LNG containment system and hull structure. Sloshing of LNG in a partially filled tank has been an active area of research with numerous experimental and numerical investigations over the past decade. In order to accurately predict the sloshing impact load, it is necessary to develop advanced numerical simulation tools which can provide accurate resolution of local flow phenomena including wave breaking, jet formation, gas entrapping and liquid-gas interactions. In the present study, a new numerical method is developed for the simulation of violent sloshing flow inside a three-dimensional LNG tank considering wave breaking and liquid-gas interaction. The sloshing flow inside a membrane-type LNG tank is simulated numerically using the Finite-Analytic Navier-Stokes (FANS) method. The governing equations for two-phase air and water flows are formulated in curvilinear coordinate system and discretized using the finite-analytic method on a non-staggered grid. Simulations were performed for LNG tank in transverse and longitudinal motions including horizontal, vertical, and rotational motions. The predicted impact pressures were compared with the corresponding experimental data. The validation results clearly illustrate the capability of the present two-phase FANS method for accurate prediction of impact pressure in sloshing LNG tank including violent free surface motion, three-dimensional instability and air trapping effects.

2017 ◽  
Vol 2017 ◽  
pp. 1-9 ◽  
Author(s):  
Zhengzhi Wang ◽  
Chunling Zhu

In view of the rotor icing problems, the influence of centrifugal force on rotor blade icing is investigated. A numerical simulation method of three-dimensional rotor blade icing is presented. Body-fitted grids around the rotor blade are generated using overlapping grid technology and rotor flow field characteristics are obtained by solving N-S equations. According to Eulerian two-phase flow, the droplet trajectories are calculated and droplet impingement characteristics are obtained. The mass and energy conservation equations of ice accretion model are established and a new calculation method of runback water mass based on shear stress and centrifugal force is proposed to simulate water flow and ice shape. The calculation results are compared with available experimental results in order to verify the correctness of the numerical simulation method. The influence of centrifugal force on rotor icing is calculated. The results show that the flow direction and distribution of liquid water on rotor surfaces change under the action of centrifugal force, which lead to the increasing of icing at the stagnation point and the decreasing of icing on both frozen limitations.


2018 ◽  
Vol 857 ◽  
pp. 270-290 ◽  
Author(s):  
Josef Hasslberger ◽  
Markus Klein ◽  
Nilanjan Chakraborty

This paper presents a detailed investigation of flow topologies in bubble-induced two-phase turbulence. Two freely moving and deforming air bubbles that have been suspended in liquid water under counterflow conditions have been considered for this analysis. The direct numerical simulation data considered here are based on the one-fluid formulation of the two-phase flow governing equations. To study the development of coherent structures, a local flow topology analysis is performed. Using the invariants of the velocity gradient tensor, all possible small-scale flow structures can be categorized into two nodal and two focal topologies for incompressible turbulent flows. The volume fraction of focal topologies in the gaseous phase is consistently higher than in the surrounding liquid phase. This observation has been argued to be linked to a strong vorticity production at the regions of simultaneous high fluid velocity and high interface curvature. Depending on the regime (steady/laminar or unsteady/turbulent), additional effects related to the density and viscosity jump at the interface influence the behaviour. The analysis also points to a specific term of the vorticity transport equation as being responsible for the induction of vortical motion at the interface. Besides the known mechanisms, this term, related to surface tension and gradients of interface curvature, represents another potential source of turbulence production that lends itself to further investigation.


2014 ◽  
Vol 541-542 ◽  
pp. 1288-1291
Author(s):  
Zhi Feng Dong ◽  
Quan Jin Kuang ◽  
Yong Zheng Gu ◽  
Rong Yao ◽  
Hong Wei Wang

Calculation fluid dynamics software Fluent was used to conduct three-dimensional numerical simulation on gas-liquid two-phase flow field in a wet flue gas desulfurization scrubber. The k-ε model and SIMPLE computing were adopted in the analysis. The numerical simulation results show that the different gas entrance angles lead to internal changes of gas-liquid two-phase flow field, which provides references for reasonable parameter design of entrance angle in the scrubber.


2013 ◽  
Vol 734-737 ◽  
pp. 1488-1492
Author(s):  
Zhen Yu Liu ◽  
Li Hong Yao ◽  
Hu Zhen Wang ◽  
Cui Cui Ye

The fractures after artificial steering fracturing appear in shades of curved surface. Aiming at the problem of steering fracture, in the paper, numerical simulation method under the condition of three-dimensional two-phase flow is presented based on finite element method. In this method, of steering fracture was achieved by adopting surface elements fractures and tetrahedron elements to describe formation. By numerical simulation, the change rule of oil and water production performance of steering fractures can be calculated, and then the steering fracture parameters can be optimized before fracturing. A new method was supplied for the numerical simulation of artificial fractured well.


Metals ◽  
2019 ◽  
Vol 9 (1) ◽  
pp. 40 ◽  
Author(s):  
Md Irfanul Haque Siddiqui ◽  
Man-Hoe Kim

The sequential casting of slabs is a major trend in the steel industry where steel quality is the most important factor. The operating parameters have the most influence on mixing phenomenon apart from the design and shapes of the tundish and its furniture. Moreover, in industrial practice, the bath height in tundish varied with time when the ladle is changed. In the present work, the numerical simulation has been carried out to study the effect of residual volume and outflow (throughput) rate on the mixing phenomenon inside the tundish. A transient, three-dimensional two-phase model using the Volume of Fluid (VOF) method and Level Set interface tracking method has been used to investigate the intermixed grade steel formation. A comparison of the two interface tracking schemes, i.e., Geo-reconstruct and Modified HRIC (High-Resolution Interface Capturing Scheme) has also been presented. The results obtained through numerical simulation has been compared with experimental results. In a later section, the results showed that residual volume has a significant effect on the grade mixing. The mixing phenomenon in tundish is considerably influenced by the advance-pouring box (APB). Further, the outflow rate of tundish has little impact on the grade intermixing phenomenon.


Sign in / Sign up

Export Citation Format

Share Document