• Journal of the Korean Society of Visualization
• /
• v.18 no.3
• /
• pp.42-51
• /
• 2020

This study performs a numerical simulation of lid driven rectangular cavity flow with different aspect ratios of k = 0.5 to 4 under Reynolds 100, 1,000, 10,000 by using multi-relaxation time (MRT) Lattice Boltzmann Method (LBM). In order to achieve better convergence, well-posed boundary conditions in the domain should be defined such as no-slip condition on side and bottom solid wall surfaces and uniform horizontal velocity on the top of the cavity. This study focuses on the flow inside different shape of rectangular cavity with the aim to observe the effect of the Reynolds number and aspect ratio on the flow characteristics and primary/secondary vortex formation. In order to validate the study, the results have been compared with existing works. The result shows that the Reynolds number and the aspect ratio both has substantial effects on the flow inside the lid-driven rectangular cavity.

• Journal of Ocean Engineering and Technology
• /
• v.36 no.1
• /
• pp.11-20
• /
• 2022

Generally, tsunamis are generated by the rapid crustal movements of the ocean floor. Other factors of tsunami generation include landslides on coastal and ocean floor slopes, glacier collapses, and meteorite collisions. In this study, two numerical analyses were conducted to examine the formation, propagation, and deformation properties of landslide tsunamis. First, LS-DYNA was adopted to simulate the formation and propagation processes of tsunamis generated by dropping rigid bodies. The generated tsunamis had smaller wave heights and wider waveforms during their propagation, and their waveforms and flow velocities resembled those of theoretical solitary waves after a certain distance. Second, after the formation of the landslide tsunami, a tsunami based on the solitary wave approximation theory was generated in a numerical wave tank (NWT) with a computational domain that considered the stability/steady phase. The comparison of two numerical analysis results over a certain distance indicated that the waveform and flow velocity were approximately equal, and the maximum wave pressures acting on the upright wall also exhibited similar distributions. Therefore, an effective numerical model such as LS-DYNA was necessary to analyze the formation and initial deformations of the landslide tsunami, while an NWT with the wave generation method based on the solitary wave approximation theory was sufficient above a certain distance.

• Journal of the Society of Naval Architects of Korea
• /
• v.59 no.3
• /
• pp.141-148
• /
• 2022

A hybrid turbulence model has developed by combining a sub-grid scale model using dynamic k equation in LES with k-𝜔 SST model of RANS equation. To ascertain potential applicability of the hybrid turbulence model, fully developed turbulent channel flows at Re_𝜏=180 have been simulated of which computational domain has a top wall with coarse cells and a bottom wall with fine cells. The streamwise mean velocity and turbulent intensity profiles showed a good agreement with DNS data when using the hybrid model rather than using a single model in k-𝜔 SST or dynamic k equation models. Computational simulations of turbulent flows around KVLCC2 with a pre-swirl duct have been mainly performed using the hybrid turbulence model. Compared to the results obtained from RANS simulation with k-𝜔 SST model as well as LES with dynamic k equation SGS model, turbulent wakes of the duct in the present simulation using the hybrid turbulence model were very similar to that of LES. Also, the resistances acting on hull, rudder and duct in hybrid turbulence model were similar to those in RANS simulation whereas the viscous forces acting on the hull in LES had a significant error due to coarse cells inappropriate to the sub-grid scale model.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.5
• /
• pp.545-565
• /
• 2018

A numerical code is developed based on potential flow theory to investigate nonlinear sloshing in rectangular Liquefied Natural Gas (LNG) tanks under forced excitation. Using this code, internal free-surface elevation and sloshing loads on liquid tanks can be obtained both in time domain and frequency domain. In the mathematical model, acceleration potential is solved in the calculation of pressure on tanks and the artificial damping model is adopted to account for energy dissipation during sloshing. The Boundary Element Method (BEM) is used to solve boundary value problems of both velocity potential and acceleration potential. Numerical calculation results are compared with published results to determine the efficiency and accuracy of the numerical code. Sloshing properties in partially filled rectangular and membrane tank under translational and rotational excitations are investigated. It is found that sloshing under horizontal and rotational excitations share similar properties. The first resonant mode and excitation frequency are the dominant response frequencies. Resonant sloshing will be excited when vertical excitation lies in the instability region. For liquid tank under rotational excitation, sloshing responses including amplitude and phase are sensitive to the location of the center of rotation. Moreover, experimental tests were conducted to analyze viscous effects on sloshing and to validate the feasibility of artificial damping models. The results show that the artificial damping model with modifying wall boundary conditions has better applicability in simulating sloshing under different fill levels and excitations.

• Proceedings of the KSME Conference
• /
• 2001.06d
• /
• pp.97-103
• /
• 2001

The effects of the electrohydrodynamic (EHD) flow and turbulent diffusion on the collection efficiency of a model ESP composed of the plates with a cavity were studied through numerical computation. The electric field and ion space charge density were calculated by the Poisson equation of the electrical potential and the current continuity equation. The EHD flow field was solved by the continuity and momentum equations of the gas phase including the electrical body force induced by the movement of ions under the electric field. The RNG $k-{\varepsilon}$ model was used to analyze the turbulent flow. The particle concentration distribution was calculated from the convective diffusion equation of the particle phase. As the ion space charge increased, the particulate collection efficiency increased because the electrical potential increased over the entire domain in the ESP. The collection efficiency decreased and then increased, i.e. had a minimum value, as the EHD circulating flow became stronger when the electrical migration velocity of the charged particle was low. However, the collection efficiency decreased with the stronger EHD flow when the electrical migration of the particle was higher relatively. The collection efficiency of the model ESP increased as the turbulent diffusivity of the particle increased when the electrical migration velocity of the particle was low. However, the collection efficiency decreased for increasing the turbulent diffusivity when the electrical migration of the particle was higher relatively.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.3
• /
• pp.167-174
• /
• 2021

Recent investigation into the integrity of nuclear containment buildings has highlighted the importance of developing an elaborate diagnostic method to evaluate the distribution and size of cavities inside concrete walls. As part of developing such a method, this paper presents a finite element approach to modeling elastic waves propagating in the containment building walls of a nuclear power plant. We introduce a perfectly matched layer (PML) wave-absorbing boundary to limit the large-scale nuclear containment wall to the region of interest. The formulation results in a semi-discrete form with symmetric damping and stiffness matrices. The transient elastic wave equations for a mixed unsplit-field PML were solved for displacement and stresses in the time domain. Numerical results show that the sensitivity of displacement, velocity, acceleration, and stresses is large depending on the size and location of the cavity. The dynamic response of the wall slightly differs depending on the existence of the containment liner plate. The results of this study can be applied to a full-waveform inversion approach for characterizing cavities inside a containment wall.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.35 no.5
• /
• pp.299-308
• /
• 2022

This study was conducted to predict the tendency for heat exchange and boil-off gas (BOG) in a liquefied hydrogen tank under sloshing excitation. First, athe fluid domain excited by sloshing was modeled using a multiphase-thermal flow domain in which liquid hydrogen and hydrogen gas are in the saturated state. Both the the volume of fluid (VOF) and Eulerian-based multi-phase flow methods were applied to validate the accuracy of the pressure prediction. Second, it was indirectly shown that the fluid velocity prediction could be accurate by comparing the free surface and impact pressure from the computational fluid dynamics with those from the experimental results. Thereafter, the heat ingress from the external convective heat flux was reflected on the outer surfaces of the hydrogen tank. Eulerian-based multiphase-heat flow analysis was performed for a two-dimensional Type-C cylindrical hydrogen tank under rotational sloshing motion, and an inflation technique was applied to transform the fluid domain into a computational grid model. The heat exchange and heat flux in the hydrogen liquid-gas mixture were calculated throughout the analysis,, whereas the mass transfer and vaporization models were excluded to account for the pure heat exchange between the liquid and gas in the saturated state. In addition, forced convective heat transfer by sloshing on the inner wall of the tank was not reflected so that the heat exchange in the multiphase flow of liquid and gas could only be considered. Finally, the effect of sloshing on the amount of heat exchange between liquid and gas hydrogen was discussed. Considering the heat ingress into liquid hydrogen according to the presence/absence of a sloshing excitation, the amount of heat flux and BOG were discussed for each filling ratio.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.499-504
• /
• 2001

A validation of a 3D CFD model for predicting local subcooling of moderator in the vicinity of calandria tubes in a CANDU reactor is performed. The small scale moderator experiments performed at Sheridan Park Experimental Laboratory(SPEL) in Ontario, Canada[1] is used for the validation. Also a comparison is made between previous CFD analyses based on 2DMOTH and PHOENICS, and the current model analysis for the same SPEL experiment. For the current model, a set of grid structures for the same geometry as the experimental test section is generated and the momentum, heat and continuity equations are solved by CFX-4.3, a CFD code developed by AEA technology. The matrix of calandria tubes is simplified by the porous media approach. The standard $k-\varepsilon$ turbulence model associated with logarithmic wall treatment and SIMPLEC algorithm on the body fitted grid are used and buoyancy effects are accounted for by the Boussinesq approximation. For the test conditions simulated in this study, the flow pattern identified is a buoyancy-dominated flow, which is generated by the interaction between the dominant buoyancy force by heating and inertial momentum forces by the inlet jets. As a result, the current CFD moderator analysis model predicts the moderator temperature reasonably, and the maximum error against the experimental data is kept at less than $2.0^{\circ}C$ over the whole domain. The simulated velocity field matches with the visualization of SPEL experiments quite well.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.3
• /
• pp.232-241
• /
• 2011

The hybrid rocket shows interesting characteristics of complicated mixing layer developed by the interaction between turbulent oxidizer flow and injected surface mass flow from fuel vaporization. In this study, the compressible LES was conducted to explore the physical phenomena of surface oscillatory flow induced by the flow interferences in a duct domain. From the numerical results, the wall injection generates the stronger streamwise vorticites and the negative components of axial velocity accompanied with the azimuthal vorticity near the surface. And the vortex shedding with a certain time scale was found to be developed by hydrodynamic instability in the mixing layer. The pressure fluctuations in this calculation exhibit a peculiar peak at a specific angular frequency($\omega$=8.8) representing intrinsic oscillation due to the injection.

• Journal of Energy Engineering
• /
• v.8 no.2
• /
• pp.224-232
• /
• 1999

Latent heat storage system using micro-encapsuled phase change material is effective method for floor heating of house and building. The temperature profile in capsule block and flow rate of hot water are important parameters for the development of heat storage system. In the present study, a mathematical model based on 3-D, non-steady state, Navier-Stokes equations, scalar conservation equations and turbulence model ($\kappa$-$\varepsilon$), is used to predict the temperature profiles in capsule and the velocity vectors in hot water pipe. The multi-block grids and fine grids embedding are used to join the circle in hot water pipe and square in capsule block. The phase change process of the capsule is quite complex not only because the size of phase change material is very small, but also because phase change material is mixed with the cement to form thermal storage block. In calculation, it's assumed that the phenomena of phase change is limited only the thermal properties of phase change material and the change of boundary is not happened in capsule. The purpose of this study is to calculate the temperature profiles in capsule block and velocity vectors in hot water pipe using the numerical calculation. Two kinds of thermal boundary condition were considered, the first (case 1) is the adiabatic condition for the both outside surfaces of the wall, the second (case 2) is the case in which one surface is natural convection with atmosphere and another surface is adaibatic. Calculation results are shown that the temperature profile in capsule block for case 1 is higher than that for case 2 due to less heat loss in adaibatic surface. Specially, in the domain of near Y=0, the difference of temperature is greater in case 1 than in case 2. The detailed experimental data of capsule block on the temperature profile and the thermal properties such as specific heat and coefficient of heat transfer with the various temperature are required to predict more exact phenomena of heat transfer.