• Journal of Korea Water Resources Association
• /
• v.54 no.11
• /
• pp.943-953
• /
• 2021

The propagation of impact wave induced by landslide and debris flow occurred on the slope of lake, reservoir and bays is a three-dimensional natural phenomenon associated with strong interaction of debris flow and water flow in complex geometrical environments. We carried out 3D numerical modeling of such impact wave in a bay using a multiphase turbulence flow model and a rheology model for non-Newtonian debris flow. Numerical results are compared with previous experimental result to evaluate the performance of present numerical approach. The results underscore that the reasonable predictions of both thickness and speed of debris flow head penetrating below the water surface are crucial to accurately reproduce the maximum peak height and free surface profiles of impact wave. Two predictions computed using different initial debris flow thicknesses become different from the instant when the peaks of impact waves fall due to the gravity. Numerical modeling using relatively thick initial debris flow thickness appears to well reproduce the water surface profile of impact wave propagating across the bay as well as wave run-up on the opposite slope. The results show that the maximum run-up height on the opposite slope is not sensitive to the initial thickness of debris flows of same total volume. Meanwhile, appropriate rheology model for debris flow consisting of inviscid particle only should be employed to more accurately reproduce the debris flow propagating along the channel bottom.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.7
• /
• pp.723-733
• /
• 2011

A numerical method of the CSF(Continuum Surface Force) model is presented for the calculation of the surface tension force and implemented in an in-house solution code(PowerCFD). The present method(code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with volume capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The application of the present method to a 2-D liquid drop problem is illustrated by an equilibrium and nonequilibrium oscillating drop calculation. It is found that the present method simulates efficiently and accurately surface tension-dominant multiphase flows.

• Nuclear Engineering and Technology
• /
• v.53 no.5
• /
• pp.1540-1555
• /
• 2021

The plate-type fuel assembly adopted in nuclear research reactor suffers from complicated effect induced by non-uniform irradiation, which might affect its stress conditions, mechanical behavior and thermal-hydraulic performance. A reliable numerical method is of great importance to reveal the complex evolution of mechanical deformation, flow redistribution and temperature field for the plate-type fuel assembly under non-uniform irradiation. This paper is the first part of a two-part study developing the numerical methodology for the thermal-fluid-structure coupling behaviors of plate-type fuel assembly under irradiation. In this paper, the thermal-fluid-structure coupling methodology has been developed for plate-type fuel assembly under non-uniform irradiation condition by exchanging thermal-hydraulic and mechanical deformation parameters between Finite Element Model (FEM) software and Computational Fluid Dynamic (CFD) software with Mesh-based parallel Code Coupling Interface (MpCCI), which has been validated with experimental results. Based on the established methodology, the effects of non-uniform irradiation and fluid were discussed, which demonstrated that the maximum mechanical deformation with irradiation was dozens of times larger than that without irradiation and the hydraulic load on fuel plates due to differential pressure played a dominant role in the mechanical deformation.

• Nuclear Engineering and Technology
• /
• v.55 no.8
• /
• pp.3017-3029
• /
• 2023

In the late in-vessel phase of a nuclear reactor severe accident, the internal heat transfer and crust evolution during the debris bed melting process have important effects on the thermal load distribution along the vessel wall, and further affect the reactor pressure vessel (RPV) failure mode and the state of melt during leakage. This study coupled the phase change model and large eddy simulation to investigate the variations of the temperature, melt liquid fraction, crust and heat flux distributions during the debris bed melting process in the hypothetical severe accident of HPR1000. The results indicated that the heat flow towards the vessel wall and upper surface were similar at the beginning stage of debris melting, but the upward heat flow increased significantly as the development of the molten pool. The maximum heat flux towards the vessel wall reached 0.4 MW/m². The thickness of lower crust decreased as the debris melting. It was much thicker at the bottom region with the azimuthal angle below 20° and decreased rapidly at the azimuthal angle around 20-50°. The maximum and minimum thicknesses were 2 and 90 mm, respectively. By contrast, the distribution of upper crust was uniform and reached stable state much earlier than the lower crust, with the thickness of about 10 mm. Moreover, the sensitivity analysis of initial condition indicated that as the decrease of time interval from reactor scram to debris bed dried-out, the maximum debris temperature and melt fraction became larger, the lower crust thickness became thinner, but the upper crust had no significant change. The sensitivity analysis of in-vessel retention (IVR) strategies indicated that the passive and active external reactor vessel cooling (ERVC) had little effect on the internal heat transfer and crust evolution. In the case not considering the internal reactor vessel cooling (IRVC), the upper crust was not obvious.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.19 no.2
• /
• pp.55-64
• /
• 2015

In this study, a 3-D flame cooling analysis is conducted to examine thermal safety for the flame deflector of the 75 tonf class propulsion test facility, and the safe discharge of the exhaust gas is assessed by using numerical results. The Mixture multiphase model is adopted for the simulation of heat transfer and phase exchange process between flame and cooling water, and the computational study using the single species unreacted model for the exhaust plume is carried out for the flame cooling. Numerical analysis predicts maximum temperature on the flame deflector wall for different water flow rates, and evaluates the safe minimum flow rate of water corresponding to the fire-resistant temperature for concrete.

• Journal of Power Electronics
• /
• v.16 no.2
• /
• pp.553-563
• /
• 2016

This paper presents a digital hardware implementation of a real-time simulator for a multiphase drive using a field-programmable gate array (FPGA) device. The simulator was developed with a modular and hierarchical design using very high-speed integrated circuit hardware description language (VHDL). Hence, this simulator is flexible and portable. A state-space representation model suitable for FPGA implementations was proposed for a dual three-phase induction machine (DTPIM). The simulator also models a two-level 12-pulse insulated-gate bipolar transistor (IGBT)-based voltage-source converter (VSC), a pulse-width modulation scheme, and a measurement system. Real-time simulation outputs (stator currents and rotor speed) were validated under steady-state and transient conditions using as reference an experimental test bench based on a DTPIM with 15 kW-rated power. The accuracy of the proposed digital hardware implementation was evaluated according to the simulation and experimental results. Finally, statistical performance parameters were provided to analyze the efficiency of the proposed DTPIM hardware implementation method.

• Advances in nano research
• /
• v.6 no.4
• /
• pp.299-321
• /
• 2018

A layerwise shear deformation theory is applied in this paper for buckling analysis of piezoelectric truncated conical shell. The core is a multiphase nanocomposite reinforced by carbon nanotubes (CNTs) and carbon fibers. The top and bottom face sheets are piezoelectric subjected to 3D electric field and external voltage. The Halpin-Tsai model is used for obtaining the effective moisture and temperature dependent material properties of the core. The proposed layerwise theory is based on Mindlin's first-order shear deformation theory in each layer and results for a laminated truncated conical shell with three layers considering the continuity boundary condition. Applying energy method, the coupled motion equations are derived and analyzed using differential quadrature method (DQM) for different boundary conditions. The influences of some parameters such as boundary conditions, CNTs weight percent, cone semi vertex angle, geometrical parameters, moisture and temperature changes and external voltage are investigated on the buckling load of the smart structure. The results show that enhancing the CNTs weight percent, the buckling load increases. Furthermore, increasing the moisture and temperature changes decreases the buckling load.

• Geomechanics and Engineering
• /
• v.13 no.1
• /
• pp.1-23
• /
• 2017

A coupled liquid-gas-solid three-phase model, linking two numerical codes (TOUGH2/EOS3 and $FLAC^{3D}$), was firstly established and validated by simulating an in-situ air flow test in Essen. Then the coupled model was employed to investigate responses of multiphase flow and soil skeleton deformation to compressed air or freshwater injection using the same simulation conditions in an aquifer of Tianjin, China. The simulation results show that with injecting pressurized fluids, the vertical effective stress in some area decreases owing to the pore pressure increasing, an expansion of soil skeleton appears, and land uplift occurs due to support actions from lower deformed soils. After fluids injection stops, soil deformation decreases overall due to injecting fluids dissipating. With the same applied pressure, changes in multiphase flow and geo-mechanical deformation caused by compressed air injection are relatively greater than those by freshwater injection. Furthermore, the expansion of soil skeleton induced by compressed air injection transfers upward and laterally continuously with time, while during and after freshwater injection, this expansion reaches rapidly a quasi-steady state. These differences induced by two fluids injection are mainly because air could spread upward and laterally easily for its lower density and phase state transition appears for compressed air injection.

• Coupled systems mechanics
• /
• v.8 no.4
• /
• pp.339-350
• /
• 2019

A numerical simulation of the incompressible multiphase hydraulic jump flow was performed to compare the interface prediction through the use of the three RANS turbulence models: $k-{\varepsilon}$, $RNGk-{\varepsilon}$ and SST $k-{\omega}$. A three dimensional no submerged hydraulic jump and a two dimensional submerged hydraulic jump were modeled. Both the geometry and the mesh were created using the open source Gmsh code. The project's geometry consists of a rectangular channel with length and height differences between the two dimensional and three dimensional simulations. Uniform hexahedral cells were used for the mesh. Three refining meshes were constructed to allow to verify simulation convergence. The Volume of Fluid (abbr. VOF) method was used for treatment of the air-water surface. The turbulence models were evaluated in three distinct set up configurations to provide a greater accuracy in the flow representation. In the two-dimensional analysis of a submerged hydraulic jump simulation, the turbulence model RNG RNG $k-{\varepsilon}$ provided a better interface adjust with the experimental results than the model $k-{\varepsilon}$ and SST $k-{\omega}$. In the three-dimensional simulation of a no-submerged hydraulic jump the k-# showed better results than the SST $k-{\omega}$ and RNG $k-{\varepsilon}$ capturing the height and length of the ledge with a better fit with the experimental results.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.1
• /
• pp.1-10
• /
• 2020

The application of multiphase flows is increasingly being applied to analyze phenomena such as single phase flows where the fluid boundary changes continuously over time or the problem of mixing a liquid phase and a gas phase. In particular, multiphase flow models that take into account incompressible Newtonian fluids for liquid and gas are often applied to solve the problems of the free water surface such as wave fields. In general, multi-phase flow models require time-based the surface tracking of each fluid's phase boundary, which determines the accuracy of the final calculation of the model. This study evaluates the advection performance of representative VOF-type boundary tracking techniques applied to various CFD numerical codes. The effectiveness of the FCT method to control the numerical flux to minimize the numerical diffusion in the conventional VOF-type boundary tracking method and advection calculation was mainly evaluated. In addition, the possibility of tracking performance of free surface using CIP method (Yabe and Aoki, 1991) was also investigated. Numerical results show that the FCT-VOF method introducing an anti-diffusive flux to precent excessive diffusion is superior to other methods under the confined conditions in this study. The results from this study are expected to be used as an important basic data in selecting free surface tracking techniques applied to various numerical codes.