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

Title: A locally conservative stress recovery technique for continuous Galerkin FEM in linear elasticity

; ;
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
FE0004832; SC0004982; WYDEQ49811GNTG
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Computer Methods in Applied Mechanics and Engineering
Additional Journal Information:
Journal Volume: 286; Journal Issue: C; Related Information: CHORUS Timestamp: 2017-06-23 15:28:08; Journal ID: ISSN 0045-7825
Country of Publication:

Citation Formats

Bush, Lawrence, Deng, Quanling, and Ginting, Victor. A locally conservative stress recovery technique for continuous Galerkin FEM in linear elasticity. Netherlands: N. p., 2015. Web. doi:10.1016/j.cma.2015.01.002.
Bush, Lawrence, Deng, Quanling, & Ginting, Victor. A locally conservative stress recovery technique for continuous Galerkin FEM in linear elasticity. Netherlands. doi:10.1016/j.cma.2015.01.002.
Bush, Lawrence, Deng, Quanling, and Ginting, Victor. Wed . "A locally conservative stress recovery technique for continuous Galerkin FEM in linear elasticity". Netherlands. doi:10.1016/j.cma.2015.01.002.
title = {A locally conservative stress recovery technique for continuous Galerkin FEM in linear elasticity},
author = {Bush, Lawrence and Deng, Quanling and Ginting, Victor},
abstractNote = {},
doi = {10.1016/j.cma.2015.01.002},
journal = {Computer Methods in Applied Mechanics and Engineering},
number = C,
volume = 286,
place = {Netherlands},
year = {Wed Apr 01 00:00:00 EDT 2015},
month = {Wed Apr 01 00:00:00 EDT 2015}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/j.cma.2015.01.002

Citation Metrics:
Cited by: 4works
Citation information provided by
Web of Science

Save / Share:
  • A discontinuous Galerkin finite element method (DGFEM) has been developed and tested for the linear, three-dimensional, rotating incompressible Euler equations. These equations admit complicated wave solutions, which poses numerical challenges. These challenges concern: (i) discretisation of a divergence-free velocity field; (ii) discretisation of geostrophic boundary conditions combined with no-normal flow at solid walls; (iii) discretisation of the conserved, Hamiltonian dynamics of the inertial-waves; and, (iv) large-scale computational demands owing to the three-dimensional nature of inertial-wave dynamics and possibly its narrow zones of chaotic attraction. These issues have been resolved, for example: (i) by employing Dirac’s method of constrained Hamiltonian dynamicsmore » to our DGFEM for linear, compressible flows, thus enforcing the incompressibility constraints; (ii) by enforcing no-normal flow at solid walls in a weak form and geostrophic tangential flow along the wall; and, (iii) by applying a symplectic time discretisation. We compared our simulations with exact solutions of three-dimensional incompressible flows, in (non) rotating periodic and partly periodic cuboids (Poincaré waves). Additional verifications concerned semi-analytical eigenmode solutions in rotating cuboids with solid walls. Finally, a simulation in a tilted rotating tank, yielding more complicated wave dynamics, demonstrates the potential of our new method.« less
  • A linear discontinuous Galerkin treatment is developed for the continuous slowing down operator. Multigroup Legendre coefficients are derived that allow this treatment to be implemented in standard S/sub ETA/ codes. A synthetic acceleration scheme is developed to accelerate the outer iterations required by this treatment. A computational comparison between the diamond difference treatment and the linear discontinuous treatment is given.
  • The effect of the nonlinear elasticity of carbon fibers on the axial Young's modulus of a uniaxial composite is reported. Previous work on a metal matrix composite yielded an equation to predict the axial Young's modulus and its change with temperature. The equation takes into account the thermal strains generated in the fiber and the matrix by the mismatch of thermal expansion, the nonlinear elasticity of the constituents, and the normal effects of temperature and applied stress. Nonlinear elasticity data for pitch-based carbon fibers are used together with data for an aluminum matrix and the predictive equation to calculate Young'smore » modulus. The results are compared with experimental data and the simpler Rule-of-Mixtures (ROM) predictions which neglect the effects of thermally generated strains coupled with nonlinear elasticity. The results display fairly close agreement with measured values of composite modulus, but large differences from the simpler ROM predictions. Although the validity of the predictive equation has not been fully verified, the variation of composite modulus with temperature and applied stress which it predicts is quite significant. Under certain conditions, a composite can be designed with a modulus that is independent of temperature or applied stress. 17 refs.« less
  • A general framework for implementing numerical boundary conditions, based upon the rigorous application of characteristic theory, has been developed for the two-step Taylor–Galerkin FEM scheme. The method consists of solving the compatibility equations for the temporal change in characteristic variables using the first step of the two-step Taylor–Galerkin FEM. This application of characteristic boundary conditions is consistent with the spatial and temporal discretization of the two-step Taylor–Galerkin FEM. It is ideal for domains discretized with linear unstructured finite elements as time and space extrapolation from interior elements is not required. Boundary conditions are constructed from the characteristic solution for solidmore » wall and symmetry boundaries, flow exits, and flow inlets. Two simulations are shown highlighting the performance of these boundary conditions.« less