• Journal of the Earthquake Engineering Society of Korea
• /
• v.7 no.2
• /
• pp.21-27
• /
• 2003

In general, simple fluid added mass method is used for the seismic and vibration analysis of the immersed structure to consider the fluid-structure interaction effect. Actually, the structural response of the immersed structure can be affected by both the fluid added mass and damping caused by the fluid viscosity. These variables appeared as a consistent matrix form with the coupling terms. In this paper, finite element formula for the inviscid fluid case and viscous fluid case are derived from the linearized Navier Stoke's equations. Using the finite element program developed in this paper, the analyses of fluid added mass and damping for the hexagon core structure of the liquid metal reactor are carried out to investigate the effect of fluid viscosity with variation of the fluid gap and Reynolds number. From the analysis results, it is verified that the viscosity significantly affects the fluid added mass and damping as the fluid gap size decrease. From the analysis results of eccentricity effect on the fluid added mass and damping of the concentric cylinders, the fluid added mass increase as the eccentricity increases, however the fluid damping increases only when the eccentricity is very severe.

• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.23-36
• /
• 2015

The $2^{nd}$ generation ice accretion analysis program has been developed and validated for various icing conditions. The essential feature of the $2^{nd}$ generation code lies in its capability of handling general 3-D geometry and improved accuracy. The entire velocity fields are obtained based on Navier-Stokes equations in order to take the massively separated flow field into account. Unlike $1^{st}$ generation code, the droplet trajectories are calculated using Eulerian approach, which is adopted to yield appropriate collection efficiency even in the shadow region. For improved thermodynamic analysis on the surfaces, water film model and modified Messinger model are newly included in the present analysis. The ice shape for a given time step is obtained by considering the exact amount of ice accreted on the surface. Each module of the icing analysis code has been seamlessly integrated on the OpenFOAM platform. The developed code was validated against available experimental data for 2D airfoils and 3D DLR-F4. Due to the lack of experimental data, the computed results of DLR-F4 were compared with those obtained from FENSAP-ICE, which is state-of-the-art 3D icing analysis code. It was clearly shown that the present code produces comparable results to those of FENSAP-ICE, in terms of prediction accuracy and the capability of handling general 3-D geometries.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.4
• /
• pp.274-280
• /
• 2017

Generally, global performance analysis in offshore platforms is performed using potential-based numerical tools, which neglect hydrodynamic viscous effects. In comparison with the potential theory, computational fluid dynamics (CFD) methods can take into account the viscous effects by solving the Navier-Stokes equation using the finite-volume method. The open-source field operation and manipulation (OpenFOAM) C++ libraries are employed for a finite volume method (FVM) numerical analysis. In this study, in order to apply CFD to the global performance analysis of a hull-mooring coupled system, we developed a numerical wave basin to analyze the global performance problem of a floating body with a catenary mooring system under regular wave conditions. The mooring system was modeled using a catenary equation and solved in a quasi-static condition, which excluded the dynamics of the mooring lines such as the inertia and drag effects. To demonstrate the capability of the numerical basin, the global performance of a barge with four mooring lines was simulated under regular wave conditions. The simulation results were compared to the analysis results from a commercial mooring analysis program, Orcaflex. The comparison included the motion of the barge, catenary shape, and tension in the mooring lines. The study found good agreement between the results from the developed CFD-based numerical calculation and commercial software.

• Journal of KIISE:Computer Systems and Theory
• /
• v.32 no.1
• /
• pp.29-38
• /
• 2005

The fluid animation procedure consists of physical simulation and visual rendering. In the physical simulation of fluids, the most frequently used practices are the numerical simulation of fluid particles using particle dynamics equations and the continuum analysis of flow via Wavier-Stokes equation. Particle dynamics method is fast in calculation, but the resulting fluid motion is conditionally unrealistic The method using Wavier-Stokes equation, on the contrary, yields lifelike fluid motion when properly conditioned, yet the complexity of calculation restrains this method from being used in real-time applications. Global illumination is generally successful in producing premium-Duality rendered images, but is also excessively slow for real-time applications. In this paper, we propose a rapid fluid animation method incorporating enhanced particle dynamics simulation method and pre-integrated volume rendering technique. The particle dynamics simulation of fluid flow was conducted in real-time using Lennard-Jones model, and the computation efficiency was enhanced such that a small number of particles can represent a significant volume. For real-time rendering, pre-integrated volume rendering method was used so that fewer slices than ever can construct seamless inter-laminar shades. The proposed method could successfully simulate and render the fluid motion in real time at an acceptable speed and visual quality.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.13-18
• /
• 2009

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought about by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 2.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}106$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental data in the low aspect ratio cavity (L/D = ~ 4.5). In the large aspect ratio cavity, however, there are other low dominant frequencies due to the leading edge shear layer with the dominant frequencies of the feedback mechanism. The characteristics of the acoustic wave propagation are analyzed using the Correlation of Pressure Distribution (CPD).

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.7-13
• /
• 2010

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 5.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}10^6$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental datum in the low aspect ratio cavity (L/D = ~4.5). In the high aspect ratio cavity, however, there are other low dominant frequencies of the leading edge shear layer with the dominant frequencies of the feedback mechanism.

• Transactions of the Korean hydrogen and new energy society
• /
• v.27 no.6
• /
• pp.764-772
• /
• 2016

This paper presents on the evaluation of numerical analysis model of a ball valve used for a gas pipeline. The ball valve has important role to control the gas flow of the pipeline as well as safety operation to prevent gas explosion at the emergency. For the validation of numerical simulation, the computational domains are introduced three different types: a hexahedron chamber connected to a pipeline outlet without considering the geometry of pressure tubes, a pipeline only considered the geometry of pressure tubes, and a pipeline connected both of the a hexahedron chamber and pressure tubes. The commercial code, SC/Tetra, is introduced to solve the three-dimensional steady-state Reynolds-averaged Navier-Stokes analysis in the present study. The valve flow coefficient and valve loss coefficient with respect to the valve opening rate of 30%, 50%, and 70% are compared with experimental results. Throughout the numerical analysis for the three analysis domains, pressure computed along the pipeline is affected by computational domains. It is noted pressure obtained by the computational model considering both of the a hexahedron chamber and pressure tubes has a relatively good agreement to the experimental data.

• Journal of the Korean Solar Energy Society
• /
• v.35 no.1
• /
• pp.35-43
• /
• 2015

A numerical simulation on the wake flow of a wind turbine which is a scaled version of a multi-megawatt wind turbine has been performed. Two different inlet conditions of averaged wind speed including one below and one above the rated wind speed were used in the simulation. Steady-state pitch angles of the blade associated with the two averaged wind speeds were imposed for the simulation. The steady state analysis based on the Reynolds averaged Navier-Stokes equations with the method of frame motion were used for the simulation to find the torque of the rotor and the wake field behind the wind turbine. The simulation results were compared with the results obtained from the wind tunnel testing. From comparisons, it was found that the simulation results on the turbine power are pretty close to the experimental values. Also, the wake results were relatively close to the experimental results but there existed some discrepancy in the shape of velocity deficit. The reason for the discrepancy is considered due to the steady state solution with the frame motion method used in the simulation. However, the method is considered useful for solutions with much reduced calculation time and reasonably good accuracy compared to the transient analysis.

• Nuclear Engineering and Technology
• /
• v.56 no.1
• /
• pp.9-18
• /
• 2024

In-vessel retention through external reactor vessel cooling (IVR-ERVC) is a severe accident management (SAM) strategy that has been adopted and used in many nuclear reactors such as AP1000, APR1400, and light water reactor etc. Some reactor accidents have raised concerns about nuclear reactors among residents, leading to a decrease in residents' acceptability and many studies on SAM are being conducted. Experiments on IVR-ERVC are almost impossible due to its specificity, so fluid characteristics are analyzed through BALI experiments with similar condition. In this study, computational fluid dynamics (CFD) via Reynolds-averaged Navier-Stokes (RANS) and large eddy simulation (LES) for BALI experiments were performed. Steady-state CFD analysis was performed on three turbulence models, and SST k-ω model was in good agreement with the experimental measurement temperature within the maximum error range of 1.9%. LES CFD analysis was performed based on the RANS analysis results and it was confirmed that the temperature and wall heat flux for depth was consistent within an error range of 1.0% with BALI experiment. The LES CFD analysis results were compared with those of the Lagrangian-based solver. LES matched the temperature distribution better than SOPHIA, but SOPHIA calculated the position of boundary between stratified layer and convective layer more accurately. On the other hand, Lagrangian-based solver predicted several small eddy behaviors of the convective layer and LES predicted large vortex behavior. The vibration characteristics near the cooling part of the BALI experimental device were confirmed through Fast Fourier Transform (FFT) investigation. It was found that the power spectral density for pressure at least 10 times higher near the side cooling than near the top cooling.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.3
• /
• pp.301-307
• /
• 2011

Turbo blowers are mainly used in refuse collection systems. We discuss blower performance in relation to the angle of the inlet vane installed at the upstream of the blower. The flow characteristics of the components are analyzed by three-dimensional Navier-Stokes analysis and compared to experimental results. A two-stage serially connected turbo blower is introduced to analyze the performance experimentally. Throughout the experimental measurements and the numerical simulation, the distorted inlet velocity generated in the small vane angle reduces the performance of the blower, because of the local leading-edge separation and the resulting non-uniform blade loading. We also perform a detailed flow analysis using the results obtained in the numerical simulation.