• 대한의용생체공학회:의공학회지
• /
• 제13권2호
• /
• pp.87-96
• /
• 1992

In thIns paper, a numerical simulation of steady laminar and turbulent flow in a two dimensional model for the total artificial heart is'presented. A trlleaflet polyurethane valve was simulated at the outflow orifice while the Inflow orifice had a trileaflet or a flap valve. The finite analytic numerical method was employed to obtain solutions to the governing equations in the Cartesian coordinates. The closure for turbulence model was achieved by employing the k-$\varepsilon$-E model. The SIMPLER algo rithm was used to solve the problem in primitive variables. The numerical solutions of the slulated model show that regions of relative stasis and trapped vortices were smaller within the ventricular chamber with the flap valve at the Inflow orifice than that with the trileaflet valve. The predicted Reynolds stresses distal to the inflow valve within the ventricular chamber were also found to be smaller wlth the flap valve than with the trlleaflet valve. These resu1ts also suggest a correlation be- tween high turbulent stresses and the presence of thrombus In the vicinity of the valves in the total artificial hearts. The computed velocity vectors and trubulent stresses were comparable with previ ously reported in vitro measurements in artificial heart chambers. Analysis of the numerical solo talons suggests that geometries similar to the flap valve(or a tilting disc valve) results in a better flow dynamics within the total artificial heart chamber compared to a trileaflet valve.

• 동력기계공학회지
• /
• 제22권1호
• /
• pp.11-17
• /
• 2018

The steady-state, incompressible and three-dimensional numerical analysis was carried out to evaluate the internal flow fields characteristics of wind tunnel contractions made by Morel's curve equations. The turbulence model used in this study is a realizable ${\kappa}-{\varepsilon}$ well known to be excellent for predicting the performance of the flow separation and recirculation flow as well as the boundary layer with rotation and strong back pressure gradient. As a results, when the flow passes through the interior space of the analytical models, the flow resistance at the inlet of the plenum chamber is the largest at $Z_m=300$, 400 mm, but the smallest at $Z_m=700mm$. The maximum turbulence intensity in the test section is about 2.5% when calculated by the homogeneous flow, so it is improved by about 75% compared to the 10% turbulence intensity at the inlet of the plenum chamber due to the contraction.

• Wind and Structures
• /
• 제9권1호
• /
• pp.37-58
• /
• 2006

The paper describes a study about effects of upstream hills on design wind loads using two mathematical approaches: Computational Fluid Dynamics (CFD) and Artificial Neural Network (NN for short). For this purpose CFD and NN tools have been developed using an object-oriented approach and C++ programming language. The CFD tool consists of solving the Reynolds time-averaged Navier-Stokes equations and $k-{\varepsilon}$ turbulence model using body-fitted nearly-orthogonal coordinate system. Subsequently, design wind load parameters such as speed-up ratio values have been generated for a wide spectrum of two-dimensional hill geometries that includes isolated and multiple steep and shallow hills. Ground roughness effect has also been considered. Such CFD solutions, however, normally require among other things ample computational time, background knowledge and high-capacity hardware. To assist the enduser, an easier, faster and more inexpensive NN model trained with the CFD-generated data is proposed in this paper. Prior to using the CFD data for training purposes, extensive validation work has been carried out by comparing with boundary layer wind tunnel (BLWT) data. The CFD trained NN (CFD-NN) has produced speed-up ratio values for cases such as multiple hills that are not covered by wind design standards such as the Commentaries of the National Building Code of Canada (1995). The CFD-NN results compare well with BLWT data available in literature and the proposed approach requires fewer resources compared to running BLWT experiments.

• 한국전산유체공학회지
• /
• 제13권4호
• /
• pp.1-7
• /
• 2008

In the process of continuous hot-dip galvanizing, it is well known that the gas wiping through an air knife system is most effective because of its uniformity in coating thickness, possibility of thin coating, workability in high speed, and simplicity of control. However, gas wiping used in the galvanizing process brings about a problem of splashing at the strip edge above a certain high speed of process. It is also known that the problem of edge splashing is more harmful than that at the mid strip surface. For a given liquid(of a certain viscosity and surface tension), the onset of splashing mainly depends upon the strip velocity, the gas-jet pressure, and the nozzle's stand-off distance. In these connections in the present study, we proposed three kinds of air knife system having nozzles of constant expansion rate, and compared the jet structures issuing from newly proposed nozzle systems with the result by a conventional one. In numerical analysis, the governing equations are consisted of two-dimensional time dependent Navier-Stokes equations, and the standard k-${\varepsilon}$ turbulence model is employed to solve turbulence stress and so on. As the result, it is found that we had better use the constant expansion-rate nozzle which can be interpreted from the point view of the energy saving for the same coating thickness. Also, we better reduce the size of separation bubble and enhance the cutting ability at the strip surface, by using an air-knife having constant expansion-rate nozzle.

• Wind and Structures
• /
• 제14권5호
• /
• pp.435-447
• /
• 2011

A Computational Fluid Dynamics model is presented in this study for the simulation of the complex fluid flows with free surfaces inside the Tuned Liquid Column Dampers in horizontal motion. The characteristics of the fluid model of the TLCD in horizontal motion include the free surface of the multiphase flow and the horizontal moving frame. In this study, the time depend unsteady Standard ${\kappa}-{\varepsilon}$ turbulent model based on Navier-Stokes equations is chosen. The volume of fluid (VOF) method and sliding mesh technique are adopted to track the free surface of water inside the vertical columns of TLCD and treat the moving boundary of the walls of TLCD in horizontal motion. Several model solution parameters comprising different time steps, mesh sizes, convergence criteria and discretization schemes are examined to establish model parametric independency results. The simulation results are compared with the experimental data in the dimensionless amplitude of the water column in four different configured groups of TLCDs with four different orifice areas. The predicted natural frequencies and the head loss coefficient of TLCDs from CFD model are also compared with the experimental data. The predicted numerical results agree well with the available experimental data.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 추계학술대회
• /
• pp.1482-1487
• /
• 2004

In general the swirl jet is generated by the injected flow that is forced to the tangential direction. A spiral nozzle which is composed of an annular slit and a convergent nozzle, is released the spiral jet that is generated by the radial flow injection through an annular slit. The objective of the present study is to investigate the additional study that is studied a changed the convergent nozzle angle and nozzle length. In the present computation, a finite volume scheme is used to solve three dimensional Navier-Stokes equations with RNG $k-{\varepsilon}$ turbulent model. The convergent nozzle angle and the nozzle length of the spiral nozzle are varied to obtain different spiral flows inside the conical convergent nozzle. The present computational results are compared with the previous experimental data. The results obtained show that the convergent nozzle angle and the nozzle length of the spiral jet strongly influence the characteristics of the spiral jets, such as a tangential and a jet width.

• 대한기계학회논문집B
• /
• 제24권10호
• /
• pp.1368-1379
• /
• 2000

A numerical study on swirling pulverized coal combustion in an axisymmetric enclosure is carried out by applying the 2-phase weighted sum of gray gases model (WSGGM) approach with the discrete ordinate method (DOM) to model the radiative heat transfer equation. In the radiative transfer equation, the same polynomial equation and coefficients for weighting factors as those for gas are adopted for the coal/char particles as a function of partial pressure and particle temperature. The Eulerian balance equations for mass, momentum, energy, and species mass fractions are adopted with the standard and RNG k-${\varepsilon}$ turbulence model, whereas the Lagrangian approach is used for the particulate phase. The eddy-dissipation model is employed for the reaction rate for gaseous mixture, and the single-step and two-step first-order reaction model for the devolatilization process for coal. Special attention is given to establish the thermal boundary conditions on radiative transfer equation By comparing the numerical results with experimental ones, the radiation model used here is confirmed and found to provide an alternative for simulating the radiative transfer.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2007년도 춘계학술대회논문집
• /
• pp.275-282
• /
• 2007

In this study, a fluid-structure interaction (FSI) analysis system has been developed in order to evaluate the turbine cascade performance with blade structural deformation effect. Relative movement of the rotor with respect to stator is reflected by modeling independent two computational domains. To consider the deformed position of rotor airfoil, dynamic moving grid method is applied. Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation SST $k-{\varepsilon}$ turbulence models are solved to predict unsteady fluid dynamic loads. A fully implicit time marching scheme based on the Newmark direct integration method with high artificial damping is used to compute the fluid-structure interaction problem. Cascade performance evaluations for different elastic axis positions are presented and compared each other. It is importantly shown that the predicted aerodynamic performance considering structural deformation effect of blade can show some deviations compared to the data generally computed from rigid blade configurations and the position of elastic axis also tend to give sensitive effect.

• 한국자동차공학회논문집
• /
• 제11권3호
• /
• pp.115-123
• /
• 2003

With the advent of high performance computers and more efficient numerical algorithms, computational fluid dynamics(CFD) has come out as a modem alternative for reducing the use of wind tunnels test in automotive engineering. However, in spite of the fact that many competent researchers have made all their talents in developing turbulence model over since the past dozen or more years, it has been an important impediment in using the CFD effectively to design machinery and to diagnose or to improve engineering problems in the industry since the turbulence model has been acting as the Achilles' tendon in aspect of the reliability even to this time. In this study, Reynolds-averaged Wavier-Stokes equations were solved to simulate an incompressible turbulent flow around a bus-like bluff body near ground plane. In order to investigate the effect of the discretisation schemes and turbulence model on the aerodynamic forces several turbulence models with five convective difference schemes are adopted. From the results of this study, it is clear that choice of turbulence model and discretisation scheme profoundly affects the computational outcome. The results also show that the adoption of RNG $k-\varepsilon$ turbulence model and nonlinear quadratic turbulence model with the second order accurate discretisation scheme predicts fairly well the aerodynamic coefficients.

• 대한기계학회논문집
• /
• 제18권7호
• /
• pp.1851-1865
• /
• 1994

The objective of the present study is to predict the film cooling effectiveness by a row of holes at various injection ratios and injection angles. Numerical calculations have been performed to investigate the characteristics of flow and temperature distributions in a region near the down-stream of injection hole including the region of adverse pressure gradient. The elliptic turbulent 3-dimensional governing equations with variable thermal properties using the low-Reynolds number k-$\bar{varepsilon}$ model was solved by SIMPLE algorithm. The results showed that the presence of adverse pressure gradient and secondary vortex in the region near the downstream of injection hole induces large temperature gradent. The $45^{\circ}$ injection has higher averaged film cooling effectiveness than $60^{\circ}$ injection. But neverthless the $90^{\circ}$ injection has greater deviation from a flat plate than $45^{\circ}$ and $60^{\circ}$ injection, the $90^{\circ}$ injection has higher averaged film cooling effectiveness than $45^{\circ}$ and $60^{\circ}$ injection in the region near the downstream of injection hole.