• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제25권3호
• /
• pp.66-81
• /
• 2021

In this paper, we review the lowest order staggered discontinuous Galerkin methods on polygonal meshes in 2D. The proposed method offers many desirable features including easy implementation, geometrical flexibility, robustness with respect to mesh distortion and low degrees of freedom. Discrete function spaces for locally H¹ and H(div) spaces are considered. We introduce special properties of a sub-mesh from a given star-shaped polygonal mesh which can be utilized in the construction of discrete spaces and implementation of the staggered discontinuous Galerkin method. For demonstration purposes, we consider the lowest case for the Poisson equation. We emphasize its efficient computational implementation using only geometrical properties of the underlying mesh.

• Journal of Mechanical Science and Technology
• /
• 제14권9호
• /
• pp.985-996
• /
• 2000

The issue of quench is related to safety operation of large-scale superconducting magnet system fabricated by cable-in-conduit conductor. A numerical method is presented to simulate the thermal hydraulic quench characteristics in the superconducting Tokamak magnet system, One-dimensional fluid dynamic equations for supercritical helium and the equation of heat conduction for the conduit are used to describe the thermal hydraulic characteristics in the cable-in-conduit conductor. The high heat transfer approximation between supercritical helium and superconducting strands is taken into account due to strong heating induced flow of supercritical helium. The fully implicit time integration of upwind scheme for finite volume method is utilized to discretize the equations on the staggered mesh. The scheme of a new adaptive mesh is proposed for the moving boundary problem and the time term is discretized by the-implicit scheme. It remarkably reduces the CPU time by local linearization of coefficient and the compressible storage of the large sparse matrix of discretized equations. The discretized equations are solved by the IMSL. The numerical implement is discussed in detail. The validation of this method is demonstrated by comparison of the numerical results with those of the SARUMAN and the QUENCHER and experimental measurements.

• Interaction and multiscale mechanics
• /
• 제5권3호
• /
• pp.287-318
• /
• 2012

We apply a partitioned-solution (iterative-staggered) coupling method based on a fixed Eulerian mesh with the level set function to a large-deformation fluid-structure interaction (FSI) problem where a large-deformable thin structure moves in a high-speed flow field, as an airbag does during deployment. This method combines advanced fluid and structure solvers-specifically, the constrained interpolation profile finite element method (CIP-FEM) for fluid Eulerian mesh and large-deformable structural elements for Lagrangian structural mesh. We express the large-deformable interface as a zero isosurface by the level set function, and introduce virtual nodes with level sets and structural normal velocities to generate the level set function according to the large-deformable interfacial geometry and enforce the kinematic condition at the interface. The virtual nodes are located in the direction normal to the structural mesh. It is confirmed that application of the method to unfolded airbag deployment simulation shows the adequacy of the method.

• Nuclear Engineering and Technology
• /
• 제16권2호
• /
• pp.47-57
• /
• 1984

2상유동을 해석하기 위한 3차원 코드인 THERMIT-6S의 미분 방정식을 세우기 위해, 수학적으로 정확하게 유도된 시간과 공간에 대해 평균한 보존 방정식을 단순화했다. 미분 방정식을 불연속화(discretization)하여 THERMIT-6S의 차분방정식을 얻는다. First-order spatial scheme, donor cell method, 그리고, staggered mesh layout을 써서 공간에 대한 불연속화를 한다. 그리고 시간에 대한 불연속화는 first-order semi-implicit scheme으로써, sonic terms와 국부적인 전달 현상에 관계되는 항들은 implicit하게 그리고 대류 전달 항들은 explicit하게 취급한다. 이렇게 얻어진 방정식들은Newton-Raphson 방법으로 선형화된다. 축소된 압력 방정식을 만들기 위해 모든 변수들이 mesh cells사이에서 단지 압력 변수를 통해서만 결부되도록, 선형화된 방정식들을 처리한다. OPERA-15 실험을 수치해석적으로 모의실험하여 본 결과, THERMIT-6S가 flow coastdown, 역류, 유체진동(flow oscillation) 등을 포함하고, sodium boiling 후의 원자로내의 변화를 예측하는데 매우 유효하다는 것이 밝혀졌다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.68-74
• /
• 2011

Recently, thanks to the advanced computational power and numerical methods, it is made possible to analyze the flow around moving bodies using computational fluid dynamics techniques. In those simulations, moving mesh techniques should be able to represent both the body motion and boundary deformation, which are frequently encountered in fluid-structure interaction and/or six degree-of-freedom problems. In the present study, the staggered loosely coupling algorithm was used for fluid-structure interaction and the Laplacian operator based technique was used for moving mesh. For the verification of the developed computational method, the flow around a two-dimensional cylinder was simulated and analyzed.

• Journal of Advanced Marine Engineering and Technology
• /
• 제32권1호
• /
• pp.82-93
• /
• 2008

The flowfield behind two cylinders and flow-induced noise generated from the cylinders in various arrangement are numerically investigated based on the finite difference lattice Boltzmann model with 21 velocity bits. which is introduced a flexible specific heat ${\gamma}$ to simulate diatomic gases like air. In an isolated cylinder with two type of mesh. some flow parameters such as Strouhal number $S_t$ and acoustic pressure ${\Delta}p$ simulated from the solution are given and quantitatively compared with those provided the previous works. The effects of the center-to-center pitch ratio $L_{cc}/d=2.0$ in staggered circular cylinders as shown in Fig. 1 and angles of incidence ${\alpha}=30^{\circ}(T_{cc}/d=0.5)$, $45^{\circ}(T_{cc}/d =0.707)$ and $60^{\circ}\;(T_{cc}/d=0.866)$, respectively, are studied. Our analysis focuses on the small-scale instabilities of vortex shedding, which occurs in staggered arrangement. With the results of drag $C_d$ and lift $C_l$ coefficients and vorticity contours. the mechanisms of the interference phenomenon and its interaction with the two-dimensional vortical structures are present in the flowfields under $Re\;{\le}\;200$. The results show that we successively capture very small pressure fluctuations, with the same frequency of vortex shedding, much smaller than the whole pressure fluctuation around pairs of circular cylinders. The upstream cylinder behaves like an isolated single cylinder, while the downstream one experiences wake-induced flutter. It is expected that, therefore, the relative position of the downstream cylinder has significant effects on the flow-induce noise, hydrodynamic force and vortex shedding characteristics of the cylinders.

• 한국전산유체공학회지
• /
• 제16권1호
• /
• pp.73-82
• /
• 2011

In this paper, the vortex-induced vibration of circular cylinders is studied using the immersed boundary method on the Cartesian mesh. The Reynolds numbers considered is from 100 to 200. Using the configuration of tendemly arranged multiple circular cylinders, the vortex shedding behind of the cylinders and their flow-induced motion are investigated. The staggered MAC grid arrangement, which is the typical grid system for the incompressible flow on the Cartesian meshes, is utilized. Pressure correction method is applied for solving the divergence-free incompressible velocity field. The body motion is described by immersed boundary technique that has advantages for moving object on the fixed computational domain. It is also discussed for the computational noise in hydrodynamic forces when body motion is represented by the immersed boundary method. The Predictor/Corrector method is used for simulating the nonlinear response of the elastically mounted cylinder excited by vortex-shedding.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2007년도 추계 학술대회논문집
• /
• pp.243-248
• /
• 2007

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed. A two-fluid three-field model was adopted for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. The hydrodynamic solver is for the 3D component of a nuclear system code and the component-scale analysis tools for transient two-phase flows. The finite volume method and unstructured grid are adopted, which are useful for the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been adapted to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing problems well.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.326-334
• /
• 2010

In this study, we provide a comprehensive review of the CIP(Constrained Interpolation Pro file/Cubic Interpolated Propagation) method with a pressure-based algorithm that is known as a general numerical solver for soled liquid, gas and plasmas. And also we introduce a body-fitted grid system(Soroban grid) for computation of strongly nonlinear marine hydrodynamic problems such as slamming water on deck, wave impact by green water. This grid system can keep the third-order accuracy in time and space with the help of the CIP method. The grid system consists of the straight lines and grid points. In the 2-dimensional grid case, each grid points moving in these lines like abacus - Soroban in Japanese. The length of each line can be different and the number of grid points in each line can be different. Mesh generation and searching of upstream departure point are very simple and possible to mesh-free treatment. To optimize computation of free-surface and multi-fluid flows, We adopt the C-CUP method. In most of the earlier computations, the C-CUP method was used with a staggered-grid approach. Here, because of the mesh free nature of the Soroban grid, we use the C-CUP method with a collocated-grid approach.

• 에너지공학
• /
• 제17권4호
• /
• pp.218-226
• /
• 2008

가압 경수로의 주요 기기에서 발생할 수 있는 과도 2상 유동(Two-phase flow) 현상에 대한 해석을 수행하기 위해 원자로 기기 열수력 해석 코드를 개발 중에 있다. 개발중인 기기 열수력 해석 코드는 지배 방정식으로 Two-phase, three-field model을 사용하고 있으며, 복잡한 기하학적 형상의 원자로 기기를 모사하기 위해 비정렬 격자계(Unstructured grid)를 활용하고 있다. 수치해석 기법으로는, 원자로 계통 해석코드 REIAP5가 사용 중이며 대부분의 원자로 내 2상 유동 조건에서 안정적이며 정확하다고 알려진 Semi-implicit 방법을 적용하였다. 그러나 기존의 Semi-implicit 방법은 1차원, 엇갈림격자(Staggered grid)에 대해 개발되었기 때문에 이를 다차원, 비정렬, 비엇갈림 격자(Non-staggered grid)에 적용하기 위해 기존의 Semi-implicit 방법을 수정하였다. 본 논문에서는 수정된 Semi-implicit 방법을 소개하고 이를 이용해 수행한 예비 계산결과를 수록하였다.