• Nuclear Engineering and Technology
• /
• v.20 no.1
• /
• pp.1-8
• /
• 1988

During power plant operation, the flow rate from the CVCS makeup system is estimated using the continuity equation and mass balance equation, when the primary loop boron concentration change is required due to the power transient. For this purpose, primary loop, pressurizer and VCT(volume control tank)(in CVCS) are modeled by three control volumes which contain each mass and boron concentration. Connecting pipes between primary loop, pressurizer and CVCS are also modeled by time delay. Calculation for 14-2-6-2 (power 100-50-100) load follow case (at EOL, for KNU 7) is made using these models.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.20-24
• /
• 2011

In this paper, numerical simulation of cavitation flow for modified NACA66 hydrofoil was made by using the multi-phase RANS equation based on pseudo-compressibility. The Homogeneous mixture model comprised of the mixture continuity, mixture momentum and liquid volume fraction equations was utilized. A vertex-centered finite-volume method was used in conjunction 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing The Spalart-Allmaras one equation model was employed for the closure of turbulence. Reasonable agreements were obtained between the calculation results and the experiment for pressure coefficients on the hydrofoil surface.

• Journal of Mechanical Science and Technology
• /
• v.14 no.7
• /
• pp.789-795
• /
• 2000

Free-surface flows with an arbitrary deformation, induced by a submerged hydrofoil, are simulated numerically, considering two-fluid flows of both water and air. The computation is performed by a finite volume method using unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell wise local mesh refinement. The integration in space is of second order, based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels. The linear equations are solved by conjugate gradient type solvers, and the non-linearity of equations is accounted for through Picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations, the continuity equation, the conservation equation of one species, and the equations for two turbulence quantities. Finally, a comparison is quantitatively made at the same speed between the computation and experiment in which the grid sensitivity is numerically checked.

• Journal of Advanced Marine Engineering and Technology
• /
• v.31 no.8
• /
• pp.954-960
• /
• 2007

The purpose of this paper is to obtain design method of air-distribution system. Air-distribution system is composed of blower, duct, diffusers and measuring equipment. The air-flow rate from each diffuser is not equal. The air-flow rate is calculated with the combined equations which are Bernoulli's equation, continuity equation and minor loss equations. Inlet condition and outlet condition are adapted in each duct system. Then square difference between function of maximum air-flow rate and minimum air-flow rate is used as an object function. Area of diffuser and velocity are established as constraints. To minimize the object function, the optimization method is used. After optimization the design variables are selected under satisfaction of constraints. The air-distribution system is calculated again with the result of optimized design variable. It is shown that the air-distribution system has the equal air-flow rate from diffusers.

• Proceedings of the KIEE Conference
• /
• 1999.11b
• /
• pp.101-103
• /
• 1999

This paper presents an analysis of 2-dimensional(2-D) multi-regions problem using boundary integral equation method(BIEM). When compared with finite element method(FEM), there are only a few unknown variables in BIEM because it implements numerical analysis only for the surface or boundary of a model. As a result, a lot of computational memory and time can be saved. Procedure to analyze 2-D multi-regions problem using potentials and its derivatives in a boundary as unknown variables, first, numerical analysis is performed for each of subregions. And then interface continuity condition is applied to the interface between them and Gauss Quadrature Formula are adopted to solve singular integral in a boundary in this paper.

• Water for future
• /
• v.24 no.1
• /
• pp.109-117
• /
• 1991

The governing equation of flow in porous media is developed on the bases of the continuity equation of fluid for transient flow through a saturated-unsaturated zone, and substitution of Darcy's law. The numerical solutions are obtained by finite element method based on the Galerkin principles weighted residuals. The analysis are carried out by using the unsteady storm data observed and rainfall intensities which are obtained by using the rainfall excess model in considering of the initial losses. The functional relationships between the hydraulic conductivity, capillary pressure head and volumetric water content are applied to the flow of water through unsaturated soil varied with changes of water content.

• Korean Journal of Hydrosciences
• /
• v.2
• /
• pp.13-24
• /
• 1991

The governing partial differential equation of flow in porous media is developed on the bases of the continuity equation of fluid for transient flow through a saturated-unsaturated zone, and substitution of Dercy's law. The numerical solution is obtained by the Galerkin finite element method based on the principle of weighted residuals. The analysis is carried out by using the unsteady storm data observed and the functional relationships between the hydraulic conductivities, capillary pressure heads, and volumetric water contents under saturated-unsaturated conditions. As the results the hydraulic conductivities, rates of change of storage and initial moisture conditions are significantly influened on the responses of subsurface flow on a hillslope.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 1998.04a
• /
• pp.98-105
• /
• 1998

In this study, the boundary element analysis in dynamics for the multilayered semi-infinite plane is developed using the fundamental solution for moving loads. Also the indirect method and superposition method are introduced to consider the multilayered systems and moving loads. At each layer the fundamental solution can be obtained by solving the governing equation which is transformed by the Fourier transform. The governing equation can be solved by three conditions; continuity conditions of displacement and stress, the traction free condition at the surface and the radiation condition at the surface and the radiation condition at the infinite distance. To verify the solution and the developed algorithm, the theoretical solution for the homogeneous layer and commercial FEM program is compared with the results of this study.

• Proceedings of the KIEE Conference
• /
• 1994.07b
• /
• pp.1533-1535
• /
• 1994

The Monte-Carlo simulation was used to define the physical mechanisms of the initiation phase of pseudo-spark discharge. The pseudo-spark discharge employing the hollow cathode geometry is accompanied by very fast current rising and intense charged particle beams. In this model, time-dependent continuity equation for the electrons and ions were solved consistently with Poisson's equation for the electric field in a two-dimensional, sysmmetrically cylinderical geometry. From the simulation, a sequence of physical mechanisms that cause the rapid current rise associated with the onset of pseudo-spark discharge mode were identified.

• Journal of electromagnetic engineering and science
• /
• v.14 no.3
• /
• pp.273-277
• /
• 2014

In this research, integral transform technique for electromagnetic scattering formulation is reviewed. Electromagnetic boundary-value problems are presented to demonstrate how the integral transforms are utilized in electromagnetic propagation, antennas, and electromagnetic interference/compatibility. Various canonical structures of slotted conductors are used for illustration; moreover, Fourier transform, Hankel transform, Mellin transform, Kontorovich-Lebedev transform, and Weber transform are presented. Starting from each integral transform definition, the general procedures for solving Helmholtz's equation or Laplace's equation for the potentials in the unbounded region are reviewed. The boundary conditions of field continuity are incorporated into particular formulations. Salient features of each integral transform technique are discussed.