• Nuclear Engineering and Technology
• /
• v.51 no.2
• /
• pp.573-578
• /
• 2019

Anisotropic distribution of the turbulent kinetic energy and the near-field excitations are the main causes of the steady state Flow-Induced Vibration (FIV) which could lead to fretting wear damage in vertically arranged supported slender rods. In this article, a combined Computational Fluid Dynamics (CFD) and Computational Structural Mechanic (CSM) approach named two-way Fluid-Structure Interaction (FSI) is used to investigate the modal characteristics of a typical rod's vibration. Performance of an Unsteady Reynolds-Average Navier-Stokes (URANS) and Large Eddy Simulation (LES) turbulence models on asymmetric fluctuations of the flow field are investigated. Using the LES turbulence model, any large deformation damps into a weak oscillation which remains in the system. However, it is challenging to use LES in two-way FSI problems from fluid domain discretization point of view which is investigated in this article as the innovation. It is concluded that the near-wall meshes whiten the viscous sub-layer is of great importance to estimate the Root Mean Square (RMS) of FIV amplitude correctly as a significant fretting wear parameter otherwise it merely computes the frequency of FIV.

• Computers and Concrete
• /
• v.27 no.2
• /
• pp.85-97
• /
• 2021

In this study, the effect of steel fiber utilization, boundary conditions, different beam cross-section, and length parameter are investigated on the free vibration behavior of fiber reinforced self-compacting concrete beam on elastic foundation. In the analysis of the beam model recommended by Euler-Bernoulli, a method utilizing Stokes transformations and Fourier Sine series were used. For this purpose, in addition to the control beam containing no fiber, three SCC beam elements were prepared by utilization of steel fiber as 0.6% by volume. The time-dependent fresh properties and some mechanical properties of self-compacting concrete mixtures were investigated. In the modelled beam, four different beam specimens produced with 0.6% by volume of steel fiber reinforced and pure (containing no fiber) SCC were analyzed depending on different boundary conditions, different beam cross-sections, and lengths. For this aim, the effect of elasticity of the foundation, cross-sectional dimensions, beam length, boundary conditions, and steel fiber on natural frequency and frequency parameters were investigated. As a result, it was observed that there is a noticeable effect of fiber reinforcement on the dynamic behavior of the modelled beam.

• Journal of The Korean Astronomical Society
• /
• v.56 no.1
• /
• pp.23-33
• /
• 2023

We investigate the diffusion process of Thomson-scattered line photons in both real space and frequency space through a Monte Carlo approach. The emission source is assumed to be monochromatic and point-like embedded at the center of a free electron region in the form of a sphere and a slab. In the case of a spherical region, the line profiles emergent at a location of Thomson optical depth τ_Th from the source exhibit the full width of the half maximum σ_λ ≃ τ^1.5_Th. In the slab case, we focus on the polarization behavior where the polarization direction flips from the normal direction of the slab to the parallel as the slab optical depth τ_Th increases from τ_Th ≪ 1 to τ_Th ≫ 1. We propose that the polarization flip to the parallel direction to the slab surface in optically thick slabs is attributed to the robustness of the Stokes parameter Q along the vertical axis with respect to the observer's line of sight whereas randomization dominates the remaining region as τ_Th increases. A brief discussion on the importance of our study is presented.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.21 no.3
• /
• pp.10-16
• /
• 2020

Computational simulations of adsorption columns were carried out to investigate the removal characteristics of VOCs from a laundry shop. n-Decane was selected as the representative component among the VOCs emitted, and the activity of the adsorbents, such as activated carbon, was evaluated using commercial CFD code. The mathematical framework was composed of continuity and Navier-stokes equations, and the simulation was performed using the Matlab program. The adsorption isotherms of LDF, Freundlich, and Langmuir were evaluated, and the adsorption amount of the adsorption isotherms with the adsorption parameter was compared. The simulation was carried out using a particle porosity, dispersion coefficient, particle density, bed diameter, and bed length of 0.79, 42.4 ㎠/min, 485 g/L, 2.0 cm, and 2.5 cm, respectively. The effect of the gas velocity, dispersion coefficient, and voidage on the adsorption amount was compared in the Langmuir adsorption isotherm. The simulation was carried out in the velocity range of 50 to 200 cm/min, dispersion coefficient range of 100 to 400 ㎠/min, and particle porosity range of 0.66 to 0.79. The simulation results of activated carbon with benzene coincided with the Langmuir isotherm. Three types of adsorption isotherm were compared under similar conditions, and the simulation results showed the efficient adsorption condition for hydrocarbons.

• Transactions of the Korean hydrogen and new energy society
• /
• v.22 no.3
• /
• pp.386-401
• /
• 2011

Mathematical models for various steps in coal gasification reactions were developed and applied to investigate the effects of operation parameters on dynamic behavior of gasification process. Chemical reactions considered in these models were pyrolysis, volatile combustion, water shift reaction, steam-methane reformation, and char gasification. Kinetics of heterogeneous reactions between char and gaseous agents was based on Random pore model. Momentum balance and Stokes' law were used to estimate the residence time of solid particles (char) in an up-flow reactor. The effects of operation parameters on syngas composition, reaction temperature, carbon conversion were verified. Parameters considered here for this purpose were $O_2$-to-coal mass ratio, pressure of reactor, composition of coal, diameter of char particle. On the basis of these parametric studies some quantitative parameter-response relationships were established from both dynamic and steady-state point of view. Without depending on steady state approximation, the present model can describe both transient and long-time limit behavior of the gasification system and accordingly serve as a proto-type dynamic simulator of coal gasification process. Incorporation of heat transfer through heterogenous boundaries, slag formation and steam generation is under progress and additional refinement of mathematical models to reflect the actual design of commercial gasifiers will be made in the near futureK.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.4
• /
• pp.165-173
• /
• 1998

Analytical solutions to predict the movement rate of submerged mound are derived using the convection coefficient and the joint distribution function of wave heights and periods. Assuming that the sediment is moved onshore due to the velocity asymmetry of Stokes' second order nonlinear wave theory, the micro-scale bedload transport equation is applied to the sediment conservation. The nonlinear convection-diffusion equation can then be obtained which governs the migration of submerged mound. The movement rate decreases exponentially with increasing the water depth, but the movement rate tends to increase as the spectral width parameter, $ u$ increases. In comparison of the analytical solution with the measured data, it is found that the analytical solution overestimates the movement rate. However, the agreement between the analytical solution and the measured data is encouraging since this over-estimation may be due to the inaccuracy of input data and the limitation of sediment transport model. In particular, the movement rates with respect to the water depth predicted by the analytical solution are in very good agreement with the estimated result using the discritization technique with the hindcast wave data.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.2
• /
• pp.262-270
• /
• 1987

A numerical method is presented which can solve the unsteady momentum and thermal boundary layers, coupled through the agency of buoyancy force, over a heated circular cylinder impulsively started from rest. By linearizing the nonlinear finite difference equations without sacrificing accuracy, numerical solutions are obtained at each time step without iteration. To get rid of the requirement of excessive number of grid points in the region of reversed flow, special form of transformed variables are used, by which the computational boundary layer thickness is maintained almost constant. These numerical properties enable the method to easily handle the region of reversed flow and how the singularity develops in the interior of the boundary layer. In order to investigated the thermal effects on the skin friction, heat flux, displacement thickness and on the separation, we have successfully solved three different cases of the buoyancy parameter .alpha.(Gr/Re$^{2}$).

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.130-135
• /
• 2008

An aerodynamic optimization design process of multistage axial turbine is presented in this article: first, applying quasi-three dimensional(Q3D) design methods to conduct preliminary design and then adopting modern optimization design methods to implement multistage local optimization. Quasi-three dimensional(Q3D) design methods, which mainly refer to S2 flow surface direct problem calculation, adopt the S2 flow surface direct problem calculation program of Harbin Institute of Technology. Multistage local optimization adopts the software of Numeca/Design3D, which jointly adopts genetic algorithm and artificial neural network. The major principle of the methodology is that the successive design evaluation is performed by using an artificial neural network instead of a flow solver and the genetic algorithms may be used in an efficient way. Flow computation applies three-dimensional viscosity Navier Stokes(N-S) equation solver. Such optimization process has three features: (i) local optimization based on aerodynamic performance of every cascade; (ii) several times of optimizations being performed to every cascade; and (iii) alternate use of coarse grid and fine grid. Such process was applied to optimize a three-stage axial turbine. During the optimization, blade shape and meridional channel were respectively optimized. Through optimization, the total efficiency increased 1.3% and total power increased 2.4% while total flow rate only slightly changed. Therefore, the total performance was improved and the design objective was achieved. The preliminary design makes use of quasi-three dimensional(Q3D) design methods to achieve most reasonable parameter distribution so as to preliminarily enhance total performance. Then total performance will be further improved by adopting multistage local optimization design. Thus the design objective will be successfully achieved without huge expenditure of manpower and calculation time. Therefore, such optimization design process may be efficiently applied to the aerodynamic design optimization of multistage axial turbine.

• 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.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.26 no.7
• /
• pp.931-941
• /
• 2020

There is a growing interest this paper for ocean sensing where autonomous vehicles can play an essential role in assisting engineers, researchers, and scientists with environmental monitoring and collecting oceanographic data. This study was conducted to develop a rigid sail for the autonomous sailing drone. Our study aims to numerically analyze the aerodynamic characteristics of curvy twin sail and compare it with wing sail. Because racing regulations limit the sail shape, only the two-dimensional geometry (2D) was open for an optimization. Therefore, the first objective was to identify the aerodynamic performance of such curvy twin sails. The secondary objective was to estimate the effect of the sail's spacing and shapes. A viscous Navier-Stokes flow solver was used for the numerical aerodynamic analysis. The 2D aerodynamic investigation is a preliminary evaluation. The results indicated that the curvy twin sail designs have improved lift, drag, and driving force coefficient compared to the wing sails. The spacing between the port and starboard sails of curvy twin sail was an important parameter. The spacing is 0.035 L, 0.07 L, and 0.14 L shows the lift coefficient reduction because of dramatically stall effect, while flow separation is improved with spacing is 0.21 L, 0.28 L, and 0.35 L. Significantly, the spacing 0.28 L shows the maximum high pressure at the lower area and the small low pressure area at leading edges. Therefore, the highest lift was generated.