• Korean Chemical Engineering Research
• /
• v.46 no.5
• /
• pp.931-937
• /
• 2008

Rigorous multiscale modelling and simulation of the MTR for WGSR was carried out to accurately predict the behavior of process variables and the reactor performance. The MTR consists of 4 fixed bed tube reactors packed with heterogeneous catalysts, as well as surrounding shell part for the cooling purpose. Considering that fluid flow field and reaction kinetics give a great influence on the reactor performance, employing multiscale methodology encompassing Computational Fluid Dynamics (CFD) and process modeling was natural and, in a sense, inevitable conclusion. Inlet and outlet temperature of the reactant fluid at the tube side was $345^{\circ}C$ and $390^{\circ}C$, respectively and the CO conversion at the exit of the tube side with these conditions approached to about 0.89. At the shell side, the inlet and outlet temperature of the cooling fluid, which flows counter-currently to tube flow, was $190^{\circ}C$ and $240^{\circ}C$. From this heat exchange, the energy saving was achieved for the flow at shell side and temperature of the tube side was properly controlled to obtain high CO conversion. The simulation results from this research were accurately comparable to the experimental data from various papers.

• Journal of the Korean Chemical Society
• /
• v.55 no.2
• /
• pp.208-217
• /
• 2011

The effects of pH and temperature on the removal of two dyes (neutral red; NR and malachite green oxalates; MG) from aqueous effluents using Maghnia montmorillonite clay in a batch adsorption process were investigated. The results showed the stability of the optical properties of MG in aqueous solution and adsorbed onto clay under wide range of pH 3-9. However, the interaction of NR dye with clay is accompanied by a red shift of the main absorption bands of monomer cations under pH range of 3-5, whereas, those of neutral form remains nearly constant over the pH range of 8-12. The optimal pH for favorable adsorption of the dyes, i.e. ${\geq}$90% has been achieved in aqueous solutions at 6 and 7 for NR and VM respectively. The most suitable adsorption temperatures were 298 and 318 K with maximum adsorption capacities of 465.13mg/g for NR and 459.89 mg/g for MG. The adsorption equilibrium results for both dyes follow Langmuir, Freundlich isotherms. The numerical values of the mean free energy $E_a$ of 4.472-5.559 kj/mol and 2.000-2.886 kj/mol for NR and MG respectively indicated physical adsorption. Various thermodynamic parameters, such as ${\Delta}H^{\circ}$, ${\Delta}S^{\circ}$, ${\Delta}G^{\circ}$ and Ea have been calculated. The data showed that the adsorption process is spontaneous and endothermic. The sticking probability model was further used to assess the potential feasibility of the clay mineral as an alternative adsorbent for organic ion pollutants in aqueous solution.

• Clean Technology
• /
• v.20 no.1
• /
• pp.13-21
• /
• 2014

Because of the water scarcity caused by the increase of salinity in the underground water, seawater desalination stands out as one of the most promising solution. As there are so much energy costs in operating desalination plants, new hybrid process which is more effective should be researched. A geothermal VMD (vacuum membrane distillation) hybrid process is a competitive alternative for seawater desalination. Because geothermal energy has significant characteristics of high capacity factor to operate the power plant at full capacity for 24 hour per day, it can be a priority heat source of VMD superior to any other renewable energies such as solar and wind power. In this study, we design a geothermal VMD hybrid process, analyze it economically and finally compare the result with a case of conventional VMD process. Geothermal VMD hybrid process generates $23,822,409 of NPV (net present value) more than the conventional VMD process in case of 5% discount rate. The break-even point between these processes is 5.36 year. Sensitivity analysis indicates that steam cost is the most decisive influence variable to the economic outcome.

• Clean Technology
• /
• v.13 no.2
• /
• pp.151-158
• /
• 2007

An experimental study to investigate the condensing heat transfer characteristics of small diameter horizontal double pipe heat exchangers with R-22 and R-410A was performed. Experimental facility was constructed to calculate and observe HTC(heat transfer coefficients), flow patterns and pressure drop. The main components include a liquid pump, an evaporator, a condenser(test section), a sight-glass, pressure taps and measurement apparatus. Two pipes of different diameters are tested; One 5.35 mm ID 0.5 mm thick, the other 3.36 mm ID 0.7 mm thick. The mass flow rate ranged from 200 to $500\;ks/m^2{\cdot}s$ and heating capacity were form 1.0 to 2.4 kW. The flow patterns of R-22 and R-410A were observed with a high speed camera through the sight-glass. The tests revealed that HTC of R-410A was higher than that of R-22 by maximum 5%. Annular pattern was observed for the most cases but stratified flow was also detected when x<0.2. The pressure drop in 3.36 mm ID pipe was higher than that of 5.35 mm by $30{\sim}50%$. Comparing with previous correlations such as Shah, Fujii and Soliman's, Fujii' showed the best good agreement with my data with a maximum deviation of 40%.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.491-499
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed tot the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted tot in the sense of yon Karman hypothesis. The material model is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized lot anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed P-version finite element model is demonstrated through several comparative points of iew in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic tone.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.1 s.71
• /
• pp.93-103
• /
• 2006

It is absolutely essential that safety assessment of the containment buildings during service life because containment buildings are last barrier to protect radioactive substance due to the accidents. Therefore, this study describes an enhanced degenerated shell finite element(FE) which has been developed for nonlinear FE analysis of reinforced concrete(RC) containment buildings with elasto-plastic material model. For the purpose of the material nonlinear analysis, Drucker-Prager failure criteria is adapted in compression region and material parameters which determine the shape of the failure envelop are derived from biaxial stress tests. Reissner-Mindlin(RM) assumptions are adopted to develop the degenerated shell FE so that transverse shear deformation effects is considered. However, it is found that there are serious defects such as locking phenomena in RM degenerated shell FE since the stiffness matrix has been overestimated in some situations. Therefore, shell formulation is provided in this paper with emphasis on the terms related to the stiffness matrix based on assumed strain method. Finally, the performance of the present shell element to analysis RC containment buildings is tested and demonstrated with several numerical examples. From the numerical tests, the present results show a good agreement with experimental data or other numerical results.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.2
• /
• pp.145-151
• /
• 2008

The research so far has primarily analyzed efficiency improvement but in this research, it analyzes the characteristics of earthquake behavior, with changed pier heights, through ordinary and seismic analysis. For this, the kind of bridge bearing has been changed against PSC I-shaped bridge, which is mostly used in practice, and at all times earthquake analysis has been performed with through height of pier. Especially considering sectional power resulting from earthquake analysis, displacement of PSC I-shaped bridge bearing, diameter of pier pillar by earthquake load, and upper spare gap have been analyzed. In case of high-pear, seismic isolated device is decided as proper for cars' driving and for management of bridge since it decreases movement of upper structure, than elastic bearing, reducing size of elastic connect device, and it's been analyzed it is effective for improvement of fine view and economic efficiency reducing section of lower bridge structure. Finally, when design PSC I-shaped bridge bearing, for the proper structure and high-pier side, applying seismic isolated device through precise inner analysis is proper than applying equal elastic bearing.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.6
• /
• pp.635-642
• /
• 2014

In this paper, an efficient technique is developed to predict failure probability of three failure modes(case rupture, fracture and bolt breakage) related to solid rocket motor case due to the inner pressure during the mission flight. The overall procedure consists of the steps: 1) design parameters affecting the case failure are identified and their uncertainties are modelled by probability distribution, 2) combustion analysis in the interior of the case is carried out to obtain maximum expected operating pressure(MEOP), 3) stress and other structural performances are evaluated by finite element analysis(FEA), and 4) failure probabilities are calculated for the above mentioned failure modes. Axi-symmetric assumption for FEA is employed for simplification while contact between bolted joint is accounted for. Efficient procedure is developed to evaluate failure probability which consists of finding first an Most Probable Failure Point(MPP) using First-Order Reliability Method(FORM), next making a response surface model around the MPP using Latin Hypercube Sampling(LHS), and finally calculating failure probability by employing Importance Sampling.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.2
• /
• pp.63-70
• /
• 2014

Design sensitivity analysis and topology design optimization for a piezoelectric crystal resonator are developed. The piezoelectric crystal resonator is deformed mechanically when subjected to electric charge on the electrodes, or vice versa. The Mindlin plate theory with higher-order interpolations along thickness direction is employed for analyzing the thickness-shear vibrations of the crystal resonator. Thin electrode plates are masked on the top and bottom layers of the crystal plate in order to enforce to vibrate it or detect electric signals. Although the electrode is very thin, its weight and shape could change the performance of the resonators. Thus, the design variables are the bulk material densities corresponding to the mass of masking electrode plates. An optimization problem is formulated to find the optimal topology of electrodes, maximizing the thickness-shear contribution of strain energy at the desired motion and restricting the allowable volume and area of masking plates. The necessary design gradients for the thickness-shear frequency(eigenvalue) and the corresponding mode shape(eigenvector) are computed very efficiently and accurately using the analytical design sensitivity analysis method using the eigenvector expansion concept. Through some demonstrative numerical examples, the design sensitivity analysis method is verified to be very efficient and accurate by comparing with the finite difference method. It is also observed that the optimal electrode design yields an improved mode shape and thickness-shear energy.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.6
• /
• pp.521-527
• /
• 2016

This paper is concerned with the numerical analysis of dynamic response of floating offshore wind turbine subject to underwater explosion using an effective non-reflecting technique. An infinite sea water domain was truncated into a finite domain, and the non-reflecting technique called the perfectly matched layer(PML) was applied to the boundary of truncated finite domain to absorb the inherent reflection of out-going impact wave at the boundary. The generalized transport equations that govern the inviscid compressible water flow was split into three PML equations by introducing the direction-wise absorption coefficients and state variables. The fluid-structure interaction problem that is composed of the wind turbine and the sea water flow was solved by the iterative coupled Eulerian FVM and Largangian FEM. And, the explosion-induced hydrodynamic pressure was calculated by JWL(Jones-Wilkins-Lee) equation of state. Through the numerical experiment, the hydrodynamic pressure and the structural dynamic response were investigated. It has been confirmed that the case using PML technique provides more reliable numerical results than the case without using PML technique.