• Journal of the Korean Society for Marine Environment & Energy
• /
• v.11 no.1
• /
• pp.35-41
• /
• 2008

Overtopping wave energy convertor is an offshore wave energy convertor for collecting the overtopping waves converting the water pressure head into electric power through the hydro turbines. This paper presents a numerical wave tank based on the commercial CFD code Fluent. The Reynolds Averaged Naiver-Stokes and VOF model is utilized to generate the 2D numerical linear propagating waves, which has been validated by the analytical solutions. Several incident wave conditions and shape parameters are calculated in the optimal designing investigation of the overtopping characteristics and discharge for the overtopping wave energy convertor.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.2
• /
• pp.94-100
• /
• 2007

An axisymmetric supersonic jet screech in the Mach number range from 1.07 to 1.2 is numerically simulated. The axisymmetric mode is the dominant screech mode for an axisymmetric jet. The Reynolds-averaged Navier-Stokes equations in the conjunction with a modified Spalart-Allmaras turbulence model are employed. A high resolution finite volume essentially non-oscillatory(ENO) schemes are used along with nonreflecting characteristic boundary conditions that are crucial to screech tone computations to accurately capture the sound waves, shock-cell structures and large-scale instability waves.

• International Journal of Fluid Machinery and Systems
• /
• v.9 no.3
• /
• pp.229-236
• /
• 2016

Exergetic analysis was introduced in optimization of a rotating equilateral triangular internal cooling channel with staggered square ribs to maximize the net exergy gain. The objective function was defined as the net exergy gain considering the exergy gain by heat transfer and exergy losses by friction and heat transfer process. The flow field and heat transfer in the channel were analysed using three-dimensional Reynolds-averaged Navier-Stokes equations under the uniform temperature condition. Shear stress transport turbulence model has been selected as a turbulence closure through the turbulence model test. Computational results for the area-averaged Nusselt number were validated compared to the experimental data. Three design variables, i.e., the angle of rib, the rib pitch-to-hydraulic diameter ratio and the rib width-to-hydraulic diameter ratio, were selected for the optimization. The optimization was performed at Reynolds number, 20,000. Twenty-two design points were selected by Latin hypercube sampling, and the values of the objective function were evaluated by the RANS analysis at these points. Through optimization, the objective function value was improved by 22.6% compared to that of the reference geometry. Effects of the Reynolds number, rotation number, and buoyancy parameter on the heat transfer performance of the optimum design were also discussed.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2157-2161
• /
• 2003

In this paper, the response surface method using three-dimensional Navier-Stokes analysis to optimize the shape of a forward-curved blades centrifugal fan, is described. For numerical analysis, Reynolds-averaged Navier-Stokes equations with standard k-e turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in forward-curved blades centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. 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. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. It was found that the optimization process provides reliable design of this kind of fans with reasonable computing time

• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2002.09a
• /
• pp.216-218
• /
• 2002

Fluence and LET spectrum for 290,400 MeV/u $\^$12/C and 400 MeV/u $\^$20/Ne beams have been measured by a $\Delta$E-E counter telescope. Total charge-changing cross sections are deduced from measured fluence. The measured cross sections agree with previous measurements, however, they are disagreement with a model calculation. To dose-averaged LETs, the model calculation can reproduce the measured LETs except for peak LETs at Bragg peak region.

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

VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.24 no.7
• /
• pp.569-574
• /
• 2014

In order to investigate low frequency swishing noise of wind turbines, acoustic source model using a singularity in circular motion is introduced to derive analytic solution of Lowson acoustic analogy in time domain. Results in time and frequency domains computed by the solution show apparent modulation of amplitude and frequency. The solution indicates that time histories of acoustic pressure at receiver points varied significantly according to receiver's directional location, even when the retarded time distributions are similar. However, the corresponding time-averaged spectra of sound pressure at the receiver locations where the retarded time distributions are almost same are not significantly different. It can be inferred from these results that the time-averaged sound pressure spectra which cannot take into account the detailed difference in the time-variation of wind turbine noise may not represent the sound quality of wind turbines due to its swishing. Finally, as an introduction of procedure to quantify low frequency swishing noise level, relative variation of overall sound pressure level is obtained using tonal low frequency noise model.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2007.04a
• /
• pp.329-334
• /
• 2007

An axisymmetric supersonic jet screech in the Mach number range from 1.07 to 1.2 is numerically simulated. The axisymmetric mode is the dominant screech mode for an axisymmetric jet. The Reynolds-averaged Navier-Stokes equations in the conjunction with modified Spalart-Allmaras turbulence model are employed. A high resolution finite volume essentially non-oscillatory(ENO) schemes are used along with nonreflecting characteristic boundary conditions that are crucial to screech tone computations to accurately capture the sound waves, shock-cell structures, unsteady shock motions and large-scale instability waves.

• Journal of Environmental Science International
• /
• v.2 no.1
• /
• pp.17-26
• /
• 1993

The layer that is directly influenced by ground surface is called the atmospheric boutsdary layer in comparison with the free atmosphere of higher layer. In the boundary layer, the changes of wind, temperature and coefficient of turbulent diffusion in altitude are large and have great influences an atmospheric diffusion. The purpose of this paper is to express the structure and characteristics of development of mixed layer by using laboratory experiment and numerical simulation. Laboratory experiment using water tank are performed that closely simulate the process of break up of nocturnal surface inversion above heated surface and its phenomena are analyzed by the use of horizontally averaged temperature which is observed. The result obtained from the laboratory experiment is compared with theoretical ones from ; \textsc{k}-\varepsilon numerical model. The results are summarized as follows. 1) The horizontally averaged temperature was found to vary smoothly with height and the mixed layer developed obviously being affected by the convection. 2) The mean height of mixed layer may be predicted as a function of time, knowing the mean initial temperature gradient. The experimental values are associated well with the theoretical values computed for value of the universal constant $C_r$= 0.16, our $C_r$ value is little smaller than the value found by Townsend and Deardoru et al.

• 유체기계공업학회:학술대회논문집
• /
• 2003.12a
• /
• pp.385-390
• /
• 2003

Shape of a multi-blades centrifugal fan is optimized by response surface method based on three-dimensional Navier-Stokes analysis. For numerical analysis, Reynolds-averaged Wavier-Stokes equations with standard $k-{\varepsilon}$ turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in this centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Optimizations with and without constraints are carried out. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. The correlation of efficiency with relative size of inactive zone at the exit of impeller is discussed as well as with average momentum fluxes in the scroll.