• Journal of the Institute of Convergence Signal Processing
• /
• v.24 no.2
• /
• pp.97-102
• /
• 2023

Numerical simulations have been performed to investigate the hydrodynamic characteristics of the rudder since they play an important role in naval architecture fields. Although some values such as hydrodynamics forces can be measured easily in the towing tanks, it is difficult to obtain the detailed information of the flow fields such as pressure distribution, velocity distribution, vortex generation from experiments. In the present study, the effects of hydrodynamic coefficients and Reynolds number acting on the rudder were studied by using Computational Fluid Dynamics(CFD). Ansys fluent, one of commercial CFD solvers, solves the Navier-Stokes equations and the k-epsilon turbulence model is selected for the viscous model to solve RANS equations. At first, drag coefficients and lift coefficient for different angle of attack are obtained by using a CFD commercial code for KCS rudder. Secondly, the 2-D lift coefficients and drag coefficients are compared with 3-D coefficients at the same conditions. Thirdly, the effects of Reynolds number on the hydrodynamic forces are investigated.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.11-23
• /
• 1997

This paper contains a new approach to blade section design method for marine propellers. The hydrodynamic characteristics of 2-D section are highly influenced by its geometrical parameters i.e., thickness and camber distributions and leading edge radius etc. To consider fully turbulent flow field near 2-D section. the finite volume method with k-${\varepsilon}$ turbulent model which solve Reynolds time averaged Navier-Stokes(RANS) equation is applied. In this study, O-type grid system that can provide many calculation points on blade surface is used. The results were compared with those of the experiment of NACA0012 to confirm the accuracy of the developed codes. The goal of this study is the development of a blade section with high efficiency and low drag. To achieve this, we carried out the tests of lift, drag and cavitation characteristics in cavitation tunnel. The results of experiment were compared with numerical results in order to validate the proposed blades design method. By comparing the numerical results with the experiments, we found that the new blade section, KH28 allows superior performance in efficiency and cavitation avoidance characteristics. We further investigated the blade section design method and an application study of this section, KH28 to apply to the marine propeller. In order to improve the accuracy of numerical results on prediction of lift and drag, we conclude here that the 2-layer boundary model must be used.

• Journal of the Korean Electrochemical Society
• /
• v.14 no.4
• /
• pp.231-238
• /
• 2011

In this work, the numerical analysis for the zincate behavior at a zinc electrode with an electrolyte flow was carried out for a ZAFC. The Nernst-Planck equation with a boundary condition of Butler-Volmer type was adopted to describe electrochemical effects of mass transfer, migration, kinetics of electrode. The Navier-Stokes equation, coupling to the Nernst-Planck equation, is also applied to describe the internal electrolyte flow fields. The validity of the numerical model is proved through the comparative analysis between numerical and experimental results. The concentration of zincate and the current density were also investigated at a zinc anode according to various electrolyte velocities. We have found the concentration of zincate decreased and the current density increased with an increase in the electrolyte velocity.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.41 no.3
• /
• pp.173-184
• /
• 2013

A multi-level design optimization framework for aerodynamic design of rotary wing such as propeller and helicopter rotor blades is presented in this study. Strategy of the proposed framework is to enhance aerodynamic performance by sequentially applying the planform and sectional design optimization. In the first level of a planform design, we used a genetic algorithm and blade element momentum theory (BEMT) based on two-dimensional aerodynamic database to find optimal planform variables. After an initial planform design, local flow conditions of blade sections are analyzed using high-fidelity CFD methods. During the next level, a sectional design optimization is conducted using two dimensional Navier-Stokes analysis and a gradient based optimization algorithm. When optimal airfoil shape is determined at the several spanwise locations, a planform design is performed again. Through this iterative design process, not only an optimal flow condition but also an optimal shape of an EAV propeller blade is obtained. To validate the optimized propeller-blade design, it is tested in wind-tunnel facility with different flow conditions. An efficiency, which is slightly less than the expected improvement of 7% predicted by our proposed design framework but is still satisfactory to enhance the aerodynamic performance of EAV system.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.5
• /
• pp.323-334
• /
• 2020

Since boundary layer transition has a significant impact on the aero-thermodynamic performance of hypersonic flight vehicles, capability of accurate prediction of transition location is essential for design and performance analysis. In this study, γ-Re_θt model is improved to predict transition of hypersonic boundary layers and validated. A coefficient in the production term of the intermittency transport equation that affects the transition onset location is constructed and applied as a function of Mach number, wall temperature, and freestream stagnation temperature based on the similarity numerical solution of compressible boundary layer. To take into account a Mach number dependency of transition onset momentum thickness Reynolds number and transition length, additional correlation equations are determined as function of Mach number and applied to Re_θc and F_length correlations of the baseline model. The suggested model is implemented to a commercial CFD code in consideration of practical use. Analysis of hypersonic flat plate and circular cone boundary layers is carried out by using the model for validation purpose. An improvement of prediction capability with respect to variation of Mach number and unit Reynolds number is identified from the comparison with experimental data.

• Wind and Structures
• /
• v.20 no.4
• /
• pp.549-564
• /
• 2015

The three-dimensional unsteady incompressible Reynolds-averaged Navier-Stokes equations and k-${\varepsilon}$ double equations turbulent model were used to investigate the effect on the measurements of anemometers due to a passing high-speed train. Sliding mesh technology in Fluent was utilized to treat the moving boundary problem. The high-speed train considered in this paper was with bogies and inter-carriage gaps. Combined with the results of the wind tunnel test in a published paper, the accuracy of the present numerical method was validated to be used for further study. In addition, the difference of slipstream between three-car and eight-car grouping models was analyzed, and a series of numerical simulations were carried out to study the influences of the anemometer heights, the train speeds, the crosswind speeds and the directions of the induced slipstream on the measurements of the anemometers. The results show that the influence factors of the train-induced slipstream are the passing head car and tail car. Using the three-car grouping model to analyze the train-induced flow is reasonable. The maxima of horizontal slipstream velocity tend to reduce as the height of the anemometer increases. With the train speed increasing, the relationship between $V_{train}$ and $V_{induced\;slipstream}$ can be expressed with linear increment. In the absence of natural wind conditions, from the head car arriving to the tail car leaving, the induced wind direction changes about $330^{\circ}$, while under the crosswind condition the wind direction fluctuates around $-90^{\circ}$. With the crosswind speed increasing, the peaks of $V_X,{\mid}V_{XY}-V_{wind}{\mid}$ of the head car and that of $V_X$ of the tail car tend to enlarge. Thus, when anemometers are installed along high-speed railways, it is important to study the effect on the measurements of anemometers due to the train-induced slipstream.

• Ocean Systems Engineering
• /
• v.6 no.1
• /
• pp.61-97
• /
• 2016

This is the second of two papers on the 3D numerical modeling of nearshore hydro- and morphodynamics. In Part I, the focus was on surf and swash zone hydrodynamics in the cross-shore and longshore directions. Here, we consider nearshore processes with an emphasis on the effects of oceanic forcing and beach characteristics on sediment transport in the cross- and longshore directions, as well as on foreshore bathymetry changes. The Delft3D and XBeach models were used with four turbulence closures (viz., ${\kappa}-{\varepsilon}$, ${\kappa}-L$, ATM and H-LES) to solve the 3D Navier-Stokes equations for incompressible flow as well as the beach morphology. The sediment transport module simulates both bed load and suspended load transport of non-cohesive sediments. Twenty sets of numerical experiments combining nine control parameters under a range of bed characteristics and incident wave and tidal conditions were simulated. For each case, the general morphological response in shore-normal and shore-parallel directions was presented. Numerical results showed that the ${\kappa}-{\varepsilon}$ and H-LES closure models yield similar results that are in better agreement with existing morphodynamic observations than the results of the other turbulence models. The simulations showed that wave forcing drives a sediment circulation pattern that results in bar and berm formation. However, together with wave forcing, tides modulate the predicted nearshore sediment dynamics. The combination of tides and wave action has a notable effect on longshore suspended sediment transport fluxes, relative to wave action alone. The model's ability to predict sediment transport under propagation of obliquely incident wave conditions underscores its potential for understanding the evolution of beach morphology at field scale. For example, the results of the model confirmed that the wave characteristics have a considerable effect on the cumulative erosion/deposition, cross-shore distribution of longshore sediment transport and transport rate across and along the beach face. In addition, for the same type of oceanic forcing, the beach morphology exhibits different erosive characteristics depending on grain size (e.g., foreshore profile evolution is erosive or accretive on fine or coarse sand beaches, respectively). Decreasing wave height increases the proportion of onshore to offshore fluxes, almost reaching a neutral net balance. The sediment movement increases with wave height, which is the dominant factor controlling the beach face shape.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.24 no.2
• /
• pp.187-207
• /
• 2019

In order to develop a high performance ocean model, we used Julia, a Just-In-Time compile language, and to obtain the solution of the momentum equation, we made the code to solve the Poisson equation by the Successive Over-Relaxation method. And then we made two models to test Julia calculation codes. First, a simple channel form is modeled to test constant source/sink conditions. Second, the simplified Yellow Sea was modeled to test tidal forcing, Coriolis forces, and the effect of vertical eddy diffusivity coefficients. The model has been tested with a total of eight cases in the two scenarios. As a result of the test, the depth-averaged current speed of the three cases in Scenario 1 converged perfectly to the theoretical value, and that showed well a vertical flow velocity gradient due to the bottom friction. Also, the result of Scenario 2 represented well the amphidromic points of Yellow Sea and the tidal characteristics of mid-western and southwestern coast of Korea. Therefore, it is considered that the ocean model using Julia language has developed successfully, this suggests that the ocean model has come to the stage of successful transition from a classical compile language to a Just-In-Time compile language.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.1B
• /
• pp.95-109
• /
• 2008

We numerically analyzed the nonlinear shoaling, a plunging breaker and its accompanying energetic suspension of sediment at a bed, and a redistribution of suspended sediments by a down rush of preceding waves and the following plunger using SPH with a Gaussian kernel function, Lagrangian Dynamic Smagorinsky model (LDS), Van Rijn's pick up function. In that process, we came to the conclusion that the conventional model for the tractive force at a bottom like a quadratic law can not accurately describe the rapidly accelerating flow over a swash zone, and propose new methodology to accurately estimate the bottom tractive force. Using newly proposed wave model in this study, we can successfully duplicate severely deformed water surface profile, free falling water particles, a queuing splash after the landing of water particles on the free surface and a wave finger due to the structured vortex on a rear side of wave crest (Narayanaswamy and Dalrymple, 2002), a circulation of suspended sediments over a swash zone, net transfer of sediments clouds suspended over a swash zone toward the offshore, which so far have been regarded very difficult features to mimic in the computational fluid mechanics.