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

• Transactions of the Korean Society of Automotive Engineers
• /
• v.12 no.2
• /
• pp.54-61
• /
• 2004

The engine-out flow is highly transient and hot, and may place tremendous thermal and inertial loads on a closed-coupled catalyst. Therefore, time-dependent and detailed flow and thermal field simulation may be crucial. The aim of this study is to develop combined chemical reaction and multi-dimensional fluid dynamic mathematical model and to study the effect of unsteady pulsating thermal and flow characteristics on thermal reliability of closed-coupled catalyst. The effect of cell density on the conversion performance under real running condition is also investigated. Unlike previous studies, the present study focuses on coupling between the problems of pulsating flow pattern and catalyst thermal response and conversion efficiency. The results are expressed in terms of temporal evolution of flow, pollutant and temperature distribution as well as transient characteristics of conversion efficiency. Fundamental understanding of the flow and thermal phenomena of closed-coupled catalyst under real running condition is presented. It is shown that instants of significantly low values of flow uniformity and conversion efficiency exist during exhaust blowdown and the temporal varaition of flow uniformity is very similar in pattern to one of conversion efficiency. It is also found that the location of hot spot in monolith is directly affected by transient flow pattern in closed-coupled catalyst.

• Journal of the Microelectronics and Packaging Society
• /
• v.16 no.4
• /
• pp.5-8
• /
• 2009

In this paper, a high-quality low-temperature co-fired ceramic (LTCC) solenoid inductor using a solder ball and an air cavity on a silicon wafer for three-dimensional (3-D) system-in-package (SiP) is proposed. The LTCC multi-layer solenoid inductor is attached using Ag paste and solder ball on a silicon wafer with the air cavity structure. The air cavity is formed on a silicon wafer through an anisotropic wet-etching technology and is able to isolate the LTCC dielectric loss which is equivalent to a low k material effect. The electrical coupling between the metal layer and the LTCC dielectric layer is decreased by adopting the air cavity. The LTCC solenoid inductor using the solder ball and the air cavity on silicon wafer has an improved Q factor and self-resonant frequency (SRF) by reducing the LTCC dielectric resistance and parasitic capacitance. Also, 3-D device stacking technologies provide an effective path to the miniaturization of electronic systems.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.20 no.4
• /
• pp.397-404
• /
• 2010

The use of equivalent circuit models of piezoelectric energy harvesters is inevitable when power circuitry including rectifying and smoothing circuit elements is connected to them for evaluating DC electrical outputs. This is because it is difficult to incorporate the electro-mechanical coupling resulting from the additional circuitry into the conventional finite element analysis. Motivated by this observation, we propose a method to accurately extract the equivalent circuit parameters by using commercially available FEM software such as ANSYS which provides three-dimensional AC piezoelectric analysis. Then the equivalent circuit can be analyzed by circuit simulators such as $SimPowerSystems^{TM}$ of MATLAB. While the previous works have estimated the circuit parameters by experimental measurements or by analytical solutions developed only for limited geometries and boundary conditions, the proposed method has no such limitation because piezoelectric energy harvesters of any shapes and boundary conditions can be treated in FEM software. For the verification of the proposed method, multi-modal AC electrical power output by using a corresponding equivalent circuit is compared with that by ANSYS. The proposed method is then shown to be very useful in the subsequent evaluation of DC electrical output which is obtained by attaching a bridge diode and a storage capacitor to a piezoelectric energy harvester.

• Nuclear Engineering and Technology
• /
• v.53 no.8
• /
• pp.2477-2487
• /
• 2021

To improve safety analysis technology for a nuclear reactor containment considering an interaction between a reactor coolant system (RCS) and containment, this study aims at an experimental investigation on the integrated simulation of the RCS and containment, with an integral effect test facility, ATLAS-CUBE. For a realistic simulation of a pressure and temperature (P/T) transient, the containment simulation vessel was designed to preserve a volumetric scale equivalently to the RCS volume scale of ATLAS. Three test cases for a steam line break (SLB) transient were conducted with variation of the initial condition of the passive heat sink or the steam flow direction. The test results indicated a stratified behavior of the steam-gas mixture in the containment following a high-temperature steam injection in prior to the spray injection. The test case with a reduced heat transfer on the passive heat sink showed a faster increase of the P/T inside the containment. The effect of the steam flow direction was also investigated with respect to a multi-dimensional distribution of the local heat transfer on the passive heat sink. The integral effect test data obtained in this study will contribute to validating the evaluation methodology for mass and energy (M/E) and P/T transient of the containment.

• Journal of Drive and Control
• /
• v.17 no.4
• /
• pp.113-124
• /
• 2020

Most of the work of construction equipment and agricultural machinery is done in off-road conditions. Autonomous driving required in these conditions uses GPS sensors, and PID controllers to control their speed and position. The hydrostatic steering, which is composed of a PSU, hydraulic hoses, and cylinders, rather than a mechanical coupling is used in these equipments. The PSU plays a key role in hydrostatic steering. Precise control of the position under various conditions requires detailed behavioral analysis of the basic components and operation. Two Gerotor PSU is now a commonly used safer option. The components of the PSU can be divided into mechanical and hydraulic actuating elements by its behavior. Since the system is combined by mechanical and hydraulic elements, the modelings are performed using Amesim, which is one of the most effective for the multi-domain dynamic system analysis. To confirm the validity of the model, input torque and pressures are checked with varying steering speed. The opening and the steering speed of normal and newly designed control valve set is investigated with the effect of centering spring force and friction. Finally, simulation results with fully detailed model with two gerotors are analyzed and compared with simple model.

• 한국신재생에너지학회:학술대회논문집
• /
• 2008.05a
• /
• pp.33-36
• /
• 2008

The transient behavior of a passive air breathing direct methanol fuel cell (DMFC) operated on vapor-feeding mode is studied in this paper. It generally takes 30 minutes after starting for the cell response to come to its steady-state and the response is sometimes unstable. A mathematical dynamic one-dimensional model for simulating transient response of the DMFC is presented. In this model a DMFC is decomposed into its subsystems using lumped model and divided into five layers, namely the anodic diffusion layer, the anodic catalyst layer, the proton exchange membrane (PEM), the cathodic catalyst layer and the cathodic diffusion layer. All layers are considered to have finite thickness, and within every one of them a set of differential-algebraic governing equations are given to represent multi-components mass balance, such as methanol, water, oxygen and carbon dioxide, charge balance, the electrochemical reaction and mass transport phenomena. A one-dimensional, isothermal and mass transport model is developed that captures the coupling between water generation and transport, oxygen consumption and natural convection. The single cell is supplied by pure methanol vapor from a methanol reservoir at the anode, and the oxygen is supplied via natural air-breathing at the cathode. The water is not supplied from external source because the cell uses the water created at the cathode using water back diffusion through nafion membrane. As a result of simulation strong effects of water transport were found out. The model analysis provides several conclusions. The performance drop after peak point is caused by insufficiency of water at the anode. The excess water at the cathode makes performance recovery impossible. The undesired crossover of the reactant methanol through the PEM causes overpotential at the cathode and limits the feeding methanol concentration.

• Nuclear Engineering and Technology
• /
• v.36 no.6
• /
• pp.528-539
• /
• 2004

It is essential in commercial reactors that the safety limits imposed on the fuel pellets and fuel clad barriers, such as the linear power density (LPD) and the departure from nucleate boiling ratio (DNBR), are not violated during reactor operations. In order to accurately monitor the safety limits of current reactor states, a detailed three-dimensional (3D) core power distribution should be estimated from the in-core detector signals. In this paper, we propose a calculation methodology for detailed 3D core power distribution, using in-core detector signals and core monitoring constants such as the 3D Coupling Coefficients (3DCC), node power fraction, and pin-to-node factors. Also, the calculation method for several core safety parameters is introduced. The core monitoring constants for the real core state are promptly provided by the core design code and on-line MASTER (Multi-purpose Analyzer for Static and Transient Effects of Reactors), coupled with the core monitoring program. through the plant computer, core state variables, which include reactor thermal power, control rod bank position, boron concentration, inlet moderator temperature, and flow rate, are supplied as input data for MASTER. MASTER performs the core calculation based on the neutron balance equation and generates several core monitoring constants corresponding to the real core state in addition to the expected core power distribution. The accuracy of the developed method is verified through a comparison with the current CECOR method. Because in all the verification calculation cases the proposed method shows a more conservative value than the best estimated value and a less conservative one than the current CECOR and COLSS methods, it is also confirmed that this method secures a greater operating margin through the simulation of the YGN-3 Cycle-1 core from the viewpoint of the power peaking factor for the LPD and the pseudo hot pin axial power distribution for the DNBR calculation.

• Geomechanics and Engineering
• /
• v.11 no.1
• /
• pp.77-94
• /
• 2016

A forked tunnel, as a special complicated underground structure, is composed of big-arch tunnel, multi-arch tunnel, neighborhood tunnels and separate tunnels according to the different distances between two separate tunnels. Due to the complicated process of design and construction, surrounding jointed rock mass stability of the big-arch tunnel which belongs to the forked tunnel during excavation is a hot issue that needs special attentions. In this paper, an elasto-plastic damage constitutive model for jointed rock mass is proposed based on the coupling method considering elasto-plastic and damage theories, and the irreversible thermodynamics theory. Based on this elasto-plastic damage constitutive model, a three dimensional elasto-plastic damage finite element code (D-FEM) is implemented using Visual Fortran language, which can numerically simulate the whole excavation process of underground project and perform the structural stability of the surrounding rock mass. Comparing with a popular commercial computer code, three dimensional fast Lagrangian analysis of continua (FLAC3D), this D-FEM has advantages in terms of rapid computing process, element grouping function and providing more material models. After that, FLAC3D and D-FEM are simultaneously used to perform the structural stability analysis of the surrounding rock mass in the forked tunnel considering three different computing schemes. The final numerical results behave almost consistent using both FLAC3D and D-FEM. But from the point of numerically obtained damage softening areas, the numerical results obtained by D-FEM more closely approach the practical behaviors of in-situ surrounding rock mass.

• Smart Structures and Systems
• /
• v.19 no.2
• /
• pp.181-194
• /
• 2017

In order to investigate the interface debonding defects detection mechanism between steel tube and concrete core of concrete-filled steel tubes (CFSTs), multi-physical fields coupling finite element models constituted of a surface mounted Piezoceramic Lead Zirconate Titanate (PZT) actuator, an embedded PZT sensor and a circular cross section of CFST column are established. The stress wave initiation and propagation induced by the PZT actuator under sinusoidal and sweep frequency excitations are simulated with a two dimensional (2D) plain strain analysis and the difference of stress wave fields close to the interface debonding defect and within the cross section of the CFST members without and with debonding defects are compared in time domain. The linearity and stability of the embedded PZT response under sinusoidal signals with different frequencies and amplitudes are validated. The relationship between the amplitudes of stress wave and the measurement distances in a healthy CFST cross section is also studied. Meanwhile, the responses of PZT sensor under both sinusoidal and sweep frequency excitations are compared and the influence of debonding defect depth and length on the output voltage is also illustrated. The results show the output voltage signal amplitude and head wave arriving time are affected significantly by debonding defects. Moreover, the measurement of PZT sensor is sensitive to the initiation of interface debonding defects. Furthermore, wavelet packet analysis on the voltage signal under sweep frequency excitations is carried out and a normalized wavelet packet energy index (NWPEI) is defined to identify the interfacial debonding. The value of NWPEI attenuates with the increase in the dimension of debonding defects. The results help understand the debonding defects detection mechanism for circular CFST members with PZT technique.