• Environmental Engineering Research
• /
• v.11 no.5
• /
• pp.241-249
• /
• 2006

A mathematical model is described for the prediction of convective upward transport of an organic solvent driven by evaporation at the surface, which is known as the major transport mechanism in the in-situ photolysis of a soil contaminated with 2,3,7,8-tetrachlorodibenzo-p-dioxin(TCDD). A finite-element model was proposed to incorporate the effects of multiphase flow on the distribution of each fluid, gravity as a driving force, and the use of hysteretic models for more accurate description of k-S-p relations. Extensive numerical calculations were performed to study fluid flow through three types of soils under different water table conditions. Predictions of relative permeability-saturation-pressure (k-S-p) relations and fluids distribution for an illustrative soil indicate that hysteresis effects may be quite substantial. This result emphasizes the need to use hysteretic models in performing flow simulations including reversals of flow paths. Results of additional calculations accounting for hysteresis on the one-dimensional unsaturated soil columns show that gravity affects significantly on the flow of each fluid during gravity drainage, solvent injection, and evaporation, especially for highly permeable soils. The rate and duration of solvent injection also have a profound influence on the fluid saturation profile and the amount of evaporated solvent. Key factors influencing water drainage and solvent evaporation in soils also include hydraulic conductivity and water table configuration.

• Journal of Industrial and Engineering Chemistry
• /
• v.68
• /
• pp.274-281
• /
• 2018

Experimental studies of the hydrodynamic performance of the draft tube conical spout-fluid bed (DCSF) were conducted using pharmaceutical pellets. The experiments were carried out in a DCSF consisted of two sections: (a) a conical section with the cross section of $120mm{\times}250mm$ and the height of 270 mm, (b) a cylindrical section with the diameter of 250 mm and the height of 600 mm. The flow characteristics of solids were investigated with a high speed camera and a pezoresistive absolute pressure transducer simultaneously. These characteristics revealed different flow regimes in the DCSF: packed bed at low gas velocities, fluidized bed in draft tube at higher gas velocities until minimum spouting, and spouted bed. The stable spouting was identified by the presence of two dominant frequencies of the power spectrum density of pressure fluctuation signature: (i) the frequency band 6-9 Hz and (ii) the frequency band 12-15 Hz. The pressure drops across the draft tube as well as the annulus measured in order to better recognize the flow structure in the DCSF. It was observed that the pressure drop across the draft tube, the pressure drop across the annulus, and the minimum spouting velocity increase with the increase in the height of draft tube and distance of the entrainment zone, but with the decrease in the distributor hole pitch. Finally, this study provided novel insight into the hydrodynamic of DCSF, particularly minimum spouting and stable spouting in the DCSF which contains valuable information for process design and scale-up of spouted bed equipment.

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

• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.4
• /
• pp.361-372
• /
• 2019

The flow pattern of air layers and skin-friction drag reduction by air injection are investigated to find the suitable multiphase flow model using unstructured finite-volume CFD solver for the Reynolds-averaged Navier-Stokes equations. In the present computations, two different multiphase flow modeling approaches, such as the Volume of Fluid (VOF) and the Eulerian Multi-Phase (EMP), are adopted to investigate their performances in resolving the two-phase flow pattern and in estimating the frictional drag reduction. First of all, the formation pattern of air layers generated by air injection through a circular opening on the bottom of a flat plate are investigated. These results are then compared with those of MMkiharju's experimental results. Subsequently, the quantitative ratios of skin-friction drag reduction including the behavior of air layers, within turbulent boundary layers in large scale and at high Reynolds number conditions, are investigated under the same conditions as the model test that has been conducted in the US Navy's William B. Morgan Large Cavitation Channel (LCC). From these results, it is found that both VOF and EMP models have similar capability and accuracy in capturing the topology of ventilated air cavities so called'air pockets and branches'. However, EMP model is more favorable in predicting quantitatively the percentage of frictional drag reduction by air injection.

• Journal of Korea Society of Waste Management
• /
• v.34 no.7
• /
• pp.685-696
• /
• 2017

This study focuses on computational particle fluid dynamics (CPFD) modeling for the fast pyrolysis of biomass in a conical spouted bed reactor. The CPFD simulation was conducted to understand the hydrodynamics, heat transfer, and biomass fast pyrolysis reaction of the conical spouted bed reactor and the multiphase-particle in cell (MP-PIC) model was used to investigate the fast pyrolysis of biomass in a conical spouted bed reactor. A two-stage semi-global kinetics model was applied to model the fast pyrolysis reaction of biomass and the commercial code (Barracuda) was used in simulations. The temperature of solid particles in a conical spouted bed reactor showed a uniform temperature distribution along the reactor height. The yield of fast pyrolysis products from the simulation was compared with the experimental data; the yield of fast pyrolysis products was 74.1wt.% tar, 17.4wt.% gas, and 8.5wt.% char. The comparison of experimental measurements and model predictions shows the model's accuracy. The CPFD simulation results had great potential to aid the future design and optimization of the fast pyrolysis process for biomass.

• Elastomers and Composites
• /
• v.51 no.4
• /
• pp.263-268
• /
• 2016

In the polymer shaping process that uses molds, the quality of the molded products is determined not only by the flow of the (molten) polymer but also by the air venting in the cavity. Inadequate air venting in the cavity can cause defects in the product, such as voids, short shot, or black streaks. As it is critical to consider the location and size of the vents for proper venting of the air in the cavity, a method that predicts the flow of air and material is required. The venting of air by the flow of rubber inside the cavity was simulated by using a multi-phase computational fluid dynamics method. Through computer simulation, the interface of rubber and air over time was predicted. Then, the velocity and pressure distribution of the venting air were observed. Our research proposes a fundamental method for analyzing the multi-phase flow of polymer materials and air inside the cavity of a mold.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.1
• /
• pp.35-44
• /
• 2017

Supercavitation is one of the most attractive technologies to achieve high speed for underwater vehicles. However, the multiphase flow with high-speed around the supercavitating vehicle (SCV) is difficult to simulate accurately. In this paper, we use modified the turbulent viscosity formula in the Standard K-Epsilon (SKE) turbulent model to simulate the supercavitating flow. The numerical results of flow over several typical cavitators are in agreement with the experimental data and theoretical prediction. In the last part, a flying SCV was studied by unsteady numerical simulation. The selected computation setup corresponds to an outdoor supercavitating experiment. Only very limited experimental data was recorded due to the difficulties under the circumstance of high-speed underwater condition. However, the numerical simulation recovers the whole scenario, the results are qualitatively reasonable by comparing to the experimental observations. The drag reduction capacity of supercavitation is evaluated by comparing with a moving vehicle launching at the same speed but without supercavitation. The results show that the supercavitation reduces the drag of the vehicle dramatically.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.1
• /
• pp.495-509
• /
• 2019

The objective of this study is to analyze the compressibility effects of multiphase cavitating flow during the water-entry process. For this purpose, the water-entry of a projectile at transonic speed is investigated computationally. A temperature-adjusted Tait equation is used to describe the compressibility effects in water, and air and vapor are treated as ideal gases. First, the computational methodology is validated by comparing the simulation results with the experimental measurements of drag coefficient and the theoretical results of cavity shape. Second, based on the computational methodology, the hydrodynamic characteristics of flow are investigated. After analyzing the cavitating flow in compressible and incompressible fluids, the characteristics under compressible conditions are focused upon. The results show that the compressibility effects play a significant role in the development of cavitation and the pressure inside the cavity. More specifically, the drag coefficient and cavity size tend to be larger in the compressible case than those in the incompressible case. Furthermore, the influence of entry velocities on the hydrodynamic characteristics is investigated to provide an insight into the compressibility effects on cavitating flow. The results show that the drag coefficient and the impact pressure vary with the entry velocity, and the prediction formulas for drag coefficient and impact pressure are established respectively in the present study.

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

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