• 대한수학회지
• /
• 제59권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.

• Structural Engineering and Mechanics
• /
• 제55권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.

• 설비공학논문집
• /
• 제10권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.

• 대한수학회지
• /
• 제58권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.

• 대한기계학회논문집B
• /
• 제22권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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 춘계학술대회논문집A
• /
• 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.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2005년도 학술발표회 논문집
• /
• 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.

• Nuclear Engineering and Technology
• /
• 제55권12호
• /
• pp.4395-4407
• /
• 2023

Bubble columns are widely encountered in several industries, especially in the field of nuclear safety. The Eulerian-Eulerian and the Eulerian-Lagrangian methods are commonly used to investigate bubble columns. Eulerian approaches require additional tasks such as strict volume conservation at the interface and a predefined well-structured grid. In contrast, the Lagrangian approach can be easily implemented. Hence, we introduce a fully Lagrangian approach for the simulation of bubble columns using the discrete bubble model (DBM) and moving particle semi-implicit (MPS) methods. Additionally, we propose a rigorous method to estimate the volume fraction accurately, and verified it through experimental data and analytical results. The MPS method was compared with the experimental data of Dambreak. The DBM was verified by analyzing the terminal velocity of a single bubble for each bubble size. It agreed with the analytical results for each of the four drag correlations. Additionally, the improved method for calculating the volume fraction showed agreement with the Ergun equation for the pressure drop in a packed bed. The implemented MPS-DBM was used to simulate the bubble column, and the results were compared with the experimental results. We demonstrated that the MPS-DBM was in quantitative agreement with the experimental data.

• 천문학회지
• /
• 제34권4호
• /
• pp.237-241
• /
• 2001

Our goal is to present a simple volume-of-fluid type interface-tracking algorithm to compressible two-phase flow in two space dimensions. The algorithm uses a uniform underlying Cartesian grid with some cells cut by the tracked interfaces into two subcells. A volume-moving procedure that consists of two basic steps: (1) the update of volume fractions in each grid cell at the end of the time step, and (2) the reconstruction of interfaces from discrete set of volume fractions, is employed to follow the dynamical behavior of the interface motion. As in the previous work with a surface-tracking procedure for general front tracking (LeVeque & Shyue 1995, 1996), a high resolution finite volume method is then applied on the resulting slightly nonuniform grid to update all the cell values, while the stability of the method is maintained by using a large time step wave propagation approach even in the presence of small cells and the use of a time step with respect to the uniform grid cells. A sample preliminary numerical result for an underwater explosion problem is shown to demonstrate the feasibility of the algorithm for practical problems.

• 한국전산유체공학회지
• /
• 제9권2호
• /
• pp.43-51
• /
• 2004

Interaction between fluid flow and thermal radiation has received considerable attention due to its numerous applications in engineering field. In this case the thermofluid properties of radiating fluid vary with the variation of temperature field caused by absorption and emission of radiant heat. To analyze the radiation heat transfer in radiating fluid, the simultaneous solution of the radiative transfer equation (RTE) and the fluid dynamics equations is required. This means that the numerical procedure used for the RTE must be computationally efficient to permit its inclusion in the other submodels, and must be compatible with the other transport equations. The finite volume method (FVM) and the discrete ordinates method (DOM) are usually employed to simulate radiation problems in generalized coordinates. These two representative methods are examined and compared, especially in view of the numerical integration of the radiation intensity over solid angle. The FVM shows better accuracy than the DOM owing to less constraints of the selection of control angle.