• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.258-270
• /
• 2020

A new task of using the Jacobian-Free-Newton-Krylov (JFNK) method for the PWR core transient simulations involving neutronics, thermal hydraulics and mechanics is conducted. For the transient scenario of PWR, normally the Picard iteration of the coupled coarse-mesh nodal equations and parallel channel TH equations is performed to get the transient solution. In order to solve the coupled equations faster and more stable, the Newton Krylov (NK) method based on the explicit matrix was studied. However, the NK method is hard to be extended to the cases with more physics phenomenon coupled, thus the JFNK based iteration scheme is developed for the nodal method and parallel-channel TH method. The local gap conductance is sensitive to the gap width and will influence the temperature distribution in the fuel rod significantly. To further consider the local gap conductance during the transient scenario, a 1D mechanics model is coupled into the JFNK scheme to account for the fuel thermal expansion effect. To improve the efficiency, the physics-based precondition and scaling technique are developed for the JFNK iteration. Numerical tests show good convergence behavior of the iterations and demonstrate the influence of the fuel thermal expansion effect during the rod ejection problems.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1651-1656
• /
• 2003

The present study has been conducted to analytically investigate the thermal control performance of variable conductance heat pipe(YCHP) with axial grooves. The condenser port of the YCHP is occupied by a inert gas in which the concentration of gas is varied with the operation temperature and the heat transport capacity is thus varied with the operating temperature due to the variation of inert gas concentration. In this study, numerical evaluation for the thermal control of the YCHP with axial grooves is made from the 1st order diffusion model that considers the diffusive expansion of inert gas by concentration gradient. Ammonia is used as a working fluid and Nitrogen as a control gas in the Aluminum tube. As a result, the thermal performance of YCHP based on diffusion model has been compared with that of YCHP from flat front model. Additionally, it is found that the concentration of inert gas is distributed in the condenser region of YCHP with axial grooves.

• Nuclear Engineering and Technology
• /
• v.54 no.6
• /
• pp.2094-2106
• /
• 2022

The solid-state core of a heat pipe cooled reactor operates at high temperatures over 1000 K with thermal and irradiation-induced expansion during burnup. The expansion changes the gap thickness between the solid components and the material properties, and may even cause the gap closure, which then significantly influences the thermal and mechanical characteristics of the reactor core. This study developed an irradiation behavior model for HPRTRAN, a heat pipe reactor system analysis code, to introduce the irradiation effects such as swelling and creep. The megawatt heat pipe reactor MegaPower was chosen as an application case. The coupled irradiation-thermal-mechanical model was developed to simulate the irradiation effects on the heat transfer and stresses of the whole reactor core. The results show that the irradiation deformation effect is significant, with the irradiation-induced strains up to 2.82% for fuel and 0.30% for monolith at the end of the reactor lifetime. The peak temperatures during the lifetime are 1027:3 K for the fuel and 956:2 K for monolith. The gap closure enhances the heat transfer but caused high stresses exceeding the yield strength in the monolith.

• Korean Journal of Materials Research
• /
• v.25 no.12
• /
• pp.678-682
• /
• 2015

The electrical property of polymer matrix composites with added carbon powder is studied based on the temperature dependency of the conduction mechanism. The temperature coefficient of the resistance of the polymer matrix composites below the percolation threshold (x) changed from negative to positive at 0.20 < x < 0.21; this trend decreased with increasing of the percolation threshold. The temperature dependence of the electrical property(resistivity) of the polymer matrix composites below the percolation threshold can be explained by using a tunneling conduction model that incorporates the effect of the thermal expansion of the polymer matrix composites into the tunneling gap. The temperature coefficient of the resistance of the polymer matrix composites above the percolation threshold has a positive value; its absolute value increased with increasing volume fraction of carbon powder. By assuming that the electrical conduction through the percolating paths is a thermally activated process and by incorporating the effect of thermal expansion into the volume fraction of the carbon power, the temperature dependency of the resistivity above the percolation threshold can be well explained without violating the universal law of conductivity.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2010.11a
• /
• pp.359-360
• /
• 2010

• Journal of the Korean Society for Precision Engineering
• /
• v.17 no.6
• /
• pp.199-209
• /
• 2000

Particles reinforced MMCs have higher specific modulus, higher specific strength, better properties at elevated temperatures and better wear resistance than monolithic metals. But the coefficient of thermal expansion(CTE) of Al6061 is 5 times larger than that of SiCp. The discrepancy of CTE makes some residual stresses inside of MMCs. This work investigates Si$C_p$/Al6061 composites at high temperatures in the microscopic view by three-dimensional elasto-plastic finite element analyses and compares the analytical results with the experimental ones. The theoretical model is not able to consider the nonuniform shape of particle. So the shape of particle is assumed to be perfect global shape. And also particle distribution is not homogeneous in experimental specimen. It is assumed to be homogeneous in simulation model. The type of particle distribution is face-centered cubic array(FCC array). Furthermore, non-homogeneous distribution is modeled by combination of several volume fractions.

• The Bulletin of The Korean Astronomical Society
• /
• v.36 no.1
• /
• pp.54.1-54.1
• /
• 2011

Formation of self-gravitating gas clouds and hence stars in galaxies is a consequence of both thermal and dynamical evolution of a gaseous medium. Using hydrodynamics simulations including cooling and heating explicitly, we follow simultaneously thermal and dynamical evolution of galactic gas disks to study dynamics and structures of galactic spiral shocks with thermal instability and regulation of the star formation rates (SFRs). We first perform one-dimensional simulations in direction perpendicular to spiral arms. The multiphase gas flows across the arm soon achieve a quasi-steady state characterized by transitions from warm to cold phases at the shock and from cold to warm phases in the postshock expansion zone, producing a substantial fraction of intermediate-temperature gas. Next, we allow a vertical degree of freedom to model vertically stratified disks. The shock front experiences unsteady flapping motions, driving a significant amount of random gas motions, and self-gravity promotes formation of bound clouds inside spiral arms. Finally, we include the star formation feedback in both mechanical (due to supernova explosion) and radiative (due to FUV heating by young stars) forms in the absence of spiral arms. At saturation, gravitationally bound clouds form via thermal and gravitational instabilities, which are compensated by disruption via supernova explosions. We find that the FUV heating regulates the SFRs when gas surface density is low, confirming the prediction of the thermal and dynamical equilibrium model of Ostriker et al. (2010) for star formation regulation.

• Journal of the Korea Concrete Institute
• /
• v.23 no.5
• /
• pp.667-674
• /
• 2011

This paper discusses the analysis method of prestressed concrete girder integral abutment bridges for a 75-year bridge life and the development of prediction models for abutment displacements under thermal loading due to annual temperature fluctuation and time-dependent loading. The developed nonlinear numerical modeling methodologies considered soil-structure interaction between supporting piles and surrounding soils and between abutment and backfills. Material nonlinearity was also considered to simulate differential rotation in construction joints between abutment and backwall. Based on the numerical modeling methodologies, a parametric study of 243 analysis cases, considering five parameters: (1) thermal expansion coefficient, (2) bridge length, (3) backfill height, (4) backfill stiffness, and (5) pile soil stiffness, was performed to established prediction models for abutment displacements over a bridge life. The parametric study results revealed that thermal expansion coefficient, bridge length, and pile-soil stiffness significantly influenced the abutment displacement. Bridge length parameter significantly influenced the abutment top displacement at the centroid of the superstructure, which is similar to the free expansion analysis results. Developed prediction model can be used for a preliminary design of integral abutment bridges.

• Journal of the Korea Convergence Society
• /
• v.11 no.8
• /
• pp.153-158
• /
• 2020

In this study, the thermal and stress analyses were performed by applying a temperature condition of 100℃ at CPU cooler model. The value of heat flux value is shown to be the most at the lower rod area. The upper part becomes, the smaller the heat flow rate. The highest temperature is shown at the bottom of the CPU cooler model. Overall, the upper part becomes, the smaller the temperature becomes. Based on the temperature analysis, the thermal deformation caused by expansion, the deformation becomes smaller as the upper part of the overlapping plates. The great deformation happens at the bent area of the small rod as the lower part of model and the least deformation is shown at the lowest floor of model. In addition, the maximum thermal stress of 570.63 MPa happens at the floor below model. The stress is shown to decrease as the upper part of the overlapping plates becomes. But the stress is shown to increase somewhat at the middle part of model. By applying the study result on the thermal and stress analyses of CPU cooler, this study is seen to be suitable for the aesthetic convergence.

• Journal of the Korean Society for Precision Engineering
• /
• v.28 no.1
• /
• pp.96-101
• /
• 2011

Under extreme test or operation condition, refrigerator endures complicated stresses state and thermal bowing deformation arises on the sidewall. Shelf rails designed in the inner case provide increased surface area to permit expansion without bowing, and also increase structural rigidity to resist bowing. In this study, we designed six different shelf patterns of refrigerator model and studied the bead on refrigerator deformation using finite element method (FEM). Analysis result shows that increasing the numbers of beads properly in refrigerator is more helpful to reduce thermal bowing deformation. In addition, the beads would decrease stress on refrigerator sidewall. However, vertical beads have no effect to reduce thermal deformation of the bowing.