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

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.05a
• /
• pp.309-313
• /
• 2004

Prediction of green water loads acting an the bow deck is an essential part for the design of bow structures against the green water impact. Proper technique of the green water simulation is highly required for the prediction of green water loads. In this paper, a new numerical method for green water simulation, which is based an predictor-corrector-upwind finite difference scheme of the 2nd kind, is introduced. Through the comparisons between computed' results and experimental measurements, it is verified that the present numerical tool is adequate as a practical calculation tool for the green water problem.

• International Journal of Automotive Technology
• /
• v.1 no.1
• /
• pp.35-41
• /
• 2000

The crashworthiness of vehicles with finite element methods depends on the geometry modeling and the material properties. The vehicle body structures are generally composed of various members such as frames, stamped panels and deep-drawn parts from sheet metals. In order to ensure the impact characteristics of auto-body structures, the dynamic behavior of sheet metals must be examined to provide the appropriate constitutive relation. In this paper, high strain-rate tensile tests have been carried out with a tension type split Hopkinson bar apparatus specially designed for sheet metals. Experimental results from both static and dynamic tests with the tension split Hopkinson bar apparatus are interpolated to construct the Johnson-Cook and a modified Johnson-Cook equation as the constitutive relation, that should be applied to simulation of the dynamic behavior of auto-body structures. Simulation of auto-body structures has been carried out with an elasto-plastic finite element method with explicit time integration. The stress integration scheme with the plastic predictor-elastic corrector method is adopted in order to accurately keep track of the stress-strain relation for the rate-dependent model accurately. The crashworthiness of the structure with quasi-static constitutive relation is compared to the one with the rate-dependent constitutive model. Numerical simulation has been carried out for frontal frames and a hood of an automobile. Deformed shapes and the Impact energy absorption of the structure are investigated with the variation of the strain rate.

• Wind and Structures
• /
• v.26 no.1
• /
• pp.25-34
• /
• 2018

This paper presents a beam finite element model of a vibrate wind blade in large elastic deformation subjected to the aerodynamic, centrifugal, gyroscopic and gravity loads. The gyroscopic loads applied to the blade are induced by her simultaneous vibration and rotation. The proposed beam finite element model is based on a simplex interpolation method and it is mainly intended to the numerical analysis of wind blades vibration in large elastic deformation. For this purpose, the theory of the sheared beams and the finite element method are combined to develop the algebraic equations system governing the three-dimensional motion of blade vibration. The applicability of the theoretical approach is elucidated through an original case study. Also, the static deformation of the used wind blade is assessed by appropriate software using a solid finite element model in order to show the effectiveness of the obtained results. To simulate the nonlinear dynamic response of wind blade, the predictor-corrector Newmark scheme is applied and the stability of numerical process is approved during a large time of blade functioning. Finally, the influence of the modified geometrical stiffness on the amplitudes and frequencies of the wind blade vibration induced by the sinusoidal excitation of gravity is analyzed.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.8 no.3
• /
• pp.116-121
• /
• 2005

In the present paper, a new calculation algorithm far solving large scale environmental or geophysical flows with free surface is proposed where the non-hydrostatic pressure component is taken into consideration. Predictor-corrector fractional step approach with explicit, forward time marching scheme in the sigma coordinate system is employed. In order to validate the present calculation algorithm and to estimate the effects of non-hydrostatic pressure on resultant flow and free surface movements, example calculations are carried out for typical steady and unsteady flow problems. Present method can be applied to the meso-scale free surface flows with complex bottom topography where MAC-like 3-d hydrodynamic calculations are quite ineffective and uneconomic.

• Ocean Systems Engineering
• /
• v.4 no.2
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.8
• /
• pp.1087-1096
• /
• 2001

The dynamic structures of unsteady CH$_4$/Air jet diffusion flame with a flame-vortex interaction were numerically investigated. A timed-dependent, axisymmetric computational model and a low mach number approximation were employed in the present calculation. A two-step global reaction mechanism which considers 6 species, was used to calculate the reaction rates. The predicted results including the gravitational effect show that the large outer vortices and the small inner vortices can be well simulated without any additional disturbances near nozzle tip. It was found that the temperature and species concentrations have deviated values even for the same mixture fraction in the flame-vortex interaction region. It was also shown that the flame surface is not deformed by the inner vortex in upstream region, while in downstream region, the flame surface is compressed or stretched by the outer vortex roll-up. The present unsteady jet flame configuration accompanying a flame-vortex interaction is expected to give good implications for the unsteady structures of turbulent flames.

• Korean Journal of Soil Science and Fertilizer
• /
• v.36 no.4
• /
• pp.181-192
• /
• 2003

We developed a mathematical simulation model to portray the vertical distribution of soil water from the measured weather data and the known soil hydraulic properties, and then compared simulation results with the periodically measured soil water profiles obtained on Jungdong sandy loam to verify the model, In this model, we solved potential-based Richards' equation by the implicit finite difference method superimposed on the predictor-corrector scheme. We presumed that: soil hydraulic properties are homogeneous; soil water flows isothermally; hysteresis is not considered; no vapor flows; no heat transfers into the soil profiles; and water added to soil surface is distributed along the soil profile following partial displacement principle. The input data were broadly classified into two groups: (1) daily weather data such as rainfall, maximum and minimum air temperatures, relative humidity and solar radiation and (2) soil hydraulic data to approximate unsaturated hydraulic conductivity and water retention. Each hydraulic polynomial function approximated using the Chebyshev polynomial and least square difference technique in tandem showed a fairly good fit of the given set of data. Vertical distribution of soil water as approximations to the Richards' equation subject to changing surface and phreatic boundaries was solved numerically during 53 days with a comparatively large time increment, and this pattern agreed well with field neutron scattering data, except for the surface 0.1 m slab.

• Ocean Systems Engineering
• /
• v.10 no.2
• /
• pp.201-226
• /
• 2020

A systematic numerical comparative study of the performance of semicircular and rectangular submerged breakwaters interacting with solitary waves is the basis of this paper. To accomplish this task, Nwogu's extended Boussinesq model equations are employed to simulate the interaction of the wave with breakwaters. The finite difference technique has been used to discretize the spatial terms while a fourth-order predictor-corrector method is employed for time discretization in our numerical model. The proposed computational scheme uses a staggered-grid system where the first-order spatial derivatives have been discretized with fourth-order accuracy. For validation purposes, five test cases are considered and numerical results have been successfully compared with the existing analytical and experimental results. The performances of the rectangular and semicircular breakwaters have been examined in terms of the wave reflection, transmission, and dissipation coefficients (RTD coefficients) denoted by K_R, K_T, K_D. The latter coefficient K_D emerges due to the non-energy conserving K_R and K_T. Our computational results and graphical illustrations show that the rectangular breakwater has higher reflection coefficients than semicircular breakwater for a fixed crest height, but as the wave height increases, the two reflection coefficients approach each other. un the other hand, the rectangular breakwater has larger dissipation coefficients compared to that of the semicircular breakwater and the difference between them increases as the height of the crest increases. However, the transmission coefficient for the semicircular breakwater is greater than that of the rectangular breakwater and the difference in their transmission coefficients increases with the crest height. Quantitatively, for rectangular breakwaters the reflection coefficients K_R are 5-15% higher while the diffusion coefficients K_D are 3-23% higher than that for the semicircular breakwaters, respectively. The transmission coefficients K_T for rectangular breakwater shows the better performance up to 2.47% than that for the semicircular breakwaters. Based on our computational results, one may conclude that the rectangular breakwater has a better overall performance than the semicircular breakwater. Although the model equations are non-dissipative, the non-energy conserving transmission and reflection coefficients due to wave-breakwater interactions lead to dissipation type contribution.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.10
• /
• pp.28-37
• /
• 2004

A two-dimensional direct numerical simulation was performed to investigate the evolution and vortical structure of a single vortex in reacting and non-reacting jet flow fields. A predictor-corrector-type numerical scheme with a low Mach number approximation was used, and a two-step global reaction mechanism was adopted as the combustion model. Through the comparisons of single vortex behaviors in reacting and non-reacting jet flow fields, it was found that the evolution characteristics and vortical structure of the single vortex were significantly influenced by a outer vortex that was generated from the buoyance effect as well as the chemical heat release. Furthermore, it was also identified that the differences of the vortical structure in reacting and non-reacting jet flow fields were mainly attributed to the thermal expansion, Baroclinic torque and buoyance effect.