• Nuclear Engineering and Technology
• /
• v.41 no.8
• /
• pp.1045-1052
• /
• 2009

A calculation model is developed to predict the transient free surface flow on the containment floor following a loss-of-coolant accident (LOCA) of pressurized water reactors (PWR) for the use of debris transport evaluation. The model solves the two-dimensional Shallow Water Equation (SWE) using a finite volume method (FVM) with unstructured triangular meshes. The numerical scheme is based on a fully explicit predictor-corrector method to achieve a fast-running capability and numerical accuracy. The Harten-Lax-van Leer (HLL) scheme is used to reserve a shock-capturing capability in determining the convective flux term at the cell interface where the dry-to-wet changing proceeds. An experiment simulating a sudden break of a water reservoir with L-shape open channel is calculated for validation of the present model. It is shown that the present model agrees well with the experiment data, thus it can be justified for the free surface flow with accuracy. From the calculation of flow field over the simplified containment floor of APR1400, the important phenomena of free surface flow including propagations and interactions of waves generated by local water level distribution and reflection with a solid wall are found and the transient flow rates entering the Holdup Volume Tank (HVT) are obtained within a practical computational resource.

• 제어로봇시스템학회:학술대회논문집
• /
• 1994.10a
• /
• pp.665-671
• /
• 1994

Some extended results in the study of two-position alignment for strapdown inertial navigation system are presented. In [1], an observability analysis for two-position alignment was done by analytic rank test of the stripped observability matrix and numerical calculation of the error covariance propagation using ten-state error model. In this paper, it is done by an analytic approach which utilizes the nonsingular condition of the determinant of simplified stripped observability matrix and by numerical calculation of the error covariance propagation accomplished in more cases than [1], and the twelve-state error model including vertical channel is used instead of ten-state error model. In addition, it is confirmed that this approach more clearly produces the same result as shown in the original work in terms of complete observability and there exist some better two-position configurations than [1] using the twelve-state error model.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.2
• /
• pp.470-483
• /
• 1990

Turbulent jet in a crossflow, issuing from a row of holes on a convex surface of 90 .deg. bend duct, is predicted by a 3-dimensional numerical method. The Cartesian coordinate system in adopted in upstream and downstream tangents and the cylindrical polar coordinate system in curved region. The Reynolds stresses and heat fluxes are obtained from a standard k-e model in the core region and van Driest model in the vicinity of the wall. The governing equations are discretized by a finite volume method and solutions are obtained by a locally elliptic calculation procedure. Pressure and convective terms are treated by SIMPLE algorithm and hybrid scheme respectively. A vortex initially induced by the injected jet has been built up due to the interaction with the secondary flow caused by pressure gradient and centrifugal force. The vortex structure has a strong influence on the wall cooling effectiveness. Another vortex like horseshoe is formed in the vicinity of the injection hole and its strength is getting weak as it moves downward.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.2
• /
• pp.95-103
• /
• 2013

Prediction of the running posture is important to evaluate the resistance by the numerical calculation for a high speed vessel. Especially for a planing craft having a large variation of running attitude it becomes more essential, but it can not be obtained easily because the running posture and the hydrodynamic forces including the resistance are interacted with each other. So iterative calculation to obtain the dynamic forces according to the changes in attitude is necessary, in this study, considering the calculated hydrodynamic force at the assumed draft as the additional buoyancy the corrected draft is calculated through satisfying the equilibrium between the buoyancy and the hull weight. To verify the derived method three kinds of hull forms were used with the results of model tests, R/V ATHENA and 150 tons class guide vessel for middle-speed semi-planing crafts, 28 feet fast boat for a high-speed planing boat. For all cases with several iterations the converged value of draft can be obtained, lastly the resistance and flow around hull were simulated by using VOF method.

• Journal of Korean Society for Atmospheric Environment
• /
• v.15 no.6
• /
• pp.757-765
• /
• 1999

Lagrangian particle dispersion model(LPDM) is an effective tool to calculate the dispersion from a point source since it dose not induce numerical diffusion errors in solving the pollutant dispersion equation. Fictitious particles are released to the atmosphere from the emission source and they are then transported by the mean velocity and diffused by the turbulent eddy motion in the LPDM. The concentration distribution from the dispersed particles in the calculation domain are finally estimated by applying a particle count method or a Gaussian kernel method. The two methods for calculating concentration profiles were compared each other and tested against the analytic solution and the tracer experiment to find the strength and weakness of each method and to choose computationally time saving method for the LPDM. The calculated concentrations from the particle count method was heavily dependent on the number of the particles released at the emission source. It requires lots fo particle emission to reach the converged concentration field. And resulting concentrations were also dependent on the size of numerical grid. The concentration field by the Gaussian kernel method, however, converged with a low particle emission rate at the source and was in good agreement with the analytic solution and the tracer experiment. The results showed that Gaussian kernel method was more effective method to calculate the concentrations in the LPDM.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.10
• /
• pp.1384-1391
• /
• 2001

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67) and to evaluate the performances of Abid's low-Reynolds-number k-$\omega$ and Baldwin-Lomax turbulence models. A finite volume method is used fur spatial discretization. The equations are solved implicitly in time by the use of approximate factorization. The upwind difference scheme is used for inviscid terms and viscous terms are approximated with central difference. The flux-difference-splitting method of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. The results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, it is concluded that Abid'k-$\omega$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost the same.

• Journal of KSNVE
• /
• v.9 no.5
• /
• pp.955-965
• /
• 1999

Centrifugal fans are widely used and the noise generated by the these machines causes one of the most serious problems. In general, the centrifugal fan noise is often dominated by tones at BPF(blade passage frequency) and its higher harmonics. This is a consequence of the strong interaction between the flow discharged form the impeller and the cutoff in the casing. However, only a few researches have been carried out on predicting the noise because of the difficulty in obtaining detailed information about the flow field and casing effects on noise radiation. The objective of this study is to understand the generation mechanism of sound and to develop a prediction method for the unsteady flow field and the acoustic pressure field of a centrifugal fan. We assume that the impeller rotates with a constant angular velocity and the flow field of the impeller is incompressible and inviscid. So, a discrete vortex method (DVM) is used to model the centrifugal by the unsteady Bernoulli equation. Lowson's method is used to predict the acoustic source. A centrifugal impeller and wedge introduced by Weidemann are used in the numerical calculation and the results are compared with the experimental data. Reasonable results are obtained not only for the peak frequencies but also for the amplitudes of the tonal sound.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.3
• /
• pp.71-79
• /
• 1992

In order to obtain the wave-making resistance of a ship, so-called the Neumann-Kelvin problem is solved numerically. For computing the Havelock source, which is the Green's function of the problem, we adopted the methods given by Newman(1987) for the term representing the local disturbance, and Baar and Price(1988) for the wave disturbance, respectively. In the numerical code we developed, the source strength is assumed as bilinear on each panel and continuous throughout the hull surface. The wave-making resistance is calculated using the algorithm of de Sendagorta and erases(1988), which makes use of the wave amplitude far downstream. The Wigley hull was chosen for the sample calculation, and our results showed a good agreement with other existing experimental and numerical results.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.8
• /
• pp.47-54
• /
• 2006

In this paper, numerical algorithms applied in two-step eulerian scheme are investigated and implemented. Element quantities are advected with donor or van Leer algorithm. Nodal quantities are advected with Simplifed ALE [SALE] algorithm. Material interfaces are determined with Simple Line Interface Calculation [SLIC] algorithm. Practical aspects considered for code development are addressed in detail. The results of developed two-step Eulerian code are verified by comparing with those from pure lagrangian scheme and commercial code.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.3 no.6
• /
• pp.8-22
• /
• 1995

Three dimensional numerical study of non-evaporating and evaporating spray characteristics was performed in a quiescent and motoring condition of direct injection diesel engine. The calculation parameter was breakup model. The breakup models used were Reitz & Diwakar model and TAB model. The modified k-${\varepsilon}$ turbulence model considering the compressibility effect due to the compression and expansion of piston was used. The calculation results of the spray tip penetration and tip velocity using the TAB model showed similar trends comparing with the experimental data. Although the evaporation rate was not nearly affected with the breakup model at the higher injection pressure, in the low injection case, the evaporation rate result using the TAB model became higher than that of R&D model. The evaporation rate was increased with the injection pressure due to the vigorous interaction with the gas field.