• Journal of the Korean Society of Propulsion Engineers
• /
• v.18 no.4
• /
• pp.1-8
• /
• 2014

Overexpanded flow of an axisymmetric thruster nozzle is numerically simulated to investigate effects of nozzle pressure ratio (NPR) on the shock structure and thrust performance. The Reynolds-averaged Navier-Stokes equations with k-${\omega}$ SST turbulence model are solved utilizing FLUENT solver. As the NPR is raised, thrust performance monotonically increases with the shock structure and flow-separation point being pushed toward the nozzle exit. It is also discussed that the flow structure at nozzle-exit plane which is immediately affected by a position of nozzle-interior shocks and expansion waves, has strong influence upon the thrust performance of thruster nozzle.

• Structural Engineering and Mechanics
• /
• v.9 no.6
• /
• pp.615-636
• /
• 2000

Some elevated water tanks have failed due to torsional vibrations in past earthquakes. The overall axisymmetric structural geometry and mass distribution of such structures may leave only a small accidental eccentricity between centre of stiffness and centre of mass. Such a small accidental eccentricity is not expected to cause a torsional failure. This paper studies the possibility of amplified torsional behaviour of elevated water tanks due to such small accidental eccentricity in the elastic as well as inelastic range; using two simple idealized systems with two coupled lateral-torsional degrees of freedom. The systems are capable of retaining the characteristics of two extreme categories of water tanks namely, a) tanks on staging with less number of columns and panels and b) tanks on staging with large number of columns and panels. The study shows that the presence of a small eccentricity may lead to large displacement of the staging edge in the elastic range, if the torsional-to-lateral time period ratio $({\tau})$ of the elevated tanks lies within a critical range of 0.7< ${\tau}$ <1.25. Inelastic behaviour study reveals that such excessive displacement in some of the reinforced concrete staging elements may cause unsymmetric yielding. This may lead to progressive strength deterioration through successive yielding in same elements under cyclic loading during earthquakes. Such localized strength drop progressively develop large strength eccentricity resulting in large localized inelastic displacement and ductility demand, leading to failure. So, elevated water tanks should have ${\tau}$ outside the said critical range to avoid amplified torsional response. The tanks supported on staging with less number of columns and panels are found to have greater torsional vulnerability. Tanks located near faults seem to have torsional vulnerability for large ${\tau}$.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.1
• /
• pp.147-158
• /
• 2000

For the large scale reinforced concrete cooling tower shells, the shape imperfection can be introduced due not only to mistakes in the process of construction but also to the long term behavior of concrete. The shape imperfection evokes the additional responses such as displacements and stresses in addition to the design values. In this study, the statistical behavior of the RC cooling tower shell due to the shape imperfection is investigated using the Monte Carlo simulation. The radius of cooling tower and the shell thickness are adopted as the parameters which cause the shape imperfection. The shape imperfection is modeled as a stochastic field rather than the local one of axisymmetric or bulge type of imperfection. The randomness in the radius is shown to be more affecting the structural responses than the randomness in the shell thickness. In addition to the geometrical randomness, the effect of randomness in the modulus of elasticity on the structural response is also investigated and compared with that of the geometrical ones.

• Journal of the Society of Naval Architects of Korea
• /
• v.52 no.3
• /
• pp.236-247
• /
• 2015

In this paper, the experimental equipments for ventilated supercavitation in cavitation tunnel is constructed and the basic data of ventilated supercavitation regard to the entrainment coefficient and Froude number is fulfilled. The experiments are conducted for the disk cavitator with injecting air and the pressure inside cavity and the shape of cavity are measured. As the entrainment coefficient increases while the Froude number is kept constant, the ventilated cavitation number decreases to a minimum value which decreases no more even with increasing the air entrainment. The minimum value of ventilated cavitation number, caused by the blockage effect, decreases according to increasing the diameter ratio of test section to cavitator. The cavity length is rapidly enlarged near the minimum cavitation number. In low Froude numbers, the cavity tail is floating up due to buoyancy and the air inside the cavity is evacuated from its rear end with twin-vortex hollow tubes. However, in high Froude numbers, the buoyancy effect is almost negligible and there is no more twin-vortex tubes so that the cavity shape becomes close to axisymmetric. In order to measure the cavity length and width, the two methods, which are to be based on the cavity shapes and the maximum width of cavity, are applied. As the entrainment coefficient increases after the ventilated cavitation number gets down to the minimum cavitation number, the cavity length still increases gradually. These phenomenon can be confirmed by the measurement using the method based on the cavity shapes. On the other hand, when the method based on the maximum width of cavity is used, the length and width of the cavity agree well with a semi-empirical formular of natural cavity. So the method based on the maximum width of cavity can be a valid method for cavitator design.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.38 no.3
• /
• pp.82-91
• /
• 1996

This is a basic study for the static and dynamic analysis on the elasto-plastic and elasto-viscoplastic of an axi-symmetric shell. The objective of this study was to investigate the mechanical characteristics of a nuclear reinforced concrete containment structure, which was selected as a model, by a numerical analysis using a finite element method. The structure was modeled with discrete ring elements of 8-noded isoparametric element rotating against the symmetrical axis, and the interaction between the foundation and the structure was modeled by Winkler's model. Also, the meridional tendon was modeled with 2-node truss elements, and the hoop tendon was done with point elements in two degrees of freedom. The effect of the tendon was considered without the increasement in total degree of freedom as the stiffness matrix of modeled tendon elements was assembled on the stiffness matrix of ring elements linked with the tendon. The results obtained from the analysis of an example were summarized as follows : 1. The stresses in the hoop direction on the interior and exterior surfaces of the structure were shown in changes of similar trend, and high stresses appeared on the structure wall 2. The stresses in the meridional direction on the interior and exterior surfaces were shown in change of different trend. Especially, the stresses at the junctions between the dome and the wall and between the wall and the bottom plate of the structure were very high, compared with those at other parts of the structure. 3. The stress changes in the direction of thickness on the crown of the dome were much linearly distributed. However, as the amount of tendon increased, the stresses in the upper and lower parts of the wall established with the tendon were shown stress concentration. 4. The stress changes in the direction of thickness on the center of the structure wall was linearly distributed in the all cases, and special stress due to the use of the tendon was not shown.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.9 no.1
• /
• pp.271-281
• /
• 2005

The effects of radiation damping is used to compensate the truncated boundary which is relatively close to the structure-fluid interface in the fluid element surrounding the submerged structures. An efficient ring element is presented to model the shell and fluid element which fully utilizes the characteristics of the axisymmetry. The computational model uses the technique which separate the meridional shape and circumferential wave mode and gets similar result with the exact solution in the eigenvalues and the earthquake analysis. The fluid-structure interaction techniques is developed in the finite element analysis of two dimensional problems using the relations between pressure, nodal unknown acceleration and added mass assuming the fluid to be invicid, incompressible and irrotational. The effectiveness and efficiency of the technique is demonstrated by analyzing the free vibration and seismic analysis using the added mass matrix considering the structural deformation effect.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.1
• /
• pp.1-8
• /
• 1996

A potential-based panel method has been developed for numerical computation of wave resistance on a hybrid hydrofoil. Hybrid hydrofoil is composed of a main body, two struts and two hydrofoils. The main body, which is assumed to be an axisymmetric body for the present analysis, is normally used to support displacement of a body with its buoyancy. Normal dipoles and the sources are distributed on the body(main body, struts, hydrofoils) and the sources are distributed on the free surface. Linearized free surface and the radiation conditions are satisfied using the fourth order finite difference operator and the semi-linear pressure Kutta condition is used for the numerical computation of the hydrofoils. Poisson type free surface condition has been used for the numerical computation and hyperboloidal panel method has been used for better numerical accuracy. To verify this numeric method, model tests are performed in circulation water channel. From the comparison of experimental results with numeric ones, the present method can be used as a useful tool for the design of high speed vessels.

• Journal of Ocean Engineering and Technology
• /
• v.30 no.6
• /
• pp.474-483
• /
• 2016

Accidental events such as collisions, groundings, and hydrocarbon explosions in marine structures can cause catastrophic damage. Thus, it is extremely important to predict the extent of such damage, which determines the total amount of oil spills and the residual hull girder strength. Punching fracture tests were conducted by Choung (2009b), where various sizes of indenters and circular unstiffened steel plates with different thicknesses were used to quasi-statically realize damage extents. A three-dimensional fracture strain surface was developed based on a reference (Choung et al., 2015b), where the average stress triaxiality and average normalized Lode angle were used as the parameters governing the fracture of ductile steels. In this study, new numerical analyses were performed using very fine axisymmetric elements in combination with an Abaqus user-subroutine to implement the three-dimensional fracture strain surface. Conventional numerical analyses were also conducted for the tests to identify the best fit fracture strain values by changing the fracture strains. Based on the phenomenon of the average normalized Lode angle starting out positive and then becoming slightly negative, it was inferred that the shear stress primarily dominates in determining the fractures locations, with a partial contribution from the compressive stress. It should be stated that the three-dimensional fracture surface effectively predicted at least the shear stress-dominant fracture behavior of a mild steel.

• The Bulletin of The Korean Astronomical Society
• /
• v.41 no.1
• /
• pp.40.2-41
• /
• 2016

To understand the observed kinematics in the narrow-line region (NLR) of type 2 AGNs, we construct a model of 3-D biconical outflow combined with a thin dust plane. The model consists of two identical cones whose apex is located at the nucleus, and the cones are axisymmetric with respect to the bicone axis. After we define the properties of the bicone and the dust plane, we calculate a spatially integrated velocity and velocity dispersion along the line-of-sight using various physical parameters. As we test the effect of model parameters, we find three key parameters determining the integrated kinematics: intrinsic outflow velocity, bicone inclination, and the amount of dust extinction. The velocity dispersion increases as the intrinsic outflow velocity or the bicone inclination increases, while the velocity shift increases as the amount of dust extinction increases. We confirm that the integrated velocity dispersion can be a good indicator of the intrinsic outflow velocity unless dust extinction is not very strong (>~80%), while the effect of dust extinction can be alleviated by combining the integrated velocity and the velocity dispersion. Based on the simulated velocity distributions using the 3-D models, the variety of the observed [O $_{III}$] line profiles of type 2 AGNs can be well reproduced. In addition, we perform Monte Carlo simulations based on the different sets of model parameters. By comparing the model results with the observed [O $_{III}$] kinematics of ~39,000 SDSS type 2 AGNs (Woo et al. 2016), we find that the observed [O $_{III}$] velocity-velocity dispersion distribution is well reproduced by the biconical outflow model, enabling us to constrain the intrinsic physical parameters of outflows.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.10a
• /
• pp.133-138
• /
• 2002

A 3-D axisymmetric computer program is developed to predict the NO behavior in SNCR system for the stoker incinerator with the waste treatment capacity, 200ton/day. To this end a turbulent reacting flow field calculation is made using proper assumption and empiricism. The stoker bed is assumed to be a homogeneous waste-volatilized gaseous state. The initial composition or reactants are assumed based on the data of the ultimate analysis. Turbulent is resolved by k-e model and turbulent reaction is handled by eddy-breakup model harmonized with empirical chemistry data for gaseous combustion, NO and urea reaction. The liquid droplet is traced by Lagrangian method incorporated by aerodynamic drag, Coriolis and crntrifugal forces. Radiation is treated by sensible heat loss model. Calculation results are in good agreement with experimental data at the outlet of post combustion chamber in Daejon 4th industrial complex. The flue gas shows the temperature range of $900\sim1000^{\circ}C$, velocity of 5m/s and NO concentration of 140ppm at the exit while the measured temperature, flue gas velocity and NO concentration are $967^{\circ}C$, $3\sim4m/s$ and $100\sim200ppm$respectively. Using the developed computer program a parametric study has been made with the variation of heat content of waste, castable length and SNCR variables for the determination of proper injector location. In general, the calculated results are consistent and physically acceptable.