• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.5-10
• /
• 2003

In order to reduce the excessive numerical dissipation, the new spatial discretization scheme is introduced. The present method in this paper has the formula that has an additional procedure of defining transferred properties at a cell-interface, based on AUSMPW+. The newly defined transferred property could eliminate numerical dissipation effectively in non-flow aligned grid system. In addition, the present method guarantees the monotonic characteristic in capturing a discontinuity. Through a stationary or moving contact discontinuity and a stationary or moving shock discontinuity, a vortex discontinuity and shock wave/ boundary layer interaction, it is verified that the accuracy of the present method is improved.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.2
• /
• pp.69-80
• /
• 2002

Sequential process of roll-bending and line heating has been used to deform the curved hull-plates in shipyards. A growing interest for the mechanization or automation of the line heating process has been noted. Relations between heating conditions and residual deformations are important components needed for the mechanization. The residual deformations are investigated by using a thermal elastic-plastic analysis based on the finite element analysis(FEA). Several experiments are also performed to examine the validity of the results of FEA. The input parameters of line heating are suggested by dimensional analysis of line heating. The dimensional analysis can extract the primary input-parameters of line heating. The relations between the heating conditions and the residual deformations are set up by multi-variate analysis and multiple-regression method. This study suggests a method for the relation between the heating conditions and the deformations lying under the line heating.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.290-297
• /
• 2009

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second order scheme.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2012.05a
• /
• pp.278-281
• /
• 2012

Numerical studies were performed to investigate the three-dimensional front structure and dynamics features of detonation wave propagating in a circular tube such as Pulse Detonation Engine (PDE). By carrying out a series of parametric study using one step irreversible Arrhenius kinetics model, mechanisms of the three-dimensional front structure were investigated for two-, three-, four and six-cell mode detonations. A comparison with two-dimensional results, the effects of slapping transverse waves in radial direction were confirmed. In the all muti-cell modes, the detonation front structures and smoked-records on the wall are formed by the propagation of transverse waves along the wall in clockwise and counter-clockwise while the slapping move in radial direction. And the strength of reflected waves on the curved wall is changed by the multi-dimensional confinement effect.

• Geotechnical Engineering
• /
• v.1 no.1
• /
• pp.59-72
• /
• 1985

This paper deals with the numerical analysis by the (mite element method introducing Biot's theory of consolidation and the modified Cambridge model proposed by Roscoe school of Cambridge University as constitutive equation and using Christian-Boehner's technique. Especially, time interval and division of elements are investigated in vies of stability and economics. In order to check the validity of author's program, the program was tested with one-dimensional consolidation case followed by Terzaghi's exact solution and with the results of the Magnan's analysis for existing banking carried out for study at Cubzac-les-ports in France. The main conclusions obtained are summarized as follows: 1. In the case of one-dimensional consolidation, the more divided the elements are near the surface of the foundation, the higher the accuracy of the numerical analysis is. 2. For the time interval, it is stable to divide 20 times per 1-lg cycle. 3. At the element which has long drain distance, the Mandel-fryer effect appears due to time lag. 4. Lateral displacement at an initial loading stage predicted by author's program, in which the load was assumed as not concentrative. but rather in grid form, is well consistent with the value of observation. 5. The pore water pressure predicted by author's program has a better accordance with the value of observation compared with Magnan's results. 6. Optimum construction control by Matsuo-Kawamura's method is possible with the predicted lateral displacement and settlement by the program.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.10 no.6
• /
• pp.155-163
• /
• 2010

The past researches on flood inundation simulation mainly focused on development of numerical models based on unstructured mesh networks to improve model performances. However, despite the accurate simulation results, such models are not suitable for real-time flood inundation forecasting due to a huge computational burden in terms of geographic data processing. In addition, even though various types of vector and raster data are available to be compatible with flood inundation models for post-processes such as flood hazard mapping and flood inundation risk analysis, the unstructured mesh-based models are not effective to fully use such information due to data incommensurability. Therefore, this study aims to develop a raster-based two-dimensional inundation model; it guarantees computational efficiency because of direct application of DEM for flood inundation modeling and also has a good compatibility with various types of raster data, compared to a commercial model such as FLUMEN. We applied the model to simulate the BaekSan levee break in the Nam river during a flood period from August 10 to 13, 2002. The simulation results showed a good agreement with the field-surveyed inundation area and were also very similar with results from the FLUMEN. Moreover, the model provided physically-acceptable velocity vectors with respect to inundating and returning flows due to the difference of water level between channel and lowland.

• Journal of Energy Engineering
• /
• v.25 no.1
• /
• pp.69-75
• /
• 2016

Bubble size is a key parameter for an accurate prediction of bubble behaviours in the multi-dimensional two-phase flow. In the current STAR CCM+ CFD code, a mechanistic bubble size model $S{\gamma}$ is available for the prediction of bubble size in the flow channel. As another model, Yun model is developed based on DEBORA that is subcooled boiling data in high pressure. In this study, numerical simulation for the gas-liquid two-phase flow was conducted to validate and confirm the performance of $S{\gamma}$ model and Yun model, using the commercial CFD code STAR CCM+ ver. 10.02. For this, local bubble models was evaluated against the air-water data from DEDALE experiments (1995) and Hibiki et al. (2001) in the vertical pipe. All numerical results of $S{\gamma}$ model predicted reasonably the two-phase flow parameters and Yun model is needed to be improved for the prediction of air-water flow under low pressure condition.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.43 no.4
• /
• pp.469-476
• /
• 2023

Soil plasticity models for cyclic and dynamic loads are essential in non-linear numerical analysis of geotechnical structures. While a single yield surface model shows a linear behavior for cyclic loads, J₂-bounding surface plasticity model with zero elastic region can effectively simulate a nonlinearity of the ground response with the same material properties. The radius of the yield surface inside the boundary surface converged to 0 to make the elastic region disappear, and plastic hardening modulus and dilatancy define plastic strain increment. This paper presents the stress-strain incremental equation of the developed model, and derives plastic hardening modulus for the hyperbolic model. The comparative analyses of the triaxial compression test and the shallow foundation under the cyclic load can show stable numerical convergence, consistency with the theoretical solution, and hysteresis behavior. In addition, plastic hardening modulus for the modified hyperbolic function is presented, and a methodology to estimate model variables conforming 1D equivalent linear model is proposed for numerical modeling of the multi-dimensional behavior of the ground.

• Journal of Soil and Groundwater Environment
• /
• v.13 no.5
• /
• pp.33-46
• /
• 2008

A series of three-dimensional numerical simulations using a generalized multidimensional hydrodynamic dispersion numerical model is performed to simulate effectively and to evaluate quantitatively impacts of fault existence on densitydependent groundwater flow and salt transport in coastal aquifer systems. A series of steady-state numerical simulations with calibration is performed first for an actual coastal aquifer system which contains a major fault. A series of steadystate numerical simulations is then performed for a corresponding coastal aquifer system which does not have such a major fault. Finally, the results of both numerical simulations are compared with each other and analyzed. The results of the numerical simulations show that the major fault produces hydrogeologically significant heterogeneity and true anisotropy in the actual coastal aquifer system, and density-dependent groundwater flow, salt transport, and seawater intrusion patterns in the coastal aquifer systems are intensively and extensively dependent upon the existence or absence of such a major fault. Especially, the major fault may act as a pathway for groundwater flow and salt transport along the direction parallel to its plane, while it may also behave as a barrier against groundwater flow and salt transport along the direction perpendicular to its plane.

• Nuclear Engineering and Technology
• /
• v.24 no.3
• /
• pp.274-284
• /
• 1992

A 1/15-scale CREARE experiment, which simulates the thermal-hydraulic behavior in the reactor pressure vessel of a PWR during a hypothetical Loss Of Coolant Accident, has been analyzed using CATHARE code for the associated model assessment to represent the phenomenon. The key parameters examined in the CREARE experiment were known as ECC water injection rate. ECC water subcooling, system pressure, and steam flow rate coming out from the core bottom. The present CATHARE simulation, however, has been mainly focused on qualitative analysis of a countercurrent flow in the downcomer. The discrepancy of the simulation results with the experimental data is considered arising primarily from an inadequate numerical representation as well as an interfacial friction model. Accordingly it is suggested from the sensitivity studies that either multidimensional approach or further examination of momentum equations at a junction near a volume element in CATHARE be necessary in order to represent the phenomenon more realistically.