• Journal of Ocean Engineering and Technology
• /
• v.28 no.6
• /
• pp.517-526
• /
• 2014

The leakage of cryogenic LNG through cracks in the insulation membrane of an LNG carrier causes the hull structure to experience a cold spot as a result of the heat transfer from the LNG. The hull structure will become brittle at this cold spot and the evaporated natural gas may potentially lead to a hazard because of its flammability. This paper presents a computational model for the LNG flow and heat diffusion in an LNG insulation panel subject to leakage. The temperature distribution in the insulation panel and the speed of gas diffusion through it are simulated to assess the safety level of an LNG carrier subject that experiences a leak. The behavior of the leaked LNG is modeled using a multiphase flow that considers the mixture of liquid and gas. The simulation model considers the phase change of the LNG, gas-liquid multiphase interactions in the porous media, and accompanying rates of heat transfer. It is assumed that the NO96-GW membrane storage is composed of glass wool and plywood for the numerical simulation. In the numerical simulation, the seepage, heat diffusion, and evaporation of the LNG are investigated. It is found that the diffusion speed of the leakage is very high to accelerate the evaporation of the LNG.

• Journal of computational fluids engineering
• /
• v.17 no.4
• /
• pp.41-48
• /
• 2012

A numerical analysis of thermal stratification in the upper plenum of the MONJU fast breeder reactor was performed. Calculations were performed for a 1/6 simplified model of the MONJU reactor using the commercial code, CFX-13. To better resolve the geometrically complex upper core structure of the MONJU reactor, the porous media approach was adopted for the simulation. First, a steady state solution was obtained and the transient solutions were then obtained for the turbine trip test conducted in December 1995. The time dependent inlet conditions for the mass flow rate and temperature were provided by JAEA. Good agreement with the experimental data was observed for steady state solution. The numerical solution of the transient analysis shows the formation of thermal stratification within the upper plenum of the reactor vessel during the turbine trip test. The temporal variations of temperature were predicted accurately by the present method in the initial rapid coastdown period (~300 seconds). However, transient numerical solutions show a faster thermal mixing than that observed in the experiment after the initial coastdown period. A nearly homogenization of the temperature field in the upper plenum is predicted after about 900 seconds, which is a much shorter-term thermal stratification than the experimental data indicates. This discrepancy is due to the shortcoming of the turbulence models available in the CFX-13 code for a natural convection flow with thermal stratification.

• Journal of Energy Engineering
• /
• v.21 no.4
• /
• pp.385-392
• /
• 2012

The purpose of this study is investigated on the combustion and the thermal characteristics of porous media burners which are many using for a condensing boiler recently. In addition, results of this study will be used the fundamental information to decide the burner type which will be applied to the future development of EGR(Exhaust gas recirculation) condensing boiler. Two flat type of burners made of a the metal fiber(MF) and the ceramic(CM) were selected and examined, experimentally. As experimental results, the emitted CO concentration of CM was higher than that of MF. However, the NO concentration of MF was higher than that of CM. The efficiencies of both burners were increased as increasing the burner capacity. While the efficiency of MF was higher than that of CM, regardless of the burner capacity. In the experimental range, MF is appropriated for the burner material and 0.8 of equivalence ratio is an optimal operation condition, regarding of the proportional control, the thermal efficiency and emitted NO and CO concentration based on the regulations of KS B standard and EN 677 standard.

• Structural Engineering and Mechanics
• /
• v.66 no.4
• /
• pp.505-513
• /
• 2018

Steam flooding is widely used in heavy oil reservoir with coupling effects among the formation temperature change, fluid flow and solid deformation. The effective stress, porosity and permeability in this process can be affected by the multi-physical coupling of thermal, hydraulic and mechanical processes (THM), resulting in a complex interaction of geomechanical effects and multiphase flow in the porous media. Quantification of the state of deformation and stress in the reservoir is therefore essential for the correct prediction of reservoir efficiency and productivity. This paper presents a coupled fluid flow, thermal and geomechanical model employing a program (MATLAB interface code), which was developed to couple conventional reservoir (ECLIPSE) and geomechanical (ABAQUS) simulators for coupled THM processes in multiphase reservoir modeling. In each simulation cycle, time dependent reservoir pressure and temperature fields obtained from three dimensional compositional reservoir models were transferred into finite element reservoir geomechanical models in ABAQUS as multi-phase flow in deforming reservoirs cannot be performed within ABAQUS and new porosity and permeability are obtained using volumetric strains for the next analysis step. Finally, the proposed approach is illustrated on a complex coupled problem related to steam flooding in an oil reservoir. The reservoir coupled study showed that permeability and porosity increase during the injection scenario and increasing rate around injection wells exceed those of other similar comparable cases. Also, during injection, the uplift occurred very fast just above the injection wells resulting in plastic deformation.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.5
• /
• pp.109-122
• /
• 1993

ln this paper, a finite element analysis procedure is proposed for the incremental multi-step excavations in a fluid-saturated porous medium such as saturated soil ground. As the basis of derivation, Biot's equation was used. The proposed procedure was applied to some one- and two-dimensional problems under incremental excavations. Unsaturated cases as well as saturated cases were considered for comparison. Through numerical tests, the effects of permeability and excavation speed on the deformation history was investigated. Results showed that pore pressure built up during incremental excavation has a significant effect on the deformation and stresses of solid skeleton and validated the use of the present procedure for the analysis of multi-step excavations in fluid-saturated media such as in saturated shallow ground.

• 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 Computational Structural Engineering Institute of Korea
• /
• v.20 no.2
• /
• pp.191-206
• /
• 2007

The general theoretical solutions for the wavespeed and damping derived in Part 1 of this work, are incorporated into the computer code. In this paper the code is used in a parametric study of the influence of excitation frequency and variations in material properties on propagation velocity and damping. Compressional wave velocity for waves of the first kind is shown to vary as a function of the frequency-permeability product, with a zone where wavespeed transitions from a lower bound value to a higher bound value with increasing values of the product. Damping is seen to be a maximum where the rate of change in wavespeed is greatest. Waves of the second kind also show a transition in wavespeed from near zero at low values of the frequency-permeability product to an upper bound value at higher values of the product.

• Journal of Soil and Groundwater Environment
• /
• v.23 no.1
• /
• pp.1-13
• /
• 2018

A series of experiments using transparent micromodels with an artificial pore network etched on glass plates was performed to investigate the effects of flow rate on the migration and distribution of resident wetting porewater (deionized water) and injecting non-wetting fluid (n-hexane). Multicolored images transformed from real RGB images were used to distinguish n-hexane from porewater and pore structure. Hexane flooding followed by immiscible displacement with porewater, migration through capillary fingering, preferential flow and bypassing were observed during injection experiments. The areal displacement efficiency increases as the injection of n-hexane continues until the equilibrium reaches. Experimental results showed that the areal displacement efficiency at equilibrium increases as the flow rate increases. Close observation reveals that preferential flowpaths through larger pore bodies and throats and clusters of entrapped porewater were frequently created at lower flow rate. At higher flow rate, randomly oriented forward and lateral flowpaths of n-hexane displaces more porewater at an efficiency close to stable displacement. It may resulted from that the pore pressure of n-hexane, at higher flow rate, increases fast enough to overcome capillary pressure acting on smaller pore throats as well larger ones. Experimental results in this study may provide fundamental information on migration and distribution of immiscible fluids in subsurface porous media.

• Journal of Hydrogen and New Energy
• /
• v.28 no.6
• /
• pp.610-617
• /
• 2017

This paper presents a study of the effect of thickness of porous Al-Ni electrodes, on the Hydrogen Evolution Reaction (HER) in alkaline media. As varying deposition time at 300 W DC sputtering power, the thickness of the Al-Ni electrodes was controlled from 1 to $20{\mu}m$. The heat treatment was carried out in $610^{\circ}C$, followed by selective leaching of the Al-rich phase. XRD studies confirmed the presence of $Al_3Ni_2$ intermetallic compounds after the heat treatment, indicating the diffusion of Ni from the Ni-rich phase to Al-rich phase. The porous structure of the Al-Ni electrodes after the selective leaching of Al was also confirmed in SEM-EDS analysis. The double layer capacitance ($C_{dl}$) and roughness factor ($R_f$) of the electrodes were increased for the thicker Al-Ni electrodes. As opposed to the general results in above, there were no further improvements of the HER activity in the case of the electrode thickness above $10{\mu}m$. This result may indicate that the $R_f$ is not the primary factor for the HER activity in alkaline media.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.197-215
• /
• 2015

Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front. A one-dimensional finite element formulation is developed for the numerical simulation of variably saturated flow systems. First order backward Euler implicit and second order Crank-Nicolson time discretization schemes are adopted as a solution strategy in this formulation based on Picard and Newton iterative techniques. Five examples are used to investigate the numerical performance of two approaches and the different factors are highlighted that can affect their convergence and efficiency. The first test case deals with sharp moisture front that infiltrates into the soil column. It shows the capability of providing a mass-conservative behavior. Saturated conditions are not developed in the second test case. Involving of dry initial condition and steep wetting front are the main numerical complexity of the third test example. Fourth test case is a rapid infiltration of water from the surface, followed by a period of redistribution of the water due to the dynamic boundary condition. The last one-dimensional test case involves flow into a layered soil with variable initial conditions. The numerical results indicate that the Crank-Nicolson scheme is inefficient compared to fully implicit backward Euler scheme for the layered soil problem but offers same accuracy for the other homogeneous soil cases.