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

• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.2
• /
• pp.155-174
• /
• 2006

Parabolized stability equations for compressible flows in general curvilinear coordinate system are derived to deal with a broad range of transition prediction problems on complex geometry. A highly accurate finite difference PSE code has been developed using an implicit marching procedure. Compressible and incompressible flat plate flow stability under two-dimensional and three¬dimensional disturbances has been investigated to test the present code. Results of the present computation are found to be in good agreement with the multiple scale analysis and DNS data. Stability calculation results by the present PSE code for compressible boundary layer at Mach numbers ranging from 0.02 to 1.5 are also presented and are again seen to be as accurate as the spectral method.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.134-140
• /
• 2007

Nonlinear parabolized stability equations for compressible flow in general curvilinear coordinate system are derived to deal with a broad range of transition prediction problems on complex geometry. A highly accurate finite difference PSE code has been developed using an implicit marching procedure. Blasius flow is tested. The results of the present computation show good agreement with DNS data. Nonlinear interaction can make the T-S fundamental wave more unstable and the onset of its amplitude decay is shifted downstream relative to linear case. For nonlinear calculations, rather small difference in initial amplitude can produce large change during nonlinear region. Compressible secondary instability at Mach number 1.6 is also simulated and showed that 1.1% initial amplitude for primary mode is enough to trigger the secondary growth.

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

• Journal of Advanced Marine Engineering and Technology
• /
• v.18 no.3
• /
• pp.94-101
• /
• 1994

To develope the new simulator for the analysis of fluid flow information, the influence of various convective difference schemes were evaluated. General curvilinear coordinate system with nonorthogonal grids was adopted for the successful analysis of various complex geometries. Computation results show that if we can not obtain full grid numbers within available computational environment, we need to use higher order finite difference schemes to keep the prediction accuracy.

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

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

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.1 no.1
• /
• pp.22-30
• /
• 1989

Wave propagation is changed by the effect of shoaling, current-depth refraction and shelter-ing etc. To solve these problems. numerous models have been developed. In the present study, a coordinate system is proposed based on the wave ray equation with the wave number equation including diffraction effects . The governing equation for the study was derived from the mild slope wave equation in non-steady state, including current effects (Kirby, 1986a) and trans-formed into an orthogonal curvilinear coordinate system on the basis of the wave ray equation. To obtain a numerical solution, an explicit finite difference scheme was used, and solved by the relaxation method. This model was tested for various cases: Firstly a submersed circular shoal and a constant unit depth. Secondly a submerged elliptic shoal on a slope, and finally a breakwater harbour with obliquely incident waves on a slope. The model was found to simulate the experimental results and other theoretical results in wave height and wave angle fairy well, and the applicability of the model around an arbitrary shaped coastal structure was also verified. To demonstrate the general usefullness of the present approach , the model is to be applied to a field situation with a complex bed topography.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.3
• /
• pp.158-165
• /
• 2005

In order to analyze incompressible two-phase flow, Level Set method was applied on pseudo-compressibility formulation. Level Set function is defined as a signed distance function from the phase interface, and gives the information of the each phase location and the geometric data to the flow. In this study, Level Set function transport equation was coupled with flow conservation equations, and owing to pseudo-compressibility technique we could solve the resultant vector equation iteratively. Two-phase flow analysis code was developed on general curvilinear coordinate, and numerical tests of bubble dynamics and surging wave problems demonstrate its capability successfully.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.110-119
• /
• 2006

An accurate and cost efficient method PSE is used for the stability analysis of 2D or 3D compressible boundary layers. A highly accurate finite difference PSE code has been developed at a general curvilinear coordinate system using an implicit marching procedure to deal with a broad range of transition predictions problems. Evolution of disturbances in compressible flat plate boundary layers are studied for free-stream Mach numbers ranging from 0 to 1.5. The effect of mean-flow nonparallelism is found to be weak on two dimensional waves and strong on three dimensional waves. The maximum amplification rate increases monotonically with Mach number. The present PSE solutions are compared with previous numerical investigations and experimental results and are found to be in good agreement.