Fuel Dispersion and Flow Blockage Analyses for MYRRHA-FASTEF Reactor by SIMMER Code

Author(s):  
Giacomino Bandini ◽  
Marica Eboli ◽  
Nicola Forgione

This work illustrates the 3D set up model and the results concerning the recent analysis of fuel dispersion in the MYRRHA-FASTEF reactor performed with SIMMER code within the EU-FP7 SEARCH Project. Under severe accidental conditions, the release of fuel in the primary system can occur in case of fuel rod clad failure and degradation. Two cases were therefore taken into account, an imposed fuel release to study key parameters which influence the dispersion phenomenon and a coolant flow blockage in a fuel assembly. The reactor was simulated by a 3D Cartesian geometry with 65×63×42 cell mesh. Steady-state and transient analyses were performed by SIMMER-IV. Steady-state analysis was performed in order to assess the correct operability of the code and of the model. The results were compared with the design values. The most significant results obtained for temperature trends and profiles, velocity and mass flow rate trends are reported. Transient results were also analysed, i.e. fuel dispersion transients were simulated, comparing the effect of fuel porosity on the fuel dispersion inside the pool. In addition, the effects of the release position and the fuel particle dimension on the dispersion phenomenon were also investigated. The final section of the paper describes the effects of a flow blockage on the core degradation and dispersion of fuel particles in the primary circuit of the MYRRHA reactor. This simulation, with fuel porosity equal to 5%, started after a preliminary steady state condition. The mass flow rate in one of the inner fuel assemblies was then reduced to about 10% of the initial value. The results show that the SIMMER-IV code is capable of predicting steady-state results in good agreement with the nominal values, also confirming the correctness of the set up model.

Author(s):  
Xue-Nong Chen ◽  
Danilo D’Andrea ◽  
Claudia Matzerath Boccaccini ◽  
Werner Maschek

Safety analyses for the XT-ADS were performed with the reactor safety code SIMMER-III. Besides a brief description of the numerical model, three typical transients are presented in this paper, namely, the unprotected loss of flow (ULOF), unprotected transient over-current (UTOC), and the unprotected coolant flow blockage accident (UBA). Because of the important phenomenon of mass flow rate undershooting in the ULOF case, an integral equation model was set up for a further theoretical study of ULOF. The model confirms the numerical simulation results for various cases and gives a deeper understanding of this phenomenon. The faster the pump shut down, the larger is the undershooting of the mass flow rate. On the other hand a larger coolant cold leg area leads to a weaker undershooting. The stability analysis shows that the natural convection state is in the region of the damped oscillation for the current XT-ADS design.


Author(s):  
Xiaofeng Yang ◽  
Zhaohui Chen ◽  
Tang-Wei Kuo

Steady-state port flow simulations were carried out with a commercial three dimensional (3D) Computational Fluid Dynamics (CFD) code using Cartesian mesh with cut cells to study the prediction accuracy. The accuracy is assessed by comparing predicted and measured mass-flow rate and swirl and tumble torques at various valve lifts using different boundary condition setup and mesh topology relative to port orientation. The measured data is taken from standard steady-state flow bench tests of a production intake port. The predicted mass-flow rates agree to within 1% with the measured data between the intermediate and high valve lifts. At low valve lifts, slight over prediction in mass-flow rate can be observed. The predicted swirl and tumble torques are within 25% of the flow bench measurements. Several meshing parameters were examined in this study. These include: inlet plenum shape and outlet plenum/extension size, embedded sphere with varying minimum mesh size, finer meshes on port and valve surface, orientation of valve and port centerline relative to the mesh lines. For all model orientations examined, only the mesh topology with the valve axis aligned closely with the mesh lines can capture the mass-flow rate drop for very high valve lifts due to flow separation. This study further demonstrated that it is possible to perform 3D CFD flow analyses to adequately simulate steady-state flow bench tests.


Author(s):  
Tao Fang ◽  
Satbir Singh

Steady-state port-flow simulations with static valve lift are often utilized to optimize the performance of intake system of an internal combustion engine. Generally, increase in valve lift results in higher mass flow rate through the valve. But in certain cases, mass flow rate can actually decrease with increased valve lift, caused by separation of turbulent flow at the valve-seat. Prediction of this phenomenon using computational fluid dynamics (CFD) models is not trivial. It is found that the computational mesh significantly influences the simulation results. A series of steady-state port flow simulation are carried out using a commercial CFD code. Several mesh topologies are applied for the simulations. The predicted results are compared with available experimental data from flow bench measurements. It is found that the flow separation and reduction in mass flow rate with increased valve lift can be predicted when high mesh density is used in the proximity of the valve seat and the walls of the intake port. Higher mesh density also gives better predictions of mass flow rate compared to the experiments, but only for high valve lifts. For low valve lifts, the error in predicted flow rate is close to 13%.


Author(s):  
Basant K. Jha ◽  
Taiwo S. Yusuf

This paper presents MHD transient flow in an infinite vertical concentric annulus when the fluid is set in motion by free convection current occurring in the annulus as a result of application of isothermal heating on the inner surface of the outer cylinder while the outer surface of the inner cylinder is thermally insulated. The solution of the governing equations are obtained using the well-known Laplace transform technique while the Riemann-sum approximation method has been used to invert the solution from Laplace domain to time domain. The numerical values obtained using Riemann-sum approximation approach is validated by presenting a comparison with the values obtained using the implicit finite difference method as well as the steady-state solution. These comparisons with the steady state solution shows a remarkable agreement at large value of time. The effect of the governing parameters on the velocity field, temperature field, mass flow rate as well as the skin-friction on both surfaces of the annulus have been analysed and presented with the aid of line graph. Generally, we observed that the mass flow rate and skin friction at the isothermally heated surface increases with increase in radius ratio. However, the reverse is seen at the thermally insulated surface as the skin-friction decreases with increase in radius ratio.


Author(s):  
Alex Yeung ◽  
Nagabhushana Rao Vadlamani ◽  
Tom Hynes ◽  
Sumit Sarvankar

This paper studies the computational modelling of the flow separation over the engine nacelle lips under the off-design condition of significant crosswind. A numerical framework is set up to reproduce the general flow characteristics under crosswinds with increasing engine mass flow rate, which include: low-speed separation, attached flow and high speed shock-induced separation. A quasi-3D (Q3D) duct extraction method from the full 3D (F3D) simulations has been developed. Results obtained from the Q3D simulations are shown to largely reproduce the trends observed (isentropic Mach number variations and high-speed separation behaviour) in the 3D intake, substantially reducing the simulation time by a factor of 50. The agreement between the F3D and Q3D simulations is encouraging when the flow either fully attached or with modest levels of separation but degrades when the flow fully detaches. Results are shown to deviate beyond this limit since the captured streamtube shape (and hence the corresponding Q3D duct shape) changes with the mass flow rate. Interestingly, the drooped intake investigated in the current study is prone to earlier separation under crosswinds when compared to an axisymmetric intake. Implications of these results on the industrial nacelle lip design are also discussed.


2018 ◽  
Vol 61 (3) ◽  
pp. 887-895
Author(s):  
Michael P. Sama ◽  
Austin M. Weiss ◽  
Emma K. Benedict

Abstract. Validation of agricultural sprayers is important for quantifying as-applied coverage rates under field conditions. The complexity of modern sprayer control systems presents a challenge for precise field validation due to the use of nozzle control technologies, such as pulse width modulation, to meter chemical flow rates at individual nozzles. Non-uniform flow over time may result in local variations at high spatial resolutions that are ignored when estimating as-applied coverage rates across a field. The purpose of this study was to test several methods for estimating the mass of water applied to a water-sensitive paper spray card target using steady-state and instantaneous measurement techniques. The steady-state method consisted of a spray patternator table used to quantify the mass flow rate distribution across the nozzle width at varying nozzle pressures. The mass flow rate was then projected onto a two-dimensional area traveling across the spray width to calculate the mass of water that was deposited in the area. Two instantaneous sampling methods were used. The first method directly measured the mass of the spray card and water for 5 min after exposure to model the evaporation rate and solve for the initial mass at the time of exposure. The second method indirectly used the percent coverage of the exposed spray card by droplets. Results showed that the error between the calculated mass of water from the mass flow rate and the estimated initial mass of water from the evaporation rate varied between 2% and 8%. The relationships between the calculated and estimated initial mass of water methods and the spray card percent coverage were highly linear (R2 > 0.98). Both instantaneous methods produced results with higher variability between replications than the steady-state method, but the number of replications resulted in acceptably small differences between average mass measurements. These results show the potential for using evaporation rates for laboratory validation and percent coverage for laboratory or field validation of as-applied coverage rates. Keywords: Evaporation rate, Flow measurement, Precision agriculture, Sprayers, Water-sensitive paper.


Author(s):  
Xiaofeng Yang ◽  
Zhaohui Chen ◽  
Tang-Wei Kuo

Steady-state port flow simulations were carried out with a commercial three-dimensional (3D) computational fluid dynamics (CFD) code using Cartesian mesh with cut cells to study the prediction accuracy. The accuracy is assessed by comparing predicted and measured mass-flow rate and swirl and tumble torques at various valve lifts using different boundary condition setup and mesh topology relative to port orientation. The measured data are taken from standard steady-state flow bench tests of a production intake port. The predicted mass-flow rates agree to within 1% with the measured data between the intermediate and high valve lifts. At low valve lifts, slight overprediction in mass-flow rate can be observed. The predicted swirl and tumble torques are within 25% of the flow bench measurements. Several meshing parameters were examined in this study. These include: inlet plenum shape and outlet plenum/extension size, embedded sphere with varying minimum mesh size, finer meshes on port and valve surface, orientation of valve, and port centerline relative to the mesh lines. For all model orientations examined, only the mesh topology with the valve axis aligned closely with the mesh lines can capture the mass-flow rate drop for very high valve lifts due to flow separation. This study further demonstrated that it is possible to perform 3D CFD flow analyses to adequately simulate steady-state flow bench tests.


Author(s):  
Tao Fang ◽  
Satbir Singh

Steady-state port-flow simulations with static valve lift are often utilized to optimize the performance of intake system of an internal combustion engine. Generally, increase in valve lift results in higher mass flow rate through the valve. But in certain cases, mass flow rate can actually decrease with increased valve lift, caused by separation of turbulent flow at the valve seat. Prediction of this phenomenon using computational fluid dynamics (CFD) models is not trivial. It is found that the computational mesh significantly influences the simulation results. A series of steady-state port-flow simulations are carried out using a commercial CFD code. Several mesh topologies are applied for the simulations. The predicted results are compared with available experimental data from flow bench measurements. It is found that the flow separation and reduction in mass flow rate with increased valve lift can be predicted when high mesh density is used in the proximity of the valve seat and the walls of the intake port. Higher mesh density also gives better predictions of mass flow rate compared to the experiments, but only for high valve lifts. For low valve lifts, the error in predicted flow rate is close to 13%.


Author(s):  
I. Di Piazza ◽  
M. Tarantino ◽  
F. Magugliani ◽  
A. Alemberti

A CFD study has been carried out on fluid flow and heat transfer in the HLM-cooled Fuel Pin Bundle of the ALFRED LFR DEMO. In the context of GEN-IV Heavy Liquid Metal-cooled reactors safety studies, the flow blockage in a Fuel sub-assembly is considered one of the main issues to be addressed and the most important and realistic accident for LFR Fuel Assembly. The present paper is a first step towards a detailed analysis of such phenomena, and a CFD model and approach is presented to have a detailed thermo-fluid dynamic picture in the case of blockage. The closed hexagonal, grid-spaced fuel assembly of the LFR ALFRED has been modeled and computed. At this stage, the details of the spacer grids have not been included, but a conservative analysis has been carried out based on the current main geometrical and physical features. Reactivity feedback, as well as axial power profile, have not been included in this analysis. Results indicate that critical conditions, with clad temperatures around ∼900°C, are reached with blockage larger than 30% in terms of area fraction. Two main effects can be distinguished: a local effect in the wake/recirculation region downstream the blockage and a global effect due to the lower mass flow rate in the blocked subchannels; the former effect gives rise to a temperature peak behind the blockage and it is dominant for large blockages (>20%), while the latter effect determines a temperature peak at the end of the active region and it is dominant for small blockages (<10%). The blockage area has been placed at the beginning of the active region, so that both over-mentioned phenomena can fully take place. The mass flow rate at the different degree of blockage has been imposed from preliminary system code simulations. Transient analyses with fully resolved SST-ω turbulence model have been carried out and results indicate that a blockage of ∼15% (in terms of blocked area) leads to a maximum clad temperature around 800 °C, and this condition is reached in a characteristic time of 3–4 s without overshoot. Local clad temperatures around 1000 °C can be reached for blockages of 30% or more. CFD simulations indicate that Blockages >15% could be detected by putting some thermocouples in the plenum region of the FA.


Author(s):  
G. Belforte ◽  
T. Raparelli ◽  
V. Viktorov ◽  
A. Trivella

In porous resistances, Darcy’s law provides a good approximation of mass flow rate when the differences between upstream and downstream pressures are sufficiently small. In this range, the mass flow rates are proportional to the porous resistance’s permeability. For gas bearings, the pressure difference is normally higher, and it is known experimentally that the mass flow rates are lower than would result from Darcy’s law. Forchheimer’s law adds an inertial term to Darcy’s law and, when an appropriate coefficient is selected for this term, provides a good approximation of flow rates for the same applications even with the highest pressure differences. This paper presents an experimental and theoretical investigation of porous resistances used in gas bearing supply systems. Cylindrical sintered bronze inserts featuring lengths, diameters and particle sizes commonly used in gas bearings and thrust pads were examined. The paper describes the test set-up and experimental results obtained for: a) Mass flow rate through single porous resistances at different upstream and downstream pressures; and b) Mass flow rate and pressure distribution on a pneumatic pad featuring the same porous resistances. The theoretical permeability of the chosen porous resistances was calculated, and the results from set-up a) were then used to obtain experimental permeability and to determine the inertial coefficients. The results, which are expressed as a function of the Reynolds number, confirmed the validity of using Forchheimer’s law. The mass flow rates from set-up b) were compared with those from set-up a) at the same pressure differentials across the insert.


Sign in / Sign up

Export Citation Format

Share Document