• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.84-92
• /
• 1997

In this paper a new structured grid generation technique is introduced. This new technique utilizes predictor-corrector approach, and is a marching scheme in the global sense as the hyperbolic scheme is. In the predictor step, one layer of grid cells is obtained by using Modified Advancing Front Method which generates a collection of quadrilateral cells simultaneously. In the corrector step, the layer of grid cells that is calculated in the predictor step is adjusted by solving Laplace equations to prevent grid lines from skewing and overlapping in highly curved configurations. It is shown that the resultant algorithm, named a MAP scheme, which combines the Modified Advancing Front Method as a Predictor with an elliptic scheme as a corrector can be used to generate globally smooth and locally near-orthogonal grids for external flow fields even for highly curved configurations. Examples of grid generations for external flow fields about several configurations by use of the present approach are given, and its applicability and flexibility have been demonstrated and discussed.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.117-128
• /
• 1997

Effects of dam and weir on the fluid flow and behavior of inclusions in the continuous casting tundish have been studied using the CFD (Computational Fluid Dynamics) technique. Inclusions affecting the quality of steel products have been considered to be passive: the fluid flow has been obtained for unstaggered grid points defined on body-fitted generalized cuvilinear coordinates with no attention on inclusions, and the spatial propagation of inclusions has been determined by using the flow field data. The result show that the dam and weir direct the flow to the free surface and increase the residence time of inclusions significantly, and thereby that inclusions have much more chance to be floated to the free surface of the tundish where it is eliminated. It is also found that they offer more margin on the geometric design of exit nozzles connected to moulds. This finding is particularly important for twin casting operations where the quality of steel products from the two moulds be kept uniformly.

• Journal of computational fluids engineering
• /
• v.20 no.1
• /
• pp.77-83
• /
• 2015

Flow characteristics of synthetic jet based flow supplying devices have been computationally investigated for different device shapes. Jet momentum was produced by the volume change of a cavity by two piezoelectric-driven diaphragms. The devices have additional flow path compared with the original synthetic jet actuator, and these flow path changes the flow characteristics of synthetic jet actuator. Four non-dimensional parameters, which were functions of the shapes of the additional flow path, were considered as the most critical parameters in jet performance. Comparative studies were conducted to compare volume flow rate and jet velocity. Computed results were solved by 2-D incompressible Navier-Stokes solver with k-w SST turbulence model. Detailed computations revealed that the additional flow path diminishes suction strength of the synthetic jet actuator. In addition, the cross section area of the flow path has more influence over the jet performances than the length of the flow path. Based on the computational results, the synthetic jet based flow supplying devices could be improved by applying suitable shape of the flow path.

• Journal of computational fluids engineering
• /
• v.20 no.1
• /
• pp.47-56
• /
• 2015

In this paper, we study the effect of various turbulence models by comparing the aerodynamic characteristics and the flow patterns computed for aircraft models. An in-house CFD solver, MSAPv, that solves the three dimensional RANS equations with the turbulence model equations is used. The turbulence models used in this study are the Spalart-Allmaras model, Menter's $k-{\omega}$ SST model, Coakley's $q-{\omega}$ model, and Huang and Coakley's $k-{\varepsilon}$ model. DLR-F6 WB and WBNP configurations are selected for the study. We concentrate on the separated flow pattern variations with the turbulence models at the wing-body junction and the wing-pylon junction as well as drag polar curves.

• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.54-60
• /
• 2012

Transition prediction results are validated with experimental data obtained from a transonic wind tunnel for the KU109C airfoil. A Reynolds-Averaged Navier-Stokes code is simultaneously coupled with the transition transport model of Langtry and Menter and applied to the numerical prediction of aerodynamic performance of the KU109C airfoil. Drag coefficients from the experiment are better correlated to the numerical prediction results using a transition transport model rather than the fully turbulent simulation results. Maximum lift coefficient and drag divergence at the zero-lift condition with Mach number are investigated. Through the present validation procedure, the accuracy and usefulness of both the experiment and the numerical prediction are assessed.

• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.92-98
• /
• 2010

In a modern high speed drawing process of optical fibers, it is necessary to use helium as a cooling gas in a glass fiber cooling unit in order to sufficiently cool down the fast moving glass fiber freshly drawn from the heated silica preform in the furnace. Since the air is entrained unavoidably when the glass fiber passes through the cooling unit, the helium is needed to be injected constantly into the cooling unit. The present numerical study investigates and analyzes the air entrainment using an axisymmetric geometry of glass fiber cooling unit. The effects of helium injection rate and direction on the air entrainment rate are discussed in terms of helium purity of cooling gas inside the cooling unit. For a given rate of helium injection, it is found that there exists a certain drawing speed that results in sudden increase in the air entrainment rate, which leads to the decreasing helium purity and therefore the cooling performance of the glass fiber cooling unit. Also, the helium injection in aiding direction is found to be more advantageous than the injection in opposing direction.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.78-84
• /
• 2013

In this paper, numerical simulations of both low- and high-frequency buzz phenomena at the throttle ratios (T.R.) in Nagashima's experiment are performed. The dominant frequencies of the low-and high-frequency buzz in the experiment are about 109 Hz with T.R.=0.97 and 376 Hz with T.R.=0.55, respectively. An axisymmetric solver with the S-A turbulence model is used for the simulations, and DFT(Discrete Fourier Transform) on pressure histories is conducted for the buzz frequency analysis. In the present simulations, the free-stream Mach number and the Reynolds number based on the inlet diameter are 2 and $10^7$, respectively. Both the low- and high-frequency buzz phenomena are accomplished without the changes in the grid topology. The dominant frequency of the simulation is about 125 Hz with T.R.=0.97, while it is 399 Hz with T.R.=0.55.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.49-55
• /
• 2013

The goal of the research is to evaluate the open source code of OpenFOAM for the use of nuclear plant flow simulation objectively. Of the various incompressible flow solvers, simpleFoam, pimpelFoam are then tested under three validated cases (backward facing step, flow over circular cylinder and turbulent round jet flow). For the evaluation of steady state incompressible laminar flow simulation, low reynolds number of backward facing step flow was solved by simpleFoam. The resultant of the reattached lengths turned out to be similar with the other experimental and simulation results. For transient flow simulation, flow over circular cylinder and turbulent round jet flow were solved by pimpleFoam. The simulation accuracy was evaluated by comparing the resultant flow patterns with the description of the characteristics of the flow over the circular cylinder. The quantitative accuracy was evaluated for no more than 85% by comparing it to the decaying constants of the turbulent round jet velocity.

• Journal of computational fluids engineering
• /
• v.5 no.1
• /
• pp.33-42
• /
• 2000

The LiBr-H₂O absorption process on a horizontal tube has been analyzed numerically. The flow field, which was calculated in the authors' previous study by solving the fully elliptic Navier-Stokes equations with accurate free-surface-tracking method, is used to solve the temperature and concentration distributions in the absorption film. With the assumption that the absorbent is linear, calculations have been made for various inlet temperature and flow-rate conditions. For low inlet temperature, the absorption rate is large in the upstream region but the mean temperature also increases and as a result the absorption decreases as the film flows to downstream while high-inlet-temperature case does the opposite. The difference in the absorption rate due to the inlet temperature change becomes smaller in the downstream than that in the upstream. For large flow rate, the heat transfer to the wall becomes poor due to the thick film and so does the absorption rate. The analyses have also been carried out for multiple tube arrangement and the results show that the absorption rate converges after a few tube rows.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.24-30
• /
• 2016

In the water purification plant, the mixture of water and chemical from the mixing basin enters the flocculation basin. The rotating flocculators are generally used for the efficient flocculation of dregs. In this paper, the performance of flocculators of a vertical paddle type, widely used in the typical flocculation basins, and a hydro-foil type, recently disseminated in the field, are compared with each other by use of the numerical method. Also the characteristics and the efficiency are analyzed with CFD techniques. The strain rate and the eddy viscosity are compared for two types to predict the mixing efficiency, and the maximum speed and its location are pursued from the computed data. The hydrofoil type shows that the eddy viscosity is enhanced 1.66 to 3.03 times larger than that of vertical paddle type, and also produced 1.87 to 1.95 times larger flocs for each stage. However, the rapid rotation of hydrofoil may chop the floc to small size due to the higher turbulence intensity. From the result of computation, the strong and weak points of each type have been analyzed for the decision making.