• Proceedings of the Korean Geotechical Society Conference
• /
• 2010.09a
• /
• pp.808-812
• /
• 2010

A great deal of attention is focused on coupled Thermo-Hydro-Mechanical (THM) behavior of multiphase porous media in diverse geo-mechanical and geo-environmental areas. This paper presents general governing equations for coupled THM processes in unsaturated porous media. Coupled partial differential equations are derived from 3 mass balances equations (solid, water, and air), energy balance equation, and force equilibrium equation. Finite element code is developed from the Galerkin formulation and time integration of these governing equations for 4 main variables (displacement $\underline{u}$, gas pressure $P_g$, liquid pressure $P_l$), and temperature T). The code is validated with theoretical solutions for linear material with simple boundary conditions.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.26-27
• /
• 2006

Characteristics of spatial structure of particle clusters are investigated by using the flow field data obtained from three-dimensional numerical simulations. Eulerian/Lagrangian approach with two-way coupling is applied and individual particle-particle collisions are taken into account by using the hard-sphere model. More than 16 million particles are traced in the maximum case. The results show that the cluster is consisted from the multiple-spatial scale components while the low wave-number, hence the large-scale structure, is dominant. Three-dimensional structure reconstructed from the low-pass filtered data enables us to investigate the essential dynamics of particle clusters in detail.

• Transactions of Materials Processing
• /
• v.20 no.5
• /
• pp.377-386
• /
• 2011

During stamping processes, the air trapped between sheet metal and the die cavity can be highly compressed and ultimately reduce the shape accuracy of formed panels. To prevent this problem, vent holes and passages are sometimes drilled into the based on expert experience and know-how. CAE can be also used for analyzing the air behavior in die cavity during stamping process, incorporating both elasto-plastic behavior of sheet metal and the fluid dynamic behavior of air. This study presents sheet metal forming simulation combined simultaneously with simulation of air behavior in the die cavity. There are three approaches in modeling of air behavior. One is a simple assumption of the bulk modulus having a constant pressure depending on volume change. The next is the use of the ideal gas law having uniform pressure and temperature in air domain. The third is FPM (Finite point method) having non-uniform pressure in air domain. This approach enables direct coupling of mechanical behavior of solid sheet metal and the fluid behavior of air in sheet metal forming simulation, and its result provides the first-hand idea for the location, size and number of the vent holes. In this study, commercial software, PAM-$STAMP^{TM}$ and PAM-$SAFE^{TM}$, were used.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.388-391
• /
• 2017

In this study, heat transfer analysis of variable thrust control system have been conducted by using commercial CFD code and FEM code. We Carried out computational fluid dynamics analysis to obtain the temperature and convective heat transfer coefficient of hot gas of variable thrust control system. Data are used as boundary condition for heat transfer analysis using mapping method. Temperature of O-ring for sealing was predicted

• Advances in Energy Research
• /
• v.5 no.2
• /
• pp.91-105
• /
• 2017

Thermal hydraulic (TH) analysis of nuclear power reactors is utmost important. In this way, the numerical codes that preparing TH data in reactor core are essential. In this paper, a subchannel analysis of a Russian pressurized water reactor (WWER1000) core with enhanced numerical code is carried out. For this, in fluid domain, the mass, axial and lateral momentum and energy conservation equations for desired control volume are solved, numerically. In the solid domain, the cylindrical heat transfer equation for calculation of radial temperature profile in fuel, gap and clad with finite difference and finite element solvers are considered. The dependence of material properties to fuel burnup with Calza-Bini fuel-gap model is implemented. This model is coupled with Isotope Generation and Depletion Code (ORIGEN2.1). The possibility of central hole consideration in fuel pellet is another advantage of this work. In addition, subchannel to subchannel and subchannel to rod connection data in hexagonal fuel assembly geometry could be prepared, automatically. For a demonstration of code capability, the steady state TH analysis of a the WWER1000 core is compromised with Thermal-hydraulic analysis code (COBRA-EN). By thermal hydraulic parameters averaging Fuel Assembly-to-Fuel Assembly method, the one sixth (symmetry) of the Boushehr Nuclear Power Plant (BNPP) core with regular subchannels are modeled. Comparison between the results of the work and COBRA-EN demonstrates some advantages of the presented code. Using the code the thermal modeling of the fuel rods with considering the fission gas generation would be possible. In addition, this code is compatible with neutronic codes for coupling. This method is faster and more accurate for symmetrical simulation of the core with acceptable results.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.3
• /
• pp.75-83
• /
• 2011

Emerging issues related with fully coupled Thermo-Hydro-Mechanical (THM) behavior of unsaturated soil demand the development of a numerical tool in diverse geo-mechanical and geo-environmental areas. This paper presents general governing equations for coupled THM processes in unsaturated porous media. Coupled partial differential equations are derived from three mass balances equations (solid, water, and air), energy balance equation, and force equilibrium equation. With Galerkin formulation and time integration of these governing equations, finite element code is developed to find nonlinear solution of four main variables (displacement-u, gas pressure-$P_g$), liquid pressure-$P_1$), and temperature-T) using Newton's iterative scheme. Three cases of numerical simulations are conducted and discussed: one-dimensional drainage experiments (u-$P_g-P_1$), thermal consolidation (u-$P_1$-T), and effect of pile on surrounding soil due to surface temperature variation (u-$P_1$-T).

• Earthquakes and Structures
• /
• v.26 no.4
• /
• pp.261-267
• /
• 2024

With the increase in the exploitation depth of offshore oil and gas, it is possible to control the global buckling of deep-sea pipelines by the snake lay method. Previous studies mainly focused on the analysis of critical buckling force and critical temperature of pipelines under the snake-like laying method, and pipelines often suffer structural failure due to seismic disasters during operation. Therefore, seismic action is a necessary factor in the design and analysis of submarine pipelines. In this paper, the seismic action of steel pipes in the operation stage after global buckling has occurred under the active control method is analyzed. Firstly, we have established a simplified finite element model for the entire process cycle and found that this modeling method is accurate and efficient, solving the problem of difficult convergence of seismic wave and soil coupling in previous solid analysis, and improving the efficiency of calculations. Secondly, through parameter analysis, it was found that under seismic action, the pipe diameter mainly affects the stress amplitude of the pipeline. When the pipe wall thickness increases from 0.05 m to 0.09 m, the critical buckling force increases by 150%, and the maximum axial stress decreases by 56%. In the pipe soil interaction, the greater the soil viscosity, the greater the pipe soil interaction force, the greater the soil constraint on the pipeline, and the safer the pipeline. Finally, the pipeline failure determination formula was obtained through dimensionless analysis and verified, and it was found that the formula was accurate.

• Journal of the Korean Chemical Society
• /
• v.42 no.6
• /
• pp.618-628
• /
• 1998

Pure CaO, Mn-doped CaO, Mn/CaO, and K/CaO catalysts were prepared and tested as catalysts for the oxidative coupling of methane in the temperature range of 600 to 800$^{\circ}C$ to investigate the effects of Mn- and K-addition on the catalytic activity of calcium oxide. To characterize the catalysts, X-ray powder diffraction(XRD), XPS, SEM, DSC, and TG analyses were performed. The catalytic reaction was carried out in a single-pass flow reactor using on-line gas chromatography system. Normalized reaction conditions were generally $p(CH_4)/p(O_2)=250$ Torr/50 Torr, total feed flow rate=30 mL/min, and 1 atm of total pressure with He being used as diluent gas. Among the catalysts tested, 6.3 mol% Mn-doped CaO catalyst showed the best $C_2$ yield of 8.0% with a selectivity of 43.2% at 775$^{\circ}C$. The $C_2$ selectivity increased on lightly doped CaO catalysts, while decreased on heavily doped CaO([Mn] > 6.3 mol%) catalysts. 6 wt.% Mn/CaO and 6 wt.% K/CaO catalysts showed the $C_2$ selectivities of 13.2% and 30.9%, respectively, for the reaction. Electrical conductivities of CaO and Mn-doped CaO were measured in the temperature range of 500 to 1000$^{\circ}C$ at Po2's of $10^{-3}\; to\;10^{-1}\;atm.$ The electrical conductivity was decreased with Mn-doping and increased with increasing $P0_2$in the range of $10^{-3}\;to\;10^{-1}\;atm,$ indicating the specimens to be p-type semiconductors. It was suggested that the interstitial oxygen ions formed near the surface can activate methane and the formation of interstitial oxygen ions was discussed on the basis of solid-state chemistry.