• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.15-25
• /
• 2001

Numerical prediction of the diffusion controlled transition in a turbine gas pass is important because it can change the local heat transfer rate over a turbine blade as much as three times. In this study, the gas flow over turbine blade is simplified to the flat plate boundary layer, and an adaptive grid scheme redistributing grid points within the computation domain is proposed with a great emphasis on the construction of the grid control function. The function is sensitized to the second invariant of the mean strain tensor, its spatial gradient, and the interaction of pressure gradient and flow deformation. The transition process is assumed to be described with a κ-ε turbulence model. An elliptic solver is employed to integrate governing equations. Numerical results show that the proposed adaptive grid scheme is very effective in obtaining grid independent numerical solution with a very low grid number. It is expected that present scheme is helpful in predicting actual flow within a turbine to improve computation efficiency.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.33-38
• /
• 2000

A Numerical study of laminar vortex-shedding past a circular cylinder has been performed widely by many researchers. Many factors, such as numerical technique and domain size, number and shape of grid, affected predicting vortex shedding and Strouhal number. In the present study, the effect of convection scheme, time discretization methods and grid dependence were investigated. The present paper presents the finite volume solution of unsteady flow past circular cylinder at Re=200, 400. The Strouhal number was predicted using UDS, CDS, Hybrid, Power-law, LUDS, QUICK scheme for convection term, implicit and crank-nicolson methods for time discretization. The grid dependence was investigated using H-type mesh and O-type mesh. It also studied that the effect of mesh size of the nearest adjacent grid of circular cylinder. The effect of convection scheme is greater than the effect of time discretization on predicting Strouhal. It has been found that the predicted Strouhal number changed with mesh size and shape.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.3
• /
• pp.269-277
• /
• 2014

A numerical scheme is proposed to control the BBMB (Benjamin-Bona-Mahony-Burgers) equation, and the scheme consists of three steps. Firstly, BBMB equation is converted to a finite set of nonlinear ordinary differential equations by the quadratic B-spline finite element method in spatial. Secondly, the controller is designed based on the linear quadratic regulator (LQR) theory; Finally, the system of the closed loop compensator obtained on the basis of the previous two steps is solved by the backward Euler method. The controlled numerical solutions are obtained for various values of parameters and different initial conditions. Numerical simulations show that the scheme is efficient and feasible.

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

In the aeroacoustic application of computational fluid dynamics, the physical phenomena like the crackle in the unsteady compressible jets should be based on very time-accurate numerical solution. The accuracy of the present numerical scheme is extended to the fifth order, using the WENO filter to the sixth-order central difference computation. However, the computational capacity is very restricted by the environment of computational power, so therefore the quadrilateral adaptive grids technique is introduced for this high-order accuracy scheme. The first problem is the multi-dimensional interpolation between fine and coarse grids. Some general benchmark problems are solved to show the effectiveness of this method.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.36-42
• /
• 2006

The existing numerical approximations of convection flux, especially the spatial higher-order difference schemes, in unstructured cell-centered finite volume methods are examined in detail with each other and evaluated with respect to the accuracy through their application to a 2-D benchmark problem. Six higher-order schemes are examined, which include two second-order upwind schemes, two central difference schemes and two hybrid schemes. It is found that the 2nd-order upwind scheme by Mathur and Murthy(1997) and the central difference scheme by Demirdzic and Muzaferija(1995) have more accurate prediction performance than the other higher-order schemes used in unstructured cell-centered finite volume methods.

• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.44-53
• /
• 2007

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed for a 3D component of a nuclear system code and a component-scale analysis tool. A two-fluid three-field model is used for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. An unstructured grid is adopted for realistic simulations of the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been applied to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the modified numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing very well.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.35 no.2
• /
• pp.57-64
• /
• 1993

A modified kinematic wave model of the overland flow in vegetated filter strips was developed. The model can predict both flow depth and hydraulic radius of the flow. Existing models can predict only mean flow depth. By using the hydraulic radius, erosion, deposition and flow's transport capacity can be more rationally computed. Spacing hydraulic radius was used to compute flow's hydraulic radius. Numerical solution of the model was accomplished by using both a second-order nonlinear scheme and a linear solution scheme. The nonlinear portion of the model ensures convergence and the linear portion of the model provides rapid computations. This second-order nonlinear scheme minimizes numerical computation errors that may be caused by linearization of a nonlinear model. This model can also be applied to golf courses, parks, no-till fields to route runoff and production and attenuation of many nonpoint source pollutants.

• Coupled systems mechanics
• /
• v.7 no.6
• /
• pp.669-690
• /
• 2018

This paper deals with numerical modeling of dynamic failure phenomena in rate-sensitive brittle and/or ductile materials. To this end, a two-dimensional continuum viscodamage-embedded discontinuity model, which is based on our previous work (see Do et al. 2017), is developed. More specifically, the pre-peak nonlinear and rate-sensitive hardening response of the material behavior, representing the fracture-process zone creation, is described by a rate-dependent continuum damage model. Meanwhile, an embedded displacement discontinuity model is used to formulate the post-peak response, involving the macro-crack creation accompanied by exponential softening. The numerical implementation in the context of the finite element method exploiting the second-order mid-point scheme is discussed in detail. In order to show the performance of the model several numerical examples are included.

• Bulletin of the Korean Mathematical Society
• /
• v.47 no.6
• /
• pp.1225-1234
• /
• 2010

Numerical solutions for the fractional differential dispersion equations with nonlinear forcing terms are considered. The backward Euler finite difference scheme is applied in order to obtain numerical solutions for the equation. Existence and stability of the approximate solutions are carried out by using the right shifted Grunwald formula for the fractional derivative term in the spatial direction. Error estimate of order $O({\Delta}x+{\Delta}t)$ is obtained in the discrete $L_2$ norm. The method is applied to a linear fractional dispersion equations in order to see the theoretical order of convergence. Numerical results for a nonlinear problem show that the numerical solution approach the solution of classical diffusion equation as fractional order approaches 2.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.1
• /
• pp.1-10
• /
• 2010

In the present paper, a numerical method based on the constraint interpolation profile (CIP) method is applied for simulating two-dimensional violent sloshing problems. The free surface boundary value problem is considered as a multiphase problem which includes water and air. A stationary Cartesian grid system is adopted, and an interface capturing method is used to trace the shape of free surface profile. The CIP combined unified procedure (CCUP) scheme is applied for flow solver, and the tangent of hyperbola for interface capturing (THINC) scheme is used for interface capturing. Numerical simulations have been carried out for partially-filled 2D tanks under forced sway and roll motions at various filling depths and frequencies. The computational results are compared with experiments and/or the other numerical results to validate the present numerical method.