• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.10a
• /
• pp.229-236
• /
• 2004

A new three-node triangular shell element based on higher order zig-zag theory is developed for laminated composite shells with multiple delaminations. The present higher order zig-zag shell theory is described in a general curvilinear coordinate system and in general tensor notation. All the complicated curvatures of surface including twisting curvatures can be described in an exact manner in the present shell element because this element is based on geometrically exact surface representation. The displacement field of the proposed finite element includes slope of deflection. which requires continuity between element interfaces. Thus the nonconforming shape function of Specht's three-node triangular plate bending element is employed to interpolate out-of-plane displacement. The present element passes the bending and twisting patch tests in flat surface configurations. The developed element is evaluated through the buckling problems of composite cylindrical shells with multiple delaminations. Through the numerical examples it is demonstrated that the proposed shell element is efficient because it has minimal degrees of freedom per node. The accuracy of the present element is demonstrated in the prediction of buckling loads and buckling modes of shells with multiple delaminations. The present shell element should serve as a powerful tool in the prediction of buckling loads and modes of multi-layered thick laminated shell structures with arbitrary-shaped multiple delaminations.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2004.10a
• /
• pp.69-72
• /
• 2004

A new three-node triangular shell element based on higher order zig-zag theory is developed for laminated composite shells with multiple delaminations. The present higher order zig-zag shell theory is described in a general curvilinear coordinate system and in general tensor notation. All the complicated curvatures of surface including twisting curvatures can be described in an exact manner in the present shell element because this element is based on geometrically exact surface representation. The displacement field of the proposed finite element includes slope of deflection, which requires continuity between element interfaces. Thus the nonconforming shape function of Specht's three-node triangular plate bending element is employed to interpolate out-of-plane displacement. The present element passes the bending and twisting patch tests in flat surface configurations. The developed element is evaluated through the eigenvalue problems of composite cylindrical shells with multiple delaminations. Through the numerical examples it is demonstrated that the proposed shell element is efficient because it has minimal degrees of freedom per node. The present shell element should serve as a powerful tool in the prediction of natural frequency and modes of multi-layered thick laminated shell structures with arbitrary-shaped multiple delaminations.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2157-2161
• /
• 2003

In this paper, the response surface method using three-dimensional Navier-Stokes analysis to optimize the shape of a forward-curved blades centrifugal fan, is described. For numerical analysis, Reynolds-averaged Navier-Stokes equations with standard k-e turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in forward-curved blades centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Linear Upwind Differencing Scheme(LUDS) is used to approximate the convection terms in the governing equations. SIMPLEC algorithm is used as a velocity-pressure correction procedure. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. It was found that the optimization process provides reliable design of this kind of fans with reasonable computing time

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.46-53
• /
• 1997

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. A H-H type of multi-block grid system is generated around a three dimensional minivan-like body. Turbulent flows have been modeled by the Baldwin-Lomax turbulent model. To validate present procedure, the flows around the Ahmed body with 12.5° of slant angle are simulated. A good agreement with other numerical results is achived. After code validation, the flows around a mimivan-like body are simulated. The simulation shows three dimensional vortex-pair just behind body. The flow separation is also observed on the rear of the body. It has concluded that the results of present study properly agreed with physical flow phenomena.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.11
• /
• pp.3573-3588
• /
• 1996

The flow through a centrifugal compressor impeller was calculated using the 3-dimensional Navier-Stokes solution method. A control volume method based on a rotating curvilinear coordinate system was used to solve the time-averaged Navier-Stokes equations, and a standard k-.epsilon. model was used to obtain eddy viscosity. Numerical results and experimental data were compared for the overall performance of the impeller, the pressure distributions along the shroud wall and the detailed flowfields at the design and off-design conditions, which showed good coincidence. The flow through the impeller is complex with the curvature of the streamlines and rotation. The development of secondary flows and the jet-wake flow characteristics, which is the main source of flow loss, was discussed. Calculation results show quite different patterns as the flow rate changes.

• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.13-21
• /
• 2013

Cavitating flow is usually formed on the surface of a high speed underwater object. When a object moves near a free surface at very high speed, the cavity signature becomes one of the major factors to be overcome by sensors of military satellite. The present work was to study the free surface effect on the ventilated cavitation. The governing equations were Navier-Stokes equations based on a homogeneous mixture model. The multiphase flow solver used an implicit preconditioning method in the curvilinear coordinate system. The cavitation model used here was the one first presented by Merkle et al.(2006) and redeveloped by Park & Ha(2009). Computations considered the free surface effects were carried out with a NACA0012 hydrofoil and the corresponding results were compared with the experimental data to have a good agreement. Calculations were then performed considering the ventilated cavitation, including the effect of non-condensable gas under the free surface effects.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.1
• /
• pp.189-205
• /
• 1996

The twin-roll type strip continuous casting process(or direct rolling process) of steel materials is characterized by two rotating water cooled rolls receiving a steady supply of molten metal which solidifies onto the rolls. A solidification analysis of molten metal considering phase transformation and thermofluid is performed using finite diffefence method with curvilinear coordinate to reduce computing time and molten region analysis with arbitrary shape. An enthalpy-specific heat method is used to determine the temperatures inthe roll and the steel. The temperature distribution of cooling roll is calculated using two dimensional finite element method, because of complex roll shape due to cooling hole in rolls and improvemnt accuracy of calculation result. The energy equaiton of cooling roll is solved simultanuously with the conservation equaiton of molten metal in order to consider heat transfer through the cooling roll. The calculated roll temperature is compared to experimental results and the heat transfer coefficient between cooling roll surface and rolling material(steel) is also determined from comparison of measured roll temperature and calculated temperature.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.10
• /
• pp.2617-2629
• /
• 1995

The flow through a centrifugal compressor rotor was calculated using the quasi-3-dimensional and fully 3-dimensional Navier-Stokes solution methods. The calculated results, obtained during the development of the computer codes for both methods are discussed. In the inviscid quasi 3-dimensional analysis, stream function formulation was used for the blade to blade (B-B) plane calculations, and the streamline curvature method was used for the meridional (H-S) plane calculations. In the viscous 3-dimensional flow analysis, a control volume method based on a general rotating curvilinear coordinate system was used to solve the time-averaged Navier-Stokes equations, and a standard k-.epsilon. model was used to obtain eddy viscosity. The quasi-3-dimensional analysis reasonably predicts the pressure distributions and requires much less computation time in the region where viscous effects are not strong; however, it fails to predict velocity field and loss mechanism through the impeller passage. The viscous 3-dimensional flow analysis shows reasonable pressure distributions and typical jet-wake flow field through the impeller passage. Secondary flow and total pressure distributions on cross-sectional planes explain the loss mechanisms through the impeller.

• Journal of Bio-Environment Control
• /
• v.3 no.1
• /
• pp.10-19
• /
• 1994

Wind load is known to be one of major forces to influence the stability of agricultural structures. General flow fields were calculated to determine flow characteristics over the envelop of the following three types of greenhouses with arched roof : single span, twin span greenhouses, and two single span greenhouses apart 3m inbetween. Pressure coefficients along the envelop of greenhouse were numerically calculated by the k-$\varepsilon$ turbulence model, which lead to determine wind forces on it. Curvilinear coordinate for an arched roof and the upwind scheme were adopted for the study. The calculated pressure coefficients were validated with the avaliable data of Japanese Standard and NGAM Standard. The Magnitude of calculated forces over the envelop was not in good accordance with data except the windward wall. Even tile data of Japanese and NGAM Standard for validation deviated a lot from each other in quantity and quality. Such discrepancy may be attributed to different geometric and/or flow configuration conditions for experiments, or the insenstivity of the k-$\varepsilon$ turbulence model to recirculation flow.

• Ocean Systems Engineering
• /
• v.1 no.1
• /
• pp.31-57
• /
• 2011

Impact pressure due to sloshing is of great concern for the ship owners, designers and builders of the LNG carriers regarding the safety of LNG containment system and hull structure. Sloshing of LNG in partially filled tank has been an active area of research with numerous experimental and numerical investigations over the past decade. In order to accurately predict the sloshing impact load, a new numerical method was developed for accurate resolution of violent sloshing flow inside a three-dimensional LNG tank including wave breaking, jet formation, gas entrapping and liquid-gas interaction. The sloshing flow inside a membrane-type LNG tank is simulated numerically using the Finite-Analytic Navier-Stokes (FANS) method. The governing equations for two-phase air and water flows are formulated in curvilinear coordinate system and discretized using the finite-analytic method on a non-staggered grid. Simulations were performed for LNG tank in transverse and longitudinal motions including horizontal, vertical, and rotational motions. The predicted impact pressures were compared with the corresponding experimental data. The validation results clearly illustrate the capability of the present two-phase FANS method for accurate prediction of impact pressure in sloshing LNG tank including violent free surface motion, three-dimensional instability and air trapping effects.