• Journal of the Korean Society of Combustion
• /
• v.6 no.2
• /
• pp.15-22
• /
• 2001

A pressure-based, unstructured finite volume method has been applied to couple the chemical kinetics and fluid dynamics and to capture effectively and accurately the steep gradient flame field. The pressure-velocity coupling is handled by two methodologies including the pressure-correction algorithm and the projection scheme. A stiff, operator-split projection scheme for the detailed nonequilibrium chemistry has been employed to treat the stiff reaction source terms. The conservative form of the governing equations are integrated over a cell-centered control volume with collocated storage for all transport variables. Computations using detailed chemistry and variable transport properties were performed for two laminar reacting flows: a counterflow hydrogen-air diffusion flame and a lifted methane-air triple flame. Numerical results favorably agree with measurements in terms of the detailed flame structure.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.39-47
• /
• 2006

In this paper, the authors developed a multidimensional interpolation method inside a finite volume cell in the computation of high-order accurate numerical flux such as the fifth order WEND (weighted essentially non-oscillatory) scheme. This numerical method starts from a simple Taylor series expansion in a proper spatial order of accuracy, and the WEND filter is used for the reconstruction of sharp nonlinear waves like shocks in the compressible flow. Two kinds of interpolations are developed: one is for the cell-averaged values of conservative variables divided in one mother cell (Type 1), and the other is for the vertex values in the individual cells (Type 2). The result of the present study can be directly used to the cell refinement as well as the convective flux between finer and coarser cells in the Cartesian adaptive grid system (Type 1) and to the post-processing as well as the viscous flux in the Navier-Stokes equations on any types of structured and unstructured grids (Type 2).

• Journal of Mechanical Science and Technology
• /
• v.14 no.7
• /
• pp.789-795
• /
• 2000

Free-surface flows with an arbitrary deformation, induced by a submerged hydrofoil, are simulated numerically, considering two-fluid flows of both water and air. The computation is performed by a finite volume method using unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell wise local mesh refinement. The integration in space is of second order, based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels. The linear equations are solved by conjugate gradient type solvers, and the non-linearity of equations is accounted for through Picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations, the continuity equation, the conservation equation of one species, and the equations for two turbulence quantities. Finally, a comparison is quantitatively made at the same speed between the computation and experiment in which the grid sensitivity is numerically checked.

• Journal of applied mathematics & informatics
• /
• v.12 no.1_2
• /
• pp.335-349
• /
• 2003

Almost all business are affected by the weather so that weather derivatives has been traded to hedge weather risk. Since the weather itself is not an asset with a market price, some analysts believe that the Black-Scholes equation could not be used appropriately to price weather derivative options. But some weather derivatives can be considered as an Asian option, we revisit the Black-scholes model. Numerical solution of the Black-Scholes equation has a significant error at the money option or around the money option, it is necessary to adopt adaptive mesh near to the strike value. Here we propose a numerical method with an adaptive grid refinement.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.771-776
• /
• 2000

Fine grid calculations are reported for the developing turbulent flow in a straight duct and a curved duct of square cross-section with a radius of curvature to hydraulic diameter ratio ${\delta}=R_c/H_H=3.357$ and a bend angle of 180 deg. A sequence of modeling refinements is introduced; the replacement of wall function by a fine mesh across the sublayer and a low Reynolds number second moment closure up to the near wall sublayer in which the non-linear return to isotropy model and the cubic-quasi-isotropy model for the pressure strain are adopted; and the introduction of a multiple source model for the exact dissipation rate equation. Each refinement is shown to lead to an appreciable improvement in the agreement between measurement and computation.

• Journal of Korea Water Resources Association
• /
• v.46 no.6
• /
• pp.655-666
• /
• 2013

Lately, intensity and frequency of natural disasters such as flood are increasing because of abnormal climate. Casualties and property damages due to large-scale floods such as Typhoon Rusa in 2002 and Typhoon Maemi in 2003 rapidly increased, and these show the limits of the existing disaster prevention measures and flood forecasting systems regarding irregular climate changes. In order to efficiently respond to extraordinary flood, it is important to provide effective countermeasures through an inundation model that can accurately simulate flood inundation patterns. However, the existing flood inundation analysis model has problems such as excessive take of analysis time and accuracy of the analyzed results. Therefore, this study conducted a flood inundation analysis by using the Gerris flow solver that uses quadtree grid, targeting the Baeksan Levee in the Nakdong River Basin that collapsed because of a concentrated torrential rainfall in August, 2002. Through comparisons with the FLUMEN model that uses unstructured grid among the existing flood inundation models and the actual flooded areas, it determined the applicability and efficiency of the quadtree grid-based flood inundation model of the Gerris flow solver.

• Journal of Software Engineering Society
• /
• v.23 no.4
• /
• pp.151-163
• /
• 2010

In the defense software domain where large-scale software products in various application areas need to be built, reusing software is regarded as one of the important practices to build software products efficiently and economically. There have been many efforts to apply various methods to support software reuse in the defense software domain. However, developers in the defense software domain still experience many difficulties and face obstacles in reusing software assets. In this paper, we analyze practical problems of software reuse in the defense software domain, and define core requirements to solve those problems. To meet these requirements, we are currently developing the Component Grid system, a reuse-support system that provides a developer-centric software reuse environment. We have designed an architecture of Component Grid, and defined essential elements of the architecture. We have also developed the core approaches for developing the Component Grid system: a semantic-tagging-based requirement tracing method, a reuse-knowledge representation model, a social-network-based asset search method, a web-based asset management environment, and a wiki-based collaborative and participative knowledge construction and refinement method. We expect that the Component Grid system will contribute to increase the reusability of software assets in the defense software domain by providing the environment that supports transparent and efficient sharing and reuse of software assets.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.4
• /
• pp.313-325
• /
• 2011

Finite element grid system using h-refinement on the Yellow Sea was constructed based on previous study (Suh, 1999b) from 14 K to 210 K and special attention was concentrated on refining the coastal zone. In grid generation, depth change between adjacent points and non-dimensional tidal wave length ratio were considered. As a result approximately a quarter of the total nodes are located nearby 5 m of shallow area. Accurate bathymetry data using 30's and ETOPO1 with open boundary conditions of 8 major tidal constituents extracted automatically from FES2004 have been applied. In tidal simulation a 3-dimensional nonlinear harmonic model was setup and tidal amplification due to changes in vertical turbulent and bottom friction were simulated. In this study not only 8 major tidal constituents but also nonlinear shallow tides $M_4,$, $MS_4$ and long period $M_f,$, $M_{sf}$ were reproduced. It is found that implication of spatial variation of friction coefficient plays a very important role in reproduction of astronomical and shallow tides which are computed by iterative computation of nonlinear terms. Also it should be considered differently with respect to tidal periods. To understand the distribution of tidal asymmetry, amplitude ratio of $M_4/M_2$ and phase differences $2g(M_2)-g(M_4)$ were calculated. Tidal distortion ratio marks up to 0.2 on the west coast showing shallow coastal characteristics and somewhat wide range of ebb-dominances in front of Mokpo area are reproduced.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.4
• /
• pp.244-252
• /
• 2018

This paper provides quantification of the effects of the turbulence model and grid refinement on the analysis of tip vortex flows by using the RANS(Reynolds averaged Navier-Stokes) method. Numerical simulations of the tip vortex flows of the NACA $66_2$-415 elliptic hydrofoil were conducted, and two turbulence models for RANS closure were tested, i.e., the Realizable $k-{\varepsilon}$ model and the Reynolds stress transport model. Numerical results were compared with available experimental data, and it was shown that the data for the Reynolds stress transport model that were computed on the finest grid system had better agreement in reproducing the development and propagation of the tip vortex. The Realizable $k-{\varepsilon}$ model overestimated the turbulence level in the vortex core and showed a diffusive behavior of the tip vortex. The tip vortex cavitation on the hydrofoil and its trajectory also showed good agreement between the current numerical results that were obtained using the Reynolds stress transport model and the results observed in the experiment.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.207-219
• /
• 1993

Tsunamis along the east coast of Korea accompanying the 1983 East Sea central region earthquake is hindcasted with numerical models for tsunami propagation and inundation. Both linear and nonlinear models were used to compute propagation and elevation of tsunami waves on the coastal area of Korea. For the mesh refinement, grid system was divided into two sub-regions in Korean coastal area with final 10m grid resolution at interior area where serious inundation was observed. Calculated tsunami height distribution showed a general agreement with coastal observation. With interior detailed mesh system at mid-east coast region, the inundatin at the port of Imwon were qualitatively well reproduced by inundation and runup model.