• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.18-24
• /
• 2003

In the present study, solution algorithms for the connotation of unsteady flows on an unstructured mesh me presented Dual time stepping is incorporated to achieve the 2-nd order temporal accuracy while reducing the linearization and the factorization errors associated with a linear solver. Hence, any time step can be used by only considering physical phenomena. Gauss-Seidel scheme is used to solve linear system of equations. Rigid motion and suing analogy method for moving mesh are all considered and compared. Special treatments of suing analogy for high aspect ratio cells are presented. Finally, numerical results for oscillating ing are compared with experimental data.

• The KSFM Journal of Fluid Machinery
• /
• v.14 no.3
• /
• pp.10-17
• /
• 2011

Free energy based lattice Boltzmann method (LBM) has been used to simulate the contact angle and the bubble necking with large density ratio. LBM with the proper contact angle model is able to reduce the spurious currents and eliminate the singularity in the contact lines. The numerical results of the contact angles are satisfied with the Youngs law. For bubble necking flows, simulations are executed for various viscosities and contact angles. The phenomena of the bubble necking are simulated successfully and the subsequent results are presented. The present method is also applicable to the nucleate boiling flows.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.279-282
• /
• 2002

Recently, gas turbines for power generation adopt multistage DLN(Dry Low NOx) type combustion, where diffusion combustion is applied at low load and, with increase in load, the combustion mode is changed to lean premixed combustion to reduce NOx emissive concentration. However, during the mode changeover from diffusion to premixed flame, unfavorable phenomena, such as flashback, high amplitude combustion oscillations, or thermal damage of combustor parts could frequently occur. In the present study, to apply for the analysis of such unfavorable phenomena, three-dimensional CFD investigations are carried out to compare the detailed flow characteristics and temperature distribution inside the gas turbine combustor before and after combustion mode changeover. The fuel considered here is pure methane gas. A standard $k-{\varepsilon}$ turbulence model with wall function and a P-N type radiation heat transfer model, have been utilized. To analyze the complex geometric effects of combustor parts on combustion characteristics, fuel nozzles, a swirl vane f3r fuel-air mixing, and cooling air holes on the combustor liner wall, are included in this simulation.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.4B
• /
• pp.339-346
• /
• 2009

In this study, the flow characteristics of the sea region, where tidal current and shock waves are combined, are identified using a three-dimensional numerical model (Princeton Ocean Model, POM). The model is adopted and applied for simulating the flows of the sea region near the open sections during the seadike closure work of Sihwa Seadike which was closed in 1994. The simulation results show that the shock waves with high velocities propagate through the sections toward the inside and outside of the seadike during the periods of the spring and ebb tides, respectively. It is found that the phenomena of flow separation occur near the shock waves; as the shock waves extend to wider zones after passing the sections, their effects on the tidal current become weak. In addition, the longitudinal velocity profiles of the flows are revealed to be affected by the shock waves. For all the simulations, at the ebb tide, the drawdown of the water levels occurs in front of the open section, respectively, especially, hydraulic jump occurs when simulating the case of maximum difference in water level between the inside and outside of the seadike. As a result, it is thought that the flow characteristics of the sea region dominated by shock waves need to be identified employing three-dimensional analysis approach, which is expected to provide the information for ocean engineering works and facility management.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.1
• /
• pp.179-192
• /
• 1993

Three-dimensional cohesive sediment transport model, COSETM-3, is develpoed using a finite difference method. The model results are compared with the physical experimental results for the relative concentration with time at the mid-depth of the recirculating flume and are found to be in good agreement. This model is applied to Suyoung Bay in Pusan of Korea to verify the field applicability of the model and to investigate on the SS (suspended solids) diffusion phenomena at the bay. Behaviors of discharging SS from Suyoung River at normal river flow and flood river flow are predicted. The numerical results appear to be reasonable and qualitative agreement with field data. The influence of settling velocity on the concentration distribution of SS is also investigated. In case of not considering settling velocity, SS concentration at surface layer is higher than that at lower layer, but in case of considering settling velocity, SS concentration at lower layer is higher than that at surface layer. The fluctuation of SS concentration at surface layer is large due to the strong mixing, but the fluctuation of the concentration at lower layer is small due to the weak mixing. SS diffusion patterns at flood river flow are similar to those at normal river flow, while the concentration at that flow is so much higher than that at this flow. SS concentration increases with time until the peak discharge occurs, but the concentration decreases with time with decreasing river flow after the peak discharge.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.32-33
• /
• 2003

Mechanism of a periodic oscillation of shock-induced combustion over a two- dimensional wedges and axi-symmetric cones were investigated through a series of numerical simulations at off-attaching condition of oblique detonation waves(ODW). A same computational domain over 40 degree half-angle was considered for two-dimensional and axi-symmetric shock-induced combustion phenomena. For two-dimensional shock-induced combustion, a 2H2+02+17N2 mixture was considered at Mach number was 5.85with initial temperature 292 K and initial pressureof 12 KPa. The Rankine-Hugoniot relation has solution of attached waves at this condition. For axi-symmetric shock-induced combustion, a H2+2O2+2Ar mixture was considered at Mach number was 5.0 with initial temperature 288 K and initial pressure of 200 mmHg. The flow conditions were based on the conditions of similar experiments and numerical studies.[1, 3]Numerical simulation was carried out with a compressible fluid dynamics code with a detailed hydrogen-oxygen combustion mechanism.[4, 5] A series of calculations were carried out by changing the fluid dynamic time scale. The length wedge is varied as a simplest way of changing the fluid dynamic time scale. Result reveals that there is a chemical kinetic limit of the detached overdriven detonation wave, in addition to the theoretical limit predicted by Rankine-Hugoniot theory with equilibrium chemistry. At the off-attaching condition of ODW the shock and reaction waves still attach at a wedge as a periodically oscillating oblique shock-induced combustion, if the Rankine-Hugoniot limit of detachment isbut the chemical kinetic limit is not.Mechanism of the periodic oscillation is considered as interactions between shock and reaction waves coupled with chemical kinetic effects. There were various regimes of the periodicmotion depending on the fluid dynamic time scales. The difference between the two-dimensional and axi-symmetric simulations were distinct because the flow path is parallel and uniform behind the oblique shock waves, but is not behind the conical shock waves. The shock-induced combustion behind the conical shockwaves showed much more violent and irregular characteristics.From the investigation of characteristic chemical time, condition of the periodic instability is identified as follows; at the detaching condition of Rankine-Hugoniot theory, (1) flow residence time is smaller than the chemical characteristic time, behind the detached shock wave with heat addition, (2) flow residence time should be greater than the chemical characteristic time, behind an oblique shock wave without heat addition.

• Nuclear Engineering and Technology
• /
• v.40 no.3
• /
• pp.213-224
• /
• 2008

The pumping of coolant in a liquid metal fast reactor may be performed with an annular linear induction electro-magnetic (EM) pump. Linear induction pumps use a traveling magnetic field wave created by poly-phase currents, and the induced currents and their associated magnetic field generate a Lorentz force, whose effect can be the pumping of the liquid metal. The flow behaviors in the pump are very complex, including a time-varying Lorentz force and pressure pulsation, because an induction EM pump has time-varying magnetic fields and the induced convective currents that originate from the flow of the liquid metal. These phenomena lead to an instability problem in the pump arising from the changes of the generated Lorentz forces along the pump's geometry. Therefore, a magneto-hydro-dynamics (MHD) analysis is required for the design and operation of a linear induction EM pump. We have developed a time-harmonic 2-dimensional axisymmetry MHD analysis method based on the Maxwell equations. This paper describes the analysis and numerical method for obtaining solutions for some MHD parameters in an induction EM pump. Experimental test results obtained from an induction EM pump of CLIP-150 at the STC "Sintez," D.V. Efremov Institute of Electro-physical Apparatus in St. Petersburg were used to validate the method. In addition, we investigated some characteristics of a linear induction EM pump, such as the effect of the convective current and the double supply frequency (DSF) pressure pulsation. This simple model overestimated the convective eddy current generated from the sodium flow in the pump channel; however, it had a similar tendency for the measured data of the pump performance through a comparison with the experimental data. Considering its simplicity, it could be a base model for designing an EM pump and for evaluating the MHD flow in an EM pump.

• Journal of computational fluids engineering
• /
• v.9 no.2
• /
• pp.30-35
• /
• 2004

A numerical study of underexpanded jet and impingement on a wall mounted at various distances from the nozzle exit is presented. The 3-dimensional Wavier-Stokes equations and κ-ω turbulence equations are solved. The grids are constructed as overlapped grid systems to examine the distance effect. The DADI method is applied to obtain steady-state solutions. To avoid numerical instability such as the carbuncle phenomena that sometimes accompany approximate Riemann solver, the HLLE+ scheme is employed for the inviscid flux at the cell interfaces. A goal of this work is to apply a number of two-equation turbulence models based on the w equation to the impinging jet problem.

• 한국해양학회지
• /
• v.29 no.2
• /
• pp.95-106
• /
• 1994

Tidal and wind-driven currents in Chinhae Bay are investigated using a three-dimensional numerical model developed by Kim et al. (1993). The simulations indicate that the flow patterns in the bay are predominated by the bathymetry, wind and river inflow, and the effects of wind on the flow pattern in the inner bay are much stronger than those in the entrance channel. Computed tidal currents coincide with the field measurements. The horizontal and vertical velocities of tidal and residual currents are strong in the entrance channel of the bay, whereas the velocities are relatively weak in the western and northern parts of the bay. Computed velocity fields show the expected phase difference between the velocities in the surface and those in the bottom layer, and these characteristics are more remarkable during the spring tide than the neap tide. The surface currents in the bay depend strongly on the wind and river inflow, and such phenomena are more remarkable during the neap tide than the spring tide.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.3
• /
• pp.115-125
• /
• 1992

The variation of flow pattern caused by the topographical change of Kwangyang bay, is analyzed using the numerical tidal model for the depth-integrated two- dimensional long wave equation. The results of study are as follows. Due to pier construction, the area of water surface is deceased and the water inflow into the Kwangyang bay is reduced. For this result, at the outer bay of Myo island, the tidal range is slightly increased. And at the inner bay, water level is dropped generally, and especially at the time of low water tide, the phenomena of water level drop obviously appears. According to the variation pattern, flow velocities is lower than those of non-construction condition over the Kwangyang bay. But at the channel(from Kwangyang east stream) flowing into the east Kwangyang bay, for the contraction of channel profile, flow velocity is increased. The study based on the 100 year frequency design flood discharge from Sueocheon(river) and Dongcheon(river) which are flowing into the bay and Seomjin River flowing along the boundary of the bay is also performed. During the spring tide condition, the results showed the rise of water level about 1.2 m at Seomjin River Estuary and 0.3 m at inner bay is occurred.