• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.294-298
• /
• 2003

The complex, unsteady, self-sustained pressure oscillations induced by supersonic flow past a rectangular cavity is investigated using numerical simulations. The present numerical study is performed using a parallel, multiblock solver for the two-dimensional, compressible Navier­Stokes equations. Open cavities with length-to-depth (L / D) ratio in the range 0.5 - 3.3 are considered. This paper sheds light on the cavity physics, cavity oscillatory mechanism, and the organisation of vortical structures inside the cavity. The vortex shedding phenomenon, the shear layer impingement event at the aft wall and the movement of the acoustic/compression wave within the cavity are well predicted. The vortical structures· and the source of the acoustic disturbances are found to be located near the aft wall of the cavity. With the increase in the cavity length, strong recompression of the flow near the aft wall leading to a sudden jump in the cavity form drag is observed. The estimated cavity tones are in good agreement with the available semi­empirical relation. Multiple peaks are noticed in deep and long cavities. For the present free­stream Mach number 1.71, it is observed that around L/D=2.0, the cavity oscillatory mechanism changes from the transverse to longitudinal oscillatory mode. The effects of this transition on various fluid dynamics and acoustic properties are also discussed.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2015년도 학술발표회
• /
• 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.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 2004년도 유체기계 연구개발 발표회 논문집
• /
• pp.396-401
• /
• 2004

In this paper, a 1MW HAWT(FIL-1000) rotor blade has been designed by BEMT(Blade Element Momentum Theory) with Prandtl's tip loss. Also, a 3-D flow and performance analysis on the FIL-1000 rotor blade has been carried out by using the 3-D Navier-Stokes commercial solver (CFX-5.7) to provide more efficient design techniques to the large-scale HAWT engineers. The rated power and itsapproaching wind velocity at design point (TSR=7.5) are 1MW and 9.99m/s respectively. The rotor diameter is 54.5m and the rotating speed is 26.28rpm. Airfoils such as FFA W-301, DU91-W-250, DU93-W-210, NACA 63418, NACA 63415 consist of the rotor blade from hub to tip. Recent CFX version, 5.7 was adopted to simulate 3-D flow field and to analyze the performance characteristics of the rotor blade. Entire mesh node number is about 730,000 and it is generated by ICEM-CFD to achieve better mesh quality The predicted maximum power occurringat the design tip speed ratio is 931.45kW. Approaching to the root, the inflow angle becomes large, which causesthe blade to be stalled in the region. Therefore, k-$\omega$ SST turbulence model was used to predict the quantitative flow information more accurately. Application of commercial CFD code to optimum blade design and performance analysis was proved to be more effective environment to HAWT blade designers.

• 대한조선학회논문집
• /
• 제51권5호
• /
• pp.409-418
• /
• 2014

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.

• 한국자동차공학회논문집
• /
• 제6권2호
• /
• pp.12-20
• /
• 1998

Exter ballistics of a typical high-speed projectile is studied through a flow-visualization experiment and an unstructured grid Navier-Srokes computation. Experiment produced a schlieren photograph that adequately shows the characteristic features of this complex flow, namely two kinds of oblique cone shocks and turbulent wake developing into the downstream. A hybrid scheme of finite volume-element method is used to simulate the compressible Reynolds-Averaged Navier-Stok- es solution on unstructured grids. Osher's approximate Riemann solver is used to discretize the cinvection term. Higher-order spatial accuracy is obtained by MUSCL extension and van Albada ty- pe flux limiter is used to stabilize the numerical oscillation near the solution discontinuity. Accurate Gakerkin method is used to discretize the viscous term. Explict fourth-order Runge-Kutta method is used for the time-stepping, which simplifies the application of MUSCL extension. A two-layer k-$\varepsilon$ turbulence model is used to simulate the turbulent wakes accurately. Axisymmetric folw and two-dimensional flow with an angle of attack have been computed. Grid-dependency is also checked by carrying out the computation with doubled meshes. 2-D calculation shows that effect of angle of attack on the flow field is negligible. Axi-symmetric results of the computation agrees well with the flow visualization. Primary oblique shock is represented within 2-3 meshes in numerical results, and the varicose mode of the vortex shedding is clearly captured in the turbulent wake region.

• 한국해양공학회지
• /
• 제24권5호
• /
• pp.16-21
• /
• 2010

During the loading/offloading operation of a liquefied natural gas carrier (LNGC) that is moored in a side-by-side configuration with an offshore plant, sloshing that occurs due to the partially filled LNG tank and the interactive effect between the two floating bodies are important factors that affect safety and operability. Therefore, a time-domain software program, called CHARM3D, was developed to consider the interactions between sloshing and the motion of a floating body, as well as the interactions between multiple bodies using the potential-viscous hybrid method. For the simulation of a floating body in the time domain, hydrodynamic coefficients and wave forces were calculated in the frequency domain using the 3D radiation/diffraction panel program based on potential theory. The calculated values were used for the simulation of a floating body in the time domain by convolution integrals. The liquid sloshing in the inner tanks is solved by the 3D-FDM Navier-Stokes solver that includes the consideration of free-surface non-linearity through the SURF scheme. The computed sloshing forces and moments were fed into the time integration of the ship's motion, and the updated motion was, in turn, used as the excitation force for liquid sloshing, which is repeated for the ensuing time steps. For comparison, a sloshing motion coupled analysis program based on linear potential theory in the frequency domain was developed. The computer programs that were developed were applied to the side-by-side offloading operation between the offshore plant and the LNGC. The frequency-domain results reproduced the coupling effects qualitatively, but, in general, the peaks were over-predicted compared to experimental and time-domain results. The interactive effects between the sloshing liquid and the motion of the vessel can be intensified further in the case of multiple floating bodies.

• 대한토목학회논문집
• /
• 제29권3B호
• /
• pp.289-301
• /
• 2009

본 연구에서는 Navier-Stokes방정식과 자유수면추적을 위해 도입한 VOF함수의 이류방정식에 기초하고 있는 CADMAS-SURF(CDIT, 2001)에 Fenton(1972)의 3차근사고립파이론에 의한 조파이론을 새롭게 도입한 수치파동수로로부터 직립호안상에서 지진해일(고립파)의 처오름 및 작용파력을 해석한다. 기존의 연구결과와 대비하여 본 수치해석결과의 적용성을 논의하며, 최소자승법에 기초하여 직립호안상에서 최대처오름과 최대파력을 합리적으로 추정할 수 있는 회귀식을 제안한다. 또한, 육상구조물에 작용하는 지진해일(tsunami wave)의 파력산정에 CADMAS-SURF(CDIT, 2001)의 적용을 시도한다. 육상에서 지진해일의 수위변동, 유속 및 파력변동의 특성으로부터 육상구조물에 미치는 지진해일의 직접적인 영향을 추론할 수 있을것이다.

• Ocean Systems Engineering
• /
• 제6권1호
• /
• pp.23-60
• /
• 2016

The major objective of this study was to develop further understanding of 3D nearshore hydrodynamics under a variety of wave and tidal forcing conditions. The main tool used was a comprehensive 3D numerical model - combining the flow module of Delft3D with the WAVE solver of XBeach - of nearshore hydro- and morphodynamics that can simulate flow, sediment transport, and morphological evolution. Surf-swash zone hydrodynamics were modeled using the 3D Navier-Stokes equations, combined with various turbulence models (${\kappa}-{\varepsilon}$, ${\kappa}-L$, ATM and H-LES). Sediment transport and resulting foreshore profile changes were approximated using different sediment transport relations that consider both bed- and suspended-load transport of non-cohesive sediments. The numerical set-up was tested against field data, with good agreement found. Different numerical experiments under a range of bed characteristics and incident wave and tidal conditions were run to test the model's capability to reproduce 3D flow, wave propagation, sediment transport and morphodynamics in the nearshore at the field scale. The results were interpreted according to existing understanding of surf and swash zone processes. Our numerical experiments confirm that the angle between the crest line of the approaching wave and the shoreline defines the direction and strength of the longshore current, while the longshore current velocity varies across the nearshore zone. The model simulates the undertow, hydraulic cell and rip-current patterns generated by radiation stresses and longshore variability in wave heights. Numerical results show that a non-uniform seabed is crucial for generation of rip currents in the nearshore (when bed slope is uniform, rips are not generated). Increasing the wave height increases the peaks of eddy viscosity and TKE (turbulent kinetic energy), while increasing the tidal amplitude reduces these peaks. Wave and tide interaction has most striking effects on the foreshore profile with the formation of the intertidal bar. High values of eddy viscosity, TKE and wave set-up are spread offshore for coarser grain sizes. Beach profile steepness modifies the nearshore circulation pattern, significantly enhancing the vertical component of the flow. The local recirculation within the longshore current in the inshore region causes a transient offshore shift and strengthening of the longshore current. Overall, the analysis shows that, with reasonable hypotheses, it is possible to simulate the nearshore hydrodynamics subjected to oceanic forcing, consistent with existing understanding of this area. Part II of this work presents 3D nearshore morphodynamics induced by the tides and waves.

• Ocean Systems Engineering
• /
• 제7권4호
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.