Flow instability is investigated in a two-dimensional channel with thin baffles placed symmetrically in the vertical direction and periodically in the streamwise dircetion. At low Reynolds numbers, the flow is steady and symmetric. Above a critical Reynolds number, the steady flow undergoes a Hopf bifurcation leading to unsteady periodic flow. As Reynolds number further increases, we observe the onset of secondary instability. At high Reynolds numbers, the two-dimensional periodic flow becomes three dimmensional. To identify the onset of secondary instability, we carry out Floquet stability analysis. We obseved the transition to 3D flow at a Reynolds number of about 125. Also, we computed dominant spanwise wavenumbers near the critical Reynolds number, and visualized vortical structures associated with the most unstable spanwise wave.

The fluid-structure interaction of a towed body is simulated using a developed code, which is based on the flux-difference splitting scheme on the hybrid Cartesian/immersed boundary method. To improve the stability in the coupling between the fluid and structure domains, a scheme is used, in which the effects of structure deformation are treated implicitly. The developed code is validated for the fluid-structure interaction problem through comparisons with other results on the vortex-induced vibration of elastically mounted cylinders. To simulate behavior of a towed body, an asymmetric tension modelling for a towing cable is suggested. In the suggested model, the tension is proportional to the elongation of the cable, but the cable has no effect on the body motion whenever the distance between the endpoints of the cable is smaller than the original length of the cable. The fluid-structure interactions of a towed body are simulated on the basis of different parameters of the towing cables. It is observed that the suggested tension model predicts the snapping for a shorter towing cable, which is in accordance with the reported results.

Super-cavitating flows around under-water bodies are being studied for drag reduction and dramatic speed increase. In this paper, high speed super-cavitating flow around a two-dimensional symmetric wedge-shaped body were studied using an unsteady Reynolds-averaged Navier-Stokes equations solver based on a cell-centered finite volume method. To verify the computational method, flow over a hemispherical head-form body was simulated and validated against existing experimental data. Various computational conditions, such as different wedge angles and caviation numbers, were considered for the super-cavitating flow around the wedge-shaped body. Super-cavity begins to form in the low pressure region and propagates along the wedge body. The computed cavity lengths and velocities on the cavity boundary with varying cavitation number were validated by comparing with analytic solution.

In this study, the glass fiber drawing from a silica preform in the furnace for the optical fiber manufacturing process is numerically simulated by considering the radiative heating of cylindrically shaped preform. The one-dimensional governing equations of the mass, momentum, and energy conservation for the heated and softened preform are solved as a set of the boundary value problems along with the radiative transfer approximation between the muffle tube and the deformed preform shape, while the furnace heating is modeled by prescribing the temperature distribution of muffle tube. The temperature-dependent viscosity of silica plays an important role in formation of preform neck-down profile when the glass fiber is drawn at high speed. The calculated neck-down profile of preform and the draw tension are found to be reasonable and comparable to the actual results observed in the optical fiber industry. This paper also presents the effects of key operating parameters such as the muffle tube temperature distribution and the fiber drawing speed on the preform neck-down profile and the draw tension. Draw tension varies drastically even with the small change of furnace heating conditions such as maximum heating temperature and heating width, and the fine adjustment of furnace heating is required in order to maintain the appropriate draw tension of 100~200 g.

In this paper, our effort to apply 3-D Virtual Reality system for stereoscopic visualization of mesh data on the web is briefly described. This study is an extension of our previous and on-going research efforts to develop an automatic grid generation program specialized for wing mesh, named as eGWing. The program is developed by using JAVA programming language, and it can be used either as an application program on a local computer or as an applet in the network environment. In this research advancing layer method(ALM) augmented by elliptic smoothing method is used for the structured grid generation. And to achieve a stereoscopic viewing capability, two graphic windows are used to render its own viewing image for the left and right eye respectively. These two windows are merged into one image using 3D monitor and the viewers can see the mesh data visualization results with stereoscopic depth effects by using polarizing glasses. In this paper three dimensional mesh data visualization with stereoscopic technique combined with 3D monitor is demonstrated, and the current achievement would be a good start-up for further development of low-cost high-quality stereoscopic mesh data visualization system which can be shared by many users through the web.

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.

In this paper, various operating parameters of stream reforming process from methane in stream reformer and preconverter for MCFC is studied by numerical method. Commercial code is used to simulated the porous catalyst with user subroutine to model three dominant chemical reactions which are Stream Reforming(SR), Water-Gas Shift(WGS), and Direct Stram Reforming(DSR). The hydrogen production is tested with different wall temperature and different reactor shapes. The calculated results of the concentration of hydrogen in stream reformer are very well consistent with experimental results. This numerical study gives the design reactor wall temperature condition and size of reactor to satisfy the required fuel conversion.

Aerodynamic analysis was done for a fuselage and wing configuration of a mid-sized aircraft using unsteady 3-dimensional Navier-Stokes solver. Various turbulent models including a transitional SST were used to observe a dynamic stall as well as cruise characteristics. Also, different mesh moving methods were evaluated. Flow hysteresis which causes dynamic stall was investigated through flow field investigations.

In this study, steady and unsteady aerodynamic analyses of a huge rocket configuration have been conducted in a transonic flow region. The launch vehicle structural response are coupled with the transonic flow state transitions at the nose of the payload fairing. Before performing the coupled fluid-structure transonic aeroealstic simulations transonic aerodynamic characteristics are investigated for the pitching motions of the rocket at finite angle-of-attack. An unsteady CFD analysis method with a moving grid technique based on the Reynolds-averaged Navier-Stokes equations with the k-w SST transition turbulence model is applied to accurately predict the transonic loads of the rocket at pitching motion. It is shown that the fluctuating amplitude of the lateral aerodynamic loads imposed on the rocket due to the pitching motion can be significantly increased in the transonic flow region.

Engineering interests of submerged bodies and turbomachinery has led researchers to study various cavitation models for decades. The governing equations used for the present work are the two-phase Navier-Stokes equations with homogeneous mixture model. The solver employed on implicit dual time preconditioning algorithm in curvilinear coordinates. Three different cavitation models were applied to two axisymmetric cylinders and compared with experiments. It is concluded that the Merkle's new cavitation model has successfully accounted for cavitating flows and well captured the re-entrant jet phenomenon over the 0-caliber cylinder.

HRSG (Heat Recovery Steam Generator) is a boiler to recover heat from the exhaust gas of an engine and to generate steam for more power generation or process. For the HRSG, water-tube type boiler is commonly adopted to accommodate the working pressure or capacity requirement of the system. The water-tube type boiler has a steam drum to separate steam from the water-steam mixture supplied from the evaporator tube (riser). The drum should be sized properly to separate the steam by the gravity and auxiliary internals, such as a demister, which are installed to filter the steam. To size the steam drum and to estimate the filter efficiency of drum internals, the velocity distribution inside the drum needs to be identified. In the present study, a series of CFD has been conducted to find the velocity distributions inside steam drums for conventional HRSGs and water-tube type industrial boilers. The velocity distributions obtained from the simulation have been normalized and a correlation to predict them has been found. The correlation is applied to the steam drum design by determining a proper position of a demister to show proper separation performance.

In this paper, the vortex-induced vibration of circular cylinders is studied using the immersed boundary method on the Cartesian mesh. The Reynolds numbers considered is from 100 to 200. Using the configuration of tendemly arranged multiple circular cylinders, the vortex shedding behind of the cylinders and their flow-induced motion are investigated. The staggered MAC grid arrangement, which is the typical grid system for the incompressible flow on the Cartesian meshes, is utilized. Pressure correction method is applied for solving the divergence-free incompressible velocity field. The body motion is described by immersed boundary technique that has advantages for moving object on the fixed computational domain. It is also discussed for the computational noise in hydrodynamic forces when body motion is represented by the immersed boundary method. The Predictor/Corrector method is used for simulating the nonlinear response of the elastically mounted cylinder excited by vortex-shedding.

Laminar flow over a cube near a plane wall is numerically investigated in order to understand the effects of the cube-wall gap on the flow characteristics as well as the drag and lift coefficients. The main focus is placed on the three-dimensional vortical structures and its relation to the lift force applied on the cube. Numerical simulations are performed for the Reynolds numbers between 100 and 300, covering several different flow regimes. Without a wall nearby, the flow at Re=100 is planar symmetric with no vortical structure in the wake. However, when the wall is located close to the cube, a pair of streamwise vortices is induced behind the cube. At Re=250, the wall strengthens the existing streamwise vortices and elongates them in the streamwise direction. As a result, the lift coefficients at Re=100 and 250 increase as the cube-wall gap decreases. On the other hand, without a wall, vortex shedding takes place at Re=300 in the form of a hairpin vortex whose strength changes in time. The head of hairpin vortex or loop vortex, which is closely related to the lift force, seems to disappear due to the nearby wall. Therefore, unlike at Re=100 and 250, the lift coefficient tends to decrease more or less as the cube approaches the wall.

The Navier-Stokes equations are solved to study the flow characteristics of a micro viscous pump. The viscous micropump is consisted of a stationary disk with a spiral shaped channel and a rotating disk. A simple geometrical model for the tip clearance is proposed and validated by comparing computed flow rate with corresponding experimental data. Present numerical solutions show satisfactory agreement with the corresponding experimental data. The tip clearance effect is found to become significant as the rotational speed increases. As the pressure load increases, a reversed flow region is seen to form near the stationary disk. The height of the channel is shown to be optimized in terms of the flow rate for a given rotational speed and pressure load. The optimal height of the channel becomes small as the rotational speed decreases or the pressure load increases. The flow rate of the pump is found to be in proportion to the width of channel.