• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.145-150
• /
• 2005

Generally, Compressible Navier-Stokes codes are used to solve high mach number flows. But, Most of high mach number flows embrace low mach number flows. This phenomenon results in low convergence rate and non-physical solution in CFD analysis. So Many researchers developed preconditioning technique to solve these problems. This Study presents how to modify previous compressible N-S computer code with little changes of structure into preconditioned compressible N-S code applying Roe's Approximate Riemann Solver. And this study show developed preconditioning code is very well operated at all mach number flows.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.48-54
• /
• 1998

Three-dimensional incompressible Navier-Stokes equations have been solved by the 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 is used for time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetrahedra, prisms, pyramids, hexahedra, or mixed-element grid. The numerical efficiency and accuracy of the present method is critically evaluated for several example problems.

• 한국HCI학회:학술대회논문집
• /
• 2009.02a
• /
• pp.463-468
• /
• 2009

This paper proposes a method for extending simulating fluid on mobile device, which was only possible on desktop PC. Fluid simulation is done by solving Navier-Stokes equation numerically, and previous research were mainly focused on numerical stability [1], and realism [2]. However, such methods assume rich computational resources, which is not available on mobile devices. On the other hand, rigid-body solver is the mostly used physically-based technique [3], and only simple height field-based method is released for fluid simulation [4]. To overcome these problems, we proposes a modified incompressible fluid dynamics solver for the mobile device, and also we propose a technique for visualizing fluids on the mobile device.

• The Journal of the Acoustical Society of Korea
• /
• v.40 no.4
• /
• pp.261-269
• /
• 2021

Without any validation of the incompressible assumption, most of previous studies on cavitation flow and its noise have utilized numerical methods based on the incompressible Reynolds Average Navier-Stokes (RANS) equations because of advantage of its efficiency. In this study, to investigate the effects of the flow compressibility on the Tip Vortex Cavitation (TVC) flow and noise, both the incompressible and compressible simulations are performed to simulate the TVC flow, and the Ffowcs Williams and Hawkings (FW-H) integral equation is utilized to predict the TVC noise. The DARPA Suboff submarine body with an underwater propeller of a skew angle of 17 degree is targeted to account for the effects of upstream disturbance. The computation domain is set to be same as the test-section of the large cavitation tunnel in Korea Research Institute of Ships and Ocean Engineering to compare the prediction results with the measured ones. To predict the TVC accurately, the Delayed Detached Eddy Simulation (DDES) technique is used in combination with the adaptive grid techniques. The acoustic spectrum obtained using the compressible flow solver shows closer agreement with the measured one.

• 한국전산유체공학회:학술대회논문집
• /
• 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 Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.8
• /
• pp.10-20
• /
• 2002

The present paper investigates methods to control dynamic stall using an optimal approach. An unsteady aerodynamic sensitivity analysis code is developed by a direct differentiation method from a two-dimensional unsteady compressible Navier-Stokes solver including a two-equation turbulence model. Dynamic stall control is conducted by minimizing an objective function defined at an instant instead of integrating for a period of time. Unsteady sensitivity derivatives of the objective function are calculated by the sensitivity code, and optimization is carried out using a linear line search method at every physical time step. Numerous examples of dynamic stall control using control parameters such as nose radius, maximum thickness of airfoil, or suction show satisfactory results.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.6
• /
• pp.8-16
• /
• 2003

Aerodynamic characteristics over Micro Aerial Vehicle(MAV) in low Reynolds number regime are numerically studied using 3-D unsteady, incompressible Navier-Stokes flow solver with single partitioning method for multi-block grid. For more efficient computation of unsteady flows, this flow solver is parallel-implemented with MPl(Message Passing Interface) programming method. Firstly, MAV wing with not complex geometry is considered and then, we analyze aerodynamic characteristics over full MAV configuration varying the angle of attack. Present computational results show a better agreement with the experimental data by MACDL(Micro Aerodynamic Control and Design Lab.), Seoul National University. We can also find the conceptually designed MAV by MACDL has the static stability.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.17-29
• /
• 2004

Internal waves are computed using a ghost fluid method on an unstructured grid. Discontinuities in density and dynamic pressure are captured in one cell without smearing or oscillations along a multimaterial interface. A time-accurate incompressible Navier-Stokes/Euler solver is developed based on a three-point backward difference formula for the physical time marching. Artificial compressibility is introduced with respect to pseudotime and an implicit method is used for the pseudotime iteration. To track evolution of an interface, a level set function is coupled with the governing equations. Roe's flux difference splitting method is used to calculate numerical fluxes of the coupled equations. To get higher order accuracy, dependent variables are reconstructed based on gradients which are calculated using Gauss theorem. For each edge crossing an interface, dynamic pressure is assigned for a ghost node to enforce the continuity of total pressure along the interface. Solitary internal waves are computed and the results are compared with other computational and experimental results.

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

The flowfield around transonic wing-body configuration was simulated using in-house CFD code and compared with the experimental data to understand the influence of several features of CFD(Computational Fluid Dynamics) ; grid dependency, turbulence models, spatial discretization, and viscosity. The wing-body configuration consists of a simple planform RAE Wing 'A' with an RAE 101 airfoil section and an axisymmetric body. The in-house CFD code is a compressible Euler/Navier-Stokes solver based on unstructured grid. For the turbulence model, the $k-{\omega}$ model, the Spalart-Allmaras model, and the $k-{\omega}$ SST model were applied. For the spatial discretization method, the central differencing scheme with Jameson's artificial viscosity and Roe's upwind differencing scheme were applied. The results calculated were generally in good agreement with experimental data. However, it was shown that the pressure distribution and shock-wave position were slightly affected by the turbulence models and the spatial discretization methods. It was known that the turbulent viscous effect should be considered in order to predict the accurate shock wave position.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.12
• /
• pp.1-8
• /
• 2005

A parametric study has been accomplished to figure out the effects of the elliptic cylinder thickness, angle of attack, and Reynolds number on the lift and drag forces exerted on the elliptic cylinder. A two-dimensional incompressible Navier-Stokes flow solver is developed using SIMPLER method to analyze the unsteady viscous flow over elliptic cylinder. Thickness-to-chord ratios of 0.2, 0.4, and 0.6 elliptic cylinders are simulated at different Reynolds numbers of 400 and 600, and angles of attack of $10^{\circ}$, $20^{\circ}$, and $30^{\circ}$. Through this study, it is observed that the elliptic cylinder thickness, angle of attack, and Reynolds number affect significantly not only the time-mean values and the amplitudes of the drag and lift forces but also the frequencies of the force oscillations.