• Proceedings of the KIEE Conference
• /
• 1999.07a
• /
• pp.366-368
• /
• 1999

In this paper, the analysis of the gas flow in the GCB is presented by using the modified FLIC method and the arc model. The modified FLIC method has two step calculation procedure. And it adopts the upwind scheme, which results in the stability.[1] The arc model used in this paper makes arc a heating source in the energy equation. The heating source is composed of ohmic heating and radiation energy transfer. And the type of the GCB in this paper is a auto-expansion type.[2]

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.90-97
• /
• 1999

An implicit turbulent flow solver is developed for 2-D unstructured hybrid meshes. Spatial discretization is accomplished by a cell-centered finite volume formulation using an upwind flux differencing. Time is advanced by an implicit backward Euler time stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one equation model, which is coupled with wall function. The numerical method is applied for flows on a flat plate, the NACA 0012 airfoil, and the Douglas 3 element airfoil. The results are compared with experimental data.

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

The performance characteristics of partial admission supersonic turbines are analyzed by using the commercial CFD program FLUENT6.0. The governing equations were discretized with Euler implicit method in time and 2nd-order upwind scheme of FVM in space. The k-$\varepsilon$ turbulence model was utilized to describe the turbulent flow field. In order to investigate the nozzle--rotor interactions and the effect of partial admission, the flows in supersonic turbine rotor cascades with a nozzle are computed. Extensive computations of partial admission supersonic turbines provide the shock structures and flow patterns in the nozzle and rotor. It is clearly shown that the nozzle flow is highly affected by the shocks or expansion waves propagated from the rotor leading edge. And the rotor flow is also affected by the shocks or wakes originated from the nozzle.

• Wind and Structures
• /
• v.1 no.2
• /
• pp.127-144
• /
• 1998

In this paper, the wake structures behind a square cylinder at the Reynolds number of 22,000 are simulated using the large eddy simulation, and the main features of the wake structure associated with unsteady vortex-shedding are investigated. The Smagorinsky model is used for parametrization of the subgrid scales. The finite element method with isoparametric linear elements is employed in the computations. Unsteady computations are performed using the explicit method with streamline upwind scheme for the advection term. The time integration incorporates a subcycling strategy. No-slip condition is enforced on the wall surface. A comparative study between two-and three-dimensional computations puts a stress on the three-dimensional effects in turbulent flow simulations. Simulated three-dimensional wake structures are compared with numerical and experimental results reported by other researchers. The results include time-averaged, phase-averaged flow fields and numerically visualized vortex-shedding pattern using streaklines. The results show that dynamics of the vortex-shedding phenomenon are numerically well reproduced using the present method of finite element implementation of large eddy simulation.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.1 no.2
• /
• pp.128-137
• /
• 1989

Numerical solutions for two-dimensional, steady, free convection are presented for a cylinder filled with saturated porous media. An annulus is bounded by inner wall with constant heat flux and two adiabatic horizontal walls with outer wall isothermally cooled. Governing equations are numerically solved for the range of Aspect Ratio 1 to 20, Radius Ratio, 1 to 20, and Rayleigh number, 50 to $10^4$ by Finite Difference method utilizing upwind scheme. Results are presented in terms of stream lines and isotherms, temperature distributions and local Nusselt numbers at the heated wall. Average Nusselt numbers are also presented for the comparisons.

• Journal of Hydrospace Technology
• /
• v.2 no.2
• /
• pp.36-43
• /
• 1996

In this work, a panel method is described, which cart solve the flow field round a surface-piercing body that experiences lift and wave resistance. As the body boundary condition, a Dirichlet type is employed, and as the free surface boundary condition the Poisson type is implemented, while in its discretization Dawson's 4-point upwind difference scheme is utilized, and as the Kutta condition a Morino-Kuo type is chosen. As to the type of singularity, source panels are distributed on the free surface, and source and dipole panels on the body surface, and dipole panels on the wake surface. For a sample run, a catamaran of the parabolic Wigley hull is chosen, for which experimental data are available, and the predictions by the numerical means and by the experiment are compared for a wide range of parameters.

• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• pp.179-186
• /
• 1999

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67), and to evaluate the performances of k-$\epsilon$ and Baldwin-Lomax turbulence models. A finite volume method is used for spatial discretization. And, the equations are solved implicitly in time with the use of approximate factorization. Upwind difference scheme is used for inviscid terms, but viscous terms are centrally differenced. The flux-difference-splitting of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. And, the results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, we conclude that k-$\epsilon$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost same.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.49 no.3
• /
• pp.172-179
• /
• 2000

In this paper, the analysis of the arc-flow in the GCB is presented by using the modified FLIC method and the arc model. The modified FLIC method adopts the upwind scheme and requires short calculation time. The arc model used in this paper treats the arc as a energy source in the energy equation. The energy source is composed of the ohmic heating and the radiation energy transfer. At each step, the movement of electrode is simulated. From the simulation, reasonable results can be obtained.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.260-265
• /
• 2008

A numerical scheme for solid propellant rocket has been studied using preconditioning method to research unsteady combustion processes for the double-base propellant with a converging-diverging nozzle. The Navier-Stokes equation is solved by dualtime stepping method with finite volume method. The turbulence model uses a shear stress transport modeling. The species equation follows up the method of Xinping WI, Mridul Kumar and Kenneth K. Kuo. A preconditioned algorithm is applied to solve incompressible regime inside the combustor and compressible flow at nozzle. Mass flux was evaluated using modified advective upwind splitting method. The simulated result the comparison a fully coupled implicit method and a semi implicit method in terms of accuracy and efficiency. This report shows the result of solid rocket propellant combustion.

• Journal of Korea Water Resources Association
• /
• v.34 no.1
• /
• pp.19-30
• /
• 2001

A laterally averaged 2-dimensional hydrodynamic model is newly developed. The coordinate system is first transformed to minimize the effects of irregularity of bottom and surface. The advection terms of the governing equations are then discretized by an upwind scheme. By employing an explicit scheme for longitudinal direction and an implicit scheme for vertical direction, the model is free from restriction of temporal step size caused by a relatively small grid ratio. To demonstrate the applicability of the model, calculated time histories of free surface displacements and distributions of velocity and salinity are compared with the field measurements of the Keum River Estuary before construction of the estuary dam. A reasonable agreement is observed between them.