scholarly journals Numerical Simulation of the Solid Particle Sedimentation and Bed Formation Behaviors Using a Hybrid Method

Energies ◽  
2020 ◽  
Vol 13 (19) ◽  
pp. 5018
Author(s):  
Md Abdur Rob Sheikh ◽  
Xiaoxing Liu ◽  
Tatsuya Matsumoto ◽  
Koji Morita ◽  
Liancheng Guo ◽  
...  

In the safety analysis of sodium-cooled fast reactors, numerical simulations of various thermal-hydraulic phenomena with multicomponent and multiphase flows in core disruptive accidents (CDAs) are regarded as particularly difficult. In the material relocation phase of CDAs, core debris settle down on a core support structure and/or an in-vessel retention device and form a debris bed. The bed’s shape is crucial for the subsequent relocation of the molten core and heat removal capability as well as re-criticality. In this study, a hybrid numerical simulation method, coupling the multi-fluid model of the three-dimensional fast reactor safety analysis code SIMMER-IV with the discrete element method (DEM), was applied to analyze the sedimentation and bed formation behaviors of core debris. Three-dimensional simulations were performed and compared with results obtained in a series of particle sedimentation experiments. The present simulation predicts the sedimentation behavior of mixed particles with different properties as well as homogeneous particles. The simulation results on bed shapes and particle distribution in the bed agree well with experimental measurements. They demonstrate the practicality of the present hybrid method to solid particle sedimentation and bed formation behaviors of mixed as well as homogeneous particles.

Author(s):  
Liancheng Guo ◽  
Koji Morita ◽  
Hirotaka Tagami ◽  
Yoshiharu Tobita

The postulated core disruptive accidents (CDAs) are regarded as particular difficulties in the safety analysis of liquid-metal fast reactors (LMFRs). In the CDAs, the self-leveling behavior of debris bed is a crucial issue to the relocation of molten core and heat-removal capability of the debris bed. The fast reactor safety analysis code, SIMMER-III, which is a 2D, multi-velocity-field, multiphase, multicomponent, Eulerian, fluid dynamics code coupled with a fuel-pin model and a space- and energy-dependent neutron kinetics model, was successfully applied to a series of CDA assessments. However, strong interactions among rich solid particles as well as particle characteristics in multiphase flows were not taken into consideration for fluid-dynamics models of SIMMER-III. In this article, a developed hybrid method, by coupling the discrete element method (DEM) with the multi-fluid model of SIMMER-III, is applied in the numerical simulation of self-leveling behavior in debris bed. In the coupling algorithm the motions of gas and liquid phases are solved by a time-factorization (time-splitting) method. For particles, contact forces among particles and interactions between particles and fluid phases are considered through DEM. The applicability of the method in such complicate three phase flow is validated by taking the simulation of a simplified self-leveling experiment in literature. Reasonable agreement between simulation results and corresponding experimental data shows that the present method could provide a promising means for the analysis of self-leveling behavior of debris bed in CDAs.


2020 ◽  
Vol 12 (1) ◽  
pp. 703-717
Author(s):  
Yin Wei ◽  
Wang Jiaqi ◽  
Bai Xiaomin ◽  
Sun Wenjie ◽  
Zhou Zheyuan

AbstractThis article analyzes the technical difficulties in full-section backfill mining and briefly introduces the technical principle and advantages of backfilling combined with caving fully mechanized mining (BCCFM). To reveal the strata behavior law of the BCCFM workface, this work establishes a three-dimensional numerical model and designs a simulation method by dynamically updating the modulus parameter of the filling body. By the analysis of numerical simulation, the following conclusions about strata behavior of the BCCFM workface were drawn. (1) The strata behavior of the BCCFM workface shows significant nonsymmetrical characteristics, and the pressure in the caving section is higher than that in the backfilling section. φ has the greatest influence on the backfilling section and the least influence on the caving section. C has a significant influence on the range of abutment pressure in the backfilling section. (2) There exits the transition area with strong mine pressure of the BCCFM workface. φ and C have significant effect on the degree of pressure concentration but little effect on the influence range of strong mine pressure in the transition area. (3) Under different conditions, the influence range of strong mine pressure is all less than 6 m. This article puts forward a control strategy of mine pressure in the transition area, which is appropriately improving the strength of the transition hydraulic support within the influence range (6 m) in the transition area according to the pressure concentration coefficient. The field measurement value of Ji15-31010 workface was consistent with numerical simulation, which verifies the reliability of control strategy of the BCCFM workface.


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.


2014 ◽  
Vol 1008-1009 ◽  
pp. 850-860 ◽  
Author(s):  
Zhou Wei Zhang ◽  
Jia Xing Xue ◽  
Ya Hong Wang

A calculation method for counter-current type coil-wound heat exchanger is presented for heat exchange process. The numerical simulation method is applied to determine the basic physical parameters of wound bundles. By controlling the inlet fluid velocity varying in coil-wound heat exchanger to program and calculate the iterative process. The calculation data is analyzed by comparison of numerical result and the unit three dimensional pipe bundle model was built. Studies show that the introduction of numerical simulation can simplify the pipe winding process and accelerate the calculation and design of overall configuration in coil-wound heat exchanger. This method can be applied to the physical modeling and heat transfer calculation of pipe bundles in coil wound heat exchanger, program to calculate the complex heat transfer changing with velocity and other parameters, and optimize the overall design and calculation of spiral bundles.


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.


2019 ◽  
Vol 9 (5) ◽  
pp. 847
Author(s):  
Lide Wei ◽  
Changfu Wei ◽  
Sugang Sui

This paper suggests a large-scale three-dimensional numerical simulation method to investigate the fluorine pollution near a slag yard. The large-scale three-dimensional numerical simulation method included an experimental investigation, laboratory studies of solute transport during absorption of water by soil, and large-scale three-dimensional numerical simulations of solute transport. The experimental results showed that the concentrations of fluorine from smelting slag and construction waste soil were well over the discharge limit of 0.1 kg/m3 recommended by Chinese guidelines. The key parameters of the materials used for large-scale three-dimensional numerical simulations were determined based on an experimental investigation, laboratory studies, and soil saturation of survey results and back analyses. A large-scale three-dimensional numerical simulation of solute transport was performed, and its results were compared to the experiment results. The simulation results showed that the clay near the slag had a high saturation of approximately 0.9, consistent with the survey results. Comparison of the results showed that the results of the numerical simulation of solute transport and the test results were nearly identical, and that the numerical simulation results could be used as the basis for groundwater environmental evaluation.


2012 ◽  
Vol 446-449 ◽  
pp. 3875-3878
Author(s):  
Bai Feng Ji ◽  
Wei Lian Qu

Thunderstorm microbursts, which are sources of extreme wind loadings in nature, have caused numerous structural failures, especially collapses of transmission tower around the world. Numerical simulation using computational fluid dynamics (CFD) has recently made significant progress in simulating downbursts. In this paper, transient simulation of a three-dimensional moving downburst was studied using computational fluid dynamics simulation method. Transient simulation of a three-dimensional moving downburst was conducted using time-filtered Reynolds Averaged Navier-Stokes (RANS) numerical simulation method. The three-dimensional transient wind field characteristics in a moving downburst were studied in detail. The results indicate that transient wind field characteristics in a moving downburst present quite different characteristics compared with stationary downburst at different heights and radial positions.


Sign in / Sign up

Export Citation Format

Share Document