• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.4
• /
• pp.609-620
• /
• 1992

A numerical solution of temperature and thermally induced stress in a wafer during rapid thermal processing (R.T.P) is obtained, and an analysis of onset and propagation of slip is performed and compared with experiment. In order to calculate temperature distribution of a wafer in R.T.P system, heat conduction equation that incorporated with radiative and convective heat transfer model is solved, and the solution of the equation is calculated numerically using alternating direction implicit (A.D.I) method. In dealing with radiative heat transfer, a partially transparent body that absorbs the radiation energy is assumed and this transparent body undergoes multiple internal reflections and absorptions. Two dimensional (assuming plane stress) thermoelastic constitutive equation is used to calculate thermal stress induced in a wafer and finite element method is employed to solve the equation numerically. The stress resolved in the slip directions on the slip planes of silicon is compared with the yield stress of silicon in order to predict the slip. The result of the analysis shows that the wafer temperature at which slip occurs is affected by the heating rate of the R.T.P system. It is observed that once slip occurs in the wafer, the slip grows.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.7
• /
• pp.1-8
• /
• 2005

For study on the unsteady wall interference effect, flows around a circular cylinder in closed test-section wind tunnels have been numerically investigated by solving compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with the Roe's flux-difference splitting and an implicit time-integration method coupled with dual time-step sub-iteration. The computed results showed that the unsteady pressure gradient over the cylinder is enhanced by the wall interference, and as a result the fluctuations of lift and drag are augmented. The drag is further increased because of the lower base pressure. The vortex shedding frequency is also increased by the wall interference. The pressure on the test section wall shows the harmonics having the shedding frequency contained in the wall effect.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1487-1503
• /
• 2019

The methods and performance of a pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT employing the drift-flux model are presented. This code aims at providing an accurate yet fast core thermal-hydraulics solution capability to high-fidelity multiphysics core analysis systems targeting massively parallel computing platforms. The four equation drift-flux model is adopted for two-phase calculations, and numerical solutions are obtained by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models involving turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analyses. ESCOT is parallelized by a domain decomposition scheme that involves both radial and axial decomposition to enable highly parallelized execution. The ESCOT solutions are validated through the applications to various experiments which include CNEN $4{\times}4$, Weiss et al. two assemblies, PNNL $2{\times}6$, RPI $2{\times}2$ air-water, and PSBT covering single/two-phase and unheated/heated conditions. The parameters of interest for validation include various flow characteristics such as turbulent mixing, spacer grid pressure drop, cross-flow, reverse flow, buoyancy effect, void drift, and bubble generation. For all the validation tests, ESCOT shows good agreements with measured data in the extent comparable to those of other subchannel-scale codes: COBRA-TF, MATRA and/or CUPID. The execution performance is examined with a mini-sized whole core consisting of 89 fuel assemblies and for an OPR1000 core. It turns out that it is about 1.5 times faster than a subchannel code based on the two-fluid three field model and the axial domain decomposition scheme works as well as the radial one yielding a steady-state solution for the OPR1000 core within 30 s with 104 processors.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.8 no.3
• /
• pp.23-27
• /
• 2008

The new conjunctive surface-subsurface flow model at a large scale was developed by using a 1-D Diffusion Wave (DW) model for surface flow interacting with the 3-D Volume Averaged Soil-moisture Transport (VAST) model for subsurface flow for the comprehensive terrestrial water and energy predictions in Land Surface Models (LSMs). A selection of numerical implementation schemes is employed for each flow component. The 3-D VAST model is implemented using a time splitting scheme applying an explicit method for lateral flow after a fully implicit method for vertical flow. The 1-D DW model is then solved by MacCormack finite difference scheme. This new conjunctive flow model is substituted for the existing 1-D hydrologic scheme in Common Land Model (CLM), one of the state-of-the-art LSMs. The new conjunctive flow model coupled to CLM is tested for a study domain around the Ohio Valley. The simulation results show that the interaction between surface flow and subsurface flow associated with the flow routing scheme matches the runoff prediction with the observations more closely in the new coupled CLM simulations. This improved terrestrial hydrologic module will be coupled to the Climate extension of the next-generation Weather Research and Forecasting (CWRF) model for advanced regional, continental, and global hydroclimatological studies and the prevention of disasters caused by climate changes.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.9 no.1
• /
• pp.9-23
• /
• 1997

A numerical study of a two-layer, stratified flow is investigated, using the implicit finite difference method in one dimension. The results of computational method have been tested and, in case of lock exchange flow, compared with the results of experimental data. The results of model experiments with various interfacial, bottom friction coefficients along with various time weighting factor of numerical scheme and dissipative interface are shown and discussed. Two-layer model experiment has been also carried out to investigate the generation and propagation characteristics of internal tidal wave over the steep bottom topography under stratified condition. The internal wave seems to well radiate through the downstream boundary under the experiments adopting radiation conditions both at two layers and only at upper layer, confirming the applicability of radiational boundary condition in stratified flows. It is also shown that the internal wave through the downstream boundary propagates more actively with increasing thickness of lower layer in the downstream. This implies that the potential tidal energy in the interface will depend upon the thickness of lower layer for the constant thickness of upper layer.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.7 no.4
• /
• pp.185-191
• /
• 2004

A mooring system can be applied to keep the position of a floating structures. In this study, the structural analysis is carried out to analyze the dynamic characteristics of a mooring line for a floating breakwater. A three-dimensional equations of motion for a submerged chain are derived. Bending stiffness is considered for the necessary restoring force in the regions of zero tension. A fortran program is to be developed by employing finite difference method. In the algorithm, an implicit time integration and Newton-Raphson iteration are adopted. The results of simulation show good agreement in tension response pattern with the experimental results of a reference. The results of this study can contribute for the design of mooring system for a floating breakwater.

• Journal of the Korea Computer Graphics Society
• /
• v.8 no.1
• /
• pp.37-45
• /
• 2002

We present a cloth simulation technique that is very stable yet also responsive. The stability of the technique allows the use of a large fixed time step when simulating various types of fabrics and character motions. The animations generated using this technique are strikingly realistic. Wrinkles form and disappear in a quite natural way, which is the feature that most distinguishes textile fabrics from other sheet materials. Significant improvements in both the stability and realism were made possible by overcoming the post-buckling instability as well as the numerical instability. The instability caused by buckling arises from a structural instability and therefore cannot be avoided by simply employing an implicit method. Addition of a damping force may help to avoid instabilities; however, it can significantly degrade the realism of the cloth motion. In this paper, a new buckling model based on immediate buckling assumption is proposed. A cloth element is assumed to reach a stable configuration immediately once it begins to buckle. This assumption makes it possible to simulate the fabric buckling stably without introducing any fictitious damping force. Consequently, it produces highly responsive cloth motion as well as improves the stability by modeling the fabric-specific buckling property adequately.

• Journal of Korea Water Resources Association
• /
• v.42 no.7
• /
• pp.571-578
• /
• 2009

Levee-break Analysis model is developed to predict the variation of breach width according to time and to estimate inundation area and depth in protected lowland. This Model calculate flood depth using 4 point implicit finite difference method in river channel and analyze breach flow based on physical theory introducing soil transport equation and erosion process. Breach analysis model and channel flood model are combined into Levee-Break Model and this model is applied to actual levee break case. Then, this model can simulate reasonably many levee-break parameters such as river stage, breach width, breach formation and so on. If the applicability of this model is proved through applications to more various actual levee-break cases, the suggested model is expected to do more accurate flood analyses on levee break site.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.1
• /
• pp.17-23
• /
• 2014

Graphics processing units (GPUs) are ideal for solving problems involving parallel data computations. In this study, the GPU is used for effectively carrying out a multi-body dynamic simulation with particle dynamics. The Hilber-Hushes-Taylor (HHT) implicit integration algorithm is used to solve the integral equations. For detecting collisions among particles, the spatial subdivision algorithm and discrete-element methods (DEM) are employed. The developed program is verified by comparing its results with those of ADAMS. The numerical efficiencies of the serial program using the CPU and the parallel program using the GPU are compared in terms of the number of particles, and it is observed that when the number of particles is greater, more computing time is saved by using the GPU. In the present example, when the number of particles is 1,300, the computational speed of the parallel analysis program is about 5 times faster than that of the serial analysis program.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.31 no.1
• /
• pp.117-130
• /
• 1989

In this study, an implicit one-dimensional model, DWRM(Dynamic Wave Routing Model) was developed by using the four-point weighted difference method. By applying the developed model to the Keum River, the parameters were calibrated and the model applicability was tested through the comparison between observed and computed water levels. In addition, the effects of the construction of an estuary dam to the flood wave were estimated as a result of the model application. The results of the study can be summarized as follows; 1. The roughness coefficients were evaluated by comparison between observed and computed water level at Jindu, Gyuam and Ganggyeung station in 1985. The Root Mean Squares for water level differences between observed and computed values were 0.10, 0.11, 0. 29m and the differences of peak flood levels were 0.07, 0.02, 0. 07m at each station. Since the evaluated roughness coefficients were within the range of 0.029-0.041 showing the realistic value for the general condition of rivers, it can be concluded that the calibration has been completed. 2. By the application of model using the calibrated roughness coefficients, the R. M. S. for water level differences were 0.16, 0.24, 0. 24m and the differences of peak flood level were 0.17, 0.13,0.08 m at each station. The arrival time of peak flood at each station and the stage-discharge relationship at Gongju station agreed well with the observed values. Therefore, it was concluded that the model could be applied to the Keum River. 3. The model was applied under conditions before and after the construction of the estuary dam. The 50-year frequency flood which had 7, 800m$^3$/sec of peak flood was used as the upstream condition, and the spring tide and the neap tide were used as the downstream condition. As the results of the application, no change of the peak flood level was showed in the upper reaches of 19.2km upstream from the estuary dam. For areas near 9.6km upstream from the estuary dam, the change of the peak flood level under the condition before and after the construction was 0. 2m. However considering the assumptions for the boundary conditions of downstream, the change of peak flood level would be decreased.