• Journal of Korea Water Resources Association
• /
• v.31 no.5
• /
• pp.599-612
• /
• 1998

Numerical analysis of convection-dominated transport problems are challenging because of dual characteristics of the governing equation. In the finite element method, a strategy is to modify the test function to weight more in the upwind direction. This is called as the Petrov-Galerkin method. In this paper, both N+1 and N+2 Petrov-Galerkin methods are applied to transport problems at high grid Peclet number. Frequency fitting algorithm is used to obtain optimal levels of N+2 upwinding, and the results are discussed. Also, a new Petrov-Galerkin method, named as "Optimal Test Function Petrov-Galerkin Method," is proposed in this paper. The test function of this numerical method changes its shape depending upon relative strength of the convection to the diffusion. A numerical experiment is carried out to demonstrate the performance of the proposed method.

• Water for future
• /
• v.29 no.6
• /
• pp.155-166
• /
• 1996

A finite element model for simulating gradually and rapidly varied unsteady flow in open channel is developed based on dynamic wave equation using Petrov-Galerkin method. A matrix stability analysis shows the selective damping of short wave lengths and excellent phase accuracies achived by Petrov-Galerkin method. Whereas the Preissmann scheme displays less selective damping and poor phase accuracies, and Bubnov-Galerkin method shows nondissipative characteristics whicn causes a divergence problem in short wave length. The analysis also shows that the Petrov-Galerkin method displays the desirable combination of selective damping of high frequency progressive waves over a wide range of Courant number and good phase accuracy at low Courant number. Therefore, the Petrov-Galerkin can be effectively applied to gradually and rapidly varied unsteady flow.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.103-111
• /
• 2005

In this paper, a new meshfree technique which improves the numerical integration accuracy is introduced. This new method called thc Petrov-Galerkin natural clement method(PG-NEM) by authors is based on the Voronoi diagram and the Delaunay triangulation which is based on the same concept used lot conventional natural clement method called the Bubnov-Galerkin natural element method(BG-NEM). But, unlike the BG-NEM, the test basis function is differently chosen, based on the concept of Petrov-Galerkin, such that its support coincides exactly with a regular integration region in background mesh. Therefore, it is expected that the proposed technique ensures the remarkably improved numerical integration accuracy in comparison with the BG-NEM.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.6 s.237
• /
• pp.703-710
• /
• 2005

A finite element discretization of the advection and redistancing equations of level set method has been studied. It has been shown that Galerkin spatial discretization combined with Crank-Nicolson temporal discretization of the advection equation of level set yields a good result and that consistent streamline upwind Petrov-Galerkin(CSUPG) discretization of the redistancing equation gives satisfactory solutions for two test problems while the solutions of streamline upwind Petrov-Galerkin(SUPG) discretization are dissipated by the numerical diffusion added for the stability of a hyperbolic system. Furthermore, it has been found that the solutions obtained by CSUPG method are comparable to those by second order ENO 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.10
• /
• pp.2667-2677
• /
• 1995

A new numerical algorithm using equal order linear finite element and fractional step method has been developed which is capable of analyzing unsteady fluid flow and heat transfer problems. Streamline Upwind Petrov-Galerkin (SUPG) method is used for the weighted residual formulation of the Navier-Stokes equations. It is shown that fractional step method, in which pressure term is splitted from the momentum equation, reduces computer memory and computing time. In addition, since pressure equation is derived without any approximation procedure unlike in the previously developed SIMPLE algorithm based FEM codes, the present numerical algorithm gives more accurate results than them. The present algorithm has been applied preferentially to the well known bench mark problems associated with steady flow and heat transfer, and proves to be more efficient and accurate.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1999.05a
• /
• pp.547-552
• /
• 1999

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.113-121
• /
• 2005

In order to resolve a common numerical integration inaccuracy of meshfree methods, we introduce an improved natural clement method called Petrov-Galerkin natural element method(PG-NEM). While Laplace basis function is being taken for the trial shape function, the test shape function in the present method is differently defined such that its support becomes a union of Delaunay triangles. This approach eliminates the inconsistency of tile support of integrand function with the regular integration domain, and which preserves both simplicity and accuracy in the numerical integration. In this paper, the validity of the PG-NEM is verified through the representative benchmark problems in 2-d linear elasticity. For the comparison, we also analyze the problems using the conventional Bubnov-Galerkin natural element method(BG-NEM) and constant strain finite clement method(CS-FEM). From the patch test and assessment on convergence rate, we can confirm the superiority of the proposed meshfree method.

• 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).

• Journal of Korea Water Resources Association
• /
• v.32 no.3
• /
• pp.237-246
• /
• 1999

A finite element model is studied to simulate unsteady free surface flow based on dynamic wave equation and collocation method. The collocation method is used in conjunction with Hermite polynomials, and resulting matrix equations are solved by skyline method. The model is verified by applying to hydraulic jump, nonlinear disturbance propagation and dam-break flow in a horizontal frictionless channel. The computed results are compared with those by Bubnov-Galerkin and Petrov-Galerkin methods. It is also applied to the North Han River to simulate the floodwave propagation. The computed results have good agreements with those of DWOPER model in terms of discharge hydrographs. The suggested model has proven to be one of the promising scheme for simulating the gradually and rapidly varied unsteady flow in open channels.