• Journal of the Korea Institute of Military Science and Technology
• /
• v.4 no.2
• /
• pp.202-214
• /
• 2001

In this paper the upwind flux difference splitting Navier-Stokes method has been applied to study quasi one-dimensional nozzle flow and axisymmetric sphere-cone($5^{\circ}$) flow for the perfect gas, the equilibrium and the nonequilibrium chemically reacting hypersonic flow. The effective gamma(${ \tilde{\gamma}}$), enthalpy to internal energy ratio was used to couple chemistry with the fluid mechanics for equilibrium chemically reacting air. The influences of the various altitude(30km, 50km) at Mach number(15.0, 20.0) were investigated. The equilibrium shock position was located farthest downstream when compared with those of perfect gas in a quasi one-dimensional nozzle. The equilibrium shock thickness over the blunt body region was much thinner than that of perfect gas shock.

• Journal of computational fluids engineering
• /
• v.2 no.2
• /
• pp.26-34
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geon unstructured triangular meshes. The flux terms are discretized based on a cell-centered formulation with the Roe's flux-difference splitting. The solution is advanced in time us backward-Euler time-stepping scheme. At each time step, the linear system of equation approximately solved wi th the Gauss-Seidel relaxation scheme. The effect of turbulence is with a standard k-ε two-equation model which is solved separately from the mean flow equation the same backward-Euler time integration scheme. The triangular meshes are generated advancing-front/layer technique. Validations are made for flows over the NACA 0012 airfoil. Douglas 3-element airfoil. Good agreements are obtained between the numerical result experiment.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.156-157
• /
• 2003

A simple method is proposed to split the flux vector of the Euler equations by introducing two artificial wave speeds. The direction of wave propagation can be adjusted by these two wave speeds. This idea greatly simplifies the upwinding, and leads to a new family of upwind schemes. Numerical flux function for multi-dimensional Euler equations is formulated for any grid system, structured or unstructured. A remarkable simplicity of the scheme is that it successfully achieves one-sided approximation for all waves without recourse to any matrix operation. Moreover, its accuracy is comparable with the exact Riemann solver. For 1-D Euler equations, the scheme actually surpasses the exact solver in avoiding expansion shocks without any additional entropy fix. The scheme can exactly resolve stationary contact discontinuities, and it is also freed of the carbuncle problem in multi­dimensional computations.

• Journal of Mechanical Science and Technology
• /
• v.18 no.12
• /
• pp.2190-2203
• /
• 2004

A combined procedure for two-dimensional Delaunay mesh generation algorithm and an adaptive remeshing technique with higher-order compressible flow solver is presented. A pseudo-code procedure is described for the adaptive remeshing technique. The flux-difference splitting scheme with a modified multidimensional dissipation for high-speed compressible flow analysis on unstructured meshes is proposed. The scheme eliminates nonphysical flow solutions such as the spurious bump of the carbuncle phenomenon observed from the bow shock of the flow over a blunt body and the oscillation in the odd-even grid perturbation in a straight duct for the Quirk's odd-even decoupling test. The proposed scheme is further extended to achieve higher-order spatial and temporal solution accuracy. The performance of the combined procedure is evaluated on unstructured triangular meshes by solving several steady-state and transient high-speed compressible flow problems.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.29-37
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geometries on unstructured triangular meshes. The flux terms are discretized based on a cell-centered finite-volume formulation with the Roe's flux-difference splitting. The solution is advanced in time using an implicit backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with the Gauss-Seidel relaxation scheme. The effect of turbulence effects is approximated with a standard $k-{\varepsilon}$ two-equation model which is solved separately from the mean flow equations using the same backward-Euler time integration scheme. The triangular meshes are generated using an advancing-front/layer technique. Validations are made for flows over the NACA0012 airfoil and the Douglas 3-element airfoil. Good agreements are obtained between the numerical results and the experiment.

• Journal of electromagnetic engineering and science
• /
• v.9 no.1
• /
• pp.12-16
• /
• 2009

This paper presents a new iterative locally one-dimensional [mite-difference time-domain(LOD-FDTD) method that has a simpler formula than the original iterative LOD-FDTD formula[l]. There are fewer arithmetic operations than in the original LOD-FDTD scheme. This leads to a reduction of CPU time compared to the original LOD-FDTD method while the new method exhibits the same numerical accuracy as the iterative ADI-FDTD scheme. The number of arithmetic operations shows that the efficiency of this method has been improved approximately 20 % over the original iterative LOD-FDTD method.

• Journal of Magnetics
• /
• v.16 no.4
• /
• pp.408-412
• /
• 2011

We investigated theoretically the magnetic field dependence of the quantum optical transition of qusi 2-Dimensional Landau splitting system, in GaN and ZnO. We apply the Quantum Transport theory (QTR) to the system in the confinement of electrons by square well confinement potential. We use the projected Liouville equation method with Equilibrium Average Projection Scheme (EAPS). Through the analysis of this work, we found the increasing properties of the optical Quantum Transition Line Shapes(QTLSs) which show the absorption power and the Quantum Transition Line Widths(QTLWs) with the magnetic-field in GaN and ZnO. We also found that QTLW, ${\gamma}(B)_{total}$ of GaN < ${\gamma}(B)_{total}$ of ZnO in the magnetic field region B < 25 Tesla.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.3 no.2
• /
• pp.23-34
• /
• 2000

Three dimensional, compressible, mass weighted averaging of Favre, Navier-Stokes system with k-$\varepsilon$ turbulence, is numerically discretized to compute three dimensional multiple jet interaction flow fields for a hybrid projectile containing three rocket motors in the ogive section. Numerical flow field computations have been made for angled nose jets and rockets at supersonic speed using multiblock structured grid. The jet conditions include very high jet to free stream pressure ratio and high temperature. It is shown that the strength of nozzle stagnation pressure affects the flow field near the side nozzle and the high stagnation pressure increases total amount of drag by a few percent. However, minor drag loss due to the pressure drag might be fully overcomed by an additional axial thrust. The results of present study can be applied for the design of future hybrid projectile.

• Journal of The Korean Astronomical Society
• /
• v.28 no.2
• /
• pp.223-243
• /
• 1995

We describe the implementation of a multi-dimensional numerical code to solve the equations for idea! magnetohydrodynamics (MHD) in cylindrical geometry. It is based on an explicit finite difference scheme on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. Multiple spatial dimensions are treated through a Strang-type operator splitting. Curvature and source terms are included in a way to insure the formal accuracy of the code to be second order. The constraint of a divergence-free magnetic field is enforced exactly by adding a correction, which involves solving a Poisson equation. The Fourier Analysis and Cyclic Reduction (FACR) method is employed to solve it. Results from a set of tests show that the code handles flows in cylindrical geometry successfully and resolves strong shocks within two to four computational cells. The advantages and limitations of the code are discussed.

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

A three-dimensional CSCM upwind flux difference splitting Navier-stokes code with two-equation turbulence models was developed to predict the transonic flows in centrifugal compressor diffuser. The k-$\epsilon$ model of Abe et al. performed well in predicting the pressure distribution in the shock wave/turbulent boundary-layer interaction. Three turbulence models predicted the similar distribution of static pressure through the diffuser and showed a good agreement with the experimental results. The secondary flows in the corner were predicted well by these turbulence models. The pressure increase before the throat of the diffuser vane is important for the overall pressure recovery. As the mass flow rate increased the blockage decreased at the throat. The pressure coefficient distribution through the diffuser depended on the throat blockage not on the rotational speed of the impeller.