• Carbon letters
• /
• v.25
• /
• pp.25-32
• /
• 2018

In this study, a thermal-gradient chemical vapor infiltration (TG-CVI) process was numerically studied in order to enhance the deposition uniformity within the preform. The computational fluid dynamics technique was used to solve the governing equations for heat transfer and gas flow during the TG-CVI process for two- and three-dimensional (2-D and 3-D) models. The temperature profiles in the 2-D and 3-D models showed good agreement with each other and with the experimental results. The densification process was investigated in a 2-D axisymmetric model. Computation results showed the distribution of the SiC deposition rate within the preform. The results also showed that using two-zone heater gave better deposition uniformity.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2005.05b
• /
• pp.21-24
• /
• 2005

A computational analysis of hygro-thermal and mechanical behaviour of concrete column at high temperature is presented. The objective of this study is to develop a finite difference model that simulates coupled heat and transport phenomena in reinforced concrete structures exposed to rapid heating conditions such as fires. The theoretical basis for the integrated finite difference method is presented to describe a powerful numerical technique for solving of fluid flow in porous media. The numerical results predict the phenomena of 'moisture clog' and the explosive spalling of concrete under fire. The investigations show that high-strength concrete(HSC) and normal-strength concrete(NSC) exposed to high temperature have different pore pressure buildup dependent on porosity, permeability and moisture contents. HSC has more possibility than NSC on spalling.

• The Journal of Engineering Geology
• /
• v.23 no.4
• /
• pp.447-455
• /
• 2013

A micro-mechanical pipe network model with the shape of a cube was developed to simulate the behavior of fluid flow through a porous medium. The fluid-flow mechanism through the cubic pipe network channels was defined mainly by introducing a well-known percolation theory (Stauffer and Aharony, 1994). A non-uniform flow generally appeared because all of the pipe diameters were allocated individually in a stochastic manner based on a given pore-size distribution curve and porosity. Fluid was supplied to one surface of the pipe network under a certain driving pressure head and allowed to percolate through the pipe networks. A percolation condition defined by capillary pressure with respect to each pipe diameter was applied first to all of the network pipes. That is, depending on pipe diameter, the fluid may or may not penetrate a specific pipe. Once pore pressures had reached equilibrium and steady-state flow had been attained throughout the network system, Darcy's law was used to compute the resultant permeability. This study investigated the sensitivity of network size to permeability calculations in order to find out the optimum network size which would be used for all the network modelling in this study. Mean pore size and pore size distribution curve obtained from field are used to define each of pipe sizes as being representative of actual oil sites. The calculated and measured permeabilities are in good agreement.

• Nuclear Engineering and Technology
• /
• v.20 no.2
• /
• pp.132-136
• /
• 1988

A strategic question in groundwater contaminant transport modeling is whether we need to treat waste packages or drums as individual, discrete sources or as approximately lumped sources. In this paper we present analyses of array sources in porous media. We analyze a planar array of sources in porous media with groundwater flow. We compare the concentration field predicted by a detailed model of individual point sources to concentration fields predicted by an infinite plane source and a single point source, all of the same equivalent strength. From this study we identified three regions: (1) a region close to the sources where the effects of adjacent sources are significant and individual source models should be used, (2) a region extending from a few meters to hundreds to thousands of meters downstream, where an equivalent source of infinite extent gives accurate results, and (3) a far-field region, where in an equivalent source of finite extent gives accurate results.

• Nuclear Engineering and Technology
• /
• v.50 no.2
• /
• pp.297-305
• /
• 2018

Fuel assembly (FA) bow in pressurized water reactor (PWR) cores is considered to be a complex process with a large number of influencing mechanisms and several unknowns. Uncertainty and sensitivity analyses are a common way to assess the predictability of such complex phenomena. To perform such analyses, a structural model of a row of 15 FAs in the reactor core is implemented with the finite-element code ANSYS Mechanical APDL. The distribution of lateral hydraulic forces within the core row is estimated based on a two-dimensional Computational Fluid Dynamics model with porous media, assuming symmetric or asymmetric core inlet and outlet flow profiles. The influence of the creep rate on the bow amplitude is tested based on different creep models for guide tubes and fuel rods. Different FA initial states are considered: fresh FAs or FAs with higher burnup, which may be initially straight or exhibit an initial bow from previous cycles. The simulation results over one reactor cycle demonstrate that changes in the creep rate and the hydraulic conditions may have a considerable impact on the bow amplitudes and the bow patterns. A good knowledge of the specific creep behavior and the hydraulic conditions is therefore crucial for making reliable predictions.

• Nuclear Engineering and Technology
• /
• v.36 no.6
• /
• pp.559-570
• /
• 2004

A computational fluid dynamics model for predicting moderator circulation inside the Canada deuterium uranium (CANDU) reactor vessel has been developed to estimate the local subcooling of the moderator in the vicinity of the calandria tubes. The buoyancy effect induced by the internal heating is accounted for by the Boussinesq approximation. The standard $k-{\varepsilon}$ turbulence model with logarithmic wall treatment is applied to predict the turbulent jet flows from the inlet nozzles. The matrix of the calandria tubes in the core region is simplified to a porous media in which the anisotropic hydraulic impedance is modeled using an empirical correlation of pressure loss. The governing equations are solved by DFX-4.4, a commercial CFD code developed by AEA technology. The resultant flow patterns of the constant-z slices containing the inlet nozzles and the outlet port are "mined-type", as observed in the former 2-dimensional experimental investigations. With 103% full power for conservatism, the maximum temperature of the moderator is $82.9^{\circ}C$ at the top of the core region. Considering the hydrostatic pressure change, the minimum subcooling is $24.8^{\circ}C$.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.9
• /
• pp.923-932
• /
• 2011

Because of the instability of a flow pattern in the inlet transition square duct (hereinafter referred to as "transition duct") of a heat recovery steam generator (hereinafter referred to as "HRSG") in a combined cycle power plant, the Reynolds number in the first row of a tube bank is differs sharply from that in the sectional area of the transition duct. This causes differences in the heat flux in each tube in the tube bank. The computational fluid dynamics (CFD) predictions provide three-dimensional results for velocity, temperature, and other flow parameters over the entire domain of the duct and HRSG. A renormalization group theory (RNG) based k-${\epsilon}$�� turbulent model is used for obtaining the results cited in this study. A porous media option is used for modeling the tube banks and the number of transfer units method is used for determining the heat transfer characteristics. This study describes a comparison between the numerical simulation results and actual design output.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.8
• /
• pp.65-73
• /
• 2013

Ground-coupled heat pump system has attracted attention as a promising renewable energy technology due to its improving energy efficiency and eco-friendly mechanism for space cooling and heating. Pipes buried in the ground play a role of direct thermal interaction between circulating fluid inside the pipe and surrounding soils in the geothermal exchange system. However, both complexities of turbulent flow coupling thermal-hydraulic phenomena and very long aspect ratio of the pipe make it difficult to model the heat exchange system directly. Energy balance for fluid flow inside the pipe was derived to model thermal-hydraulic phenomena, and one-dimensional pipe element was proposed through Galerkin formation and time integration of the equation. Developed element is combined to pre-developed FEM code for THM phenomena in porous media. Numerical results of Thermal Response Test showed that line-source model overestimates equivalent thermal conductivity of surrounding soils due to thermal interaction between adjacent pipes and finite length of the pipe. Thus, inverse analysis for the TRT simulation was conducted to present optimal transformation matrix with utmost convergence.

• Journal of Soil and Groundwater Environment
• /
• v.6 no.2
• /
• pp.69-81
• /
• 2001

The constitutive relation among capillary pressure, saturation and relative permeability should be predetermined in order to simulate immiscible water-gas flow in porous media. The relation between saturation and relative permeability becomes more important when the capillary force can be disregarded and viscous friction force governs the flow. In this study, a 2-dimensional finite difference numerical model was developed, in which the variation of viscosity with pressure and that of relative permeability with water saturation can be treated. Seven cases of parallel plate tests were performed in order to obtain the characteristic equation of relative permeability which would be used in. the developed numerical model. It was not possible, however, to match the curves of relative permeability from the plate tests with the existing emperical models. Consequently a logistic equation was proposed as a new emperical model. As this model was composed of the parameter involving aperture size, any aperture size of fracture can be applied to the model. For the purpose of verification, the characteristic equation of relative permeability was applied to the developed numerical model and the computed results were compared with those of plate test. As a result of application of numerical model, in order to check the field applicability, to single fracture surrounding an underground storage cavern, the simultaneous flow of water and propane gas was able to be simulated properly by the model.

• Journal of Soil and Groundwater Environment
• /
• v.9 no.1
• /
• pp.1-11
• /
• 2004

Foam reduces the mobility of gas phase in porous media to overcome gravity override and to divert acid into desired layers in the petroleum industry and to enhance the efficiency of environmental remediation. Recent experimental studies on foam show that foam exhibits a remarkably different flow rheology depending on the flow regime. This study, for the first time, focuses on the issues of foam diversion process under the conditions relevant to groundwater remediation, combining results from laboratory linear-flow experiments and a simple numerical model with permeability contrasts. Linear flow tests performed at two different permeabilities (k = 9.1 and 30.4 darcy) confirmed that two flow regimes of steady-state strong foams were also observed within the permeability range of shallow geological formations. Foam exhibited a shear-thinning behavior in a low-quality regime and near Newtonian rheology in a high-quality regime. Data taken from linear flow tests were incorporated into a simple numerical model to evaluate the efficiency of foam diversion process in the presence of permeability contrasts. The simple model illustrated that foam in the high-quality regime exhibited a successful diversion but foam in the low-quality regime resulted in anti-diversion, implying that only foam in the high-quality regime would be applicable to the diversion process. Sensitivity study proved that the success of diversion process using foam in the high-quality regime was primarily controlled by the limiting capillary pressures (${P_c}{^*}$) of the two layers of interest. Limitations and implications are also discussed and included.