• Journal of the Korean Society of Visualization
• /
• v.7 no.2
• /
• pp.35-41
• /
• 2010

Supersonic microjets acquire considerable research interest from a fundamental fluid dynamics perspective, in part because the combination of highly compressible flow at low-to-moderate Reynolds number is not very common, and in part due to the complex nature of the flow itself. In addition, microjets have a great variety engineering applications such as micro-propulsion, MEMS(Micro-Electro Mechanical Systems) components, microjet actuators and fine particle deposition and removal. Numerical simulations have been carried out at moderate nozzle pressure ratios and for different nozzle exit diameters to investigate and to understand in-depth of aerodynamic characteristics of supersonic microjets.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2396-2400
• /
• 2008

The motion of small heavy particles in homogeneous isotropic turbulence in the present of gravity is investigated using Direct Numerical Simulations (DNS) at moderate Reynolds number. The Lagrangian velocity and acceleration statistics of particles and of flow for a wide range of Stokes number, defined as the ratio of the particle response time to Kolmogorov time scale of turbulence, were obtained for the direction of the gravity and normal direction, respectively. It is found that particles lose their correction faster than the case without gravity. Then, a significant increase in the average settling velocity was observed for a certain range of Stokes number. Our focus is placed on gravitational effect on very small particles. Our simulations show that as the Stokes number reduces to zero, their mean settling velocity approaches the terminal velocity in still fluid.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.11
• /
• pp.895-903
• /
• 2007

CMBT(Curved Moving Boundary Treatment) is a newly developed scheme for the treatment of a no slip condition on the curved solid wall of moving obstacle in a flow field. In our research CMBT was used to perform LBM simulation of a flow over a moving circular cylinder to determine the flow feature and aerodynamics characteristic of the cylinder. To ascertain the applicability of CMBT on the complex shape of the obstacle, it was first simulated for the case of the flow over a fixed circular cylinder in a channel and the results were compared against the solution of Navier-Stokes equation with deforming mesh technique. The simulations were performed in a moderate range of reynolds number at each moving cylinder to identify the flow feature and aerodynamic characteristics of circular cylinder in a channel. The drag coefficients of the cylinder were calculated from the simulation results. We have numerically confirmed that the critical reynolds number for vortex shedding is ar Re=250 and the result is the same as the case of fixed cylinder. As the cylinder approaching to one wall, the 2nd vortex is developed by interacting with the wall boundary-layer vorticity. As the velocity ratio increase the third vortex are generated by interacting with the 2nd vortexes developed on the upper and lower wall boundary layer. The resultant $C_d$ decrease as reynolds number increasing and the Cd approached to a value when Re>1000.

• Journal of computational fluids engineering
• /
• v.11 no.3 s.34
• /
• pp.28-36
• /
• 2006

NSBT(No Slip Boundary Treatment) is a newly developed scheme for the treatment of a no slip condition on the solid wall of obstacle in a flow field. In our research, NSBT was used to perform LBM simulation of a flow over a circular cylinder to determine the flow feature and aerodynamics characteristic of the cylinder. To ascertain the applicability of NSBT on the complex shape of the obstacle, it was first simulated for the case of the flow over a circular and square cylinder in a channel and the results were compared against the solution of Navier-Stokes equation. The simulations were performed in a moderate range of Reynolds number at each cylinder position to identify the flow feature and aerodynamic characteristics of circular cylinder in a channel. The drag coefficients of the cylinder were calculated from the simulation results. We have numerically confirmed that the critical reynolds number for vortex shedding is in the range of 200$\sim$250. For the gap parameter $\gamma$ = 2 cases at Re > 240, the vortex shedding were symmetric and it resembled the Karmann vortex. As the cylinder approached to one wall, the vorticity significantly reduced in length while the vorticity on the other side elongated and the vorticity combined with the wall boundary-layer vorticity. The resultant $C_d$ by LBM concurred with the results of DNS simulation performed by previous researchers.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.12
• /
• pp.977-982
• /
• 2009

Direct numerical simulation (DNS) of a turbulent boundary layer with moderate Reynolds number was performed to scrutinize streamwise-coherence of hairpin packet motions. The Reynolds number based on the momentum thickness (${\theta}_{in}$) and free-stream velocity (U${\infty}$) was varied in the range $Re_{\theta}$=1410${\sim}$2540 which was higher than the previous numerical simulations in the turbulent boundary layer. In order to include the groups of hairpin packets existing in the outer layer, large computational domain was used (more than 50${\delta}_o$, where ${\theta}_o$ is the boundary layer thickness at the inlet in the streamwise domain). Characteristics of packet motions were investigated by using instantaneous flow fields, two-point correlation and conditional average flow fields in xy-plane. The present results showed that a train of hairpin packet motions was propagating coherently along the downstream and these structures induced the very large-scale motions in the turbulent boundary layer.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.13 no.5
• /
• pp.717-724
• /
• 2010

The initial trajectory of a spheroid configuration bobbin for precision guidance has been investigated by analyzing its aerodynamic load and six-degree-of-freedom motion. The effects of changes in the spheroidal head configuration, flow angle and lateral center-of-gravity offset are numerically studied using the commercial software "FLUENT". A wind tunnel test is also conducted to validate the numerical scheme and to examine effect of the Reynolds number on the flow around the bobbin. It is shown that the size of the separation bubble formed on the surface decreases significantly when the Reynolds number is varied between 110,000 and 140,000. At a zero flow angle, an oblate spheroidal head shows relatively moderate rotation while a prolate spheroidal head shows rapid rotation. The bobbin with a spherical head shape has little effect on the flow direction; however, the oblate bobbin is sensitive to the flow angle. The roll motion of the bobbin is greatly influenced by the lateral center-of-gravity offset and maximum dispersion is observed at half of the radius.

• Korea-Australia Rheology Journal
• /
• v.17 no.3
• /
• pp.131-140
• /
• 2005

In this work we show the results of our most recent Direct Numerical Simulations (DNS) of turbulent viscoelastic channel flow using spectral spatial approximations and a stabilizing artificial diffusion in the viscoelastic constitutive model. The Finite-Elasticity Non-Linear Elastic Dumbbell model with the Peterlin approximation (FENE-P) is used to represent the effect of polymer molecules in solution, The corresponding rheological parameters are chosen so that to get closer to the conditions corresponding to maximum drag reduction: A high extensibility parameter (60) and a moderate solvent viscosity ratio (0.8) are used with two different friction Weissenberg numbers (50 and 100). We then first find that the corresponding achieved drag reduction, in the range of friction Reynolds numbers used in this work (180-590), is insensitive to the Reynolds number (in accordance to previous work). The obtained drag reduction is at the level of $49\%\;and\;63\%$, for the friction Weissenberg numbers 50 and 100, respectively. The largest value is substantially higher than any of our previous simulations, performed at more moderate levels of viscoelasticity (i.e. higher viscosity ratio and smaller extensibility parameter values). Therefore, the maximum extensional viscosity exhibited by the modeled system and the friction Weissenberg number can still be considered as the dominant factors determining the levels of drag reduction. These can reach high values, even for of dilute polymer solution (the system modeled by the FENE-P model), provided the flow viscoelasticity is high, corresponding to a high polymer molecular weight (which translates to a high extensibility parameter) and a high friction Weissenberg number. Based on that and the changes observed in the turbulent structure and in the most prevalent statistics, as presented in this work, we can still rationalize for an increasing extensional resistance-based drag reduction mechanism as the most prevalent mechanism for drag reduction, the same one evidenced in our previous work: As the polymer elasticity increases, so does the resistance offered to extensional deformation. That, in turn, changes the structure of the most energy-containing turbulent eddies (they become wider, more well correlated, and weaker in intensity) so that they become less efficient in transferring momentum, thus leading to drag reduction. Such a continuum, rheology-based, mechanism has first been proposed in the early 70s independently by Metzner and Lamley and is to be contrasted against any molecularly based explanations.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.6
• /
• pp.1-10
• /
• 1996

In the present study, the structure and the characteristics of gaseous premixed flame impinging normal to the flat plate have been investigated experimentally. For the examination of the heat transfer and combustion characteristics, measurements of temperature, direct and schlieren photography were performed. The results of present study show that the length of inner flame becomes smaller as distance from nozzle exit to plate decrease. The width of flame becomes larger as air-fuel ratio decreases. The smaller Reynolds number at nozzle exit and the smaller distance from nozzle exit to plate lead to the higher heat transfer rate in the region of center of plate. As the air-fuel ratio decreases, the heat transfer at plate with moderate rate occurs on wide region.

• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.107-112
• /
• 2012

We explored the dynamic motions and the lateral equilibrium positions of an elastic capsule in channel flow at moderate Reynolds number varying Re, aspect ratio, size ratio, membrane stretching and bending coefficient. The transition of tank-treading/swinging to tumbling motion was observed in the simulations and the transition of dynamic motions for capsules resulted in different trend of the variation in the lateral equilibrium positions. Though other conditions were similar, the capsule with tumbling motion migrated closer to the wall than that with tank-treading motion.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.163-166
• /
• 2006

We have introduced a new version of the 3-D lid-driven cavity problem that leads to more complicated fluid parcel trajectories and thus, enhanced mixing, but at the same time weakens corner singularities. We employed an advanced form of LES to solve this problem and presented preliminary results that show very complicated streamline structures on both large and small scales, despite a relatively low Reynolds number. Finally, we demonstrated moderate speedups via parallelization. Ongoing tests are expected to resolve the questions raised regarding possible sources of the rather poor parallel performance compared with that seen in earlier studies with the same code. Because it is expected that findings may be significant for parallel performance in general, we plan to emphasize this aspect in the oral presentation the Parrel (CFD 2006 Conference.