• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.10
• /
• pp.1384-1391
• /
• 2001

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67) and to evaluate the performances of Abid's low-Reynolds-number k-$\omega$ and Baldwin-Lomax turbulence models. A finite volume method is used fur spatial discretization. The equations are solved implicitly in time by the use of approximate factorization. The upwind difference scheme is used for inviscid terms and viscous terms are approximated with central difference. The flux-difference-splitting method of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. 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, it is concluded that Abid'k-$\omega$ 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 the same.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.3
• /
• pp.201-209
• /
• 2011

An unsteady RANS analysis study of the 3-D flow around the NREL Phase VI horizontal axis wind turbine(HAWT) was performed using sliding mesh approach. Two different analysis models such as rotor-only and rotor with tower/nacelle were constructed to investigate the blade/tower interaction. Analysis results for the rotor with tower/nacelle were compared with the corresponding NREL's experimental data which produced fairly good validation of the present CFD model. Comparison of flows around those two models also clearly showed the blade/tower interaction even it was small for upwind configuration. Other visualization results and integrated aerodynamic loads including torque of the blade demonstrated the effective unsteady flow simulation capability of the present CFD model.

• Journal of Bio-Environment Control
• /
• v.3 no.1
• /
• pp.10-19
• /
• 1994

Wind load is known to be one of major forces to influence the stability of agricultural structures. General flow fields were calculated to determine flow characteristics over the envelop of the following three types of greenhouses with arched roof : single span, twin span greenhouses, and two single span greenhouses apart 3m inbetween. Pressure coefficients along the envelop of greenhouse were numerically calculated by the k-$\varepsilon$ turbulence model, which lead to determine wind forces on it. Curvilinear coordinate for an arched roof and the upwind scheme were adopted for the study. The calculated pressure coefficients were validated with the avaliable data of Japanese Standard and NGAM Standard. The Magnitude of calculated forces over the envelop was not in good accordance with data except the windward wall. Even tile data of Japanese and NGAM Standard for validation deviated a lot from each other in quantity and quality. Such discrepancy may be attributed to different geometric and/or flow configuration conditions for experiments, or the insenstivity of the k-$\varepsilon$ turbulence model to recirculation flow.

• Coupled systems mechanics
• /
• v.8 no.2
• /
• pp.129-146
• /
• 2019

The main goal of this work is to develop a numerical simulator to study an isothermal single-phase two-component flow in a naturally fractured oil reservoir, taking into account advection and diffusion effects. We use the Peng-Robinson equation of state with a volume translation to evaluate the properties of the components, and the discretization of the governing partial differential equations is carried out using the Finite Difference Method, along with implicit and first-order upwind schemes. This process leads to a coupled non-linear algebraic system for the unknowns pressure and molar fractions. After a linearization and the use of an operator splitting, the Conjugate Gradient and Bi-conjugated Gradient Stabilized methods are then used to solve two algebraic subsystems, one for the pressure and another for the molar fraction. We studied the effects of fractures in both the flow field and mass transport, as well as in computing time, and the results show that the fractures affect, as expected, the flow creating a thin preferential path for the mass transport.

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

Computations of the mean and turbulence flows over three-dimensional hill of conical shape have implemented. Beside the standard ${\kappa}-{\varepsilon}$ , two other modifications proposed by Detering & Etling and Duynkerke for atmospheric applications were also considered. These predictions were compared with the data of a wind tunnel experiment. From the comparison, it was concluded that all three models predict the mean flow velocities equally well while only the Duynkerke's model accurately predicts the turbulence data statistics. It also concluded that there are large discrepancies between model predictions and the measurements near the ground surface. The flow field, which was obtained by using the Duynkerke's modification, was used to simulate gas dispersion from an upwind source. The calculation results are verified based on the measurement data. Modifications of the turbulent Schmidt number were carried out in order to match the measured results. The code was used to investigate the influence of the recirculation zone behind a building of cubical shape on the transport and dispersion of pollutant. For a stack behind and near the obstacle, some conclusions about the effect of the stack height and stack location were derived.

• Journal of computational fluids engineering
• /
• v.5 no.2
• /
• pp.38-46
• /
• 2000

It has been indicated that the rotor/stator interaction has distinct causes of unsteadiness, such as the viscous vortex shedding, wake/stator interaction and potential rotor/stator interaction. In this paper, the mechanism of unsteady potential interaction in one stage axial compressor is numerically investigated for blade row ratio 1:1 and 2:3 at design point and for blade row ratio 2:3 at off-design point in two-dimensional view point. The numerical technique used is the upwind scheme of Van-Leer's Flux Vector Splitting(FVS) and Cubic spline interpolation is applied on zonal interface. In this study the flow unsteadiness due to potential interaction are found to be larger in blade row ratio 2:3 than in 1:1. The total pressure rise in blade row ratio 2:3 is closer to the real value in design point than that in 1:1. The change of unsteady pressure amplitude according to the variation of stator exit pressure is very small.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.11
• /
• pp.1518-1526
• /
• 2001

In this paper, the mechanism of unsteady potential interaction and wake interaction in one stage axial turbine is numerically investigated at design point in two-dimensional viewpoint. The numerical technique used is the upwind scheme of Van-Leer's Flux Vector Splitting (FVS) and Cubic spline interpolation is applied on zonal interface between stator and rotor. The inviscid analysis is used to embody the influence of potential interaction only and viscous analysis is used to embody the influences of both potential interaction and wake interaction at the same time. The potential-flow disturbance from the stator into a rotor passage and the periodic blockage effect of rotor produce the unsteady pressure on the blade surface in inviscid analysis. After the wake is cut by rotor, two counterrotating votical patterns flanking the wake centerline in the passage are generated. So, these phenomena magnify the unsteady pressure in viscous analysis than that in inviscid analysis. The resulting unsteady forces on the rotor, generated by the combined interaction of the two effects by potential and wake interaction, are discussed.

• Proceedings of the SAREK Conference
• /
• 2008.06a
• /
• pp.1366-1371
• /
• 2008

Four simulation models of plastic plate heat exchangers are designed and simulated. The flat plate type heat exchanger is designed as the reference model in order to evaluate how much thermal performance increases. The turbulence promoter type heat exchanger is fabricated with cylindrical-type vortex generators and rib-type turbulence promoters. The corrugate type is obtained from the conventional stainless steel compact heat exchangers, which are called the herringbone-type compact heat exchangers. The dimple type heat exchanger has a number of dimples on its surface. In this study, the flow and heat transfer characteristics of the plastic plate heat exchanger are investigated using numerical simulation and compared with experimental results. The flows are assumed as a three-dimensional, incompressible and turbulent model. The standard k-$\varepsilon$ model is used as the turbulent flow modeling, the SIMPLE algorithm is used to treat the coupling between pressure and velocity, and first order upwind scheme is used for discretization of momentum, turbulent and energy. The computational analysis and experimental results both show that the friction coefficient and Nu number is highest in the corrugate type.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.144-151
• /
• 2003

Computations of the mean and turbulence flows over three-dimensional hill of conical shape have implemented. Beside the standard $\kappa-\epsilon$, two other modifications proposed by Detering & Etling and Duynkerke for atmospheric applications were also considered. These predictions were compared with the data of a wind tunnel experiment. From the comparison, it was concluded that all three models predict the mean flow velocities equally well while only the Duynkerke's model accurately predicts the turbulence data statistics. It also concluded that there are large discrepancies between model predictions and the measurements near the ground surface. The flow field, which was obtained by using the Duynkerke's modification, was used to simulate gas dispersion from an upwind source. The calculation results are verified based on the measurement data. Modifications of the turbulent Schmidt number were carried out in order to match the measured results. The code was used to investigate the influence of the recirculation zone behind a building of cubical shape on the transport and dispersion of pollutant. For a stack behind and near the obstacle, some conclusions about the effect of the stack height and stack location were derive

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.10
• /
• pp.754-760
• /
• 2008

Computation of moving interface by the level set method typically requires the reinitialization of level set function. An inaccurate estimation of level set function $\phi$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, an accurate and robust reinitialization process is essential to the simulation of free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates level set function directly using a normal vector on the interface without solving there-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1 splitting/SUPG (Streamline Upwind Petrov-Galerkin) FEM are adopted to discretize advection equation of the level set function and the incompressible Navier-Stokes equation, respectively. Advection equation and re-initialization process of free surface capturing are validated with benchmark problems, i.e., a broken dam flow and timereversed single vortex flow. The simulation results are in good agreement with the existing results.