• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.3
• /
• pp.31-39
• /
• 1999

In this paper, numerical experiment is attempted to analyze and compare the combustion efficiency of the burning sprays due to OFO, FOF triplet / FOOF split doublet injectors. Preconditioned Wavier-Stokes equation system with low Reynolds number $\kappa$-$\varepsilon$ model for turbulence closure, is LU-SGS time-integrated. Spray processes are modeled by DSF analysis with experimentally determined injection characteristics. n-heptane/air global reaction model approximates the combustion for simplicity, and the influence of turbulence on the chemical reaction is included using eddy dissipation model. The results showed the FOF triplet injector of highest combustion efficiency, whereas the OFO type of poet performance. It was also observed that the droplet mean diameter and the average gas temperature due to the mixing efficiency, are the representative parameters for the performance design of combustion.

• Journal of the Korean Ceramic Society
• /
• v.37 no.2
• /
• pp.111-116
• /
• 2000

700nm red emission(3F3longrightarrow3H6) in Tm3+ ion with 800 nm(3H6longrightarrow3H4) excitation via upconversion process has been reported only in host materials which have low phonon energies such as halide crystals. However, we observed 700nm and 480nm(1G4longrightarrow3H6) upconverted emission with 800nm excitation in several oxide glasses which has never reported. With spectroscopic analyses and lifetime measurements of each nergy level of Tm3+ ion doped in various oxide glasses, following mechanisms are suggested. For red upconversion, upconversion mechanism changed with Tm3+ concentration. While direct excitation up to 3F3 level via anti-Stokes excitation was dominated at low concentration, two-step excitation via 3H6longrightarrow3H4 and 3F4longrightarrow3F3 transitions was dominated at high concentration. For blue upconversion, two step excitation mechanism up to 1G4 level was suggested as follows : electrons are exciated up to 3H5 with direct excitation with pumping light up to 3H4 followed by multiphonon relaxation, and then additional reabsorption of pumping light excites electrons up to 1G4.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.127-138
• /
• 1999

In this study steady state solutions of cavity flows driven by two moving walls are studied. The north and east walls of the cavity are movable where as the remaining two walls are fixed in space. Numerical experiments for three different driving schemes for moving walls are done at two different Reynolds numbers of Re=40 and 400. The first scheme is to accelerate north and east walls simultaneously. In the second one, the north wall is started first and the east wall is accelerated later. In the third one the east wall starts first. It is usually expected that all these three cases yield the same steady state solution after sufficiently long time. However, present numerical experiments show that such a usual belief is valid only when the Reynolds number is low enough (Re=40). At higher Reynolds number (Re=400), the flow develops to three different steady states depending on the history of the boundary condition change.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.12
• /
• pp.4013-4026
• /
• 1996

A numerical simulation has been carried out for the jet impinging on a flat plate and a semi-circular concave surface. In this computation finite volume method was employed to solve the full Navier-Stokes equation based on a non-orthogonal coordinate with non staggered variable arrangement. The standard k-.epsilon. turbulent model and low Reynolds number k-.epsilon. model(Launder-Sharmar model) with Yap's correction were adapted. The accuracy of the numerical calculations were compared with various experimental data reported in the literature and showed good predictions of centerline velocity decay, wall pressure distribution and skin friction. For the jet impingement on a semi-circular concave surface, potential core length was calculated for two different nozzle(round edged nozzle and rectangular edged nozzle) to consider effects of the nozzle shape. The result showed that round edged nozzle had longer potential core length than rectangular edged nozzle for the same condition. Heat transfer rate along the concave surface with constant heat flux was calculated for various nozzle exit to surface distance(H/B) in the condition of same jet velocity. The maximum local Nusselt number at the stagnation point occurred at H/B = 8 where the centerline turbulent intensity had maximum value. The predicted Nusselt number showed good agreement with the experimental data at the stagnation point. However heat transfer predictions along the downstream were underestimated. This results suggest that the improved turbulence modeling is required.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

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

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.12
• /
• pp.1738-1746
• /
• 2002

In the FSI (Fluid-Structure Interaction) problems, two different governing equations are to be solved together. One is fur the fluid and the other for the structure. Furthermore, a kinematic constraint should be imposed along the boundary between the fluid and the structure. We use the combined formulation, which incorporates both the fluid and structure equations of motion into a single coupled variational equation so that it is not necessary to calculate the fluid force on the surface of structure explicitly when solving the equations of motion of the structure. A two-dimensional channel flow divided by a Bernoulli-Euler beam is considered and the dynamic response of the beam under the influence of channel flow is studied. The Navier-Stokes equations are solved using a P2P1 Galerkin finite element method with ALE (Arbitrary Lagrangian-Eulerian) algorithm. The internal structural damping effect is not considered in this study and numerical results are compared with a previous work fer steady case. In addition to the Reynolds number, two non-dimensional parameters, which govern this fluid-structure system, are proposed. It is found that the larger the dynamic viscosity and density of the fluid are, the larger the damping of the beam is. Also, the added mass is found to be linearly proportional to the density of the fluid.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2010.08a
• /
• pp.74-74
• /
• 2010

The aqua-plasma is the non-thermal plasma in electrical conductive electrolyte by generates the vapor film layer on the immersed metal electrode surface. This plasma can generate the hydroxyl radical by dissociate the water molecule with the plasma electron. To develop the plasma discharge device for high efficiency in the hydroxyl radical generation, proper model for estimation of plasma power is necessary. In this work, the 1-D spherical model was developed, considering temperature dependence material constants. The relation between the plasma power and hydroxyl generation was also studied by the comparison between the optical emission intensity from the hydroxyl radical using monochromator and estimated plasma power. First, the thickness of vapor layer thickness was estimated using the Navier-Stokes fluid equation in order to calculate the discharge E-field inside vapor layer. Using the E-field magnitude and power balance on the plasma generation, it was possible to estimate the plasma power. The plasma power was assumed to uniformly fill the vapor layer and the temperature of vapor layer was fixed in the boiling temperature of electrolyte, 375K. In the experiment, the aqua-plasma was discharged in the saline by applied the voltage on the bipolar electrode. The range of applied voltage was 234 to 280V-rms in the frequency of 380 kHz. Two type electrodes were produced with two ${\Phi}0.2$ tungsten. The plasma power was estimated from the V-I signal from the two high voltage probes and current probe. The estimated plasma power agreed with the profile of emission intensity when the plasma discharged between the metal electrode and vapor layer surface. However, when the plasma discharged between the metal electrodes, the increasing rate of emission intensity was lower than the increase of plasma power. It implies that the surface reaction is more sufficient rather than the volume reaction in the radical generation, due to the high density of water molecule in the liquid.

• 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 the Society of Naval Architects of Korea
• /
• v.48 no.3
• /
• pp.222-232
• /
• 2011

This article examines the scale effect of the flow characteristics, resistance and propulsion performance on a 317k VLCC. The turbulent flows around a ship in both towing and self-propulsion conditions are analyzed by solving the Reynolds-averaged Navier-Stokes equation together with the application of Reynolds stress turbulence model. The computations are carried out in both model- and full-scale. A double-body model is applied for the treatment of free surface. An asymmetric body-force propeller is used. The speed performances including resistance and propulsion factors are obtained from two kinds of methods. One is to analyze the computational results in model scale through the revised ITTC' 78 method. The other is directly to analyze the computational results in full scale. Based on the computational predictions, scale effects of the resistance and the self-propulsion factors including form factor, thrust deduction fraction, effective wake fraction and various efficiencies are investigated. Scale effects of the streamline pattern, hull pressure and local flow characteristics including x-constant sections, propeller and center plane, and transom region are also investigated. This study presents a useful tool to hull-form and propeller designers, and towing-tank experimenters to take the scale effect into consideration.