• Nuclear Engineering and Technology
• /
• v.56 no.3
• /
• pp.772-784
• /
• 2024

The hexagonal nodal code RENUS has been enhanced to handle irregularly deformed hexagonal assemblies. The underlying RENUS methods involving triangle-based polynomial expansion nodal (T-PEN) and corner point balance (CPB) were extended in a way to use line and surface integrals of polynomials in a deformed hexagonal geometry. The nodal calculation is accelerated by the coarse mesh finite difference (CMFD) formulation extended to unstructured geometry. The accuracy of the unstructured nodal solution was evaluated for a group of 2D SFR core problems in which the assembly corner points are arbitrarily displaced. The RENUS results for the change in nuclear characteristics resulting from fuel deformation were compared with those of the reference McCARD Monte Carlo code. It turned out that the two solutions agree within 18 pcm in reactivity change and 0.46% in assembly power distribution change. These results demonstrate that the proposed unstructured nodal method can accurately model heterogeneous thermal expansion in hexagonal fueled cores.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.12
• /
• pp.1926-1932
• /
• 2001

In this second part of the paper, the automatic mesh generation and remeshing algorithm using bubble packing method is applied to the nonlinear problem. The remeshing/refinement procedure is necessary in the large deformation process especially because the mesh distortion deteriorates the convergence and accuracy. To perform the nonliear analysis, the transfer of state variables such as displacement and strain is added to the algorithm of Part 1. The equilibrium equation based on total Lagrangian formulation and elasto-viscoplastic model is used. For the numerical experiment, the upsetting process including the contact constraint condition is analyzed by two refinement criteria. And from the result, it is addressed that the present algorithm can generate the refined meshes easily at the largely deformed area with high error.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2003.05a
• /
• pp.323-326
• /
• 2003

This paper presents a new algorithm that can refine hexahedral elements while maintaining the appropriate connectivity. In the algorithm, at first the regions of mesh to be refined are defined and, then, the zero-thickness element layers inserted into the interfaces between the regions. All the meshes in the regions, in which the zero-thickness layers are inserted, are to be regularized in order to improve the shape of the slender elements on the interfaces. This algorithm is applied to the analysis of plastic deformation process. The results show that the refined mesh gives smaller relative errors than the original mesh.

• Transactions of Materials Processing
• /
• v.12 no.4
• /
• pp.334-339
• /
• 2003

This paper presents a new algorithm that enables the refinement of hexahedral elements while maintaining the appropriate connectivity. In the algorithm, at first the regions of mesh to be refined are defined and, then, the zero-thickness element layers are inserted into the interfaces between the regions. All the meshes in the regions, in which the zero-thickness layers are inserted, are to be regularized in order to improve the shape of the slender elements on the interfaces. This algorithm is applied to the analysis of plastic deformation process. The results show that the refined mesh gives smaller relative errors than the original mesh.

• Structural Engineering and Mechanics
• /
• v.25 no.5
• /
• pp.501-518
• /
• 2007

In this work, a new approach is developed for dynamic analysis of a composite beam with an interply crack, based on finite element solution of partial differential equations with the use of the COMSOL Multiphysics package, allowing for fast and simple change of geometric characteristics of the delaminated area. The use of COMSOL Multiphysics package facilitates automatic mesh generation, which is needed if the problem has to be solved many times with different crack lengths. In the model, a physically impossible interpenetration of the crack faces is prevented by imposing a special constraint, leading to taking account of a force of contact interaction of the crack faces and to nonlinearity of the formulated boundary value problem. The model is based on the first order shear deformation theory, i.e., the longitudinal displacement is assumed to vary linearly through the beam's thickness. The shear deformation and rotary inertia terms are included into the formulation, to achieve better accuracy. Nonlinear partial differential equations of motion with boundary conditions are developed and written in the format acceptable by the COMSOL Multiphysics package. An example problem of a clamped-free beam with a piezoelectric actuator is considered, and its finite element solution is obtained. A noticeable difference of forced vibrations of the delaminated and undelaminated beams due to the contact interaction of the crack's faces is predicted by the developed model.

• Journal of computational fluids engineering
• /
• v.18 no.3
• /
• pp.26-33
• /
• 2013

In the present study, unsteady flow simulations of a variable geometry nozzle were conducted using a two-dimensional flow solver based on hybrid unstructured meshes. The variable geometry nozzle is used to achieve efficient performances of aircraft engines at various operating conditions. To describe the motion of the variable geometry nozzle, an algebraic method based on the basis decomposition of normal edge vector was used for the deformation of viscous elements. A ball-vertex spring analogy was used for inviscid elements. The aerodynamic data were obtained for a range of nozzle pressure ratios, and the validations were made by comparing the present results with available experimental data. The unsteady nozzle flows were simulated with an oscillating diverging section and a converging-diverging section. It was found that the nozzle performances are influenced by the nozzle exit flow characteristics, mass flow rate, as well as unsteady effects. These unsteady effects are shown to behave differently depending on the frequency of the nozzle motion.

• Transactions of Materials Processing
• /
• v.15 no.6 s.87
• /
• pp.459-466
• /
• 2006

Trimming line design plays an important role in obtaining accurate edge profile after flanging. Compared to the traditional section-based method, simulation-based method can produce more accurate trimming line by considering deformation mechanics. Recently, the use of a finite element inverse method is proposed to obtain optimal trimming line. By analyzing flanging inversely from the final mesh after flanging, trimming line can be obtained from initial mesh on the drawing die surface. Initial guess generation fer finite element inverse method is obtained by developing the final mesh onto drawing tool mesh. Incremental development method is adopted to handle irregular mesh with various size and undercut. In this study, improved incremental development algorithm to handle complex shape is suggested. When developing the final mesh layer by layer, the algorithm which can define the development sequence and the position of developing nodes is thoroughly described. Flanging of front fender is analyzed to demonstrate the effectiveness of the present method. By using section-based trimming line and simulation-based trimming line, incremental finite element simulations are carried out. In comparison with experiment, it is clearly shown that the present method yields more accurate edge profile than section-based method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.9 s.180
• /
• pp.2397-2405
• /
• 2000

The remeshing is a method to replace a distorted mesh by a new mesh without interrupting the finite element calculation. The remeshing procedure in this paper refers to the rezoning, for which a sm oothing process is developed to alleviate the distortions of mesh. In the paper, an automatic finite element rezoning system with tetrahedral elements for large deformation analysis has been developed. Our smoothing process is composed of two steps, a surface smoothing and a volume smoothing. In the surface smoothing, checking the dihedral angle and projection on surface patch reduced the change of shape and nodes penetrating die. The constrained Laplacian smoothing has been employed for the volume smoothing process. The state variables are mapped from old mesh to new mesh by using volume coordinates within a tetrahedral element. All these procedures have been linked to the NIKE3D program As illustrated in the examples the overall strategy ensures a robust and efficient rezoning scheme for finite element simulation of metal-forming processes.

• Korean Journal of Optics and Photonics
• /
• v.17 no.4
• /
• pp.359-365
• /
• 2006

Wave-front error(WFE) is the main parameter that determines the optical performance of the opto-mechanical system. In the development of opto-mechanics, WFE due to the main loading conditions are set to the important specifications. The deformation of the optical surface can be exactly calculated thanks to the evolution of numerical methods such as the finite element method(FEM). To calculate WFE from the deformation results of FEM, another approximation of the optical surface deformation is required. It needs to construct additional grid or element mesh. To construct additional mesh is troublesomeand leads to transformation error. In this work, the moving least-squares approximation is used to reconstruct wave front error It has the advantage of accurate approximation with only nodal data. There is no need to construct additional mesh for approximation. The proposed method is applied to the examples of GOCI scan mirror in various loading conditions. The validity is demonstrated through examples.

• Korea-Australia Rheology Journal
• /
• v.15 no.2
• /
• pp.75-82
• /
• 2003

In this study, with the finite element method we numerically investigate the deformation of liquid drops surrounded by Newtonian or non-Newtonian viscous medium in the axisymmetric contraction flow. 1, 2 or 4 Newtonian or non-Newtonian drops are considered and the truncated power-law model is applied In order to describe non-Newtonian viscous behavior for both fluids. In this type of flow the drop exhibits considerably large deformation, and thus techniques of unstructured mesh generation and auto-remeshing are employed to accurately express the fluid mechanical behavior. We examine the deformation pattern of liquid drops with viscosity dependence different from that of the surrounding medium and also explain their interactions by comparing relative position or speed of drop front.