Development and Verification of Unstructured Adaptive Mesh Technique With Edge Compatibility

Author(s):  
Kei Ito ◽  
Tomoaki Kunugi ◽  
Hiroyuki Ohshima

In a design study of the large-sized sodium-cooled fast reactors in Japan (JSFR), one key issue to establish an economically superior design is suppression of a gas entrainment (GE) phenomenon at a free surface in the reactor vessel. However, the GE phenomenon is highly non-linear and too difficult to be evaluated theoretically. Therefore, we are developing a high-precision CFD method to evaluate the GE phenomenon accurately. The CFD method is formulated on an unstructured mesh to establish an accurate modeling for a complicated shape of the JSFR system. As a two-phase flow simulation method, a high-precision volume-of-fluid algorithm is employed in the CFD method. In addition, physically appropriate formulations at gas-liquid interfaces are introduced into the CFD method. The developed CFD method is already applied to the simulation of a GE phenomenon in a basic GE experiment and the simulation results show good agreement with experimental results. Therefore, it is confirmed that the proposed CFD method can reproduce a GE phenomenon. However, for the simulation of the GE phenomenon in the JSFR, we still have one problem on a mesh subdivision. Though a fine mesh subdivision has to be applied to the regions where the GE occurs, it is difficult to preliminarily know the regions because the GE occurrence is strongly affected by a local instant flow pattern, i.e. a vortex generation. Therefore, an adaptive mesh technique is necessary to apply a fine mesh subdivision automatically to only the local GE occurrence regions in the large-sized JSFR. In this study, as one part of an adaptive mesh development, a two-dimensional unstructured adaptive mesh technique is developed and verified. In the proposed two-dimensional adaptive mesh technique, each cell is isotropically subdivided to reduce distortions of the mesh. In addition, a connection cell is formed to eliminate the edge incompatibility between a refined and a non-refined cells. A connection cell has several subdivision patterns and one of them is selected to be compatible with adjacent cells on every cell edge. Finally, the present unstructured adaptive mesh technique is verified by solving well-known driven cavity problem. As the result, the present unstructured adaptive mesh technique succeeds in providing a high-precision solution, although we employ a poor-quality distorted mesh at the initial state. In addition, the simulation error on the unstructured adaptive mesh at the steady state is much less than the error on the structured mesh consisting of a larger number of cells.

Author(s):  
Kei Ito ◽  
Tomoaki Kunugi ◽  
Hiroyuki Ohshima

In the design study of large-sized sodium-cooled fast reactors in Japan (JSFR), the suppression of gas entrainment (GE) phenomena at a free surface in the reactor vessel is very important to establish an economically superior design. However, the GE phenomena are highly non-linear and too difficult to be evaluated theoretically. Therefore, we are developing high-precision CFD method for gas-liquid two-phase flows to evaluate the GE phenomena accurately. To reproduce the GE phenomena by CFDs, there are three key issues, i.e. geometry dependency, interfacial dynamics and locality. Former two issues are already addressed by employing unstructured mesh schemes and a high-precision simulation method for gas-liquid two-phase flow based on the PLIC (Piecewise Linear Interface Calculation) method, respectively. In fact, the simulation results of the GE phenomena in a simple GE experiment showed good agreements with experimental data. Recently, therefore, we focus on the locality of the GE phenomena. In our previous study (presented in ICONE17), the two-dimensional unstructured adaptive mesh technique for single-phase flows was developed to address the third issue. The isotropic cell refinement method was employed and the connection cell method was proposed to eliminate the edge incompatibility. The verification/validation results showed that the developed unstructured adaptive mesh technique succeeded in providing a high-precision solution, even though a poor-quality distorted mesh at the initial state was employed. In this study, the unstructured adaptive mesh technique is extended to the numerical simulations of gas-liquid two-phase flows. The redistribution methods of two-phase flow variables are newly developed to satisfy the conservations of the variables, i.e. the volumes of gas and liquid phases, the location of interfaces and the momentum of each phase. This improved unstructured adaptive mesh technique for gas-liquid two-phase flows is validated by solving the well-known slotted disk revolution and dam-break problems. As a result, the unstructured adaptive mesh technique succeeds in maintaining the slotted-disk shape after one revolution and shows more than first order accuracy (grid convergence) in the slotted-disk revolution problem. In addition, thanks to the momentum-conservative formulation, the dam-break phenomenon is well simulated by the unstructured adaptive mesh technique. Especially, wave-breaking phenomena are simulated by refined cells near the gasliquid interface. It should be noted that these simulation results are obtained by using relatively small number of cells because of the efficient mesh adaptation by the unstructured adaptive mesh technique.


Author(s):  
Guibo Li ◽  
Yongsheng Lian ◽  
Matt Mersch ◽  
Chris Omalley ◽  
Adam Hofmann

This paper focuses on the numerical simulation of two-phase flow near the outlet of flat fan nozzles under low operating pressure. The moment-of-fluid (MOF) method is used for the representation of the liquid gas interface and the directional split method is used for the advection of the interface. A variable density pressure projection algorithm is used for the fluid solver and a block structured adaptive mesh refinement (AMR) method is used to locally increase the resolution near interface. The internal geometry of the nozzle is defined by three parameters: nozzle inlet diameter D, V-cut or U-cut width W and V-cut or U-cut offset H. The effects of these three parameters on the fan exit angle and pressure loss through the nozzles are studied. Experiments are conducted to validate our simulation results. Finally, we give some conclusions about the modeling of the flat fan nozzle under low operating pressure based on our study.


2020 ◽  
Vol 125 ◽  
pp. 103360 ◽  
Author(s):  
Irfan Khan ◽  
Mingjun Wang ◽  
Yapei Zhang ◽  
Wenxi Tian ◽  
Guanghui Su ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document