• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.409-413
• /
• 2005

Even though the relative importance of length scale of flow system allow us to simplify three dimensional flow problem to one or two dimensional representation, many systems still require three dimensional analysis. The objective of this study is to develop an efficient and accurate finite element model for analyzing and predicting three dimensional flow features in natural rivers and to offend to model spreading of pollutants and transport of sediments in the future. Firstly, three dimensional Reynolds averaged Navier-Stokes equations with the hydrostatic pressure assumption in generalized curvilinear coordinates were combined with the kinematic free-surface condition. Secondly. to simulate realistic high Reynolds number flow, the model employed the Streamline Upwind/Petrov-Galerkin(SU/PG) scheme as a weighting function for the finite element method in conjunction with an appropriate turbulence model(Smagorinsky scheme for the horizontal plain and Mellor-Yamada scheme for the vertical direction). Several tests is performed for the purpose of validation and verification of the developed model. A simple rectangular channel, 5-shaped and U-shaped channel are used for tests and comparisons are made with RMA-10 model. Runs for each case is converged stably without a oscillation and calculated water-surface deformation, longitudinal and transversal velocities, and velocity vector fields are in good agreement with the results of RMA-10 model.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.3B
• /
• pp.175-183
• /
• 2012

In this study, Galerkin scheme and SU/PG scheme of Petrov-Galerkin family were applied to the shallow water equations and a finite element model for shallow water flow was developed. Numerical simulations were conducted in several flumes with convection-dominated flow condition. Flow simulation of channel with slender structure in the water course revealed that Galerkin and SU/PG schemes showed similar results under very low Fr number and Re number condition. However, when the Fr number increased up to 1.58, Galerkin scheme did not converge while SU/PG scheme produced stable solutions after 5 iterations by Newton-Raphson method. For the transcritical flow simulation in diverging channel, the present model predicted the hydraulic jump accurately in terms of the jump location, the depth slope, and the flow depth after jump, and the numerical results showed good agreements with the hydraulic experiments carried out by Khalifa(1980). In the oblique hydraulic jump simulation, in which convection-dominated supercritical flow (Fr=2.74) evolves, Galerkin scheme blew up just after the first iteration of the initial time step. However, SU/PG scheme captured the boundary of oblique hydraulic jump accurately without numerical oscillation. The maximum errors quantified with exact solutions were less than 0.2% in water depth and velocity calculations, and thereby SU/PG scheme predicted the oblique hydraulic jump phenomena more accurately compared with the previous studies (Levin et al., 2006; Ricchiuto et al., 2007).

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2005.05a
• /
• pp.43-46
• /
• 2005

Micro molding technology is a promising mass production technology for polymer based microstructures. Mass production technologies such as the micro injection/compression molding, hot embossing, and micro reaction molding are already in use. In the present study, we have developed a numerical analysis system to simulate three-dimensional non-isothermal cavity filling for hot embossing, with a special emphasis on the free surface capturing. Precise free surface capturing has been successfully accomplished with the level set method, which is solved by means of the Runge-Kutta discontinuous Galerkin (RKDG) method. The RKDG method turns out to be excellent from the viewpoint of both numerical stability and accuracy of volume conservation. The Stokes equations are solved by the stabilized finite element method using the equal order tri-linear interpolation function. To prevent possible numerical oscillation in temperature Held we employ the streamline upwind Petrov-Galerkin (SUPG) method. With the developed code we investigated the detailed change of free surface shape in time during the mold filling. In the filling simulation of a simple rectangular cavity with repeating protruded parts, we find out that filling patterns are significantly influenced by the geometric characteristics such as the thickness of base plate and the aspect ratio and pitch of repeating microstructures. The numerical analysis system enables us to understand the basic flow and material deformation taking place during the cavity filling stage in microstructure fabrications.

• Steel and Composite Structures
• /
• v.35 no.1
• /
• pp.77-92
• /
• 2020

The present paper outlined a procedure for geometrically nonlinear dynamic analysis of functionally graded graphene platelets-reinforced (GPLR-FG) nanocomposite cylinder subjected to mechanical shock loading. The governing equation of motion for large deformation problems is derived using meshless local Petrov-Galerkin (MLPG) method based on total lagrangian approach. In the MLPG method, the radial point interpolation technique is employed to construct the shape functions. A micromechanical model based on the Halpin-Tsai model and rule of mixture is used for formulation the nonlinear functionally graded distribution of GPLs in polymer matrix of composites. Energy dissipation in analyses of the structure responding to dynamic loads is considered using the Rayleigh damping. The Newmark-Newton/Raphson method which is an incremental-iterative approach is implemented to solve the nonlinear dynamic equations. The results of the proposed method for homogenous material are compared with the finite element ones. A very good agreement is achieved between the MLPG and FEM with very fine meshing. In addition, the results have demonstrated that the MLPG method is more effective method compared with the FEM for very large deformation problems due to avoiding mesh distortion issues. Finally, the effect of GPLs distribution on strength, stiffness and dynamic characteristics of the cylinder are discussed in details. The obtained results show that the distribution of GPLs changed the mechanical properties, so a classification of different types and volume fraction exponent is established. Indeed by comparing the obtained results, the best compromise of nanocomposite cylinder is determined in terms of mechanical and dynamic properties for different load patterns. All these applications have shown that the present MLPG method is very effective for geometrically nonlinear analyses of GPLR-FG nanocomposite cylinder because of vanishing mesh distortion issue in large deformation problems. In addition, since in proposed method the distributed nodes are used for discretization the problem domain (rather than the meshing), modeling the functionally graded media yields to more accurate results.

• Journal of Korea Water Resources Association
• /
• v.35 no.5
• /
• pp.463-474
• /
• 2002

A two-dimensional water quality management model, Unsteady/Uncertainty Water Quality Model(UUWQM), is developed for a hydrodynamic analysis, an advection-diffusion analysis, and a reliability analysis by using uncertainty technique. The model is applied to the 35 km reach of Sungju to Hyunpoong in the midstream of Nakdong River. 2-D hydrodynamic and water quality analyses are peformed in this reach. Important input variables are decided by sensitivity analysis and verified by Monte Carlo method. Frequency distributions of water quality concentrations are computed from MFOSM method and Monte Carlo method at several locations in this study area. A water quality management system is constructed by calculating the violation probabilities of existing water quality standards.

• Journal of Wetlands Research
• /
• v.13 no.3
• /
• pp.395-409
• /
• 2011

The objective of this study is to develop an accurate and robust two-dimensional finite element method for turbulence simulation in open channels. The model is based on Streamline Upwind/Petrov-Galerkin finite element method and Boussinesq's eddy viscosity theory. The method developed in the study is depth-averaged mixing length model which assumes anisotropic and local equilibrium state of turbulence. The model calibration and validation were performed by comparing with analytical solutions and observed data. Several numerical simulations were carried out, which examined the performance of the turbulence model for the purpose of sensitivity analysis. The uniform channels that appear horizontal flow and vertical flow were carried out. The model was also applied to the Han river was in for the applicability test. The results were compared with the observed data. The suggested model displayed reasonable flow distribution compare to the observed data in natural river flow. As a result of this study, the two-dimensional finite element model provides a reliable results for flow distribution based on the turbulence simulation in open channels.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.414-418
• /
• 2005

A finite element model is developed for the numerical simulation of bed elevation change in a curved channel. The SU/PG (Streamline-Upwind/Petrov-Galerkin) method is used to solve 2D shallow water equations and the BG (Bubnov-Galerkin) method is used for the Exner equation. For the time derivative terms, the Crank-Nicolson scheme is used. The developed model is a decoupled model in a sense that the bed elevation does not change simultaneously with the flow during the computational time step. The total load formula with is used for the sediment transport model. The slip conditions are described along the lateral boundaries. The effects of gravity force due to geometry change and the secondary flows in a curved channel are considered in the model. For the verification, the model is applied to two laboratory experiments. The first is $140^{\circ}$ bended channel data at Delft Hydraulics Laboratory and the second is $140^{\circ}$ bended channel data at Laboratory of Fluid Mechanics of the Delft University of Technology. The finite element grid is constructed with linear quadrilateral elements. It is found that the computed results are in good agreement with measured data, showing a point bar at the inner bank and a pool at the outer bank.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.10
• /
• pp.754-760
• /
• 2008

Computation of moving interface by the level set method typically requires the reinitialization of level set function. An inaccurate estimation of level set function $\phi$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, an accurate and robust reinitialization process is essential to the simulation of free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates level set function directly using a normal vector on the interface without solving there-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1 splitting/SUPG (Streamline Upwind Petrov-Galerkin) FEM are adopted to discretize advection equation of the level set function and the incompressible Navier-Stokes equation, respectively. Advection equation and re-initialization process of free surface capturing are validated with benchmark problems, i.e., a broken dam flow and timereversed single vortex flow. The simulation results are in good agreement with the existing results.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1999.03b
• /
• pp.253-256
• /
• 1999

The Streamline Upwind Petrov-Galerkin(SUPG) finite element method is used to solve the two-dimensional laminar and turbulent flow. The flow is simulated by averaged Navier-Stokes equations with a penalty function approach and the lograithmic(k-$\varepsilon$) turbulent model is employed to take into account its turbulent behavior. The near-wall viscous sub-layer model is employed to approach the dominant viscous effects in the near wall zones. To find a good-enough initial guess of the Newton-Raphson iteration solving Nonlinear Matrix the Incremental method is considered for momentum and the Incomplete logarithmic turbu-lent equations for Turbulence. The validation of our method is investigated in comparision with published experimental data.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.10a
• /
• pp.375-382
• /
• 2003

Porous media consist of physically and chemically different materials and have an extremely complicated behavior due to the different material properties of each of its constituents. In addition, the internal structure of porous media has generally a complex geometry that makes the description of its mechanical behavior quite complex. Thus, in order to describe and clarify the deformation behavior of porous media, constitutive models for deformation of porous media coupling several effects such as flow of fluids of thermodynamical change need to be developed in frame of Arbitrary Lagrangian Eulerian (ALE) description. The aim of ALE formulations is to maximize the advantages of Lagrangian and Eulerian methods, and to minimize the disadvantages. Therefore, this method is appropriate for the analysis of porous media that are considered for the behavior of solids and fluids. First of all, governing equations for saturated porous media based on ALE description are derived. Then, weak forms of these equations are obtained in order to implement numerical method using finite element method. Finally, Petrov-Galerkin method Is applied to develop finite element formulation.