• Coupled systems mechanics
• /
• v.11 no.2
• /
• pp.107-119
• /
• 2022

The aim of the work presented in this paper is development of numerical model for prediction of temperature distribution in pavement according to the measured meteorological parameters, with introduction of non-linear heat transfer coefficient which is a function of temerature difference between the air and the pavement. Developed model calculates heat radiated from the pavement back in the air, which is an important part of the heat trasfer process in the open air surfaces. Temperature of the pavement surface, heat radiation together with many meteorological parameters were measured in series during two years in order to validate the model and calibrate model parameters. Special finite element method for temperature heat transfer towards the soil together with the time integration scheme are used to solve the governing equation. It is proved that non-linear heat transfer coefficient, which is a function of time and temperature difference between the air and the pavement, is required to decribe this phenomena. Proposed model includes heat tranfer coefficient callibration for specific climate region, through the iterative inverse procedure.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.53-61
• /
• 1999

Two-dimensional, unsteady, incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. The compressible code involves a conventional upwind-differenced scheme for the convective terms and LU-SGS scheme for temporal integration. The incompressible code with pseudo-compressibility method also adopts the same schemes as the compressible code. Three two-equation turbulence models are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by predicting the flow around the RAE 2822 transonic airfoil and the NACA 4412 airfoil, respectively. In addition, both the incompressible and compressible code are used to compute the flow over the NLR 7301 airfoil with flap to study the compressible effect near the high-loaded leading edge. The grid systems are efficiently generated using Chimera overlapping grid scheme. Overall, the κ-ω SST model shows closer agreement with experiment results, especially in the prediction of adverse pressure gradient region on the suction surfaces of high-lift airfoils.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.11
• /
• pp.1371-1378
• /
• 1999

Flow through turbomachinery has a very complex structure and Is Intrinsically unsteady. In addition, trend to highly loaded turbomachinery makes the flow extremely complex due to the interaction between rotor and stator. In this study, flows through UTRC LSRR turbine are numerically analyzed using 2 dimensional Navier-Stokes equations. The convective terms of the governing equations are discretized using the Van-Leer's FVS(Flux vector splitting) with an upwind TVD scheme. The conventional central differencing is used to discretize the diffusion terms on the finite volume. The accurate unsteady motion is achieved by using a 2nd order accurate, 3-point Euler implicit scheme. The quasi-conservative zonal scheme is used for calculating the flow variables on the zonal interface between the rotor and stator. The axial gap between stator and rotor has been configured in two variations, 15% and 65% of average chord length. The analysis program is validated using experimental results and the effect of axial gap is examined. The numerical analysis results are presented by time averaged pressure coefficient and pressure magnitude coefficient and compared with experimental results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.7
• /
• pp.873-881
• /
• 1997

A numerical simulation has been carried out for three-dimensional turbulent flows around an Ahmed body. The Reynolds-averaged Navier-Stokes equation is solved with the SIMPLE method in general curvilinear coordinates system. Several k-.epsilon. turbulence models with two convective difference schemes are evaluated for the performance such as drag coefficient, velocity and pressure fields. The drag coefficient, the velocity and pressure fields are found to be changed considerably with the adopted k-.epsilon. turbulence models as well as the finite difference schemes. The results of simulation prove that the RNG k-.epsilon. model with the QUICK scheme predicts fairly well the tendency of velocity and pressure fields and gives more reliable drag coefficient. It is also demonstrated that the large difference between simulations and experiment in the drag coefficient is due to relatively high predicted values of pressure drag from vertical rear end base.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.2 no.1
• /
• pp.38-50
• /
• 1994

To improve the efficiency of internal combustion engines, it is necessary to understand mixed air-fuel in-cylinder flow processes accurately at intake and compression strokes. There is experimental and numerical methods to analyse in-cylinder flow process. In numerical method, standard $k-{\varepsilon}$ model with wall function was mostly adopted in in-cylinder flow process. But this type model was not efficiently predicted in the near wall region. Therefore in the present study, low Reynolds number $k-{\varepsilon}$ model was adopted near the cylinder wall and standard $k-{\varepsilon}$ model in other region. Also QUICK scheme was used for convective difference scheme. This study takes axisymmetric reciprocating model engine motored at 200rpm with a centrally located valve, incorporated 60 degree seat angie, and flat piston surface excluding inlet port. Because in-cylinder flow processes are undergoing unsteady and compressible, averaged cylinder pressure and inlet velocity at arbitrary crank angle are determined from thermodynamic analytic method and incylinder states at that crank angle are iteratively determined from the numerical analytic method.

• Nuclear Engineering and Technology
• /
• v.41 no.8
• /
• pp.1045-1052
• /
• 2009

A calculation model is developed to predict the transient free surface flow on the containment floor following a loss-of-coolant accident (LOCA) of pressurized water reactors (PWR) for the use of debris transport evaluation. The model solves the two-dimensional Shallow Water Equation (SWE) using a finite volume method (FVM) with unstructured triangular meshes. The numerical scheme is based on a fully explicit predictor-corrector method to achieve a fast-running capability and numerical accuracy. The Harten-Lax-van Leer (HLL) scheme is used to reserve a shock-capturing capability in determining the convective flux term at the cell interface where the dry-to-wet changing proceeds. An experiment simulating a sudden break of a water reservoir with L-shape open channel is calculated for validation of the present model. It is shown that the present model agrees well with the experiment data, thus it can be justified for the free surface flow with accuracy. From the calculation of flow field over the simplified containment floor of APR1400, the important phenomena of free surface flow including propagations and interactions of waves generated by local water level distribution and reflection with a solid wall are found and the transient flow rates entering the Holdup Volume Tank (HVT) are obtained within a practical computational resource.

• Journal of Mechanical Science and Technology
• /
• v.16 no.5
• /
• pp.683-695
• /
• 2002

The objective of the present study is to analyze the fluid flow with moving boundary using a finite element method. The algorithm uses a fractional step approach that can be used to solve low-speed flow with large density changes due to intense temperature gradients. The explicit Lax-Wendroff scheme is applied to nonlinear convective terms in the momentum equations to prevent checkerboard pressure oscillations. The ALE (Arbitrary Lagrangian Eulerian) method is adopted for moving grids. The numerical algorithm in the present study is validated for two-dimensional unsteady flow in a driven cavity and a natural convection problem. To extend the present numerical method to engine simulations, a piston-driven intake flow with moving boundary is also simulated. The density, temperature and axial velocity profiles are calculated for the three-dimensional unsteady piston-driven intake flow with density changes due to high inlet fluid temperatures using the present algorithm. The calculated results are in good agreement with other numerical and experimental ones.

• Journal of Advanced Marine Engineering and Technology
• /
• v.21 no.4
• /
• pp.430-436
• /
• 1997

Time-dependent Flow characteristics of a two-dimensional lid-driven square cavity flow of six high Reynolds numbers, $10^4 2{\times}lO^4., 3{\times}l0^4, 5{\times}lO^4, 7.5{\times}lO^4$ and $10^5$ were investigated. A convection conservative difference scheme based upon SOLA to maintain the nearly 2nd-order spatial accuracy was adopted on irregular grid formation. Irregular grid number is $80{\times}80$ and its minimum size is about 1/400 of the cavity height(H) and its maximum is about 1/53 H. The result shows that at Re= $10^4$, periodic migration of small eddies appearing in corner separation region and its temporal sinusoidal fluctuation are represented. And another critical Reynolds number which shows the commencement of flow randomness emerging from the periodic fluctuation is assumed to be around Re= $1.5{\times}10$. At five higher Reynolds numbers ($2{\times}lO^4., 3{\times}l0^4, 5{\times}lO^4, 7.5{\times}lO^4$ and $10^5$), an organizing structure of four consecutive vortices similar to a Moffat vortex at two lower corners is revealed from time-mean flow patterns.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.1
• /
• pp.1-14
• /
• 1998

A two-dimensional numerical model for tidal currents based on the depth averaged equation is developed. The mode1 employs a rectangular grid system for its simplicity in the application of complicate coastal shore lines. To raise computing efficiency, implicit approximate factorization scheme is implemented in solving governing equations. An upwind-differencing is used to discretize convective terms, which provides a numerical dissipation automatically and suppresses any oscillations caused by nonlinear instabilities. Some numerical tests are made against the analytic solutions of a linearized shallow water equation to validate the developed numerical scheme, and comparisons of the model prediction with the analytic solution are satisfactory. As a real application, the tidal currents are computed on the Inchon area where the tidal currents are important for the design of new canal which is under construction.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.296-296
• /
• 2015

The aim of this paper is to numerically explore the feasibility of designing a Mini-Hydro turbine. The interest for this kind of horizontal axis turbine relies on its versatility. For instance, in the field of renewable energy, this kind of turbine may be considered for different applications, such as: tidal power, run-of-the-river hydroelectricity, wave energy conversion. It is fundamental to improve the turbine performance and to decrease the equipment costs for achievement of "environmental friendly" solutions and maximization of the "cost-advantage". In the present work, the commercial CFD code ANSYS is used to perform 3D simulations, solving the incompressible Unsteady Reynolds-Averaged Navier-Stokes (U-RANS) equations discretized by means of a finite volume approach. The implicit segregated version of the solver is employed. The pressure-velocity coupling is achieved by means of the SIMPLE algorithm. The convective terms are discretized using a second order accurate upwind scheme, and pressure and viscous terms are discretized by a second-order-accurate centered scheme. A second order implicit time formulation is also used. Turbulence closure is provided by the realizable k - turbulence model. In this study, a mini hydro turbine (3kW) has been considered for utilization of horizontal axis impeller. The turbine performance and flow behavior have been evaluated by means of numerical simulations. Moreover, the performance of the impeller varied in the pressure distribution, torque, rotational speed and power generated by the different number of blades and angles. The model has been validated, comparing numerical results with available experimental data.