A Numerical Model for Coupling of Neutron Diffusion and Thermomechanics in Fast Burst Reactors

We develop a numerical model for coupling of neutron diffusion adn termomechanics in order to stimulate transient behavior of a fast burst reactor. The problem involves solving a set of non-linear different equations which approximate neutron diffusion, temperature change, and material behavior. With this equation set we will model the transition from a supercritical to subcritical state and possible mechanical vibration

A Comparative Study of the Harmonic and Arithmetic Averaging of Diffusion Coefficients for Non-linear Heat Conduction Problems

We perform a comparative study for the harmonic versus arithmetic averaging of the heat conduction coefficient when solving non-linear heat transfer problems. In literature, the harmonic average is the method of choice, because it is widely believed that the harmonic average is more accurate model. However, our analysis reveals that this is not necessarily true. For instance, we show a case in which the harmonic average is less accurate when a coarser mesh is used. More importantly, we demonstrated that if the boundary layers are finely resolved, then the harmonic and arithmetic averaging techniques are identical in the truncation error sense. Our analysis further reveals that the accuracy of these two techniques depends on how the physical problem is modeled

A Novel Hyperbolization Procedure for The Two-Phase Six-Equation Flow Model

We introduce a novel approach for the hyperbolization of the well-known two-phase six equation flow model. The six-equation model has been frequently used in many two-phase flow applications such as bubbly fluid flows in nuclear reactors. One major drawback of this model is that it can be arbitrarily non-hyperbolic resulting in difficulties such as numerical instability issues. Non-hyperbolic behavior can be associated with complex eigenvalues that correspond to characteristic matrix of the system. Complex eigenvalues are often due to certain flow parameter choices such as the definition of inter-facial pressure terms. In our method, we prevent the characteristic matrix receiving complex eigenvalues by fine tuning the inter-facial pressure terms with an iterative procedure. In this way, the characteristic matrix possesses all real eigenvalues meaning that the characteristic wave speeds are all real therefore the overall two-phase flowmodel becomes hyperbolic. The main advantage of this is that one can apply less diffusive highly accurate high resolution numerical schemes that often rely on explicit calculations of real eigenvalues. We note that existing non-hyperbolic models are discretized mainly based on low order highly dissipative numerical techniques in order to avoid stability issues

A Novel Hyperbolization Procedure for The Two-Phase Six-Equation Flow Model

We introduce a novel approach for the hyperbolization of the well-known two-phase six equation flow model. The six-equation model has been frequently used in many two-phase flow applications such as bubbly fluid flows in nuclear reactors. One major drawback of this model is that it can be arbitrarily non-hyperbolic resulting in difficulties such as numerical instability issues. Non-hyperbolic behavior can be associated with complex eigenvalues that correspond to characteristic matrix of the system. Complex eigenvalues are often due to certain flow parameter choices such as the definition of inter-facial pressure terms. In our method, we prevent the characteristic matrix receiving complex eigenvalues by fine tuning the inter-facial pressure terms with an iterative procedure. In this way, the characteristic matrix possesses all real eigenvalues meaning that the characteristic wave speeds are all real therefore the overall two-phase flowmodel becomes hyperbolic. The main advantage of this is that one can apply less diffusive highly accurate high resolution numerical schemes that often rely on explicit calculations of real eigenvalues. We note that existing non-hyperbolic models are discretized mainly based on low order highly dissipative numerical techniques in order to avoid stability issues

Analysis and Numerical Solution for Multi-Physics Coupling of Neutron Diffusion and Thermomechanics in Spherical Fast Burst Reactors

Coupling neutronics to thermomechanics is important for the analysis of fast burst reactors, because the criticality and safety study of fast burst reactors heavily depends on the thermomechanical behavior of fuel materials. For instance, the shut down mechanism or the transition between super and sub-critical states are driven by the fuel material expansion or contraction. The material expansion or contraction is due to temperature gradient which results from fission power. In this paper, we introduce a numerical model for coupling of neutron diffusion and thermomechanics in fast burst reactors. We also provide some analysis of the coupled system. We studied material behaviors corresponding to different levels of power pulses