• Korea-Australia Rheology Journal
• /
• v.16 no.4
• /
• pp.183-191
• /
• 2004

High Deborah or Weissenberg number problems in viscoelastic flow modeling have been known formidably difficult even in the inertialess limit. There exists almost no result that shows satisfactory accuracy and proper mesh convergence at the same time. However recently, quite a breakthrough seems to have been made in this field of computational rheology. So called matrix-logarithm (here we name it tensor-logarithm) formulation of the viscoelastic constitutive equations originally written in terms of the conformation tensor has been suggested by Fattal and Kupferman (2004) and its finite element implementation has been first presented by Hulsen (2004). Both the works have reported almost unbounded convergence limit in solving two benchmark problems. This new formulation incorporates proper polynomial interpolations of the log­arithm for the variables that exhibit steep exponential dependence near stagnation points, and it also strictly preserves the positive definiteness of the conformation tensor. In this study, we present an alternative pro­cedure for deriving the tensor-logarithmic representation of the differential constitutive equations and pro­vide a numerical example with the Leonov model in 4:1 planar contraction flows. Dramatic improvement of the computational algorithm with stable convergence has been demonstrated and it seems that there exists appropriate mesh convergence even though this conclusion requires further study. It is thought that this new formalism will work only for a few differential constitutive equations proven globally stable. Thus the math­ematical stability criteria perhaps play an important role on the choice and development of the suitable con­stitutive equations. In this respect, the Leonov viscoelastic model is quite feasible and becomes more essential since it has been proven globally stable and it offers the simplest form in the tensor-logarithmic formulation.

• Nuclear Engineering and Technology
• /
• v.52 no.8
• /
• pp.1862-1870
• /
• 2020

The purpose of this study is to investigate the effect of welding residual stress on operating stress in designing a nuclear turbine welded rotor. A two-dimensional axisymmetric finite element model is employed to calculate the residual stress before and after post weld heat treatment (PWHT), and then the superposition of residual stress after PWHT and operating stress at normal speed and overspeed were discussed. The investigated results show that operating stress can be affected significantly by welding residual stress, and the distribution trend of superposition stress at the weld area is mainly determined by welding residual stress. The superposition of residual stress and operating stress is linear superposition, and the hoop stress distribution of superposition stress is similar with the distribution of residual stress. With the increasing overspeed, the distribution pattern of the hoop superimposed stress remains almost unchanged, while the stress level increases.

• Applied Science and Convergence Technology
• /
• v.24 no.6
• /
• pp.250-253
• /
• 2015

High contrast grating (HCG) is the structure made up of the sub-wavelength grating of high-index and the surrounding layer of low-index, which reveals high contrast between two materials. Its advantages include high reflectivity over a broad bandwidth, polarization and wavelength selectivity, optical high-Q resonator, and phase modulation. In this work, the HCG structure comprising of indium tin oxide (ITO) and Silicon (Si), for the surrounding layer and the grating layer respectively, was studied. Its theoretical model was established, and transmittance, phase and optical behavior were calculated by rigorous coupled-wave analysis and finite element method. Furthermore, the established structure was fabricated to validate its feasibility. The fabricated structure shows the focusing capability whose length is about $10{\mu}m$, and the feasibility of the structure was demonstrated. It is also meaningful that ITO layer can contribute to the fabrication of the HCG structure, leading to enable the structure to be electrical-driven.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.6
• /
• pp.688-696
• /
• 2004

Nonlinear dynamic behavior of diffusive-thermal instability in diluted CH$_4$/O$_2$ diffusion flames is numerically investigated by adopting detailed chemistry and transport. Counterflow diffusion flame is adopted as a model flamelet. Particular attention is focused on the pulsating-instability regime, which arises for Lewis numbers greater than unity, and the instability occurs at high strain rate near extinction condition in this flame configuration. Once a steady flame structure is obtained for a prescribed value of initial strain rate, transient solution of the flame is calculated after a finite amount of strain-rate perturbation is imposed on the steady flame. Transient evolution of the flame depends on the initial strain rate and the amount of perturbed strain rate. Basically, the dynamic behaviors can be classified into two types, namely non-oscillatory decaying solution and diverging solution leading to extinction. The peculiar oscillatory solution, which has been found in the previous study adopting one-step chemistry and constant Lewis numbers, is net observed in this study, which is attributed to both convective flow and preferential diffusion effects.

• Membrane and Water Treatment
• /
• v.7 no.4
• /
• pp.297-311
• /
• 2016

Modeling and simulation of air humidification by hollow fiber membrane contactors are investigated in the current study. A computational fluid dynamic model was developed by solving the k-epsilon turbulence 2D Navier-Stokes equations as well as mass conservation equations for steady-state conditions in membrane contactors. Finite element method is used for the study of the air humidification under different operating conditions, with a focus on the humidity density, total mass transfer flux and velocity field. There has been good agreement between simulation results and experimental data obtained from literature. It is found that the enhancement of air stream decreases the outlet humidity from 0.392 to 0.340 (module 1) and from 0.467 to 0.337 (module 2). The results also indicated that there has been an increase in air velocity in the narrow space of shell side compared with air velocity wide space of shell side. Also, irregular arrangement has lower dead zones than regular arrangement which leads to higher water flux.

• 한국연소학회:학술대회논문집
• /
• 2004.06a
• /
• pp.95-101
• /
• 2004

Dynamic behavior of diffusive-thermal instability in diluted $CH_4/O_2$ diffusion flames is numerically investigated by adopting detailed chemistry and transport. Counterflow diffusion flame is adopted as a model flamelet. Particular attention is focused on the pulsating-instability regime, which arises for Lewis numbers greater than unity, and the instability occurs at high strain rate near extinction condition in this flame configuration. Once a steady flame structure is obtained for a prescribed value of initial strain rate. transient solution of the flame is calculated after a finite amount of strain-rate perturbation is imposed Oil the steady flame. Transient evolution of the flame depends on the initial strain rate and the amount of perturbed strain rate. Basically, the dynamic behaviors can be classified into two types, namely non-oscillatory decaying solution and diverging solution leading to extinction. The peculiar oscillatory solution. which has been found in the previous study adopting one-step chemistry and constant Lewis numbers, is not observed in this study, which is attributed to both convective flow and preferential diffusion effects.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.201-206
• /
• 2005

In order to realistically represent the complex turbulence-chemistry interaction at the partially premixed turbulent lifted flames encountered in the gas turbine combustors, the combined conserved-scalar/level-set flamelet approach has been adopted. The parallel unstructured-grid finite-volume method has been developed to maintain the geometric flexibility and computational efficiency for the solution of the physically and geometrically complex flows. Special emphasis is given to the swirl effects on the combustion characteristics of the lean-premixed gas turbine combustor. Numerical results suggest that the present approach is capable of realistically simulating the combustion characteristics for the lean-premixed gas turbine engines and the lifted turbulent jet flame with a vitiated coflow.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.10a
• /
• pp.204-207
• /
• 2003

The present paper reports high-fidelity simulation of direct initiation processes of cylindrical detonation waves by concentrated energy deposition. The goal is to understand the underpinning mechanisms in failed or successful detonation initiation processes. We employed the Space-Time CESE method to solve the reacting flow equations, including realistic finite-rate chemistry model of the nine species and twenty-four reactions for H$_2$-O$_2$-Ar mixtures. Detailed results of sub-critical, critical. and supercritical initiation process are reported.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.1
• /
• pp.29-35
• /
• 2001

Radiation-induced oscillatory instability in CH$_4$/Air diffusion flames is numerically investigated by adopting detailed chemistry. Counterflow diffusion flame is employed as a model flamelet and optically thin gas-phase radiation is assumed. Attention is focused on the extinction regime induced by radiative heat loss, which occurs at low strain rate. Once a steady flame structure is obtained for a prescribed value of initial strain rate, transient solution of the flame is calculated after a finite amount of strain-rate perturbation is imposed on the steady flame. Depending on the initial strain rate and the amount of perturbed strain rate, transient evolution of the flame exhibits various types of flame-evolution behaviors. Basically, the dynamic behaviors can be classified into two types, namely oscillatory decaying solution and diverging solution leading to extinction.

• Proceedings of the Korea Concrete Institute Conference
• /
• 1997.04a
• /
• pp.254-260
• /
• 1997

Chloride-induced corrosion of reinforcement is one of the main factors which cause the deterioration of concrete structures. Durability and service lives of the concrete sturctures should be predicted in order to minimize the risk of corrosion of reinforcement. The objective of this study is to suggest the basis of analytical methods of predicting the corrosion threhold time of concrete structures. Based on the chemistry and physics of chloride ion transport and corrosion process, chloride intrusion with various exposure conditions, variability of diffusivity and transport of pore water in concrete are taken into consideration in applying finite element formulation to the predicion of corrosion threhold time. The effects of main factors on the prediction of chloride intrusion and corrosion threhold time are examined. In addition, after chloride diffusivities of several mixture proportions with different parameters are measured by chloride diffusion test, the exemplary anayses of corrosion threhold time of those mixture proportions are carried out.