• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.345-350
• /
• 2000

Adaptive finite element analysis, which its solution error meets with the user defined allowable error, is recently used far improving reliability of finite element analysis results. This adaptive analysis is composed of two procedures; one is the error estimation of an analysis result and another is the reconstruction of finite elements. In the rp-method, an element size is controlled by relocating of nodal positions(r-method) and the order of an element shape function is determined by the hierarchical polynomial(p-method) corresponding to the element solution error. In order to show the effectiveness and accuracy of the suggested rp-method, various numerical examples were analyzed and these analysis results were examined by comparing with those obtained by the existed methods. As a result of this study, following conclusions are obtained. (1) rp-method is more accurate and effective than the r- and p-method. (2) The solution convergency of the rp-method is controlled by means of the iterative calculation numbers of the r- and p- method each other.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.12
• /
• pp.2298-2305
• /
• 1992

To obtain the convenient solution of the multibody dynamic systems composed of flexible beams, linear finite element technique is adopted and the nodal coordinates are interpolated in the global inertia frame. Mass matrix becomes an extremely simple constant matrix and the force vector also becomes extremely simple because Coriolis acceleration and centrifugal force are not required. And the elastic force is also simply computed from the moving frame attached to the material. To solve the global differential algebraic euation. an ODE technique is adopted after Lagrange multiplier is computed by the accelerated iterative technique, and the time demanding procedures such as Newton-Raphson iterations and decomposition of the big matrix are not required. The accuracy of the present solution is checked by a well-known example problem.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.1-10
• /
• 2004

Adaptive finite element analysis, in which its solution error meets with the user defined allowable error, is recently used to improve the reliability of finite element analysis results. This adaptive analysis is composed of two procedures; one is the error estimation of an analysis result and the other is the reconstruction of finite elements. In the (p-method, an element size is controlled by relocating of nodal positions (r-method) and the order of an element shape function is determined by the hierarchical polynomial (p-method) corresponding to the clement solution error by the enhanced SPR. In order to show the effectiveness and the accuracy of the suggested rp-method, various numerical examples were analyzed and these analysis results were examined by comparing with those obtained by the existed methods.

• Nuclear Engineering and Technology
• /
• v.45 no.6
• /
• pp.707-720
• /
• 2013

We analyze piecewise homogenization with flux-weighted cross sections and preservation of averaged currents at the boundary of the homogenized domain. Introduction of a set of flux discontinuity ratios (FDR) that preserve reference interface currents leads to preservation of averaged region reaction rates and fluxes. We consider the class of numerical discretizations with one degree of freedom per volume and per surface and prove that when the homogenization and computing meshes are equal there is a unique solution for the FDRs which exactly preserve interface currents. For diffusion submeshing we introduce a Jacobian-Free Newton-Krylov method and for all cases considered obtain an 'exact' numerical solution (eight digits for the interface currents). The homogenization is completed by extending the familiar full assembly homogenization via flux discontinuity factors to the sides of regions laying on the boundary of the piecewise homogenized domain. Finally, for the familiar nodal discretization we numerically find that the FDRs obtained with no submesh (nearly at no cost) can be effectively used for whole-core diffusion calculations with submesh. This is not the case, however, for cell-centered finite differences.

• Water for future
• /
• v.27 no.1
• /
• pp.111-121
• /
• 1994

The purpose of this paper is to construct a network model for the optimal allocation of limited water resources to the nodal system with given priorities. The solution technique for the model is based on the out-of-kilter algorithm(OKA). For the verification and application of the theoretical methodology and computer programs, the Geum river system is selected. Using release of Daecheong dam and water demand in Geum river basin, optimal allocation of water resources is accomplished for 4 cases(case 1 - case 4) which consider priority numbers in the demand nodes. The results of the application show that the model can reasonably represent the physical system, and water shortage at the demand nodes with high priority numbers is reduced. Its system solution was verified with that by the revised simplex algorithm.

• Nuclear Engineering and Technology
• /
• v.55 no.12
• /
• pp.4491-4503
• /
• 2023

This study aimed to verify and validate the transient simulation capability of the hybrid code system RAST-F for fast reactor analysis. For this purpose, control rod (CR) drop experiments involving eight separate CRs and six CR groups in the China Experimental Fast Reactor (CEFR) start-up tests were utilized to simulate the CR drop transient. The RAST-F numerical solution, including the neutron population, time-dependent reactivity, and CR worth, was compared against the measurement values obtained from two out-of-core detectors. Moreover, the time-dependent reactivity and CR worth from RAST-F were verified against the results obtained by the Monte Carlo code Serpent using continuous energy nuclear data. A code-to-code comparison between Serpent and RAST-F showed good agreement in terms of time-dependent reactivity and CR worth. The discrepancy was less than 160 pcm for reactivity and less than 110 pcm for CR worth. RAST-F solution was almost identical to the measurement data in terms of neutron population and reactivity. All the calculated CR worth results agreed with experimental results within two standard deviations of experimental uncertainty for all CRs and CR groups. This work demonstrates that the RAST-F code system can be a potential tool for analyzing time-dependent phenomena in fast reactors.

• Journal of Bio-Environment Control
• /
• v.14 no.4
• /
• pp.298-301
• /
• 2005

Potato (Solamum tuberosum 'Dejima') plantlets were investigated on culture type and initial quantity of inoculation in bioreactor and survival rate by hydroponics for mass production. rode stems (1 to 1.5cm in length) of potato plantlets multiplied in vitro were grown for 3 weeks in liquid Murashige and Skoog (MS) medium with sucrose $30 g\; L^{-1}$. When plantlets (80-node inoculation) were raised in 10L balloon type bubble (BB) bioreactor, the healthiest growth of plantlets was obtained from explants cultured in ebb & flow culture with medium supplied periodically 12 times per day. The suitable inoculation quantity of 20L BB bioreactor was 120 pieces of stem segments (mean 2.2g fresh weight) in ebb & flow culture. Number of nodal shoot was eight on the average. In controlled culture room, survival rate of plantlets at 7 days after stem cutting was above $70\%$ when they were acclimatized by hydroponics grown in deep flow and solid medium culture. The highest survival rate of the stem cutting plantlets was in nutrient solution adjusted to EC $1.4dS{\cdot}m^{-1}$. Stem cutting plantlets through one culture could be obtained $670\~900$, when plantlets were grown in ebb & flow culture during 3 weeks using a 20L bioreactor with initial 120 pieces of nodal segments. 11 is possible In do mass production of seedlings cultured in bioreactor and hydroponics.

• Journal of Korean Association for Spatial Structures
• /
• v.3 no.2 s.8
• /
• pp.85-90
• /
• 2003

If we use a third order approximation for the displacement function of beam element in finite element methods, finite element solutions of beams yield nodal displacement values matching to beam theory results to have no connection with the number increasing of elements of beams. It is assumed that, as the member displacement value at beam nodes are correct, the calculation procedure of beam element stiffness matrix have no numerical errors. A the member forces are calculated by the equations of $\frac{-M}{EI}=\frac{{d^2}{\omega}}{dx^2}\;and\;\frac{dM}{dx}=V$, the member forces at nodes of beams have errors in a moment and a shear magnitudes in the case of smaller number of element. The nodal displacement value of plate subject to the lateral load converge to the exact values according to the increase of the number of the element. So it is assumed that the procedures of plate element stiffness matrix calculations has a error in the fundamental assumptions. The beam methods for the high accuracy ratio solution Is also applied to the plate analysis. The method of reducing a error ratio of member forces and element stiffness matrix in the finite element methods is studied. Results of study were as follows. 1. The matrixes of EI[B] and [K] in the equations of M(x)=EI[B]{q} and M(x) = [K]{q}+{Q} of beams are same. 2. The equations of $\frac{-M}{EI}=\frac{{d^2}{\omega}}{dx^2}\;and\;\frac{dM}{dx}=V$ for the member forces have a error ratio in a finite element method of uniformly loaded structures, so equilibrium node loads {Q} must be substituted in the equation of member forces as the numerical examples of this paper revealed.

• Structural Engineering and Mechanics
• /
• v.8 no.4
• /
• pp.325-341
• /
• 1999

A procedure for the computation of the load carrying capacity of perfectly plastic plates in bending is presented. The approach, based on the kinematic theorem of limit analysis, requires the evaluation of the minimum of a convex, but non-smooth, function under linear equality constraints. A systematic solution procedure is devised, which detects and eliminates the finite elements which are predicted as rigid in the collapse mechanism, thus reducing the problem to the search for the minimum of a smooth and essentially unconstrained function of nodal velocities. Both Kirchhoff and Mindlin plate models are considered. The effectiveness of the approach is illustrated by means of some examples.

• Structural Engineering and Mechanics
• /
• v.26 no.2
• /
• pp.151-162
• /
• 2007

Rigorous analytical solutions are obtained for the plane stress problem of a rectangular plate subjected to non-linearly distributed bending loads on two opposite edges. They are then used in a Galerkin type solution to obtain the corresponding convergent buckling loads. It is shown that the critical bending moment depends significantly on the actual edge load distribution and further the number of nodal lines of the buckled configuration can also be different from that corresponding to a linear antisymmetric distribution of the bending stresses. Results are tabulated for future use while judging approximate numerical solutions.