• Journal of Korea Water Resources Association
• /
• v.41 no.8
• /
• pp.773-784
• /
• 2008

The present study was carried out to examine flow properties in linear drainage channels such as road surface drainage facilities. The finite difference formulation for the varied flow analysis was solved for flow profiles in the channels. Starting the first step at the control section, the Newton-Raphson method was applied for producing numerical solutions of the equation. We considered two types of linear drainage channels, a channel with one outlet at downstream end and a channel with two outlets at both ends. Moreover, the flow analysis for various channel slopes was performed. However, we considered channels with the two outlets of slopes satisfying the condition that the both ends are the control section. The maximum of those slopes was decided from the relation between the channel slope and the location of control section. The flow of a channel with one outlet was calculated upward and downward from the control section existing in channel or upward from the control section at downstream end. The flow of a channel with two outlets at both ends were calculated for upstream and downstream channel segments divided by the water dividend, respectively and the flow analysis was completed when the water depth at the water dividend calculated from upstream end was equal to that calculated from downstream end. If the slope was larger than the critical slope, the channel with two outlets was likely to behave like the channel with one outlet. The maximum water depth was investigated and compared with that calculated additionally from the uniform flow analysis. The uniform flow analysis was likely to lead a excessive design of a drainage channel with mild slope.

• Clean Technology
• /
• v.10 no.1
• /
• pp.9-17
• /
• 2004

Direct methanol fuel cell(DMFC) with proton exchange membrane (PEM) has advantages over the conventional power source (e.g. vehicle). DMFC, however, has a problem to be solved such as methanol crossover, high anodic overpotential and limiting current density, etc. The physicochemical phenomena in DMFC can be described by coupled PDEs (partial differential equations), which can be solved by a PDE solver. In this paper, we utilized a commercial software FEMLAB to solve the PDEs. The FEMLAB is one of the software programs available which are developed as a solver for building physics problems based on PDEs and is designed to simulate systems of coupled PDEs which may be 1D, 2D, 3D, non-liner and time dependent. We performed simulation using the Tafel equation as an electrochemical reaction model to analyze methanol concentration profile in DMFC system. We confirm that the rapid decrease of methanol concentration at anodic catalyst layer with the increase of the current density is a main reason of the low performance in DMFC through simulation results.

• Computational Structural Engineering
• /
• v.4 no.2
• /
• pp.87-94
• /
• 1991

In this paper, a mixed variational principle applicable to the linear elasticity of inhomogeneous anisotropic materials is presented. For derivation of the general variational principle, a systematic procedure for the variational formulation of linear coupled boundary value problems developed by Sandhu et al. is employed. Consistency condition of the field operators with the boundary operators results in explicit inclusion of boundary conditions in the governing functional. Extensions of admissible state function spaces and specialization to a certain relation in the general governing functional lead to the desired mixed variational principle. In the physical sense, the present variational principle is analogous to the Reissner's recent formulation obtained by applying Lagrange multiplier technique followed by partial Legendre transform to the classical minimum potential energy principle. However, the present one is more advantageous for the application to the general anisotropic materials since Reissner's principle contains an implicit function which is not easily converted to an explicit form.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.1B
• /
• pp.21-27
• /
• 2012

The HLLL scheme, proposed by T. Linde, determines all the wave speeds from the initial states because the middle wave is evaluated by the introduction of a generalized entropy function. The scheme is considered a genuine successor to the original HLL scheme because it is completely separated form the Roe's linearization scheme unlike the HLLE scheme and does not rely on the exact solution unlike the HLLC scheme. In this study, a numerical model was configured by the HLLL scheme with the total energy as a generalized entropy function to solve governing equations, which are the one-dimensional shallow water equations without source terms and with an additional conserved variable relating a concentration. Despite the limitations of the first order solutions, results to three cases with the exact solutions were generally accurate. The HLLL scheme appeared to be superior in comparison with the other HLL-type schemes. In particular, the scheme gave fairly accurate results in capturing the front of wetting and drying. However, it revealed shortcomings of more time-consuming calculations compared to the other schemes.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.2
• /
• pp.13-29
• /
• 1996

This paper presents calculated results of the free-surface flow around a Series 60($C_B=0.6$) model. Three-dimensional Navier-Sotkes equations are solved and Baldwin-Lomax algebraic turbulence model is adopted to simulate the high Reynolds-number flow. To reduce computational efforts, velocity components near the wall are extrapolated with a the solved by using the Implicit Approximate Factorization method[2]. The successive-over-relaxation method is used for solving pressure-Poisson equation when obtaining the pressure field projecting the divergence-free velocity field. To simulate the free-surface flows more precisely, the numerical scheme solving the equation for the kinematic boundary condition is very important. In this paper, there numerical schemes are employed and the results are compared with the available experimental data.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.4
• /
• pp.251-258
• /
• 1993

The main purpose of this paper is to present both the natural frequencies and the buckling loads of beam-columns on Winkler-type foundations. The ordinary differential equations governing the free vibrations and the buckling loads of beam-columns on Winkler-type foundation are derived as nondimensional forms. The Runge-Kutta method and Determinant Search method are used to perform the integration of the differential equations and to determine the eigenvalues(natural frequencies and buckling loads), respectively. Hinged-hinged and damped-clamped end constraints are applied in numerical examples. The relation between frequency parameter and elastic foundation parameter is presented in figure. The effects of axial loads on the natural frequencies of beam-columns on elastic foundations are investigated and the relation between buckling load parameter and elastic foundation parameter is also analyzed. The relation between foundation rested ratio and frequency parameter, buckling load parameter are investigated. The beam-columns on non-homogeneous elastic foundation are analyzed and typical mode shapes are also presented.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.5 no.6
• /
• pp.55-64
• /
• 2001

Industrial machines are sometimes exposed to the danger of earthquake. In the design of a mechanical system, this factor should be accounted for from the viewpoint of reliability to analyze a complex nonlinear structure system under random excitation is proposed. First, the actual random excitation, such as earthquake, is approximated to the corresponding Gaussian process for the statistical analysis. The modal equations of overall system are expanded sequentially. Then, the perturbed equations are synthesized into the overall system and solved in probabilistic way. Several statistical properties of a random process that are of interest in random vibration are evaluated in each substructure. Comparing with the results of the numerical simulation proved the efficiency of the proposed method.

• Journal of the Korean GEO-environmental Society
• /
• v.15 no.3
• /
• pp.41-48
• /
• 2014

Experimental study of sedimentation and self-weight consolidation has been primary research area in dredged soil. However, good quality of the dredged soil and minimum water pollution caused by the pumping of reclaimed soil require intensive study of the flow characteristics of dredged material due to dumping. In this study, continuity and the equilibrium equations for mass flow assuming single phase was derived to simulate mass flow in dredged containment area. To optimize computation and modeling time for three dimensional geometry and boundary conditions, depth integration is applied to governing equations to consider three dimensional topography of the site. Petrov-Galerkin formulation is applied in spatial discretization of governing equations. Generalized trapezoidal rule is used for time integration, and Newton iteration process approximated the solution. DG and CDG technique were used for weighting matrix in discontinuous test function in dredged flow analysis, and numerical stability was evaluated by performed a square slump simulation. A comparative analysis for numerical methods showed that DG method applied to SU / PG formulation gives minimal pseudo oscillation and reliable numerical results.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.4
• /
• pp.319-327
• /
• 2017

In this study, characteristics of seismic behaviors of offshore wind turbine systems using concrete-suction-type supporting structures are investigated. Applying hydrodynamic pressure from the surrounding sea water and interaction forces from the underlying soil to the structural system which is composed of RNA, the tower, and the supporting structure, a governing equation of the system is derived and its earthquake responses are obtained. It can be observed from the analysis results that the responses are significantly influenced by soil-structure interaction because dynamic responses for higher natural vibration modes are increased due to the flexibility of soil. Therefore, the soil-structure interaction must be taken into consideration for accurate assessment of dynamic behaviors of offshore wind turbine systems using concrete-suction-type supporting structures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.6
• /
• pp.331-337
• /
• 2018

This paper presents dynamic algorithm of the MLS(moving least squares) difference method using first order differential Approximation. The governing equations are only discretized by the first order MLS derivative approximation. The system equation consists of an assembly of the approximate function, so the shape of system equation is similar to FEM(finite element method). The CDM(central difference method) is used for time integration of dynamic equilibrium equation. The natural frequency analyses of the MLS difference method and FEM are performed, and two analysis results are compared. Also, the accuracy of the proposed numerical method is verified by displaying the dynamic analysis results together with the results by the existing second order differential approximation. In the process of assembling the first order MLS derivative approximation, the oscillation error was suppressed and the stress distribution was interpreted as relatively uniform.