• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.4
• /
• pp.601-613
• /
• 2017

An adjoint-based error estimation and mesh adaptation study is conducted for two-dimensional viscous flows on unstructured hybrid meshes. The error in an integral output functional of interest is estimated by a dot product of the residual vector and adjoint variable vector. Regions for the mesh to be adapted are selected based on the amount of local error at each nodal point. Triangular cells in the adaptive regions are refined by regular refinement, and quadrangular cells near viscous walls are bisected accordingly. The present procedure is applied to single-element airfoils such as the RAE2822 at a transonic regime and a diamond-shaped airfoil at a supersonic regime. Then the 30P30N multi-element airfoil at a low subsonic regime with a high incidence angle (${\alpha}=21deg.$) is analyzed. The same level of prediction accuracy for lift and drag is achieved with much less mesh points than the uniform mesh refinement approach. The detailed procedure of the adjoint-based mesh refinement for the multi-element airfoil case show that the basic flow features around the airfoil should be resolved so that the adjoint method can accurately estimate an output error.

• Ocean Systems Engineering
• /
• v.1 no.4
• /
• pp.317-336
• /
• 2011

A time domain finite element based method is employed to analyze wave radiation by multiple cylinders. The nonlinear free surface and body surface boundary conditions are satisfied based on the perturbation method up to the second order. The first- and second-order velocity potential problems at each time step are solved through a finite element method (FEM). The matrix equation of the FEM is solved through an iteration and the initial solution is obtained from the result at the previous time step. The three-dimensional (3D) mesh required is generated based on a two-dimensional (2D) hybrid mesh on a horizontal plane and its extension in the vertical direction. The hybrid mesh is generated by combining an unstructured grid away from cylinders and two structured grids near the cylinder and the artificial boundary, respectively. The fluid velocity on the free surface and the cylinder surface are calculated by using a differential method. Results for various configurations including two-cylinder and four-cylinder cases are provided to show the mutual influence due to cylinders on the first and second waves and forces.

• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.3
• /
• pp.403-413
• /
• 2017

Dual-point aerodynamic design optimization is conducted for DLR-F6 wing-body-nacelle-pylon configuration adopting an efficient surface mesh movement method for complex junction geometries. A three-dimensional unstructured Euler solver and its discrete adjoint code are utilized for flow and sensitivity analysis, respectively. Considered design conditions are a low-lift condition and a cruise condition in a transonic regime. Design objective is to minimize drag and reduce shock strength at both flow conditions. Shape deformation is made by variation of the section shapes of inboard wing and pylon, nacelle vertical location and nacelle pitch angle. Hicks-Henne shape functions are employed for deformation of the section shapes of wing and pylon. By the design optimization, drag coefficients were remarkably reduced at both design conditions retaining specified lift coefficient and satisfying other constraints. Two-point design results show mixed features of the one-point design results at low-lift condition and cruise conditions.

• Journal of the Korea Computer Graphics Society
• /
• v.26 no.1
• /
• pp.17-25
• /
• 2020

We present a novel multi-core CPU based parallel algorithm for the cell-connectivity information extraction algorithm, which is one of the preprocessing steps for volume rendering of unstructured grid data. We first check the synchronization issues when parallelizing the prior serial algorithm naively. Then, we propose a 3-step parallel algorithm that achieves high parallelization efficiency by removing synchronization in each step. Also, our 3-step algorithm improves the cache utilization efficiency by increasing the spatial locality for the duplicated triangle test process, which is the core operation of building cell-connectivity information. We further improve the efficiency of our parallel algorithm by employing a memory pool for each thread. To check the benefit of our approach, we implemented our method on a system consisting of two octa-core CPUs and measured the performance. As a result, our method shows continuous performance improvement as we add threads. Also, it achieves up to 82.9 times higher performance compared with the prior serial algorithm when we use thirty-two threads (sixteen physical cores). These results demonstrate the high parallelization efficiency and high cache utilization efficiency of our method. Also, it validates the suitability of our algorithm for large-scale unstructured data.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.2
• /
• pp.41-51
• /
• 2013

The tidal creek dependent mass transport characteristic in Gyeong-Gi Bay (west coast of Korea) was studied using field measured data and numerical model. Gyeong-Gi Bay consists of 3 main tidal channels and contains a well-developed vast tidal flat. This region is famous for its large tidal difference and strong current. We aim to study the effect of tidal creek in the tidal flat on the mass exchange between the estuary and the ocean. For numerical application, the application of unstructured grid feature is essential, since the tidal creek has complicated shape and form. For this purpose, the FVCOM is applied to the study area and simulation is performed for 2 different cases. In case A, geographic characteristics of the tidal creek is ignored in the numerical grid and in case B, the tidal creek are constructed using unstructured grid. And these 2 cases are compared with the field measured cross-channel mass transport data. The cross-channel mass transport at the Yeomha waterway mouth and Incheon harbor was measured in June, 9~10 (Spring tide) and 17~18 (Neap tide), 2009. CTD casting and ADCP cross-channel transect was conducted 13 times in one tidal cycle. The observation data analysis results showed that mass transport has characteristic of the ebb dominance Line 1 (Yeomha waterway mouth), on the other hand, a flood dominant characteristic is shown in Line 2 (Incheon harbor front). By comparing the numerical model (case A & B) with observation data, we found that the case B results show much better agreement with measurement data than case A. It is showed that the geographic feature of tidal creek should be considered in grid design of numerical model in order to understand the mass transport characteristics over large tidal flat area.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.11
• /
• pp.759-765
• /
• 2017

In this study, a three-dimensional least-square, level-set-based two-phase flow code was developed for the simulation of three-dimensional sloshing problems using finite element discretization. The code was validated by solving some benchmark problems. The proposed method was found to provide improved results against other existing methods, by using a coarser mesh. The results of the numerical experiments conducted during the course of this study showed that the proposed method was both robust and accurate for the simulation of three-dimensional sloshing problems. Using a substantially coarse grid, historical results of the dynamic pressure at a selected position corresponded with existing experimental data. The pressure history with a finer grid was similar to that of a coarse grid; however, a fine grid provided higher peak pressures. The present method could be extended to the analysis of a sloshing problem in a complex geometrical configuration using unstructured meshes owing to the features of FEM.

• The Bulletin of The Korean Astronomical Society
• /
• v.40 no.1
• /
• pp.65.2-65.2
• /
• 2015

We developed a three-dimensional magnetohydrodynamic (MHD) code to reproduce the structure of a solar wind, the properties of a coronal mass ejection (CME) and the interaction between them. This MHD code is based on the finite volume method incorporating total variation diminishing (TVD) scheme with an unstructured grid system. In particular, this grid system can avoid the singularity at the north and south poles and relax tight CFL conditions around the poles, both of which would arise in a spherical coordinate system (Tanaka 1994). In this model, we first apply an MHD tomographic method (Hayashi et al. 2003) to interplanetary scintillation (IPS) observational data and derive a solar wind from the physical values obtained at 50 solar radii away from the Sun. By comparing the properties of this solar wind to observational data obtained near the Earth orbit, we confirmed that our model captures the velocity, temperature and density profiles of a solar wind near the Earth orbit. We then insert a spheromak-type CME (Kataoka et al. 2009) into the solar wind to reproduce an actual CME event. This has been done by introducing a time-dependent boundary condition to the inner boundary of our simulation domain. On the basis of a comparison between a simulated CME and observations near the Earth, we discuss the physics involved in an ICME interacting with a solar wind.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.71-76
• /
• 2009

The research to predict the resistance performance of the ship using the CFD analysis is increasing. For the CFD numerical analysis the computational mesh, which is proper to computational model, has to be made before the analysis is begun. In the parametric study, even though the deformation of each case is not very sharp, the whole computational mesh should be regenerated according to the conventional way. Hence, lots of effort is needed to repeated mesh generation work. To solve these problems, the automatic mesh generation method using the macro function of commercial CAD program and mesh generation program is introduced in this study. First, in the CAD program, by using the macro function and putting the deformation rate of bow and stern in lengthwise, the repeated modeling work is performed automatically. Next, the generated geometries are read by the mesh generation program and the proper mesh for the geometry is created automatically also using the macro function. The hybrid mesh which has unstructured grid near the bow and stern and structured grid in the remaining part of domain is used. The verification of the developed method is done by applying the method to predict the resistance performance of the podded propulsion cruise ship of the Daewoo Shipbuilding & Marine Engineering (DSME) in the cases of different length of bow and stern and pod set in different position. The author believes that the introduced method can help to make the database to optimize the resistance performance of the ship in various cases can be constructed without difficulty.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.11a
• /
• pp.848-854
• /
• 2006

In this study, a computational analysis system has been developed in order to investigate flow-induced vibration(FIV) phenomenon for general stator-rotor cascade configurations. Relative movement of the rotor with respect to stator is reflected by modeling independent two computational domains. Fluid domains are modeled using the unstructured grid system with dynamic moving and local deforming methods. Unsteady, Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation SST $k-\omega$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used flow computing the coupled governing equations of the fluid-structure interaction problem. Detailed FIV responses for different flow conditions are presented with respect to time and vibration characteristics are also physically investigated in the time domain.

• The KSFM Journal of Fluid Machinery
• /
• v.19 no.5
• /
• pp.42-49
• /
• 2016

Aerodynamic analysis of a low-voltage electric motor has been performed with various inlet vent shapes. Effects of inlet vent shape on aerodynamic performance of a motor cooling fan have been investigated numerically using three-dimensional Reynolds-averaged Navier-Stokes equations. The k-${\varepsilon}$ turbulence model was used for the analysis of turbulence. The finite volume method and unstructured tetrahedral grids were used in the numerical analysis. Optimal grid system in the computational domain was selected through a grid-dependency test. From the results of the flow analysis, considerable energy loss by flow separation was observed in the flow passage. It was found that mass flow rate through the cooling fan in the low-voltage motor can be increased by modifying the inlet vent shape. And, some inlet vent shapes were suggested to improve the aerodynamic performance of the motor cooling fan.