• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.1
• /
• pp.13-26
• /
• 2006

When the stagnation temperature of a perfect gas increases, the specific heat for constant pressure and ratio of the specefic heats do not remain constant any more and start to vary with this temperature. The gas remains perfect: its state equation remains always valid, with exception that it will be named by calorically imperfect gas. The aim of this research is to develop the relations of the necessary thermodynamics and geometrical ratios. and to study the supersonic flow at high temperature. lower than the threshold of dissociation. The results are found by the resolution of nonlinear algebraic equations and integration of complex analytical functions where the exact calculation is impossible. The dichotomy method is used to solve the nonlinear equation. and the Simpson algorithm for the numerical integration of the found integrals. A condensation of the nodes is used. Since. the functions to be integrated have a high gradient at the extremity of the interval of integration. The comparison is made with the calorifcally perfect gas to determine the error made by this last. The application is made for the air in a supersonic nozzle.

• Journal of Ocean Engineering and Technology
• /
• v.27 no.2
• /
• pp.39-46
• /
• 2013

In this study, the impact loads on tank walls by sloshing phenomena and on a tall structure in a three-dimensional rectangular tank were predicted using multiphase flow simulations. The solver was based on the CIP/CCUP (Constraint interpolation CIP/CIP combined unified procedure) method, and the THINC-WLIC (Tangent hyperbola for interface capturing-weighted line interface calculation) scheme was used to capture the air-water interface. For the convection terms of the Navier-Stokes equations, the USCIP (Unsplit semi-lagrangian CIP) method was adopted. The results of simulations were compared with those of experiments. Overall, the comparisons were reasonably good.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.186-191
• /
• 2004

The reactive flowfield of the transverse injecting combustor has been studied using Euler-Lagrange method in order to develop an efficient solution procedure for the understanding of liquid spray combustion in the transverse injecting combustor which has been widely used in ramjets and turbojet afterburners. The unsteady two-dimensional gas-phase equations have been represented in Eulerian coordinates and the liquid-phase equations have been formulated in Lagrangian coordinates. The gas-phase equations based on the conservation of mass, momentum, and energy have been supplemented by combustion. The vaporization model takes into account the transient effects associated with the droplet heating and the liquid-phase internal circulation. The droplet trajectories have been determined by the integration of the Lagrangian equation in the flow field obtained from the separate calculation without considering the iterative effect between liquid and gas phases. The reported droplet trajectories had been found to deviate from the initial conical path toward the flow direction in the very end of its lifetime when the droplet size had become small due to evaporation. The integration scheme has been based on the TEACH algorithm for gas-phase equation, the second order Runge-Kutta method for liquid-phase equations and the linear interpolation between the two coordinate systems. The calculation results has shown that the characteristics of the droplet penetration and recirculation have been strongly influenced by the interaction between gas and liquid phases in such a way that most of the vaporization process has been confined to the wake region of the injector, thereby improving the flame stabilization properties of the flowfield.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.6
• /
• pp.725-733
• /
• 1998

The pressure fluctuation on the surface of a submerged body has been recognized as a dominant noise source. There have been many studies concerning the flow induced noise on a flat plate. However, the noise over an axisymmetric body has not been well reported. This paper addresses the way in which we have investigated the mechanism of noise generation due to an axisymmetric body. The associated experiments and signal processing methods are introduced. A 3-dimensional axisymmetric body whose length and diameter were 2 m and 10.4 cm, was prepared as a test specimen. The wall pressure on the surface of the body was measured in a large scale low noise wind tunnel at KIMM(Korea Institute of Machinery and Metals). To measure the wall pressure, we used two microphone arrays which were tangential and normal to the flow. Based on the measured signal, frequency-wavenumber spectrum which explains the structure of turbulence noise, was estimated. Tangential to the flow, there exists convective ridge at a relatively higher wavenumber region; this can cause spatial aliasing. To circumvent this problem, the cross spectrum was interpolated. The interpolation has been performed by unwrapping the phase and smoothing the cross spectrum. The phase unwrapping was done based on the Corcos model; the phase of cross spectrum decreases linearly with the distance between microphones. Aforementioned signal processings are possible by employing the experimental results that the estimated wavenumber spectrum quite resembles the Corcos model. We try to modify the Corcos model which is applicable to the flat plate, by altering the magnitude of cross spectrum to fit the experimental data more accurately. We proposed that this wavenumber spectrum model is suitable for the 3-dimensional axisymmetric body. Normal to the flow, there exists a little correlation between signals of different microphones. The circumferential wavenumber spectrum contains uniform power along the wavenumbers.

• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2218-2229
• /
• 2006

A conservative finite-volume numerical method for unstructured grids with the cell-centered method has been developed for computing flow and heat transfer by combining the attractive features of the existing pressure-based procedures with the advances made in unstructured grid techniques. This method uses an integral form of governing equations for arbitrary convex polyhedra. Care is taken in the discretization and solution procedure to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. For both convective and diffusive fluxes the forms superior to both accuracy and stability are particularly adopted and formulated through a systematic study on the existing approximation ones. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are computed by using a linear reconstruction based on the divergence theorem. Momentum interpolation is used to prevent the pressure checkerboarding and a segregated solution strategy is adopted to minimize the storage requirements with the pressure-velocity coupling by the SIMPLE algorithm. An algebraic solver using iterative preconditioned conjugate gradient method is used for the solution of linearized equations. The flow analysis code (PowerCFD) developed by the present method is evaluated for its application to several 2-D structured-mesh benchmark problems using a variety of unstructured quadrilateral and triangular meshes. The present flow analysis code by using unstructured grids with the cell-centered method clearly demonstrate the same accuracy and robustness as that for a typical structured mesh.

• Journal of computational fluids engineering
• /
• v.19 no.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.

• Spatial Information Research
• /
• v.20 no.5
• /
• pp.1-13
• /
• 2012

This study intends to evaluate the seasonal flow direction of carbon dioxide in Northeast Asia by using GOSAT, the first Greenhouse Observing SATellite, in an attempt to overcome costly, laborious and time consuming ground observation which has been frequently pointed out in existing studies. For this purpose, missing values were supplemented by applying the Kriging interpolation and the overall flow direction of carbon dioxide was determined through anisotoropy semi-variogram. As a result, it was found that the overall spatial distribution of carbon dioxide in Northeast Asia varies depending on the latitude, and that carbon dioxide mainly flows southeast or east in spring, autumn and winter, but northeast or north in summer. Similar to the flow of monsoons in Northeast Asia, these results show that carbon dioxide flows mainly from the west to the east, which proves that carbon dioxide discharged from China is influencing even the Korean Peninsula and Japan. However, as the flow of carbon dioxide varies depending on a variety of factors such as artificial sources, plant respiration, and the absorption and discharge of the ocean, follow-up studies are requested to evaluate such variables and the correlations.

• Journal of the Society of Naval Architects of Korea
• /
• v.43 no.5 s.149
• /
• pp.538-549
• /
• 2006

A code is developed to simulate a viscous flow field around a deformable body using the hybrid Cartesian/immersed boundary method. In this method, the immersed boundary(IB) nodes are defined near the body boundary then velocities at the IB nodes are reconstructed based on the interpolation along the normal direction to the body surface. A new method is suggested to define the IB nodes so that a closed fluid domain is guaranteed by a set of IB nodes and the method is applicable to a zero-thickness body such as a sail. To validate the developed code, the vorticity fields are compared with other recent calculations where a cylinder orbits and moves into its own wake. It is shown the code can handle a sharp trailing edge at Reynolds number of $10^5$ under moderate requirements on girds. Finally the developed code is applied to simulate the vortex shedding behind a deforming foil with flapping tail like a fish. It is shown that the acceleration of fluids near the flapping tail contributes to the generation of the thrust for propulsion.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.3
• /
• pp.29-36
• /
• 2002

A numerical analysis based on two-dimensional and three-dimensional incompressible Wavier-Stokes equations has been carried out for double-circular-arc compressor cascades and the results are compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE algorithm adopt pressure weighted interpolation method for non-staggered grid and hybrid scheme for the convective terms. Turbulence modeling is very important for prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. Considering computation times, $\kappa$-$\varepsilon$ turbulence model with wall function is used.

• Journal of Biomedical Engineering Research
• /
• v.20 no.1
• /
• pp.29-36
• /
• 1999

The objective of this study is to analyse the flow field in the branching duct by visualizing the flow phenomena using the PIV system. A bifurcation model is fabricated with transparent acrylic resin to visualize the whole flow field with the PIV system. Water was used as the working fluid and the conifer powder as the tracer particles. The single-frame and two-frame methods of the PIV system and 2-frame of the grey level correlation method are applied to obtain the velocity vectors from the images captured in the flow filed. The velocity distributions in a lid-driven cavity flow are compared with the so-called standard experimental data, which was obtained from by 4-frame method in order to validate experimental results of the PIV measurements. The flow patterns of a Newtonian fluid in a branching duct were successfully visualized by using the PIV system and the sub-pixel and the area interpolation method were used to obtain the final velocity vectors. The velocity vectors obtained from the PIV system are in good agreement with the numerical results of the 3-dimensional branch flow. The results of numerical analyses and the PIV experiments for the three-dimensional flows in the branch ing duct show the recirculation zone distal to the branching point and the sizes of the recirculation length and height of the tow different methods are in good agreement.