• 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 Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.60-68
• /
• 1991

The updated Lagrangian Finite Element Method is introduced to analyse rigid body-fluid impact problem which is characterized by incompressible Navier-Stokes equations and impact-contact conditions between free surface and rigid body. For the convenience of numerical computation, velocity fields are splinted into vicous and pressure parts, and then the governing equations and boundary conditions are decomposed in accordance with the decomposition. However, Viscous stresses acting an the solid boundaries are neglected on the assumption that very small velocity gradients may occur during extremely small time interval of the impact. Four coded quadrilateral elements are used to discretize the space domain and the fully explicit time-marching algorithm is employed with a reasonably small time step. At the beginning of each time step, contact velocity of the rigid body is computed from the momentum balance between the body and the fluid. The velocity field is then computed to satisfy the discretized equations of motions and incompressibility and contact constraints as well as an exact free surface boundary condition. At the end of each time step, the fluid domain is updated from the velocity field. In the present time stepping numerical analysis, behaviour of the free surface near the body can be observed without any difficulty which is very important in the water impact problem. The applicability of the algorithm is illustrated by a wedge type falling body problem. The numerical solutions for time-varying pressure distributions and impact loadings acting ion the surface are obtained.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.3
• /
• pp.180-198
• /
• 2019

In recent years, countries like Europe and Japan have been involved in many researches on the Low-Crested Structure (LCS) which is the method to protect beach erosion and it is regarded as an alternative to the submerged breakwaters, and compiled its results and released the design manual. In the past, studies on LCS have focused on two-dimensional wave transmission and calculating required weight of armor units, and these were mainly examined and discussed based on experiments. In this study, three-dimensional numerical analysis is performed on permeable LCS. The open-source CFD code olaFlow based on the Navier-Stokes momentum equations is applied to the numerical analysis, which is a strongly nonlinear analysis method that enables breaking and turbulence analysis. As a result, the distribution characteristics of the LCS such as water level, water flow, and turbulent kinetic energy were examined and discussed, then they were carefully compared and examined in the case of submerged breakwaters. The study results indicate that there is a difference between the flow patterns of longshore current near the shoreline, the spatial distribution of longshore and on-offshore directions of mean turbulent kinetic energy in case of submerged breakwaters and LCS. It is predicted that the difference in these results leads to the difference in sand movement.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.47 no.9
• /
• pp.649-657
• /
• 2019

Hydrazine, which is used as a representative monopropellant, is an extremely poisonous substance and has a disadvantage that it is harmful to the human body and is very difficult to handle. In recent years, research on the development of non-toxic and environmentally friendly propellants has attracted much attention. Ammonium dinitramide(ADN) based propellant developed by Swedish Space Corporation has superior performance to hydrazine and has been commercialized through performance verification in space environment. On the other hand, the exhaust gas from a thruster nozzle collides with a satellite while it is spreading in the vacuum space, thermal load and surface contamination may occur and may reduce the performance and lifetime of the satellite. However, a study on the effect of the exhaust gas of the green propellant thruster on the satellite has not been conducted in earnest yet. Therefore, the exhaust gas behavior in space was analyzed in this study for the ADN based green monopropellant using Navier-Stokes equations and the DSMC method. As a result, it can be expected to be used as design validation data in the development of satellite when using the ADN based green monopropellant.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.127-136
• /
• 2022

In this work a multi-fidelity non-intrusive polynomial chaos (MF-NIPC) has been applied to a structural wind engineering problem in architectural design for the first time. In architectural design it is important to design structures that are safe in a range of wind directions and speeds. For this reason, the computational models used to design buildings and bridges must account for the uncertainties associated with the interaction between the structure and wind. In order to use the numerical simulations for the design, the numerical models must be validated by experi-mental data, and uncertainties contained in the experiments should also be taken into account. Uncertainty Quantifi-cation has been increasingly used for CFD simulations to consider such uncertainties. Typically, CFD simulations are computationally expensive, motivating the increased interest in multi-fidelity methods due to their ability to lev-erage limited data sets of high-fidelity data with evaluations of more computationally inexpensive models. Previous-ly, the multi-fidelity framework has been applied to CFD simulations for the purposes of optimization, rather than for the statistical assessment of candidate design. In this paper MF-NIPC method is applied to flow around a rectan-gular 5:1 cylinder, which has been thoroughly investigated for architectural design. The purpose of UQ is validation of numerical simulation results with experimental data, therefore the radius of curvature of the rectangular cylinder corners and the angle of attack are considered to be random variables, which are known to contain uncertainties when wind tunnel tests are carried out. Computational Fluid Dynamics (CFD) simulations are solved by a solver that employs the Finite Element Method (FEM) for two turbulence modeling approaches of the incompressible Navier-Stokes equations: Unsteady Reynolds Averaged Navier Stokes (URANS) and the Large Eddy simulation (LES). The results of the uncertainty analysis with CFD are compared to experimental data in terms of time-averaged pressure coefficients and bulk parameters. In addition, the accuracy and efficiency of the multi-fidelity framework is demonstrated through a comparison with the results of the high-fidelity model.

• The Journal of the Acoustical Society of Korea
• /
• v.40 no.6
• /
• pp.555-565
• /
• 2021

The goal of this study is to improve flow and noise performances of existing backward-curved blade centrifugal fan system used for circulating cold air in a refrigerator freezer by optimally designing airfoil shape. The unique characteristics of the system is to drive cold airflow with two volute tongues in combination with duct system in a back side of a refrigerator without scroll housing generally used in a typical centrifugal fan system. First, flow and noise performances of existing fan system were evaluated experimentally. A P-Q curve was obtained using a fan performance tester in the flow experiment, and noise spectrum was measured in an anechoic chamber in the noise experiment. Then, flow characteristics were numerically analyzed by solving the three-dimensional unsteady Navier-Stokes equations and noise analysis was performed by solving the Ffowcs Williams and Hawkins equation with input from the flow simulation results. The validity of numerical results was confirmed by comparing them with the measured ones. Based on the verified numerical method, blade inlet and outlet angles were optimized for maximum flow rate using the two-factor central composite design of the response surface method. Finally, the flow and noise performances of a prototype manufactured with the optimum design were experimentally evaluated, which showed the improvement in flow and noise performance.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.3B
• /
• pp.261-270
• /
• 2008

Even though there is a great deal of progress in a numerical method of high caliber like SPH, it is very rarely deployed in a water resources community. Despite the great stride in computing environment, depth averaged approach like a nonlinear shallow equation is still efficient tool for flood routing in large watershed, but it can give some misleading information like the inundation height of flood. In this rationale, we numerically simulate the flow into the dry channel, dry channel with an obstacle triggered by the collapse of a two dimensional water column using SPH (Smoothed Particle Hydrodynamics) in order to boost the application of numerical method of high caliber like SPH in a water resources community. As a most severe test of the robustness of SPH, we also carry out the simulation of the flow through a clearance into the wet channel driven by the rapid removal of a water gate. As a hydrodynamic model, we used the Navier-Stokes equation, a numerical integration of which was carried out using SPH. To verify the validity of newly proposed numerical model, we compare the numerically simulated flow with the others in the literature mainly from VOF and MAC, and hydraulic experiments by Martin and Moyce (1952), Koshizuka et al. (1995) and Janosi et al. (2004). It was shown that agreements between the numerical results in this study and hydraulic experiments are remarkable.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.2
• /
• pp.147-157
• /
• 2011

A RANS(Reynolds averaged Navier-Stokes) based numerical method is developed for the evaluation of ship resistance and self-propulsion performances. In the usability aspect of CFD for the hull form design, the field grid around practical hull forms is generated by solving a grid Poisson equation based on the hull surface grid generated from station offsets and centerline profile. A body force technique is introduced to model the effects of the propeller in which the propeller loads are obtained from potential flow analysis using an unsteady lifting surface method. The free surface is captured by using a two-phase level-set method and the realizable $k-{\varepsilon}$ model is used for turbulence closure. The hull attitude in vertical plane, i.e., trim and sinkage, is calculated by using a quasi-steady method and then considered in the computation by translating and rotating the grid system according to the values. For the validation of the proposed method, the numerical results of resistance tests for KCS, KLNG, and KVLCC1 and of self-propulsion test for KCS are compared with experimental data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.6
• /
• pp.429-434
• /
• 2017

A molecular-continuum hybrid method was developed to simulate microscale and nanoscale fluids where continuum fluidics cannot be used to predict Couette flow. Molecular dynamics simulation is used near the solid surface where the flow cannot be predicted by continuum fluidics, and Navier-Stokes equations are used in the other regions. Numerical simulation of Couette flow was performed using the hybrid method to investigate the effect of solid-liquid interaction and surface roughness in a nanochannel. It was found that the solid-liquid interaction and surface roughness influence the boundary condition. When the surface energy is low, slippage occurs near the solid surface, and the magnitude of slippage decreases with increase in surface energy. When the surface energy is high, a locking boundary condition is formed. The roughness disturbs slippage near the solid surface and promotes the locking boundary condition.

• Steel and Composite Structures
• /
• v.19 no.4
• /
• pp.861-876
• /
• 2015

Composite sections design consists on checking that the point defined by axial load and bending moment keeps included within the surface enclosed by the section interaction curve. Eurocode 4 suggests a method for tracing this diagram based on the plastic stress distribution method. However curves obtained according to this criterion overvalue concrete encased sections bearing capacity, especially when axial force comes with high bending moment values, so a correction factor is required. This article proposes a method for tracing this diagram based on the strain compatibility method. When stresses on the section are integrated by considering the Navier hypothesis, the use of the materials nonlinear constitutive equations provides curves much more adjusted to reality. This process requires the use of rather complex software which might reveal as too complex for practitioners. Preserving the same criteria of an elastic-plastic stress distribution, this article presents alternative expressions to obtain the failure internal forces in five significant points of the interaction diagram having considered five different positions of the neutral axis. These expressions are simply enough for their practical application. Concordance of curves traced strictly relying on these five points with those obtained by computer assisted stress integration considering the strain compatibility method and even with Eurocode 4 weighted curves will be presented for three different cross-sections and two different concrete strengths, revealing very good results.