skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Generalization of mixed multiscale finite element methods with applications

Abstract

Many science and engineering problems exhibit scale disparity and high contrast. The small scale features cannot be omitted in the physical models because they can affect the macroscopic behavior of the problems. However, resolving all the scales in these problems can be prohibitively expensive. As a consequence, some types of model reduction techniques are required to design efficient solution algorithms. For practical purpose, we are interested in mixed finite element problems as they produce solutions with certain conservative properties. Existing multiscale methods for such problems include the mixed multiscale finite element methods. We show that for complicated problems, the mixed multiscale finite element methods may not be able to produce reliable approximations. This motivates the need of enrichment for coarse spaces. Two enrichment approaches are proposed, one is based on generalized multiscale finte element metthods (GMsFEM), while the other is based on spectral element-based algebraic multigrid (rAMGe). The former one, which is called mixed GMsFEM, is developed for both Darcy’s flow and linear elasticity. Application of the algorithm in two-phase flow simulations are demonstrated. For linear elasticity, the algorithm is subtly modified due to the symmetry requirement of the stress tensor. The latter enrichment approach is based on rAMGe. Themore » algorithm differs from GMsFEM in that both of the velocity and pressure spaces are coarsened. Due the multigrid nature of the algorithm, recursive application is available, which results in an efficient multilevel construction of the coarse spaces. Stability, convergence analysis, and exhaustive numerical experiments are carried out to validate the proposed enrichment approaches. iii« less

Authors:
 [1]
  1. Texas A & M Univ., College Station, TX (United States)
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1260492
Report Number(s):
LLNL-TH-696037
DOE Contract Number:
AC52-07NA27344
Resource Type:
Thesis/Dissertation
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE

Citation Formats

Lee, C S. Generalization of mixed multiscale finite element methods with applications. United States: N. p., 2016. Web. doi:10.2172/1260492.
Lee, C S. Generalization of mixed multiscale finite element methods with applications. United States. doi:10.2172/1260492.
Lee, C S. Mon . "Generalization of mixed multiscale finite element methods with applications". United States. doi:10.2172/1260492. https://www.osti.gov/servlets/purl/1260492.
@article{osti_1260492,
title = {Generalization of mixed multiscale finite element methods with applications},
author = {Lee, C S},
abstractNote = {Many science and engineering problems exhibit scale disparity and high contrast. The small scale features cannot be omitted in the physical models because they can affect the macroscopic behavior of the problems. However, resolving all the scales in these problems can be prohibitively expensive. As a consequence, some types of model reduction techniques are required to design efficient solution algorithms. For practical purpose, we are interested in mixed finite element problems as they produce solutions with certain conservative properties. Existing multiscale methods for such problems include the mixed multiscale finite element methods. We show that for complicated problems, the mixed multiscale finite element methods may not be able to produce reliable approximations. This motivates the need of enrichment for coarse spaces. Two enrichment approaches are proposed, one is based on generalized multiscale finte element metthods (GMsFEM), while the other is based on spectral element-based algebraic multigrid (rAMGe). The former one, which is called mixed GMsFEM, is developed for both Darcy’s flow and linear elasticity. Application of the algorithm in two-phase flow simulations are demonstrated. For linear elasticity, the algorithm is subtly modified due to the symmetry requirement of the stress tensor. The latter enrichment approach is based on rAMGe. The algorithm differs from GMsFEM in that both of the velocity and pressure spaces are coarsened. Due the multigrid nature of the algorithm, recursive application is available, which results in an efficient multilevel construction of the coarse spaces. Stability, convergence analysis, and exhaustive numerical experiments are carried out to validate the proposed enrichment approaches. iii},
doi = {10.2172/1260492},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Mon Aug 01 00:00:00 EDT 2016},
month = {Mon Aug 01 00:00:00 EDT 2016}
}

Thesis/Dissertation:
Other availability
Please see Document Availability for additional information on obtaining the full-text document. Library patrons may search WorldCat to identify libraries that hold this thesis or dissertation.

Save / Share:
  • New finite-element methods are proposed for mixed variational formulations. The methods are constructed by adding to the classical Galerkin method various least-squares like terms. The additional terms involve integrals over element interiors, and include mesh-parameter dependent coefficients. The methods are designed to enhance stability. Consistency is achieved in the sense that exact solutions identically satisfy the variational equations.Applied to several problems, simple finite-element interpolations are rendered convergent, including convenient equal-order interpolations generally unstable within the Galerkin approach. The methods are subdivided into two classes according to the manner in which stability is attained: (1) circumventing Babuska-Brezzi condition methods; (2) satisfyingmore » Babuska-Brezzi condition methods. Convergence is established for each class of methods. Applications of the first class of methods to Stokes flow and compressible linear elasticity are presented. The second class of methods is applied to the Poisson, Timoshenko beam and incompressible elasticity problems. Numerical results demonstrate the good stability and accuracy of the methods, and confirm the error estimates.« less
  • This dissertation contains the basic theory and a systematic procedure for constructing a least-squares finite element method for linear and nonlinear first-order systems of partial differential equations in engineering and mathematical physics. Error estimates, a condition number bound, and the optimal choice for weighting factors for the least squares finite element approximation are presented. The procedure is implemented in conjunction with an element-by-element preconditioned conjugate-gradient method. A new preconditioner is formulated and analyzed. An adaptive refinement strategy is developed and implemented in a series of numerical studies. Applications to practical problems in incompressible and compressible flow analysis are presented.
  • The author considers an approximate scheme for treating a class of two dimensional, singularly perturbed, exterior boundary value problems. The scheme is a combination of boundary element and finite element methods. By introducing a factitious boundary the region is divided into two regions, an inner and outer region. In the inner region, the finite element method is employed for solving the field equations under consideration with proper natural boundary conditions. These natural boundary conditions are obtained by solving the corresponding boundary integral equations. In the outer regions, solutions are expressed in terms of simple layer potentials which are matched tomore » the finite element solution on the factitious boundary. Error estimates and numerical experiments are included.« less
  • A mixed method, which combines a nonlinear conjugate gradient and Newton's methods, has proved effective for solution of three-dimensional nonlinear finite-element models. The class of problems to which the methods have been applied consists of finite deformations of nearly incompressible hyperelastic solids. Incrementation of a penalty parameter allows the effective use of conjugate-gradient methods, in contrast to the failure of search methods when Lagrange multipliers are employed. The global nature of the nonlinear conjugate-gradient scheme has proven effective in approaching solution points far removed from initial estimates of the solution. On the other hand, the quadratic convergence rate of Newton'smore » method is utilized, to good effect, once the solution has been closely approximated by the conjugate-gradient calculation. The use of vector computers (e.g., Cray X-MP, Cyber 205) was considered from the beginning in the design of the algorithms and the computer code. Comparisons between global matrix and element by element algorithms show that the element-by-element algorithm is better than the global-matrix algorithm when quadratic or higher-order elements are used. Comparison of the method presented here with other iterative schemes, made using ITPACK, demonstrate the effectiveness of this method.« less
  • Residual microstresses arising from the thermal expansion mismatch between phases in WC-Ni composites with particle volume fractions of 0.11, 0.21, 0.41 and 0.60 are studied using neutron diffraction and nonlinear plane stress finite element (FE) analysis. A method to interpret strain distributions from measurements of line broadening using FE results is developed. This method utilizes diffraction data to validate FE models which are then used to interpret the microstructural details of the stress distribution. The variation in composite microstress is confirmed as a major source of diffraction peak broadening in WC composites. The FE analysis indicates that models based onmore » composite microstructures, as opposed to models consisting of regular arrays of particles, are necessary to properly model stress distributions. A means to approximate the breadth of the stress distribution from strain distributions measured with diffraction is presented. The effect of the thermal residual stress results on the stress-strain response of the material are considered.« less