• 한국신재생에너지학회:학술대회논문집
• /
• 2008.05a
• /
• pp.507-510
• /
• 2008

Characteristics of an autothermal reformer at various operating parameters have been studied in this paper. Numerical method has been used, and simulation model has been developed for the analysis. Full Combustion reaction, Steam Reforming(SR) reaction, Water-Gas Shift(WGS) reaction, and Direct Steam Reforming(DSR) reaction are assumed as dominant chemical reactions in the autothermal reformer. Simulation results are compared with experimental results for code validation. Operating parameters of the autothermal reformer are inlet temperature, Oxygen to Carbon Ratio(OCR), Steam to Carbon Ratio(SCR), and Gas Hourly Space Veolcity(GHSV). SR reaction rate decreases with low inlet temperature. If OCR is increased, $H_2$ yield is increased but optimal point is suggested. WGS reaction is activated with high SCR. When GHSV is increased, reforming efficiency is increased but pressure drop may decrease the system efficiency.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.12
• /
• pp.1501-1510
• /
• 2004

The pumping characteristics of a single-stage disk-type drag pump (DTDP) are calculated for the variation of the vertical clearance between a rotor and stator by the three-dimensional direct simulation Monte Carlo (DSMC) method. The gas flow mainly belongs to the molecular transition flow region. Spiral channels of a DTDP are cut on the both the upper and lower sides of a rotating disk, but a stationary disk is planar. The interaction between molecules is described by the variable hard-sphere model. The no time counter method is used as a collision sampling technique. The vertical clearance has a significant effect on the pumping performance. Experiments are performed under the outlet pressure range of 0.4∼533 Pa. When the numerical results are compared with the experimental data, the numerical results agree well quantitatively

• Proceedings of the Optical Society of Korea Conference
• /
• 2003.02a
• /
• pp.184-185
• /
• 2003

Unlike the conventional waveguide such as optical fiber using total internal reflection, photonic crystal waveguide(PCW), a waveguide made of a line defect in a photonic crystal(PC) structure, does not admit an analytic approach due to its complexity but requires a direct numerical approach. Here, we present numerical results of computer simulation for PCW by using the three-dimensional(3D) Finite-Difference Time -Domain(FDTD) algorithm. (omitted)

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1059-1069
• /
• 2023

Dam-break flow occurs when an elevated dam suddenly collapses, resulting in the catastrophic release of rapid and uncontrolled impounded water. This study compares laminar and turbulent closure models for simulating three-dimensional dam-break flows using OpenFOAM. The Reynolds-Averaged Navier-Stokes (RANS) model, specifically the k-ε model, is employed to capture turbulent dissipation. Two scenarios are evaluated based on a laboratory experiment and a modified multi-layered block obstacle scenario. Both models effectively represent dam-break flows, with the turbulent closure model reducing oscillations. However, excessive dissipation in turbulent models can underestimate water surface profiles. Improving numerical schemes and grid resolution enhances flow recreation, particularly near structures and during turbulence. Model stability is more significantly influenced by numerical schemes and grid refinement than the use of turbulence closure. The k-ε model's reliance on time-averaging processes poses challenges in representing dam-break profiles with pronounced discontinuities and unsteadiness. While simulating turbulence models requires extensive computational efforts, the performance improvement compared to laminar models is marginal. To achieve better representation, more advanced turbulence models like Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) are recommended, necessitating small spatial and time scales. This research provides insights into the applicability of different modeling approaches for simulating dam-break flows, emphasizing the importance of accurate representation near structures and during turbulence.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.28 no.6
• /
• pp.483-495
• /
• 2008

The paper presents a numerical study to evaluate the phase velocities and attenuations of the average longitudinal and shear ultrasonic waves resulting from multiple scattering in fiber-reinforced composites. A computational procedure developed in this work is first used to produce a random, yet largely even distribution of fibers. Both the viscoelastic epoxy matrix and lossless randomly distributed graphite fibers are modeled using the mass-spring-dashpot lattice model, with no damping for the latter. By numerically simulating ultrasonic through-transmission tests using this direct model of composites, phase velocities and attenuations of the longitudinal and shear waves through the composite are found as functions of frequency or fiber concentration. The numerical results are observed to generally agree with the corresponding results in the literature. Discrepancies found in some detail aspects, particularly in the attenuation results, are also addressed.

• Journal of the Korean Society of Visualization
• /
• v.15 no.3
• /
• pp.27-33
• /
• 2017

Direct numerical simulations of turbulent channel flows with moving wall conditions on the top wall are performed to examine the effects of the moving wall on the turbulent characteristics. The moving wall velocity only applied to the top wall with the opposite direction to the main flow is systematically varied to reveal the sustained-mechanism for turbulence. The turbulence statistics for the Couette-Poiseuille flow, such as mean velocity, root mean square of the velocity fluctuations, Reynolds shear stress and pre-multiplied energy spectra of the velocity fluctuations, are compared with those of canonical turbulent channel flows. The comparison suggests that although the turbulent activity on the top wall increases with increasing the Reynolds number, that on the bottom wall decreases, contrary to the previous finding for the canonical turbulent channel flows. The increase of the turbulent energy on the top wall is attributed to not only the increase of the Reynolds number but also elongation of the logarithmic layer due to increase of the wall layer on the top wall. However, because the logarithmic layer is shortened on the bottom wall due to the decrease of the wall layer, the turbulence energy on the bottom wall decreases despite of the increase of the Reynolds number.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.45-52
• /
• 2008

In order to investigate the initiation and propagation processes of a spherical detonation wave induced by direct initiation, numerical simulations were carried out using two-dimensional compressible Euler equations with an axisymmetric assumption and a one-step reaction model based on Arrhenius kinetics with various levels of ignition energy. By varying the amount of ignition energy, three typical initiation behaviors, which were subcritical, supercritical and critical regimes, were observed. Then, the ignition energy of more than $137.5{\times}10^6$ in non-dimensional value was required for initiating a spherical detonation wave, and the minimum ignition energy(i.e., critical energy) was less than that of the one-dimensional simulation reported by a previous numerical work. When the ignition energy was less than the critical energy, the blast wave generated from an ignition source continued to attenuate due to the separation of the blast wave and a reaction front. Therefore, detonation was not initiated in the subcrtical regime. When the ignition energy was more than the minimum initiation energy, the blast wave developed into a multiheaded detonation wave propagating spherically at CJ velocity, and then a cellular pattern radiated regularly out from the ignition center in the supercritical regime. The influence on ignition energy was observed in the cell width near the ignition center, but the cell width on the fully developed detonation remained constant during the expanding of detonation wave due to the consecutive formation of new triple points, regardless of ignition energy. When the ignition energy was equal to the critical energy, the decoupling of the blast wave and a reaction front appeared, as occurred in the subcrtical regime. After that, the detonation bubble induced by the local explosion behind the blast wave expanded and developed into the multiheaded detonation wave in the critical regime. Although few triple points were observed in the vicinity of the ignition core, the regularly located cellular pattern was generated after the onset of the multiheaded detonation. Then, the average cell width on the fully developed detonation was almost to that in the supercritical regime. These numerical results qualitatively agreed with previous experimental works regarding the initiation and propagation processes.

• Journal of the Korean Institute of Electrical and Electronic Material Engineers
• /
• v.35 no.6
• /
• pp.616-621
• /
• 2022

Numerical simulation is a good way to predict the conversion efficiency of solar cells without a direct experimentation and to achieve low cost and high efficiency through optimizing each step of solar cell fabrication. TOPCon industrial solar cells fabricated with n-type silicon wafers on a larger area have achieved a higher efficiency than p-type TOPCon solar cells. Electrical and optical losses of the front surface are the main factors limiting the efficiency of the solar cell. In this work, an optimization of boron-doped emitter surface and front electrodes through numerical simulation using "Griddler" is reported. Through the analysis of the results of simulation, it was confirmed that the emitter sheet resistance of 150 Ω/sq along the front electrodes having a finger width of 20 ㎛, and the number of finger lines ~130 for silicon wafer of M6 size is an optimized technology for the front emitter surface of the n-type TOPCon solar cells that can be developed.

• Structural Engineering and Mechanics
• /
• v.32 no.1
• /
• pp.1-20
• /
• 2009

The present contribution addresses uncertainty quantification and uncertainty propagation in structural mechanics using stochastic analysis. Presently available procedures to describe uncertainties in load and resistance within a suitable mathematical framework are shortly addressed. Monte Carlo methods are proposed for studying the variability in the structural properties and for their propagation to the response. The general applicability and versatility of Monte Carlo Simulation is demonstrated in the context with computational models that have been developed for deterministic structural analysis. After discussing Direct Monte Carlo Simulation for the assessment of the response variability, some recently developed advanced Monte Carlo methods applied for reliability assessment are described, such as Importance Sampling for linear uncertain structures subjected to Gaussian loading, Line Sampling in linear dynamics and Subset simulation. The numerical example demonstrates the applicability of Line Sampling to general linear uncertain FE systems under Gaussian distributed excitation.

• Structural Engineering and Mechanics
• /
• v.66 no.3
• /
• pp.379-386
• /
• 2018

The present study addresses the direct and indirect methods of determining the mode-I fracture toughness of concrete using experimental tests and particle flow code. The direct method used is compaction tensile test and the indirect methods are notched Brazilian disc test, semi-circular bend specimen test, and hollow center cracked disc. The experiments were carried out to determine which indirect method yields the fracture toughness closer to the one obtained by the direct method. In the numerical analysis, the PFC model was first calibrated with respect to the data obtained from the Brazilian laboratory test. The crack paths observed in the simulated tests were in reasonable accordance with experimental results. The discrete element simulations demonstrated that the macro fractures in the models are caused by microscopic tensile breakages on large numbers of bonded particles. The mode-I fracture toughness in the direct tensile test was smaller than the indirect testing results. The fracture toughness obtained from the SCB test was closer to the direct test results. Hence, the semi-circular bend test is recommended as a proper experiment for determination of mode-I fracture toughness of concrete in the absence of direct tests.