• Journal of Power System Engineering
• /
• v.20 no.2
• /
• pp.5-16
• /
• 2016

Generally, methods using transfer techniques, like the transfer matrix method and the transfer stiffness coefficient method, find natural frequencies using the sign change of frequency determinants in searching frequency region. However, these methods may omit some natural frequencies when the initial frequency interval is large. The Sylvester-transfer stiffness coefficient method ("S-TSCM") can always obtain all natural frequencies in the searching frequency region even though the initial frequency interval is large. Because the S-TSCM obtain natural frequencies using the number of natural frequencies existing under a searching frequency. In this paper, the algorithm for the free vibration analysis of axisymmetric conical shells was formulated with S-TSCM. The effectiveness of S-TSCM was verified by comparing numerical results of S-TSCM with those of other methods when analyzing free vibration in two computational models: a truncated conical shell and a complete (not truncated) conical shell.

• Nuclear Engineering and Technology
• /
• v.56 no.4
• /
• pp.1310-1319
• /
• 2024

The reactor cavity cooling system (RCCS) is a passive reactor safety system commonly present in the designs of High-Temperature Gas-cooled Reactors (HTGR) that removes heat from the reactor pressure vessel by means of natural convection and radiation. It is one of the factors responsible for ensuring that the reactor does not melt down under any plausible accident scenario. For the simulation of accident scenarios, which are transient phenomena unfolding over a span of up to several days, intermediate fidelity methods and system codes must be employed to limit the models' execution time. These models can quantify radiation heat transfer well, but heat transfer caused by natural convection must be quantified with the use of correlations for the heat transfer coefficient. It is difficult to obtain reliable correlations for HTGR RCCS heat transfer coefficients experimentally due to such a system's size. They could, however, be obtained from high-fidelity steady-state simulations of RCCSs. The Rayleigh number in RCCSs is too high for using a Direct Numerical Simulation (DNS) technique; thus, a Reynolds-Averaged Navier-Stokes (RANS) approach must be employed. There are many RANS models, each performing best under different geometry and fluid flow conditions. To find the most suitable one for simulating an RCCS, the RANS models need to be validated. This work benchmarks various RANS models against three experiments performed on the HTTR RCCS Mockup by the Japanese Atomic Energy Agency (JAEA) in 1993. This facility is a 1/6 scale model of a vessel cooling system (VCS) for the High Temperature Engineering Test Reactor (HTTR), which is operated by JAEA. Multiple RANS models were evaluated on a simplified 2d-axisymmetric geometry. They were found to reproduce the experimental temperature profiles with errors of up to 22% for the lowest temperature benchmark and 15% for the higher temperature benchmarks. The results highlight that the pragmatic turbulence models need to be validated for high Rayleigh natural convection-driven flows and improved accordingly, more publicly available experimental data of RCCS resembling experiments is needed and indicate that a 2d-axisymmetric geometry approximation is likely insufficient to capture all the relevant phenomena in RCCS simulations.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.4
• /
• pp.295-302
• /
• 2009

The flow in the draft tube cone of Francis turbines operated at partial discharge is a complex hydrodynamic phenomenon where an incoming steady axisymmetric swirling flow evolves into a three-dimensional unsteady flow field with precessing helical vortex (also called vortex rope) and associated pressure fluctuations. The paper addresses the following fundamental question: is it possible to compute the circumferentially averaged flow field induced by the precessing vortex rope by using an axisymmetric turbulent swirling flow model? In other words, instead of averaging the measured or computed 3D velocity and pressure fields we would like to solve directly the circumferentially averaged governing equations. As a result, one could use a 2D axi-symmetric model instead of the full 3D flow simulation, with huge savings in both computing time and resources. In order to answer this question we first compute the axisymmetric turbulent swirling flow using available solvers by introducing a stagnant region model (SRM), essentially enforcing a unidirectional circumferentially averaged meridian flow as suggested by the experimental data. Numerical results obtained with both models are compared against measured axial and circumferential velocity profiles, as well as for the vortex rope location. Although the circumferentially averaged flow field cannot capture the unsteadiness of the 3D flow, it can be reliably used for further stability analysis, as well as for assessing and optimizing various techniques to stabilize the swirling flow. In particular, the methodology presented and validated in this paper is particularly useful in optimizing the blade design in order to reduce the stagnant region extent, thus mitigating the vortex rope and expending the operating range for Francis turbines.

• Journal of Mechanical Science and Technology
• /
• v.18 no.1
• /
• pp.153-166
• /
• 2004

This paper presents a comparative study of a fully coupled, upwind, compressible Navier-Stokes code with three two-equation models and the Baldwin-Lomax algebraic model in predicting transonic/supersonic flow. The k-$\varepsilon$ turbulence model of Abe performed well in predicting the pressure distributions and the velocity profiles near the flow separation over the axisymmetric bump, even though there were some discrepancies with the experimental data in the shear-stress distributions. Additionally, it is noted that this model has y$\^$*/ in damping functions instead of y$\^$＋/. The turbulence model of Abe and Wilcox showed better agreements in skin friction coefficient distribution with the experimental data than the other models did for a supersonic compression ramp problem. Wilcox's model seems to be more reliable than the other models in terms of numerical stability. The two-equation models revealed that the redevelopment of the boundary layer was somewhat slow downstream of the reattachment portion.

• Steel and Composite Structures
• /
• v.44 no.1
• /
• pp.81-89
• /
• 2022

The purposes of the present work it to study the effect of shear deformation on the static post-buckling response of simply supported functionally graded (FGM) axisymmetric beams based on classical, first-order, and higher-order shear deformation theories. The behavior of postbuckling is introduced based on geometric nonlinearity. The material properties of functionally graded materials (FGM) are assumed to be graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of the constituents. The equations of motion and the boundary conditions derived using Hamilton's principle. This article compares and addresses the efficiency, the applicability, and the limits of classical models, higher order models (CLT, FSDT, and HSDT) for the static post-buckling response of an asymmetrically simply supported FGM beam. The amplitude of the static post-buckling obtained a solving the nonlinear governing equations. The results showing the variation of the maximum post-buckling amplitude with the applied axial load presented, for different theory and different parameters of material and geometry. In conclusion: The shear effect found to have a significant contribution to the post-buckling behaviors of axisymmetric beams. As well as the classical beam theory CBT, underestimate the shear effect compared to higher order shear deformation theories HSDT.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.591-596
• /
• 2000

Flow through compliant tubes with linear taper in wall thickness is numerically simulated by finite element analysis. Two models are examined: a planar two-dimensional channel, and an axisymmetric tube. For verification of the numerical method, flow through a compliant stenotic vessel is simulated and compared to existing experimental data. Computational results for an axisymmetric tube show that as cross-sectional area falls with a reduction in downstream pressure, flow rate increases and reaches a maximum when the speed index (mean velocity divided by wave speed) is near unity at the point of minimum cross-section area, indicative of wave speed flow limitation or "choking" (flow speed equals wave speed) in previous one-dimensional studies. For further reductions in downstream pressure, flow rate decreases. Cross-sectional narrowing is significant but localized. When the ratio of downstream-to-upstream wall thickness is ${\le}$ 2 the area throat is located near the downstream end; as wall taper is increased to ${\ge}$ 3 the constriction moves to the upstream end of the tube. In the planar two-dimensional channel, area reduction and flow limitation are also observed when outlet pressure is decreased. In contrast to the axisymmetric case, however, the elastic wall in the two-dimensional channel forms a smooth concave surface with the area throat located near the mid-point of the elastic wall. Though flow rate reaches a maximum and then falls, the flow does not appear to be choked.

• Journal of Power System Engineering
• /
• v.19 no.6
• /
• pp.60-67
• /
• 2015

While designing and operating machines, it is very important to understand the dynamic characteristic of the machines. Authors developed the Sylvester-transfer stiffness coefficient method in order to analyze effectively the free vibration of machines or structures. The Sylvester-transfer stiffness coefficient method was derived from the combination of the Sylvester's inertia theorem and the transfer stiffness coefficient method. In this paper, the authors formulate the computational algorithm for flexural free vibration analysis of axisymmetric annular plate using the Sylvester-transfer stiffness coefficient method. To confirm the usefulness of the Sylvester-transfer stiffness coefficient method, the natural frequencies and modes for two computational models computed using the Sylvester-transfer stiffness coefficient method are compared with those computed using the exact solution and the finite element method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.5
• /
• pp.787-797
• /
• 1986

An implicit approach is employed to obtain a general boundary integral formulation of axisymmetric elastic problems in terms of a pair of singular integral equations. The corresponding kernel functions from the solutions of Navier's equation are derived by applying a three dimensional integral and a direct axisymmetrical approach. A numerical discretization schem including the evaluation of Cauchy principal values of the singular integral is described. Finally the typical axisymmetric elastic models are analyzed, i.e. the hollow sphere, the constant thickness and the V-notched round bar.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.3
• /
• pp.321-332
• /
• 2005

In this paper, the time-domain finite element formulations for axisymmetric linear viscoelastic problems, especially for the viscoelastic hollow sphere and cylinder, under various boundary conditions are presented with the theoretical solutions of them obtained by using the elastic-viscoelastic correspondence principle. It is assumed that the viscoelastic material behaves like a standard linear solid in distortion and elastically in dilatation. Numerical examples are solved based on the spherically symmetric, axisymmetric and plane strain finite element models. Good agreements are obtained between numerical and theoretical solutions, which shows the validity and accuracy of the presented method.

• Journal of computational fluids engineering
• /
• v.15 no.1
• /
• pp.1-8
• /
• 2010

RCCS is a passive safety-related system that removes the decay heat of VHTR when normal decay heat removal systems are in failure. Understanding thermo-hydraulics of RCCS is important to design a safer VHTR. RCCS consists of 292 cooling panels, which are placed in the reactor cavity. The layout of RCCS gives an idea that, for CFD simulations, cooling panels can be assumed to be one annulus tube. This assumption can reduce significantly the computational time, especially for the unsteady simulation. To simulate RCCS in an axisymmetric manner, three models were suggested and compared. Each model has (1) the same outer radius, (2) the same cross-sectional area (3) the same pressure drop, respectively, as the RCCS cooling panels. The steady-state simulation was conducted with these three models and the DO radiation model. It is found that over 90% of the heat from the outer wall of the reactor pressure vessel is transported to the RCCS by radiative heat transfer. The simulation with the third model, which has the same pressure drop as the design, estimates the closest wall temperature profiles to a thermo-hydraulic code, GAMMA+, result.