• Advanced Composite Materials
• /
• v.18 no.1
• /
• pp.21-42
• /
• 2009

The vibration characteristics of post-buckled laminated composite doubly curved shells are investigated. The finite element method is used for the analysis of post-buckling and free vibration of post-buckled laminated shells. The geometric non-linear finite element model includes the general non-linear terms in the strain-displacement relationships. The shell geometry used in the present formulation is derived using an orthogonal curvilinear coordinate system. Based on the principle of virtual work the non-linear finite element equations are derived. Arc-length method is implemented to capture the load-displacement equilibrium curve. The vibration characteristics of post-buckled shell are performed using tangent stiffness obtained from the converged deflection. The code is first validated and then employed to generate numerical results. Parametric studies are performed to analyze the snapping and vibration characteristics. The relationship between loads and fundamental frequencies and between loads and the corresponding displacements are determined for various parameters such as thickness ratio and shallowness.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.11
• /
• pp.1500-1508
• /
• 2001

This paper represents the importance of dependent variables in non-orthogonal curvilinear coordinates just as the importance of those variables of convective scheme and turbulence model in computational fluid dynamics. Each of Cartesian, physical covariant and physical contravariant velocity components was tested as the dependent variables of momentum equations in the staggered grid system. In the flow past a circular cylinder, the results were computed to use each of three variables and compared to experimental data. In the skewed driven cavity flow, the results were computed to check the grid dependency of the variables. The results used in Cartesian and physical contravariant components of velocity in cylinder flow show the nearly same accuracy. In the case of Cartesian and contravariant component, the same number of vortex was predicted in the skewed driven cavity flow. Vortex strength of Cartesian component case has about 30% lower value than that of the other two cases.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.95-98
• /
• 2007

A general purpose program NUFLEX has been extended for two-phase flows with topologically complex interface and cavitation flows with liquid-vapor phase change caused by large pressure drop. In analysis of two-phase flow, the phase interfaces are tracked by employing a LS(Level Set) method. Compared with the VOF(Volume-of-Fluid} method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. Also, the cavitation process is computed by including the effects of evaporation and condensation for bubble formation and collapse as well as turbulence in flows. The volume-faction and continuity equations are adapted for cavitation models with phase change. The LS and cavitation formulation are implemented into a general purpose program for 3-D flows and verified through several test problems.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.3 no.3
• /
• pp.186-196
• /
• 1991

Numerical analysis has been performed to investigate the effects of the turbulence intensity and Prandtl number on the local heat transfer around a circular cylinder in crossflow. The governing equations were reformulated in a non-orthogonal coordinate system with Cartesian velocity components and discretised by the finite volume method with a non-staggered variable arrangement. For laminar flow, the calculations were performed for the Reynolds numbers 26 and 200. The results showed good agreement with the experimental results. For turbulent flow of the Reynolds number $1{\times}10^5$ and $2{\times}10^6$, the results showed that with an increase in the turbulent intensity in the main stream, the local Nusselt number increases in the front region of the circular cylinder. But the effect of turbulent intensity on the local Nusselt number diminishes in the wake region. The influence of Prandtl numbers show similar trend to that of turbulent intensity.

• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.46-51
• /
• 1996

The flow field around a three dimensional minivan-like body has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using second-order accurate schemes for the time derivatives, and third/second-order scheme for the spatial derivatives. The Marker-and-Cell concept is applied to efficiently solve continuity equation. The fourth -order artificial damping is added to the continuity equation for numerical stability. A H-H type multi-block grid system is generated around a three dimensional minivan-like body. Turbulent flows have been modeled by the Baldwin-Lomax turbulent model. The simulation shows three dimensional vortex-pair just behind body. And the flow separation is also observed the rear of the body. It has concluded that the results of present study properly agree with physical flow phenomena.

• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• pp.163-169
• /
• 1999

The numerical procedure has been developed for simulating incompressible viscous flow around a turbine stage with rotor-stator interaction. This study solves 2-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system. The Marker-and-Cell concept is applied to efficiently solve continuity equation. To impose an accurate boundary condition, O-H multiblocked grid system is generated. O-type grid and H-type grid is generated near and outer rotor-stator The cubic-spline interpolation is applied to handle a relative motion of a rotor to the stator. Turbulent flows have been modeled by the Baldwin- Lomax turbulent model. To validate present procedure, the time averaged pressure coefficients around the rotor and stator are compared with experiment and a good agreement obtained.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.3
• /
• pp.375-385
• /
• 1996

Turbulent flow and heat transfer characteristics around staggered tube banks were studied using the 3-D Navier-Stokes equations and energy equation governing a steady incompressible flow, which were reformulated in a non-orthogonal coordinate system with cartesian velocity components and discretized by the finite volume method with a non-staggered variable arrangement. The predicted turbulent kinetic energy using RNG $k-{\varepsilon}$ model was lower than that of standard $k-{\varepsilon}$ model but showed same result for mean flow field quantities. The prediction of the skin friction coefficient using RNG $k-{\varepsilon}$ model showed better trend with experimental data than standard $k-{\varepsilon}$ model result. The inclined flow showed higher velocity and skin friction coefficient than transverse flow because of extra strain rate ($\frac{{\partial}w}{{\partial}y}$). Also, this was why the inclined flow showed higher local heat transfer coefficient than the transverse flow.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.445-456
• /
• 1997

Continuous casting process is numerically analyzed using the continuum model in a non-orthogonal coordinate system. Flow damping in the mush is modeled by combining the viscosity dependence on liquid fraction in dilute mush and the permeability dependence on liquid fraction in concentrated mush. The effect of turbulence is indirectly considered by effective diffusivity determined elsewhere by experiment. The main objective is to investigate the effects of casting parameters such as casting speed and tundish superheat on the distribution of surface temperature, shell thickness, metallurgical length and centerline segregation. Some of the computed results are compared with available experiments, and reasonable agreements are obtained.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.86-91
• /
• 1997

This paper describes analysis of complex flow around Car-like body using Chimera grid technique. As a computational algorithm, Pullboat and Chaussee's Diagonal algorithm is selected to reduce computational time. Introducing hole points flag to this Diagonal algorithm, an algorithm for Chimera grid is generated easily. This study solves 3-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using second-order accurate schemes for the time derivatives, and third/second-order scheme for the spatial derivatives. The Marker-and-Cell concept is applied to efficiently solve continuity equation. The fourth-order artificial damping is added to the continuity equation for numerical stability, It has concluded that the results of present study properly agree with physical flow phenomena.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.2
• /
• pp.641-648
• /
• 1996

Thermophoretic particle deposition on a circular cylinder in a uniform laminar air flow was numerically investigated using a control volume method based on the generalized non-orthogonal coordinate system. Variation of air properties due to the change of temperature was taken into account. Effects of variable property on the distribution of heat transfer and deposition rates of particle were discussed. A new correlation of thermophoretic particle deposition on a circular cylinder was proposed in the present study.