• 한국전산유체공학회지
• /
• 제16권1호
• /
• pp.36-41
• /
• 2011

Numerical simulation of air-launched rockets from a helicopter was conducted to investigate the aerodynamic interference between air-launched rocket and helicopter. For this purpose, a three-dimensional inviscid flow solver has been developed based on unstructured meshes. An overset mesh technique was used to describe the relative motion between rocket and rocket launcher. The flow solver was coupled with six degree-of-freedom equation to predict the trajectory of free-flight rockets. For the validation, calculations were made for the impinging jet with inclined plate. The rotor downwash of helicopter was calculated and applied to simulation of air-launched rocket. It is shown that the rotor downwash has non-negligible effect on the air-launched rocket and its plume development.

• 한국전산유체공학회지
• /
• 제11권2호
• /
• pp.1-7
• /
• 2006

An unstructured overset mesh method has been developed for the simulation of unsteady viscous flow fields around multiple bodies in relative motion. For this purpose, a robust and fast search technique is proposed for both triangle and high-aspect ratio quadrilateral cell elements. The interpolation boundary is defined for data communication between grid systems and an interpolation method is suggested for viscous and inviscid cell elements. This method has been applied to calculate the flow fields around 2-D airfoils involving relative motion. Validations were made by comparing the predicted results with those of experiments or other numerical results. It was demonstrated that the present method is efficient and robust for the prediction of unsteady time-accurate flow fields involving multiple bodies in relative motion.

• 한국전산유체공학회지
• /
• 제2권1호
• /
• pp.129-137
• /
• 1997

We consider long nonlinear waves in the two-layer flow of an inviscid and incompressible fluid bounded above by a free surface and below by a rigid boundary. The flow is forced by a bump on the bottom. The derivation of the forced KdV equation fails when the density ratio h and the depth ratio ρ yields a condition 1＋hρ=(2-h)((1-h)²＋4ρh)/sup 1/2/. To overcome this difficulty we derive a forced modified KdV equation by a refined asymptotic method. Numerical solutions are given and hydraulic fall solution of a two layer fluid is expressed analytically in the case that derivation of the forced KdV(FKdV) equaition fails.

• 한국전산유체공학회지
• /
• 제20권1호
• /
• pp.32-38
• /
• 2015

A 3-dimensional compressible turbulent boundary layer solver has been developed. A time marching method is used to integrate the turbulent boundary layer equations. While the direct integration of the boundary layer equations is performed for unseparated flow regions, the inverse integration is performed for separated flow regions. The program is verified for flows that have analytical solutions or other numerical results. The solver will be merged with an Euler solver for viscous-inviscid interaction.

• 한국전산유체공학회지
• /
• 제3권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.

• 한국전산유체공학회지
• /
• 제3권2호
• /
• pp.65-72
• /
• 1998

From the Boltzmann equation with BGK approximation, a gas-kinetic BGK scheme is developed and methods for its efficient calculation, using the convergence acceleration techniques, are presented in a framework of an implicit time integration. The characteristics of the original gas-kinetic BGK scheme are improved in order for the accurate calculation of viscous and heat convection problems by considering Osher's linear subpath solutions and Prandtl number correction. Present scheme applied to various numerical tests reveals a high level of accuracy and robustness and shows advantages over flux vector splittings and Riemann solver approaches from Euler equations.

• 한국전산유체공학회지
• /
• 제14권3호
• /
• pp.33-48
• /
• 2009

In the present study, a fast grid deformation technique has been incorporated into the unsteady compressible and incompressible viscous flow solvers on unstructured hybrid meshes. An algebraic method based on the basis decomposition of normal edge vector was used for the deformation of viscous elements, and a ball-vertex spring analogy was adopted for inviscid elements among several spring analogy methods due to its robustness. The present method was validated by comparing the results obtained from the grid deformation and the rigid motion of entire grids. Fish swimming motion of an NACA0012 airfoil and flapping wing motion of a generic fighter were also simulated to demonstrate the robustness of the present grid deformation technique.

• 한국전산유체공학회지
• /
• 제19권2호
• /
• pp.86-92
• /
• 2014

In this paper, prediction of separation trajectory for Two-stage-To-Orbit space launch vehicle has been numerically simulated by using an aerodynamic database based on steady state analysis. Aerodynamic database were obtained for matrix of longitudinal and vertical positions. The steady flow simulations around the launch vehicle have been made by using a 3-D RANS flow solver based on unstructured meshes. For this purpose, a vertex-centered finite-volume method was adopted to discretize inviscid and viscous fluxes. Roe's finite difference splitting was utilized to discretize the inviscid fluxes, and the viscous fluxes were computed based on central differencing. To validate this flow solver, calculations were made for the wind-tunnel experiment model of the LGBB TSTO vehicle configuration on steady state conditions. Aerodynamic database was constructed by using flow simulations based on test matrix from the wind-tunnel experiment. ANN(Artificial Neural Network) was applied to construct interpolation function among aerodynamic variables. Separation trajectory for TSTO launch vehicle was predicted from 6-DOF equation of motion based on the interpolated function. The result of present separation trajectory calculation was compared with the trajectory using experimental database. The predicted results for the separation trajectory shows fair agreement with reference[4] solution.

• 한국전산유체공학회지
• /
• 제17권1호
• /
• pp.8-15
• /
• 2012

A code developed using the flux-difference splitting scheme on the hybrid Cartesian/immersed boundary method is applied to simulate three-dimensional internal waves. The material interface is regarded as a moving contact discontinuity and is captured on the basis of mass conservation without any additional treatment across the interface. Inviscid fluxes are estimated using the flux-difference splitting scheme for incompressible fluids of different density. The hybrid Cartesian/immersed boundary method is used to enforce the boundary condition for a moving three-dimensional body. Immersed boundary nodes are identified within an instantaneous fluid domain on the basis of edges crossing a boundary. The dependent variables are reconstructed at the immersed boundary nodes along local normal lines to provide the boundary condition for a discretized flow problem. The internal waves are simulated, which are generated by an pitching ellipsoid near an material interface. The effects of density ratio and location of the ellipsoid on internal waves are compared.

• 한국해양환경ㆍ에너지학회지
• /
• 제10권2호
• /
• pp.107-111
• /
• 2007

경계를 가지는 비점성 2유체 유동을 해석할 수 있는 수치해법을 개발하고 이를 이용하여 초기에 정지되어 있는 액체판의 끝단에서 일어나는 수축현상을 분석하였다. 경계면은 보오텍스 쉬트로 모델하였고 유동장은 와도격자법과 경계면추적법을 적용하여 계산하였다. 액체판의 끝단은 표면장력에 의해 수축되고 뭉툭한 형상을 취하게 되고, 이러한 끝단은 잘룩한 형상의 통로에 의하여 액체판에 연결되어 안쪽으로 끌려들게 된다. 이러한 현상을 비점성 조건에서 운동에너지 변화를 포함한 무차원 수의 함수로 분석하였다. 끝단의 수축속도는 로 파악되었고 표면장력파의 전파특성이 엄밀하게 조사되어야 함을 밝혔다.