• Journal of Ocean Engineering and Technology
• /
• v.29 no.2
• /
• pp.154-162
• /
• 2015

In this study, experimental and numerical methods were applied to observe sloshing impact phenomena. A two-dimensional rectangular tank filled with water and air was considered with a specific excitation condition that induced a hydrodynamic impact without an air pocket at the top corner of the tank. High-speed cameras and a pressure measurement system were synchronized, and a particle image velocimetry (PIV) technique was applied to measure the velocity field and corresponding pressure. The experimental condition was implemented in a numerical computation to solve incompressible two-phase flows using a Cartesian-grid method. The discretized solution was obtained using the finite difference and constraint-interpolation-profile (CIP) methods, which adopt a fractional step scheme for coupling the pressure and velocity. The tangent of the hyperbola for interface capturing (THINC) scheme was used with the weighed line interface calculation (WLIC) method to capture the interface between the air and water. The calculated impact pressures and velocity fields were compared with experimental data, and the relationship between the local velocity and pressure was investigated based on the computational results.

• Structural Engineering and Mechanics
• /
• v.56 no.4
• /
• pp.589-603
• /
• 2015

The mixed-mode stress intensity factors of 2-D angled cracks are evaluated by Petrov-Galerkin natural element (PG-NE) method in which Voronoi polygon-based Laplace interpolation functions and CS-FE basis functions are used for the trial and test functions respectively. The interaction integral is implemented in a frame of PG-NE method in which the weighting function defined over a crack-tip integral domain is interpolated by Laplace interpolation functions. Two Cartesian coordinate systems are employed and the displacement, strains and stresses which are solved in the grid-oriented coordinate system are transformed to the other coordinate system aligned to the angled crack. The present method is validated through the numerical experiments with the angled edge and center cracks, and the numerical accuracy is examined with respect to the grid density, crack length and angle. Also, the stress intensity factors obtained by the present method are compared with other numerical methods and the exact solution. It is observed from the numerical results that the present method successfully and accurately evaluates the mixed-mode stress intensity factors of 2-D angled cracks for various crack lengths and crack angles.

• Proceedings of the KIEE Conference
• /
• 2003.11b
• /
• pp.23-26
• /
• 2003

The vector field histogram(VFH) uses a two-dimensional Cartesian histogram grid as a world model. The VFH method subsequently employs a two-stage data-reduction process in order to compute the desired control commands for the vehicle. In the first stage the histogram grid is reduced to a one dimensional polar histogram that is constructed around the robot's momentary location. Each sector in the polar histogram contains a value representing the polar obstacle density in that direction. In the second stage, the algorithm selects the most suitable sector from among all polar histogram sectors with a low polar obstacle density, and the steering of the robot is aligned with that direction. We applied this algorithm to our four-legged robot.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2003.06a
• /
• pp.1076-1079
• /
• 2003

The vector field histogram(VFH) uses a two-dimensional Cartesian histogram grid as a world model. The VFH method subsequently employs a two-stage data-reduction process in order to compute the desired control commands for the vehicle. In the first stage the histogram grid is reduced to a one dimensional polar histogram that is constructed around the robot's momentary location. Each sector in the polar histogram contains a value representing the polar obstacle density in that direction. In the second stage, the algorithm selects the most suitable sector from among all polar histogram sectors with a low polar obstacle density, and the steering of the robot is aligned with that direction. We applied this algorithm to our simulation program and tested..

• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.49-57
• /
• 2013

Analysis of fluid-structure interaction for two nearby underwater vehicles immersed in the sea is quite challenging because simulation of flow around them is very difficult due to the complexity of underwater vehicle shapes. The conventional approach using body-fitted or unstructured grids demands much time in dynamic grid generation, and yields slow convergence of solution. Since an analysis of fluid-structure interaction must be based on accurate simulation results, a more efficient way of simulating flow around underwater vehicles, without sacrificing accuracy, is desirable. An immersed boundary method facilitates implementation of complicated underwater-vehicle shapes on a Cartesian grid system. An LES modeling is also incorporated to resolve turbulent eddies. In this paper, we will demonstrate the effectiveness of the immersed boundary method we adopted, by presenting the simulation results on the flow around a modeled high-speed underwater vehicle interacting with a modeled low-speed one.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.397-403
• /
• 2010

The Immersed boundary method(IBM) is one of CFD techniques which can simulate flow field around complex objectives using simple Cartesian grid system. In the previous studies the IBM has mostly been implemented to fractional step method based Navier-Stokes solvers. In these cases, pressure buildup near IB was found to occur when linear interpolation and stadard mass conservation is used and the interpolation scheme became complicated when higher order of interpolation is adopted. In this study, we implement the IBM to an incompressible Navier-Stokes solver which uses SIMPLE algorithm. Bi-linear and quadratic interpolation equations were formulated by using only geometric information of boundary to reconstruct velocities near IB. Flow around 2D circular cylinder at Re=40 and 100 was solved by using these formulations. It was found that the pressure buildup was not observed even when the bi-linear interpolation was adopted. The use of quadratic interpolation made the predicted aerodynamic forces in good agreement with those of previous studies.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.1
• /
• pp.1-10
• /
• 2010

In the present paper, a numerical method based on the constraint interpolation profile (CIP) method is applied for simulating two-dimensional violent sloshing problems. The free surface boundary value problem is considered as a multiphase problem which includes water and air. A stationary Cartesian grid system is adopted, and an interface capturing method is used to trace the shape of free surface profile. The CIP combined unified procedure (CCUP) scheme is applied for flow solver, and the tangent of hyperbola for interface capturing (THINC) scheme is used for interface capturing. Numerical simulations have been carried out for partially-filled 2D tanks under forced sway and roll motions at various filling depths and frequencies. The computational results are compared with experiments and/or the other numerical results to validate the present numerical method.

• The Bulletin of The Korean Astronomical Society
• /
• v.36 no.1
• /
• pp.82.1-82.1
• /
• 2011

We present a new code for solving non-LTE radiative transfer problems in a general grid (RIG). RIG develops from RATRAN code (Hogerheijde & van der Tak 2000) using the Accelerated Monte-Carlo method, and it can cope with line overlap effect among multiple molecular and atomic species. In this algorithm we make grids in arbitrary coordinates adequate to the problem, but, on the other hand, photons propagate in the Cartesian coordinates. For spherical, cylindrical and other well defined coordinate, the problem of tracing photon's path reduces to solving simple quadratic equations. For example, the outflow in the star formation have high dynamic range in scales from a few AU to ~ 0.1 pc and have also cylindrical symmetry. So, we have used (r, ${\alpha}$) coordinate system, where r is the distance from the origin and ${\alpha}$ is z/ R2 in the cylindrical coordinate of (R,z). The (r, ${\alpha}$) coordinate realizes the density - power function of r - and temperature distributions of the problems with smaller numbers of grid than the cylindrical coordinate does, and the former consumes less time to solve the problems than the latter.

• Journal of computational fluids engineering
• /
• v.20 no.4
• /
• pp.51-57
• /
• 2015

In the present study, two-dimensional numerical investigation of flow past a square cylinder beneath the free surface has been performed to identify the effects of presence of the free surface. An immersed boundary method was adopted for implementation of the cylinder cross-section in a Cartesian grid system. Also, a level-set method was used to capture the interface of two fluids. To prevent transition to three-dimensional flow, Reynolds number chosen for this simulation was 150. The cases for Froude number 0.2 and gap ratio(h/D) between 0.25 and 5.00 were examined. At the specific Reynolds number, we study the effects of gap ratio on flow characteristics around a square cylinder by computing flow fields, force coefficients and Strouhal number.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2737-2742
• /
• 2007

Strong wind flow around a building complex was numerically studied by LES. The original motivation of this work stemmed from the efforts to develop a risk assessment technique for windstorm hazards. Lagrangian-averaged scale-invariant dynamic subgrid-scale model was used for turbulence modeling, and a log-law-based wall model was employed on all the solid surfaces including the ground and the surface of buildings to replace the no-slip condition. The shape of buildings was implemented on the Cartesian grid system by an immersed boundary method. Key flow quantities for the risk assessment such as mean and RMS values of pressure on the surface of the selected buildings are presented. In addition, characteristics of the velocity field at some selected locations vital to safety of human beings is also reported.