• Nuclear Engineering and Technology
• /
• v.52 no.11
• /
• pp.2431-2441
• /
• 2020

Currently, the pin-by-pin homogenized solvers are a very active research field as they can, unlike the nodal codes, directly predict the local power, while requiring significantly less computational resources than the heterogeneous transport codes. This paper presents a recently developed pin-by-pin diffusion/SP3 solver Tortin, its spatial discretization method and the reflector treatment. Regarding the spatial discretization, it was observed that the finite difference method applied on pin-cell size mesh does not properly capture the big flux change between MOX and uranium fuel, while the nodal expansion method is more accurate but too slow. If the finite difference method is used with a finer mesh in the outer two pin rows of the fuel assembly, it increases the required computation time by only 50%, but decreases the pin power errors below 1% with respect to lattice code reference solutions. The paper further describes the coupling of Tortin with a microscopic depletion solver. Several verification tests show that the SP3 pin-by-pin solver can reproduce the heterogeneous transport solvers results with very good accuracy, even for fuel cycle depletion of very heterogeneous core employing MOX fuel or inserted control rods, while being two orders of magnitude faster.

• Journal of the Korean Wood Science and Technology
• /
• v.22 no.2
• /
• pp.61-70
• /
• 1994

This study was executed to find the applicability of press drying of tree disk by investigating the shrinkage and drying defect and to form appropriate model by comparing the actual moisture content(MC) and internal temperature in respect of drying time with calculated values based computer simulation to which was applied finite difference method. In press drying disk, heating period, constant drying rate period maintained plateau temperature at 100$^{\circ}C$ and falling drying rate period were significantly distinguished. Actual MC and internal temperature were analogous to those calculated at comparing points. Heat transfer model formed by Fourier's law using specific heat of moist wood and conduction coefficient considering fractional volume of each element of wood cell wall, bound water, free water and air showed applicability as basic data to developing heat expansion, shrinkage and drying stress during press drying. Also mass transfer model formed by Fick's diffusion law using water vapor diffusion coefficient showed applicability. Longitudinal shrinkage was developed by pressure of hot press and tangential shrinkage was restrained by hygrothermal recovery. The heart check, surface check and ring failure were occurred differently in species, but V-shaped crack didn't develop.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.2
• /
• pp.139-148
• /
• 2012

The MLS(Moving Least Squares) Difference Method is a numerical scheme that combines the MLS method of Meshfree method and Taylor expansion involving not numerical quadrature or mesh structure but only nodes. This paper presents an dynamic algorithm of MLS difference method for solving transient solid mechanics problems. The developed algorithm performs time integration by using Newmark method and directly discretizes strong forms. It is very convenient to increase the order of Taylor polynomial because derivative approximations are obtained by the Taylor series expanded by MLS method without real differentiation. The accuracy and efficiency of the dynamic algorithm are verified through numerical experiments. Numerical results converge very well to the closed-form solutions and show less oscillation and periodic error than FEM(Finite Element Method).

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.4
• /
• pp.493-499
• /
• 2007

The Taylor expansion expresses a differentiable function and its coefficients provide good approximations for the given function and its derivatives. In this study, m-th order Taylor Polynomial is constructed and the coefficients are computed by the Moving Least Squares method. The coefficients are applied to the governing partial differential equation for solid problems including crack problems. The discrete system of difference equations are set up based on the concept of point collocation. The developed method effectively overcomes the shortcomings of the finite difference method which is dependent of the grid structure and has no approximation function, and the Galerkin-based meshfree method which involves time-consuming integration of weak form and differentiation of the shape function and cumbersome treatment of essential boundary.

• Structural Engineering and Mechanics
• /
• v.7 no.2
• /
• pp.181-202
• /
• 1999

Dynamic response of axisymmetric arbitrary laminated composite cylindrical shell of finite length, using three-dimensional elasticity equations are studied. The shell is simply supported at both ends. The highly coupled partial differential equations are reduced to ordinary differential equations (ODE) with variable coefficients by means of trigonometric function expansion in axial direction. For cylindrical shell under dynamic load, the resulting differential equations are solved by Galerkin finite element method, In this solution, the continuity conditions between any two layer is satisfied. It is found that the difference between elasticity solution (ES) and higher order shear deformation theory (HSD) become higher for a symmetric laminations than their unsymmetric counterpart. That is due to the effect of bending-streching coupling. It is also found that due to the discontinuity of inplane stresses at the interface of the laminate, the slope of transverse normal and shear stresses aren't continuous across the interface. For free vibration analysis, through dividing each layer into thin laminas, the variable coefficients in ODE become constants and the resulting equations can be solved exactly. It is shown that the natural frequency of symmetric angle-ply are generally higher than their antisymmetric counterpart. Also the results are in good agreement with similar results found in literatures.

• Structural Engineering and Mechanics
• /
• v.44 no.1
• /
• pp.109-125
• /
• 2012

Post-buckling behavior of Timoshenko beams subjected to uniform temperature rising with temperature dependent physical properties are studied in this paper by using the total Lagrangian Timoshenko beam element approximation. The beam is clamped at both ends. In the case of beams with immovable ends, temperature rise causes compressible forces end therefore buckling and post-buckling phenomena occurs. It is known that post-buckling problems are geometrically nonlinear problems. Also, the material properties (Young's modulus, coefficient of thermal expansion, yield stress) are temperature dependent: That is the coefficients of the governing equations are not constant in this study. This situation suggests the physical nonlinearity of the problem. Hence, the considered problem is both geometrically and physically nonlinear. The considered highly non-linear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. The beams considered in numerical examples are made of Austenitic Stainless Steel (316). The convergence studies are made. In this study, the difference between temperature dependent and independent physical properties are investigated in detail in post-buckling case. The relationships between deflections, thermal post-buckling configuration, critical buckling temperature, maximum stresses of the beams and temperature rising are illustrated in detail in post-buckling case.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.43 no.1
• /
• pp.1-7
• /
• 2015

This paper predicts the spring-in effect of curved laminated composite structure for various stacking sequence using finite element analysis(ABAQUS). In composite manufacturing process, large temperature difference, different coefficient of thermal expansion and chemical shrinkage effect cause distortion of composite parts such as spring-in and warpage. Distortion of composite structure is important issue on quality of product, and it should be considered in manufacturing process. In finite element analysis, a CHILE(Cure Hardening Instantaneously Linear Elastic) model and chemical shrinkage effects are considered developing user subroutine in ABAQUS and some cases are simulated.

• Journal of Ocean Engineering and Technology
• /
• v.10 no.2
• /
• pp.106-119
• /
• 1996

The ceramic has various high mechanical properties such as heat, abrasion, corrosion resistance and high temperature strength compared with metal. It also has low speciffic weight, low thermal expansibillity, low thermal conductivity. However, it could not be used as structural material since it is brittle and difficult for the machining. Therefore, there have been many researches to attempt to join ceramic with metal which is full of ductillity in order to compensate the weakness of ceramic.The problem is that residual stress develops around the joint area while the ceramic/metal joint material is cooled from high joining temperature to room temperature due to remarkable difference of thermal expansion coefficients between ceramic and metal. Especially, the residual stress at both edges of the specimen reduces the strngth of joint to a large amount by forming a singular stress field. In this study, two dimensional finite element method is attempted for the thermal elastic analysis. The joint residual stress of ceramic/metal developed in the cooling process is investigated and the change of joint residual stress resulted from the repetitive heat cycle is also examined. In addition, it is attempted to clarify the joint stress distribution of the case of tensile load and of the case of superposition of residual stress and actual loading stress.

• Journal of computational fluids engineering
• /
• v.7 no.3
• /
• pp.53-59
• /
• 2002

A two-dimensional laminar flow through a channel with a pair of symmetric vertical fins is investigated. At far up- and down-stream from the fins, the plane Poiseuille flow exists in the channel. The Stokes flow for this channel is first investigated analytically and then the other laminar flows by numerical method. For analytic method, the method of eigen function expansion and collocation method are employed. In numerical solution for laminar flows, finite difference method(FDM) is used to obtain vorticity and stream function. From the results, the streamline patterns are shown and the additional pressure drop due to the attached fins and the force exerted on the fin are calculated. It is clear that the force depends on the length of fins and Reynolds number. When the Reynolds number exceeds a critical value, the flow becomes asymmetric. This critical Reynolds number Re/sub c/ depends on the length of the fins.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2002.05a
• /
• pp.820-823
• /
• 2002

MEMS structures Generally have been fabricated using surface-machining, but the interface failure between silicon substrate and evaporated thin film frequently takes place due to difference of linear coefficient of thermal expansion. Therefore this paper studied the effect of the residual stress caused by variable external loads. This study did not analyzed accurate quantity of the residual stress but trend for the effect of residual stress. Several specimens were fabricated using other material(Al, Au and Cu) and thermal load was applied. The residual stress was measured by nano-indentation using AFM. The results showed the existence of the residual stress due to thermal load. The indentation area of the thermal loaded thin film reduced about 3.5% comparing with the virgin thin film caused by residual stress. The finite element analysis results are similar to indentation test.