• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.69-69
• /
• 2011

본 논문에서는 비정상 상태의 비압축성 유동장을 해석하기 위하여 물체맞춤격자방법이 아닌 가상경계법을 사용하였다. 가상경계법은 구조격자를 사용하여 구조물 경계면에서 Momentum Forceing을 사용하여 가상의 경계를 만들어 유동장을 해석하는 방법이다. Navier-Stoke 방정식의 수치 이산화 방법으로 Kim et al(1985)이 사용한 Fractional Step Method(FSM)을 사용하였다. 시간에 대하여 semi-implicit FSM를 사용하였고, 확산항에 대해서는 2차 정확도의 Crank-Nicolson Method를 대류항은 3차 정확도의 Runge-Kutta Method를 사용 하였다. 본 연구에서는 가상경계법을 이용한 유동장 해석이 교량 단면에 대하여 수치해석이 가능한지 검토하였다. 가상경계법은 현재 많은 연구가 유선형의 구조물에 대하여 수행되어 오고 있다. 교량 단면과 같은 각 진 구조물에 대한 검토는 아직 미비한 실정이다. 가상경계법에서 다루고 있는 구조물 경계면에서의 Momentum Forcing 방법이 유선형의 구조물에 맞추어 연구가 진행되었기 때문이다. 먼저 본 연구의 프로그램을 검증하기 위하여 원형 실린더에 대하여 가상경계법을 적용한 결과 Re 수 200에서 Strouhal Number, 양력계수, 항력계수를 이전 연구 결과와 비교하였다. Williamson(1988)과 Zhang(1995)의 연구결과와 유사한 결과를 얻을 수 있다. 그리고 교량의 단면과 같은 각진 구조물(Bluff Body)에 대하여 가상경계법 적용하였다. 본 논문의 연구에서 평가 대상으로 하고 있는 2차원 교량 단면에 대하여 유동장 해석을 하였다. 본 논문에서 정량적인 유체력과 유동장에 대한 비교 및 검토가 이루어지지 못했지만 압력장과 유선의 형태가 이론적인 값을 벗어나지 않고 있는 것으로 확인 되었다. Re 수 2700에서 전산 해석을 수행하였으며, 교량 단면 주위의 압력계수와 박리현상 그리고 후류에서의 Vortex shedding 현상이 모두 적절한 분포가 나타나는 것을 확인할 수 있었다. 따라서 가상경계법을 이용하여 각진 구조물에 대한 주위 유동장해석에 대한 가능성을 확인하였으며, 풍동실험과의 결과비교를 통하여 가상경계법을 이용하여 교량 단면 주위의 유동장 해석 결과를 정량적으로 비교할 것이다.

• Geophysics and Geophysical Exploration
• /
• v.22 no.2
• /
• pp.56-61
• /
• 2019

In this paper, we studied efficient modeling strategies when we simulate the 3D time-domain acoustic wave propagation using a cell-based finite difference method which can handle the variations of both P-wave velocity and density. The standard finite difference method assigns physical properties such as velocities of elastic waves and density to grid points; on the other hand, the cell-based finite difference method assigns physical properties to cells between grid points. The cell-based finite difference method uses average physical properties of adjacent cells to calculate the finite difference equation centered at a grid point. This feature increases the computational cost of the cell-based finite difference method compared to the standard finite different method. In this study, we used additional memory to mitigate the computational overburden and thus reduced the calculation time by more than 30 %. Furthermore, we were able to enhance the performance of the modeling on several media with limited density variations by using the cell-based and standard finite difference methods together.

• Journal of Navigation and Port Research
• /
• v.31 no.8
• /
• pp.683-687
• /
• 2007

Hydroelastic deformation of pontoon type floating structure in waves is critical in structural design. Therefore, it is necessary to develop additional technology that make to dissipate the wave energy as the submerged horizontal plate. In this study, we investigate the characteristics of hydrodynamic interaction effect by the submerged plate affecting to the radiation forces on a pontoon type floating structure using numerical analysis. We have developed the numerical method based on the composite grid system that consists of moving and fixed grid to compute the radiation forces due to the heaving motion of pontoon type floating structure and submerged plate. The numerical simulations based on the finite difference method are carried out to solve the fully nonlinear free surface involving the breaking waves and compared with the experimental data to confirm the reliability of the numerical method. Then, we discuss the interaction effects on the hydrodynamic forces that could influence on the hydroelastic response of floating structure.

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

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved by successively in the moving least square sense. Some weighing functions were tested in order to investigate the up-winding effect for the convection term. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as FVM.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.2
• /
• pp.109-123
• /
• 2008

The present paper presents the CFD result for a beam wave test case. An ONR tumblehome ship model with bilge keels is used. The beam wave test is for zero forward speed and roll and heave 2DOF with wave slope $a_k=0.156$ and wavelength ${\lambda}=1.12L_{PP}$, with $L_{PP}$ the ship length. The problems is solved numerically with an unsteady Reynolds averaged Navier-Stokes approach. The free surface flow is computed using a single-phase level-set method and the motions in each time step are integrated using a predictor-corrector iteration approach which uses dynamic overset grids moving with relative ship motion. The predicted CFD results for motions and forces are compared with experimental data, showing a reasonable agreement.

• Journal of Navigation and Port Research
• /
• v.34 no.5
• /
• pp.367-373
• /
• 2010

When the numerical analysis is carried out, it is necessary to set proper elements as a feature of analysis domains for more accurate simulations. In this study, Distinct Element Method(DEM) is applied, only considering repulsive force and tensile force except for frictional force and resisting force of particle. When the filled particles with initial Quad-tree type is relocated by DEM, a blank space existing among the particles can be minimized because the shape of particle is circular. Finally, it is the effective feature that the centroidal disposion of the particles is similar to an equilateral triangle. Triangular mesh are formed by using the Delaunay triangular technique on these relocated particles, the quality of triangular mesh is more improved by carrying out Laplace interpolations. The compared result of Aspect Ratio before and after the Laplace interpolation is shown that although the quality of triangular mesh made by DEM is good, the later triangular mesh are higher quality than the formers. In this study, although the developed technique takes a longer calculational time than the previous technique to generate triangular mesh, it is considered that the applicable possibility is very high in the generation of finite element mesh about wave analysis and various numerical simulation to need a complex or reappearance of exact topography.

• Composites Research
• /
• v.30 no.6
• /
• pp.331-337
• /
• 2017

The composite lattice structures have advantages of high specific stiffness and strength and are mainly applied to the structures of launch vehicles that carry the compressive load. However, since these structures are manufactured by filament winding technology, there are some defects and voids found in the knots. For these reasons, the stiffness and strength of the lattice structures have to be compared with finite element model for predicting design load. But, the full scale test is difficult because time and space are limited and the shape of structure is complex, and hence the simple and reliable test methods for examination of stiffness are needed. In this paper, subelements of composite lattice structures were prepared and compressive and bending test were conducted for examination of stiffness of helical and hoop rib. Test methods for subelements of composite lattice structures that has curved and twisted shape were supposed and compared with finite element analysis results.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.4
• /
• pp.18-27
• /
• 2002

Aerodynamic shape optimization of two-dimensional airfoils in inviscid compressible flows is performed using a continuous adjoint formulation on unstructured meshes. Accurate evaluation of the gradient is achieved by using a reconstruction scheme based on the Laplacian averaging. A least-square method with extended stencil is used for flow gradient calculations. Proper convergence criterion is studied on Euler and adjoint equations for efficient design. The present method has been applied to RAE2822 and NACA0012 airfoils such that wave drag can be minimized by removing the shock wave. An inverse design is also performed to recover the shock wave on the designed RAE2822 airfoil.

• Journal of computational fluids engineering
• /
• v.3 no.2
• /
• pp.17-26
• /
• 1998

The three-dimensional incompressible Navier-Stokes equations have been solved by a node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method with Jacobi matrix solver is used for the time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetragedra, prisms, pyramids, hexahedra, or mixed-element grid. Inviscid bump flow is solved to check the accuracy of high order convective flux discretisation. And viscous flows around a circular cylinder and a sphere are studied to show the efficiency and accuracy of the solver.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.7
• /
• pp.716-726
• /
• 2010

In this research, design optimization of an aircraft wing has been performed using the fully automated Multidisciplinary Design Optimization (MDO) framework, which integrates aerodynamic and structural analysis considering nonlinear structural behavior. A computational fluid dynamics (CFD) mesh is generated automatically from parametric modeling using CATIA and Gambit, followed by an automatic flow analysis using FLUENT. A computational structure mechanics (CSM) mesh is generated automatically by the parametric method of the CATIA and visual basic script of NASTRAN-FX. The structure is analyzed by ABAQUS. Interaction between CFD and CSM is performed by a fully automated system. The Response Surface Method (RSM) is applied for optimization, helping to achieve the global optimum. The optimization design result demonstrates successful application of the fully automated MDO framework.