• Journal of Korea Water Resources Association
• /
• v.43 no.10
• /
• pp.911-920
• /
• 2010

The closed water bodies, such as reservoirs and lakes, could be contaminated by an inflow of pollutants in the upstream as well as a stratification caused by seasonal natural phenomena. The vertical circulation particularly plays an important role in reduction of environmental pollutants. The factors of the vertical circulation are the temperature, wind, thermal diffusivity and sunlight. The wind is probably the most significant factor among them. Thus, it is necessary to describe the validation and application of a three-dimensional numerical model of wind-induced circulation in a thermally stratified flow. In this paper, a three-dimensional numerical model for the thermally stratified flows is presented. The model is conducted in three steps to calculate the velocity components from the momentum equations in x- and y- axis directions, the elevations from the free surface equation and the temperature from the scalar transport equation. Numerical predictions are compared with available analytical solutions for the sloshing free surface movement in a rectangular basin. The numerical results generally show a reasonable agreement with analytical solutions. And the model is applied to the circulation for the wind induced flow in a thermally stratification. Consequently, the developed model is validated by two verifications and phenomena of the internal flow.

• Computational Structural Engineering
• /
• v.8 no.4
• /
• pp.87-97
• /
• 1995

For the same configuration of two-dimensional finite element models, 6-node element exhibits stiffer bending stiffness than 8-node element. This is true in the relation between 16-node element and 20-node element for three-dimensional model. This stiffening phenomenon comes from the elimination of several mid nodes from full-node elements. Therefore, this may be called 'relative stiffness stiffening phenomenon'. It seems that there are a couple of ways to correct the stiffening effect, however, we could find only one effective method-the method of modification of Gauss sampling points-which passes the patch test and does not alter other kinds of stiffness, such as extensional stiffness. The quantity of modification is a function of Poisson's ratios of the constituent materials. We could obtain two modification equations, one for plane stress case and the other for plane strain case. This method can be extended to 3-dimensional solid elements. Except the exact plane strain cases, most 3-dimensional plates could be modeled successfully with 16-node element modified by the equation for the plane stress case. The effectiveness of the modification method is checked by applying it to several examples with excellent improvements. In numerical examples, beams with various boundary conditions are subjected to static and time-dependent loads. Free and forced motion analyses of beams and plates are also tested. The beam and plate may be composed of isotropic multilayers as well as a single layer.

• Journal of Advanced Marine Engineering and Technology
• /
• v.36 no.5
• /
• pp.635-644
• /
• 2012

In order to design a wave energy generating system, a 6-DOF analysis technique is applied to the three-Dimensional CFD analysis on of a floating body and the behavior is interpreted according to the nature of the incoming wave. A wave period of 5.5s & amplitude of 0.57m from Marado is chosen. 12 case of natural pitching period from 1.25 to 2.8s has been modeled. The relation between tuning factor & pitch angle for the waves generated is compared to analyze the effects of energy absorption variables, namely mass moment of inertia, angular velocity and angular acceleration. From the results obtained, we conclude that model L is the maximum power absorbed, 6kW approximately. A maximum pitch angle of 1.91 degree was attained by Model F, and the maximum displacement of nearly 0.7m was attained by Model L among models D, F and L.

• Journal of the Korean Geotechnical Society
• /
• v.32 no.9
• /
• pp.51-62
• /
• 2016

Parametric studies for various site conditions by using 3d numerical model were carried out in order to estimate dynamic behavior of soil-pile-structure system in dry soil deposits. Proposed model was analyzed in time domain using FLAC3D which is commercial finite difference code to properly simulate nonlinear response of soil under strong earthquake. Mohr-Coulomb criterion was adopted as soil constitutive model. Soil nonlinearity was considered by adopting the hysteretic damping model, and an interface model which can simulate separation and slip between soil and pile was adopted. Simplified continuum modeling was used as boundary condition to reduce analysis time. Also, initial shear modulus and yield depth were appropriately determined for accurate simulation of system's nonlinear behavior. Parametric study was performed by varying weight of superstructure, pile length, pile head fixity, soil relative density with proposed numerical model. From the results of parametric study, it is identified that inertial force induced by superstructure is dominant on dynamic behavior of soil-pile-structure system and effect of kinematic force induced by soil movement was relatively small. Difference in dynamic behavior according to the pile length and pile head fixity was also numerically investigated.

• Journal of Korean Society of Steel Construction
• /
• v.23 no.2
• /
• pp.199-210
• /
• 2011

The cables in cable-stayed bridges are under high stress and are very sensitive to vibration due to their small section areas compared with other members. Therefore, it is reasonable to evaluate the cable impact factor by taking into account the dynamic effect due to moving-vehicle motion. In this study, the cable impact factors were evaluated via moving-vehicle-load analysis, considering the design parameters, i.e., vehicle weight, cable model, road surface roughness, vehicle speed, longitudinal distance between vehicles. For this purpose, two steel composite cable-stayed bridges with 230- and 540-m main spans were selected. The results of the analysis were then compared with those of the influence line method that is currently being used in design practice. The road surface roughness was randomly generated based on ISO 8608, and the convergence of impact factors according to the number of generated road surfaces was evaluated to improve the reliability of the results. A9-d.o.f. tractor-trailer vehicle was used, and the vehicle motion was derived from Lagrange's equation. 3D finite element models for the selected cable-stayed bridges were constructed with truss elements having equivalent moduli for the cables, and with beam elements for the girders and the pylons. The direct integration method was used for the analysis of the bridge-vehicle interaction, and the analysis was conducted iteratively until the displacement error rate of the bridge was within the specified tolerance. It was acknowledged that the influence line method, which cannot consider the dynamic effect due to moving-vehicle motion, could underestimate the impact factors of the end-cables at the side spans, unlike moving-vehicle-load analysis.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.7 no.3
• /
• pp.69-78
• /
• 2007

In this study, a time-dependent aspect of an embankment failure is considered to simulate a flood inundation map and calculate overtopping discharge induced by an embankment failure. A numerical model has been developed by solving the two dimensional nonlinear shallow water equations with a finite volume method on unstructured grids. To analyze a Riemann problem, the HLLC approximate Riemann solver and the Weighted Averaged Flux method are employed by using a TVD limiter and the source term treatment is also employed by using the operator splitting method. Firstly, the numerical model is applied to a dam break problem and a sloping seawall. Obtained numerical results show good agreements with experimental data. Secondly, the model is applied to a flow induced by an embankment failure by assuming that the width and elevation of embankment are varied with time-dependent functions. As a result of the comparison with each numerical overtopping discharge, established flood inundation discharges in the previous studies are overestimated than the result of the present numerical model.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.559-573
• /
• 2008

In this study, we newly proposed 3-D nonlinear wave driver utilizing the Navier-Stokes Eq. the numerical integration of which is carried out using SPH (Smoothed Particle Hydrodynamics), an internal wave generation with the source function of Gaussian distribution and an energy absorbing layer. For the verification of new 3-D nonlinear wave driver, we numerically simulate the sloshing problem within a parabolic water basin triggered by a Gaussian hump and uniformly inclined water surface by Thacker (1981). It turns out that the qualitative behavior of sloshing caused by relaxing the external force which makes a free surface convex or uniformly inclined is successfully simulated even though phase error is visible and an inundation height shrinks as numerical simulation more proceeds. For the more severe test, we also simulate the nonlinear shoaling and refraction over uniform beach of wedge shape. It is shown that numerically simulated waves are less refracted than the linear counterpart by Hamiltonian ray theory due to nonlinearity, energy dissipation at the bottom and side walls, energy loss induced by breaking, and the hydraulic jump occurring when breaking waves encounter a down-rush by the preceding wave.