• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2013년도 춘계학술대회 논문집
• /
• pp.528-532
• /
• 2013

본 연구에서는, 해석적 모델로 불규칙하게 분포된 질량을 가진 열탄성 댐핑을 포함하는 마이크로빔 구조물을 연구하였다. 마이크로 스케일의 기계적 공명체(mechanical resonator)에 대한 열탄성 댐핑의 중요성은 높은 Q-factor를 설계하는데 고려된다. 본 연구에서의 빔 모델은 Euler-Bernoulli 빔 이론을 기조로 한다. 빔의 고유 진동수를 결정하기 위하여, 에너지 기법이 적용되었다. 또한, 열탄성 댐핑 효과는 열전도 방정식을 사용할으로써 고려되었고, Q-factor가 결정될 수 있었다. 운동방정식의 유도에는 체계적인 무차원화를 수행하였다. 임의의 집중된 질량을 포함하는 열탄성 댐핑을 가진 마이크로빔에 대해 모델의 결과값을 입증하였고 mode shape과 Q-factor를 제시하였다.

• Nuclear Engineering and Technology
• /
• 제46권5호
• /
• pp.641-654
• /
• 2014

A new computer code, named CORONA (Core Reliable Optimization and thermo-fluid Network Analysis), was developed for the core thermo-fluid analysis of a prismatic gas cooled reactor. The CORONA code is targeted for whole-core thermo-fluid analysis of a prismatic gas cooled reactor, with fast computation and reasonable accuracy. In order to achieve this target, the development of CORONA focused on (1) an efficient numerical method, (2) efficient grid generation, and (3) parallel computation. The key idea for the efficient numerical method of CORONA is to solve a three-dimensional solid heat conduction equation combined with one-dimensional fluid flow network equations. The typical difficulties in generating computational grids for a whole core analysis were overcome by using a basic unit cell concept. A fast calculation was finally achieved by a block-wise parallel computation method. The objective of the present paper is to summarize the motivation and strategy, numerical approaches, verification and validation, parallel computation, and perspective of the CORONA code.

• 한국자기학회지
• /
• 제8권3호
• /
• pp.161-168
• /
• 1998

본 논문에서는 라플라스 변환과 유환요소법의 결합에 의하여 확산반정식의 과도해석에 적용이 가능한 알고리즘을 제안하였다. 제안한 방법은 시간항을 라플라스 변환을 이용아여 제거한 후 유한요소법을 적용하여 해를 구한다. 이렇게 주파수 영역에서 구해진 해는 라플라스 역변환을 이용하여 시간영역의 값으로 변환한다. 제안된 방법의 타당성을 검증하기 위하여 열전도문제를 해석하엿으며, 제안한 방법이 해석해와 잘 이치한다는 것을 알 수 잇었다. 제안한 방법은 시간 차분이 필요하지 않기 때문에 여러 가지 확산방정식을 해석함에 있어서 매우 유용할 것으로 사료된다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제25권4호
• /
• pp.196-220
• /
• 2021

This study formulates a new geometric parameterization scheme to effectively address numerical analysis subject to the variation of the fiber radius of a square unit cell. In particular, the proposed mesh-morphing approach may lead to a parameterized weak form whose bilinear and linear forms are affine in the geometric parameter of interest, i.e. the fiber radius. As a result, we may certify the reduced basis analysis of a square unit cell model for any parameters in a predetermined parameter domain with a rigorous a posteriori error bound. To demonstrate the utility of the proposed geometric parameterization, we consider a two-dimensional, steady-state heat conduction analysis dependent on two parameters: a fiber radius and a thermal conductivity. For rapid yet rigorous a posteriori error evaluation, we estimate a lower bound of a coercivity constant via the min-θ method as well as the successive constraint method. Compared to the corresponding finite element analysis, the constructed reduced basis analysis may yield nearly the same solution at a computational speed about 29 times faster on average. In conclusion, the proposed geometric parameterization scheme is conducive for accurate yet efficient reduced basis analysis.

• 한국재료학회지
• /
• 제29권9호
• /
• pp.547-552
• /
• 2019

In the present study, the thermal conductivity of a silicon nitride($Si_3N_4$) thin-film is evaluated using the dual-wavelength pump-probe technique. A 100-nm thick $Si_3N_4$ film is deposited on a silicon (100) wafer using the radio frequency plasma enhanced chemical vapor deposition technique and film structural characteristics are observed using the X-ray reflectivity technique. The film's thermal conductivity is measured using a pump-probe setup powered by a femtosecond laser system of which pump-beam wavelength is frequency-doubled using a beta barium borate crystal. A multilayer transient heat conduction equation is numerically solved to quantify the film property. A finite difference method based on the Crank-Nicolson scheme is employed for the computation so that the experimental data can be curve-fitted. Results show that the thermal conductivity value of the film is lower than that of its bulk status by an order of magnitude. This investigation offers an effective way to evaluate thermophysical properties of nanoscale ceramic and dielectric materials with high temporal and spatial resolutions.

• Nuclear Engineering and Technology
• /
• 제53권11호
• /
• pp.3653-3664
• /
• 2021

In PWR three-dimensional transient coupling calculation software CORCA-K, the nodal Green's function method and diagonal implicit Runge Kutta method are used to solve the spatiotemporal neutron dynamic diffusion equation, and the single-phase closed channel model and one-dimensional cylindrical heat conduction transient model are used to calculate the coolant temperature and fuel temperature. The LMW, NEACRP and PWR MOX/UO₂ benchmarks and FangJiaShan (FJS) nuclear power plant (NPP) transient control rod move cases are used to verify the CORCA-K. The effects of burnup, fuel effective temperature and ejection rate on the control rod ejection process of PWR are analyzed. The conclusions are as follows: (1) core relative power and fuel Doppler temperature are in good agreement with the results of benchmark and ADPRES, and the deviation between with the reference results is within 3.0% in LMW and NEACRP benchmarks; 2) the variation trend of FJS NPP core transient parameters is consistent with the results of SMART and ADPRES. And the core relative power is in better agreement with the SMART when weighting coefficient is 0.7. Compared with SMART, the maximum deviation is -5.08% in the rod ejection condition and while -5.09% in the control rod complex movement condition.

• Nuclear Engineering and Technology
• /
• 제55권9호
• /
• pp.3213-3228
• /
• 2023

As an important part of the digital reactor, the pin-by-pin wise fine coupling calculation is a research hotspot in the field of nuclear engineering in recent years. It provides more precise and realistic simulation results for reactor design, operation and safety evaluation. CORCA-K a nodal code is redeveloped as a robust pin-by-pin wise neutronics and thermal-hydraulic coupled calculation code for pressurized water reactor (PWR) core. The nodal green's function method (NGFM) is used to solve the three-dimensional space-time neutron dynamics equation, and the single-phase single channel model and one-dimensional heat conduction model are used to solve the fluid field and fuel temperature field. The mesh scale of reactor core simulation is raised from the nodal-wise to the pin-wise. It is verified by two benchmarks: NEACRP 3D PWR and PWR MOX/UO₂. The results show that: 1) the pin-by-pin wise coupling calculation system has good accuracy and can accurately simulate the key parameters in steady-state and transient coupling conditions, which is in good agreement with the reference results; 2) Compared with the nodal-wise coupling calculation, the pin-by-pin wise coupling calculation improves the fuel peak temperature, the range of power distribution is expanded, and the lower limit is reduced more.

• 대한기계학회논문집B
• /
• 제39권7호
• /
• pp.567-578
• /
• 2015

Leidenfrost 온도 이상으로 가열된 벽면과 충돌하는 액적의 속도가 열전달 특성에 미치는 영향에 관한 실험 연구를 수행하였다. 동기화된 초고속 가시화 카메라와 적외선 카메라를 이용하여 벽면과 충돌하는 액적의 충돌 특성과 충돌면의 온도 분포를 측정하였다. 획득한 표면온도 분포를 충돌면의 경계 조건으로 이용하여 가열 벽면의 3차원 비정상 열전도 수치해석을 통해 표면 열유속 분포를 얻었다. 수직방향 충돌속도가 증가할수록 최대 액막 직경이 증가하고 가열 벽면과 액막 사이에 존재하는 증기막의 두께가 감소하여 열전달 효율이 증가하였다. 액적은 웨버수가 30보다 작은 경우 되튐현상이 발생하였으며, 큰 경우 작은 액적들로 분쇄되어졌다. 충돌속도에 의한 열전달량의 증가 경향이 되튐영역에서 분쇄영역에서 가면서 약화되었으며, 이는 분쇄현상에 의해 유효 열전달 면적의 확대 효과가 저감되었기 때문으로 해석된다.

• Journal of the Korean Physical Society
• /
• 제73권10호
• /
• pp.1584-1588
• /
• 2018

During hyperthermia therapy, cancer cells are heated to a temperature in the range of $40{\sim}45^{\circ}C$ for a defined time period to damage these cells while keeping healthy tissues at safe temperatures. Prior to hyperthermia therapy, the amount of heat energy transferred to the cancer cells must be predicted. Among various non-invasive methods, the thermal prediction method using the specific absorption rate (SAR) is the most widely used method. The existing methods predict the thermal distribution by using a single constant for the mass density in one organ through assignment. However, because the SAR and the bio heat equation (BHE) vary with the mass density, the mass density of each organ must be accurately considered. In this study, the mass density distribution was calculated using the relationship between the Hounsfield unit and the mass density of tissues in preceding research. The SAR distribution was found using a quasi-static approximation to Maxwell's equation and was used to calculate the potential distribution and the energy distributions for capacitive RF heating. The thermal distribution during exposure to RF waves was determined by solving the BHE with consideration given to the considering contributions of heat conduction and external heating. Compared with reference data for the mass density, our results was within 1%. When the reconstructed temperature distribution was compared to the measured temperature distribution, the difference was within 3%. In this study, the density distribution and the thermal distribution were reconstructed for the agar phantom. Based on these data, we developed an algorithm that could be applied to patients.

• Tribology and Lubricants
• /
• 제30권5호
• /
• pp.303-310
• /
• 2014

In order to reduce friction and improve reliability, researchers have applied various surface texturing methods to highly sliding machine elements such as mechanical seals and piston rings. Despite extensive theoretical research on surface texturing, previous numerical results are only applicable to isothermal and iso-viscous conditions. Because the lubricant flow pattern of textured bearing surfaces is much more complicated than that for non-textured bearings, the Navier?Stokes equation is more suitable than the Reynolds equation for the former. This study carries out a thermohydrodynamic (THD) lubrication analysis to investigate the lubrication characteristics of a single micro-dimpled parallel thrust bearing cell. The analysis involves using the continuity, Navier?Stokes, energy, temperature?viscosity relation, and heat conduction equations with the commercial computational fluid dynamics (CFD) code FLUENT. This study discretizes these equations using the finite volume method and solves them using the SIMPLE algorithm. The results include finding the streamlines, pressure and temperature distributions, and variations in the friction force and leakage for various dimple radii and depths. Increasing the dimple radius and decreasing the depth causes a recirculation flow to form because of a strong vortex, and the oil temperature greatly increases compared with the non-textured case. The present numerical scheme and results are applicable to THD analysis of various surface-textured sliding bearings and can lead to further study.