• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.9
• /
• pp.615-622
• /
• 2017

A normal shock wave is initially formed in the shock tube that migrates towards the closed end of the tube, which, in turn, leads to the reflection of shock. Due to the interaction of the reflected shock with the boundary layer, bifurcation of the shock wave takes place. A shock train will be generated after the bifurcated shock wave approaches the contact surface. Until now, only a few studies have been conducted to investigate this shock train phenomenon inside the shock tube. For the present study, a CFD analysis has been performed on a two dimensional axisymmetric model of a shock tube using unsteady, compressible Navier-Stokes equations. In order to investigate the detailed characteristics of the shock train phenomenon, quantitative studies have been performed by varying shock tube length, diameter under fixed diaphragm, and pressure ratio inside a shock tube.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.42 no.5
• /
• pp.151-159
• /
• 2000

Dynamic loading of structures often causes excursions of stresses will into the inelastic range and the influence of geometry changes on the response is also significant in may cases. In general , the shell structures designed according to quasi-Static analysis may collapse under condition of dynamic loading. Therefore, for a more realistic prediction on the lad carrying capacity of these shell. both material and geometric nonlinear effects should be considered. In this study , the material nonlinearity effect on the dynamic response is formulated by the elasto-viscoplastic model highly corresponding to the real behavior of the material. Also, the geometrically nonlinear behavior is taken into account using a Total Lagrangian formulation. the reinforcing bars are modeled by the equivalent steel layer at the location of reinforcements, and Von Mises yield criteria is adopted for the steel layer behavior. Also, Drucker-Prager yield criteria is applied for the behavior of concrete. the shape imperfection of dome is assumed as 'dimple type' which can be expressed Wd1=Wd0(1-(r-a)m)n while the shape imperfection of wall is assumed as sinusoidal curve which is Wwi =Wwo sin(n $\pi$y/l). In numerical test, three cases of shape imperfection of 0.0 -5.0cm(opposite direction to loading ; inner shape imperfection)and 5cm (direction to loading : outward shape imperfection) and thickness of steel layer determined by steel ratio of 0,3, and 5% were analyzed. The effect of shape imperfection and steel ratio and behavior characteristics of perfect shape shell and imperfect shape shell are identified through analysis of above mentioned numerical test. Dynamic behaviors of dome and wall according toe combination of shape imperfection and steel ratio are also discussed in this paper.

• Journal of the Korean Magnetics Society
• /
• v.6 no.2
• /
• pp.106-111
• /
• 1996

Electric machines such as motors which have rmving parts are designed for producing mechanical force or torque. The accurate calculations of electromagnetic force and torque are important in the design these machines. Electromagnetic force calculation method using the results of Finite Element Method(FEM) has been presented variously in 2-D problems. Typically the Maxwell's Stress Tensor method and the method of virtual work are used. The former calculates forces by integrating the surface force densities which can be expressed in terms of Maxwell Stress Tensor(MST), and the latter by differentiating the electromagnetic energy with respect to the virtual dis¬placement of rigid bodies of interest. In the problems including current source, magnetic vector potentials(MVP) have rmstly been used as unknown variables for field analysis by a numerical method; e. g. FEM. This paper, thus, introduces the two both methods using MVP in 3-D case. To verify the usefulness of presented methods, a solenoid model is chosen and analyzed by 3-D and axisymmetric FEM. It is found that the force calculation results are in good agreement for several mesh schemes.

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

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

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.729-732
• /
• 2010

The present study has been motivated by the development of a reliable numerical methodology for simulation of kerosene/LOx coaxial swirl injectors. To deal with thermodynamic non-ideality and anomalies of transport properties pronounced at supercritical pressures, a set of subroutine libraries has been constructed based on the cubic equations of state, and applied to an existing flamelet analysis code. For computational efficiency, two-dimensional axisymmetric RANS formulation with swirl was adopted and validated successfully against an isothermal coaxial swirling jet. For the actual problem with high pressure combustion, however, numerical results show that the RANS models yield excessive production of turbulence probably due to high density gradient magnitude in the vicinity of mixing layer of swirling film flow, and imply strongly further improvement of the turbulence models.

• Journal of Korean Dental Science
• /
• v.5 no.1
• /
• pp.13-20
• /
• 2012

Purpose: The purpose of this study was to investigate the effect of patterning on the stress distribution in the bone tissue using the finite element analysis (FEA) model. Materials and Methods: For optimal comparison, it was assumed that the implant was axisymmetric and infinitely long. The implant was assumed to be completely embedded in the infinitely long cortical bone and to have 100% bone apposition. The implant-bone interface had completely fixed boundary conditions and received an infinitely big axial load. von Mises stress and maximal principal stress were analyzed. Conventional thread and 2 or 3 patterns on the upper and lower flank of the thread were compared. Result: The surface areas of patterned implants were increased up to 106~115%. The thread with patterns distributed stress better than conventional thread. Patterning in threads may produce more stress in the implant itself, but reduce stress in the surrounding bone. Stress patterns of von Mises stress were favorable with patterns, while the maximal principal stress was increased with patterns. Patterns in the lower flank showed favorable stress distribution. Conclusion: The patterns in implant thread reduce the stress generated in surrounding bone, but the number and position of patterns were crucial factors in stress distribution.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.45-52
• /
• 2008

In order to investigate the initiation and propagation processes of a spherical detonation wave induced by direct initiation, numerical simulations were carried out using two-dimensional compressible Euler equations with an axisymmetric assumption and a one-step reaction model based on Arrhenius kinetics with various levels of ignition energy. By varying the amount of ignition energy, three typical initiation behaviors, which were subcritical, supercritical and critical regimes, were observed. Then, the ignition energy of more than $137.5{\times}10^6$ in non-dimensional value was required for initiating a spherical detonation wave, and the minimum ignition energy(i.e., critical energy) was less than that of the one-dimensional simulation reported by a previous numerical work. When the ignition energy was less than the critical energy, the blast wave generated from an ignition source continued to attenuate due to the separation of the blast wave and a reaction front. Therefore, detonation was not initiated in the subcrtical regime. When the ignition energy was more than the minimum initiation energy, the blast wave developed into a multiheaded detonation wave propagating spherically at CJ velocity, and then a cellular pattern radiated regularly out from the ignition center in the supercritical regime. The influence on ignition energy was observed in the cell width near the ignition center, but the cell width on the fully developed detonation remained constant during the expanding of detonation wave due to the consecutive formation of new triple points, regardless of ignition energy. When the ignition energy was equal to the critical energy, the decoupling of the blast wave and a reaction front appeared, as occurred in the subcrtical regime. After that, the detonation bubble induced by the local explosion behind the blast wave expanded and developed into the multiheaded detonation wave in the critical regime. Although few triple points were observed in the vicinity of the ignition core, the regularly located cellular pattern was generated after the onset of the multiheaded detonation. Then, the average cell width on the fully developed detonation was almost to that in the supercritical regime. These numerical results qualitatively agreed with previous experimental works regarding the initiation and propagation processes.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.2
• /
• pp.52-58
• /
• 2018

Hysteresis phenomena are often encountered in a wide variety of fluid flow systems used in industrial and engineering applications. Hence, in recent years, a significant amount of research been focusing on clarifying the physics of the flow hysteresis appearing during the transient change of the pressure ratios and influencing the performance of the supersonic wind tunnel. However, investigations on the hysteresis phenomenon, particularly when it occurs inside the supersonic wind tunnel, are rare. In this study, numerical simulations were carried out to investigate the hysteresis phenomena of the shock waves encountered in a supersonic wind tunnel. The unsteady and compressible flow was analyzed with an axisymmetric model, and the N-S equations were solved by using a fully implicit finite volume scheme. The optimal pressure ratio was determined from the hysteresis curves, and the results can be utilized to operate the wind tunnel efficiently.

• Journal of the Korean Institute of Gas
• /
• v.26 no.3
• /
• pp.38-44
• /
• 2022

As the supply of hydrogen charging stations for hydrogen supply accelerates due to the hydrogen economy revitalization policy, the risk of accidents is also increasing. Since most hydrogen explosion accidents lead to major accidents, it is very important to secure safety when using hydrogen energy. In order to utilize hydrogen energy, it is essential to secure the safety of hydrogen storage containers used for production, storage, and transportation of liquid hydrogen. In this paper, in order to evaluate the structural safety of a hydrogen-filled pressure vessel, the behavioral characteristics of gas pressure were analyzed by finite element analysis. SA-372 Grade J / Class 70 was used for the material of the pressure vessel, and a hexahedral mesh was applied in the analysis model considering only the 1/4 shape because the pressure vessel is axisymmetric. A finite element analysis was performed at the maximum pressure using a hydrogen gas pressure vessel, and the von Mises stress, deformation, and strain energy density of the vessel were observed.