• Journal of the Korean Society for Heat Treatment
• /
• v.14 no.4
• /
• pp.228-234
• /
• 2001

The present study was an attempt for systematic data conversion between FDM and FEM in order to evaluate the thermal stress distribution during quenching process. It has been generally recognized that FDM is efficient in flow and temperature analysis and FEM in that of stress. But it induced difficulty and tedious work in analysis that one uses both FDM and FEM to take their advantages because of the discrepancy of nodes between analysis tools. So we proposed field data conversion procedure from FDM to FEM in 3-dimensional space, then applied this procedure to analysis of quenching process. The simulation procedure calculates the distributions of temperature and microstructure using FDM and microstructure evolution equations of diffusion and diffusionless transformation. FEM was used for predicting the distributions of thermal stress. The present numerical code includes coupled temperaturephase transformation kinetics and temperature-microstructure dependent material properties. Calculated results were compared with previous experimental data to verify the method, which showed good agreements.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.21 no.3
• /
• pp.224-229
• /
• 2009

Partial safety factors(PSFs) for Level I reliability based design of caisson type quay walls were calculated. First order reliability method(FORM) based PSFs are the functions of sensitivities of limit state function with respect to design random variables, target reliability index, characteristic values and first moment of random variables. Modified PSFs for water level and resilient water level are newly defined to keep consistency with the current design code. In the numerical example, PSFs were calculated by using a target reliability index. Seismic coefficient is defined to show extreme distribution. It was found that PSFs for seismic coefficient becomes smaller as the return period for design seismic coefficient grows longer.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.1
• /
• pp.165-176
• /
• 1996

The flow field and spray characteristics for loop scavenged type 2stroke engine having pancake shape was numerically computed using KIVA-Ⅱ code. The cylinder has 1intake port, 2side intake ports and 1exhaust port with induced flow angle 25 deg. In engine calculation, the chop techniques is used to strip or add planes of cells across the mesh adjacent to the TDC and the BDC(ports parts) for preventing the demand of exceed time during the computation, providing a control on cell height in the squish region. The modified turbulent model including the consideration of the compressibility effect due to the compression and expansion of piston was also used. The case of 25 deg.(injection angle) which is opposite to scavenging flow direction shows better the distribution of droplets and the evaporation rate of droplets compared to other cases(0 deg., - 25 deg.). When injection pressure was increased, the spray tip penetration became longer. When injection pressure was increased, the interaction between the upward gas velocity and spray droplets strongly cause. Thus the breakup of droplets is strongly occurred and the evaporation rate of droplets was found to be better.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.7 no.8
• /
• pp.150-159
• /
• 1999

Exhaust system is composed of several parts. Among, them , design of muffler system strongly influences on engine efficiency and noise reduction. So , through comprehension of flow characteristics inside muffler is necessary . In this study , three-dimensional steady and unsteady compressible flow analysis was performed to understand the flow characteristics, pressure loss and amplitude variation of pulsating pressure. The computational grid generation was carried out using commercial preprocessor ICEM CFD/CAE. And the three-dimensional fluid motion inside the muffler was analyzed by STAR-CD, the computational fluid dynamics code. RNG k-$\varepsilon$ tubulence model was applied to consider the complexity of the geometry and fluid motion. The steady and unsteady flow field inside muffler such as velocity distribution, pulsating pressure and pressure loss was examined. In case of unsteady state analysis, velocity of inlet region was converted from measured pulsating pressure. Experimental measurement of pressure and temperature was carried out to provide the boundary and initial condition for computational study under three engine operating conditions. As a result of this study, we could identify the flow characteristics inside the muffler and obtain the pressure loss, amplitude variation of pulsating exhaust gas.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.43 no.9
• /
• pp.805-813
• /
• 2015

This paper introduces the development of temperature mapping code between thermal mesh and structural mesh in KARI Satellite Design Software. Generally, temperature distribution of a satellite varies with the time by the space environment of the orbit, so thermal expansion of the structure should be analysed in design of the satellite. For the sake of the coupled thermal structural analysis, an interpolation algorithm between two finite element heterogeneous grids has been proposed by which temperature transfer is successively conducted.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.8
• /
• pp.611-620
• /
• 2020

The effect of discharge voltage on electron mean energy, electric potential, ionization rate, neutral and ion density of Hall thruster was analyzed using a two-dimensional axisymmetric hybrid model. The results of the code developed for this study such as discharge current, thrust, and plasma distribution according to discharge voltage of SPT-100ML Hall thruster were compared by experiments and calculations of other researchers for validation. The results show that the electron mean energy, the ionization rate, and the ion density are increased while the neutral density is decreased as the discharge voltage is increased. The thrust and the discharge current are proportional to the discharge voltage.

• Journal of ILASS-Korea
• /
• v.8 no.4
• /
• pp.9-16
• /
• 2003

The Purpose of this study is to investigate the relationships between the local heat release rate and CH concentration have been investigated by numerical simulations of methane-air premixed flames. And simultaneous CH and OH PLIF(Planar Laser Induced Fluorescence) measurement has been also conducted for lean premixed flame as well as for laminar flames. Numerical simulations are conducted for laminar premixed flames and turbulent ones by using PREMIX in CHEMKIN and two dimensional DNS code with GRI mechanism version 2.11, respectively. In the case of laminar premixed flame, the distance between the peak of heat release rate and that of CH concentration is under $91{\mu}m$ for all equivalence ratio calculated in present work. Even for the premixed flame in high intensity turbulence, the distribution of the heat release rate coincides with that of CH mole fraction. For CH PLIF measurements in the laminar premixed flame burner, CH fluorescence intensity as a function of equivalence ratio shows a similar trend with CH mole fraction computed by GRI mechanism. Simultaneous CH and OH PLIF measurement gave us useful information of instantaneous reaction zone. In addition, CH fluorescence can be measured even for lean conditions where CH mole fraction significantly decreases compared with that of stoichiometric condition. It was found that CH PLIF measurements can be applicable to the estimation of the spatial fluctuations of heat release rate in the engine combustion.

• Nuclear Engineering and Technology
• /
• v.31 no.1
• /
• pp.17-28
• /
• 1999

A thermal-hydraulic analysis is conducted on the loss-of-coolant-accident (LOCA) during shutdown operation of YGN Units 3/4. Based on the review of plant-specific characteristics of YGN Units 3/4 in design and operation, a set of analysis cases is determined, and predicted by the RELAP5/MOD3.2 code during LOCA in the hot-standby mode. The evaluated thermal-hydraulic phenomena are blowdown, break flow, inventory distribution, natural circulation, and core thermal response. The difference in thermal-hydraulic behavior of LOCA at shutolown condition from that of LOCA at full power is identified as depressurization rate, the delay in peak natural circulation timing and the loop seal clearing (LSC) timing. In addition, the effect of high pressure safety injection (HPSI) on plant response is also evaluated. The break spectrum analysis shows that the critical break size can be between 1% to 2% of cold leg area, and that the available operator action time for the Sl actuation and the margin in the peak clad temperature (PCT) could be reduced when considering uncertainties of the present RELAP5 calculation.

• Nuclear Engineering and Technology
• /
• v.46 no.5
• /
• pp.619-632
• /
• 2014

The Best Estimate Plus Uncertainty (BEPU) method has been widely used to evaluate the uncertainty of a best-estimate thermal hydraulic system code against a figure of merit. This uncertainty is typically evaluated based on the physical model's uncertainties determined by expert judgment. This paper introduces the application of data assimilation methodology to determine the uncertainty bands of the physical models, e.g., the mean value and standard deviation of the parameters, based upon the statistical approach rather than expert judgment. Data assimilation suggests a mathematical methodology for the best estimate bias and the uncertainties of the physical models which optimize the system response following the calibration of model parameters and responses. The mathematical approaches include deterministic and probabilistic methods of data assimilation to solve both linear and nonlinear problems with the a posteriori distribution of parameters derived based on Bayes' theorem. The inverse problem was solved analytically to obtain the mean value and standard deviation of the parameters assuming Gaussian distributions for the parameters and responses, and a sampling method was utilized to illustrate the non-Gaussian a posteriori distributions of parameters. SPACE is used to demonstrate the data assimilation method by determining the bias and the uncertainty bands of the physical models employing Bennett's heated tube test data and Becker's post critical heat flux experimental data. Based on the results of the data assimilation process, the major sources of the modeling uncertainties were identified for further model development.

• Nuclear Engineering and Technology
• /
• v.48 no.5
• /
• pp.1109-1119
• /
• 2016

Although, the direct use of spent pressurized water reactor (PWR) fuel in CANda Deuterium Uranium (CANDU) reactors (DUPIC) cycle is still under investigation, DUPIC cycle is a promising method for uranium utilization improvement, for reduction of high level nuclear waste, and for high degree of proliferation resistance. This paper focuses on the effect of DUPIC cycle on CANDU reactor safety parameters. MCNP6 was used for lattice cell simulation of a typical 3,411 MWth PWR fueled by $UO_2$ enriched to 4.5w/o U-235 to calculate the spent fuel inventories after a burnup of 51.7 MWd/kgU. The code was also used to simulate the lattice cell of CANDU-6 reactor fueled with spent fuel after its fabrication into the standard 37-element fuel bundle. It is assumed a 5-year cooling time between the spent fuel discharges from the PWR to the loading into the CANDU-6. The simulation was carried out to calculate the burnup and the effect of DUPIC fuel on: (1) the power distribution amongst the fuel elements of the bundle; (2) the coolant void reactivity; and (3) the reactor point-kinetics parameters.