• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.21-28
• /
• 1997

This paper presents a two-dimensional horizontal implicit model to general circulation in estuaries and coastal seas. The model is developed in non-orthogonal curvilinear coordinates system, using the Interpolating Matrix Method (IMM), in combination with a technique of multi-domain. In the propose model, the Saint-Venant equations are solved by a splitting-up technique, in the successive steps; convection, diffusion and wave propagation. The ability of the proposed model to deal with full scale nature is illustrated by the interpretation of a dye-tracing experiment in the Gironde estuary.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.2
• /
• pp.198-209
• /
• 1999

This study reports the selection of dependent variables for momentum equations in general curvilinear coordinates. Catesian, covariant and contravariant velocity components were examined for the dependent variable. The focus of present study is confined to staggered grid system Each dependent variable selected for momentum equations are tested for several flow fields. Results show that the selection of Cartesian and covariant velocity components intrinsically can not satisfy mass conservation of control volume unless additional converting processes ore used. Also, Cartesian component can only be used for the flow field in which main-flow direction does not change significantly. Convergence rate for the selection of covariant velocity component decreases quickly as with the increase of non-orthogonality of grid system. But the selection of contravariant velocity component reduces the total mass residual of discretized equations rapidly to the limit of machine accuracy and the solutions are insensitive to the main-flow direction.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.11a
• /
• pp.173-176
• /
• 2000

Finite element model based on the Naghdi's shell theory in the general tensor-based form is formulated in the present study. Partial mixed variational functional for assumed strain is formulated in order to avoid the severe locking troubles known as transverse shear and membrane locking. The proposed assumed strain element in general tensor Naghdi's shell model provides very accurate solutions for thin shells in benchmark problems. In additions, linear elastic constitutive equations are given in the general curvilinear coordinate system including anisotropic layered structures. Thus laminated composited shell structures are easily analyzed in the present formulation.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.516-523
• /
• 2009

The artificial compressibility (AC) method for the incompressible Navier-Stokes equations in the generalized curvilinear coordinates using the primitive form is implemented. The main advantage of the AC approach is that the resulting system of equations resembles the system of compressible N-S equations and can thus be integrated in time using standard, well-established time-marching methods. The errors, which are the odd-even oscillation, for pressure field in using the artificial compressibility can be eliminated by using the $4^{th}$ order artificial dissipation term which is explicitly included. Even though this paper focuses exclusively on 2D laminar flows to validate and assess the performance of this solver, this numerical method is general enough so that it can be readily extended to carry out 3D URANS simulation of engineering flows. This algorithm yields practically identical velocity profiles and primary vortex and secondary vortices that are in excellent overall agreement with the results of the vorticity-stream function formulation (Ghia et al., 1982). However, the grid resolution have to be required to be large enough to express the various vortices.

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

• 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.21 no.7
• /
• pp.873-881
• /
• 1997

A numerical simulation has been carried out for three-dimensional turbulent flows around an Ahmed body. The Reynolds-averaged Navier-Stokes equation is solved with the SIMPLE method in general curvilinear coordinates system. Several k-.epsilon. turbulence models with two convective difference schemes are evaluated for the performance such as drag coefficient, velocity and pressure fields. The drag coefficient, the velocity and pressure fields are found to be changed considerably with the adopted k-.epsilon. turbulence models as well as the finite difference schemes. The results of simulation prove that the RNG k-.epsilon. model with the QUICK scheme predicts fairly well the tendency of velocity and pressure fields and gives more reliable drag coefficient. It is also demonstrated that the large difference between simulations and experiment in the drag coefficient is due to relatively high predicted values of pressure drag from vertical rear end base.

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

The Reynolds-averaged Navier-Stokes equations with the equations of the k-.epsilon. turbulence model are solved numerically in a general curvilinear system for a three-dimensional turbulent flow around an Ahmed body. The simulation is especially aimed at the evaluation of three finite differencing schemes for the convection term, which include the upwind differencing scheme(UDS), the second order upwind differencing scheme(SOU scheme) and the QUICK scheme. The drag coefficient, the velocity and pressure fields are found to be changed considerably with the adopted finite differencing schemes. It is clearly demonstrated that the large difference between computation and experiment in the drag coefficient is due to relatively high predicted values of pressure drag from both front part and vertical rear end base. The results also show that the simulation with the QUICK or SOU scheme predicts fairly well the flow field and gives more accurate drag coefficient than other finite differencing scheme.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.11 no.2
• /
• pp.126-133
• /
• 2003

The aerodynamic characteristics of automobiles have received substantial interest recently. Detailed knowledge of the vehicle aerodynamics is essential to improve fuel efficiency and enhance stability at high-speed cruising. In this study, a numerical simulation has been carried out for three-dimensional turbulent flows around a commercial bus body. Also, the effect of rear-spoiler attached at rear end of bus body was investigated. The Wavier-Stokes equation is solved with SIMPLE method in general curvilinear coordinates system. RNG $k-\varepsilon$ turbulence model with the MARS scheme was used for the evaluating aerodynamic forces, velocity and pressure distribution. The results showed details of the three-dimensional wake flow in the immediate rear of bus body and the effect of rear-spoiler on the wake structure. A maximum of 14% reduction in drag coefficient was achieved for a model with a rear-spoiler.