• Korea-Australia Rheology Journal
• /
• v.18 no.2
• /
• pp.109-118
• /
• 2006

A finite-element method was employed to analyze axisymmetric unsteady motion of a deformable bubble near the wall. In the present study a deformable bubble in a Newtonian medium near the wall was considered. In solving the governing equations a structured mesh generator was used to describe the collapse of highly deformed bubbles with the Arbitrary Lagrangian Eulerian (ALE) method being employed in order to capture the transient bubble boundary effectively. In order to check the accuracy of the present FE analysis we compared the results of our FE solutions with the result of the collapse of spherical bubbles in a large body of fluid in which solutions can be obtained using a 1D FE analysis. It has been found that 1D and 2D bubble deformations are in good agreement for spherically symmetric problems confirming the validity of the numerical code. Non-spherically symmetric problems were also solved for the collapse of bubble located near a plane solid wall. We have shown that a microjet develops at the bubble boundary away from the wall as already observed experimentally. We have discussed the effect of Reynolds number and distance of the bubble center from the wall on the transient collapse pattern of bubble.

• Nuclear Engineering and Technology
• /
• v.53 no.3
• /
• pp.1029-1040
• /
• 2021

For the structural integrity evaluation of pressurized water reactor (PWR) steam generator (SG) tubes subjected to transient hydraulic loading, determination of the tube-to-tube gap velocity and static pressure distributions along the tubes is prerequisite. This paper addresses both computational fluid dynamics (CFD) and analytical approaches for predicting the tube-to-tube gap velocity and static pressure distributions during blowdown following a feedwater line break (FWLB) accident at a PWR SG. First of all, a comparative study on CFD calculations of the transient velocity and pressure distributions in the SG secondary sides for two different models having 30 or no tubes is performed. The result shows that the velocities of sub-cooled water flowing between any adjacent two tubes of a tubed SG model during blowdown can be roughly estimated by applying the specified SG secondary side porosity to those of the no-tubed SG model. Secondly, simplified analytical approximate solutions for the steady two-dimensional SG secondary flow velocity and pressure distributions under a given discharge flowrate are derived using a line sink model. The simplified analytical solutions are validated by comparing them to the CFD calculations.

• 전기의세계
• /
• v.23 no.3
• /
• pp.43-52
• /
• 1974

The matrix inversion is very inefficient for computing direct solutions of the large spare systems of linear equations that arise in many network problems as a large electrical power system. Optimally ordered triangular factorization of sparse matrices is more efficient and offers the other important computational advantages in some applications with this method. The direct solutions are computed from sparse matrix factors instead of a full inverse matrix, thereby gaining a significant advantage is speed and computer memory requirements. In this paper, it is shown that the sparse matrix method is superior to the inverse matrix method to solve the linear equations of large sparse networks. In addition, it is shown that the sparse matrix method is superior to the inverse matrix method to solve the linear equations of large sparse networks. In addition, it is shown that the solutions may be applied directly to sove the load flow in an electrical power system. The result of this study should lead to many aplications including short circuit, transient stability, network reduction, reactive optimization and others.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.6 no.3
• /
• pp.247-255
• /
• 1982

The heat diffusion equation for an annular fin is analyzed by Laplace transformation. The fin has a uniform thickness, with its end insulated, and three different temperature profiles at the base such as step change, harmonic and exponential functions. The exact solutions for the temperature and heat flux of the fins are obtained with the infinite series. The series solutions converge rapidly for large values of dimensionless time, but slowly for small values. Therefore some approximate solutions are presented here to fine the temperature distribution and heat flux for small values of dimensionless time. Furthermore a simple approximate heat flux, .OMEGA.=1.13c.tau.$^{1}$2/ is found in the range of .tau. .leg. o.1/c for the exponential function at the base.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.24 no.4
• /
• pp.9-18
• /
• 1987

The hydrodynamic radiation forces acting on a ship travelling in waves have been conventionally treated by strip theories or by direct three dimensional approaches, most of which have been formulated in frequency domain. If the forward speed of a ship varies with time, or if its path is not a straight line, conventional frequency domain analysis can no more be used, and for these cases time domain analysis may be used. In this paper, formulations are made in time domain with applications to some problems the results of which are known in frequency domain. And the results of both domains are compared to show the characteristics and validity of time domain solutions. The radiation forces acting on a three dimensional body within the framework of a linear theory. If the linearity of entire system is assumed, radiation forces due to arbitrary ship motions can be expressed by the convolution integral of the arbitrary motion velocity and the so called impulse response function. Numerical calculations are done for some bodies of simple shapes and Series-60[$C_B=0.7$] ship model. For all cases, integral equation techniques with transient Green's function are used, and velocity or acceleration potentials are obtained as the solution of the integral equations. In liner systems, time domain solutions are related with frequency domain solutions by Fourier transform. Therefore time domain solutions are Fourier transformed by suitable relations and the results are compared with various frequency domain solutions, which show good agreements.

• Journal of computational fluids engineering
• /
• v.17 no.4
• /
• pp.41-48
• /
• 2012

A numerical analysis of thermal stratification in the upper plenum of the MONJU fast breeder reactor was performed. Calculations were performed for a 1/6 simplified model of the MONJU reactor using the commercial code, CFX-13. To better resolve the geometrically complex upper core structure of the MONJU reactor, the porous media approach was adopted for the simulation. First, a steady state solution was obtained and the transient solutions were then obtained for the turbine trip test conducted in December 1995. The time dependent inlet conditions for the mass flow rate and temperature were provided by JAEA. Good agreement with the experimental data was observed for steady state solution. The numerical solution of the transient analysis shows the formation of thermal stratification within the upper plenum of the reactor vessel during the turbine trip test. The temporal variations of temperature were predicted accurately by the present method in the initial rapid coastdown period (~300 seconds). However, transient numerical solutions show a faster thermal mixing than that observed in the experiment after the initial coastdown period. A nearly homogenization of the temperature field in the upper plenum is predicted after about 900 seconds, which is a much shorter-term thermal stratification than the experimental data indicates. This discrepancy is due to the shortcoming of the turbulence models available in the CFX-13 code for a natural convection flow with thermal stratification.

• Nuclear Engineering and Technology
• /
• v.55 no.6
• /
• pp.2230-2245
• /
• 2023

KANT (KAIST Advanced Nuclear Tachygraphy) is a PWR core simulator recently developed at Korea Advance Institute of Science and Technology, which solves three-dimensional steady-state and transient multigroup neutron diffusion equations under Cartesian geometries alongside the incorporation of thermal-hydraulics feedback effect for multi-physics calculation. It utilizes the standard Nodal Expansion Method (NEM) accelerated with various Coarse Mesh Finite Difference (CMFD) methods for neutronics calculation. For thermal-hydraulics (TH) calculation, a single-phase flow model and a one-dimensional cylindrical fuel rod heat conduction model are employed. The time-dependent neutronics and TH calculations are numerically solved through an implicit Euler scheme, where a detailed coupling strategy is presented in this paper alongside a description of nodal equivalence, macroscopic depletion, and pin power reconstruction. For validation of the steady, transient, and depletion calculation with pin power reconstruction capacity of KANT, solutions for various benchmark problems are presented. The IAEA 3-D PWR and 4-group KOEBERG problems were considered for the steady-state reactor benchmark problem. For transient calculations, LMW (Lagenbuch, Maurer and Werner) LWR and NEACRP 3-D PWR benchmarks were solved, where the latter problem includes thermal-hydraulics feedback. For macroscopic depletion with pin power reconstruction, a small PWR problem modified with KAIST benchmark model was solved. For validation of the multi-physics analysis capability of KANT concerning large-sized PWRs, the BEAVRS Cycle1 benchmark has been considered. It was found that KANT solutions are accurate and consistent compared to other published works.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2005.06a
• /
• pp.965-968
• /
• 2005

Roll drafting operation causes variations in the linear density of bundles because the bundle flow cannot be controlled completely by roll pairs. Defects occurring in this operation bring about many problems successively in the next processes. In this paper, we attempt to analyze the draft dynamics and the linear density irregularity based on the governing equation of a bundle motion that has been suggested in our previous studies. For analyzing the dynamic characteristics of the roll drafting operation, it is indispensable to investigate a transient state in time domain before the bundle flux reaches a steady state. However, since governing equations of bundle flow consisting of continuity and motion equations turn out to be nonlinear, and coupled between variables, the solutions for a transient state cannot be obtained by an analytical method. Therefore, we use the Finite Difference Method(FDM), particularly, the FTBS(Forward-Time Backward-Space) difference method. Then, the total equations system yields to an algebraic equations system and is solved under given initial and boundary conditions in an iterative fashion. From the simulation results, we confirm that state variables show different behavior in the transient state; e.g., the velocity distribution in the flow field changes more quickly the linear density distribution. During a transient flow in a drafting zone, the output irregularity is influenced differently by the disturbances, e.g., the variation in input bundle thickness, the drafting speed, and the draft ratio.

• Journal of the Microelectronics and Packaging Society
• /
• v.28 no.2
• /
• pp.13-28
• /
• 2021

Since transient electronic devices can operate under harsh conditions such as electrolytic solutions or inside the body, and be removed by hydrolysis after operation, they can replace conventional electronic devices in various research areas like biomedical implantable devices. Moreover, transient electronic devices that can dissolve in water and enzymes are the focus of the new concept of green technology, which can solve electrical waste issues. However, the surroundings of transient electronic devices can deteriorate internal device components. Thus, an encapsulation strategy is introduced for stable operation in solution by shielding the outside of a device with a passive barrier. This article summarizes recent research trends in transient electronic devices, including their background, dissolution behavior, and encapsulation strategies to enhance reliability by blocking water permeation.

• Nuclear Engineering and Technology
• /
• v.43 no.6
• /
• pp.531-546
• /
• 2011

The backward differentiation formula (BDF) method is applied to a three-dimensional reactor kinetics calculation for efficient yet accurate transient analysis with adaptive time step control. The coarse mesh finite difference (CMFD) formulation is used for an efficient implementation of the BDF method that does not require excessive memory to store old information from previous time steps. An iterative scheme to update the nodal coupling coefficients through higher order local nodal solutions is established in order to make it possible to store only node average fluxes of the previous five time points. An adaptive time step control method is derived using two order solutions, the fifth and the fourth order BDF solutions, which provide an estimate of the solution error at the current time point. The performance of the BDF- and CMFD-based spatial kinetics calculation and the adaptive time step control scheme is examined with the NEACRP control rod ejection and rod withdrawal benchmark problems. The accuracy is first assessed by comparing the BDF-based results with those of the Crank-Nicholson method with an exponential transform. The effectiveness of the adaptive time step control is then assessed in terms of the possible computing time reduction in producing sufficiently accurate solutions that meet the desired solution fidelity.