• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.8
• /
• pp.677-684
• /
• 2007

In the Present Study, a 3-D viscous flow solver, based on unstructured hybrid meshses containing tetrahedra, prisms and pyramids, has been developed. A finite-volume discretization scheme is used for solving the compressible Navier-Stokes equations. A cell-vertex median dual volume is used for spatial discretization. The one-equation Spalart-Allmaras turbulence model has been adopted to evaluate the eddy viscosity. Validation were made by computing laminar and turbulent flows around a 3-D wing for steady flows and turbulent flows around an oscillating 3-D wing in harmonic motion for unsteady flows.

• International Journal of Aeronautical and Space Sciences
• /
• v.10 no.2
• /
• pp.95-105
• /
• 2009

The flow fields around the HARTII rotor were numerically investigated using a viscous flow solver on adaptive unstructured meshes. An overset mesh and a deforming mesh technique were used to handle the blade motion including blade deflection, which was obtain from the HARTII experimental data. A solution-adaptive mesh refinement technique was also used to capture the rotor wake effectively. Comparison of the sectional normal force and pitching moment at 87% radial station between the two cases, with and without the blade deflection, showed that the blade loading is significantly affected by blade torsion. It was found that as the mesh was refined, the strength of tip vortex is better preserved, and the magnitude of high frequency blade loading, caused by blade-vortex interaction (BVI), is further magnified. It was also found that a proper time step size, which corresponds to the cell size, should be used to predict unsteady solutions accurately. In general, the numerical results in terms of the unsteady blade loading and the rotor wake show good agreement with the experimental data.

• Nuclear Engineering and Technology
• /
• v.46 no.2
• /
• pp.263-272
• /
• 2014

Scattering source calculations using conventional spherical harmonic expansion may require lots of computation time to treat full-coupled three-dimensional photon-electron transport in a highly anisotropic scattering medium where their scattering cross sections should be expanded with very high order (e.g., $P_7$ or higher) Legendre expansions. In this paper, we introduce a modified scattering kernel approach to avoid the unnecessarily repeated calculations involved with the scattering source calculation, and used it with parallel computing to effectively reduce the computation time. Its computational efficiency was tested for three-dimensional full-coupled photon-electron transport problems using our computer program which solves the multi-group discrete ordinates transport equation by using the discontinuous finite element method with unstructured tetrahedral meshes for complicated geometrical problems. The numerical tests show that we can improve speed up to 17~42 times for the elapsed time per iteration using the modified scattering kernel, not only in the single CPU calculation but also in the parallel computing with several CPUs.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.2
• /
• pp.1-9
• /
• 2003

In the present study, a 3-D Parallel DSMC method in developed on unstructured meshes for the efficient simulation of rarefied gas flows. Particle tracing between cells in achieved based on a linear shape function extended to three dimensions. For high parallel efficiency, successive domain decomposition is applied to achieve load balancing between processors by accounting for the number of particles. A particle weighting technique is also adopted to handle flows containing gases of significantly dirrerent number densities in the same flow domain. Application is made for flow past a 3-D delta wing and the result is compared with that from experiment and other calculation. Flow around a rocket payload at 100km altitude is also solved and the effect of plume back flow from the nozzle in studied.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.45 no.3
• /
• pp.259-267
• /
• 2017

Flow fields around a KARI-11-180 airfoil, SDM and transonic body are numerically simulated by using an unstructured meshes based compressible flow solver developed at KAIST. RANS equations are solved to analyse the flow fields and Roe's FDS method is adopted to evaluate convective fluxes. Turbulence effect of the flow fields is modeled by a SA model, SST model and ${\gamma}-{\widetilde{Re}}_{{\theta}t}$ model. It is found that smaller drag coefficients are predicted for the KARI-11-180 airfoil when a transition phenomenon is considered and small deviations exist between CFD and EFD results. For the SDM, flow separation is observed at a leading edge and calculated aerodynamic properties show similar tendencies to experimental results. A shock wave on main wings of the transonic body is successfully captured by the present flow solver at a Mach number 0.9. Estimated pressure profiles by means of the present CFD method also agree well with those of wind tunnel results.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.15-23
• /
• 2016

In the present study, two phase flows around a projectile vertically launched from an underwater platform have been numerically investigated by using a three dimensional multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. The relative motion between the platform and projectile was described by six degrees of freedom equations of motion with Euler angles and a chimera technique. The propulsive power of the projectile was modeled as the fluid force acting on the lower surface of the body by the compressed air emitted from the underwater platform. Various flow conditions were considered to analyze the fluid-dynamics motion parameters of the projectile. The water level of platform and the current speed around the projectile were the main parametric variables. The numerical calculations were conducted up to 0.75sec in physical time scale. The dynamics tendency of the projectile was almost identical with respect to the water level variation due to the constant buoyancy term. The moving speed of the projectile along the vertical axis inside the platform decreased when the current speed increased. This is because the inflow from outside of the platform impeded development of the compressed air emitted from the floor surface of the launch platform. As a result, the fluid force acting on the lower surface of the projectile decreased, and injection time of the projectile from the platform was delayed.

• Journal of Mechanical Science and Technology
• /
• v.15 no.5
• /
• pp.592-601
• /
• 2001

An automatic mesh generation scheme with unstructured quadrilateral elements on trimmed NURBS surfaces has been developed. In this paper NURBS surface geometries in the IGES format have been employed to represent geometric models. For unstructured mesh generation with quadrilateral elements, a domain decomposition algorithm employing loop operators has been modified. As for the surface meshing, an indirect 2D approach is proposed in which both quasi-expanded planes and projection planes are employed. Sampled meshes for complex models are presented to demonstrate the robustness of the algorithm.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.73-77
• /
• 2004

In the present study, an efficient and robust implicit operator for the LU-SGS method is proposed. Numerical experiments for supersonic flow are performed to demonstrate the performance of the proposed method.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.2
• /
• pp.196-202
• /
• 2012

In this study, numerical analysis of viscous flows is carried out based on the unstructured grid. There exist some difficulties in expressing and computing numerical derivatives on the unstructured grid due to lack of the structured characteristics. The general computer algorithms are developed to perform numerical derivatives easily and extended to be applicable to various geometries composed of hybrid meshes. And the optimal method of strongly implicit procedure is newly contrived to accelerate the rate of convergence in solving the pressure Poisson equation. To verify numerical schemes, the driven cavity problems of 2 and 3 dimension are simulated. The numerical results are compared with others and our numerical schemes are shown to be valid.

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

In the present study, a CFD program is developed for the Fluid-Structure Interaction(FSI) analysis. The non-staggered, non-orthogonal, and unstructured grid system was also used to handle the complicated geometries in the program. In order to validate the capabilities of the developed CFD program, various models are investigated by using unstructured and nonorthogonal meshes. The predicted results are a good agreement with analytic solution, experimental data and commercial software. And also PISO algorithm is applied for transient flow analysis. The cyclic boundary condition and baffle cell are developed in order to improve the effectiveness of the calculation for complex geometry.