• Journal of Ship and Ocean Technology
• /
• v.8 no.2
• /
• pp.61-69
• /
• 2004

Unsteady nonlinear wave motions on the free surface over a plane beach of constant slope are numerically simulated using a finite difference method in rectangular grid system. Two-dimensional Navier-Stokes equations and the continuity equation are used for the computations. Irregular leg lengths and stars are employed near the boundaries of body and free surface to satisfy the boundary conditions. Also, the free surface which consists of markers or segments is determined every time step with the satisfaction of kinematic and dynamic free surface conditions. Moreover, marker-density method is also adopted to allow plunging jets impinging on the free surface. The second-order Stokes wave theory is employed for the generation of waves on the inflow boundary. For the simulation of wave breaking phenomena, the computations are carried out with the plane beach of constant slope in surf zone. The results are compared with other existing experimental results. Agreement between the experimental data and the computation results is good.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.4 s.247
• /
• pp.364-372
• /
• 2006

The objective of this study is to evaluate the efficiency and the prediction accuracy of developed large eddy simulation (LES) program for complex turbulent flows, such as recirculating and swirling flows. To save the computational cost, a Beowulf cluster system consisting 16 processors was constructed. The flows in backward-facing step and dump combustor were examined as representative recirculating and swirling flows. Firstly, a direct numerical simulation (DNS) for laminar backward-facing step flows was previously conducted to validate the overall performance of program. Then LES was carried out for turbulent backward-facing step flows. The results of laminar flow showed a qualitative and quantitative agreement between simulations and experiments. The simulations of the turbulent flow also showed reasonable results. Secondly, LES results for non-swirling and swirling flows in a dump combustor were compared with the results of Reynolds-averaged Navier-Stokes (RANS) using standard $k-{\varepsilon}$ model. The results show that LES has a better performance in predicting the mean axial and azimuthal velocities, comer recirculation zone (CRZ) and center toroidal recirculation zone (CTRZ) than those of RANS. Finally, it was examined the capability of LES for the description of unsteady phenomena.

• Wind and Structures
• /
• v.24 no.3
• /
• pp.223-247
• /
• 2017

By using numerical simulation, vast and detailed information and observation of the physics of flow over a train model can be obtained. However, the accuracy of the numerical results is questionable as it is affected by grid convergence error. This paper describes a systematic method of computational grid refinement for the Unsteady Reynolds Navier-Stokes (URANS) of flow around a generic model of trains using the OpenFOAM software. The sensitivity of the computed flow field on different mesh resolutions is investigated in this paper. This involves solutions on three different grid refinements, namely fine, medium, and coarse grids to investigate the effect of grid dependency. The level of grid independence is evaluated using a form of Richardson extrapolation and Grid Convergence Index (GCI). This is done by comparing the GCI results of various parameters between different levels of mesh resolutions. In this study, monotonic convergence criteria were achieved, indicating that the grid convergence error was progressively reduced. The fine grid resolution's GCI value was less than 1%. The results from a simulation of the finest grid resolution, which includes pressure coefficient, drag coefficient and flow visualization, are presented and compared to previous available data.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.1 no.1
• /
• pp.8-23
• /
• 1997

Steady and unsteady numerical simulations are conducted for the comparison with the experiments performed to investigate the ram accelerator flow field by using an expansion tube facility in Stanford University. Wavier-Stokes equations for chemically reacting flows are analyzed by fully implicit and time accurate numerical methods with Jachimowski's detailed chemistry model for hydrogen-air combustion involving 9 species and 19 reaction steps. Although the steady state numerical simulation shows a good agreement with the experimental schlieren and OH PLIF images for the case of $2H_2$＋$O_2$＋$17N_2$ fails in reproducing the combustion region behind the shock intersection point shown in the case of $2H_2$＋$O_2$＋$12N_2$ mixture. Therefore, an unsteady numerical simulation is conducted for this case and the result shows all the detailed flow stabilization process. From the result of unsteady numerical simulation, the experimental result seems to be an instantaneous state during the flow stabilization process. The combustion behind the shock intersection point is the result of a normal detonation formed by the intersection of strong oblique shocks that exist at early stage of the stabilization process. At final stage, the combustion region behind the shock intersection point disappears and the steady state result is retained. The time required for stabilization of the reacting flow in the model ram accelerator is found to be very long in comparison with the experimental test time.

• Wind and Structures
• /
• v.14 no.5
• /
• pp.435-447
• /
• 2011

A Computational Fluid Dynamics model is presented in this study for the simulation of the complex fluid flows with free surfaces inside the Tuned Liquid Column Dampers in horizontal motion. The characteristics of the fluid model of the TLCD in horizontal motion include the free surface of the multiphase flow and the horizontal moving frame. In this study, the time depend unsteady Standard ${\kappa}-{\varepsilon}$ turbulent model based on Navier-Stokes equations is chosen. The volume of fluid (VOF) method and sliding mesh technique are adopted to track the free surface of water inside the vertical columns of TLCD and treat the moving boundary of the walls of TLCD in horizontal motion. Several model solution parameters comprising different time steps, mesh sizes, convergence criteria and discretization schemes are examined to establish model parametric independency results. The simulation results are compared with the experimental data in the dimensionless amplitude of the water column in four different configured groups of TLCDs with four different orifice areas. The predicted natural frequencies and the head loss coefficient of TLCDs from CFD model are also compared with the experimental data. The predicted numerical results agree well with the available experimental data.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.5 no.4
• /
• pp.502-512
• /
• 2013

Simulations of cavitation flow and hull pressure fluctuation for a marine propeller operating behind a hull using the unsteady Reynolds-Averaged Navier-Stokes equations (RANS) are presented. A full hull body submerged under the free surface is modeled in the computational domain to simulate directly the wake field of the ship at the propeller plane. Simulations are performed in design and ballast draught conditions to study the effect of cavitation number. And two propellers with slightly different geometry are simulated to validate the detectability of the numerical simulation. All simulations are performed using a commercial CFD software FLUENT. Cavitation patterns of the simulations show good agreement with the experimental results carried out in Samsung CAvitation Tunnel (SCAT). The simulation results for the hull pressure fluctuation induced by a propeller are also compared with the experimental results showing good agreement in the tendency and amplitude, especially, for the first blade frequency.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.32-38
• /
• 2006

The objective of this study is to evaluate the prediction accuracy of development large eddy simulation(LES) program for turbulent flow behind a bluff-body. The LES solver was implemented on parallel computer consisting 16 processors. To verify the capability of LES code, the results were compared with those of Reynolds Averaged Navier-Stokes(RANS) using standard ${\kappa}-{\varepsilon}$ model as well as experimental data. The results showed that the LES and RANS qualitatively well predicted the experimental results, such as mean axial, radial velocities and turbulent kinetic energy. In the quantitative analysis, however, the LES showed a better prediction performance than RANS. Specially, the LES well described characteristics of the recirculation zones, such as air stagnation point and jet stagnation point. Finally, the unsteady phenomena on the Bluff-body, such as the transition of recirculation region and vorticity, was examined with LES methodology.

• Journal of Power System Engineering
• /
• v.19 no.6
• /
• pp.19-25
• /
• 2015

The unsteady-state, incompressible and three-dimensional large eddy simulation(LES) was carried out to analyze the aerodynamic performance of three-dimensional small-size axial fan(SSAF) with the different depth of bellmouth. The static pressure coefficients analyzed by LES predict a little bit larger than measurements except stall region regardless of the installation depth between SSAF and bellmouth. Moreover, static pressure efficiencies analyzed by LES show about maximum 30% at the actual operating point ranges, but measurements do not. Therefore, if the blades of conventional SSAF have some more rigidity and complete dynamic balance, the aerodynamic performance of SSAF will be some more improved. In consequence, LES shows the best prediction performance in comparison with any other Reynolds averaged Navier-Stokes(RANS) method.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.2
• /
• pp.162-170
• /
• 2014

A CFD(Computational Fluid Dynamics) analysis is presented to predict hydrodynamic characteristics of a marine propeller. A commercial RANS(Reynolds Averaged Navier-Stokes equation) solver, namely FLUENT, is utilized in conjunction with fully unstructured meshes around rotating propeller. Mesh generation process is greatly accelerated by using fully unstructured meshes composed of both isotropic and anisotropic tetrahedral elements. The anisotropic tetrahedral elements were used in the flow domain near the blade and shaft, where the viscous effect is important, having complex shape yet resolving the thin boundary layers. For other regions, isotropic tetrahedral elements are utilized. Two different approaches simulating rotational effect of the propeller are employed, namely Moving reference frame technique for steady simulation, and Sliding mesh technique for unsteady simulation. Both approaches are applied to the propeller open water (POW) test simulation. The current results, which are thrust and torque coefficients, are compared with available experimental data.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.12 no.6
• /
• pp.48-55
• /
• 2008

We have simulated the performance of a simple engine model with a gas turbine engine simulation program based on CFD. 2-dimensional Navier-Stokes code for the viscous flow was applied to simulate a compressor and a turbine, and the chemical equilibrium code with the lumped method was applied to simulate the combustor. Unsteady-flow phenomenon between rotor and stator of the compressor and the turbine was analyzed by steady mixing-plane method. In this way, the influence of the turbine blade pitch on the engine was investigated. It was shown that the compressor is operated at more higher pressure conditions as narrower the pitch distance of the turbine.