The numerical simulation techniques developed for the flows of non-Newtonian fluids including the viscoelastic fluids and the short fiber suspensions are introduced. And the causes of the numerical breakdown and the recent efforts to resolve them are presented in particular.

극동아시아 지역국가에서 곡수지연은 역사적으로 크게 번성했다. 이들 고대국가의 지배층은 유상곡수 둘레에 앉아서 시를 지었으며 술을 채운 술잔이 앉은 자리에 도달할 때까지 시를 짓지 못하면 벌주를 내리는 유희(遊戱)를 즐겼다. 본 논문에서는 한국, 중국, 일본에 존재하는 곡수의 역사를 간결히 고찰한다. 또한 곡수지연을 통하여 고대인물들이 즐겼을 유상곡수를 지금은 완전히 기능이 정지된 포석정에 대하여 현대적인 전산기법의 도움을 받아 약간의 근사를 써서 시뮬레이션하였다. 또한, 술잔의 흥미로운 유로(流路)를 예측하여 보았다. 본 연구를 통하여 우리의 선조들이 즐겼던 우아했던 삶의 방식하나를 자세히 들여다 보는 것이 가능해졌다.

A previously developed pressure based calculation procedure for incompressible flows was modified and applied to transonic and supersonic flows. It uses pressure as a primary variable in preference to density and body fitted coordinate and non-staggered grid system. The discretized momentum equations were rearranged as a system of equations with respect to covariant velocity components. Three different discretization schemes, QUICK hybrid and first order upwind, were used to approximate the convective terms and compared. Present approach was tested far two transonic flow and one supersonic flow problems. Comparison with previous results show that present approach can be used as a solver for compressible flows.

An incompressible flow stover based on the unstructured finite volume method has been developed. The flow domain is discretized by triangles in 2D or tetrahedra in 3D. The convective and viscous fluxes are obtained using edge connectivities of the unstructured meshes. The pressure-velocity coupling is handled by the artificial compressibility algorithm due to its computational efficiency associated with the hyperbolic nature of the resulting equations. Laminar test flow problems are computed and presented with a comparison against other numerical solutions or experimental results.

This paper discribes the viscous interaction of Hypersonic Wedge Flows using Roe FDS and AUSM+. For this purpose we developed the frozen and the equilibrium code and numerically simulated the viscous interaction by changing the surface temperature and the mach number. We used curve fitting data in NASA Reference Publication 1181, 1260 to calculate equilibrium properties. We compare the equilibrium flow with the frozen flow. We conclude that the mach number and the surface temperature are significant parameters, as the surface temperature and the mach number increase the viscous interaction becomes stronger, and we must consider high-temperature effects in hypersonic flow

The flow field around a three dimensional minivan-like body has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using second-order accurate schemes for the time derivatives, and third/second-order scheme for the spatial derivatives. The Marker-and-Cell concept is applied to efficiently solve continuity equation. The fourth -order artificial damping is added to the continuity equation for numerical stability. A H-H type multi-block grid system is generated around a three dimensional minivan-like body. Turbulent flows have been modeled by the Baldwin-Lomax turbulent model. The simulation shows three dimensional vortex-pair just behind body. And the flow separation is also observed the rear of the body. It has concluded that the results of present study properly agree with physical flow phenomena.

Mixing is most important for developing an electric washer which transforms angular momentum from rotating solid wall to laundry clothes inside it. For magnification of this mixing effect, some inventions are introduced to washing machine system, i. e., washing plate, washing rod, and even for washing cap in a model of a Korean manufacture. However, the previous efforts show dissatisfaction up till now. In this paper, a triumph to enhance mixing effects to increase washing performance is presented and verified by numerical investigation. The present model to simulate a washing tub is the simple circular cylinder with two endwall disks which is completely filled with a viscous liquid. The present improvement is to change mounting position of a bottom disk of the model cylinder. Therefore, the aim of this work just proposes a new idea, which is numerically inspected, to a producer of washing machine, In detail, this invention is alternating the mounting position of a rotating bottom disk. Actually skewed pulsator is placed in steady of a flat disk, so the two endwall disks at top and bottom are not in parallel. The angle between an inclined bottom disk and the horizontal plane is fixed as 5 degree and physical domain to consider poses a sliced cylinder. Flow fields in both a right circular cylinder and the present improved model are fully depicted by numerical integration on a body fitted nonorthogonal regular uniform grid system. Numerical data to explain flow structure are plotted for understanding of the effects of the inclined disk. Also enhanced mixing effects by the inclined rotating disk are gauged by accurate numerical data used in this work.

This paper describes some computational results of various energy and environmental systems using Patankar's SIMPLE method. The specific topics handled in this study are jet bubbling reactor for flue gas desulfurization, cyclone-type afterburner for incineration, 200m tall stack for 500 MW electric power generation, double skin and heat storage systems of building energy saving for the utilization of solar heating, finally turbulent combustion systems with liquid droplet or pulverized coal particle. A control-volume based finite-difference method with the power-law scheme is employed for discretization. The pressure-velocity coupling is resolved by the use of the revised version of SIMPLE, that is, SIMPLEC. Reynolds stresses are closed using the standard $k-{\varepsilon}$ and RNG $k-{\varepsilon}$ models. Two-phase turbulent combustion of liquid drop or pulverized coal particle is modeled using locally-homogeneous, gas-phase, eddy breakup model. However simple approximate models are incorporated for the modeling of the second phase slip and retardation of ignition without consideration of any detailed particle behavior. Some important results are presented and discussed in a brief note. Especially, in order to make uniform exit flow for the jet bubbling reactor, a well-designed structure of distributor is needed. Further, the aspect ratio in the double skin system appears to be one of important factors to give rise to the visible change of the induced air flow rate. The computational tool employed in this study, in general, appears as a viable method for the design of various engineering system of interest.

The full incompressible Navier-Stokes equations are numerically integrated to solve the unsteady channel flow using a new numerical scheme of second-order accuracy developed by the authors. It is well known that in spite of the symmetry in the boundary condition and geometry, asymmetry can develop with time-dependency in a channel with sudden expansion. The instability of the shear flow and the cross-channel pressure contribute to such asymmetric flow. In this paper, we successfully generated a channel flow in which vortex waves were propagated downstream due to the harmonically oscillating inlet flow. The structure of the eddies and wall vorticity are parametrically investigated.

Characteristic boundary conditions are discussed in conjunction with a flux-difference splitting formulation as modified from Roe's linearization. Details of how one can implement the characteristic boundary conditions which are compatible with the discrete formulation at interior points are given for different types of boundaries including subsonic outflow and adiabatic wall. The latter conditions are demonstrated through computation of supersonic ogive-cylinder flow at high angle of attack and the computed wall pressure distribution is compared with experiment.

Numerical simulations are made for the three-dimensional flow around a wing in ground effect craft haying the complex geometry. A numerical tool is developed for the primary design of hull and wing shape of practical Wing-In-Ground effect(WIG) stop. The finite-difference method is utilized to descretize the governing equations and pressure field is obtained by using Marker-And-Cell(MAC) method. The air and water flows are simultaneously simulated in the time-marching solution procedure for the Navier-Stokes equation. The porosity technique and the density function are devised for the implementation of the three-dimensional body-boundary and the free-surface conditions, respectively. In this paper, a craft is modeled simply by three blocks containing a wing mounted on a main body horizontally, with the endplate. The numerical calculations of a WIG advancing in a calm water are performed and the WIG-generated wave profiles are also obtained. In the final paper, details of the numerical methods employed for the present study and calculated results are discussed.

파랑중을 항주하는 Wigley 선형주위의 점성유동장을 수치적으로 해석하였다. Navier-Stokes 방정식을 유한차분법을 이용하여 차분화하고 압력과 속도의 연성은 MAC방법으로 처리하였다. 격자계로는 O-H헝 격자계를 사용하였으며 정수중에서 항주하는 결과를 다른 실험결과 및 계산결과와 비교하여 계산정도를 검증하였으며, 파랑이 선체주위의 유동장에 미치는 영향에 관하여 고찰하였다.

The jet pump is being used in many fields for several purposes because of its simple construction and easy operation. The characteristics of the geometrical variables, pressure gradient and velocity distribution of the jet pump are studied using the CFD technique. The flow calculations through a bended nozzle. a mixing chamber and a venturi are presented and phenomenological aspects are discussed. This study solve 3-D steady incompressible Navier-Stokes equations using the Iterative time marching scheme. The governing equations are differenced with 1st-order accurate backward difference scheme for the time derivatives and 3rd-order accurate QUICK scheme for the convective terms. The Mark-and-cell concept was applied efficiently to solve continuity equation, which is differenced 2nd-order accurate central differenced scheme. The 4th-order artificial damping is added to the continuity equation for numerical stability. A O-type of grid system is generated inside a nozzle and venturi of the jet pump. It has concluded that the results of present study properly agree with physical flow phenomena.

In the present paper, flow fields around two dimensional single and two circular cylinders are analysed by a finite difference method. Navier-Stokes and the continuity equations an solved to simulate the flow fields. A overlapped grid system(the composite of a body boundary-fitted grid system near the body and a rectangular grid system for other flow fields) is used for this calculation. In the use of overlapped grid system, it is most significant thing to exchange the physical quantities from one grid system to the other one continuously, In this research, the linear interpolations of physical quantaties are done for this purpose in the overlapped region. The numerical calculations are carried out for the flows around a circular cylinder and two cylinders to verify the accuracy of present method. The flow fields around two cylinders facing the flow with side by side and tandem arrangement are analysed. The results are compared to other experimental and computational ones done in other single grid system.

The Delaunay triangulation technique was tested for complicated shapes of computational domain. While a simple geometry, both in topology and in geometry, was discretized well into triangular elements. a complex geometry often failed in triangularization. A complex geometry should be devided into smaller sub-domains whose shape is simple both topologically and geometrically. The present study developed the data structures not only for relationships among neibering elements but also for shape information, and coupled these into the Delaunay triangulation technique. This approach was able to enhance greatly the reliability of triangularization specially in complicated shapes of computational domains. The GUI (Graphic User Interface) and OOP (Object-Oriented Programming) were used in order to develop the user-friendly and efficient computer code.

When computing the flow around complex three dimensional configurations, the generation of the grid is the most time consuming. This paper presents a survey of the grid generation technique. The technique for treating problems arising in gird generation for complex geometries is addressed. A general procedure for generating boundary conforming grids is included

A numerical analysis is conducted on heat transfer and fluid flow of a plasma spraying process under the DC-RE hybrid electromagnetic field. Plasma flow is analyzed by using Eulerian approach and the equation of particle motion is simultaneously solved using a trajectory analysis with a lumped-heat-capacity model. Axisymmetric two dimensional electromagnetic fields governed by Maxwell's equations are solved based on a vector potential concept. The effects of the RF electromagnetic field on the temperature and velocity fields of the turbulent plasma flow are clarified. Control characteristics of phase changes and dispersed features of particles by applying the RF electromagnetic field are also clarified in an attempt to improve the plasma spraying process

The numerical method for the flow field of a gas atomization process is presented. For the analysis of the compressible supersonic jet flow of a gas. an axisymmetric Navier-Stokes equations are solved using a LU-factored upwind method. The MUSCL type TVD scheme is used for the discretization of inviscid flux, whereas Steger-Warming splitting and LU factorization is applied to the implicit operator. For the validation of the present method, we computed the flow field around the simple gas atomizer proposed by Issac. The numerical results has shown excellent agreement with the experimental data.

A new wall impaction model for spray and its assessment are described in this paper. The gas phase is modelled in terms of the Eulerian continuum conservation equations of mass, momentum, energy and fuel vapour fraction. The liquid phase is modelled following the discrete droplet model approach. The droplet parcel contains many thousands of drops assumed to have the same size, temperature and velocity components. The droplet parcel equations of trajectory, momentum, mass and energy are written in Lagrangian form. The new drop-wall interaction model is proposed, which is based on experimental investigations on individual drops, and it is applied for the general non-orthogonal gird. The model is then assessed through comparison with experiments over a wide range of test conditions of sprays. The results are in good agreement with experimental data.