• Wind and Structures
• /
• v.8 no.6
• /
• pp.407-426
• /
• 2005

A hybrid RANS/LES approach, based on the Limited Numerical Scales concept, is applied to the numerical simulation of the flow around a square cylinder. The key feature of this approach is a blending between two eddy-viscosities, one given by the $k-{\varepsilon}$ RANS model and the other by the Smagorinsky LES closure. A mixed finite-element/finite-volume formulation is used for the numerical discretization on unstructured grids. The results obtained with the hybrid approach are compared with those given by RANS and LES simulations for three different grid resolutions; comparisons with experimental data and numerical results in the literature are also provided. It is shown that, if the grid resolution is adequate for LES, the hybrid model recovers the LES accuracy. For coarser grid resolutions, the blending criterion appears to be effective to improve the accuracy of the results with respect to both LES and RANS simulations.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.129-134
• /
• 2011

In this paper, the confined single slot jet impingement is investigated numerically. Although the geometry of the jet impingement is simple, the flow pattern of the jet impingement is complex and the numerical results of the jet impingement is affected much by numerical methods. The first goal of this study is to analyze the effects of Reynolds models and numerical spatial discretization schemes on the results of heat transfer performance and the flaw characteristics and to select the best method. Various versions of the low Reynolds number k-epsilon turbulence models are compared. Using the selected numerical method, the flow field and heat transfer characteristics of confined single slot jet impingement on a moving plate are analyzed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.7
• /
• pp.916-923
• /
• 2000

Numerical study of three-dimensional turbulent flow in a forward curved centrifugal fan is presented. Standard $k-{\varepsilon}$ turbulence model and non-orthogonal curvilinear coordinates arc used to consider the turbulent flow field and complex geometry. Finite Volume approach is adopted for discretization scheme and structured grid system is used to help convergence. Multiblock grid system is used for flow field and divided into five domains that are inlet, outlet, impeller, tip clearance and scroll. It is assumed that the flow field is steady and incompressible. These numerical results are compared with the experimental data inside a rotor and at the fan outlet. Most important flow features are captured through this numerical approach. Finally details of flow field inside a fan are described and analyzed.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.1
• /
• pp.60-67
• /
• 2012

This study considers the numerical analysis on parametric roll for container ships. As a method of numerical simulation, an impulse-response-function approach is applied in time domain. A systematic study is carried out for the parametric roll of two container ships, particularly observing the sensitivity of computational results to some parameters which can affect the analysis of parametric roll. The parameters to be considered are metacentric height (GM), simulation time window, and the discretization of wave spectrum. Based on the result of parametric roll simulation, numerical sensitivity and uncertainty in computational analysis are discussed.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.22 no.1
• /
• pp.75-89
• /
• 2018

This paper presents the numerical valuation of the arithmetic Asian option by using the operator-splitting method (OSM). Since there is no closed-form solution for the arithmetic Asian option, finding a good numerical algorithm to value the arithmetic Asian option is important. In this paper, we focus on a two-dimensional PDE. The OSM is famous for dealing with plural-dimensional PDE using finite difference discretization. We provide a detailed numerical algorithm and compare results with MCS method to show the performance of the method.

• Korea-Australia Rheology Journal
• /
• v.20 no.1
• /
• pp.7-14
• /
• 2008

The objective of this work is to propose a procedure to simulate the flow in the LCM (Liquid Composites Molding) processes by finite difference discretization in a curvilinear coordinate system adapted to the shape of the saturated zone. The numerical results obtained are compared with experimental results obtained by an experimental device elaborated at our laboratory. It allows to realize linear and radial injections for different porosities and to observe the flow front kinetics. Numerical and experimental results are then compared with those of the literatures and excellent agreements are noticed. Finally, we suggest a concept of the capillary number to explain the variations of the permeability obtained for pressure values lower than 0.25 Bar.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.528-531
• /
• 2004

For a simple one-dimensional ignition problem a mathematical model is described to investigate the difficulties in numerical simulations. Some computation results are obtained and comparison is made with analytical solution. Discussions are made on topics such as 1) coordinate transformation, 2) gas-phase and solid-phase analysis; (divergence form of the governing system, a finite-volume discretization, implicit time integration, upwind split flux, spatial accuracy improvement are described. Mass, reagent mass, and energy conservations are solved.), and 3) method to determine quantities on the burning surface (matching). Results obtained for small values of the non-dimensional pressure show a steady-combustion and good agreement with the analytical solution. Numerical instability appeared for larger values of the pressure, discussion on the cause of the problem is made. This effort is a part of a study of flame spread phenomena on solid propellant surface.

• Nuclear Engineering and Technology
• /
• v.54 no.4
• /
• pp.1195-1205
• /
• 2022

Recent developments on spectral diffusion algorithms, i.e., algorithms which exploit the projection of the solution on the eigenfunctions of the Laplacian operator, demonstrated their effective applicability in fast transient conditions. Nevertheless, the numerical error introduced by these algorithms, together with the uncertainties associated with model parameters, may impact the reliability of the predictions on short-lived volatile fission product release from nuclear fuel. In this work, we provide an upper bound on the numerical error introduced by the presented spectral diffusion algorithm, in both constant and time-varying conditions, depending on the number of modes and on the time discretization. The definition of this upper bound allows introducing a methodology to a priori bound the numerical error on short-lived volatile fission product retention.

• Journal of applied mathematics & informatics
• /
• v.41 no.2
• /
• pp.295-309
• /
• 2023

The Richards' equation attracts the attention of several scientific researchers due to its importance in the hydrogeology field especially porous soil. This work presents a numerical method to solve the two dimensional Richards' equation. The pressure form and the mixed form of Richards' equation are solved numerically using a bivariate diamond finite volumes scheme. Euler explicit scheme is used for the time discretization. Different test cases are done to validate the accuracy and the efficiency of our numerical model and to compare the possible numerical strategies. We started with a first simple test case of Richards' pressure form where the hydraulic capacity and the hydraulic conductivity are taken constant and then a second test case where the hydrodynamics parameters are linear variables. Finally, a third test case where the soil parameters are taken according the Van Gunchten empirical model is presented.

• Journal of the Korean Society of Safety
• /
• v.21 no.2 s.74
• /
• pp.143-149
• /
• 2006

This paper deals with the space-time finite element analysis of two-dimensional vibration problem with a single variable. The method of space-time finite elements enables the simpler solution than the usual finite element analysis with discretization in space only. We present a discretization technique in which finite element approximations are used in time and space simultaneously for a relatively large time period. The weighted residual process is used to formulate a finite element method for a space-time domain. A stability problem is described and some investigations for chosen type of rectangular space-time finite elements are carried out. Instability is caused by a too large time step of successive time steps in the traditional time-dependent problems. It has been shown that the numerical stability of time-stepping on the larger time steps is quite good. The unstructured space-time finite element not only overcomes the shortcomings of the stability in the traditional numerical methods, but it is also endowed with the features of an effective computational technique. Some numerical examples have been presented to illustrate the efficiency of the described method.