• 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 the Convergence on Culture Technology
• /
• v.9 no.4
• /
• pp.423-430
• /
• 2023

This dissertation discusses the formulation of a nonlinear storytelling board that preserves the contextual perspective of characters. Storytelling encompasses the director's creative intention by leveraging the interaction of various elements to construct a logical narrative that explores cause and effect. Its primary objective is to enhance viewers' empathy. Consequently, there is a pressing need for comprehensive research on differentiating storytelling from storyboarding. Moreover, the integrated approach to storytelling and storyboarding holds scholarly value in understanding the process of narrative composition and visualization. Thus, a study proposes a method for constructing nonlinear storytelling boards considering the discrete camera perspective and contextual scene continuity, ultimately contributing to visual complexity and correlation comprehension. This approach enables a careful and simultaneous consideration of the correlations that deepen cognition, including the physical, emotional, and event rhythms mentioned in Karen Perlman's theory.

• Journal of Korea Society of Industrial Information Systems
• /
• v.28 no.5
• /
• pp.67-75
• /
• 2023

This paper deals with the operation of air taxis, which is one of the latest innovative technologies aimed at solving the issue of traffic congestion in cities. A key challenge for the successful introduction of the technology and efficient operation is a vertiport location problem. This paper employs a discrete choice model to calculate choice probabilities of transportation modes for each route, taking into account factors such as cost and travel time associated with different modes. Based on this probability, a mathematical formulation to maximize the utilization rate for air taxi is proposed. However, the proposed model is NP-hard, effective and efficient solution methodology is required. Compared to previous studies that simply proposed the optimization models, this study presents a solution methodology using the cross-entropy algorithm and confirms the effectiveness and efficiency of the algorith through numerical experiments. In addition to the academic excellence of the algorithm, it suggests that decision-making that considers actual data and air taxi utilization plans can increase the practial usability.

• Structural Engineering and Mechanics
• /
• v.89 no.2
• /
• pp.181-197
• /
• 2024

A numerical method is presented in this paper, for buckling analysis of thin arbitrary stiffened composite cylindrical shells under axial compression. The stiffeners can be placed inside and outside of the shell. The shell and stiffeners are operated as discrete elements, and their interactions are taking place through the compatibility conditions along their intersecting lines. The governing equations of motion are obtained based on Koiter's theory and solved by utilizing the principle of the minimum potential energy. Then, the buckling load coefficient and the critical buckling load are computed by solving characteristic equations. In this formulation, the elastic and geometric stiffness matrices of a single curved strip of the shell and stiffeners can be located anywhere within the shell element and in any direction are provided. Moreover, five stiffened composite shell specimens are made and tested under axial compression loading. The reliability of the presented method is validated by comparing its numerical results with those of commercial software, experiments, and other published numerical results. In addition, by using the ANSYS code, a 3-D finite element model that takes the exact geometric arrangement and the properties of the stiffeners and the shell into consideration is built. Finally, the effects of Poisson's ratio, shell length-to-radius ratio, shell thickness, cross-sectional area, angle, eccentricity, torsional stiffness, numbers and geometric configuration of stiffeners on the buckling of stiffened composite shells with various end conditions are computed. The results gained can be used as a meaningful benchmark for researchers to validate their analytical and numerical methods.

• The Journal of The Korea Institute of Intelligent Transport Systems
• /
• v.15 no.6
• /
• pp.116-126
• /
• 2016

In this study, lane-wise traffic flow characteristics were analysed on uninterrupted flow using a new notion of "lane-wise travel speed reversal (LTSR)" which is defined as a phenomena that travel speed in the median lane is lower than other lanes. Mathematical formulation was also proposed to calculate the strength of LTSR. The experiment road site is Seoul Outer Ring Expressway (Jayuro-IC~Jangsoo-IC), and travel trajectories for each four lane were collected for weekdays (Mon. through Fri.) during morning peak. Comparing lane-wise travel speeds for entire test road section, no LTSR was observed, meaning that the travel speed in the median lane is the fastest, followed by 2nd, 3rd, and 4th lane as in order. Howerver, the result of microscopic analysis using 100-meter discrete road section based data shows that LTSR occurs many times. Especially the strength of LTSR is higher in congestion area and freeway merge and diverge segment. It is expected that these results could be used as a fundamental data when establishing lane-by-lane traffic operation strategy and developing lane-wise traffic information collection and dissemination technology.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.11 no.3
• /
• pp.179-196
• /
• 1986

In this paper, the peformance of a statistical packet voice/data multiplexer is studied. In ths study we assume that in the packet voice/data multiplexer two separate finite queues are used for voice and data traffics, and that voice traffic gets priority over data. For the performance analysis we divide the output link of the multiplexer into a sequence of time slots. The voice signal is modeled as an (M+1) - state Markov process, M being the packet generation period in slots. As for the data traffic, it is modeled by a simple Poisson process. In our discrete time domain analysis, the queueing behavior of voice traffic is little affected by the data traffic since voice signal has priority over data. Therefore, we first analyze the queueing behavior of voice traffic, and then using the result, we study the queueing behavior of data traffic. For the packet voice multiplexer, both inpur state and voice buffer occupancy are formulated by a two-dimensional Markov chain. For the integrated voice/data multiplexer we use a three-dimensional Markov chain that represents the input voice state and the buffer occupancies of voice and data. With these models, the numerical results for the performance have been obtained by the Gauss-Seidel iteration method. The analytical results have been verified by computer simylation. From the results we have found that there exist tradeoffs among the number of voice users, output link capacity, voic queue size and overflow probability for the voice traffic, and also exist tradeoffs among traffic load, data queue size and oveflow probability for the data traffic. Also, there exists a tradeoff between the performance of voice and data traffics for given inpur traffics and link capacity. In addition, it has been found that the average queueing delay of data traffic is longer than the maximum buffer size, when the gain of time assignment speech interpolation(TASI) is more than two and the number of voice users is small.

• Journal of Navigation and Port Research
• /
• v.34 no.3
• /
• pp.195-203
• /
• 2010

This paper concerns about total dynamic analysis of integrated mining system. This system consists of vertical steel pipe, intermediate buffer station, flexible pipe and self-propelled miner. The self-propelled miner and buffer are assumed as rigid-body of 6-dof. Discrete models of vertical steel pipe and flexible pipe are adopted, which are obtained by means of lumped-parameter method. The motion of mining vessel is not considered. Instead, the motion of mining vessel is taken into account in form of various boundary conditions (e.g. forced excitation in slow motion and/or fast oscillation and so on). A terramechanics model of extremely cohesive soft soil is applied to the self-propelled miner. Hinged and ball constraints are used to define the connections between sub-systems (vertical steel pipe, buffer, flexible pipe, self-propelled miner). Equations of motion of the coupled model are derived with respect to the each local coordinates system. Four Euler parameters are used to express the orientations of the sub-systems. To solve the equations of motion of the total dynamic model, an incremental-iterative formulation is employed. Newmark-${\beta}$ method is used for time-domain integration. The total dynamic responses of integrated mining system are investigated.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.3
• /
• pp.167-174
• /
• 2021

Recent investigation into the integrity of nuclear containment buildings has highlighted the importance of developing an elaborate diagnostic method to evaluate the distribution and size of cavities inside concrete walls. As part of developing such a method, this paper presents a finite element approach to modeling elastic waves propagating in the containment building walls of a nuclear power plant. We introduce a perfectly matched layer (PML) wave-absorbing boundary to limit the large-scale nuclear containment wall to the region of interest. The formulation results in a semi-discrete form with symmetric damping and stiffness matrices. The transient elastic wave equations for a mixed unsplit-field PML were solved for displacement and stresses in the time domain. Numerical results show that the sensitivity of displacement, velocity, acceleration, and stresses is large depending on the size and location of the cavity. The dynamic response of the wall slightly differs depending on the existence of the containment liner plate. The results of this study can be applied to a full-waveform inversion approach for characterizing cavities inside a containment wall.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.45-59
• /
• 2006

The successful performance of any numerical geotechnical simulation depends on the accuracy and efficiency of the numerical implementation of constitutive model used to simulate the stress-strain (constitutive) response of the soil. The corner stone of the numerical implementation of constitutive models is the numerical integration of the incremental form of soil-plasticity constitutive equations over a discrete sequence of time steps. In this paper a well known two-surface soil plasticity model is implemented using a generalized implicit return mapping algorithm to arbitrary convex yield surfaces referred to as the Closest-Point-Projection method (CPPM). The two-surface model describes the nonlinear behavior of coarse-grained materials by incorporating a bounding surface concept together with isotropic and kinematic hardening as well as fabric formulation to account for the effect of fabric formation on the unloading response. In the course of investigating the performance of the CPPM integration method, it is proven that the algorithm is an accurate, robust, and efficient integration technique useful in finite element contexts. It is also shown that the algorithm produces a consistent tangent operator $\frac{d\sigma}{d\varepsilon}$ during the iterative process with quadratic convergence rate of the global iteration process.