Abstract:As compared to the two-fluid single-pressure model, the two-fluid seven-equation two-pressure model has been proved to be unconditionally well-posed in all situations, thus existing with a wide range of industrial applications. The classical 1st-order upwind scheme is widely used in existing nuclear system analysis codes such as RELAP5, CATHARE, and TRACE. However, the 1st-order upwind scheme possesses issues of serious numerical diffusion and high truncation error, thus giving rise to the challenge of accurat… Show more
“…However, high‐order unbounded schemes often lead to unphysical oscillations for the solution of sharp gradients. As shown in our previous work, unphysical oscillation occurs near discontinuities for the numerical results of TOU, CD, and SOU schemes when simulating the water faucet problem. This oscillation can be explained by Godunov's order barrier theorem .…”
Summary
Current existing main nuclear thermal‐hydraulics (T‐H) system analysis codes, such as RALAP5, TRACE, and CATHARE, play a crucial role in the nuclear engineering field for the design and safety analysis of nuclear reactor systems. However, two‐fluid model used in these T‐H system analysis codes is ill posed, easily leading to numerical oscillations, and the classical first‐order methods for temporal and special discretization are widely employed for numerical simulations, yielding excessive numerical diffusion. Two‐fluid seven‐equation two‐pressure model is of particular interest due to the inherent well‐posed advantage. Moreover, high‐order accuracy schemes have also attracted great attention to overcome the challenge of serious numerical diffusion induced by low‐order time and space schemes for accurately simulating nuclear T‐H problems. In this paper, the semi‐implicit solution algorithm with high‐order accuracy in space and time is developed for this well‐posed two‐fluid model and the robustness and accuracy are verified and assessed against several important two‐phase flow benchmark tests in the nuclear engineering T‐H field, which include two linear advection problems, the oscillation problem of the liquid column, the Ransom water faucet problem, the reversed water faucet problem, and the two‐phase shock tube problem. The following conclusions are achieved. (1) The proposed semi‐implicit solution algorithm is robust in solving two‐phase flows, even for fast transients and discontinuous solutions. (2) High‐order schemes in both time and space could prevent excessive numerical diffusion effectively and the numerical simulation results are more accurate than those of first‐order time and space schemes, which demonstrates the advantage of using high‐order schemes.
“…However, high‐order unbounded schemes often lead to unphysical oscillations for the solution of sharp gradients. As shown in our previous work, unphysical oscillation occurs near discontinuities for the numerical results of TOU, CD, and SOU schemes when simulating the water faucet problem. This oscillation can be explained by Godunov's order barrier theorem .…”
Summary
Current existing main nuclear thermal‐hydraulics (T‐H) system analysis codes, such as RALAP5, TRACE, and CATHARE, play a crucial role in the nuclear engineering field for the design and safety analysis of nuclear reactor systems. However, two‐fluid model used in these T‐H system analysis codes is ill posed, easily leading to numerical oscillations, and the classical first‐order methods for temporal and special discretization are widely employed for numerical simulations, yielding excessive numerical diffusion. Two‐fluid seven‐equation two‐pressure model is of particular interest due to the inherent well‐posed advantage. Moreover, high‐order accuracy schemes have also attracted great attention to overcome the challenge of serious numerical diffusion induced by low‐order time and space schemes for accurately simulating nuclear T‐H problems. In this paper, the semi‐implicit solution algorithm with high‐order accuracy in space and time is developed for this well‐posed two‐fluid model and the robustness and accuracy are verified and assessed against several important two‐phase flow benchmark tests in the nuclear engineering T‐H field, which include two linear advection problems, the oscillation problem of the liquid column, the Ransom water faucet problem, the reversed water faucet problem, and the two‐phase shock tube problem. The following conclusions are achieved. (1) The proposed semi‐implicit solution algorithm is robust in solving two‐phase flows, even for fast transients and discontinuous solutions. (2) High‐order schemes in both time and space could prevent excessive numerical diffusion effectively and the numerical simulation results are more accurate than those of first‐order time and space schemes, which demonstrates the advantage of using high‐order schemes.
“…Similar to Darwish and Mukalled (2003) and Tasri (2021), Wang et al (2013) studied the performance of some limiters using advection flow test cases. Kivva (2020), Wu et al (2017) and Zang et al (2015) compared some newly developed flux limiters for a convective-diffusive equation using an advection flow test case. Govin and Nair (2022) compared several flux limiters with a new high-order slope limiter using supersonic flow through a wedge as a case study.…”
Flux limiters are widely used in numerical simulations to prevent spurious oscillation in the flow with strong property gradients. However, applying flux limiter on flow without strong property gradient such as advection-diffusion flow can cause errors. This article discusses the errors caused by several flux limiters in advection-diffusion flow solution. A method for applying one-dimensional limiters to two-dimensional unstructured mesh was also suggested. The error was measured by comparing the finite volume solution of a test case with a reference solution. The study shows that the calculation error of second-order finite volume with flux limiter was higher than that of second-order finite volume without limiter. However, the error of third-order finite volume with flux limiter is less than that of second-order without flux limiter. Among the flux limiters tested in this study, Venkatakrishnan’s flux limiter produces the highest error, followed by Van leer’s limiter, EULER and SMART limiter.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.