• Journal of Korea Water Resources Association
• /
• v.50 no.2
• /
• pp.89-97
• /
• 2017

Turbulent flow structure in the high amplitude meandering channel is complex due to secondary recirculation with helicoidal motions and shear layers formed by flow separation from the curved sidewall. In this work, the secondary flow and the superelevation of the water surface produced in the high-amplitude Kinoshita channel are reproduced by the unsteady Reynolds-averaged Navier-Stokes (RANS) computations using the VOF technique for resolving the variation of water surface elevation and three statistical turbulence models ($k-{\varepsilon}$, RNG $k-{\varepsilon}$, $k-{\omega}$ SST). The numerical results computed by a second-order accurate finite volume method are compared with an existing experimental measurement. Among applied turbulence models, $k-{\omega}$ SST model relatively well predicts overall distribution of the secondary recirculation in the Kinoshita channel, while all three models yield similar prediction of water superelevation transverse slope. The secondary recirculation driven by the radial acceleration in the upstream bend affects the flow structure in the downstream bend, which yields a pair of counter-rotating vortices at the bend apex. This complex flow pattern is reasonably well reproduced by the $k-{\omega}$ SST model. Both $k-{\varepsilon}$ based models fail to predict the clockwise-rotating vortex between a pair of counter-rotating vortices which was observed in the experiment. Regardless of applied turbulence models, the present computations using the VOF method appear to well reproduce the superelevation of water surface through the meandering channel.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.2
• /
• pp.63-70
• /
• 2014

Design sensitivity analysis and topology design optimization for a piezoelectric crystal resonator are developed. The piezoelectric crystal resonator is deformed mechanically when subjected to electric charge on the electrodes, or vice versa. The Mindlin plate theory with higher-order interpolations along thickness direction is employed for analyzing the thickness-shear vibrations of the crystal resonator. Thin electrode plates are masked on the top and bottom layers of the crystal plate in order to enforce to vibrate it or detect electric signals. Although the electrode is very thin, its weight and shape could change the performance of the resonators. Thus, the design variables are the bulk material densities corresponding to the mass of masking electrode plates. An optimization problem is formulated to find the optimal topology of electrodes, maximizing the thickness-shear contribution of strain energy at the desired motion and restricting the allowable volume and area of masking plates. The necessary design gradients for the thickness-shear frequency(eigenvalue) and the corresponding mode shape(eigenvector) are computed very efficiently and accurately using the analytical design sensitivity analysis method using the eigenvector expansion concept. Through some demonstrative numerical examples, the design sensitivity analysis method is verified to be very efficient and accurate by comparing with the finite difference method. It is also observed that the optimal electrode design yields an improved mode shape and thickness-shear energy.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.19 no.2
• /
• pp.265-281
• /
• 2017

For shallow tunnel constructions, settlement of the ground surface is a main issue. Recent technical developments in shield TBM tunneling technique have enabled a decrease in such settlements based on tunneling with ground deformation controls. For this objective, the tail void grouting is a common practice. Generally surface settlements in a soil of low permeability occur during a tunnel construction but also during a long period after completion of the tunnel. The long-term settlements occur mainly due to consolidation around the tunnel. The consolidation process is caused and determined by the tail void grouting which leads to an excess pore water pressure in the vicinity of the tunnel. Because of this, the grouting pressure has a strong effect on the long-term settlements in the shield tunneling. In order to investigate this effect, a series of coupled hydro-mechanical 3D finite element simulations have been performed. The results show that an increase in grouting pressure reduces the short-term settlements, but in many cases, it doesn't lead to a reduction of the final settlements after the completion of consolidation. Thereby, the existence of a critical grouting pressure is identified, at which the minimal settlements are expected.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.5
• /
• pp.441-449
• /
• 2015

In this paper we perform a linearized buckling analysis using the Kirchhoff plate theory and the von Karman nonlinear strain-displacement relation. Design sensitivity analysis(DSA) expressions for plane elasticity and buckling problems are derived with respect to Young's modulus and thickness. Using the design sensitivity, we can formulate the topology optimization method for minimizing the compliance and maximizing eigenvalues. We develop a topology optimization method applicable to plate buckling problems using the prestress for buckling analysis. Since the prestress is needed to assemble the stress matrix for buckling problem using the von Karman nonlinear strain, we introduced out-of-plane motion. The design variables are parameterized into normalized bulk material densities. The objective functions are the minimum compliance and the maximum eigenvalues and the constraint is the allowable volume. Through several numerical examples, the developed DSA method is verified to yield very accurate sensitivity results compared with the finite difference ones and the topology optimization yields physically meaningful results.

• 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 Society of Naval Architects of Korea
• /
• v.55 no.3
• /
• pp.252-264
• /
• 2018

The calculation of steady-state cavitating flows around Supercavitating Underwater Bodies (SUB's), which consist of a circular disk head (cavitator), a conical fore-body, a cylindrical middle-body and either a boat-tail or a flare-tail, are carried out. To calculate the axisymmetric cavitating flow, used is a commercial computational fluid dynamics code based on the finite volume method, Fluent. From the analysis of numerical results, the cavity and drag, affected by the fore-body and tail of the SUB's, are investigated. Firstly, the effect of the fore-body shape is investigated with the same disk cavitator and a cylindrical rear-body of fixed diameter. Then with the same cavitator and a fixed fore-body, the effect of the rear-body shape is investigated. Before the cavity generated by the cavitator covers the slant of fore-bodies sufficiently, the larger the cone angle of the fore-body(i.e., the shorter the slant length), the larger the drag and the slower the development of cavity. After the cavity covers the fore-body completely so that the pressure drag component of the body is vanished, the characteristics of drag-velocity curves are identical. Also, as the tail angle is bigger, the cavity generated by the cavitator is suppressed further and the drag becomes larger. The peak of the drag appears for the flare-tail, i.e., when the tail angle is positive(+). On the contrary, the trough of the drag appears for the boat-tail, i.e., when the tail angle is negative(-). When the tail angle is 5 degrees, the peak of the drag appears at the body speed of 80m/s and the value of the drag is 43% larger than that at the design speed of 100m/s. When the tail angle is -5 degrees, the trough of the total drag appears at 75m/s and that drag is 30% smaller than that of the cavitator, which means the rest of the body has a negative drag.

• Journal of the Korean Vacuum Society
• /
• v.19 no.4
• /
• pp.249-255
• /
• 2010

Methods of the characteristics evaluation of turbo-molecular pumps (TMP) are well-defined in the international measurement standards such as ISO, PNEUROP, DIN, JIS, and AVS. The Vacuum Center in the Korea Research Institute of Standards and Science (KRISS) has recently designed, constructed, and established the integrated characteristics evaluation system of TMPs based on the international documents by continuously pursuing and acquiring the reliable international credibility through measurement perfection. The measurement of TMP pumping speed is normally performed with the throughput and orifice methods dependent on the mass flow regions. However, in the UHV range of the molecular flow region, the high uncertainties of the gauges, mass flow rates, and conductance are too critical to precisely accumulate reliable data. In order to solve the uncertainty problems of pumping speeds in the UHV range, we introduced a SRG with 1% accuracy and a constant volume flow meter (CVFM) to measure the finite mass flow rates down to $10^{-1}$ Pa-L/s with 3% uncertainty for the throughput method. In this way we have performed the measurement of pumping speed down to $10^{-4}$ Pa with an uncertainty of less than 6% for a 1000 L/s TMP. In this article we suggest that the CVFM has an ability to measure the conductance of the orifice experimentally with flowing the known mass through the orifice chambers, so that we may overcome the discontinuity problem encountering during introducing two measurement methods in one pumping speed evaluation sequence.

• Journal of Korea Water Resources Association
• /
• v.45 no.1
• /
• pp.39-52
• /
• 2012

Surface-subsurface interactions are an intrinsic component of the hydrologic response within a watershed. In general, these interactions are considered to be one of the most difficult areas of the discipline, particularly for the modeler who intends simulate the dynamic relations between these two major domains of the hydrological cycle. In essence, one major complexity is the spatial and temporal variations in the dynamically interacting system behavior. The proper simulation of these variations requires the need for providing an appropriate coupling mechanism between the surface and subsurface components of the system. In this study, an approach for modelling surface-subsurface flow and transport in a fully intergrated way is presented. The model uses the 2-dimensional diffusion wave equation for sheet surface water flow, and the Boussinesq equation with the Darcy's law and Dupuit-Forchheimer's assumption for variably saturated subsurface water flow. The coupled system of equations governing surface and subsurface flows is discretized using the finite volume method with central differencing in space and the Crank-Nicolson method in time. The interactions between surface and subsurface flows are considered mass balance based on the continuity conditions of pressure head and exchange flux. The major module consists of four sub-module (SUBFA, SFA, IA and NS module) is developed.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.12
• /
• pp.5-15
• /
• 2012

In this study, a series of CFD (Computational Fluid Dynamics) numerical analyses were performed in order to evaluate the thermal performance of six full-scale closed-loop vertical ground heat exchangers constructed in a test bed located in Wonju. The circulation HDPE pipe, borehole and surrounding ground formation were modeled using FLUENT, a finite-volume method (FVM) program, for analyzing the heat transfer process of the system. Two user-defined functions (UDFs) accounting for the difference in the temperatures of the circulating inflow and outflow fluid and the variation of the surrounding ground temperature with depth were adopted in the FLUENT model. The relevant thermal properties of materials measured in laboratory were used in the numerical analyses to compare the thermal efficiency of various types of the heat exchangers installed in the test bed. The simulation results provide a verification for the in-situ thermal response test (TRT) data. The CFD numerical back-analysis with the ground thermal conductivity of 4 W/mK yielded better agreement with the in-situ thermal response tests than with the ground thermal conductivity of 3 W/mK.

• Journal of the Korea Society for Simulation
• /
• v.27 no.4
• /
• pp.19-26
• /
• 2018

The aerial launching UAV(Unmanned Aerial Vehicle) mainly uses a set of folding tandem wings to maximize flight performance and minimize the space required for mounting in a mothership. This folding tandem wing has a unique aerodynamic problem that is different from the general type of fixed wing aircraft, such as the rear wing interference problem caused by the wing of the front wing wake and vortex, and the imbalance of the pivot moment applied to the front and rear wings when the wing is deployed. In this paper, we have modeled and simulated various cases through computational fluid dynamics based on the finite volume method and analyzed various aerodynamic phenomena of the tandem wing type aircraft. We find that the front wing shall be installed higher than the rear for minimizing the wake influence and the rear wing can be deployed faster than the front because of the pivot moment due to aerodynamic forces. Also, considering the pivot moment due to aerodynamic force, the rear wing can be deployed much faster than the front wing. Therefore, it is necessary to consider it when developing the wing deploy mechanism.