• Nuclear Engineering and Technology
• /
• v.54 no.8
• /
• pp.3059-3072
• /
• 2022

A Jacobian-Free Newton Krylov Two-Nodal Coarse Mesh Finite Difference algorithm based on Nodal Expansion Method (NEM_TNCMFD_JFNK) is successfully developed and proposed to solve the three-dimensional (3D) and multi-group reactor physics models. In the NEM_TNCMFD_JFNK method, the efficient JFNK method with the Modified Incomplete LU (MILU) preconditioner is integrated and applied into the discrete systems of the NEM-based two-node CMFD method by constructing the residual functions of only the nodal average fluxes and the eigenvalue. All the nonlinear corrective nodal coupling coefficients are updated on the basis of two-nodal NEM formulation including the discontinuity factor in every few newton steps. All the expansion coefficients and interface currents of the two-node NEM need not be chosen as the solution variables to evaluate the residual functions of the NEM_TNCMFD_JFNK method, therefore, the NEM_TNCMFD_JFNK method can greatly reduce the number of solution variables and the computational cost compared with the JFNK based on the conventional NEM. Finally the NEM_TNCMFD_JFNK code is developed and then analyzed by simulating the representative PWR MOX/UO₂ core benchmark, the popular NEACRP 3D core benchmark and the complicated full-core pin-by-pin homogenous core model. Numerical solutions show that the proposed NEM_TNCMFD_JFNK method with the MILU preconditioner has the good numerical accuracy and can obtain higher computational efficiency than the NEM-based two-node CMFD algorithm with the power method in the outer iteration and the Krylov method using the MILU preconditioner in the inner iteration, which indicates the NEM_TNCMFD_JFNK method can serve as a potential and efficient numerical tool for reactor neutron diffusion analysis module in the JFNK-based multiphysics coupling application.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.38C no.10
• /
• pp.920-927
• /
• 2013

Chroma key technique which composes images by separating an object from its background in specific color has restrictions on color and space. Especially, unlike general chroma key technique, image composition for stereo 3D display requires natural image composition method in 3D space. The thesis attempted to compose images in 3D space using depth keying method which uses high resolution depth information. High resolution depth map was obtained through camera calibration between the DSLR and Kinect sensor. 3D mesh model was created by the high resolution depth information and mapped with RGB color value. Object was converted into point cloud type in 3D space after separating it from its background according to depth information. The image in which 3D virtual background and object are composed obtained and played stereo 3D images using a virtual camera.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.1 s.173
• /
• pp.144-155
• /
• 2000

초록 The validity, of a single parameter such as stress intensity, factor K or J-integral in traditional fracture mechanics depends strongly on the geometry, and loading condition. Therefore the second parameter like T-stress measuring the stress constraint is additionally needed to characterize the general crack-tip fields. While many, research works have been done to verify, the J-T description of elastic-plastic crack-tip stress fields in plane strain specimens, limited works (especially. for bimaterials) have been performed to describe the structural surface crack-front stress fields with the two parameters. On this background, via detailed three dimensional finite element analyses for surface-cracked plates and straight pipes of homogeneous materials and bimaterials under various loadings, we investigate the extended validity or limitation of the two parameter approach. We here first develop a full 3D mesh generating program for semi-elliptical surface cracks, and calculate elastic T-stress from the obtained finite element stress field. Comparing the J-T predictions to the elastic-plastic stresses from 3D finite element analyses. we then confirm the extended validity of fracture mechanics methodology based on the J-T two parameters in characterizing the surface crack-front fields of welded plates and pipes under various loadings.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.9
• /
• pp.1325-1333
• /
• 2001

Critical defects in pressure vessels and pipes are generally found in the form of a semi-elliptical surface crack, and the analysis of which is consequently an important problem in engineering fracture mechanics. Furthermore, in addition to the traditional single parameter K or J-integral, the second parameter like T-stress should be measured to quantify the constraint effect. In this work, the validity of the line-spring finite element is investigated by comparing line-spring J-T solutions to the reference 3D finite element J-T solutions. A full 3D-mesh generating program for semi-elliptical surface cracks is employed to provide such reference 3D solutions. Then some structural characteristics of the surface-cracked T and L-joints are studied by mixed mode line-spring finite element. Negative T-stresses observed in T and L-joints indicate the necessity of J-T two parameter approach for analyses of surface-cracked T and L-joints.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.64-67
• /
• 2007

In the present study, the wind pressure transients on platform screen door in side platforms caused by passing trains have been investigated numerically. The transient compressible 3-D full Navier-Stokes solution is obtained with actual operational condition of subway train and the moving mesh technique is adopted considering the train movement. To achieve more accurate analysis, the entrance and exit tunnel of platform are included in a computational domain and detailed shape of train is also modeled Numerical analyses were conducted on five operational condition which are different velocity variation of subway train, existence of stationary train and passing each other trains. The results show that pressure load on platform screen door is maximized when the two trains are passing each other. It is also seen from the computational results that the maximum pressure variation was found to be satisfactory to various foreign standards.

• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.1-6
• /
• 2007

In the present study, the wind pressure transients on platform screen door in side platforms caused by passing trains have been investigated numerically. The transient compressible 3-D full Navier-Stokes solution is used with actual operational condition of subway train by adopting the moving mesh technique considering the train movement. To achieve more accurate analysis, the entrance and exit tunnel connecting the stations are included in a computational domain with modeling the detailed shape of the train. Numerical analyses are conducted on five operational conditions which include the variation of the train speed, case with or without the train stopped in the other track, and case for two trains passing each other inside the station. The results show that pressure load on platform screen door is maximized when the two trains are passing each other. It is also seen from the computational results that the maximum pressure variation for the cases considered in the present study is found to be satisfactory to various foreign standards.

• Wind and Structures
• /
• v.29 no.6
• /
• pp.405-416
• /
• 2019

Aerodynamic characteristic of a small scale wind turbine under the influence of an incoming uniform wind field is studied using k-ω Shear Stress Transport turbulence model. Firstly, the lift and drag characteristics of the blade section consisting of S826 airfoil is studied using 2D simulations at a Reynolds number of 1×10⁵. After that, the full turbine including the rotational effects of the blade is simulated using Multiple Reference Frames (MRF) and Sliding Mesh Interface (SMI) numerical techniques. The differences between the two techniques are quantified. It is then followed by a detailed comparison of the turbine's power/thrust output and the associated wake development at three tip speeds ratios (λ = 3, 6, 10). The phenomenon of blockage effect and spatial features of the flow are explained and linked to the turbines power output. Validation of wake profiles patterns at multiple locations downstream is also performed at each λ. The present work aims to evaluate the potential of the numerical methods in reproducing wind tunnel experimental results such that the method can be applied to full-scale turbines operating under realistic conditions in which observation data is scarce or lacking.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.224-229
• /
• 2004

The objective of this work is not only to perform feasibility studies on the CFD (computational fluid dynamics) analysis for the capillary system design but also to provide an enhanced understanding of the autonomous capillary flow. The capillary flow is evaluated by means of the commercial CFD software of FLUENT, which includes the VOF (volume-of-fluid) model for multiphase flow analysis. The effect of wall adhesion at fluid interfaces in contact with rigid boundaries is considered in terms of static contact angle. Feasibility studies are first performed, including mesh-resolution influence on pressure profile, which has a sudden increase at the liquid/gas interface. Then we perform both 2D and 3D simulations and examine the transient nature of the capillary flow. Analytical solutions are also derived for simple cases and compared with numerical results. Through this work, essential information on the capillary system design is brought out. Our efforts and initial success in numerical description of the microfluidic capillary flows enhance the fundamental understanding of the autonomous capillary flow and will eventually pave the road for full-scale, computer-aided design of microfluidic networks.

• Transactions of Materials Processing
• /
• v.9 no.3
• /
• pp.233-241
• /
• 2000

It is well known that the quality and efficiency of the design of metal forming processes can be significantly improved with the aid of effective numerical simulations. In the present study, a two-and three-dimensional finite element simulation system, CAMP form, was developed for the analysis of metal forming processes in the PC environment. It is composed of a solver based on the thermo-rigid-viscoplastic approach and graphic user interface (GUI) based pre-and post-processors to be used for the effective description of forming conditions and graphic display of simulation results, respectively. In particular, in the case of CAMPform 2D (two-dimensional), as the solver contains an automatic remeshing module which determines the deformation step when remeshing is required and reconstructs the new mesh system, it is possible to carry out simulations automatically without any user intervention. Also, the forming analysis considers ductile fracture of the workpiece and wear of dies for better usage of the system. In the case of CAMPform 3D, general three-dimensional problems that involve complex die geometries and require remeshing can be analyzed, but full automation of simulations has yet to be achieved. In this paper, the overall structure and computational background of CAMPform will be briefly explained and analysis results of several forming processes will be shown. From the current results, it is construed that CAMPform can be used in providing useful information to assist the design of forming processes.

• Nuclear Engineering and Technology
• /
• v.55 no.4
• /
• pp.1428-1438
• /
• 2023

A novel Fission Diffusion Synthetic Acceleration (FDSA) method is developed and implemented as a part of a hybrid neutronics method for source convergence acceleration and variance reduction in Monte Carlo (MC) criticality calculations. The acceleration of the MC calculation stems from constructing a synthetic operator and solving a low-order problem using information obtained from previous MC calculations. By applying the P₁ approximation, two correction terms, one for the scalar flux and the other for the current, can be solved in the low-order problem and applied to the transport solution. A variety of one-dimensional (1-D) and two-dimensional (2-D) numerical tests are constructed to demonstrate the performance of FDSA in comparison with the standalone MC method and the coupled MC and Coarse Mesh Finite Difference (MC-CMFD) method on both intended purposes. The comparison results show that the acceleration by a factor of 3-10 can be expected for source convergence and the reduction in MC variance is comparable to CMFD in both slab and full core geometries, although the effectiveness of such hybrid methods is limited to systems with small dominance ratios.