• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.90-97
• /
• 1999

An implicit turbulent flow solver is developed for 2-D unstructured hybrid meshes. Spatial discretization is accomplished by a cell-centered finite volume formulation using an upwind flux differencing. Time is advanced by an implicit backward Euler time stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one equation model, which is coupled with wall function. The numerical method is applied for flows on a flat plate, the NACA 0012 airfoil, and the Douglas 3 element airfoil. The results are compared with experimental data.

• Journal of Mechanical Science and Technology
• /
• v.14 no.9
• /
• pp.985-996
• /
• 2000

The issue of quench is related to safety operation of large-scale superconducting magnet system fabricated by cable-in-conduit conductor. A numerical method is presented to simulate the thermal hydraulic quench characteristics in the superconducting Tokamak magnet system, One-dimensional fluid dynamic equations for supercritical helium and the equation of heat conduction for the conduit are used to describe the thermal hydraulic characteristics in the cable-in-conduit conductor. The high heat transfer approximation between supercritical helium and superconducting strands is taken into account due to strong heating induced flow of supercritical helium. The fully implicit time integration of upwind scheme for finite volume method is utilized to discretize the equations on the staggered mesh. The scheme of a new adaptive mesh is proposed for the moving boundary problem and the time term is discretized by the-implicit scheme. It remarkably reduces the CPU time by local linearization of coefficient and the compressible storage of the large sparse matrix of discretized equations. The discretized equations are solved by the IMSL. The numerical implement is discussed in detail. The validation of this method is demonstrated by comparison of the numerical results with those of the SARUMAN and the QUENCHER and experimental measurements.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.327-330
• /
• 2011

CFD simulation is widely used in various industries, universities and research centers. In Korea most of the researchers use foreign commercial S/W packages especially in industries. But commercial CFD packages have some problems as limit to source code and very high license foe. So from several years ago open source CFD code has been widely spread as an alternative. But in Korea there are a few users of open source code. Insufficiency of performance validation as for accuracy, robustness, convenience and parallel speed-up is important obstacles of open source code. So we tested some validation cases as to incompressible external aerodynamics and internal flaws and now are doing compressible flaws. As the first stage of compressible flow validation, we simulated Korea next generation high speed train(HEMU). It's running condition is 400km/hr and maximum Mach number reaches up to 0.4. With the high speed train we tested accuracy, robustness and parallel performance of open source CFD code OpenFOAM Because there isn't experimental data we compared results with widely used commercial code. When use $1^{st}$ order upwind scheme aerodynamic forces are very similar to commercial code. But using $2^{nd}$ order upwind scheme there was some discrepancy. The reason of the difference is not clear yet. Mesh manipulation, domain decomposition, post-processing and robustness are satisfactory. Paralle lperformance is similar to commercial code.

• Journal of Environmental Science International
• /
• v.7 no.3
• /
• pp.349-360
• /
• 1998

Effects of uniform flow on a two-dimensional mesoscale horizontal convection were investigated by using the vorticity and thermodynamic equations. For thins purpose, We simulated properties of a thermal convection m a stably stratified Boussinesq flued caused by partial heating at the center of a lower boundary If we don't consider effects of the uniform flow, the convection takes the form of aidsymmetrlc with respect to the z-alds. But when uniform flow Is strong, velocity field and temperature field consist of a sin91e cell structure which spreads upstream side of the partial heating area. The flow pattern for strong uniform flows takes the form of positive temperatue near the ground and negative temperature perturbation soft over the partial heating area, and downward motion directly over the upwind portion of the partial heating area and upward motion on the downstream side. The downstream edge of the upstream cell Is shifted in the downstream direction with the Increase of uniform flow almost linearly.

• 유체기계공업학회:학술대회논문집
• /
• 2002.12a
• /
• pp.337-345
• /
• 2002

The present study addresses a computational result of unsteady gas flow through a critical nozzle. The axisymmetric, unsteady, compressible, Wavier-Stokes equations are solved using a finite volume method that makes use of the second order upwind scheme for spatial derivatives and the multi-stage Runge-Kutta integral scheme for time derivatives. The steady solutions of the governing equation system are validated with the previous experimental data to ensure that the present computational method is valid to predict the critical nozzle flows. In order to simulate the effects of back pressure fluctuations on the critical nozzle flows, an excited pressure oscillation with an amplitude and frequency is assumed downstream of the exit of the critical nozzle. The results obtained show that for low Reynolds numbers, the unsteady effects of the pressure fluctuations can propagate upstream of the throat of critical nozzle, and thus giving rise to the applicable fluctuations in mass flow rate through the critical nozzle, while for high Reynolds numbers, the pressure signals occurring at the exit of the critical nozzle do not propagate upstream beyond the nozzle throat. For very low Reynolds number, it is found that the sonic line near the throat of the critical nozzle remarkably fluctuateswith time, providing an important mechanism for pressure signals to propagate upstream of the nozzle throat, even in choked flow conditions. The present study is the first investigation to clarify the unsteady effects on the critical nozzle flows.

• 유체기계공업학회:학술대회논문집
• /
• 1998.12a
• /
• pp.42-48
• /
• 1998

Numerical study is presented for the analysis of three-dimensional incompressible turbulent flows in multiblade centrifugal fan. Reynolds-averaged Navier-Stokes equations with standard k - $\epsilon$ turbulence model are transformed to non-orthogonal curvilinear coordinates, and are discretized with finite volume approximations. Linear Upwind Differencing Scheme(LUDS) is used to approximate the convection terms in the governing equations. SIMPLEC algorithm is used as a velocity-pressure correction procedure. The computational area is divided into three blocks; core, impeller and scroll, which are linked by multi-block method. The flow inside of the fan is regarded as steady flow, and mathematical formula established from the cascade theory and empirical coefficient are employed to simulate tile flow through the impeller. From comparisons between the computational results and the experimental data, the validity of the mathematical formula for the blade forces was examined and good results were obtained qualitatively. Hence, we can get the flow characteristics of multi-blade centrifugal fan and it will be a corner stone of the development of the multiblade centrifugal fan.

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

A FA-FD hybrid method, developed for solving three-dimensional incompressible Navier-Stokes equations, is applied to calculate three-dimensional laminar flows through a square duct with a 90° bend. The method discretizes the convective terms in the primary flow direction with 3rd-order upwind finite-differences and the convective and diffusive terms in the transverse directions with the two-dimensional finite analytic method. The non-staggered grid system is used and the pressure-velocity coupling is achieved by a global iteration procedure based on the PISO algorithm. Detailed comparisons between the computed solutions and the available experimental data are given mainly for the velocity distributions at cross-sections in a 90° bend of a square duct with both fully developed and developing entry flows. Although the computational result shows generally a good agreement with the experimental data, there are some significant discrepancies underlining the necessity of more accurate numerical methods as well as reliable experimental data for their validation.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.386-391
• /
• 2001

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.

• Wind and Structures
• /
• v.20 no.3
• /
• pp.449-468
• /
• 2015

Current wind-resistance designs of large-scale indirect dry cooling towers (IDCTs) exclude an important factor: the influence of the ventilation rate for radiator shutter on wind loads on the outer surfaces of the tower shell. More seemingly overlooked aspects are the effects of various ventilation rates on the wind pressure distribution on the tower surfaces of two IDCTs, and the feature of the flow field around them. In order to investigate the effects of the radiator shutter ventilation rates on the aerodynamic interference between IDCTs, this paper established the numerical wind tunnel model based on the Computational Fluid Dynamic (CFD) technology, and analyzed the influences of various radiator shutter ventilation rates on the aerodynamic loads acting upon a single and two extra-large IDCTs during building, installation, and operation stages. Through the comparison with the results of physical wind tunnel test and different design codes, the results indicated that: the influence of the ventilation rate on the flow field and shape coefficients on the outer surface of a single IDCT is weak, and the curve of mean shape coefficients is close to the reference curve provided by the current design code. In a two-tower combination, the ventilation rate significantly affects the downwind surface of the front tower and the upwind surface of the back tower, and the larger positive pressure shifts down along the upwind surface of the back tower as the ventilation rate increases. The ventilation rate significantly influences the drag force coefficient of the back tower in a two-tower combination, the drag force coefficient increases with the ventilation rate and reaches the maximum in a building status of full ventilation, and the maximum drag coefficient is 11% greater than that with complete closure.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.8
• /
• pp.2561-2571
• /
• 1996

The flow in the turbomachinery is very unsteady due to the stator-rotor interaction. It has been indicated that the stator-rotor interaction has three distinct causes of unsteadiness: that is, the viscous vortex shedding, wake rotor interaction and potential stator-rotor interaction. In this paper, the mechanism of unsteady potential interaction and wake interaction in the stator-rotor stage flow is numerically investigated 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. Then, the flow field of a compressor stage composed of NACA 65410 is analyzed. Flow fields are found to be simulated reasonably by this method and the sensitivity due to back-pressure variation is more stronger than rotor-velocity variation.