• Journal of Korean Society of Environmental Engineers
• /
• v.30 no.5
• /
• pp.534-543
• /
• 2008

The experimental study with the prototype provides more acceptable data than the others. But there are so many limited conditions to perform the experimental study with the prototype. So the theoretical similitude with the scaled model and the numerical study with the CFD method have been chosen alternatively to analysis the fume movement. In this study, the ventilation was estimated from the results of the numerical study based on the experimental results as the boundary conditions. The grid A and B were same size and shape with the models which was used in the experimental study and consisted with 163,839, 122,965 cells respectively. The height of the fume layer was estimated form the mole fraction of fume components and the ventilation was determined by the velocity and temperature of the fume. The results of this study showed that the fume movements estimated from the numerical study are enough to apply to the prototype if there are proper heat loss correction factors. The numerical study is easier to change study conditions and faster to get results from the study than the experimental study. So if we find some proper heat loss correction factors, it's possible to execute the various and advanced study with the numerical study.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2014.02a
• /
• pp.221-221
• /
• 2014

The present work proposes an improved numerical simulator for design and modification of large area capacitively coupled plasma (CCP) processing chamber. CCP, as notoriously well-known, demands the tremendously huge computational cost for carrying out transient analyses in realistic multi-dimensional models, because electron dissociations take place in a much smaller time scale (${\Delta}t{\approx}10-8{\sim}10-10$) than time scale of those happened between neutrals (${\Delta}t{\approx}10-1{\sim}10-3$), due to the rf drive frequencies of external electric field. And also, for spatial discretization of electron flux (Je), exponential scheme such as Scharfetter-Gummel method needs to be used in order to alleviate the numerical stiffness and resolve exponential change of spatial distribution of electron temperature (Te) and electron number density (Ne) in the vicinity of electrodes. Due to such computational intractability, it is prohibited to simulate CCP deposition in a three-dimension within acceptable calculation runtimes (<24 h). Under the situation where process conditions require thickness non-uniformity below 5%, however, detailed flow features of reactive gases induced from three-dimensional geometric effects such as gas distribution through the perforated plates (showerhead) should be considered. Without considering plasma chemistry, we therefore simulated flow, temperature and species fields in three-dimensional geometry first, and then, based on that data, boundary conditions of two-dimensional plasma discharge model are set. In the particular case of SiH4-NH3-N2-He CCP discharge to produce deposition of SiNxHy thin film, a cylindrical showerhead electrode reactor was studied by numerical modeling of mass, momentum and energy transports for charged particles in an axi-symmetric geometry. By solving transport equations of electron and radicals simultaneously, we observed that the way how source gases are consumed in the non-isothermal flow field and such consequences on active species production were outlined as playing the leading parts in the processes. As an example of application of the model for the prediction of the deposited thickness uniformity in a 300 mm wafer plasma processing chamber, the results were compared with the experimentally measured deposition profiles along the radius of the wafer varying inter-electrode gap. The simulation results were in good agreement with experimental data.

• Journal of the Korean Geosynthetics Society
• /
• v.13 no.2
• /
• pp.59-67
• /
• 2014

The amount of the contaminants that can be adsorbed on the drain was evaluated for the effective remediation of the contaminated soil, and the contaminants adsorptivity of the drain was evaluated by comparing the isothermal adsorption model after carrying out the contaminants adsorption test of the reactants coated on the surface of the drain. The reactant used in the experiment is a natural zeolite, and the contaminants are copper, lead and cadmium. The results that Freundlich and Langmuir adsorption isotherm model are compared to the adsorption amount according to the change of the initial concentration by the contaminants. As a result of the component analysis, because Si, Al and O are contained approximately 28%, 11% and 48%, respectively, it is identified that the material coated on the surface of the drain is the component of the zeolite which is the reactant for the adsorption of the heavy-metal (Cu, Pb, Cd) contaminants. The heavy-metal adsorption kinetic of the zeolite which is the reactant was decreased in order of lead, copper and cadmium. The important factor of the performance evaluation of the adsorbent is the reaction rate, and if zeolite is used as the reactant in the relationship between the maximum amount of adsorption and reaction rate, it can be utilized as the design factor that determine the removal order of the complex heavy-metal. In other words, because the maximum adsorption quantity of lead is smaller compared to copper but the reaction rate is relatively fast, it can be primarily removed, and copper can be removed after removing the lead. It was analyzed that Cadmium can be finally removed after that other heavy-metal is removed.

• Korean Chemical Engineering Research
• /
• v.43 no.2
• /
• pp.249-258
• /
• 2005

A two-column six-step pressure swing adsorption(PSA) process was to study separation of hydrogen from hydrogen and methane($60_{vol}%/40_{vol}%$) binary system onto activated carbon adsorbent. The effects of the feed gas pressure, the feed flowrate and the P/F(purge to feed) ratio on the process performance were evaluated. The cyclic steady-states of PSA process were reached to after 15 cycles. $H_2$ purity increases according as the P/F ratio and pressure increase and the feed flow rate decreases; however, $H_2$ recovery shows an opposite phenomena to the purity. PSA process simulation studied to find optimum operation condition. In the results, 22 LPM feed flowrate, 11 atm adsorption pressure and 0.10 P/F ratio might be optimal values to obtain more than 75% recovery and 99% purity hydrogen. In this study was non-isothermal and non-adiabatic model considering linear driving force(LDF) model and Langmuir-Freundlich adsorption isotherm considered to compare between prediction and experimental data.

• Journal of The Korean Astronomical Society
• /
• v.33 no.2
• /
• pp.111-121
• /
• 2000

Many models of globular cluster formation assume the presence of cold dense clouds in early universe. Here we re-examine the Fall & Rees (1985) model for formation of proto-globular cluster clouds (PGCCs) via thermal instabilities in a protogalactic halo. We first argue, based on the previous study of two-dimensional numerical simulations of thermally unstable clouds in a stratified halo of galaxy clusters by Real et al. (1991), that under the protogalactic environments only nonlinear (${\delta}{\ge}1$) density inhomogeneities can condense into PGCCs without being disrupted by the buoyancy-driven dynamical instabilities. We then carry out numerical simulations of the collapse of overdense douds in one-dimensional spherical geometry, including self-gravity and radiative cooling down to T = $10^4$ K. Since imprinting of Jeans mass at $10^4$ K is essential to this model, here we focus on the cases where external UV background radiation prevents the formation of $H_2$ molecules and so prevent the cloud from cooling below $10^4$ K. The quantitative results from these simulations can be summarized as follows: 1) Perturbations smaller than $M_{min}\~(10^{5.6}\;M{\bigodot})(nh/0.05cm^{-3})^{-2}$ cool isobarically, where nh is the unperturbed halo density, while perturbations larger than $M_{min}\~(10^8\;M{\bigodot})(nh/0.05cm^{-3})^{-2}$ cool isochorically and thermal instabilities do not operate. On the other hand, intermediate size perturbations ($M_{min} < M_{pgcc} < M_{max}$) are compressed supersonically, accompanied by strong accretion shocks. 2) For supersonically collapsing clouds, the density compression factor after they cool to $T_c = 10^4$ K range $10^{2.5} - 10^6$, while the isobaric compression factor is only $10^{2.5}$. 3) Isobarically collapsed clouds ($M < M_{min}$) are too small to be gravitationally bound. For supersonically collapsing clouds, however, the Jeans mass can be reduced to as small as $10^{5.5}\;M_{\bigodot}(nh/0.05cm^{-3})^{-1/2}$ at the maximum compression owing to the increased density compression. 4) The density profile of simulated PGCCs can be approximated by a constant core with a halo of $p{\infty} r^{-2}$ rather than a singular isothermal sphere.

• Journal of The Korean Astronomical Society
• /
• v.36 no.3
• /
• pp.89-95
• /
• 2003

Observations of dark matter dominated dwarf and low surface brightness disk galaxies favor density profiles with a flat-density core, while cold dark matter (CDM) N-body simulations form halos with central cusps, instead. This apparent discrepancy has motivated a re-examination of the microscopic nature of the dark matter in order to explain the observed halo profiles, including the suggestion that CDM has a non-gravitational self-interaction. We study the formation and evolution of self-interacting dark matter (SIDM) halos. We find analytical, fully cosmological similarity solutions for their dynamics, which take proper account of the collisional interaction of SIDM particles, based on a fluid approximation derived from the Boltzmann equation. The SIDM particles scatter each other elastically, which results in an effective thermal conductivity that heats the halo core and flattens its density profile. These similarity solutions are relevant to galactic and cluster halo formation in the CDM model. We assume that the local density maximum which serves as the progenitor of the halo has an initial mass profile ${\delta}M / M {\propto} M^{-{\epsilon}$, as in the familiar secondary infall model. If $\epsilon$ = 1/6, SIDM halos will evolve self-similarly, with a cold, supersonic infall which is terminated by a strong accretion shock. Different solutions arise for different values of the dimensionless collisionality parameter, $Q {\equiv}{\sigma}p_br_s$, where $\sigma$ is the SIDM particle scattering cross section per unit mass, $p_b$ is the cosmic mean density, and $r_s$ is the shock radius. For all these solutions, a flat-density, isothermal core is present which grows in size as a fixed fraction of $r_s$. We find two different regimes for these solutions: 1) for $Q < Q_{th}({\simeq} 7.35{\times} 10^{-4}$), the core density decreases and core size increases as Q increases; 2) for $Q > Q_{th}$, the core density increases and core size decreases as Q increases. Our similarity solutions are in good agreement with previous results of N-body simulation of SIDM halos, which correspond to the low-Q regime, for which SIDM halo profiles match the observed galactic rotation curves if $Q {\~} [8.4 {\times}10^{-4} - 4.9 {\times} 10^{-2}]Q_{th}$, or ${\sigma}{\~} [0.56 - 5.6] cm^2g{-1}$. These similarity solutions also show that, as $Q {\to}{\infty}$, the central density acquires a singular profile, in agreement with some earlier simulation results which approximated the effects of SIDM collisionality by considering an ordinary fluid without conductivity, i.e. the limit of mean free path ${\lambda}_{mfp}{\to} 0$. The intermediate regime where $Q {\~} [18.6 - 231]Q_{th}$ or ${\sigma}{\~} [1.2{\times}10^4 - 2.7{\times}10^4] cm^2g{-1}$, for which we find flat-density cores comparable to those of the low-Q solutions preferred to make SIDM halos match halo observations, has not previously been identified. Further study of this regime is warranted.

• Clean Technology
• /
• v.25 no.1
• /
• pp.56-62
• /
• 2019

The equilibrium, kinetic and thermodynamic parameters of adsorption of murexide by granular activated carbon were investigated. The experiment was carried out by batch experiment with the variables of the amount of the adsorbent, the initial concentration of the dye, the contact time and the temperature. The isothermal adsorption equilibrium was best applied to the Freundlich equation in the range of 293 ~ 313 K. From the separation factor (${\beta}$) of Freundlich equation, it was found that adsorption of murexide by granular activated carbon could be the appropriate treatment method. The adsorption energy (E) obtained from the Dubinin- Radushkevich equation shows that the adsorption process is a physical adsorption process. From the kinetic analysis of the adsorption process, pseudo second order model is more consistent than pseudo first order model. It was found that the adsorption process proceeded to a spontaneous process and an endothermic process through Gibbs free energy change ($-0.1096{\sim}-10.5348kJ\;mol^{-1}$) and enthalpy change ($+151.29kJ\;mol^{-1}$). In addition, since the Gibbs free energy change decreased with increasing temperature, adsorption reaction of murexide by granular activated carbon increased spontaneously with increasing temperature. The entropy change ($147.62J\;mol^{-1}\;K^{-1}$) represented the increasing of randomness at the solid-solution interface during the adsorption reaction of murexide by activated carbon.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.22 no.1
• /
• pp.53-66
• /
• 1980

The research work is concerned with the analytical and experimental studies on the heat transfer phenomenon around the underground concrete digester used for biogas production Systems. A mathematical and computational method was developed to estimate heat losses from underground cylindrical concrete digester used for biogas production systems. To test its feasibility and to evaluate thermal parameters of materials related, the method was applied to six physical model digesters. The cylindrical concrete digester was taken as a physical model, to which the model,atical model of heat balance can be applied. The mathematical model was transformed by means of finite element method and used to analyze temperature distribution with respect to several boundary conditions and design parameters. The design parameters of experimental digesters were selected as; three different sizes 40cm by 80cm, 80cm by 160cm and l00cm by 200cm in diameter and height; two different levels of insulation materials-plain concrete and vermiculite mixing in concrete; and two different types of installation-underground and half-exposed. In order to carry out a particular aim of this study, the liquid within the digester was substituted by water, and its temperature was controlled in five levels-35。 C, 30。 C, 25。 C, 20。C and 15。C; and the ambient air temperature and ground temperature were checked out of the system under natural winter climate conditions. The following results were drawn from the study. 1.The analytical method, by which the estimated values of temperature distribution around a cylindrical digester were obtained, was able to be generally accepted from the comparison of the estimated values with the measured. However, the difference between the estimated and measured temperature had a trend to be considerably increased when the ambient temperature was relatively low. This was mainly related variations of input parameters including the thermal conductivity of soil, applied to the numerical analysis. Consequently, the improvement of these input data for the simulated operation of the numerical analysis is expected as an approach to obtain better refined estimation. 2.The difference between estimated and measured heat losses was shown to have the similar trend to that of temperature distribution discussed above. 3.It was found that a map of isothermal lines drawn from the estimated temperature distribution was very useful for a general observation of the direction and rate of heat transfer within the boundary. From this analysis, it was interpreted that most of heat losses is passed through the triangular section bounded within 45 degrees toward the wall at the bottom edge of the digesten Therefore, any effective insulation should be considered within this region. 4.It was verified by experiment that heat loss per unit volume of liquid was reduced as the size of the digester became larger For instance, at the liquid temperature of 35˚ C, the heat loss per unit volume from the 0. 1m$^3$ digester was 1, 050 Kcal/hr m$^3$, while at for 1. 57m$^3$ digester was 150 Kcal/hr m$^3$. 5.In the light of insulation, the vermiculite concrete was consistently shown to be superior to the plain concrete. At the liquid temperature ranging from 15。 C to 350 C, the reduction of heat loss was ranged from 5% to 25% for the half-exposed digester, while from 10% to 28% for the fully underground digester. 6.In the comparison of heat loss between the half-exposed and underground digesters, the heat loss from the former was fr6m 1,6 to 2, 6 times as much as that from the latter. This leads to the evidence that the underground digester takes advantage of heat conservation during winter.

• Journal of the Korea Concrete Institute
• /
• v.13 no.6
• /
• pp.610-618
• /
• 2001

This paper investigates the relationships between dynamic elastic modulus and static elastic modulus or compressive strength according to curing temperature, aging, and cement type. Based on this investigation, the new model of the relationships we proposed. Impact echo method estimates the resonant frequency of specimens and uniaxial compression test measures the static elastic modulus and compressive strength. Type I and V cement concretes, which have the water-cement ratios of 0.40 and 0.50, are cured under the isothermal curing temperatures of 10, 23, and 50$\^{C}$ Cement type and aging have no large influence on the relationship between dynamic and static elastic modulus, but the ratio of dynamic and static elastic modulus comes close to 1 as temperature increases. Initial chord elastic modulus which is calculated at lower strain level of stress-strain curve, has the similar value to dynamic elastic modulus. The relationship between dynamic elastic modulus and compressive strength has the same tendency as the relationship between dynamic and static elastic modulus according to cement type, temperature and aging. The proposcd relationship equations between dynamic elastic modulus and static elastic modulus or compressive strength properly estimates the variation of relationships according to cement type md temperature.

• Korean Chemical Engineering Research
• /
• v.50 no.1
• /
• pp.149-154
• /
• 2012

The adsorption and desorption behavior of biofilter-medium was investigated on the performance of an adsorption column. Continuous flow-isothermal adsorption experiments were performed to treat waste air containing such a VOC as ethanol under the same condition of > 90% relative humidity as the condition of the feed to a biofilter process. In case of feeding waste air containing ethanol of 1,000 ppmv (or 2,050 mg ethanol/$m^3$) to the adsorption system at the rate of 2 L/min, the onsets of its breakthrough and reaching the state of dynamic equilibrium at the exit had been delayed 10 and 3 times, respectively, later than those at the 1st stage sampling port. Moreover, in case of 2,000 ppmv (or 4,100 mg ethanol/$m^3$), they had been delayed 9 and 3 times, respectively. Thus, regardless of feeding concentration, the ratios of delaying period were observed to be quite consistent each other at the exit of the adsorption column. With regard to the period of desorption, the ratios of delaying period were consistent each other to be 1.5 for both cases. In addition, the effect of microbial activity and sterilization-process was studied on adsorption equilibrium. The ethanol concentration in the vapor phase of vials packed with sterilized granular activated carbon (GAC) was quite consistent to that with unsterilized GAC. However, the ethanol concentrations in the vapor phase of vials packed with unsterilized compost and the unsterilized mixture of GAC and compost were higher than those with sterilized compost and the sterilized mixture of GAC and compost, respectively.