• Geotechnical Engineering
• /
• v.14 no.5
• /
• pp.191-204
• /
• 1998

The feasibility of fly ash was evaluated as an alternative liner material to the conventional clay liner of landfills through modeling and laboratory experiments. In order to consider the effect of unsaturation on water flow through the liner, analyses were made to compare flow characteristics in saturated liner with that of unsaturated one. Contaminant migration characteristics in liners were investigated by batch experiment and modeling, in which phenol was employed as a model was solved by numerical techniques of finite difference method and predictor-corrector method to deal with high non-linearity. Sequential method was used to handle the system of differential equations. Results show that the alternative liner material is more capable of cutting off water flow in unsaturated condition and in preventing phenol from passing through it. It can be seen that, under the flow conditions considered in this study, the conventional saturation approach underestimates the amount of water passing through the liner and doers the cut-off capability against phenol significantly.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.3
• /
• pp.725-733
• /
• 1990

An analysis of laminar mixed convection flow adjacent to the inclined flat surface which is subjected to a uniform temperature in a uniform free stream is performed. Nonsimilar boundary layed equation are derived by using the mixed convection parameters such that smooth transition from the purely forced convection limit to the purely free convection limit is possible. The governing equations are solved by a finite difference method using the coupled box scheme of sixth order. Numerical results are presented for prandtl numbers of 0.7 and 7 with the angle of inclination ranging from 0 to 90 degree from the vertical. The velocity distributions for the buoyancy assisting flow exhibit a significant overshoot above the free stream value in the region of intense mixed convection and the velocity field is found to be more sensitive to the buoyancy effect than the temperature field. The separation point near the wall was obtained for the buoyancy opposing flow. The local Nusselt number increases for buoyancy assisting flow and decreases for opposing flow with increasing value of the local Grashoff number in the mixed convection parameter. For large Prandtl number, the Nusselt number and the friction factor decrease significantly near the separation point. Present numerical predictions are in good agreement with recent experimental results by Ramachandran.

• Journal of Korean Society of Steel Construction
• /
• v.12 no.3 s.46
• /
• pp.291-302
• /
• 2000

In recent years the development of high modulus, high strength and low density boron and graphite fibers bonded together has brought renewed interestes in structural elements. When a plate with arbitrarily oriented layers and clamped boundary conditions is subjected to uniform loading, it is difficult to analyze and apply, compared with isotropic and orthotropic cases. Therefore the numerical methods, such as finite difference method or finite element method, should be emloyed to analyse such problems. In this study the finite difference technique is used to formulate the bending analysis of symmetric composite laminated skew plates. When this technique is used to solve the problem, it is desirable to reduce the order of the derivatives in order to minimize the number of the pivotal points involved in each equation. The 4th order partial differential equations of laminated skew plates are converted to an equivalent three of 2nd order partial differential equations with three dependant variables.

• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.5A
• /
• pp.255-266
• /
• 2012

Main cable of a suspension bridge is the important member which shows the overall structure integrity at bridge completion. Configuration of main cable is a free hanging state at cable erection completion and is different from that at bridge completion supporting the dead loads such as hanger, girder, and so on. Accordingly, the configuration control under cable erection is considerably significant because the configuration at cable erection completion has direct influence on that at bridge completion. That is performed by sag adjustments at center, side span and tension adjustments at anchor span. The former needs the sag sensitivity which represents the control quantity of strand length corresponding to that of sag. The latter requires the tension sensitivity which shows the change of strand tension according to that of strand temperature. In this study, the fundamental equations of cable were derived with the assumption of either catenary or parabola shape, the differential-related equations using chain rule on horizontal tension were drawn from those and finally the estimation methods of the sag / tension sensitivity were proposed from both those. The nonlinear numerical analysis flow charts of sag sensitivity based on the catenary equations were proposed and the sag sensitivities grounded on the differential-related equations were compared with the results using them for various parameters of sag change. Also, considering the combinations of sag change parameters, the calculation method of the final variation for the cable sag was suggested. For the real suspension bridge under construction with PPWS method, the sag/tension sensitivity were estimated considering the construction conditions like the change of PPWS length, PPWS temperature, bridge span, etc.. We hope that this study will be a systematic guideline for the configuration control under main cable erection and improved highly by field verification in the real bridge site.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.31 no.4
• /
• pp.372-378
• /
• 1995

This study describes an analysis on the sinking characteristics in purse seine. The experiment was carried out using three simplified model seines in a flume tank under still water condition. The densities($\rho$) of netting materials were 0.91g/$cm^3$ for PP seine, 1.14g/$cm^3$ for PA seine and 1.38g/$cm^3$ for PES seine. Differential equations were derived from the conservation of momenta of the model seines and used to determine the sinking speeds of the depths of leadline and the other portions of seines. An analysis carried out by simultaneous differential equations for numerical method by subroutine Runge-Kutta-Gill. The results obtained were as follows: 1) Sinking speed of net margin was fastest for PP seine, followed by PA and PES seines. 2) The coefficient of resistance for netting of seine was estimated to be $K_D=0.061({\frac{\rho}{{\rho}_w}})^4$. 3) The coefficient of resistance for netting bundle of seine was estimated to be $C_R=0.91({\frac{\rho}{{\rho}_w}}$. 4) In all seines, the calculated depths of leadline closely agreed with the measured ones(meas.=0.99cal.).

• Journal of the Korean Institute of Gas
• /
• v.7 no.4 s.21
• /
• pp.44-52
• /
• 2003

In this work dynamic heat transfer in a CPFS (cable penetration fire stop) system built in the firewall of nuclear power plants is three-dimensionally investigated to develop a test-simulator that can be used to verify effectiveness of the sealant. Dynamic heat transfer in the fire stop system is formulated in a parabolic PDE (partial differential equation) subjected to a set of initial and boundary conditions. First, the PDE model is divided into two parts; one corresponding to heat transfer in the axial direction and the other corresponding to heat transfer on the vertical planes. The first PDE is converted to a series of ODEs (ordinary differential equations) at finite discrete axial points for applying the numerical method of SOR (successive over-relaxation) to the problem. The ODEs are solved by using an ODE solver In such manner, the axial heat flux can be calculated at least at the finite discrete points. After that, all the planes are separated into finite elements, where the time and spatial functions are assumed to be of orthogonal collocation state at each element. The initial condition of each finite element can be obtained from the above solution. The heat fluxes on the vertical planes are calculated by the Galerkin FEM (finite element method). The CPFS system was modeled, simulated, and analyzed here. The simulation results were illustrated in three-dimensional graphics. Through simulation, it was shown clearly that the temperature distribution was influenced very much by the number, position, and temperature of the cable stream, and that dynamic heat transfer through the cable stream was one of the most dominant factors, and that the feature of heat conduction could be understood as an unsteady-state process.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.2
• /
• pp.153-163
• /
• 2010

In this paper a frequency response analysis using Krylov subspace-based model reduction and its design sensitivity analysis with respect to design variables are presented. Since the frequency response and its design sensitivity information are necessary for a gradient-based optimization, problems of high computational cost and resource may occur in the case that frequency response of a large sized finite element model is involved in the optimization iterations. In the suggested method model order reduction of finite element models are used to calculate both frequency response and frequency response sensitivity, therefore one can maximize the speed of numerical computation for the frequency response and its design sensitivity. As numerical examples, a semi-monocoque shell and an array-type $4{\times}4$ MEMS resonator are adopted to show the accuracy and efficiency of the suggested approach in calculating the FRF and its design sensitivity. The frequency response sensitivity through the model reduction shows a great time reduction in numerical computation and a good agreement with that from the initial full finite element model.

• Geophysics and Geophysical Exploration
• /
• v.14 no.2
• /
• pp.164-175
• /
• 2011

Using natural electromagnetic (EM) fields at low frequencies, magnetotelluric (MT) surveys can investigate conductivity structures of the deep subsurface and thus are used to explore geothermal energy resources and investigate proper sites for not only geological $CO_2$ sequestration but also enhanced geothermal system (EGS). Moreover, marine MT data can be used for better interpretation of marine controlled-source EM data. In the interpretation of MT data, MT modeling schemes are important. This study improves a three dimensional (3D) MT modeling algorithm which uses edge finite elements. The algorithm computes magnetic fields by solving an integral form of Faraday's law of induction based on a finite difference (FD) strategy. However, the FD strategy limits the algorithm in computing vertical magnetic fields for a topographic model. The improved algorithm solves the differential form of Faraday's law of induction by making derivatives of electric fields, which are represented as a sum of basis functions multiplied by corresponding weightings. In numerical tests, vertical magnetic fields for topographic models using the improved algorithm overcome the limitation of the old algorithm. This study recomputes induction vectors and tippers for a 3D hill and valley model which were used for computation of the responses using the old algorithm.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.18 no.4
• /
• pp.308-320
• /
• 2006

An analytic solution to the mild-slope equation is derived for waves propagating over an axi-symmetric pit. The water depth inside the pit varies in proportion to a power of radial distance from the pit center. The governing equation is transformed into ordinary differential equations by using separation of variables, and the coefficients of the equations are transformed into explicit forms by using Hunt's (1979) approximate solution. Finally, by using the Frobenius series, the analytic solution is derived. Due to the feature of Hunt's equation, the present analytic solution is accurate in shallow and deep waters, while it is less accurate in intermediate depth water. The validity of the analytic solution is demonstrated by comparison with numerical solutions. The analytic solution is also used to examine the effects of pit geometry and relative depth on wave transformation.