• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.6
• /
• pp.92-100
• /
• 2002

In this study, transonic/supersonic nonlinear flutter analysis system of a complete aircraft including forced harmonic motion pf control surfaces has been effectively developed using the modified transonic small disturbance (TSD) equation. To consider the nonlinear effects, the coupled time marching method (CTM) combining computational structural dynamics (CFD) has been directly applied for aeroelastic computations. The grid system for a complex full aircraft configuration is effectively generated by the developed inhouse code. Intransonic and supersonic flight regimes, the characteristics of static and dynamic aeroelastic effect has been investigated for a complete aircraft model. Also, nonlinear flutter flight simulations for the forced harmonic motion of control surfaces are practically presented in detail.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.28 no.6
• /
• pp.332-349
• /
• 2016

OLAFOAM is the powerful CFD code and is an expanded version of $OpenFOAM^{(R)}$, for wave mechanics simulation. The $OpenFOAM^{(R)}$ does provide many solvers to correspond to each object of the numerical calculation in a variety of fields. OLAFOAM's governing equation bases on VARANS (Volume-Averaged Reynolds-Averaged Navier-Stokes) equation, and the finite volume method is applied to numerical techniques. The program is coded in C++ and run on the Linux operating system. First of all, in this study, OLAFOAM was validated for 1) wave transformation inside porous structure under bore and regular wave conditions, 2) wave transformation by submerged breakwater under regular wave condition, and 3) regular wave transformation and resultant vertical velocity distribution under current by comparison with existing laboratory measurements. Hereafter, this study, which is almost no examination carried out until now, analyzed closely variation characteristics of water surface level, wave height, frequency spectrum, breaking waves, averaged velocity and turbulent kinetic energy around porous submerged breakwater in the wave and current coexisting field for the case of permeable or impermeable rear beach. It was revealed that the wave height fluctuation according to current direction(following or opposing) was closely related to the turbulent kinetic energy, and others.

• Nuclear Engineering and Technology
• /
• v.54 no.3
• /
• pp.842-848
• /
• 2022

Parametric studies of heat transfer and fluid flow are very important research of interest because the design and operation of fluid flow and heat transfer systems are guided by these parametric studies. The safety of the system operation and system optimization can be determined by decreasing or increasing particular fluid flow and heat transfer parameter while keeping other parameters constant. The parameters that can be varied in order to determine safe and optimized system include system pressure, mass flow rate, heat flux and coolant inlet temperature among other parameters. The fluid flow and heat transfer systems can also be enhanced by the presence of or without the presence of particular effects including gravity effect among others. The advanced Generation IV reactors to be deployed for large electricity production, have proven to be more thermally efficient (approximately 45% thermal efficiency) than the current light water reactors with a thermal efficiency of approximately 33 ℃. SCWR is one of the Generation IV reactors intended for electricity generation. High Performance Light Water Reactor (HPLWR) is a SCWR type which is under consideration in this study. One-eighth of a proposed fuel assembly design for HPLWR consisting of 7 fuel/rod bundles with 9 coolant sub-channels was the geometry considered in this study to examine the effects of system pressure and mass flow rate on wall and fluid temperatures. Gravity effect on wall and fluid temperatures were also examined on this one-eighth fuel assembly geometry. Computational Fluid Dynamics (CFD) code, STAR-CCM+, was used to obtain the results of the numerical simulations. Based on the parametric analysis carried out, sub-channel 4 performed better in terms of heat transfer because temperatures predicted in sub-channel 9 (corner subchannel) were higher than the ones obtained in sub-channel 4 (central sub-channel). The influence of system mass flow rate, pressure and gravity seem similar in both sub-channels 4 and 9 with temperature distributions higher in sub-channel 9 than in sub-channel 4. In most of the cases considered, temperature distributions (for both fluid and wall) obtained at 25 MPa are higher than those obtained at 23 MPa, temperature distributions obtained at 601.2 kg/h are higher than those obtained at 561.2 kg/h, and temperature distributions obtained without gravity effect are higher than those obtained with gravity effect. The results show that effects of system pressure, mass flowrate and gravity on fluid flow and heat transfer are significant and therefore parametric studies need to be performed to determine safe and optimum operating conditions of fluid flow and heat transfer systems.

• Journal of Ship and Ocean Technology
• /
• v.10 no.4
• /
• pp.12-23
• /
• 2006

In this paper, a numerical study is carried out to investigate the turbulent flow around a twin-skeg container ship model with rudders including propeller effects. A commercial CFD code, FLUENT is used with body forces distributed on the propeller disk to simulate the ship stem and wake flows with the propeller in operation. A multi-block, matching, structured grid system has been generated for the container ship hull with twin-skegs in consideration of rudders and body-force propeller disks. The RANS equations for incompressible fluid flows are solved numerically by using a finite volume method. For the turbulence closure, a Reynolds stress model is used in conjunction with a wall function. Computations are carried out for the bare hull as well as the hull with appendages of a twin-skeg container ship model. For the bare hull, the computational results are compared with experimental data and show generally a good agreement. For the hull with appendages, the changes of the stem flow by the rudders and the propellers have been analyzed based on the computed result since there is no experimental data available for comparison. It is found the flow incoming to the rudders has an angle of attack due to the influence of the skegs and thereby the hull surface pressure and the limiting streamlines are changed slightly by the rudders. The axial velocity of the propeller disk is found to be accelerated overall by about 35% due to the propeller operation with the rudders. The area and the magnitude of low pressure on the hull surface enlarge with the flow acceleration caused by the propeller. The propellers are found to have an effect on up to the position where the skeg begins. The propeller slipstream is disturbed strongly by the rudders and the flow is accelerated further and the transverse velocity vectors are weakened due to the flow rectifying effect of the rudder.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.6
• /
• pp.461-468
• /
• 2012

In this study, a Cartesian-grid method based on finite volume approach is applied to simulate the ship motions in large amplitude waves. Fractional step method is applied for pressure-velocity coupling and TVD limiter is used to interpolate the cell face value for the discretization of convective term. Water, air, and solid phases are identified by using the concept of volume-fraction function for each phase. In order to capture the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with weighed line interface calculation (WLIC) method which considers multidimensional information. The volume fraction of solid body embedded in the Cartesian grid system is calculated using a level-set based algorithm, and the body boundary condition is imposed by a volume weighted formula. Numerical simulations for the two-dimensional barge type model and Wigley hull in linear waves have been carried out to validate the newly developed code. To demonstrate the applicability for highly nonlinear wave-body interactions such as green water on the deck, numerical analysis on the large-amplitude motion of S175 containership is conducted and all computational results are compared with experimental data.

• Wind and Structures
• /
• v.23 no.6
• /
• pp.505-521
• /
• 2016

The objective of this study is to investigate numerically the effect of building roof shaps on wind flow and pollutant dispersion in a street canyon with one row of trees of pore volume, $P_{vol}=96%$. A three-dimensional computational fluid dynamics (CFD) model is used to evaluate air flow and pollutant dispersion within an urban street canyon using Reynolds-averaged Navier-Stokes (RANS) equations and the Explicit Algebraic Reynolds Stress Models (EARSM) based on k-${\varepsilon}$ turbulence model to close the equation system. The numerical model is performed with ANSYS-CFX code. Vehicle emissions were simulated as double line sources along the street. The numerical model was validated by the wind tunnel experiment results. Having established this, the wind flow and pollutant dispersion in urban street canyons (with six roof shapes buildings) are simulated. The numerical simulation results agree reasonably with the wind tunnel data. The results obtained in this work, indicate that the flow in 3D domain is more complicated; this complexity is increased with the presence of trees and variability of the roof shapes. The results also indicated that the largest pollutant concentration level for two walls (leeward and windward wall) is observed with the upwind wedge-shaped roof. But the smallest pollutant concentration level is observed with the dome roof-shaped.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.4
• /
• pp.363-374
• /
• 2009

The 2D flow around 13 similar stay-vane profiles with different trailing edge geometries is investigated to determinate the main characteristics of the excitation forces for each one of them and their respective dynamic behaviors when modeled as a free-oscillating system. The main goal is avoid problems with cracks of hydraulic turbines components. A stay vane profile with a history of cracks was selected as the basis for this work. The commercial finite-volume code $FLUENT^{(R)}$ was employed in the simulations of the stationary profiles and, then, modified to take into account the transversal motion of elastically mounted profiles with equivalent structural stiffness and damping. The k-$\omega$ SST turbulence model is employed in all simulations and a deforming mesh technique used for models with profile motion. The static-model simulations were carried out for each one of the 13 geometries using a constant far field flow velocity value in order to determine the lift force oscillating frequency and amplitude as a function of the geometry. The free-oscillating stay-vane simulations were run with a low mass-damping parameter ($m^*{\xi}=0.0072$) and a single mean flow velocity value (5m/s). The structural bending stiffness of the stay-vane is defined by the Reduced Velocity parameter (Vr). The dynamic analyses were divided into two sets. The first set of simulations was carried out only for one profile with $2{\leq}Vr{\leq}12$. The second set of simulations focused on determining the behavior of each one of the 13 profiles in resonance.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.9
• /
• pp.34-40
• /
• 2005

In present study, a two dimensional unsteady Navier-Stokes solver has been developed using the Diagonalized ADI (DADI) method and implicit dual time stepping method. The jacobian matrices in steady state Navier-Stokes equations are introduced from inviscid flux terms. The implicit treatment of artificial dissipation terms results in a block penta-diagonal matrix system and it becomes a scalar penta-diagonal matrix by diagonalization. In steady state equations about fictitious time, a new residual including a real time derivative term is introduced. From a converged solution about fictitious time, a real time unsteady solution can be obtained, which is called 'implicit dual time stepping method'. For code validation, an oscillating flat plate, a regular Karman vortices past a circular cylinder and shock buffeting around a bicircular airfoil problems are numerically solved. And they are compared with a theoretical solution, experiments and other researcher's computations.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.47 no.12
• /
• pp.831-840
• /
• 2019

Flush Air Data Sensing system (FADS) estimates air data states using pressure data measured at the surface of flight vehicles. The FADS system does not require intrusive probes, so it is suitable for high performance aircrafts, stealth vehicles, and hypersonic flight vehicles. In this study, calibration procedures and solution algorithms of the FADS for a sphere-cone shape vehicle are presented for the prediction of air data from subsonic to supersonic flights. Five flush pressure ports are arranged on the surface of nose section in order to measure surface pressure data. The algorithm selects the concept of separation for the prediction of flow angles and the prediction of pressure related variables, and it uses the pressure model which combines the potential flow solution for a subsonic flow with the modified Newtonian flow theory for a hypersonic flow. The CFD code which solves Euler equations is developed and used for the construction of calibration pressure data in the Mach number range of 0.5~3.0. Tests are conducted with various flight conditions for flight Mach numbers in the range of 0.6~3.0 and flow angles in the range of -10°~+10°. Air data such as angle of attack, angle of sideslip, Mach number, and freestream static pressure are predicted and their accuracies are analyzed by comparing predicted data with reference data.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.5
• /
• pp.324-333
• /
• 2013

Recently, rudder erosion due to cavitation has been frequently reported on a semi-spade rudder of a high-speed large ship. This problem raises economic and safety issues when operating ships. The semi-spade rudders have a gap between the horn/pintle and the movable wing part. Due to this gap, a discontinuous surface, cavitation phenomenon arises and results in unresolved problems such as rudder erosion. In this study, we made a rudder model for 2-D experiments using the NACA0020 and also manufactured gap flow blocking devices to insert to the gap of the model. In order to study the gap flow characteristics at various rudder deflection angles($5^{\circ}$, $10^{\circ}$, $35^{\circ}$) and the effect of the gap flow blocking devices, we carried out the velocity measurements using PIV(Particle Image Velocimetry) techniques and cavitation observation using high speed camera in Seoul National University cavitation tunnel. To observe the gap cavitation on a semi-spade rudder, we slowly lowered the inside pressure of the cavitation tunnel until cavitation occurred near the gap and then captured it using high-speed camera with the frame rate of 4300 fps(frame per second). During this procedure, cavitation numbers and the generated location were recorded, and these experimental data were compared with CFD results calculated by commercial code, Fluent. When we use gap flow blocking device to block the gap, it showed a different flow character compared with previous observation without the device. With the device blocking the gap, the flow velocity increases on the suction side, while it decreases on the pressure side. Therefore, we can conclude that the gap flow blocking device results in a high lift-force effect. And we can also observe that the cavitation inception is delayed.