• Journal of the Korean Mathematical Society
• /
• v.59 no.3
• /
• pp.519-548
• /
• 2022

In this paper, a stabilized-penalized collocated finite volume (SPCFV) scheme is developed and studied for the stationary generalized Navier-Stokes equations with mixed Dirichlet-traction boundary conditions modelling an incompressible biological fluid flow. This method is based on the lowest order approximation (piecewise constants) for both velocity and pressure unknowns. The stabilization-penalization is performed by adding discrete pressure terms to the approximate formulation. These simultaneously involve discrete jump pressures through the interior volume-boundaries and discrete pressures of volumes on the domain boundary. Stability, existence and uniqueness of discrete solutions are established. Moreover, a convergence analysis of the nonlinear solver is also provided. Numerical results from model tests are performed to demonstrate the stability, optimal convergence in the usual L² and discrete H¹ norms as well as robustness of the proposed scheme with respect to the choice of the given traction vector.

• Journal of the Chungcheong Mathematical Society
• /
• v.36 no.4
• /
• pp.195-213
• /
• 2023

Discrete volumes of poset polytopes for a variant of up-down posets introduced in [4] were studied. We obtained the generating functions for the discrete volumes of poset polytopes for a variant of up-down poset using the characteristic matrices.

• Structural Engineering and Mechanics
• /
• v.55 no.2
• /
• pp.315-334
• /
• 2015

This study's primary aim is to check the existence of a representative volume element for granular materials and determine the link between the properties (responses) of macro structures and the size of the discrete particle assembly used to represent a constitutive relation in a two-scale model. In our two-scale method the boundary value problem on the macro level was solved using finite element method, based on the Cosserat continuum; the macro stresses and modulus were obtained using a solution of discrete particle assemblies at certain element integration points. Meanwhile, discrete particle assemblies were solved using discrete element method under boundary conditions provided by the macro deformation. Our investigations focused largely on the size effects of the discrete particle assembly and the radius of the particle on macro properties, such as deformation stiffness, bearing capacity and the residual strength of the granular structure. According to the numerical results, we suggest fitting formulas linking the values of different macro properties (responses) and size of discrete particle assemblies. In addition, this study also concerns the configuration and displacement fluctuation of discrete particle assemblies on the micro level, accompanied with the evolution of bearing capacity and deformation on the macro level.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.399-404
• /
• 2000

Functionally graded materials(FGMs) involve dual-phase graded layers in which two different constituents are mixed continuously and functionally according to a given volume fraction. For the analysis of their thermo-mechanical response, conventional homogenized methods have been widely employed in order to estimate equivalent material properties of the graded layer. However, such overall estimations are insufficient to accurately predict the local behavior. In this paper, we compare the thermo-elastic behaviors predicted by several overall material-property estimation techniques with those obtained by discrete analysis models utilizing the finite element method, for various volume fractions and loading conditions.

• Journal of the Korean Mathematical Society
• /
• v.58 no.3
• /
• pp.553-569
• /
• 2021

In this paper, we present and analyze a fully discrete numerical method for solving the time-fractional diffusion wave equation: ∂^β_tu - div(a∇u) = f, 1 < β < 2. We first construct a difference formula to approximate ∂^β_tu by using an interpolation of derivative type. The truncation error of this formula is of O(△t^2+δ-β)-order if function u(t) ∈ C^2,δ[0, T] where 0 ≤ δ ≤ 1 is the Hölder continuity index. This error order can come up to O(△t^3-β) if u(t) ∈ C³ [0, T]. Then, in combinination with the linear finite volume discretization on spatial domain, we give a fully discrete scheme for the fractional wave equation. We prove that the fully discrete scheme is unconditionally stable and the discrete solution admits the optimal error estimates in the H¹-norm and L₂-norm, respectively. Numerical examples are provided to verify the effectiveness of the proposed numerical method.

• Nuclear Engineering and Technology
• /
• v.44 no.6
• /
• pp.639-652
• /
• 2012

A three-dimensional discrete fracture network model was developed in order to simulate the hydraulic characteristics of a granitic rock mass at Korea Atomic Energy Research Institute (KAERI) Underground Research Tunnel (KURT). The model used a three-dimensional discrete fracture network (DFN), assuming a correlation between the length and aperture of the fractures, and a trapezoid flow path in the fractures. These assumptions that previous studies have not considered could make the developed model more practical and reasonable. The geologic and hydraulic data of the fractures were obtained in the rock mass at the KURT. Then, these data were applied to the developed fracture discrete network model. The model was applied in estimating the representative elementary volume (REV), the equivalent hydraulic conductivity tensors, and the amount of groundwater inflow into the tunnel. The developed discrete fracture network model can determine the REV size for the rock mass with respect to the hydraulic behavior and estimate the groundwater flow into the tunnel at the KURT. Therefore, the assumptions that the fracture length is correlated to the fracture aperture and the flow in a fracture occurs in a trapezoid shape appear to be effective in the DFN analysis used to estimate the hydraulic behavior of the fractured rock mass.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.10 no.1
• /
• pp.95-107
• /
• 1998

Natural convection and radiation heat transfer in a square enclosure containing absorbing, emitting, and isotopically scattering(participating) media is studied numerically using the finite volume method. Various numerical methods are employed to analyze the radiative heat transfer. However, it is very difficult to choose the proper method. In present study, a finite volume method(FVM) and a discrete ordinates method(DOM) are compared in rectangular enclosure. The SIMPLER algorithm is used to solve the momentum and energy equations. Thermal and flow characteristics are investigated according to the variation of radiation parameters such as optical thickness and scattering albedo. The result shows that the accuracy and the computing time of FVM are better than those of DOM in regular geometry.

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

The numerical simulation of a two-phase flow In a discrete fracture network model is presented in this paper, The purpose of this work is to consider density-driven flows induced by the density difference between hot autochthonous heavy brines and injected cold water. Mechanical consequences of high pressure waves on the fracture permeability and heat exchanges between fluids and rock matrix are neglected in this study. The finite volume method is employed to discretize spatially and the system is solved by using an IMPES(Implicit Pressure-Explicit Saturation) scheme. In order to solve the strong non-linearity of the system, the Newton-Raphson algorithm is used. The well-known Buckeley-Leverett problem is adapted to validate results calculated from the model and a relatively good agreement is obtained.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.2
• /
• pp.213-220
• /
• 1998

The conventional discrete-ordinates method (DOM) is modified and developed for the analysis of two-dimensional axisymmetric cylindrical enclosure with curved wall. The objective of the present work is to extend the capability of the conventional DOM into a general axisymmetric geometry like nozzle-shaped enclosure, by adopting the arbitrary control angle as was done in the finite-volume method (FVM), while keeping the same two-dimensional solution procedure as in the conventional DOM. The present method is validated by applying it to three different benchmark problems of axisymmetric enclosure containing absorbing, emitting and scattering medium. Results presented in this work not only support the solution accuracy, but also moderate efficiency in the numerical calculation of axisymmetric radiation problem.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.7
• /
• pp.593-602
• /
• 2016

Numerical simulations on the close-coupled supersonic gas atomizer for metallic powder production were performed in this study. A proper turbulence model was chosen and then VOF(Volume of Fluid) and DPM(Discrete Phase Model) methods were sequentially applied for the simulations of primary and secondary break-up processes of liquid metal. Diameters of parent droplets were calculated by analyzing Level-Set function contour from the VOF result. Finally, the distribution of particle diameter was obtained from the DPM result at exit of the computational domain.