• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.93-98
• /
• 2001

An adaptive Cartesian grid method having the best elements of structured, unstructured, and Cartesian grids is developed to solve the steady two-dimensional Euler equations. The solver is based on a cell-centered finite-volume method with Roe's flux-difference splitting and implicit point Gauss-seidel time integration method. Calculations of several compressible flows are carried out to show the efficiency of the developed computer code. The results were generally in good agreements with existing data in the literature and the developed code has the good ability to capture important feature of the flows.

• Korean Journal of Mathematics
• /
• v.25 no.2
• /
• pp.279-301
• /
• 2017

An L(3, 2, 1)-labeling for the graph G = (V, E) is an assignment f of a label to each vertices of G such that ${\mid}f(u)-f({\upsilon}){\mid}{\geq}4-k$ when $dist(u,{\upsilon})=k{\leq}3$. The L(3, 2, 1)-labeling number, denoted by ${\lambda}_{3,2,1}(G)$, for G is the smallest number N such that there is an L(3, 2, 1)-labeling for G with span N. In this paper, we compute the L(3, 2, 1)-labeling number ${\lambda}_{3,2,1}(G)$ when G is a cylindrical grid, which is the cartesian product $P_m{\Box}C_n$ of the path and the cycle, when $m{\geq}4$ and $n{\geq}138$. Especially when n is a multiple of 4, or m = 4 and n is a multiple of 6, then we have ${\lambda}_{3,2,1}(G)=11$. Otherwise ${\lambda}_{3,2,1}(G)=12$.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.3
• /
• pp.260-266
• /
• 2011

For long slender offshore structures, such as cables and pipe lines, their interaction with surrounding fluid flow becomes an important issue for global design of ocean systems. We employ a long circular cylinder as a representative case of slender offshore structure. A flexibly mounted cylinder in cross-flow generates complex vortex shedding and results in oscillation of the structure. In this paper, flow behind a circular cylinder at Re=100 is simulated. The vortex shedding pattern and flow induced motion are examined in the cross flow configuration as well as with various yaw-angled configurations. The "Lock-in" phenomenon is also observed when reduced velocity is approximately 4.0. The MAC Grid system, which is the typical grid system for Cartesian mesh and pressure correction methods, are used for solving the incompressible Navier-Stokes equations. Predictor/Corrector method is applied for obtaining a non-linear response of structure at the flexibly mounted. The existance and motion of the body is represented by the immersed boundary technique.

• Kyungpook Mathematical Journal
• /
• v.59 no.3
• /
• pp.515-523
• /
• 2019

Consider a graph G of order n and maximum degree ${\Delta}$. Let $f:V(G){\rightarrow}\{0,1,{\cdots},{\lceil}{\frac{{\Delta}}{2}}{\rceil}+1\}$ be a function that labels the vertices of G. Let $B_0=\{v{\in}V(G):f(v)=0\}$. The function f is a strong Roman dominating function for G if every $v{\in}B_0$ has a neighbor w such that $f(w){\geq}1+{\lceil}{\frac{1}{2}}{\mid}N(w){\cap}B_0{\mid}{\rceil}$. In this paper, we study the bounds on strong Roman domination numbers of the Cartesian product $P_m{\square}P_k$ of paths $P_m$ and paths $P_k$. We compute the exact values for the strong Roman domination number of the Cartesian product $P_2{\square}P_k$ and $P_3{\square}P_k$. We also show that the strong Roman domination number of the Cartesian product $P_4{\square}P_k$ is between ${\lceil}{\frac{1}{3}}(8k-{\lfloor}{\frac{k}{8}}{\rfloor}+1){\rceil}$ and ${\lceil}{\frac{8k}{3}}{\rceil}$ for $k{\geq}8$, and that both bounds are sharp bounds.

• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.1-8
• /
• 2007

A code is developed using the hybrid Cartesian/immersed boundary method and it is applied to simulate flows around a three-dimensional deforming body. A new criterion is suggested to distribute the immersed boundary nodes based on edges crossing a body boundary. Velocities are reconstructed at the immersed boundary nodes using the interpolation along a local normal line to the boundary. Reconstruction of the pressure at the immersed boundary node is avoided using the hybrid staggered/non-staggered grid method. The developed code is validated through comparisons with other experimental and numerical results for the velocity profiles around a circular cylinder under the forced in-line oscillation and the pressure coefficient distribution on a sphere. The code is applied to simulate the flow fields around a plate whose tail is periodically flapping under a translation. The effects of the velocity and acceleration due to the deformation on the periodic shedding of pairs of tip vortices are investigated.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.170-175
• /
• 2009

Shock focusing is related with explosive release of shock wave energy on a narrow spot in a short duration of time triggering a spontaneous high pressure near the focal point. It is well known that reflection of planar incident shock wave from the metallic concave mirror such as ellipsoidal, paraboloidal or hemispherical cavities will focus on a focal point. We intend to improve the computational results using a wave propagation algorithm and to resolve the mushroom-like structure. For computation of the concave cavity flow, it is not easy to use a single-block mesh because of the many singular points in geometry and coordinates. We have employed a uniform Cartesian-grid method for the wave propagation algorithm.

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

In this paper the wave run-up around a circular column in regular waves is numerically calculated to investigate the applicability of the Modified Marker-Density (MMD) method to prediction of wave run-up around an offshore platform. The MMD method is one of the methods to define the highly nonlinear free surface. The governing equations are the Navier-Stokes equations and the continuity equation which are computed in Cartesian grid system. To validate incident waves generated by numerical simulation, those are compared with the solutions of the Stokes $5^{th}$ order wave theory. The wave run-up simulations are performed varying the steepness and period of incident waves as referred experimental data. The numerical results are compared to the experimental data and the results show good agreements.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.409-418
• /
• 2014

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.

• Journal of Ocean Engineering and Technology
• /
• v.26 no.6
• /
• pp.7-13
• /
• 2012

This paper presents a numerical study of the wave loads acting on offshore structures using a Cartesian-grid-based flow simulation method. Finite volume discretization with a volume-of-fluid (VOF) method is adopted to solve two-phase Navier-Stokes equations. Among the many variations of the VOF method, the CICSAM scheme is applied. The body boundary conditions are satisfied using a porosity function, and wave generation is carried out by using transient (wave or damping) zone approaches. In order to validate the present numerical method, three different basic offshore structures, including a sphere, Pinkster barge, and Wigley model, are numerically investigated. First, diffraction and radiation problems are solved using the present numerical method. The wave exciting and drift forces from the diffraction problems are compared with potential-based solutions. The added mass and wave damping forces from the radiation problems are also compared with the potential results. Next, the wave-induced motion responses of the structures are calculated and compared with the existing experimental data. The comparison results are fairly good, showing the validity of the present numerical method.

• Nuclear Engineering and Technology
• /
• v.55 no.2
• /
• pp.756-768
• /
• 2023

Neutron transport calculations are extremely challenging due to the high computational cost of large and complex problems. A multilevel octree grid algorithm (MLTG) of discrete ordinates method was developed to improve the modeling accuracy and simulation efficiency on 3-D Cartesian grids. The Balakovo-3 VVER-1000 neutron dosimetry benchmark is calculated to verify and validate this numerical technique. A simplified S₂ synthetic acceleration is used in the MLTG calculation method to improve the convergence of the source iterations. For the triangularly arranged fuel pins, we adopt a source projection algorithm to generate pin-by-pin source distributions of hexagonal assemblies. MLTG provides accurate geometric modeling and flexible fixed source description at a lower cost than traditional Cartesian grids. The total number of meshes is reduced to 1.9 million from the initial 9.5 million for the Balakovo-3 model. The numerical comparisons show that the MLTG results are in satisfactory agreement with the conventional S_N method and experimental data, within the root-mean-square errors of about 4% and 10%, respectively. Compared to uniform fine meshing, approximately 70% of the computational cost can be saved using the MLTG algorithm for the Balakovo-3 computational model.