• Proceedings of the Korea Association of Crystal Growth Conference
• /
• 1998.06a
• /
• pp.35-38
• /
• 1998

The momentum, heat and mass transport in the melt of several sizes of crucibles are calculated using a three dimensional numerical simulation technique with and without the k-$\varepsilon$ turbulent model. When turbulent model is not used, non-axisymmetric profiles of velocity, temperature and oxygen concentration appear in the melt of all sizes of crucibles. Axisymmetric profiles are obtained when the k-$\varepsilon$ model is adopted.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.9
• /
• pp.2463-2476
• /
• 1994

Local mean fluctuating velocity components were measured in the separating and reattaching axisymmetrc region of turbulent boundary layer over the wall of convex cylinders placed in a water tunnel by using 2-color 4-beam fiber optics laser Doppler velocimetry. Measurements were made with three different diameters of cylinders with four different diameters of cylinders with four different diameter of the obstructions. The range of Reynolds number based on step height was between 5,000 to 25,200. The study demonstrates that the reattachment length decreases with decreasing cylinder radius and is always shorter than that for the two-dimensional backward-facing step flow at the condition of the same step height. It was also observed that the turbulent kinetic energy in the recirculating region increases with an increases in the radius of convex curvature. The measured velocity field suggests that the transverse curvature can effect definitely the formation of corner eddy.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.831-838
• /
• 2004

A fluid-solid conjugate solver has been newly developed and applied to an actual engine disk system. Most of the currently available conjugate solvers lack the special thermal modeling for turbomachinery disk system applications. In the present new code, these special models are implemented to expand the applicability of the conjugate method and to reduce the required computational resources. Most of the conjugate analysis work so far are limited to the axisymmetric framework. However, the actual disk system includes several non-axisymmetric components which inevitably affect the local heat transfer phenomena. Also the previous work devoted to this area usually concentrate their efforts on the steady-state thermal field, although the one in the transient condition is more critical to the engine components. This paper presents full 3D conjugate analysis of a single stage high pressure turbine rotor-stator disk system to assess the three-dimensional effects (Fig. 1). The analysis is carried out not only in the steady-state but also in the engine accelerating transient condition. The predicted temperatures shows good agreement with measured data.

• 연구논문집
• /
• s.24
• /
• pp.81-96
• /
• 1994

To investigate failure mechanisms observed in carbon-phenolic thermal insulators, differential equations which govern the decomposition process in a deformable anisotropic porous solid are derived for three-dimensional axisymmetric constructions. The governing equations not only couple the material deformation with pore pressure, but also couple pressure and temperature, which means that heat convected by the pyrolysis gases is properly accounted for. Then the Bubnov-Galerkin finite element method is applied to these equations to transform them into a semidescrete finite element system. A thermal insulation liner in the cowl region under typical operating conditions is analyzed to find a mechanism for plylift. The results from the structural analysis show across-ply failure in the cowl zone. The mechanism for plylift is hypothesized as a sequential procedure : 1) the across-ply failure which is the precursor to plylift and 2) the local fiber buckling caused by generation of excessive in-plane compressive stress. To prevent plylift, the across-ply stress can be reduced by using appropriate material ply angles in cowl zone design.

• 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.

• Structural Engineering and Mechanics
• /
• v.67 no.2
• /
• pp.185-206
• /
• 2018

This paper studies the non-axisymmetric 3D problem on the dynamics of the moving load acting in the interior of the hollow cylinder surrounded with elastic medium and this study is made by utilizing the exact equations of elastodynamics. It is assumed that in the interior of the cylinder the point located with respect to the cylinder axis moving forces act and the distribution of these forces is non-axisymmetric and is located within a certain central angle. The solution to the problem is based on employing the moving coordinate method, on the Fourier transform with respect to the spatial coordinate indicated by the distance of the point on the cylinder axis from the point at which the moving load acts, and on the Fourier series presentation of the Fourier transforms of the sought values. Numerical results on the critical moving velocity and on the distribution of the interface normal and shear stresses are presented and discussed. In particular, it is established that the non-axisymmetricity of the moving load can decrease significantly the values of the critical velocity.

• Journal of the Society of Naval Architects of Korea
• /
• v.60 no.4
• /
• pp.240-247
• /
• 2023

In this study, morphological and hydrodynamic characteristics of the non-axisymmetric supercavity generated behind a disk-shaped cavitator were examined. By extending the previous study on axisymmetric supercavitating flow based on a boundary element method, hydrodynamic forces acting under the angle of attack condition of 0 to 30 ° and shape characteristics of the supercavity were analyzed. The results revealed that increasing the angle of attack by 30 ° reduced the length and width of the cavity by about 15% and the volume by about 40 %. An empirical formula that can quantitatively estimate the geometrical characteristics and change of the cavity was derived. It is expected that this method can be used to evaluate the shape information and force characteristics of the supercavity for the control of the vehicle in a very short time compared to the viscous analysis in the initial design stage of the supercavity underwater vehicle.

• 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 the Computational Structural Engineering Institute of Korea
• /
• v.23 no.2
• /
• pp.235-243
• /
• 2010

Three-dimensional structural analysis system for nuclear containment building is presented in this paper. This system includes high-performance plate/shell elements as finite element library. It also adopts numerical modeling technique for unbonded tendon as well as bonded tendon in prestressed concrete structures. This system is constructed by connecting several in-house program to a commercial program DIANA, and then is capable of performing nonlinear analysis for ultimate pressure capacity of nuclear containment building. Finally, three-dimensional structural analysis of CANDU-type containment building is carried out in order to test the reliability of this system. These numerical results are compared with reference values, which obtained from axisymmetric structural analysis.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.1 s.256
• /
• pp.68-75
• /
• 2007

A three-dimensional computation is conducted to simulate a three-dimensional rotating stall in a low speed axial compressor. It is generally known that a tip leakage flow has an important role on a stall inception. However, almost of researchers have taken no interest in a role of the hub-comer-stall on the rotating stall even though it is a common feature of the flow in an axial compressor operating near stall and it has a large effect on the flows and loss characteristics. Using a time-accurate unsteady simulation, it is found that the hub-comer-stall may be a trigger to collapse the axisymmetric flows under high loads. An asymmetric disturbance is initially originated in the hub-comer-stall because separations are naturally unstable flow phenomena. Then this disturbance is transferred to the tip leakage flows from the hub-comer-stall and grows to be stationary stall cells, which adheres to blade passage and rotate at the same speed as the rotor. When stationary stall cells reach a critical size, these cells then move along the blade row and become a short-length-scale rotating stall. The rotational speed of stall cells quickly comes down to 79 percent of rotor so they rotate in the opposite direction to the rotor blades in the rotating frame.