Galerkin Weak Form (galerkin + weak_form)

Distribution by Scientific Domains


Selected Abstracts


Some numerical issues using element-free Galerkin mesh-less method for coupled hydro-mechanical problems

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 7 2009
Mohammad Norouz Oliaei
Abstract A new formulation of the element-free Galerkin (EFG) method is developed for solving coupled hydro-mechanical problems. The numerical approach is based on solving the two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Spatial variables in the weak form, i.e. displacement increment and pore water pressure increment, are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on a penalty method. Numerical stability of the developed formulation is examined in order to achieve appropriate accuracy of the EFG solution for coupled hydro-mechanical problems. Examples are studied and compared with closed-form or finite element method solutions to demonstrate the validity of the developed model and its capabilities. The results indicate that the EFG method is capable of handling coupled problems in saturated porous media and can predict well both the soil deformation and variation of pore water pressure over time. Some guidelines are proposed to guarantee the accuracy of the EFG solution for coupled hydro-mechanical problems. Copyright © 2008 John Wiley & Sons, Ltd. [source]


Conserving Galerkin weak formulations for computational fracture mechanics

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 12 2002
Shaofan Li
Abstract In this paper, a notion of invariant Galerkin-variational weak forms is proposed. Two specific invariant variational weak forms, the J-invariant and the L-invariant, are constructed based on the corresponding conservation laws in elasticity, one of which is the conservation of Eshelby's energy-momentum (Eshelby. Philos. Trans. Roy. Soc. 1951; 87: 12; In Solid State Physics, Setitz F, Turnbull D (eds). Academic Press: New York, 1956; 331; Rice, J. Appl. Mech. 1968; 35: 379). It is shown that the finite element solution obtained from the invariant Galerkin weak formulations proposed here can conserve the value of J-integral, or L-integral exactly. In other words, the J and L integrals of the Galerkin finite element solutions are path independent in the discrete sense. It is argued that by using the J-invariant Galerkin weak form to compute near crack-tip field in an elastic solid, one may accurately calculate the crack extension energy release rate and subsequently the stress intensity factors in numerical computations, because the flux of the energy-momentum is conserved in discrete computations. This may provide an alternative means to accurately simulate crack growth and propagation. Copyright © 2002 John Wiley & Sons, Ltd. [source]


Higher-order XFEM for curved strong and weak discontinuities

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 5 2010
Kwok Wah Cheng
Abstract The extended finite element method (XFEM) enables the accurate approximation of solutions with jumps or kinks within elements. Optimal convergence rates have frequently been achieved for linear elements and piecewise planar interfaces. Higher-order convergence for arbitrary curved interfaces relies on two major issues: (i) an accurate quadrature of the Galerkin weak form for the cut elements and (ii) a careful formulation of the enrichment, which should preclude any problems in the blending elements. For (i), we employ a strategy of subdividing the elements into subcells with only one curved side. Reference elements that are higher-order on only one side are then used to map the integration points to the real element. For (ii), we find that enrichments for strong discontinuities are easily extended to higher-order accuracy. In contrast, problems in blending elements may hinder optimal convergence for weak discontinuities. Different formulations are investigated, including the corrected XFEM. Numerical results for several test cases involving strong or weak curved discontinuities are presented. Quadratic and cubic approximations are investigated. Optimal convergence rates are achieved using the standard XFEM for the case of a strong discontinuity. Close-to-optimal convergence rates for the case of a weak discontinuity are achieved using the corrected XFEM. Copyright © 2009 John Wiley & Sons, Ltd. [source]


Wavelet Galerkin method in multi-scale homogenization of heterogeneous media

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 3 2006
Shafigh Mehraeen
Abstract The hierarchical properties of scaling functions and wavelets can be utilized as effective means for multi-scale homogenization of heterogeneous materials under Galerkin framework. It is shown in this work, however, when the scaling functions are used as the shape functions in the multi-scale wavelet Galerkin approximation, the linear dependency in the scaling functions renders improper zero energy modes in the discrete differential operator (stiffness matrix) if integration by parts is invoked in the Galerkin weak form. An effort is made to obtain the analytical expression of the improper zero energy modes in the wavelet Galerkin differential operator, and the improper nullity of the discrete differential operator is then removed by an eigenvalue shifting approach. A unique property of multi-scale wavelet Galerkin approximation is that the discrete differential operator at any scale can be effectively obtained. This property is particularly useful in problems where the multi-scale solution cannot be obtained simply by a wavelet projection of the finest scale solution without utilizing the multi-scale discrete differential operator, for example, the multi-scale analysis of an eigenvalue problem with oscillating coefficients. Copyright © 2005 John Wiley & Sons, Ltd. [source]


Strong and weak arbitrary discontinuities in spectral finite elements

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 8 2005
A. Legay
Abstract Methods for constructing arbitrary discontinuities within spectral finite elements are described and studied. We use the concept of the eXtended Finite Element Method (XFEM), which introduces the discontinuity through a local partition of unity, so there is no requirement for the mesh to be aligned with the discontinuities. A key aspect of the implementation of this method is the treatment of the blending elements adjacent to the local partition of unity. We found that a partition constructed from spectral functions one order lower than the continuous approximation is optimal and no special treatment is needed for higher order elements. For the quadrature of the Galerkin weak form, since the integrand is discontinuous, we use a strategy of subdividing the discontinuous elements into 6- and 10-node triangles; the order of the element depends on the order of the spectral method for curved discontinuities. Several numerical examples are solved to examine the accuracy of the methods. For straight discontinuities, we achieved the optimal convergence rate of the spectral element. For the curved discontinuity, the convergence rate in the energy norm error is suboptimal. We attribute the suboptimality to the approximations in the quadrature scheme. We also found that modification of the adjacent elements is only needed for lower order spectral elements. Copyright © 2005 John Wiley & Sons, Ltd. [source]


Meshfree weak,strong (MWS) form method and its application to incompressible flow problems

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 10 2004
G. R. Liu
Abstract A meshfree weak,strong (MWS) form method has been proposed by the authors' group for linear solid mechanics problems based on a combined weak and strong form of governing equations. This paper formulates the MWS method for the incompressible Navier,Stokes equations that is non-linear in nature. In this method, the meshfree collocation method based on strong form equations is applied to the interior nodes and the nodes on the essential boundaries; the local Petrov,Galerkin weak form is applied only to the nodes on the natural boundaries of the problem domain. The MWS method is then applied to simulate the steady problem of natural convection in an enclosed domain and the unsteady problem of viscous flow around a circular cylinder using both regular and irregular nodal distributions. The simulation results are validated by comparing with those of other numerical methods as well as experimental data. It is demonstrated that the MWS method has very good efficiency and accuracy for fluid flow problems. It works perfectly well for irregular nodes using only local quadrature cells for nodes on the natural boundary, which can be generated without any difficulty. Copyright © 2004 John Wiley & Sons, Ltd. [source]