• 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.

• Computers and Concrete
• /
• v.32 no.1
• /
• pp.61-74
• /
• 2023

This article investigates the wave propagation analysis of the imperfect functionally graded (FG) sandwich plates based on a novel simple four-variable integral quasi-3D higher-order shear deformation theory (HSDT). The thickness stretching effect is considered in the transverse displacement component. The presented formulation ensures a parabolic variation of the transverse shear stresses with zero-stresses at the top and the bottom surfaces without requiring any shear correction factors. The studied sandwich plates can be used in several sectors as areas of aircraft, construction, naval/marine, aerospace and wind energy systems, the sandwich structure is composed from three layers (two FG face sheets and isotropic core). The material properties in the FG faces sheet are computed according to a modified power law function with considering the porosity which may appear during the manufacturing process in the form of micro-voids in the layer body. The Hamilton principle is utilized to determine the four governing differential equations for wave propagation in FG plates which is reduced in terms of computation time and cost compared to the other conventional quasi-3D models. An eigenvalue equation is formulated for the analytical solution using a generalized displacements' solution form for wave propagation. The effects of porosity, temperature, moisture concentration, core thickness, and the material exponent on the plates' dispersion relations are examined by considering the thickness stretching influence.

• Current Research on Agriculture and Life Sciences
• /
• v.14
• /
• pp.49-60
• /
• 1996

In order to evaluate soil factors affecting the growth of Italian poplar, 23 areas planted with Italian poplar were surveyed. These 23 areas were classified into 3 categories, river-side, fallow-land and hill-side. The growth performance and soil factors for each area were investigated. The growth of Italian poplar at river-side was shown to be superior to that of fallow-land and fill-side. The rates of growth for fallow-land and hill-side are decreased by 8% and 21% compared to those of river-side, respectively. This suggests that plantation of Italian poplar at hill-side would not be profitable. Soil conditions of high productive area appeared liquid phase 20%, porosity 45%, water holding capacity 35 - 40%, soil hardness $1kg/cm^3$. pH 6 and rich in organic matter and total nitrogen. The results of factor analysis for soil factors affecting to Italian poplar growth that showed eigenvalue over 1 and communality value over 70% explained factor 1 : liquid phase, porosity and water holding capacity, factor 2 : pH and calcium, and factor 3 : soil hardness. This suggests that physical characteristics of soil is more important than chemical characteristics for Italian poplar growth. Multiregerssion analysis was conducted between diameter growth and soil hardness, liquid phase and calcium. The t-values for each independent variables showed significance at 1 - 10% level, but water holding capacity and pH are not significant. It is supposed that sites suitable to Italian poplar were alluvial plain of sandy loam or part of banking soil, well-ventilating soil, lower soil hardness, apposite soil moisture absorbing with about 100cm of ground water level, plentiful organic matters and total nitrogen and little acidity soil.

• Korean Journal of Remote Sensing
• /
• v.30 no.2
• /
• pp.259-274
• /
• 2014

On Aug 18, 2013, Sakurajima volcano in Japan erupted on a relatively large-scale. Geostationary Ocean Color Imager (GOCI) had used to detect volcanic ash in the surrounding area on the next day of this eruption. The geomagnetic variation has been analyzed using geomagnetic data from Cheongyang observatory in Korea and several geomagnetic observatories in Japan. First, we reconstruct geomagnetic data by principal component analysis and conduct semblance analysis by wavelet transform. Secondly, we minimize the error of solar effect by using wavelet based semblance filtering with Kp index. As a result of this study, we could confirm that the geomagnetic variation usually occur at the moment of Sakurajima volcano eruption. However, we cannot rule out the possibilities that it could have been impacted by other factors besides volcanic eruption in other variation's cases. This research is an exceptional study to analyze geomagnetic variation related with abroad volcanic eruption uncommonly in Korea. Moreover, we expect that it can help to develop further study of geomagnetic variation involved in earthquake and volcanic eruption.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.4
• /
• pp.419-429
• /
• 2003

A three-dimensional (3-D) method of analysis is presented for determining the free vibration frequencies of thick, hyperboloidal shells of revolution. Unlike conventional shell theories, which are mathematically two-dimensional (2-D), the present method is based upon the 3-D dynamic equations of elasticity. Displacement components u/sub r/, u/sub θ/, u/sub z/ in the radial, circumferential, and axial directions, respectively, we taken to be sinusoidal in time, periodic in θ, and algebraic polynomials in the r and z directions. Potential(strain) and kinetic energies of the hyperboloidal shells are formulated, and the Ritz method is used to solve the eigenvalue problem, thus yielding upper bound values of the frequencies by minimizing the frequencies. As the degree of the polynomials is increased, frequencies converge to the exact values. Convergence to four digit exactitude is demonstrated for the first five frequencies of the hyperboloidal shells of revolution. Numerical results are tabulated for eighteen configurations of completely free hyperboloidal shells of revolution having two different shell thickness ratios, three variant axis ratios, and three types of shell height ratios. Poisson's ratio (ν) is fixed at 0.3. Comparisons we made among the frequencies for these hyperboloidal shells and ones which ate cylindrical or nearly cylindrical( small meridional curvature. ) The method is applicable to thin hyperboloidal shells, as well as thick and very thick ones.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.6
• /
• pp.64-72
• /
• 2017

The longitudinal dynamic instability which can occur in the fueling process of a space launch vehicle is called pogo. It is caused by coupling between the fuselage and propulsion system and they would be formed as a closed-loop system. so that the amplitude of the response may increase or decrease. In this paper, a mathematical model which is applicable to the systematic pogo analysis of a general launch vehicle is developed for an example of space shuttle. The formulations are composed of the linearized second-order differential equation for the propulsion system, and of the pressure, weight displacement, and generalized displacement. Those are important parameters for pogo analysis, are derived through eigenvalue analysis. By the formulation suggested in this paper, it is expected that mathematical modeling method of the pogo system can be obtained and systematic pogo stability analysis for any launch vehicle will be enabled.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.23 no.5
• /
• pp.96-103
• /
• 2019

Middle-rise reinforced concrete residential buildings account for a large portion of the Korea, and structural performance analysis are needed for problems that could occur during the construction of such buildings. Thus, a middle-rise reinforced concrete residential building with 25 stories are selected as a sample model for structural performance analysis. The structural analyses are performed by dividing a sample model into the construction stage models of the 5th, 10th, 15th, 20th and 25th floors and the completion stage models with the design completed. For the comparisons of structural performances, Eigenvalue analysis results and lateral-load-resisting capabilities and structural design performances of structural members are analyzed. As a result of analyses, it was confirmed that both the construction and completion stage do not exceed KBC criteria limits at the lateral displacement and story drift ratio, and structural design performances of structural members confirm structural safety in all components except for some members of the wall. Therefore, it was concluded that if structural stability is obtained during the completion stage of a middle-rise reinforced concrete residential building, structural stability is secured under construction.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.30 no.2A
• /
• pp.81-92
• /
• 2010

In design of steel frames, it is generally believed that elastic system buckling analysis cannot predict real behaviors of structures, while inelastic system buckling analysis can give informative buckling behaviors of individual members considering inelastic material behavior. However, the use of Euler buckling equation with these system buckling analyses have the inherent problem that the methods evaluate unexpectedly large effective lengths of members having relatively small axial forces. This paper proposes a new method of obtaining elastic and inelastic effective lengths of all members in steel frames. Considering a fictitious axial force factor for each story of frames, the proposed method determines the effective lengths using the inelastic stiffness reduction factor and the iterative eigenvalue analysis. In order to verify the validity of the proposed method, the effective lengths of example frames by the proposed method were compared to those of previously established methods. As a result, the proposed method gives reasonable effective lengths of all members in steel frames. The effect of inelastic material behavior on the effective lengths of members was also discussed.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.197-206
• /
• 2003

A three-dimensional (3-D) method of analysis is presented for determining the free vibration frequencies and mode shapes of solid and hollow hemispherical shells of revolution of arbitrary wall thickness having arbitrary constraints on their boundaries. Unlike conventional shell theories, which are mathematically two-dimensional (2-D), the present method is based upon the 3-D dynamic equations of elasticity. Displacement components μ/sub Φ/, μ/sub z/, and μ/sub θ/ in the meridional, normal, and circumferential directions, respectively, are taken to be sinusoidal in time, periodic in θ, and algebraic polynomials in the Φ and z directions. Potential (strain) and kinetic energies of the hemispherical shells are formulated, and the Ritz method is used to solve the eigenvalue problem, thus yielding upper bound values of the frequencies obtained by minimizing the frequencies. As the degree of the polynomials is increased, frequencies converge to the exact values. Novel numerical results are presented for solid and hollow hemispheres with linear thickness variation. The effect on frequencies of a small axial conical hole is also discussed. Comparisons are made for the frequencies of completely free, thick hemispherical shells with uniform thickness from the present 3-D Ritz solutions and other 3-D finite element ones.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.3 no.4
• /
• pp.95-107
• /
• 1983

The general purpose programs are in their fixed algorithm and theory of mechanics which can not be altered without painful program modifications. Users are usually guided by user's manual for data input. The several symbolic manipulation programs for structural analysis are introduced recently. These programs allow users to include a wide class of solution algorithm and to specify, by means of some symbolic manipulation, a combination of analytical steps to suit a particular problem. As they can solve a single domain problem, a large computer is usually needed. The scope of this study is to develop an efficient symbolic manipulation program with space beam element, plate bending element and eigen value routines. The incorporated Substructure capability and generation capability of finite element characteristic arrays (e.g., stiffness matrix, mass matrix) enables users to analyse multidomain problem with small computer. The program consists of modulized independent processors, each having its own specific function and is easily modified, eliminated and added. The processors are efficiently handling data by the Data base approach which is the concept of integrated program network(IPN).