• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.12
• /
• pp.4013-4026
• /
• 1996

A numerical simulation has been carried out for the jet impinging on a flat plate and a semi-circular concave surface. In this computation finite volume method was employed to solve the full Navier-Stokes equation based on a non-orthogonal coordinate with non staggered variable arrangement. The standard k-.epsilon. turbulent model and low Reynolds number k-.epsilon. model(Launder-Sharmar model) with Yap's correction were adapted. The accuracy of the numerical calculations were compared with various experimental data reported in the literature and showed good predictions of centerline velocity decay, wall pressure distribution and skin friction. For the jet impingement on a semi-circular concave surface, potential core length was calculated for two different nozzle(round edged nozzle and rectangular edged nozzle) to consider effects of the nozzle shape. The result showed that round edged nozzle had longer potential core length than rectangular edged nozzle for the same condition. Heat transfer rate along the concave surface with constant heat flux was calculated for various nozzle exit to surface distance(H/B) in the condition of same jet velocity. The maximum local Nusselt number at the stagnation point occurred at H/B = 8 where the centerline turbulent intensity had maximum value. The predicted Nusselt number showed good agreement with the experimental data at the stagnation point. However heat transfer predictions along the downstream were underestimated. This results suggest that the improved turbulence modeling is required.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.28 no.1
• /
• pp.13-26
• /
• 2016

Sea water level variations are generally influenced by a variety of factors such as tides, meteorological forces, water temperature, salinity, wave, and topography, etc. Among non-tidal conditions, atmospheric pressure is one of the major factors causing water level changes. In the East Sea, due to small tidal range which is opposite to large tidal range of the Yellow Sea, it is difficult to predict water level changes using a numerical model, which consider tidal forcing only. This study focuses on the effects of atmospheric pressure variations on sea level predictions along the eastern coast of Korea. Telemac-2D model is simulated with the Inverted Barometer Effect(IBE), and then its results are analyzed. In comparison between observed data and predictions, the correlation of prediction with IBE and tide is better than that of tide-only case. Therefore, IBE is strongly suggested to be considered for the numerical simulations of sea level changes in the East Sea.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.8 no.1
• /
• pp.110-116
• /
• 2016

In this study an attempt has been made to study the hydrodynamic performance of pumpjet propulsor. Numerical investigation based on the Reynolds Averaged NaviereStokes (RANS) computational fluid dynamics (CFD) method has been carried out. The structured grid and SST ${\kappa}-{\omega}$ turbulence model have been applied. The numerical simulations of open water performance of marine propeller E779A are carried out with different advance ratios to verify the numerical simulation method. Results show that the thrust and the torque are in good agreements with experimental data. The grid independent inspection is applied to verify accuracy of numerical simulation grid. The numerical predictions of hydrodynamic performance of pumpjet propulsor are carried out with different advance ratios. Results indicate that the rotor provides the main thrust of propulsor and the balance performance of propulsor is generally satisfactory. Additionally, the curve of propulsor efficiency is in good agreement with experimental data. Furthermore, the pressure distributions around rotor and stator blades are reasonable. Beyond that, the existence of tip clearance accounts for the appearance of tip vortex that leads to a further loss in efficiency and a probability of cavitation phenomenon.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.2 no.1
• /
• pp.39-44
• /
• 2010

Recently underwater systems moving at high speed such as a super-cavitating torpedo have been studied for their practical advantage of the dramatic drag reduction. In this study we are focusing our attention on super-cavitating flows around axisymmetric cavitators. A numerical method based on inviscid flow is developed and the results for several shapes of the cavitator are presented. First using a potential based boundary element method, we find the shape of the cavtiator yielding a sufficiently large enough cavity to surround the body. Second, numerical predictions of supercavity are validated by comparing, with experimental observations carried out in a high speed cavitation tunnel at Chungnam National University (CNU CT).

• Architectural research
• /
• v.12 no.2
• /
• pp.93-98
• /
• 2010

This paper describes a numerical method for estimating the elastic moduli of cement paste. The cement paste is modeled as a unit cell which consists of three components: the unhydrated cement grain, the gel, and the capillary pore. In the unit cell, the volume fractions of the constituents are quantified using a single kinetic function calculating the degree of hydration. The elastic moduli of cement paste are calculated from the total displacements of constituents when a uniform pressure is applied to the gel contact area. The cement paste is assumed to be a homogenous isotropic matrix. Numerical simulations were conducted through the finite element analysis of the three-dimensional periodic unit cell. The model predictions are compared with experimental results. The predicted trends are in good agreement with experimental observations. This approach and some of the results might also be relevant for other technical applications.

• Wind and Structures
• /
• v.9 no.6
• /
• pp.473-492
• /
• 2006

This paper presents numerical results of pedestrian-level wind environment around the base of a row of tall buildings by CFD. Four configurations of building arrangement are computed including a single square tall building. Computed results of pedestrian-level wind flow patterns and wind speeds are compared to previous wind tunnel measurement data to enable an assessment of CFD predictions. The CFD model uses the finite-volume method with RNG $k-{\varepsilon}$ model for turbulence closure. It is found that the numerical results can reproduce key features of pedestrian-level wind environment such as corner streams around corners of upwind building, sheltered zones behind buildings and channeled high-speed flow through a building gap. However, there are some differences between CFD results and wind tunnel data in the wind speed distribution and locations of highest wind speeds inside the corner streams. In locations of high ground-level wind speeds, CFD values match wind tunnel data within ${\pm}10%$.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.298-302
• /
• 2009

Numerical analysis of two dimensional incompressible laminar free jet flow was carried out by using finite difference SMAC scheme. Flow characteristics of free jet flow such as jet width, similarity of jet velocity and hypothetical origin were investigated for different Reynolds numbers of Re=30 and 100. The reliability of predictions were confirmed by comparison with exact solution. Non-dimensional velocity distribution showed similarity of jet flow velocity after initial region. In the region of laminar flow, the location of hypothetical origin becomes more distant with Reynolds number.

• Journal of Mechanical Science and Technology
• /
• v.16 no.10
• /
• pp.1327-1335
• /
• 2002

The shock wave process represents an abrupt change in fluid properties, in which finite variations in pressure, temperature, and density occur over the shock thickness which is comparable to the mean free path of the gas molecules involved. This shock wave fluid phenomenon is simulated by using the finite difference lattice Boltzmann method (FDLBM). In this paper, a new model is proposed using the lattice BGK compressible fluid model in FDLBM for the purpose of speeding up the calculation as well as stabilizing the numerical scheme. The numerical results of the proposed model show good agreement with the theoretical predictions.

• 한국연소학회:학술대회논문집
• /
• 2001.11a
• /
• pp.35-41
• /
• 2001

A second-order conditional moment closure(CMC) model is applied to the prediction of local extinction in a turbulent hydrocarbon diffusion flame and compared with direct numerical simulation(DNS) results for the flame. Combustion of a hydrocarbon fuel is described by a simple two-step mechanism. A second-order correction for conditional mean reaction rate terms is made by the assumed pdf method. The results show that the second-order closure is necessary for accurate prediction of intermediate species, while first-order CMC gives good predictions for fuel, oxidant, product and temperature. Conditional variances and covariances are well predicted during an extinction process while they are overpredicted during a reignition process.

• Structural Engineering and Mechanics
• /
• v.5 no.4
• /
• pp.369-383
• /
• 1997

The classical theory of thin-walled members is unable to reflect the shear lag phenomenon since it is based on the assumption of no shearing strains in the middle surface of the walls. In this paper, an energy equation for the lateral buckling of thin-walled members has been derived which includes the effects of torsion, warping and, especially, the shearing strains which reflect the shear lag phenomenon. A numerical analysis for the lateral buckling of thin-walled members with openings by using Galerkin's method of weighted residuals has been presented. The proposed numerical values and the predictions by experiment for the lateral buckling loads are to agree closely in the paper. The results from these comparisons show that the proposed method here is capable of predicting the lateral buckling of thin-walled members with openings. The fast convergence of the results indicates the numerical stability of the method. By the study, a very complex practical eigenvalue problem is transformed into a very simple one of solving only a linear equation with one variable.