• Steel and Composite Structures
• /
• v.24 no.2
• /
• pp.151-176
• /
• 2017

This paper deals with general equations of motion for free vibration analysis response of thick three-layer doubly curved sandwich panels (DCSP) under simply supported boundary conditions (BCs) using higher order shear deformation theory. In this model, the face sheets are orthotropic laminated composite that follow the first order shear deformation theory (FSDT) based on Rissners-Mindlin (RM) kinematics field. The core is made of orthotropic material and its in-plane transverse displacements are modeled using the third order of the Taylor's series extension. It provides the potentiality for considering both compressible and incompressible cores. To find these equations and boundary conditions, Hamilton's principle is used. Also, the effect of trapezoidal shape factor for cross-section of curved panel element ($1{\pm}z/R$) is considered. The natural frequency parameters of DCSP are obtained using Galerkin Method. Convergence studies are performed with the appropriate formulas in general form for three-layer sandwich plate, cylindrical and spherical shells (both deep and shallow). The influences of core stiffness, ratio of core to face sheets thickness and radii of curvatures are investigated. Finally, for the first time, an optimum range for the core to face sheet stiffness ratio by considering the existence of in-plane stress which significantly affects the natural frequencies of DCSP are presented.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.239-244
• /
• 2009

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved by successively in the moving least square sense. Some weighing functions were tested in order to investigate the up-winding effect for the convection term. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as FVM.

• Journal of Ocean Engineering and Technology
• /
• v.21 no.4
• /
• pp.8-14
• /
• 2007

A particle method recognized as one of gridless methods has been developed to investigate incompressible viscous flaw. The method is more feasible and effective than conventional grid methods for solving the flaw field with complicated boundary shapes or multiple bodies. The method is consists of particle interaction models representing pressure gradient, diffusion, incompressibility and the boundary conditions. In the present study, the models in case of various simulation condition were checked with the analytic solution, and applied to the two-dimensional Poiseuille flow in order to validate the developed method.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

• Computational Structural Engineering
• /
• v.8 no.4
• /
• pp.189-198
• /
• 1995

The conventional finite element formulations for incompressible materials show pressure oscillations or pressure modes in four-node quadrilateral elements of commonly used displacement and pressure interpolations. The criterion for the stability in the pressure solution is the so-called Babugka-Brezzi stability condition, and the above elements do not satisfy this condition. In this study, a pressure continuity residual based on the pressure discontinuity at element interfaces is used to study the stabilization of pressure solutions in bilinear displacement-constant pressure four-node quadrilateral elements. This pressure residual is implemented in Q1P0 element derived from the conventional incompressible elasticity. The pressure solutions can be stable with the pressure residual though they exhibit sensitivity to the stabilization parameters. Parametric study for the solution stabilization is also discussed.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.71-80
• /
• 1996

An iterative time marching procedure for solving incompressible turbulent flow has been applied to the flows around a high speed train including cross-wind effects. This procedure solves three-dimensional unsteady incompressible Reynolds-averaged Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using first-order accurate schemes for the time derivatives and third/second-order accurate schemes for the spatial derivatives. Turbulent flows have been modeled by Baldwin-Lomax turbulent model. To validate present procedure, the flow around a high speed train at zero yaw angle was simulated and compared with experimental data. Generally good agreement with experiments was achieved. The flow fields around the high speed train at 9.2°, 16.7°, and 45° of yaw angle were also simulated.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.1
• /
• pp.108-116
• /
• 2001

The vortex tube is a simple device which separates fluid stream into a cold stream and a hot stream without any chemical reaction. The process of energy separation in the vortex tube has caused a great deal of interest. Although many studies on energy separation in the vortex tube using air as the working fluid have been made so far, few experimental studies treated energy separation for incompressible fluid. So, an experimental study for the energy separation in the vortex tube using the water which is essentially an incompressible fluid is presented. When working fluid is the water, the best geometric values of nozzle area ratio and number of nozzle holes are 0.155, 6 respectively. These geometric values are showed by the similar values which are presented by compressible fluid as working fluid. But hot side mass fraction of which maximum temperature drop is happened are different from compressible fluid.

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

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.49-55
• /
• 2013

The goal of the research is to evaluate the open source code of OpenFOAM for the use of nuclear plant flow simulation objectively. Of the various incompressible flow solvers, simpleFoam, pimpelFoam are then tested under three validated cases (backward facing step, flow over circular cylinder and turbulent round jet flow). For the evaluation of steady state incompressible laminar flow simulation, low reynolds number of backward facing step flow was solved by simpleFoam. The resultant of the reattached lengths turned out to be similar with the other experimental and simulation results. For transient flow simulation, flow over circular cylinder and turbulent round jet flow were solved by pimpleFoam. The simulation accuracy was evaluated by comparing the resultant flow patterns with the description of the characteristics of the flow over the circular cylinder. The quantitative accuracy was evaluated for no more than 85% by comparing it to the decaying constants of the turbulent round jet velocity.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.5
• /
• pp.755-761
• /
• 1987

Two-Dimensional incompressible laminar boundary layer with the reversed flow region is computed using the parially parabolized Navier-Stokes equations in primitive variables. The velocities and the pressure are explicity coupled in the difference equation and the resulting penta-diagonal matrix equations are solved by a streamwise marching technique. The test calculations for the trailing edge region of a finite flat plate and Howarth's linearly retarding flows demonstrate that the method is accurate, efficient and capable of predicting the reversed flow region.