• Journal of computational fluids engineering
• /
• v.22 no.1
• /
• pp.67-80
• /
• 2017

The high-speed bypass effect on the heat exchanger performance has been investigated numerically. The plate-fin type heat exchanger was modeled using two-dimensional porous approximation for the fin region. Governing equations of mass, momentum, and energy equations for compressible turbulent flow were solved using ideal-gas assumption for the air flow. Various bypass-channel height were considered for Mach numbers ranging 0.25-0.65. Due to the existence of the fin in the bypass channel, the main flow tends to turn into the core region of the channel, which results in the distorted velocity profile downstream of the fin region. The boundary layer thickness, displacement thickness, and the momentum thickness showed the variation of mass flow through the fin region. The mass flow variation along the fin region was also shown for various bypass heights and Mach numbers. The volumetric entropy generation was used to assess the loss mechanism inside the bypass duct and the fin region. Finally, the correlations of the friction factor and the Colburn j-factor are summarized.

• Korean Chemical Engineering Research
• /
• v.59 no.1
• /
• pp.138-151
• /
• 2021

The effect of a reactant ratio on the growth of a buoyancy-driven instability in an irreversible A+B→C reaction system is analyzed theoretically and numerically. Taking a non-stoichiometric reactant ratio into account, new linear stability equations are derived without the quasi-steady state assumption (QSSA) and solved analytically. It is found that the main parameters to explain the present system are the Damköhler number, the dimensionless density difference of chemical species and the ratio of reactants. The present initial grow rate analysis without QSSA shows that the system is initially unconditionally stable regardless of the parameter values; however, the previous initial growth rate analysis based on the QSSA predicted the system is unstable if the system is physically unstable. For time evolving cases, the present growth rates obtained from the spectral analysis and pseudo-spectral method support each other, but quite differently from that obtained under the conventional QSSA. Adopting the result of the linear stability analysis as an initial condition, fully nonlinear direct numerical simulations are conducted. Both the linear analysis and the nonlinear simulation show that the reactant ratio plays an important role in the onset and the growth of the instability motion.

• Structural Engineering and Mechanics
• /
• v.30 no.6
• /
• pp.713-732
• /
• 2008

Serviceability and durability of the concrete members can be seriously affected by the corrosion of steel rebar. Carbonation front and or chloride ingress can destroy the passive film on rebar and may set the corrosion (oxidation process). Depending on the level of oxidation (expansive corrosion products/rust) damage to the cover concrete takes place in the form of expansion, cracking and spalling or delamination. This makes the concrete unable to develop forces through bond and also become unprotected against further degradation from corrosion; and thus marks the end of service life for corrosion-affected structures. This paper presents an analytical model that predicts the weight loss of steel rebar and the corresponding time from onset of corrosion for the known corrosion rate and thus can be used for the determination of time to cover cracking in corrosion affected RC member. This model uses fully the thick-walled cylinder approach. The gradual crack propagation in radial directions (from inside) is considered when the circumferential tensile stresses at the inner surface of intact concrete have reached the tensile strength of concrete. The analysis is done separately with and without considering the stiffness of reinforcing steel and rust combine along with the assumption of zero residual strength of cracked concrete. The model accounts for the time required for corrosion products to fill a porous zone before they start inducing expansive pressure on the concrete surrounding the steel rebar. The capability of the model to produce the experimental trends is demonstrated by comparing the model's predictions with the results of experimental data published in the literature. The effect of considering the corroded reinforcing steel bar stiffness is demonstrated. A sensitivity analysis has also been carried out to show the influence of the various parameters. It has been found that material properties and their inter-relations significantly influence weight loss of rebar. Time to cover cracking from onset of corrosion for the same weight loss is influenced by corrosion rate and state of oxidation of corrosion product formed. Time to cover cracking from onset of corrosion is useful in making certain decisions pertaining to inspection, repair, rehabilitation, replacement and demolition of RC member/structure in corrosive environment.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.11
• /
• pp.905-911
• /
• 2015

In this study, to examine the effect of a hole size change(smaller hole diameter) in the outer region of the lower-support-structure-bottom plate(LSSBP) on the reactor core-inlet flow-distribution, simulations were conducted with the commercial CFD software, ANSYS CFX R.15. The predicted results were compared with those of the original LSSBP. Through these comparisons, it was concluded that a more uniform distribution of the mass flow rate at the core-inlet plane could be obtained by reducing the hole size in the outer region of the LSSBP. Therefore, from the nuclear regulatory perspective, design change of the hole pattern in the outer region of the LSSBP may be desirable in terms of improving both the mechanical integrity of the fuel assembly and the core thermal margin.

• Journal of the Korea Organic Resources Recycling Association
• /
• v.17 no.1
• /
• pp.73-79
• /
• 2009

The hydraulic conductivity of porous media is the most important property in soil characteristics. The hydraulic conductivity is determined by outdoor and indoor methods. Indoor methods normally use soil columns for flow test. Assumption of the column test is that fluid one-dimensionally flows through the column. However, fluids may flow toward the wall of the column, resulting in "boundary flow". This study investigated the effect of boundary flow on the hydraulic conductivity by using a permeameter excluding boundary flow. The results showed that the hydraulic conductivity excluding boundary flow was much smaller than the hydraulic conductivity employing the conventional determination method. This study also investigated the effects of particle size and surfactant on the hydraulic conductivity. As the particle size increased, the hydraulic conductivity was increased. The hydraulic conductivity was reduced by increasing surfactant concentration. The result showed that the viscosity of fluid significantly affected the determination of hydraulic conductivity.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.14 no.1
• /
• pp.19-33
• /
• 2002

A numerical method to efficiently secure necessary design informations of the hydraulic characteristics of rubble mound breakwater was attempted here. The method combines the exterior wave field with the interior wave field which is formulated incorporating porous media flow inside the breakwaters. An approximate method based on the long wave assumption was used for the exterior wave field while a boundary element method was used for the interior wave field. A hydraulic experiment was also performed to verify the validity of the numerical analysis. The numerical results were compared with experimental data and results from existing formulae. They generally agreed in both reflection and transmission coefficients. The calculated pore pressures also showed a similar pattern with experimental data, even if they gave some significant differences in their values fur some cases. The main cause of such differences can be attributed to the strongly nonlinear wave field occurring on the frontal slope of the breakwater. The direct input of dynamic pressures(measured from hydraulic experiment) into the numerical method was suggested as a promising method to enhance the predictability of pore pressures.

• Economic and Environmental Geology
• /
• v.28 no.5
• /
• pp.503-512
• /
• 1995

Groundwater flow in fracture networks is simulated using a discrete fracture (DF) model which assume that groundwater flows only through the fracture network. This assumption is available if the permeability of rock matrix is very low. It is almost impossible to describe fracture networks perfectly, so a stochastic approach is used. The stochastic approach assumes that the characteristic parameters in fracture network have special distribution patterns. The stochastic model generates fracture networks with some characteristic parameters. The finite element method is used to compute fracture flows. One-dimensional line element is the element type of the finite elements. The simulation results are shown by dominant flow paths in the fracture network. The dominant flow path can be found from the simulated groundwater flow field. The model developed in this study provides the tool to estimate the influences of characteristic parameters on groundwater flow in fracture networks. The influences of some characteristic parameters on the frcture flow are estimated by the Monte Carlo simulation based on 30 realizations.

• Nuclear Engineering and Technology
• /
• v.53 no.5
• /
• pp.1676-1685
• /
• 2021

We develop a new high-fidelity multiphysics model to simulate boron chemistry in the porous Chalk River Unidentified Deposit (CRUD) deposits. Heat transfer, capillary flow, solute transport, and chemical reactions are fully coupled. The evaporation of coolant in the deposits is included in governing equations modified by the volume-averaged assumption of wick boiling. The axial offset anomaly (AOA) of the Seabrook nuclear power plant is simulated. The new model reasonably predicts the distributions of temperature, pressure, velocity, volumetric boiling heat density, and chemical concentrations. In the thicker CRUD regions, 60% of the total heat is removed by evaporative heat transfer, causing boron species accumulation. The new model successfully shows the quantitative effect of coolant evaporation on the local distributions of boron. The total amount of boron in the CRUD layer increases by a factor of 1.21 when an evaporation-driven increase of soluble and precipitated boron concentrations is reflected. In addition, the concentrations of B(OH)₃ and LiBO₂ are estimated according to various conditions such as different CRUD thickness and porosity. At the end of the cycle in the AOA case, the total mass of boron incorporated in CRUD deposits of a reference single fuel rod is estimated to be about 0.5 mg.