• Nuclear Engineering and Technology
• /
• v.49 no.5
• /
• pp.953-967
• /
• 2017

In this paper, a complete station blackout (SBO) or complete loss of electrical power supplies is simulated and analyzed in a downward cooling 5-MW pool-type Material Testing Reactor (MTR). The scenario is traced in the absence of active cooling systems and operators. The code nodalization is successfully benchmarked against experimental data of the reactor's operating parameters. The passive heat removal system includes downward water cooling after pump breakdown by the force of gravity (where the coolant streams down to the unfilled portion of the holdup tank), safety flapper opening, flow reversal from a downward to an upward cooling direction, and then the upward free convection heat removal throughout the flapper safety valve, lower plenum, and fuel assemblies. Both short-term and long-term natural core cooling conditions are simulated and investigated using the RELAP5 code. Short-term analyses focus on the safety flapper valve operation and flow reversal mode. Long-term analyses include simulation of both complete SBO and long-term operation of the free convection mode. Results are promising for pool-type MTRs because this allows operators to investigate RELAP code abilities for MTR thermal-hydraulic simulations without any oscillation; moreover, the Tehran Research Reactor is conservatively safe against the complete SBO and long-term free convection operation.

• East Asian mathematical journal
• /
• v.35 no.5
• /
• pp.569-587
• /
• 2019

In this paper, we consider a split least-squares characteristic mixed element method for Sobolev equations with a convection term. First, to manipulate both convection term and time derivative term efficiently, we apply a characteristic mixed element method to get the system of equations in the primal unknown and the flux unknown and then get a least-squares minimization problem and a least-squares characteristic mixed element scheme. Finally, we obtain a split least-squares characteristic mixed element scheme for the given problem whose system is uncoupled in the unknowns. We prove the optimal order in $L^2$ and $H^1$ normed spaces for the primal unknown and the suboptimal order in $L^2$ normed space for the flux unknown.

• Communications of the Korean Mathematical Society
• /
• v.16 no.3
• /
• pp.405-413
• /
• 2001

The inverse boundary value problem for the steady state heat equation with convection term is considered in a simply connected bounded domain with smooth boundary. Taking the Dirichlet to Neumann map which maps the temperature on the boundary to the that flux on the boundary as an observation data, the global uniqueness for identifying the convection term from the Dirichlet to Neumann map is proved.

• Journal of applied mathematics & informatics
• /
• v.27 no.5_6
• /
• pp.1279-1292
• /
• 2009

In this paper, a numerical method that uses standard finite difference scheme defined on Shishkin mesh for a weakly coupled system of two singularly perturbed convection-diffusion second order ordinary differential equations with a discontinuous source term is presented. An error estimate is derived to show that the method is uniformly convergent with respect to the singular perturbation parameter. Numerical results are presented to illustrate the theoretical results.

• Journal of applied mathematics & informatics
• /
• v.34 no.1_2
• /
• pp.19-34
• /
• 2016

In this paper, we present a split least-squares characteristic mixed finite element method(MFEM) to get the approximate solutions of the convection dominated Sobolev equations. First, to manage both convection term and time derivative term efficiently, we apply a least-squares characteristic MFEM to get the system of equations in the primal unknown and the flux unknown. Then, we obtain a split least-squares characteristic MFEM to convert the coupled system in two unknowns derived from the least-squares characteristic MFEM into two uncoupled systems in the unknowns. We theoretically prove that the approximations constructed by the split least-squares characteristic MFEM converge with the optimal order in L² and H¹ normed spaces for the primal unknown and with the optimal order in L² normed space for the flux unknown. And we provide some numerical results to confirm the validity of our theoretical results.

• Journal of applied mathematics & informatics
• /
• v.37 no.5_6
• /
• pp.357-382
• /
• 2019

In this paper, we consider boundary value problem for a weakly coupled system of two singularly perturbed differential equations of convection diffusion type with discontinuous source term. In general, solution of this type of problems exhibits interior and boundary layers. A numerical method based on streamline diffusiom finite element and Shishkin meshes is presented. We derive an error estimate of order $O(N^{-2}\;{\ln}^2\;N$) in the maximum norm with respect to the perturbation parameters. Numerical experiments are also presented to support our theoritical results.

• Kyungpook Mathematical Journal
• /
• v.64 no.1
• /
• pp.113-131
• /
• 2024

In this paper we give conditions for the existence of, and describe the asymtotic behavior of, radial positive solutions of the nonlinear elliptic equation of Emden-Fowler type with convection term ∆_p u + 𝛼|u|^q-1u + 𝛽x.∇(|u|^q-1u) = 0 for x ∈ ℝ^N, where p > 2, q > 1, N ≥ 1, 𝛼 > 0, 𝛽 > 0 and ∆_p is the p-Laplacian operator. In particular, we determine ${\lim}_{r{\rightarrow}}{\infty}\,r^{\frac{p}{q+1-p}}\,u(r)$ when $\frac{{\alpha}}{{\beta}}$ > N > p and $q\,{\geq}\,{\frac{N(p-1)+p}{N-p}}$.

• Journal of applied mathematics & informatics
• /
• v.33 no.5_6
• /
• pp.635-648
• /
• 2015

We consider a weakly coupled system of singularly perturbed convection-diffusion equations with discontinuous source term. The diffusion term of each equation is associated with a small positive parameter of different magnitude. Presence of discontinuity and different parameters creates boundary and weak interior layers that overlap and interact. A numerical method is constructed for this problem which involves an appropriate piecewise uniform Shishkin mesh. The numerical approximations are proved to converge to the continuous solutions uniformly with respect to the singular perturbation parameters. Numerical results are presented which illustrates the theoretical results.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2005.06a
• /
• pp.210-215
• /
• 2005

A three-dimensional numerical simulation is performed to investigate on a low Reynolds number mixed convection in a horizontal rectangular channel with the upper part cooled and the lower part heated uniformly. The three-dimensional governing equations are solved using a finite volume method. For convective term, the central differencing scheme is used and for the pressure correction, the PISO algorithm is used. Solutions are obtained for A=4, Pr=0.72, 10, 909, the Reynolds number ranging from $2.1{\times}10^{-2}$ to $1.2{\times}10^{-1}$, the Rayleigh number is $3.5{\times}10^4$. It is found that vortex roll structures of mixed convection in horizontal rectangular channel can be classified into three roll structures which affected by Prandtl number and Reynolds number.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.5
• /
• pp.619-626
• /
• 2001

The present experimental and numerical study investigate flow and natural convection heat transfer characteristics of a composed rectangular-parallelogrammic enclosure with a guide vane. The governing equations for the two-dimensional, laminar, natural convection process in an enclosure are discretized by the control volume approach which insures the conservative characteristics to be satisfied in the calculation domain, and solved by a modified SIMPLE algorithm. The momentum and energy equations are coupled through the buoyancy term. In this results of experimental study, the natural convection heat transfer characteristics was well coincided with conclusions of other earlier experimental researches and numerical analysis.