• Journal of Korea Soil Environment Society
• /
• v.1 no.1
• /
• pp.89-102
• /
• 1996

An integrated model is presented to describe underground flow and mass transport, using a multicomponent multiphase approach. The comprehensive governing equation is derived considering mass and force balances of chemical species over four phases(water, oil, air, and soil) in a schematic elementary volume. Compact and systemati notations of relevant variables and equations are introduced to facilitate the inclusion of complex migration and transformation processes, and variable spatial dimensions. The resulting nonlinear system is solved by a multidimensional finite element code. The developed code with dynamic array allocation, is sufficiently flexible to work across a wide spectrum of computers, including an IBM ES 9000/900 vector facility, SP2 cluster machine, Unix workstations and PCs, for one-, two and three-dimensional problems. To reduce the computation time and storage requirements, the system equations are decoupled and solved using a banded global matrix solver, with the vector and parallel processing on the IBM 9000. To avoide the numerical oscillations of the nonlinear problems in the case of convective dominant transport, the techniques of upstream weighting, mass lumping, and elementary-wise parameter evaluation are applied. The instability and convergence criteria of the nonlinear problems are studied for the one-dimensional analogue of FEM and FDM. Modeling capacity is presented in the simulation of three dimensional composite multiphase TCE migration. Comprehesive simulation feature of the code is presented in a companion paper of this issue for the specific groundwater or flow and contamination problems.

• Journal of the Korea Concrete Institute
• /
• v.26 no.2
• /
• pp.211-218
• /
• 2014

In the design of reinforced rectangular concrete beam structure, the minimum amount of flexural reinforcement is required to avoid brittle failure. KCI code is based on concept of ultimate strength and usually used as a model code. But bridge design code enacted by Ministry of land, transportation and maritime affairs in 2012 is based on concept of limit state and similar to Euro code EN 1992-2. This means that the minimum reinforcement presented in both design codes has different origination and safety margin. When rectangular concrete beams with minimum reinforcement are designed according to EN and KCI codes, the amount of minimum reinforcement specified in EN code is only 76% of that in KCI code. This makes the design engineers to be confused. In this study, flexural tests were conducted on nine beams with the two different minimum reinforcement specified in KCI and EN design codes. In results, the measured ratios of nominal strength to crack strength from the test were about 25% greater than those evaluated from the equations presented in KCI and EN codes. The EN beams having only 76% of the minimum reinforcement for the KCI beams were fractured by rupture of steel reinforcement but in ductile manner. It is confirmed that the minimum reinforcement concrete beams designed according to both codes have enough safety margin in flexural capacity and moreover in ductility.

• Journal of the Korea Concrete Institute
• /
• v.16 no.3 s.81
• /
• pp.335-344
• /
• 2004

A general composite section of precast and cast-in-place concrete with prestressed and nonprestressed reinforcements was analyzed to calculate residual stresses and loss of prestressing force caused by internal constraints of concrete long-term deformation. From the analytical results, equations to design additional prestressing force and to evaluate resistant moment of the composite section were proposed. The equations shows that the additional prestressing force can be over-estimated if the loss rate of the first prestressing force is over-estimated from the lumped sum of a design code. The analytical procedure with the proposed equations has been applied to a composite section using the AASHTO Type 5 girder. The loss rates of the additional prestressing force appling to the precast concrete girder was less than those appling to the composite girder. However, the resistant moment of the additional prestressing force on the composite girder was much larger than that on the precast concrete girder. The additional prestressing force appling to the composite section was very effective for strengthening of the prestressed concrete composite girder.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.22 no.1
• /
• pp.47-57
• /
• 2010

Subgrid scale stabilization method is applied to Woo and Liu(2004)'s extended Boussinesq FEM numerical model to eliminate the 2dx wiggles. In order to optimize the computational efficiency, Hessian operator is introduced and the matrix of velocity vector is combined to one matrix for solving matrix equations. The mass lumping technique is also applied to the matrix equations of auxiliary variables. The newly developed code is applied to simulate Vincent and Briggs(1989)' wave transformation experiments and the results show that the numerical solution is almost wiggle-free and it matches very well with experimental data. Due to improvement of computational efficiency and wiggle reduction, it is plausible to apply this model to a realistic problem such as harbor oscillation problems.

• The Journal of the Institute of Internet, Broadcasting and Communication
• /
• v.9 no.5
• /
• pp.139-146
• /
• 2009

Parametric array for audio applications is analyzed by numerical modeling and analytical approximation. The nonlinear wave equations are used to provide design guidelines for the audio parametric array. A time domain finite difference code that accurately solves the KZK (Khokhlov-Zabolotskaya-Kuznetsov) nonlinear parabolic wave equation is used to predict the response of the parametric array. The time domain code relates the source size and the carrier frequency to the audible signal response including the output level and beamwidth to considering the implementation issues for audio applications of the parametric array, the emphasis is given to the frequency response and distortion. We use the time domain code to find out the optimal parameters that will help produce the parametric array with highest achievable output in terms of the average power within the demodulated signal. Parameters such as primary input frequency, audio source radius and the modulation method are given utmost importance. The output effect of those parameters are demonstrated through the numerical simulation.

• Journal of the Korean Society of Visualization
• /
• v.15 no.3
• /
• pp.20-26
• /
• 2017

The two-fluid equations are widely used to simulate two-phase flows in a nuclear reactor. For the two-fluid momentum equation, the wall and interfacial drag terms play an important role in predicting a two-phase flow behavior. Since the bubble density is much smaller than the water density, the bubble accelerates faster than the liquid in a nozzle. As a result, the bubble phase becomes faster than the liquid phase in the nozzle. In contrast, the opposite phenomena occur in the diffuser. The purpose of our study is to experimentally show these behaviors in an area-varying channel such as nozzle and diffuser. Experiments were made of turbulent bubbly flows in an area-varying horizontal channel. The velocities of the bubble and liquid phases were measured by the PIV technique. It was shown that the two-phase velocities were no longer close to each other in the area-varying regions. The bubble was faster than the liquid in the nozzle; in contrast, the bubble was slower than the liquid in the diffuser. Code simulations were also performed using the MARS code. By replacing the original wall drag model in the MARS code with Kim (1)'s wall drag partition model, we obtained the simulation results being consistent with experimental observations.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.5
• /
• pp.407-419
• /
• 2008

In this study, the transition Reynolds number of compressible axi-symmetric sharp cone boundary layer is predicted by using a linear stability theory and the -method. The compressible linear stability equation for sharp cone boundary layer was derived from the governing equations on the body-intrinsic axi-symmetric coordinate system. The numerical analysis code for the stability equation was developed based on a second-order accurate finite-difference method. Stability characteristics and amplification rate of two-dimensional second mode disturbance for the sharp cone boundary layer were calculated from the analysis code and the numerical code was validated by comparing the results with experimental data. Transition prediction was performed by application of the -method with N=10. From comparison with wind tunnel experiments and flight tests data, capability of the transition prediction of this study is confirmed for the sharp cone boundary layers which have an edge Mach number between 4 and 8. In addition, effect of wall cooling on the stability of disturbance in the boundary layer and transition position is investigated.

• Journal of the Korea Concrete Institute
• /
• v.14 no.6
• /
• pp.1010-1021
• /
• 2002

The current ACI design code does not take into account the contribution of concrete for the torsional moment of reinforced concrete(RC) beams subjected to pure torsion. This code is not capable of evaluating the inter-effects between concrete and torsional reinforcement on the torsional resistance of the RC beams. Some test results indicated that the current ACI code was not successful in predicting the observed torsional moment of the RC beams with reasonable accuracy. The research reported in this paper provides an evaluation equation to predict the torsional moment of the RC beams subjected to pure torsion. The proposed equation is derived from the equilibrium as well as compatibility equations of the truss model for the cracked RC beams. Comparisons between the observed and calculated torsional moments of the 66 tested beams, showed reasonable agreement.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2008.11a
• /
• pp.839-842
• /
• 2008

The ACI 318-05 code requires the maximum amount of shear reinforcement in reinforced concrete (RC) beams to prevent possible sudden shear failure due to over reinforcement. The design equations of the maximum amount of shear reinforcement provided by the current four design codes, ACI 318-05, CSA-04, EC2-02, and JCI-99, differ substantially from one another. The ACI 318-05, CSA-04, and EC2-02 codes provide an expression for the maximum amount of shear reinforcement ratio as a function of the concrete compressive strength, but Japanese code does not take into account the influence of the concrete compressive strength. For high strength concrete, the maximum amount of shear reinforcement calculated by the EC2-02 and CSA-04 is much greater than that calculated by the ACI 318-05. This paper presents the effects of shear reinforcement ratio and compressive strength of concrete on the maximum shear reinforcement in reinforced concrete beams. Ten RC beams having various shear reinforcement ratio were tested. Although the test beams were designed to have much more amount of shear reinforcement than that required in the ACI 318-05 code, all beams failed due to web concrete crushing after the stirrups reached the yield strain.

• Magazine of the Korea Concrete Institute
• /
• v.7 no.1
• /
• pp.117-125
• /
• 1995

Bond test was carried out to assess the effect of several variables on bond characteristics between reinforcing bar and concrete. Key variables are concrete compressive strength($f_c$'=340, 460, 6513, 904kg/$cm^2$), concrete cover (25, 38, 51, 105, IlOmm), and bar diameter(Dl3 and D22). Confining effect and bar spacing are not taken into account. Thirty-two specimens subjected to uniaxial tension were tested under hypothesis uniform bond stress distribution along the reinforcing bar embeded in concrete. Failure mode was examined and local bond stress versus slip relationship diagrams were represented to show effect of the above variables, also test results(u1timate bond stress) were compared with bond and development provisions of the ACI Building Code(AC1 318-89) and proposed equations from previous research. According to analysis, borld stress and ultimate bond stress increased although compressive strength increased beyond the ACI Building Code upper limit. Therefore in calculating development length. compressive strength effect(exceeding 700kg /$cm^2$) should be accounted.