• Journal of the Korea Computer Graphics Society
• /
• v.18 no.1
• /
• pp.29-34
• /
• 2012

High-quality tetrahedral meshes are crucial for FEM-based simulation of large elasto-plastic deformation and tetrahedral-mesh-based simulation of fluid flow. This paper proposes a volume meshing method that exploits an octree-based adaptive signed distance field to fill the inside of a polygonal object with tetrahedra, of which dihedral angles are good. The suggested method utilizes an octree structure to reduce the total number of tetrahedra by space-efficiently filling an object with graded tetrahedra. To obtain a high-quality mesh with good dihedral angles, we restrict the octree in such a way that any pair of neighboring cells only differs by one level. In octree-based tetrahedral meshing, the signed distance computation of a point to the surface of a given object is a very important and frequently-called operation. To accelerate this operation, we develop a method that computes a signed distance field directly on the vertices of the octree cells while constructing the octree using a top-down approach. This is the main focus of the paper. The suggested tetrahedral meshing method is fast, stable and easy to implement.

• Journal of computational fluids engineering
• /
• v.9 no.4
• /
• pp.55-63
• /
• 2004

The three-dimensional Stokes flow within a staggered screw channel is obtained by using a finite volume method. The geometry is intended to mimic the single screw extruder having staggered arrangement of flights. The flow solution is then subjected to the analysis of the stirring performance. In the analysis of the stirring performance, the stretching-mapping method developed by the author is employed for calculating the materials' stretching exponents, which are to be used in quantification of the mixing effect. The numerical results Indicate that the staggered geometry gives indeed far much better stirring-performance than the standard (nonstaggered) flight geometry. It was also shown that care must be given to the selection of the basis planes for evaluating the local stretching rate, and it turns out that the best method (H-method) has its basis plane just on the half way between the past and future evolution of fluid particles subjected to the defromation. In evaluating the stretching exponent, the expansion ratio must be considered which is one of the characteristic differences of the actual three-dimensional flows from the two-dimensionmal counterparts. The larger axial pressure-difference causes in general the smaller stirring performance while the flow rate is increased. The smaller channel length also increases the stirring performance.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.568-576
• /
• 2005

A 3-D finite element model combined with a volume tracking method is presented in this work to simulate the mold filling for casting processes. Especially, the analysis involves an adaptive grid method that is created under a criterion of element categorization of filling states and locations in the total region at each time step. By using an adaptive grid wherein the elements, finer than those in internal and external regions, are distributed at the surface region through refinement and coarsening procedures, a more efficient analysis of transient fluid flow with free surface is achieved. Adaptive grid based on VOF method is developed in tetrahedral element system. Through a 3-D analysis of the benchmark test of the casting process, the efficiency of the proposed adaptive grid method is verified. Developed FE code is applied to a typical industrial part of the casting process such as aluminum road wheel.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.86-94
• /
• 2010

A thermal-hydraulic code, named CUPID, has been developed for the analysis of transient two-phase flows in nuclear reactor components. A two-fluid three-field model was used for steam-water two-phase flows. To obtain numerical solutions, the finite volume method was applied over unstructured cell-centered meshes. In steam-water two-phase flows, a phase change, i.e., evaporation or condensation, results in a great change in the flow field because of substantial density difference between liquid and vapor phases. Thus, two-phase flows are very sensitive to the local pressure distribution that determines the phase change. This in turn puts emphasis on the accurate evaluation of local pressure gradient. This paper presents a new reconstruction method to evaluate the pressure gradient at cell centers on unstructured meshes. The results of the new scheme for a simple test function, a gravity-driven cavity, and a wall boiling two-phase flow are compared with those of the previous schemes in the CUPID code.

• Journal of the Korea Computer Graphics Society
• /
• v.15 no.1
• /
• pp.1-6
• /
• 2009

This paper presents a new CIP method for unstructured mesh to reduce the numerical dissipation. To reflect precise physical characteristics, CIP method updates both the physical quantity and the derivative information. The proposed method uses the Finite Volume Method(FVM) to solve the non-advection term of CIP equation. And we performed several experiments to improve the accuracy of third-order interpolation. Our result shows that our algorithm has less numerical dissipation than that of linear advection solver.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.7 no.3
• /
• pp.580-594
• /
• 2015

Liquid sloshing in two-dimensional (2-D) and three-dimensional (3-D) rectangular tanks is simulated by using a level set method based on the finite volume method. In order to examine the effect of natural frequency modes on liquid sloshing, we considered a wide range of frequency ratios ($0.5{\leq}fr{\leq}3.2$). The frequency ratio is defined by the ratio of the excitation frequency to the natural frequency of the fluid, and covers natural frequency modes from 1 to 5. When fr = 1, which corresponds to the first mode of the natural frequency, strong liquid sloshing reveals roof impact, and significant forces are generated by the liquid in the tank. The liquid flows are mainly unidirectional. Thus, the strong bulk motion of the fluid contributes to a higher elevation of the free surface. However, at fr = 2, the sloshing is considerably suppressed, resulting in a calm wave with relatively lower elevation of the free surface, since the waves undergo destructive interference. At fr = 2, the lower peak of the free surface elevation occurs. At higher modes of $fr_3$, $fr_4$, and $fr_5$, the free surface reveals irregular deformation with nonlinear waves in every case. However, the deformation of the free surface becomes weaker at higher natural frequency modes. Finally, 3-D simulations confirm our 2-D results.

• Agricultural and Biosystems Engineering
• /
• v.3 no.2
• /
• pp.65-72
• /
• 2002

Temperature is the most influential environment parameter which affects the quality change of agricultural products in cold storage. Therefore, it is essential to keep the uniform temperature distribution in the storage room. This study was performed to analyze the air movement and temperature distribution in the forced recirculating cold storage facility and to simulate optimum storage method of green groceries using 3-D CFD(three dimensional computational fluid dynamics) computer simulation which applied the standard $textsc{k}$-$\varepsilon$ turbulence model and FVM(finite volume method). The simulation was validated by the experimental results for onion storage and the simulation model was used to simulate the temperature and velocity distribution in the storage room with reference to the change of storage method such as location of storage, no stores, bulk storage, and pallet storage. In case of no stores, internal airflow was circulated without stagnation and consequently air movement and temperature distribution were uniform. In case of bulk storage, air movement was stagnated so much and temperature distribution of onion was not uniform. Furthermore, the inner temperature of onion roses more than the initial temperature of storage. In case of pallet storage, air movement and temperature distribution of onion were so uniform that the danger of quality change was decreased.

• Ocean Systems Engineering
• /
• v.12 no.3
• /
• pp.359-384
• /
• 2022

We develop in this work a new well-balanced preserving-positivity path-conservative central-upwind scheme for Saint-Venant-Exner (SVE) model. The SVE system (SVEs) under some considerations, is a nonconservative hyperbolic system of nonlinear partial differential equations. This model is widely used in coastal engineering to simulate the interaction of fluid flow with sediment beds. It is well known that SVEs requires a robust treatment of nonconservative terms. Some efficient numerical schemes have been proposed to overcome the difficulties related to these terms. However, the main drawbacks of these schemes are what follows: (i) Lack of robustness, (ii) Generation of non-physical diffusions, (iii) Presence of instabilities within numerical solutions. This collection of drawbacks weakens the efficiency of most numerical methods proposed in the literature. To overcome these drawbacks a reformulation of the central-upwind scheme for SVEs (CU-SVEs for short) in a path-conservative version is presented in this work. We first develop a finite-volume method of the first order and then extend it to the second order via the averaging essentially non oscillatory (AENO) framework. Our numerical approach is shown to be well-balanced positivity-preserving and shock-capturing. The resulting scheme could be seen as a predictor-corrector method. The accuracy and robustness of the proposed scheme are assessed through a carefully selected suite of tests.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.4
• /
• pp.525-532
• /
• 2010

In this paper, wave breakers which occur in two dimensional coasts are simulated using a SPH(Smoothed Particle Hydrodynamics) method which represents the movement of fluidic physical volume with particles. As continuative fluid is approximated to the particles, the simulations are performed using fully Lagrangian method without any grid system. Two-dimensional Navier-Stokes equations and continuity equation are used for the numerical simulations. To generate incident waves, a piston type wavemaker is employed. The accuracy of the wave which is numerically generated by the wavemaker is verified by comparing with analytical results. The computations are carried out with various wave heights and slopes. The wave patterns generated through the numerical simulations are compared with several existing experimental and computational results. Agreement between the experimental data and the computation results is comparatively good. Also, the breaker depth index and the breaker height index from the present calculations are compared with the existing experimental results, and the tendency is very similar.

• Journal of Ocean Engineering and Technology
• /
• v.26 no.1
• /
• pp.9-16
• /
• 2012

This paper presents the numerical results of a simulation of the free surface flow around a blunt bow ship model and focuses on the validation of the proposed method with a brief investigation of the relation between the resistance and free surface behavior. A finite volume method based on the Reynolds Averaged Navier-Stokes (RANS) approach is used to solve the governing flow equations, where the free surface, including wave breaking,is captured by using a two-phase Level-Set (LS) method. For turbulence closure, a two equation k-${\varepsilon}$ model with the standard wall function technique is used. Finally, the numerical results are compared with the available experimental data, showing good agreement.