• 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 Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.10a
• /
• pp.396-401
• /
• 1997

Dynamic stability of cylindrical shells subjected to follower forces is analyzed in this paper. Motion of shells is formulated in curvilinear coordinates that is consistent with assumptions made in the Timoshenko beam and the Mindlin plate. Using the finite element method, the induced equations are reduced to an equation with finite degrees of freedom. The 9-node Lagrangian element is used, and reduced integration is used to avoid shear and membrane locking. The effects of thickness ratio on the dynamic stability of cylindrical shells are studied.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.5
• /
• pp.620-627
• /
• 2003

Radiative heat transfer in an axisymmetric enclosure with absorbing, emitting, and scattering medium is studied here by using the different methods such as MDOM, FVM, and FVM2 with emphasis on the treatment of angular derivative term, which appears in a curvilinear coordinates due to angular redistribution. After final discretization equation for FVM2 is introduced by using the step scheme and directional weights, present approach is validated by applying it to three different benchmarking problems with absorbing, emitting, and scattering medium.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.3
• /
• pp.368-375
• /
• 1987

Effective thermal conductivity of the porous nuclear fuel has been investigated numerically. Difficulties associated with irregular shape of pore have been overcome by using the Body Fitted Coordinate Systems. A computer code has been developed to solve the governing equation with appropriate boundary conditions by transforming from the Cartesian coordinates to the nonorthogonal curvilinear coordinates. The effects of the porosity have been investigated. For a convenient use of the result, a correlation equation was suggested under the assumption of circular pore. The computation results by the assumption of randomly oriented elliptic pore has been agreed more closely to existing experimental result than that by the assumption of circular pore.

• 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
• /
• v.19 no.8
• /
• pp.1950-1963
• /
• 1995

A modification of the approximate-factorization method is made to accelerate the convergency rate and to take sufficiently large Courant number without loss of accuracy. And a stable implicit finite-difference scheme for solving the incompressible Navier-Stokes equations employed above modified method is developed. In the present implicit scheme, the volume fluxes with contravariant velocity components and the pressure formulation in curvilinear coordinates is adopted. In order to satisfy the continuity condition completely and to remove spurious errors for the pressure, the Navier-Stokes equations are solved by a modified SMAC scheme using a staggered gird. The upstream-difference scheme such as the QUICK scheme is also employed to the right hand side. The implicit scheme is unconditionally stable and satisfies a diagonally dominant condition for scalar diagonal linear systems of implicit operator on the left hand side. Numerical results for some test calculations of the two-dimensional flow in a square cavity and over a backward-facing step are obtained using both usual approximate-factorization method and the modified one, and compared with each other. It is shown that the present scheme allows a sufficiently large Courant number of O(10$^{2}$) and reduces the computing time.

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

• Water for future
• /
• v.23 no.3
• /
• pp.319-328
• /
• 1990

The Parameter Identification of 2-demensional estuarine model was carried out using new output ADI-FDM numerical semi-implicit schem transformed in boundary fitted(BF) - coordinate. The hydrodynamic equations which is coupled with the transport equations were used as basic equations in the model. Thompson's equations were used to transform governing equations into rectangular plane equations and his elliptic grid generation scheme was used to generate curvilinear grid system. in BF - coordinates. The parameters to be identified are friction coefficient and disperse coefficient embedded in the governing equations. The numerical output scheme is tidally averaged salinity model in BF - coordinates. The algorithm to optimize norm of error between observations and calculations is the influence coefficinet algorithm associated with least square criterion. The lumped model is conssidered in identification. This paper was concetrated on checking whether the new output scheme might be useful to identify parameters in estuarine salinity model or not. The proposed method was tested through experimental application with hypothetical simple model. The result of the test shows that the proposed method can be used for parameter identification in estuarine model.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.55-61
• /
• 2010

A dependence of weighting parameter in preconditioning method for solving low Mach number flow with incompressible flow nature is investigated. The present preconditioning method employs a finite-difference method applied Roe‘s flux difference splitting approximation with the MUSCL-TVD scheme and 4th-order Runge-Kutta method in curvilinear coordinates. From the computational results of benchmark flows through a 2-D backward-facing step duct it is confirmed that there exists a suitable value of the weighting parameter for accurate and stable computation. A useful method to determine the weighting parameter is introduced. With this method, high accuracy and stable computational results were obtained for the flow with low Mach number in the range of Mach number less than 0.3.

• Journal of computational fluids engineering
• /
• v.16 no.1
• /
• pp.60-66
• /
• 2011

Engineering interests of submerged bodies and turbomachinery has led researchers to study various cavitation models for decades. The governing equations used for the present work are the two-phase Navier-Stokes equations with homogeneous mixture model. The solver employed on implicit dual time preconditioning algorithm in curvilinear coordinates. Three different cavitation models were applied to two axisymmetric cylinders and compared with experiments. It is concluded that the Merkle's new cavitation model has successfully accounted for cavitating flows and well captured the re-entrant jet phenomenon over the 0-caliber cylinder.