• 한국전산유체공학회:학술대회논문집
• /
• 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.

• International Journal of Fluid Machinery and Systems
• /
• v.3 no.2
• /
• pp.150-159
• /
• 2010

Laidback fan shaped film-cooling hole is formulated numerically and optimized with the help of three-dimensional numerical analysis, surrogate methods, and the multi-objective evolutionary algorithm. As Pareto optimal front produces a set of optimal solutions, the trends of objective functions with design variables are predicted by hybrid multi-objective evolutionary algorithm. The problem is defined by four geometric design variables, the injection angle of the hole, the lateral expansion angle of the diffuser, the forward expansion angle of the hole, and the ratio of the length to the diameter of the hole, to maximize the film-cooling effectiveness compromising with the aerodynamic loss. The objective function values are numerically evaluated through Reynolds- averaged Navier-Stokes analysis at the designs that are selected through the Latin hypercube sampling method. Using these numerical simulation results, the Response Surface Approximation model are constructed for each objective function and a hybrid multi-objective evolutionary algorithm is applied to obtain the Pareto optimal front. The clustered points from Pareto optimal front were evaluated by flow analysis. These designs give enhanced objective function values in comparison with the experimental designs.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.19-26
• /
• 1999

This paper describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows with free-surface. The Navier-Stokes equations governing the flows have been discretized by means of finite-difference approximations, and the resulting equations have been solved via the SIMPLE-C algorithm. The free-surface is defined by the motion of a set of marker particles and the interface behaviour was investigated by means of a "Lagrangian" technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.20-28
• /
• 2014

An efficient solution algorithm for simulating free surface problem is presented. Navier-Stokes equations for variable density incompressible flow are employed as the governing equation on Cartesian meshes. In order to describe the free surface motion efficiently, VOF(Volume Of Fluid) method utilizing THINC(Tangent of Hyperbola for Interface Capturing) scheme is employed. The most time-consuming part of the current free surface flow simulations is the solution step of the linear system, derived by the pressure Poisson equation. To solve a pressure Poisson equation efficiently, the PCG(Preconditioned Conjugate Gradient) method is utilized. This study showed that the proper application of the preconditioner is the key for the efficient solution of the free surface flow when its pressure Poisson equation is solved by the CG method. To demonstrate the efficiency of the current approach, we compared the convergence histories of different algorithms for solving the pressure Poisson equation.

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

The violent free-surface motions and the corresponding impact loads are numerically simulated by using the Moving Particle Semi-implicit (MPS) method, which was originally proposed by Koshizuka and Oka (1996) for incompressible flows. In the original MPS method, there were several shortcoming including non-optimal source term, gradient and collision models, and search of free-surface particles, which led to less-accurate fluid motions and non-physical pressure fluctuations. In the present study, how those defects can be remedied is illustrated by step-by-step improvements in respective processes of the revised MPS method. The improvement of each step is explained and numerically demonstrated. The numerical results are also compared with the experimental results of Martin and Moyce (1952) for dam-breaking problem. The current numerical results for violent free-surface motions and impact pressures are in good agreement with their experimental data.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.409-413
• /
• 2005

Even though the relative importance of length scale of flow system allow us to simplify three dimensional flow problem to one or two dimensional representation, many systems still require three dimensional analysis. The objective of this study is to develop an efficient and accurate finite element model for analyzing and predicting three dimensional flow features in natural rivers and to offend to model spreading of pollutants and transport of sediments in the future. Firstly, three dimensional Reynolds averaged Navier-Stokes equations with the hydrostatic pressure assumption in generalized curvilinear coordinates were combined with the kinematic free-surface condition. Secondly. to simulate realistic high Reynolds number flow, the model employed the Streamline Upwind/Petrov-Galerkin(SU/PG) scheme as a weighting function for the finite element method in conjunction with an appropriate turbulence model(Smagorinsky scheme for the horizontal plain and Mellor-Yamada scheme for the vertical direction). Several tests is performed for the purpose of validation and verification of the developed model. A simple rectangular channel, 5-shaped and U-shaped channel are used for tests and comparisons are made with RMA-10 model. Runs for each case is converged stably without a oscillation and calculated water-surface deformation, longitudinal and transversal velocities, and velocity vector fields are in good agreement with the results of RMA-10 model.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.11a
• /
• pp.848-854
• /
• 2006

In this study, a computational analysis system has been developed in order to investigate flow-induced vibration(FIV) phenomenon for general stator-rotor cascade configurations. Relative movement of the rotor with respect to stator is reflected by modeling independent two computational domains. Fluid domains are modeled using the unstructured grid system with dynamic moving and local deforming methods. Unsteady, Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation SST $k-\omega$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used flow computing the coupled governing equations of the fluid-structure interaction problem. Detailed FIV responses for different flow conditions are presented with respect to time and vibration characteristics are also physically investigated in the time domain.

• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.1-10
• /
• 2012

This paper presents the numerical simulation results of heat transfer and friction loss for a rotating two-pass duct with the airfoil-guide vanes in the turning region. The Kriging model is used as an optimization technique with Reynolds-averaged Navier-Stokes analysis of flow field and heat transfer with shear stress transport turbulent model. To improve the heat transfer performance, angle and location of the airfoil-guide vanes have been selected as design variables. The optimization problem has been defined as a minimization of the objective function, which is defined as a linear combination of heat transfer related term and friction loss related term with a weight factor. The airfoil-guide vanes in the turning region keep the high level of heat transfer while the friction loss has a low value. By comparing the presence or absence of airfoil-guide vanes, it is shown that the airfoil-guide vanes exhibited the best heat transfer performance to improve the blade cooling except the first passage.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.275-282
• /
• 2007

In this study, a fluid-structure interaction (FSI) analysis system has been developed in order to evaluate the turbine cascade performance with blade structural deformation effect. Relative movement of the rotor with respect to stator is reflected by modeling independent two computational domains. To consider the deformed position of rotor airfoil, dynamic moving grid method is applied. Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation SST $k-{\varepsilon}$ turbulence models are solved to predict unsteady fluid dynamic loads. A fully implicit time marching scheme based on the Newmark direct integration method with high artificial damping is used to compute the fluid-structure interaction problem. Cascade performance evaluations for different elastic axis positions are presented and compared each other. It is importantly shown that the predicted aerodynamic performance considering structural deformation effect of blade can show some deviations compared to the data generally computed from rigid blade configurations and the position of elastic axis also tend to give sensitive effect.

• Computers and Concrete
• /
• v.26 no.2
• /
• pp.137-150
• /
• 2020

In this paper, the bending behavior of single-walled carbon nanotube-reinforced composite (CNTRC) laminated plates is studied using various shear deformation plate theories. Several types of reinforcement material distributions, a uniform distribution (UD) and three functionally graded distributions (FG), are inspected. A generalized higher-order deformation plate theory is utilized to derive the field equations of the CNTRC laminated plates where an analytical technique based on Navier's series is utilized to solve the static problem for simply-supported boundary conditions. A detailed numerical analysis is carried out to examine the influence of carbon nanotube volume fraction, laminated composite structure, side-to-thickness, and aspect ratios on stresses and deflection of the CNTRC laminated plates.